Activation of the beta 2-adrenergic receptor involves disruption of an ionic lock between the cytoplasmic ends of transmembrane segments 3 and 6

Comparison of the function of two novel human dopamine D2 receptor variants identifies a likely mechanism for their pathogenicity

Two recently discovered DRD2 mutations, c.634A > T, p.Ile212Phe and c.1121T > G, p.Met374Arg, cause hyperkinetic movement disorders that have overlapping features but apparently differ in severity. The two known carriers of the Met374Arg variant had early childhood disease onset and more severe motor, cognitive, and neuropsychiatric deficits than any known carriers of the Ile212Phe variant, whose symptoms were first apparent in adolescence. Here, we evaluated if differences in the function of the two variants in cultured cells could explain differing pathogenicity. Both variants were expressed less abundantly than the wild type receptor and exhibited loss of agonist-induced arrestin binding, but differences in expression and arrestin binding between the variants were minor. Basal and agonist-induced activation of heterotrimeric Gi/o/z proteins, however, showed clear differences; agonists were generally more potent at Met374Arg than at the Ile212Phe or wild type variants. Furthermore, all Gα subtypes tested were constitutively activated more by Met374Arg than by Ile212Phe. Met374Arg produced greater constitutive inhibition of cyclic AMP accumulation than Ile212Phe or the wild type D2 receptor. Met374Arg and Ile212Phe were more sensitive to thermal inactivation than the wild type D2 receptor, as reported for other constitutively active receptors, but Ile212Phe was affected more than Met374Arg. Additional pharmacological characterization suggested that the mutations differentially affect the shape of the agonist binding pocket and the potency of dopamine, norepinephrine, and tyramine. Molecular dynamics simulations provided a structural rationale for enhanced constitutive activation and agonist potency. Enhanced constitutive and agonist-induced G protein-mediated signaling likely contributes to the pathogenicity of these novel variants.

Pharmacological characterization of seven human histamine H3 receptor isoforms

The histamine H3 receptor (H3R) regulates as a presynaptic G protein-coupled receptor the release of histamine and other neurotransmitters in the brain, and is consequently a potential therapeutic target for neuronal disorders. The human H3R encodes for seven splice variants that vary in the length of intracellular loop 3 and/or the C-terminal tail but are all able to induce heterotrimeric Gi protein signaling. The last two decades H3R drug discovery and lead optimization has been exclusively focused on the 445 amino acids-long reference isoform H3R-445. In this study, we pharmacologically characterized for the first time all seven H3R isoforms by determining their binding affinities for reference histamine H3 receptor agonists and inverse agonists. The H3R-453, H3R-415, and H3R-413 isoforms display similar binding affinities for all ligands as the H3R-445. However, increased agonist binding affinities were observed for the three shorter isoforms H3R-329, H3R-365, and H3R-373, whereas inverse agonists such as the approved anti-narcolepsy drug pitolisant (Wakix®) displayed significantly decreased binding affinities for the latter two isoforms. This opposite change in binding affinity of agonist versus inverse agonists on H3R-365 and H3R-373 is associated with their higher constitutive activity in a cAMP biosensor assay as compared to the other five isoforms. The observed differences in pharmacology between longer and shorter H3R isoforms should be considered in future drug discovery programs.

Unveiling Switching Function of Amino Acids in Proteins Using a Machine Learning Approach

Dynamics of individual amino acids play key roles in the overall properties of proteins. However, the knowledge of protein structural features at the residue level is limited due to the current resolutions of experimental and computational techniques. To address this issue, we designed a novel machine-learning (ML) framework that uses Molecular Dynamics (MD) trajectories to identify the major conformational states of individual amino acids, classify amino acids switching between two distinct modes, and evaluate their degree of dynamic stability. The Random Forest model achieved 96.94% classification accuracy in identifying switch residues within proteins. Additionally, our framework distinguishes between the stable switch (SS) residues, which remain stable in one angular state and jump once to another state during protein dynamics, and unstable switch (US) residues, which constantly fluctuate between the two angular states. This study also illustrates the correlation between the dynamics of SS residues and the protein's global properties.

GPCR molecular dynamics forecasting using recurrent neural networks

G protein-coupled receptors (GPCRs) are a large superfamily of cell membrane proteins that play an important physiological role as transmitters of extracellular signals. Signal transmission through the cell membrane depends on conformational changes in the transmembrane region of the receptor, which makes the investigation of the dynamics in these regions particularly relevant. Molecular dynamics (MD) simulations provide a wealth of data about the structure, dynamics, and physiological function of biological macromolecules by modelling the interactions between their atomic constituents. In this study, a Recurrent and Convolutional Neural Network (RNN) model, namely Long Short-Term Memory (LSTM), is used to predict the dynamics of two GPCR states and three specific simulations of each one, through their activation path and focussing on specific receptor regions. Active and inactive states of the GPCRs are analysed in six scenarios involving APO, Full Agonist (BI 167107) and Partial Inverse Agonist (carazolol) of the receptor. Four Machine Learning models with increasing complexity in terms of neural network architecture are evaluated, and their results discussed. The best method achieves an overall RMSD lower than 0.139 Å and the transmembrane helices are the regions showing the minimum prediction errors and minimum relative movements of the protein.

A review of the molecular mechanisms of acaricide resistance in mites and ticks

The Arachnida subclass of Acari comprises many harmful pests that threaten agriculture as well as animal health, including herbivorous spider mites, the bee parasite Varroa, the poultry mite Dermanyssus and several species of ticks. Especially in agriculture, acaricides are often used intensively to minimize the damage they inflict, promoting the development of resistance. Beneficial predatory mites used in biological control are also subjected to acaricide selection in the field. The development and use of new genetic and genomic tools such as genome and transcriptome sequencing, bulked segregant analysis (QTL mapping), and reverse genetics via RNAi or CRISPR/Cas9, have greatly increased our understanding of the molecular genetic mechanisms of resistance in Acari, especially in the spider mite Tetranychus urticae which emerged as a model species. These new techniques allowed to uncover and validate new resistance mutations in a larger range of species. In addition, they provided an impetus to start elucidating more challenging questions on mechanisms of gene regulation of detoxification associated with resistance.

The dynamics of agonist-β2-adrenergic receptor activation induced by binding of GDP-bound Gs protein

There is considerable uncertainty about the mechanism by which the β₂-adrenergic receptor (β₂AR) is activated. Here we use molecular metadynamics computations to predict the mechanism by which an agonist induces the activation of the β₂AR and its cognate Gs protein. We found that binding agonist alone to the inactive β₂AR does not break the ionic lock and hence does not drive the β₂AR towards the activated conformation. However, we found that attaching the inactive Gs protein to the agonist-bound inactive β₂AR (containing the ionic lock) leads to partial insertion of Gαs-α5 into the core of β₂AR, which breaks the ionic lock, leading to activation of the Gs protein coupled to β₂AR. Upon activation, the Gαs protein undergoes a remarkable opening of the GDP binding pocket, making the GDP available for exchange or release. Concomitantly, Gαs-α5 undergoes a remarkable expansion in the β₂AR cytoplasmic region after the ionic lock is broken, inducing TM6 to displace outward by ~5 Å from TM3.

Molecular Biophysics of Class A G Protein Coupled Receptors-Lipids Interactome at a Glance-Highlights from the A2A Adenosine Receptor

G protein-coupled receptors (GPCRs) are embedded in phospholipid membrane bilayers with cholesterol representing 34% of the total lipid content in mammalian plasma membranes. Membrane lipids interact with GPCRs structures and modulate their function and drug-stimulated signaling through conformational selection. It has been shown that anionic phospholipids form strong interactions between positively charged residues in the G protein and the TM5-TM6-TM 7 cytoplasmic interface of class A GPCRs stabilizing the signaling GPCR-G complex. Cholesterol with a high content in plasma membranes can be identified in more specific sites in the transmembrane region of GPCRs, such as the Cholesterol Consensus Motif (CCM) and Cholesterol Recognition Amino Acid Consensus (CRAC) motifs and other receptor dependent and receptor state dependent sites. Experimental biophysical methods, atomistic (AA) MD simulations and coarse-grained (CG) molecular dynamics simulations have been applied to investigate these interactions. We emphasized here the impact of phosphatidyl inositol-4,5-bisphosphate (PtdIns(4,5)P₂ or PIP₂), a minor phospholipid component and of cholesterol on the function-related conformational equilibria of the human A_2A adenosine receptor (A_2AR), a representative receptor in class A GPCR. Several GPCRs of class A interacted with PIP₂ and cholesterol and in many cases the mechanism of the modulation of their function remains unknown. This review provides a helpful comprehensive overview for biophysics that enter the field of GPCRs-lipid systems.

Molecular mechanism of fatty acid activation of FFAR1

FFAR1 is a G-protein-coupled receptor (GPCR) that responds to circulating free fatty acids to enhance glucose-stimulated insulin secretion and release of incretin hormones. Due to the glucose-lowering effect of FFAR1 activation, potent agonists for this receptor have been developed for the treatment of diabetes. Previous structural and biochemical studies of FFAR1 showed multiple sites of ligand binding to the inactive state but left the mechanism of fatty acid interaction and receptor activation unknown. We used cryo-electron microscopy to elucidate structures of activated FFAR1 bound to a Gq mimetic, which were induced either by the endogenous FFA ligand docosahexaenoic acid or γ-linolenic acid and the agonist drug TAK-875. Our data identify the orthosteric pocket for fatty acids and show how both endogenous hormones and synthetic agonists induce changes in helical packing along the outside of the receptor that propagate to exposure of the G-protein-coupling site. These structures show how FFAR1 functions without the highly conserved "DRY" and "NPXXY" motifs of class A GPCRs and also illustrate how the orthosteric site of a receptor can be bypassed by membrane-embedded drugs to confer full activation of G protein signaling.

Pharmacological hallmarks of allostery at the M4 muscarinic receptor elucidated through structure and dynamics

Allosteric modulation of G protein-coupled receptors (GPCRs) is a major paradigm in drug discovery. Despite decades of research, a molecular-level understanding of the general principles that govern the myriad pharmacological effects exerted by GPCR allosteric modulators remains limited. The M₄ muscarinic acetylcholine receptor (M₄ mAChR) is a validated and clinically relevant allosteric drug target for several major psychiatric and cognitive disorders. In this study, we rigorously quantified the affinity, efficacy, and magnitude of modulation of two different positive allosteric modulators, LY2033298 (LY298) and VU0467154 (VU154), combined with the endogenous agonist acetylcholine (ACh) or the high-affinity agonist iperoxo (Ipx), at the human M₄ mAChR. By determining the cryo-electron microscopy structures of the M₄ mAChR, bound to a cognate G_i1 protein and in complex with ACh, Ipx, LY298-Ipx, and VU154-Ipx, and applying molecular dynamics simulations, we determine key molecular mechanisms underlying allosteric pharmacology. In addition to delineating the contribution of spatially distinct binding sites on observed pharmacology, our findings also revealed a vital role for orthosteric and allosteric ligand-receptor-transducer complex stability, mediated by conformational dynamics between these sites, in the ultimate determination of affinity, efficacy, cooperativity, probe dependence, and species variability. There results provide a holistic framework for further GPCR mechanistic studies and can aid in the discovery and design of future allosteric drugs.

Intermediate-state-trapped mutants pinpoint G protein-coupled receptor conformational allostery

Understanding the roles of intermediate states in signaling is pivotal to unraveling the activation processes of G protein-coupled receptors (GPCRs). However, the field is still struggling to define these conformational states with sufficient resolution to study their individual functions. Here, we demonstrate the feasibility of enriching the populations of discrete states via conformation-biased mutants. These mutants adopt distinct distributions among five states that lie along the activation pathway of adenosine A2A receptor (A2AR), a class A GPCR. Our study reveals a structurally conserved cation-π lock between transmembrane helix VI (TM6) and Helix8 that regulates cytoplasmic cavity opening as a "gatekeeper" for G protein penetration. A GPCR activation process based on the well-discerned conformational states is thus proposed, allosterically micro-modulated by the cation-π lock and a previously well-defined ionic interaction between TM3 and TM6. Intermediate-state-trapped mutants will also provide useful information in relation to receptor-G protein signal transduction.

Gait Abnormalities and Aberrant D2 Receptor Expression and Signaling in Mice Carrying the Human Pathogenic Mutation DRD2I212F

A dopamine D2 receptor mutation was recently identified in a family with a novel hyperkinetic movement disorder. That allelic variant D2-I212F is a constitutively active and G protein-biased receptor. We now describe mice engineered using CRISPR-Cas9-mediated gene editing technology to carry the D2-I212F variant. Drd2I212F mice exhibited gait abnormalities resembling those in other mouse models of chorea and/or dystonia and had striatal D2 receptor expression that was decreased approximately 30% per Drd2I212F allele. Electrically evoked inhibitory postsynaptic conductances in midbrain dopamine neurons and striatum from Drd2I212F mice, caused by G protein activation of potassium channels, exhibited slow kinetics (e.g., approximately four- to sixfold slower decay) compared with Drd2 +/+ mice. Current decay initiated by photolytic release of the D2 antagonist sulpiride from CyHQ-sulpiride was also ∼fourfold slower in midbrain slices from Drd2I212F mice than Drd2 +/+ mice. Furthermore, in contrast to Drd2 +/+ mice, in which dopamine is several-fold more potent at neurons in the nucleus accumbens than in the dorsal striatum, reflecting activation of Gα o versus Gα i, dopamine had similar potencies in those two brain regions of Drd2I212F mice. Repeated cocaine treatment, which decreases dopamine potency in the nucleus accumbens of Drd2 +/+ mice, had no effect on dopamine potency in Drd2 I212F mice. The results demonstrate the pathogenicity of the D2-I212F mutation and the utility of this mouse model for investigating the role of pathogenic DRD2 variants in early-onset hyperkinetic movement disorders. SIGNIFICANCE STATEMENT: The first dopamine receptor mutation to cause a movement disorder, D2-I212F, was recently identified. The mutation makes receptor activation of G protein-mediated signaling more efficient. To confirm the pathogenesis of D2-I212F, this study reports that mice carrying this mutation have gait abnormalities consistent with the clinical phenotype. The mutation also profoundly alters D2 receptor expression and function in vivo. This mouse model will be useful for further characterization of the mutant receptor and for evaluation of potential therapeutic drugs.

Mandibulofacial dysostosis with alopecia results from ETAR gain-of-function mutations via allosteric effects on ligand binding

Mutations of G protein-coupled receptors (GPCRs) cause various human diseases, but the mechanistic details are limited. Here, we establish p.E303K in the gene encoding the endothelin receptor type A (ETAR/EDNRA) as a recurrent mutation causing mandibulofacial dysostosis with alopecia (MFDA), with craniofacial changes similar to those caused by p.Y129F. Mouse models carrying either of these missense mutations exhibited a partial maxillary-to-mandibular transformation, which was rescued by deleting the ligand endothelin 3 (ET3/EDN3). Pharmacological experiments confirmed the causative ETAR mutations as gain of function, dependent on ET3. To elucidate how an amino acid substitution far from the ligand binding site can increase ligand affinity, we used molecular dynamics (MD) simulations. E303 is located at the intracellular end of transmembrane domain 6, and its replacement by a lysine increased flexibility of this portion of the helix, thus favoring G protein binding and leading to G protein-mediated enhancement of agonist affinity. The Y129F mutation located under the ligand binding pocket reduced the sodium-water network, thereby affecting the extracellular portion of helices in favor of ET3 binding. These findings provide insight into the pathogenesis of MFDA and into allosteric mechanisms regulating GPCR function, which may provide the basis for drug design targeting GPCRs.

Structural elucidation of ETHR-A and ETHR-B from Plutella xylostella and insight into non-conservative mutations in transmembrane helix-6

The development of Diamondback moth (DBM) depends on the ecdysis triggering hormone receptor (ETHR); a neuronal membrane G-protein coupled receptor (GPCR) connected to the metamorphosis cascade. Lepidopteran insect DBM is an infamous pest of cruciferous plants. This study examined the full-length coding sequences (CDS) of PxETHR-A and PxETHR-B from the DBM genome. The three-dimensional (3 D) models of both receptors were generated in an inactive state. The behaviour and stability of receptors were examined using molecular dynamics simulations in a lipid membrane system for 300 ns and established a GPCR family-based view. Secondary interactions within receptors were studied to know more about factors contributing to their stability. Multiple sequence alignment revealed conserved features of insect ETHRs those compared to the GPCR family proteins. These features were helpful during the evaluation of the molecular models of both receptors. Side-chain orientation of conserved residues, non-conserved and conserved hydrogen-bond networks (HBN) and hydrophobic clusters were examined in the structures of both receptors. The non-conserved residues L6.35, T6.39, C/S6.43, and L6.48, are present in a conserved position on the transmembrane helix-6 (TM6) of ETHRs. In TM6, PxETHR-A and PxETHR-B differ at positions C/S6.43 and Y/F6.51, both being part of the HBN.Communicated by Ramaswamy H. Sarma.

Superconserved receptors expressed in the brain: Expression, function, motifs and evolution of an orphan receptor family

GPR27, GPR85 and GPR173 constitute a small family of G protein-coupled receptors (GPCR) that share the distinctive characteristics of being highly conserved throughout vertebrate evolution and predominantly expressed in the brain. Accordingly, they have been coined as "Superconserved Receptors Expressed in the Brain" (SREB), although their expression profile is more complex than what was originally thought. SREBs have no known validated endogenous ligands and are thus labeled as "orphan" receptors. The investigation of this particular category of uncharacterized receptors holds great promise both in terms of physiology and drug development. In the largest GPCR family, the Rhodopsin-like or Class A, around 100 receptors are considered orphans. Because GPCRs are the most successful source of drug targets, the discovery of a novel function or ligand most likely will lead to significant breakthroughs for the discovery of innovative therapies. The high level of conservation is one of the characteristic features of the SREBs. We propose herein a detailed analysis of the putative evolutionary origin of this family. We highlight the properties that distinguish SREBs from other rhodopsin-like GPCRs. We present the current evidence for these receptors downstream signaling pathways and functions. We discuss the pharmacological challenge for the identification of natural or synthetic ligands of orphan receptors like SREBs. The different SREB-related scientific questions are presented with a highlight on what should be addressed in the near future, including the confirmation of published evidence and their validation as drug targets. In particular, we discuss in which pathological conditions these receptors may be of great relevance to solve unmet medical needs.

Dynamic and thermodynamic impact of L94A, W100A, and W100L mutations on the D2 dopamine receptor bound to risperidone

DRD2 is an important receptor in the mediation of antipsychotic drugs but also in Parkinson medication, hyperprolactinemia, nausea and vomiting. Recently, crystallographic studies of the DRD2-risperidone complex have provided important information about risperidone recognition in wild-type and different stabilizing DRD2-risperidone residues. Using the crystallographic structure of the DRD2-risperidone complex as a starting point, we undertook molecular dynamics (MD) simulations to investigate the structural and thermodynamic basis of molecular recognition by risperidone at the ligand-binding sites of wild-type and mutant DRD2. A solvated phospholipid bilayer was used to construct DRD2-risperidone complexes, which were then subjected to several microsecond (μs) MD simulations in order to obtain realistic receptor-ligand conformations under the equilibrated simulation time. Risperidone had a higher affinity for wild-type and L94A mutant DRD2 than the W100L and W100A mutants, according to binding free energy calculations using the Molecular Mechanics Generalized-Born Surface Area (MMGBSA) method, explaining the experimental differences in ligand residence times. Principal component (PC) analysis revealed important conformational mobility upon molecular recognition of risperidone for the L94A mutant compared to the wild type, indicating an unfavorable entropic component that may contribute to improving risperidone affinity in the L94A DRD2 mutant.

Class A and C GPCR Dimers in Neurodegenerative Diseases

Neurodegenerative diseases affect over 30 million people worldwide with an ascending trend. Most individuals suffering from these irreversible brain damages belong to the elderly population, with onset between 50 and 60 years. Although the pathophysiology of such diseases is partially known, it remains unclear upon which point a disease turns degenerative. Moreover, current therapeutics can treat some of the symptoms but often have severe side effects and become less effective in long-term treatment. For many neurodegenerative diseases, the involvement of G proteincoupled receptors (GPCRs), which are key players of neuronal transmission and plasticity, has become clearer and holds great promise in elucidating their biological mechanism. With this review, we introduce and summarize class A and class C GPCRs, known to form heterodimers or oligomers to increase their signalling repertoire. Additionally, the examples discussed here were shown to display relevant alterations in brain signalling and had already been associated with the pathophysiology of certain neurodegenerative diseases. Lastly, we classified the heterodimers into two categories of crosstalk, positive or negative, for which there is known evidence.

The mechanism for ligand activation of the GPCR-G protein complex

G protein–coupled receptors (GPCRs) activate cellular responses ranging from odorants to neurotransmitters. Binding an agonist leads to activation of a heterotrimeric G protein (GP) that stimulates external signaling. Unfortunately, the mechanism remains unknown. We show for 15 class A GPCRs, including opioids, adrenergics, adenosines, chemokines, muscarinics, cannabinoids, serotonins, and dopamines, that interaction of an inactive GP, including Gs, Gi, Go, G11, and Gq, to the inactive GPCR, containing the intracellular ionic lock between transmembrane (TM) helices 3 and 6, evolves exothermically to form a precoupled GPCR-GP complex with an opened TM3-TM6 and the GP-α5 helix partially inserted into the GPCR but not activated. We show that binding of agonist to this precoupled GPCR-GP complex causes the Gα protein to open into its active form, with the guanosine diphosphate exposed for signaling. This GP-first paradigm provides a strategy for developing selective agonists for GPCRs since it is the pharmacophore for the precoupled GPCR-GP complex that should be used to design drugs.

True-atomic-resolution insights into the structure and functional role of linear chains and low-barrier hydrogen bonds in proteins

Hydrogen bonds are fundamental to the structure and function of biological macromolecules and have been explored in detail. The chains of hydrogen bonds (CHBs) and low-barrier hydrogen bonds (LBHBs) were proposed to play essential roles in enzyme catalysis and proton transport. However, high-resolution structural data from CHBs and LBHBs is limited. The challenge is that their 'visualization' requires ultrahigh-resolution structures of the ground and functionally important intermediate states to identify proton translocation events and perform their structural assignment. Our true-atomic-resolution structures of the light-driven proton pump bacteriorhodopsin, a model in studies of proton transport, show that CHBs and LBHBs not only serve as proton pathways, but also are indispensable for long-range communications, signaling and proton storage in proteins. The complete picture of CHBs and LBHBs discloses their multifunctional roles in providing protein functions and presents a consistent picture of proton transport and storage resolving long-standing debates and controversies.

The Many Faces of G Protein-Coupled Receptor 143, an Atypical Intracellular Receptor

GPCRs transform extracellular stimuli into a physiological response by activating an intracellular signaling cascade initiated via binding to G proteins. Orphan G protein-coupled receptors (GPCRs) hold the potential to pave the way for development of new, innovative therapeutic strategies. In this review we will introduce G protein-coupled receptor 143 (GPR143), an enigmatic receptor in terms of classification within the GPCR superfamily and localization. GPR143 has not been assigned to any of the GPCR families due to the lack of common structural motifs. Hence we will describe the most important motifs of classes A and B and compare them to the protein sequence of GPR143. While a precise function for the receptor has yet to be determined, the protein is expressed abundantly in pigment producing cells. Many GPR143 mutations cause X-linked Ocular Albinism Type 1 (OA1, Nettleship-Falls OA), which results in hypopigmentation of the eyes and loss of visual acuity due to disrupted visual system development and function. In pigment cells of the skin, loss of functional GPR143 results in abnormally large melanosomes (organelles in which pigment is produced). Studies have shown that the receptor is localized internally, including at the melanosomal membrane, where it may function to regulate melanosome size and/or facilitate protein trafficking to the melanosome through the endolysosomal system. Numerous additional roles have been proposed for GPR143 in determining cancer predisposition, regulation of blood pressure, development of macular degeneration and signaling in the brain, which we will briefly describe as well as potential ligands that have been identified. Furthermore, GPR143 is a promiscuous receptor that has been shown to interact with multiple other melanosomal proteins and GPCRs, which strongly suggests that this orphan receptor is likely involved in many different physiological actions.

Computational Methods for Understanding the Selectivity and Signal Transduction Mechanism of Aminomethyl Tetrahydronaphthalene to Opioid Receptors

Opioid receptors are members of the group of G protein-couple receptors, which have been proven to be effective targets for treating severe pain. The interactions between the opioid receptors and corresponding ligands and the receptor’s activation by different agonists have been among the most important fields in opioid research. In this study, with compound M1, an active metabolite of tramadol, as the clue compound, several aminomethyl tetrahydronaphthalenes were designed, synthesized and assayed upon opioid receptors. With the resultant compounds FW-AII-OH-1 (Ki = 141.2 nM for the κ opioid receptor), FW-AII-OH-2 (Ki = 4.64 nM for the δ opioid receptor), FW-DI-OH-2 (Ki = 8.65 nM for the δ opioid receptor) and FW-DIII-OH-2 (Ki = 228.45 nM for the δ opioid receptor) as probe molecules, the structural determinants responsible for the subtype selectivity and activation mechanisms were further investigated by molecular modeling and molecular dynamics simulations. It was shown that Y^7.43 was a key residue in determining the selectivity of the three opioid receptors, and W^6.58 was essential for the selectivity of the δ opioid receptor. A detailed stepwise discovered agonist-induced signal transduction mechanism of three opioid receptors by aminomethyl tetrahydronaphthalene compounds was proposed: the 3–7 lock between TM3 and TM7, the DRG lock between TM3 and TM6 and rearrangement of I^3.40, P^5.50 and F^6.44, which resulted in the cooperative movement in 7 TMs. Then, the structural relaxation left room for the binding of the G protein at the intracellular site, and finally the opioid receptors were activated.

Ligand-Dependent Modulation of the Dynamics of Intracellular Loops Dictates Functional Selectivity of 5-HT2AR

The serotonin 2A receptor (5-HT_2AR) subtype of the G protein-coupled receptor (GPCR) family is involved in a plethora of neuromodulatory functions (e.g., neurogenesis, sleep, and cognitive processes). 5-HT_2AR is the target of pharmacologically distinct classes of ligands, binding of which either activate or inactivate the receptor. Although high-resolution structures of 5-HT_2AR as well as several other 5-HT GPCRs provided snapshots of both active and inactive conformational states, these structures, representing a truncated form of the receptor, cannot fully explain the mechanism of conformational transitions during their function. Importantly, biochemical studies have suggested the importance of intracellular loops in receptor functions. In our previous study, a model of the ligand-free form of 5-HT_2AR with the third intracellular loop (ICL3) has been meticulously built. Here, we have investigated the functional regulation of 5-HT_2AR with intact intracellular loops in ligand-free and five distinct ligand-bound configurations using unbiased atomistic molecular dynamics (MD) simulations. The selected ligands belong to either of the full, partial, or inverse agonist classes, which exert distinct pharmacological responses. We have observed significant structural, dynamic, and thermodynamic differences within ligand-bound complexes. Our results revealed, for the first time, that either activation or inactivation of the receptor upon specific ligand binding is primarily achieved through conformational transitions of its second and third intracellular loops (ICL2 and ICL3). A remarkable allosteric cross-talk between the ligand-binding site and the distal intracellular parts of the receptor, where binding of a specific ligand thermodynamically controls (either stabilizes or destabilizes) the intracellular region, consisting of crucial dynamic elements ICL2 and ICL3, and differential conformational transitions of these loops determine ligand-dependent functional selectivity.

Molecular insights into the allosteric coupling mechanism between an agonist and two different transducers for μ-opioid receptors

G protein-coupled receptors (GPCRs) as the most important class of pharmacological targets regulate G-protein and β-arrestin-mediated signaling through allosteric interplay, which are responsible for different biochemical and physiological actions like therapeutic efficacy and side effects. However, the allosteric mechanism underlying preferentially recruiting one transducer versus the other has been poorly understood, limiting drug design. Motivated by this issue, we utilize accelerated molecular dynamics simulation coupled with potential of mean force (PMF), molecular mechanics Poisson Boltzmann surface area (MM/PBSA) and protein structure network (PSN) to study two ternary complex systems of a representative class A GPCR (μ-opioid receptor (μOR)) bound by an agonist and one specific transducer (G-protein and β-arrestin). The results show that no significant difference exists in the whole structure of μOR between two transducer couplings, but displays transducer-dependent changes in the intracellular binding region of μOR, where the β-arrestin coupling results in a narrower crevice with TM7 inward movement compared with the G-protein. In addition, both the G-protein and β-arrestin coupling can increase the binding affinity of the agonist to the receptor. However, the interactions between the agonist and μOR also exhibit transducer-specific changes, in particular for the interaction with ECL2 that plays an important role in recruiting β-arrestin. The allosteric network analysis further indicates that Y148^3.33, F152^3.37, F156^3.41, N191^4.49, T160^3.45, Y106^2.42, W293^6.48, F289^6.44, I248^5.54 and Y252^5.58 play important roles in equally activating G-protein and β-arrestin. In contrast, M161^3.46 and R165^3.50 devote important contributions to preferentially recruit G-protein while D164^3.49 and R179^ICL2 are revealed to be important for selectively activating β-arrestin. The observations provide useful information for understanding the biased activation mechanism.

Dynamics and Mechanistic Underpinnings to Pharmacology of Class A GPCRs - An NMR Perspective

One-third of current pharmaceuticals target G protein-coupled receptors (GPCRs), the largest receptor superfamily in humans and mediators of diverse physiological processes. This review summarizes the recent progress in GPCR structural dynamics, focusing on class A receptors and insights derived from nuclear magnetic resonance (NMR) and other spectroscopic techniques. We describe the structural aspects of GPCR activation and the various pharmacological models that capture aspects of receptor signaling behaviour. Spectroscopic studies revealed that receptors and their signaling complexes are dynamic allosteric systems that sample multiple functional states under basal conditions. The distribution of states within the conformational ensemble and the kinetics of transitions between states are regulated through the binding of ligands, allosteric modulators, and the membrane environment. This ensemble view of GPCRs provides a mechanistic framework for understanding many of the pharmacological phenomena associated with receptor signaling, such as basal activity, efficacy, and functional bias.

Allosteric modulation of dopamine D2L receptor in complex with Gi1 and Gi2 proteins: the effect of subtle structural and stereochemical ligand modifications

Background

Allosteric modulation of G protein-coupled receptors (GPCRs) is nowadays one of the hot topics in drug discovery. In particular, allosteric modulators of D₂ receptor have been proposed as potential modern therapeutics to treat schizophrenia and Parkinson’s disease.

Methods

To address some subtle structural and stereochemical aspects of allosteric modulation of D₂ receptor, we performed extensive in silico studies of both enantiomers of two compounds (compound 1 and compound 2), and one of them (compound 2) was synthesized as a racemate in-house and studied in vitro.

Results

Our molecular dynamics simulations confirmed literature reports that the R enantiomer of compound 1 is a positive allosteric modulator of the D_2L receptor, while its S enantiomer is a negative allosteric modulator. Moreover, based on the principal component analysis (PCA), we hypothesized that both enantiomers of compound 2 behave as silent allosteric modulators, in line with our in vitro studies. PCA calculations suggest that the most pronounced modulator-induced receptor rearrangements occur at the transmembrane helix 7 (TM7). In particular, TM7 bending at the conserved P7.50 and G7.42 was observed. The latter resides next to the Y7.43, which is a significant part of the orthosteric binding site. Moreover, the W7.40 conformation seems to be affected by the presence of the positive allosteric modulator.

Conclusions

Our work reveals that allosteric modulation of the D_2L receptor can be affected by subtle ligand modifications. A change in configuration of a chiral carbon and/or minor structural modulator modifications are solely responsible for the functional outcome of the allosteric modulator.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s43440-021-00352-x.

Novel interaction between neurotrophic factor-α1/carboxypeptidase E and serotonin receptor, 5-HTR1E, protects human neurons against oxidative/neuroexcitotoxic stress via β-arrestin/ERK signaling

Protecting neurons from death during oxidative and neuroexcitotoxic stress is key for preventing cognitive dysfunction. We uncovered a novel neuroprotective mechanism involving interaction between neurotrophic factor-α1 (NF-α1/carboxypeptidase E, CPE) and human 5-HTR1E, a G protein-coupled serotonin receptor with no previously known neurological function. Co-immunoprecipitation and pull-down assays confirmed interaction between NFα1/CPE and 5-HTR1E and 125I NF-α1/CPE-binding studies demonstrated saturable, high-affinity binding to 5-HTR1E in stably transfected HEK293 cells (Kd = 13.82 nM). Treatment of 5-HTR1E stable cells with NF-α1/CPE increased pERK 1/2 and pCREB levels which prevented a decrease in pro-survival protein, BCL2, during H2O2-induced oxidative stress. Cell survival assay in β-arrestin Knockout HEK293 cells showed that the NF-α1/CPE-5-HTR1E-mediated protection against oxidative stress was β-arrestin-dependent. Molecular dynamics studies revealed that NF-α1/CPE interacts with 5-HTR1E via 3 salt bridges, stabilized by several hydrogen bonds, independent of the serotonin pocket. Furthermore, after phosphorylating the C-terminal tail and intracellular loop 3 (ICL3) of NF-α1/CPE-5-HTR1E, it recruited β-arrestin1 by forming numerous salt bridges and hydrogen bonds to ICL2 and ICL3, leading to activation of β-arrestin1. Immunofluorescence studies showed 5-HTR1E and NF-α1/CPE are highly expressed and co-localized on cell surface of human hippocampal neurons. Importantly, knock-down of 5-HTR1E in human primary neurons diminished the NF-α1/CPE-mediated protection of these neurons against oxidative stress and glutamate neurotoxicity-induced cell death. Thus, NF-α1/CPE uniquely interacts with serotonin receptor 5-HTR1E to activate the β-arrestin/ERK/CREB/BCL2 pathway to mediate stress-induced neuroprotection.

Naturally occurring mutations in G protein-coupled receptors associated with obesity and type 2 diabetes mellitus

G protein-coupled receptors (GPCRs) are the largest family of membrane receptors involved in the regulation of almost all known physiological processes. Dysfunctions of GPCR-mediated signaling have been shown to cause various diseases. The prevalence of obesity and type 2 diabetes mellitus (T2DM), two strongly associated disorders, is increasing worldwide, with tremendous economical and health burden. New safer and more efficacious drugs are required for successful weight reduction and T2DM treatment. Multiple GPCRs are involved in the regulation of energy and glucose homeostasis. Mutations in these GPCRs contribute to the development and progression of obesity and T2DM. Therefore, these receptors can be therapeutic targets for obesity and T2DM. Indeed some of these receptors, such as melanocortin-4 receptor and glucagon-like peptide 1 receptor, have provided important new drugs for treating obesity and T2DM. This review will focus on the naturally occurring mutations of several GPCRs associated with obesity and T2DM, especially incorporating recent large genomic data and insights from structure-function studies, providing leads for future investigations.

Lipopeptide Pepducins as Therapeutic Agents

Pepducins are lipidated peptides that target the intracellular loops of G protein-coupled receptors (GPCRs) in order to modulate transmembrane signaling to internally located effectors. With a wide array of potential activities ranging from partial, biased, or full agonism to antagonism, pepducins represent a versatile class of compounds that can be used to potentially treat diverse human diseases or be employed as novel tools to probe complex mechanisms of receptor activation and signaling in cells and in animals. Here, we describe a number of different pepducins including an advanced compound, PZ-128, that has successfully progressed through phase 2 clinical trials in cardiac patients demonstrating safety and efficacy in suppressing myonecrosis and arterial thrombosis.

Predicted Structure of Fully Activated Tas1R3/1R3' Homodimer Bound to G Protein and Natural Sugars: Structural Insights into G Protein Activation by a Class C Sweet Taste Homodimer with Natural Sugars

The Tas1R3 G protein-coupled receptor constitutes the main component of sweet taste sensory response in humans via forming a heterodimer with Tas1R2 or a homodimer with Tas1R3. The Tas1R3/1R3' homodimer serves as a low-affinity sweet taste receptor, stimulating gustducin G protein (G_Gust) signaling in the presence of a high concentration of natural sugars. This provides an additional means to detect the taste of natural sugars, thereby differentiating the flavors between natural sugars and artificial sweeteners. We report here the predicted 3D structure of active state Tas1R3/1R3' homodimer complexed with heterotrimeric G_Gust and sucrose. We discovered that the G_Gust makes ionic anchors to intracellular loops 1 and 2 of Tas1R3 while the Gα-α5 helix engages the cytoplasmic region extensively through salt bridge and hydrophobic interactions. We show that in the activation of this complex the Venus flytrap domains of the homodimer undergo a remarkable twist up to ∼100° rotation around the vertical axis to adopt a closed-closed conformation while the intracellular region relaxes to an open-open conformation. We find that binding of sucrose to the homodimer stabilizes a preactivated conformation with a largely open intracellular region that recruits and activates the G_Gust. Upon activation, the Gα subunit spontaneously opens up the nucleotide-binding site, making nucleotide exchange facile for signaling. This activation of G_Gust promotes the interdomain twist of the Venus flytrap domains. These structures and transformations could potentially be a basis for the design of new sweeteners with higher activity and less unpleasant flavors.

Universal Properties and Specificities of the β2-Adrenergic Receptor-Gs Protein Complex Activation Mechanism Revealed by All-Atom Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are transmembrane proteins of high pharmacological relevance. It has been proposed that their activity is linked to structurally distinct, dynamically interconverting functional states and the process of activation relies on an interconnecting network of conformational switches in the transmembrane domain. However, it is yet to be uncovered how ligands with different extents of functional effect exert their actions. According to our recent hypothesis, based on indirect observations and the literature data, the transmission of the external stimulus to the intracellular surface is accompanied by the shift of macroscopic polarization in the transmembrane domain, furnished by concerted movements of highly conserved polar motifs and the rearrangement of polar species. In this follow-up study, we have examined the β₂-adrenergic receptor (β₂AR) to see if our hypothesis drawn from an extensive study of the μ-opioid receptor (MOP) is fundamental and directly transferable to other class A GPCRs. We have found that there are some general similarities between the two receptors, in agreement with previous studies, and there are some receptor-specific differences that could be associated with different signaling pathways.

5-HT2C agonists and antagonists block different components of behavioral responses to potential, distal, and proximal threat in zebrafish

Serotonin (5-HT) receptors have been implicated in responses to aversive stimuli in mammals and fish, but its precise role is still unknown. Moreover, since at least seven families of 5-HT receptors exist in vertebrates, the role of specific receptors is still debated. Aversive stimuli can be classified as indicators of proximal, distal, or potential threat, initiating responses that are appropriate for each of these threat levels. Responses to potential threat usually involve cautious exploration and increased alertness, while responses to distal and proximal threat involve a fight-flight-freeze reaction. We exposed adult zebrafish to a conspecific alarm substance (CAS) and observed behavior during (distal threat) and after (potential threat) exposure, and treated with the 5-HT_2C receptor agonists MK-212 or WAY-161503 or with the antagonist RS-102221. The agonists blocked CAS-elicited defensive behavior (distal threat), but not post-exposure increases in defensive behavior (potential threat), suggesting inhibition of responses to distal threat. MK-212 blocked changes in freezing elicited by acute restraint stress, a model of proximal threat, while RS-102221 blocked changes in geotaxis elicited this stressor. We also found that RS-102221, a 5-HT_2C receptor antagonist, produced small effect on behavior during and after exposure to CAS. Preprint: https://www.biorxiv.org/content/10.1101/2020.10.04.324202; Data and scripts: https://github.com/lanec-unifesspa/5-HT-CAS/tree/master/data/5HT2C.

CB1 Cannabinoid Receptor Signaling and Biased Signaling

The CB1 cannabinoid receptor is a G-protein coupled receptor highly expressed throughout the central nervous system that is a promising target for the treatment of various disorders, including anxiety, pain, and neurodegeneration. Despite the wide therapeutic potential of CB1, the development of drug candidates is hindered by adverse effects, rapid tolerance development, and abuse potential. Ligands that produce biased signaling—the preferential activation of a signaling transducer in detriment of another—have been proposed as a strategy to dissociate therapeutic and adverse effects for a variety of G-protein coupled receptors. However, biased signaling at the CB1 receptor is poorly understood due to a lack of strongly biased agonists. Here, we review studies that have investigated the biased signaling profile of classical cannabinoid agonists and allosteric ligands, searching for a potential therapeutic advantage of CB1 biased signaling in different pathological states. Agonist and antagonist bound structures of CB1 and proposed mechanisms of action of biased allosteric modulators are used to discuss a putative molecular mechanism for CB1 receptor activation and biased signaling. Current studies suggest that allosteric binding sites on CB1 can be explored to yield biased ligands that favor or hinder conformational changes important for biased signaling.

The G protein-first activation mechanism of opioid receptors by Gi protein and agonists

We report the G protein-first mechanism for activation of G protein-coupled receptors (GPCR) for the three closest subtypes of the opioid receptors (OR), μOR, κOR and δOR. We find that they couple to the inactive Gi protein-bound guanosine diphosphate (GDP) prior to agonist binding. The inactive Gi protein forms anchors to the intracellular loops of the inactive apo-μOR, apo-κOR and apo-δOR, inducing opening of the cytoplasmic region to form a pre-activated state that holds Gi protein in place until agonist binds. Then, agonist binds to μOR, κOR and δOR already complexed with Gi protein, to trigger the Gαi to open up the tightly coupled GDP binding site, making GDP accessible for GTP exchange, an essential step for Gi signalling. We show that the agonist alone cannot open the intracellular region of μOR and κOR, requiring Gi protein to open the cytoplasmic region by itself. We consider that this G protein-first mechanism may apply to activation of other Class A GPCRs. However, for δOR, agonist binding can open up the intracellular region to encourage Gi protein recruitment. Thus, activation of Gi protein mediated by δOR favourably may proceed with either ligand-first or G protein-first activation mechanisms.

Implications of the simple chemical structure of the odorant molecules interacting with the olfactory receptor 1A1

G protein‒coupled receptors (GPCRs), including olfactory receptors, account for the largest group of genes in the human genome and occupy a very important position in signaling systems. Although olfactory receptors, which belong to the broader category of GPCRs, play an important role in monitoring the organism’s surroundings, their actual three-dimensional structure has not yet been determined. Therefore, the specific details of the molecular interactions between the receptor and the ligand remain unclear. In this report, the interactions between human olfactory receptor 1A1 and its odorant molecules were simulated using computational methods, and we explored how the chemically simple odorant molecules activate the olfactory receptor.

Signaling-Biased and Constitutively Active Dopamine D2 Receptor Variant

A dopamine D2 receptor mutation was recently identified in a family with a novel hyperkinetic movement disorder. Compared to the wild type D2 receptor, the novel allelic variant D2-I²¹²F activates a Gα_i1β₁γ₂ heterotrimer with higher potency and modestly enhanced basal activity in human embryonic kidney (HEK) 293 cells and has decreased capacity to recruit arrestin3. We now report that omitting overexpressed G protein-coupled receptor kinase-2 (GRK2) decreased the potency and efficacy of quinpirole for arrestin recruitment. The relative efficacy of quinpirole for arrestin recruitment to D2-I²¹²F compared to D2-WT was considerably lower without overexpressed GRK2 than with added GRK2. D2-I²¹²F exhibited higher basal activation of Gα_oA than Gα_i1 but little or no increase in the potency of quinpirole relative to D2-WT. Other signs of D2-I²¹²F constitutive activity for G protein-mediated signaling, in addition to basal activation of Gα_i/o, were enhanced basal inhibition of forskolin-stimulated cyclic AMP accumulation that was reversed by the inverse agonists sulpiride and spiperone and a ∼4-fold increase in the apparent affinity of D2-I²¹²F for quinpirole, determined from competition binding assays. In mouse midbrain slices, inhibition of tonic current by the inverse agonist sulpiride in dopamine neurons expressing D2-I²¹²F was consistent with our hypothesis of enhanced constitutive activity and sensitivity to dopamine relative to D2-WT. Molecular dynamics simulations with D2 receptor models suggested that an ionic lock between the cytoplasmic ends of the third and sixth α-helices that constrains many G protein-coupled receptors in an inactive conformation spontaneously breaks in D2-I²¹²F. Overall, these results confirm that D2-I²¹²F is a constitutively active and signaling-biased D2 receptor mutant and also suggest that the effect of the likely pathogenic variant in a given brain region will depend on the nature of G protein and GRK expression.

Analysis of Missense Variants in the Human Histamine Receptor Family Reveals Increased Constitutive Activity of E4106.30×30K Variant in the Histamine H1 Receptor

The Exome Aggregation Consortium has collected the protein-encoding DNA sequences of almost 61,000 unrelated humans. Analysis of this dataset for G protein-coupled receptor (GPCR) proteins (available at GPCRdb) revealed a total of 463 naturally occurring genetic missense variations in the histamine receptor family. In this research, we have analyzed the distribution of these missense variations in the four histamine receptor subtypes concerning structural segments and sites important for GPCR function. Four missense variants R127^3.52×52H, R139^34.57×57H, R409^6.29×29H, and E410^6.30×30K, were selected for the histamine H₁ receptor (H₁R) that were hypothesized to affect receptor activity by interfering with the interaction pattern of the highly conserved D(E)RY motif, the so-called ionic lock. The E410^6.30×30K missense variant displays higher constitutive activity in G protein signaling as compared to wild-type H₁R, whereas the opposite was observed for R127^3.52×52H, R139^34.57×57H, and R409^6.29×29H. The E410^6.30×30K missense variant displays a higher affinity for the endogenous agonist histamine than wild-type H₁R, whereas antagonist affinity was not affected. These data support the hypothesis that the E410^6.30×30K mutation shifts the equilibrium towards active conformations. The study of these selected missense variants gives additional insight into the structural basis of H₁R activation and, moreover, highlights that missense variants can result in pharmacologically different behavior as compared to wild-type receptors and should consequently be considered in the drug discovery process.

Delineating the conformational landscape of the adenosine A2A receptor during G protein coupling

G-protein-coupled receptors (GPCRs) represent a ubiquitous membrane protein family and are important drug targets. Their diverse signaling pathways are driven by complex pharmacology arising from a conformational ensemble rarely captured by structural methods. Here, fluorine nuclear magnetic resonance spectroscopy (¹⁹F NMR) is used to delineate key functional states of the adenosine A_2A receptor (A_2AR) complexed with heterotrimeric G protein (Gα_sβ₁γ₂) in a phospholipid membrane milieu. Analysis of A_2AR spectra as a function of ligand, G protein, and nucleotide identifies an ensemble represented by inactive states, a G-protein-bound activation intermediate, and distinct nucleotide-free states associated with either partial- or full-agonist-driven activation. The Gβγ subunit is found to be critical in facilitating ligand-dependent allosteric transmission, as shown by ¹⁹F NMR, biochemical, and computational studies. The results provide a mechanistic basis for understanding basal signaling, efficacy, precoupling, and allostery in GPCRs.

The Open Question of How GPCRs Interact with GPCR Kinases (GRKs)

G protein-coupled receptors (GPCRs), which regulate a vast number of eukaryotic processes, are desensitized by various mechanisms but, most importantly, by the GPCR kinases (GRKs). Ever since GRKs were first identified, investigators have sought to determine which structural features of GRKs are used to select for the agonist-bound states of GPCRs and how this binding event in turn enhances GRK catalytic activity. Despite a wealth of molecular information from high-resolution crystal structures of GRKs, the mechanisms driving activation have remained elusive, in part because the GRK N-terminus and active site tether region, previously proposed to serve as a receptor docking site and to be key to kinase domain closure, are often disordered or adopt inconsistent conformations. However, two recent studies have implicated other regions of GRKs as being involved in direct interactions with active GPCRs. Atomic resolution structures of GPCR–GRK complexes would help refine these models but are, so far, lacking. Here, we assess three distinct models for how GRKs recognize activated GPCRs, discuss limitations in the approaches used to generate them, and then experimentally test a hypothetical GPCR interaction site in GRK2 suggested by the two newest models.

Efficient G protein coupling is not required for agonist-mediated internalization and membrane reorganization of the adenosine A3 receptor

Organization of G protein‐coupled receptors at the plasma membrane has been the focus of much recent attention. Advanced microscopy techniques have shown that these receptors can be localized to discrete microdomains and reorganization upon ligand activation is crucial in orchestrating their signaling. Here, we have compared the membrane organization and downstream signaling of a mutant (R108A, R3.50A) of the adenosine A₃ receptor (A₃AR) to that of the wild‐type receptor. Fluorescence Correlation Spectroscopy (FCS) studies with a fluorescent agonist (ABEA‐X‐BY630) demonstrated that both wild‐type and mutant receptors bind agonist with high affinity but in subsequent downstream signaling assays the R108A mutation abolished agonist‐mediated inhibition of cAMP production and ERK phosphorylation. In further FCS studies, both A₃AR and A₃AR R108A underwent similar agonist‐induced increases in receptor density and molecular brightness which were accompanied by a decrease in membrane diffusion after agonist treatment. Using bimolecular fluorescence complementation, experiments showed that the R108A mutant retained the ability to recruit β‐arrestin and these receptor/arrestin complexes displayed similar membrane diffusion and organization to that observed with wild‐type receptors. These data demonstrate that effective G protein signaling is not a prerequisite for agonist‐stimulated β‐arrestin recruitment and membrane reorganization of the A₃AR.

Fungal G-Protein-Coupled Receptors: A Promising Mediator of the Impact of Extracellular Signals on Biosynthesis of Ochratoxin A

G-protein-coupled receptors (GPCRs) are transmembrane receptors involved in transducing signals from the external environment inside the cell, which enables fungi to coordinate cell transport, metabolism, and growth to promote their survival, reproduction, and virulence. There are 14 classes of GPCRs in fungi involved in sensing various ligands. In this paper, the synthesis of mycotoxins that are GPCR-mediated is discussed with respect to ligands, environmental stimuli, and intra-/interspecific communication. Despite their apparent importance in fungal biology, very little is known about the role of ochratoxin A (OTA) biosynthesis by Aspergillus ochraceus and the ligands that are involved. Fortunately, increasing evidence shows that the GPCR that involves the AF/ST (sterigmatocystin) pathway in fungi belongs to the same genus. Therefore, we speculate that GPCRs play an important role in a variety of environmental signals and downstream pathways in OTA biosynthesis. The verification of this inference will result in a more controllable GPCR target for control of fungal contamination in the future.

G protein-coupled receptor-G protein interactions: a single-molecule perspective

G protein-coupled receptors (GPCRs) regulate many cellular and physiological processes, responding to a diverse range of extracellular stimuli including hormones, neurotransmitters, odorants, and light. Decades of biochemical and pharmacological studies have provided fundamental insights into the mechanisms of GPCR signaling. Thanks to recent advances in structural biology, we now possess an atomistic understanding of receptor activation and G protein coupling. However, how GPCRs and G proteins interact in living cells to confer signaling efficiency and specificity remains insufficiently understood. The development of advanced optical methods, including single-molecule microscopy, has provided the means to study receptors and G proteins in living cells with unprecedented spatio-temporal resolution. The results of these studies reveal an unexpected level of complexity, whereby GPCRs undergo transient interactions among themselves as well as with G proteins and structural elements of the plasma membrane to form short-lived signaling nanodomains that likely confer both rapidity and specificity to GPCR signaling. These findings may provide new strategies to pharmaceutically modulate GPCR function, which might eventually pave the way to innovative drugs for common diseases such as diabetes or heart failure.

Alanine Scanning Mutagenesis of the DRYxxI Motif and Intracellular Loop 2 of Human Melanocortin-4 Receptor

The melanocortin-4 receptor (MC4R) is a member of the G-protein-coupled receptor (GPCR) superfamily, which has been extensively studied in obesity pathogenesis due to its critical role in regulating energy homeostasis. Both the Gs-cAMP and ERK1/2 cascades are known as important intracellular signaling pathways initiated by the MC4R. The DRYxxI motif at the end of transmembrane domain 3 and the intracellular loop 2 (ICL2) are thought to be crucial for receptor function in several GPCRs. To study the functions of this domain in MC4R, we performed alanine-scanning mutagenesis on seventeen residues. We showed that one residue was critical for receptor cell surface expression. Eight residues were important for ligand binding. Mutations of three residues impaired Gs-cAMP signaling without changing the binding properties. Investigation on constitutive activities of all the mutants in the cAMP pathway revealed that six residues were involved in constraining the receptor in inactive states and five residues were important for receptor activation in the absence of an agonist. In addition, mutations of four residues impaired the ligand-stimulated ERK1/2 signaling pathway without affecting the binding properties. We also showed that some mutants were biased to the Gs-cAMP or ERK1/2 signaling pathway. In summary, we demonstrated that the DRYxxI motif and ICL2 were important for MC4R function.

GPCRmd uncovers the dynamics of the 3D-GPCRome

G-protein-coupled receptors (GPCRs) are involved in numerous physiological processes and are the most frequent targets of approved drugs. The explosion in the number of new three-dimensional (3D) molecular structures of GPCRs (3D-GPCRome) over the last decade has greatly advanced the mechanistic understanding and drug design opportunities for this protein family. Molecular dynamics (MD) simulations have become a widely established technique for exploring the conformational landscape of proteins at an atomic level. However, the analysis and visualization of MD simulations require efficient storage resources and specialized software. Here we present GPCRmd (http://gpcrmd.org/), an online platform that incorporates web-based visualization capabilities as well as a comprehensive and user-friendly analysis toolbox that allows scientists from different disciplines to visualize, analyze and share GPCR MD data. GPCRmd originates from a community-driven effort to create an open, interactive and standardized database of GPCR MD simulations.

Distinctive communication networks in inactive states of β2 -adrenergic receptor: Mutual information and entropy transfer analysis

Mutual information and entropy transfer analysis employed on two inactive states of human beta-2 adrenergic receptor (β₂ -AR) unraveled distinct communication pathways. Previously, a so-called "highly" inactive state of the receptor was observed during 1.5 microsecond long molecular dynamics simulation where the largest intracellular loop (ICL3) was swiftly packed onto the G-protein binding cavity, becoming entirely inaccessible. Mutual information quantifying the degree of correspondence between backbone-C_α fluctuations was mostly shared between intra- and extra-cellular loop regions in the original inactive state, but shifted to entirely different regions in this latest inactive state. Interestingly, the largest amount of mutual information was always shared among the mobile regions. Irrespective of the conformational state, polar residues always contributed more to mutual information than hydrophobic residues, and also the number of polar-polar residue pairs shared the highest degree of mutual information compared to those incorporating hydrophobic residues. Entropy transfer, quantifying the correspondence between backbone-C_α fluctuations at different timesteps, revealed a distinctive pathway directed from the extracellular site toward intracellular portions in this recently exposed inactive state for which the direction of information flow was the reverse of that observed in the original inactive state where the mobile ICL3 and its intracellular surroundings drove the future fluctuations of extracellular regions.

Copper-mediated thiol potentiation and mutagenesis-guided modeling suggest a highly conserved copper-binding motif in human OR2M3

Sulfur-containing compounds within a physiological relevant, natural odor space, such as the key food odorants, typically constitute the group of volatiles with the lowest odor thresholds. The observation that certain metals, such as copper, potentiate the smell of sulfur-containing, metal-coordinating odorants led to the hypothesis that their cognate receptors are metalloproteins. However, experimental evidence is sparse-so far, only one human odorant receptor, OR2T11, and a few mouse receptors, have been reported to be activated by sulfur-containing odorants in a copper-dependent way, while the activation of other receptors by sulfur-containing odorants did not depend on the presence of metals. Here we identified an evolutionary conserved putative copper interaction motif CC/CSSH, comprising two copper-binding sites in TMH5 and TMH6, together with the binding pocket for 3-mercapto-2-methylpentan-1-ol in the narrowly tuned human receptor OR2M3. To characterize the copper-binding motif, we combined homology modeling, docking studies, site-directed mutagenesis, and functional expression of recombinant ORs in a cell-based, real-time luminescence assay. Ligand activation of OR2M3 was potentiated in the presence of copper. This effect of copper was mimicked by ionic and colloidal silver. In two broadly tuned receptors, OR1A1 and OR2W1, which did not reveal a putative copper interaction motif, activation by their most potent, sulfur-containing key food odorants did not depend on the presence of copper. Our results suggest a highly conserved putative copper-binding motif to be necessary for a copper-modulated and thiol-specific function of members from three subfamilies of family 2 ORs.

A Conserved Proline Hinge Mediates Helix Dynamics and Activation of Rhodopsin

Despite high-resolution crystal structures of both inactive and active G protein-coupled receptors (GPCRs), it is still not known how ligands trigger the large structural change on the intracellular side of the receptor since the conformational changes that occur within the extracellular ligand-binding region upon activation are subtle. Here, we use solid-state NMR and Fourier transform infrared spectroscopy on rhodopsin to show that Trp265^6.48 within the CWxP motif on transmembrane helix H6 constrains a proline hinge in the inactive state, suggesting that activation results in unraveling of the H6 backbone within this motif, a local change in dynamics that allows helix H6 to swing outward. Notably, Tyr301^7.48 within activation switch 2 appears to mimic the negative allosteric sodium ion found in other family A GPCRs, a finding that is broadly relevant to the mechanism of receptor activation.

GIP receptor: Expression in neuroendocrine tumours, internalization, signalling from endosomes and structure-function relationship studies

GIP is well known as a peptide regulating metabolic functions. In this review paper, we summarize a series of data on GIP receptor (GIPR). First, expression study of GIPR in human neuroendocrine tumours showed a very high incidence (nearly 100%) and a high density in both functional and non functional pancreatic tumours, ileal tumours, bronchial tumours and medullary thyroid carcinomas. Then, data on internalization of GIPR following stimulation by GIP are reported. Rapid and abundant internalization of GIPR also found in tumor pancreatic endocrine cells opens the possibility of tumor imaging and eradication using radiolabeled GIP. Interestingly, internalized GIPR continues to signal in early endosomes stimulating production of cAMP and activation of PKA, thus, supporting the view that GIPR signals from both plasma membrane and vesicles of internalization. At last, we summarize data from studies using in synergy molecular modeling and site-directed mutagenesis, which identified crucial amino acids of transmembrane domains of GIPR involved in GIPR binding site of GIP and/or in its activation and coupling to Gs protein. All together, these last molecular data may help to better understand structure-activity relationship data on GIP and GIPR.

Developmental regulation of olfactory circuit formation in mice

In mammals, odorants induce various behavioral responses that are critical to the survival of the individual and species. Binding signals of odorants to odorant receptors (ORs) expressed in the olfactory epithelia are converted to an odor map, a pattern of activated glomeruli, in the olfactory bulb (OB). This topographic map is used to identify odorants for memory‐based learned decisions. In the embryo, a coarse olfactory map is generated in the OB by a combination of dorsal‐ventral and anterior‐posterior targeting of olfactory sensory neurons (OSNs), using specific sets of axon‐guidance molecules. During the process of OSN projection, odor signals are sorted into distinct odor qualities in separate functional domains in the OB. Odor information is then conveyed by the projection neurons, mitral/tufted cells, to various regions in the olfactory cortex, particularly to the amygdala for innate olfactory decisions. Although the basic architecture of hard‐wired circuits is generated by a genetic program, innate olfactory responses are modified by neonatal odor experience in an activity‐dependent manner. Stimulus‐driven OR activity promotes post‐synaptic events and dendrite selection in the responding glomeruli making them larger. As a result, enhanced odor inputs in neonates establish imprinted olfactory memory that induces attractive responses in adults, even when the odor quality is innately aversive. In this paper, I will provide an overview of the recent progress made in the olfactory circuit formation in mice.

Structural basis of ligand recognition and self-activation of orphan GPR52

GPR52 is a class-A orphan G-protein-coupled receptor that is highly expressed in the brain and represents a promising therapeutic target for the treatment of Huntington's disease and several psychiatric disorders^1,2. Pathological malfunction of GPR52 signalling occurs primarily through the heterotrimeric G_s protein², but it is unclear how GPR52 and G_s couple for signal transduction and whether a native ligand or other activating input is required. Here we present the high-resolution structures of human GPR52 in three states: a ligand-free state, a G_s-coupled self-activation state and a potential allosteric ligand-bound state. Together, our structures reveal that extracellular loop 2 occupies the orthosteric binding pocket and operates as a built-in agonist, conferring an intrinsically high level of basal activity to GPR52³. A fully active state is achieved when G_s is coupled to GPR52 in the absence of an external agonist. The receptor also features a side pocket for ligand binding. These insights into the structure and function of GPR52 could improve our understanding of other self-activated GPCRs, enable the identification of endogenous and tool ligands, and guide drug discovery efforts that target GPR52.

GPR6 Structural Insights: Homology Model Construction and Docking Studies

GPR6 is an orphan G protein-coupled receptor that has been associated with the cannabinoid family because of its recognition of a sub-set of cannabinoid ligands. The high abundance of GPR6 in the central nervous system, along with high constitutive activity and a link to several neurodegenerative diseases make GPR6 a promising biological target. In fact, diverse research groups have demonstrated that GPR6 represents a possible target for the treatment of neurodegenerative disorders such as Parkinson's disease, Alzheimer's disease, and Huntington's disease. Several patents have claimed the use of a wide range of pyrazine derivatives as GPR6 inverse agonists for the treatment of Parkinson's disease symptoms and other dyskinesia syndromes. However, the full pharmacological importance of GPR6 has not yet been fully explored due to the lack of high potency, readily available ligands targeting GPR6. The long-term goal of the present study is to develop such ligands. In this paper, we describe our initial steps towards this goal. A human GPR6 homology model was constructed using a suite of computational techniques. This model permitted the identification of unique GPR6 structural features and the exploration of the GPR6 binding crevice. A subset of patented pyrazine analogs were docked in the resultant GPR6 inactive state model to validate the model, rationalize the structure-activity relationships from the reported patents and identify the key residues in the binding crevice for ligand recognition. We will take this structural knowledge into the next phase of GPR6 project, in which scaffold hopping will be used to design new GPR6 ligands.