Mutagenesis and modelling of the alpha(1b)-adrenergic receptor highlight the role of the helix 3/helix 6 interface in receptor activation

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.

Binding and Activation of Serotonergic G-Protein Coupled Receptors by the Multimodal Antidepressant Vortioxetine

G-protein coupled receptors (GPCRs) are important pharmacological targets. Despite substantial progress, important questions still remain concerning the details of activation: how can a ligand act as an agonist in one receptor but as an antagonist in a homologous receptor, and how can agonists activate a receptor despite lacking polar functional groups able to interact with helix 5 as is the case for the related adrenergic receptors? Studying vortioxetine (VXT), an important multimodal antidepressant drug, may elucidate both questions. Herein, we present a thorough in silico analysis of VXT binding to 5-HT_1A, 5-HT_1B, and 5-HT₇ receptors and compare it with available experimental data. We are able to rationalize the differential mode of action of VXT at different receptors, but also, in the case of the 5-HT_1A receptor, we observe the initial steps of activation that inform about an activation mechanism that does not involve polar interaction with helix 5. The results extend our current understanding of agonist and antagonist action at aminergic GPCRs.

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.

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.

Structural and functional analysis of Ccr1l1, a Rodentia-restricted eosinophil-selective chemokine receptor homologue

Mouse Ccr1l1 (Ccr1-like 1) encodes an orphan G-protein-coupled receptor (GPCR) with the highest homology to the inflammatory and highly promiscuous chemokine receptors Ccr1 and Ccr3 (70 and 50% amino acid identity, respectively). Ccr1l1 was first cloned in 1995, yet current knowledge of this putative chemokine receptor is limited to its gene organization and chromosomal localization. Here we report that Ccr1l1 is a Rodentia-specific gene selectively expressed in eosinophils. However, eosinophil phenotypes, development, and responsiveness to chemokines were all normal in naïve Ccr1l1 knockout mice. We demonstrate for the first time that recombinant Ccr1l1 is expressed on the plasma membrane of transfected cells and contains an extracellular N terminus and an intracellular C terminus, consistent with GPCR topology. Using receptor internalization, β-arrestin recruitment, calcium flux, and chemotaxis assays, we excluded all 37 available mouse chemokines, including Ccr1 ligands, and two viral chemokines as Ccr1l1 ligands, and demonstrated that mouse Ccr1, but not Ccr1l1, exhibits constitutive signaling activity. However, sequence analysis and structural modeling revealed that Ccr1l1 is well equipped to act as a classical signaling GPCR, with N-terminal sulfotyrosines as the only signaling and chemokine-binding determinant absent in Ccr1l1. Hereof, we show that a sulfatable N-terminal Ccr1 Y¹⁸ residue is essential for chemotaxis and calcium responses induced by Ccl3 and Ccl9/10, but substituting the corresponding Ccr1l1 F¹⁹ residue with tyrosine failed to confer responsiveness to Ccr1 ligands. Although Ccr1l1 remains an extreme outlier in the chemokine receptor family, our study supports that it might respond to unidentified mouse chemokine ligands in eosinophil-driven immune responses.

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.

Structural features of activated GPCR signaling complexes

G protein-coupled receptors (GPCRs) couple to diverse heterotrimeric G protein subtypes and then activate downstream signaling pathways in classical GPCR activation. It has also been found that GPCRs transduce signals through different regulatory proteins, such as arrestins. Recently, owing to the breakthroughs in cryo-electron macroscopy (Cryo-EM), numerous structures of GPCR-G protein or GPCR-arrestin complexes have been deciphered. In this review, we summarize most of reported GPCR signaling complex structures, with an emphasis on the structural features of rhodopsin-like GPCR activation and G protein-binding/arrestin-binding modes, to illustrate the activation and signaling mechanism of rhodopsin-like GPCRs.

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.

Mutations in the NPxxY motif stabilize pharmacologically distinct conformational states of the α1B- and β2-adrenoceptors

G protein-coupled receptors (GPCRs) convert extracellular stimuli to intracellular responses that regulate numerous physiological processes. Crystallographic and biophysical advances in GPCR structural analysis have aided investigations of structure-function relationships that clarify our understanding of these dynamic receptors, but the molecular mechanisms associated with activation and signaling for individual GPCRs may be more complex than was previously appreciated. Here, we investigated the proposed water-mediated, hydrogen-bonded activation switch between the conserved NPxxY motif on transmembrane helix 7 (TMH7) and a conserved tyrosine in TMH5, which contributes to α_1B-adrenoceptor (α_1B-AR) and β₂-AR activation. Disrupting this bond by mutagenesis stabilized the α_1B-AR and the β₂-AR in inactive-state conformations, which displayed decreased agonist potency for stimulating downstream IP₁ and cAMP signaling, respectively. Compared to that for wild-type receptors, agonist-mediated β-arrestin recruitment was substantially reduced or abolished for all α_1B-AR and β₂-AR inactive-state mutants. However, the inactive-state β₂-ARs exhibited decreased agonist affinity, whereas the inactive-state α_1B-ARs had enhanced agonist affinity. Conversely, antagonist affinity was unchanged for inactive-state conformations of both α_1B-AR and β₂-AR. Removing the influence of agonist affinity on agonist potency gave a measure of signaling efficacy, which was markedly decreased for the α_1B-AR mutants but little altered for the β₂-AR mutants. These findings highlight the pharmacological heterogeneity of inactive-state GPCR conformations, which may facilitate the rational design of drugs that target distinct conformational states of GPCRs.

Biased Agonist TRV027 Determinants in AT1R by Molecular Dynamics Simulations

Functional selectivity is a phenomenon observed in G protein-coupled receptors in which intermediate active-state conformations are stabilized by mutations or ligand binding, resulting in different sets of signaling pathways. Peptides capable of selectively activating β-arrestin, known as biased agonists, have already been characterized in vivo and could correspond to a new therapeutic approach for treatment of cardiovascular diseases. Despite the potential of biased agonism, the mechanism involved in selective signaling remains unclear. In this work, molecular dynamics simulations were employed to compare the conformational profile of the angiotensin II type 1 receptor (AT1R) crystal bound to angiotensin II, bound to the biased ligand TRV027, and in the apo form. Our results show that both ligands induce changes near the NPxxY motif in transmembrane domain 7 that are related to receptor activation. However, the biased ligand does not cause the rotamer toggle alternative positioning and displays an exclusive hydrogen-bonding pattern. Our work sheds light on the biased agonism mechanism and will help in the future design of novel biased agonists for AT1R.

The Role of the CXCL12/CXCR4/ACKR3 Axis in Autoimmune Diseases

Chemokine receptors are members of the G protein-coupled receptor superfamily. These receptors are intimately involved in cell movement, and thus play a critical role in several physiological and pathological situations that require the precise regulation of cell positioning. CXCR4 is one of the most studied chemokine receptors and is involved in many functions beyond leukocyte recruitment. During embryogenesis, it plays essential roles in vascular development, hematopoiesis, cardiogenesis, and nervous system organization. It has been also implicated in tumor progression and autoimmune diseases and, together with CD4, is one of the co-receptors used by the HIV-1 virus to infect immune cells. In contrast to other chemokine receptors that are characterized by ligand promiscuity, CXCR4 has a unique ligand-stromal cell-derived factor-1 (SDF1, CXCL12). However, this ligand also binds ACKR3, an atypical chemokine receptor that modulates CXCR4 functions and is overexpressed in multiple cancer types. The CXCL12/CXCR4/ACKR3 axis constitutes a potential therapeutic target for a wide variety of inflammatory diseases, not only by interfering with cell migration but also by modulating immune responses. Thus far, only one antagonist directed against the ligand-binding site of CXCR4, AMD3100, has demonstrated clinical relevance. Here, we review the role of this ligand and its receptors in different autoimmune diseases.

Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

The aminergic family of G protein-coupled receptors (GPCRs) plays an important role in various diseases and represents a major drug discovery target class. Structure determination of all major aminergic subfamilies has enabled structure-based ligand design for these receptors. Site-directed mutagenesis data provides an invaluable complementary source of information for elucidating the structural determinants of binding of different ligand chemotypes. The current study provides a comparative analysis of 6692 mutation data points on 34 aminergic GPCR subtypes, covering the chemical space of 540 unique ligands from mutagenesis experiments and information from experimentally determined structures of 52 distinct aminergic receptor-ligand complexes. The integrated analysis enables detailed investigation of structural receptor-ligand interactions and assessment of the transferability of combined binding mode and mutation data across ligand chemotypes and receptor subtypes. An overview is provided of the possibilities and limitations of using mutation data to guide the design of novel aminergic receptor ligands.

A Comprehensive Mutagenesis Screen of the Adhesion GPCR Latrophilin-1/ADGRL1

Adhesion G-protein-coupled receptors (aGPCRs) play critical roles in diverse cellular processes in neurobiology, development, immunity, and numerous diseases. The lack of molecular understanding of their activation mechanisms, especially with regard to the transmembrane domains, hampers further studies to facilitate aGPCR-targeted drug development. Latrophilin-1/ADGRL1 is a model aGPCR that regulates synapse formation and embryogenesis, and its mutations are associated with cancer and attention-deficit/hyperactivity disorder. Here, we established functional assays to monitor latrophilin-1 function and showed the activation of latrophilin-1 by its endogenous agonist peptide. Via a comprehensive mutagenesis screen, we identified transmembrane domain residues essential for latrophilin-1 basal activity and for agonist peptide response. Strikingly, a cancer-associated mutation exhibited increased basal activity and failed to rescue the embryonic developmental phenotype in transgenic worms. These results provide a mechanistic foundation for future aGPCR-targeted drug design.

Hybrid approach to structure modeling of the histamine H3 receptor: Multi-level assessment as a tool for model verification

The crucial role of G-protein coupled receptors and the significant achievements associated with a better understanding of the spatial structure of known receptors in this family encouraged us to undertake a study on the histamine H3 receptor, whose crystal structure is still unresolved. The latest literature data and availability of different software enabled us to build homology models of higher accuracy than previously published ones. The new models are expected to be closer to crystal structures; and therefore, they are much more helpful in the design of potential ligands. In this article, we describe the generation of homology models with the use of diverse tools and a hybrid assessment. Our study incorporates a hybrid assessment connecting knowledge-based scoring algorithms with a two-step ligand-based docking procedure. Knowledge-based scoring employs probability theory for global energy minimum determination based on information about native amino acid conformation from a dataset of experimentally determined protein structures. For a two-step docking procedure two programs were applied: GOLD was used in the first step and Glide in the second. Hybrid approaches offer advantages by combining various theoretical methods in one modeling algorithm. The biggest advantage of hybrid methods is their intrinsic ability to self-update and self-refine when additional structural data are acquired. Moreover, the diversity of computational methods and structural data used in hybrid approaches for structure prediction limit inaccuracies resulting from theoretical approximations or fuzziness of experimental data. The results of docking to the new H3 receptor model allowed us to analyze ligand-receptor interactions for reference compounds.

Structural determinants at the M2 muscarinic receptor modulate the RGS4-GIRK response to pilocarpine by impairment of the receptor voltage sensitivity

Membrane potential controls the response of the M2 muscarinic receptor to its ligands. Membrane hyperpolarization increases response to the full agonist acetylcholine (ACh) while decreasing response to the partial agonist pilocarpine. We previously have demonstrated that the regulator of G-protein signaling (RGS) 4 protein discriminates between the voltage-dependent responses of ACh and pilocarpine; however, the underlying mechanism remains unclear. Here we show that RGS4 is involved in the voltage-dependent behavior of the M2 muscarinic receptor-mediated signaling in response to pilocarpine. Additionally we revealed structural determinants on the M2 muscarinic receptor underlying the voltage-dependent response. By electrophysiological recording in Xenopus oocytes expressing M2 muscarinic receptor and G-protein-gated inwardly rectifying K⁺ channels, we quantified voltage-dependent desensitization of pilocarpine-induced current in the presence or absence of RGS4. Hyperpolarization-induced desensitization of the current required for RGS4, also depended on pilocarpine concentration. Mutations of charged residues in the aspartic acid-arginine-tyrosine motif of the M2 muscarinic receptor, but not intracellular loop 3, significantly impaired the voltage-dependence of RGS4 function. Thus, our results demonstrated that voltage-dependence of RGS4 modulation is derived from the M2 muscarinic receptor. These results provide novel insights into how membrane potential impacts G-protein signaling by modulating GPCR communication with downstream effectors.

The DRY motif and the four corners of the cubic ternary complex model

Recent structural data on GPCRs using a variety of spectroscopic approaches suggest that GPCRs adopt a dynamic conformational landscape, with ligands stabilizing subsets of these states to activate one or more downstream signaling effectors. A key outstanding question posed by this emerging dynamic structural model of GPCRs is what states, active, inactive, or intermediate are captured by the numerous crystal structures of GPCRs complexed with a variety of agonists, partial agonists, and antagonists. In the early nineties the discovery of inverse agonists and constitutive activity led to the idea that the active receptor state (R^⁎) is an intrinsic property of the receptor itself rather than of the RG complex, eventually leading to the formulation of the cubic ternary complex model (CTC). Here, by a careful analysis of a series of data obtained with a number of mutants of the highly conserved E/DRY motif, we show evidences for the existence of all the receptor states theorized by the CTC, four 'uncoupled (R, R^⁎ and HR and HR^⁎), and, consequently four 'coupled' (RG, R^⁎G, HRG and HR^⁎G). The E/DRY motif located at the cytosolic end of transmembrane helix III of Class A GPCRs has been widely studied and analyzed because it forms a network of interactions believed to lock receptors in the inactive conformation (R), and, thus, to play a key role in receptor activation. Our conclusions are supported by recent crystal and NMR spectra, as well as by results obtained with two prototypical GPCRs using a new FRET technology that de-couples G protein binding to the receptor from signal transduction. Thus, despite its complexity and limitations, we propose that the CTC is a useful framework to reconcile pharmacological, biochemical and structural data.

Conserved Mechanism of Conformational Stability and Dynamics in G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are transmembrane receptors involved in diverse biological functions. Despite the diversity in their amino acid sequences, class A GPCRs exhibit a conserved structural topology and possibly a common mechanism of receptor activation. To understand how this high sequence diversity translates to a conserved functional mechanism, we have compared the dynamic behavior of eight class A GPCRs comprised of six biogenic amine receptors, adenosine A_2A, and the peptide receptor protease-activated receptor 1. Starting from the crystal structures of the inactive state of these receptors bound to inverse agonists or antagonists, we have performed multiple all-atom MD simulations adding up to several microseconds of simulation. We elucidated the similarities and differences in the dynamic behavior and the conformational ensembles sampled by these eight class A GPCRs. Among the six biogenic amine receptors studied here, β₂-adrenergic receptor shows the highest level of fluctuation in the sixth and seventh transmembrane helices, possibly explaining its high basal activity. In contrast, the muscarinic acetylcholine receptors show the lowest fluctuations as well as tight packing and low hydration of the transmembrane domain. All eight GPCRs show several conserved allosteric communication pipelines from the residues in the agonist binding site with the G-protein interface. Positions of the residues along these pipelines that serve as major hubs of allosteric communication are conserved in their respective structures. These findings have important implications in understanding the dynamics and allosteric mechanism of communication in class A GPCRs and hence are useful for designing conformation-specific drugs.

Signal transmission through the CXC chemokine receptor 4 (CXCR4) transmembrane helices

The atomic-level mechanisms by which G protein-coupled receptors (GPCRs) transmit extracellular ligand binding events through their transmembrane helices to activate intracellular G proteins remain unclear. Using a comprehensive library of mutations covering all 352 residues of the GPCR CXC chemokine receptor 4 (CXCR4), we identified 41 amino acids that are required for signaling induced by the chemokine ligand CXCL12 (stromal cell-derived factor 1). CXCR4 variants with each of these mutations do not signal properly but remain folded, based on receptor surface trafficking, reactivity to conformationally sensitive monoclonal antibodies, and ligand binding. When visualized on the structure of CXCR4, the majority of these residues form a continuous intramolecular signaling chain through the transmembrane helices; this chain connects chemokine binding residues on the extracellular side of CXCR4 to G protein-coupling residues on its intracellular side. Integrated into a cohesive model of signal transmission, these CXCR4 residues cluster into five functional groups that mediate (i) chemokine engagement, (ii) signal initiation, (iii) signal propagation, (iv) microswitch activation, and (v) G protein coupling. Propagation of the signal passes through a "hydrophobic bridge" on helix VI that coordinates with nearly every known GPCR signaling motif. Our results agree with known conserved mechanisms of GPCR activation and significantly expand on understanding the structural principles of CXCR4 signaling.

Diverse activation pathways in class A GPCRs converge near the G-protein-coupling region

Class A G-protein-coupled receptors (GPCRs) are a large family of membrane proteins that mediate a wide variety of physiological functions, including vision, neurotransmission and immune responses. They are the targets of nearly one-third of all prescribed medicinal drugs such as beta blockers and antipsychotics. GPCR activation is facilitated by extracellular ligands and leads to the recruitment of intracellular G proteins. Structural rearrangements of residue contacts in the transmembrane domain serve as 'activation pathways' that connect the ligand-binding pocket to the G-protein-coupling region within the receptor. In order to investigate the similarities in activation pathways across class A GPCRs, we analysed 27 GPCRs from diverse subgroups for which structures of active, inactive or both states were available. Here we show that, despite the diversity in activation pathways between receptors, the pathways converge near the G-protein-coupling region. This convergence is mediated by a highly conserved structural rearrangement of residue contacts between transmembrane helices 3, 6 and 7 that releases G-protein-contacting residues. The convergence of activation pathways may explain how the activation steps initiated by diverse ligands enable GPCRs to bind a common repertoire of G proteins.

DIRECT-ID: An automated method to identify and quantify conformational variations--application to β2 -adrenergic GPCR

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in environment, ligand binding, and interactions with other macromolecules, among others. We present a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the norms of the differences in covariance matrices are calculated for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. For each pair of trajectories, covariance difference matrices are parsed to extract structural elements that undergo changes in conformational properties. As a demonstration of its applicability to biomacromolecular systems, the method, referred to as DIRECT-ID, was used to identify relevant ligand-modulated structural variations in the β2 -adrenergic (β2 AR) G-protein coupled receptor. Micro-second MD simulations of the β2 AR in an explicit lipid bilayer were run in the apo state and complexed with the ligands: BI-167107 (agonist), epinephrine (agonist), salbutamol (long-acting partial agonist), or carazolol (inverse agonist). Each ligand modulated the conformational dynamics of β2 AR differently and DIRECT-ID analysis of the inverse-agonist vs. agonist-modulated β2 AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand-specific microswitches across the GPCR. This study demonstrates the utility of DIRECT-ID to rapidly extract functionally relevant conformational dynamics information from extended MD simulations of large and complex macromolecular systems.

Differences in allosteric communication pipelines in the inactive and active states of a GPCR

G-protein-coupled receptors (GPCRs) are membrane proteins that allosterically transduce the signal of ligand binding in the extracellular (EC) domain to couple to proteins in the intracellular (IC) domain. However, the complete pathway of allosteric communication from the EC to the IC domain, including the role of individual amino acids in the pathway is not known. Using the correlation in torsion angle movements calculated from microseconds-long molecular-dynamics simulations, we elucidated the allosteric pathways in three different conformational states of β2-adrenergic receptor (β2AR): 1), the inverse-agonist-bound inactive state; 2), the agonist-bound intermediate state; and (3), the agonist- and G-protein-bound fully active state. The inactive state is less dynamic compared with the intermediate and active states, showing dense clusters of allosteric pathways (allosteric pipelines) connecting the EC with the IC domain. The allosteric pipelines from the EC domain to the IC domain are weakened in the intermediate state, thus decoupling the EC domain from the IC domain and making the receptor more dynamic compared with the other states. Also, the orthosteric ligand-binding site becomes the initiator region for allosteric communication in the intermediate state. This finding agrees with the paradigm that the nature of the agonist governs the specific signaling state of the receptor. These results provide an understanding of the mechanism of allosteric communication in class A GPCRs. In addition, our analysis shows that mutations that affect the ligand efficacy, but not the binding affinity, are located in the allosteric pipelines. This clarifies the role of such mutations, which has hitherto been unexplained.

Phenotypic regulation of the sphingosine 1-phosphate receptor miles apart by G protein-coupled receptor kinase 2

The evolutionarily conserved DRY motif at the end of the third helix of rhodopsin-like, class-A G protein-coupled receptors (GPCRs) is a major regulator of receptor stability, signaling activity, and β-arrestin-mediated internalization. Substitution of the DRY arginine with histidine in the human vasopressin receptor results in a loss-of-function phenotype associated with diabetes insipidus. The analogous R150H substitution of the DRY motif in zebrafish sphingosine-1 phosphate receptor 2 (S1p2) produces a mutation, miles apart m⁹³ (mil^m93), that not only disrupts signaling but also impairs heart field migration. We hypothesized that constitutive S1p2 desensitization is the underlying cause of this strong zebrafish developmental defect. We observed in cell assays that the wild-type S1p2 receptor is at the cell surface whereas in distinct contrast the S1p2 R150H receptor is found in intracellular vesicles, blocking G protein but not arrestin signaling activity. Surface S1p2 R150H expression could be restored by inhibition of G protein-coupled receptor kinase 2 (GRK2). Moreover, we observed that β-arrestin 2 and GRK2 colocalize with S1p2 in developing zebrafish embryos and depletion of GRK2 in the S1p2 R150H miles apart zebrafish partially rescued cardia bifida. The ability of reduced GRK2 activity to reverse a developmental phenotype associated with constitutive desensitization supports efforts to genetically or pharmacologically target this kinase in diseases involving biased GPCR signaling.

Unifying family A GPCR theories of activation

Several new pairs of active and inactive GPCR structures have recently been solved enabling detailed structural insight into the activation process, not only of rhodopsin but now also of the β2 adrenergic, M2 muscarinic and adenosine A2A receptors. Combined with structural analyses they have enabled us to examine the different recent theories proposed for GPCR activation and show that they are all indeed parts of the same process, and are intrinsically related through their effect on the central hydrophobic core of GPCRs. This new unifying general process of activation is consistent with the identification of known constitutively active mutants and an in-depth conservational analysis of significant residues implicated in the process.

Molecular dynamics simulations on the Tre1 G protein-coupled receptor: exploring the role of the arginine of the NRY motif in Tre1 structure

BACKGROUND

The arginine of the D/E/NRY motif in Rhodopsin family G protein-coupled receptors (GPCRs) is conserved in 96% of these proteins. In some GPCRs, this arginine in transmembrane 3 can form a salt bridge with an aspartic acid or glutamic acid in transmembrane 6. The Drosophila melanogaster GPCR Trapped in endoderm-1 (Tre1) is required for normal primordial germ cell migration. In a mutant form of the protein, Tre1sctt, eight amino acids RYILIACH are missing, resulting in a severe disruption of primordial germ cell development. The impact of the loss of these amino acids on Tre1 structure is unknown. Since the missing amino acids in Tre1sctt include the arginine that is part of the D/E/NRY motif in Tre1, molecular dynamics simulations were performed to explore the hypothesis that these amino acids are involved in salt bridge formation and help maintain Tre1 structure.

RESULTS

Structural predictions of wild type Tre1 (Tre1+) and Tre1sctt were subjected to over 250 ns of molecular dynamics simulations. The ability of the model systems to form a salt bridge between the arginine of the D/E/NRY motif and an aspartic acid residue in transmembrane 6 was analyzed. The results indicate that a stable salt bridge can form in the Tre1+ systems and a weak salt bridge or no salt bridge, using an alternative arginine, is likely in the Tre1sctt systems.

CONCLUSIONS

The weak salt bridge or lack of a salt bridge in the Tre1sctt systems could be one possible explanation for the disrupted function of Tre1sctt in primordial germ cell migration. These results provide a framework for studying the importance of the arginine of the D/E/NRY motif in the structure and function of other GPCRs that are involved in cell migration, such as CXCR4 in the mouse, zebrafish, and chicken.

Conopeptide ρ-TIA defines a new allosteric site on the extracellular surface of the α1B-adrenoceptor

The G protein-coupled receptor (GPCR) superfamily is an important drug target that includes over 1000 membrane receptors that functionally couple extracellular stimuli to intracellular effectors. Despite the potential of extracellular surface (ECS) residues in GPCRs to interact with subtype-specific allosteric modulators, few ECS pharmacophores for class A receptors have been identified. Using the turkey β(1)-adrenergic receptor crystal structure, we modeled the α(1B)-adrenoceptor (α(1B)-AR) to help identify the allosteric site for ρ-conopeptide TIA, an inverse agonist at this receptor. Combining mutational radioligand binding and inositol 1-phosphate signaling studies, together with molecular docking simulations using a refined NMR structure of ρ-TIA, we identified 14 residues on the ECS of the α(1B)-AR that influenced ρ-TIA binding. Double mutant cycle analysis and docking confirmed that ρ-TIA binding was dominated by a salt bridge and cation-π between Arg-4-ρ-TIA and Asp-327 and Phe-330, respectively, and a T-stacking-π interaction between Trp-3-ρ-TIA and Phe-330. Water-bridging hydrogen bonds between Asn-2-ρ-TIA and Val-197, Trp-3-ρ-TIA and Ser-318, and the positively charged N terminus and Glu-186, were also identified. These interactions reveal that peptide binding to the ECS on transmembrane helix 6 (TMH6) and TMH7 at the base of extracellular loop 3 (ECL3) is sufficient to allosterically inhibit agonist signaling at a GPCR. The ligand-accessible ECS residues identified provide the first view of an allosteric inhibitor pharmacophore for α(1)-adrenoceptors and mechanistic insight and a new set of structural constraints for the design of allosteric antagonists at related GPCRs.

Action of molecular switches in GPCRs--theoretical and experimental studies

G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called "molecular switches" buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homoand heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.

Understanding functional residues of the cannabinoid CB1

The brain cannabinoid (CB(1)) receptor that mediates numerous physiological processes in response to marijuana and other psychoactive compounds is a G protein coupled receptor (GPCR) and shares common structural features with many rhodopsin class GPCRs. For the rational development of therapeutic agents targeting the CB(1) receptor, understanding of the ligand-specific CB(1) receptor interactions responsible for unique G protein signals is crucial. For a more than a decade, a combination of mutagenesis and computational modeling approaches has been successfully employed to study the ligand-specific CB(1) receptor interactions. In this review, after a brief discussion about recent advances in understanding of some structural and functional features of GPCRs commonly applicable to the CB(1) receptor, the CB(1) receptor functional residues reported from mutational studies are divided into three different types, ligand binding (B), receptor stabilization (S) and receptor activation (A) residues, to delineate the nature of the binding pockets of anandamide, CP55940, WIN55212-2 and SR141716A and to describe the molecular events of the ligand-specific CB(1) receptor activation from ligand binding to G protein signaling. Taken these CB(1) receptor functional residues, some of which are unique to the CB(1) receptor, together with the biophysical knowledge accumulated for the GPCR active state, it is possible to propose the early stages of the CB(1) receptor activation process that not only provide some insights into understanding molecular mechanisms of receptor activation but also are applicable for identifying new therapeutic agents by applying the validated structure-based approaches, such as virtual high throughput screening (HTS) and fragment-based approach (FBA).

Simplified modeling approach suggests structural mechanisms for constitutive activation of the C5a receptor

Molecular modeling of conformational changes occurring in the transmembrane region of the complement factor 5a receptor (C5aR) during receptor activation was performed by comparing two constitutively active mutants (CAMs) of C5aR, NQ (I124N/L127Q), and F251A, to those of the wild-type C5aR and NQ-N296A (I124N/L127Q/N296A), which have the wild-type phenotype. Modeling involved comprehensive sampling of various rotations of TM helices aligned to the crystal template of the dark-adapted rhodopsin along their long axes. By assuming that the relative energies of the spontaneously activated states of CAMs should be lower or at least comparable to energies characteristic for the ground states, we selected the plausible models for the conformational states associated with constitutive activation in C5aR. The modeling revealed that the hydrogen bonds between the side chains of D82-N119, S85-N119, and S131-C221 characteristic for the ground state were replaced by the hydrogen bonds D82-N296, N296-Y300, and S131-R134, respectively, in the activated states. Also, conformational transitions that occurred upon activation were hindered by contacts between the side chains of L127 and F251. The results rationalize the available data of mutagenesis in C5aR and offer the first specific molecular mechanism for the loss of constitutive activity in NQ-N296A. Our results also contributed to understanding the general structural mechanisms of activation in G-protein-coupled receptors lacking the "ionic lock", R(3.50) and E/D(6.30). Importantly, these results were obtained by modeling approaches that deliberately simplify many elements in order to explore potential conformations of GPCRs involving large-scale molecular movements.

The α1-adrenergic receptors: diversity of signaling networks and regulation

The α₁-adrenergic receptor (AR) subtypes (α_1a, α_1b, and α_1d) mediate several physiological effects of epinephrineand norepinephrine. Despite several studies in recombinant systems and insightfrom genetically modified mice, our understanding of the physiological relevance and specificity of the α₁-AR subtypes is still limited. Constitutive activity and receptor oligomerization have emerged as potential features regulating receptor function. Another recent paradigm is that βarrestins and G protein-coupled receptors themselves can act as scaffolds binding a variety of proteins and this can result in growing complexity of the receptor-mediated cellular effects. The aim of this review is to summarize our current knowledge on some recently identified functional paradigms and signaling networks that might help to elucidate the functional diversity of the α₁-AR subtypes in various organs.

Allosteric rescuing of loss-of-function FFAR2 mutations

FFAR2 (GPR43) is a receptor for short-chain fatty acids (SCFAs), acetate and propionate. In the current study, we investigate the molecular determinants contributing to receptor activation by endogenous ligands. Mutational analysis revealed several important residues located in transmembrane domains (TM) 3, 4, 5, 6, and 7 for acetate binding. Interestingly, mutations that abolished acetate activity, including the mutation in the well-conserved D(E)RY motif, could be rescued by a recently identified synthetic allosteric agonist. These findings provide additional insight into agonist binding and activation which may aid in designing allosteric ligands for targeting receptor function in various diseases.

An evolutionarily conserved arginine is essential for Tre1 G protein-coupled receptor function during germ cell migration in Drosophila melanogaster

BACKGROUND

G protein-coupled receptors (GPCRs) play central roles in mediating cellular responses to environmental signals leading to changes in cell physiology and behaviors, including cell migration. Numerous clinical pathologies including metastasis, an invasive form of cell migration, have been linked to abnormal GPCR signaling. While the structures of some GPCRs have been defined, the in vivo roles of conserved amino acid residues and their relationships to receptor function are not fully understood. Trapped in endoderm 1 (Tre1) is an orphan receptor of the rhodopsin class that is necessary for primordial germ cell migration in Drosophila melanogaster embryos. In this study, we employ molecular genetic approaches to identify residues in Tre1 that are critical to its functions in germ cell migration.

METHODOLOGY/PRINCIPAL FINDINGS

First, we show that the previously reported scattershot mutation is an allele of tre1. The scattershot allele results in an in-frame deletion of 8 amino acids at the junction of the third transmembrane domain and the second intracellular loop of Tre1 that dramatically impairs the function of this GPCR in germ cell migration. To further refine the molecular basis for this phenotype, we assayed the effects of single amino acid substitutions in transgenic animals and determined that the arginine within the evolutionarily conserved E/N/DRY motif is critical for receptor function in mediating germ cell migration within an intact developing embryo.

CONCLUSIONS/SIGNIFICANCE

These structure-function studies of GPCR signaling in native contexts will inform future studies into the basic biology of this large and clinically important family of receptors.

A conserved protonation-induced switch can trigger "ionic-lock" formation in adrenergic receptors

The mechanism of signal transduction in G-protein-coupled receptors (GPCRs) is a crucial step in cell signaling. However, the molecular details of this process are still largely undetermined. Carrying out submicrosecond molecular dynamics simulations of beta-adrenergic receptors, we found that cooperation between a number of highly conserved residues is crucial to alter the equilibrium between the active state and the inactive state of diffusible ligand GPCRs. In particular, "ionic-lock" formation in beta-adrenergic receptors is directly correlated with the protonation state of a highly conserved aspartic acid residue [Asp(2.50)] even though the two sites are located more than 20 A away from each other. Internal polar residues, acting as local microswitches, cooperate to propagate the signal from Asp(2.50) to the G-protein interaction site at the helix III-helix VI interface. Evolutionarily conserved differences between opsin and non-opsin GPCRs in the surrounding of Asp(2.50) influence the acidity of this residue and can thus help in rationalizing the differences in constitutive activity of class A GPCRs.

Constitutive activity and inverse agonism at the α(₁a) and α(₁b) adrenergic receptor subtypes

The α(1b)-adrenergic receptor (AR) was, after rhodopsin, the first G protein-coupled receptor (GPCR) in which point mutations were shown to trigger constitutive (agonist-independent) activity. Constitutively activating mutations have been found in other AR subtypes as well as in several GPCRs. This chapter briefly summarizes the main findings on constitutively active mutants of the α(1a)- and α(1b)-AR subtypes and the methods used to predict activating mutations, to measure constitutive activity of Gq-coupled receptors and to investigate inverse agonism. In addition, it highlights the implications of studies on constitutively active AR mutants on elucidating the molecular mechanisms of receptor activation and drug action.

Follicle stimulating hormone receptor mutations and reproductive disorders

The follicle stimulating hormone receptor (FSHR) plays a critical role in reproductive function. In the males, FSH supports spermatogenesis, whereas in females, FSH is absolutely required for ovarian follicle growth. In females, inactivating mutations in the FSHR result in ovarian dysgenesis with amenorrhea and infertility. The few males reported with severe inactivating mutations exhibited varying spermatogenic defects, but not azoospermia. While these findings may potentially suggest that FSH action is not absolutely required for spermatogenesis, it cannot be ruled out that these individuals have some residual FSHR activity. Gain-of-function mutations in the FSHR cause spontaneous ovarian hyperstimulation syndrome in females due to the inappropriate stimulation of the mutant FSHR by human choriogonadotropin.

Lipid protein interactions couple protonation to conformation in a conserved cytosolic domain of G protein-coupled receptors

The visual photoreceptor rhodopsin is a prototypical class I (rhodopsin-like) G protein-coupled receptor. Photoisomerization of the covalently bound ligand 11-cis-retinal leads to restructuring of the cytosolic face of rhodopsin. The ensuing protonation of Glu-134 in the class-conserved D(E)RY motif at the C-terminal end of transmembrane helix-3 promotes the formation of the G protein-activating state. Using transmembrane segments derived from helix-3 of bovine rhodopsin, we show that lipid protein interactions play a key role in this cytosolic "proton switch." Infrared and fluorescence spectroscopic pK(a) determinations reveal that the D(E)RY motif is an autonomous functional module coupling side chain neutralization to conformation and helix positioning as evidenced by side chain to lipid headgroup Foerster resonance energy transfer. The free enthalpies of helix stabilization and hydrophobic burial of the neutral carboxyl shift the side chain pK(a) into the range typical of Glu-134 in photoactivated rhodopsin. The lipid-mediated coupling mechanism is independent of interhelical contacts allowing its conservation without interference with the diversity of ligand-specific interactions in class I G protein-coupled receptors.

Observation of "ionic lock" formation in molecular dynamics simulations of wild-type beta 1 and beta 2 adrenergic receptors

G protein coupled receptors (GPCRs) are a large family of integral membrane proteins involved in signal transduction pathways, making them appealing drug targets for a wide spectrum of diseases. The recently crystallized structures of two engineered adrenergic receptors have opened new avenues for the understanding of the molecular mechanisms of action of GPCRs. Taking the two crystal structures as a starting point, we carried out submicrosecond molecular dynamics simulations of wild-type beta(1) and beta(2) adrenergic receptors in a lipid bilayer under physiological conditions. These simulations give access to structural and dynamic properties of the receptors in pseudo in vivo conditions. For both systems the overall fold properties of the transmembrane region as well as the binding pocket remain close to the crystal structure of the engineered systems, thus suggesting that the ligand binding mode is not affected by the introduced modifications. Both simulations indicate the presence of one or two internal water molecules absent in both crystal structures and essential for the stabilization of the binding pocket at the interface between transmembrane helices III, IV, and V. The different interactions arising from the substitution of Tyr308 in beta(2)AR into Phe325 in beta(1)AR induce different conformations of the homologous Asn(6.55) inside the binding pockets of the two receptors, suggesting a possible origin of receptor specificity in agonist binding. The equilibrated structures of both receptors recover all of the previously suggested features of inactive GPCRs including formation of a salt bridge between the cytoplasmatic moieties of helices III and VI ("ionic lock") that is absent in the crystal structures.

Identification of two distinct inactive conformations of the beta2-adrenergic receptor reconciles structural and biochemical observations

Fully understanding the mechanisms of signaling proteins such as G protein-coupled receptors (GPCRs) will require the characterization of their conformational states and the pathways connecting those states. The recent crystal structures of the beta(2)- and beta(1)-adrenergic receptors in a nominally inactive state constituted a major advance toward this goal, but also raised new questions. Although earlier biochemical observations had suggested that these receptors possessed a set of contacts between helices 3 and 6, known as the ionic lock, which was believed to form a molecular switch for receptor activation, the crystal structures lacked these contacts. The unexpectedly broken ionic lock has raised questions about the true conformation(s) of the inactive state and the role of the ionic lock in receptor activation and signaling. To address these questions, we performed microsecond-timescale molecular dynamics simulations of the beta(2)-adrenergic receptor (beta(2)AR) in multiple wild-type and mutant forms. In wild-type simulations, the ionic lock formed reproducibly, bringing the intracellular ends of helices 3 and 6 together to adopt a conformation similar to that found in inactive rhodopsin. Our results suggest that inactive beta(2)AR exists in equilibrium between conformations with the lock formed and the lock broken, whether or not the cocrystallized ligand is present. These findings, along with the formation of several secondary structural elements in the beta(2)AR loops during our simulations, may provide a more comprehensive picture of the inactive state of the beta-adrenergic receptors, reconciling the crystal structures with biochemical studies.

Two arginine-glutamate ionic locks near the extracellular surface of FFAR1 gate receptor activation

Activation of a number of class A G protein-coupled receptors (GPCRs) is thought to involve two molecular switches, a rotamer toggle switch within the transmembrane domain and an ionic lock at the cytoplasmic surface of the receptor; however, the mechanism by which agonist binding changes these molecular interactions is not understood. Importantly, 80% of GPCRs including free fatty acid receptor 1 (FFAR1) lack the complement of amino acid residues implicated in either or both of these two switches; the mechanism of activation of these GPCRs is therefore less clear. By homology modeling, we identified two Glu residues (Glu-145 and Glu-172) in the second extracellular loop of FFAR1 that form putative interactions individually with two transmembrane Arg residues (Arg-183(5.39) and Arg-258(7.35)) to create two ionic locks. Molecular dynamics simulations showed that binding of agonists to FFAR1 leads to breakage of these Glu-Arg interactions. In mutagenesis experiments, breakage of these two putative interactions by substituting Ala for Glu-145 and Glu-172 caused constitutive receptor activation. Our results therefore reveal a molecular switch for receptor activation present on the extracellular surface of FFAR1 that is broken by agonist binding. Similar ionic locks between the transmembrane domains and the extracellular loops may constitute a mechanism common to other class A GPCRs also.

The third intracellular loop stabilizes the inactive state of the neuropeptide Y1 receptor

Constitutively active G-protein-coupled receptors (GPCRs) can signal even in the absence of ligand binding. Most Class I GPCRs are stabilized in the resting conformation by intramolecular interactions involving transmembrane domain (TM) 3 and TM6, particularly at loci 6.30 and 6.34 of TM6. Signaling by Gi/Go-coupled receptors such as the Neuropeptide Y1 receptor decreases already low basal metabolite levels. Thus, we examined constitutive activity using a biochemical assay mediated by a Gi/Gq chimeric protein and a more direct electrophysiological assay. Wild-type (WT-Y1) receptors express no measurable, agonist-independent activation, while mu-opioid receptors (MOR) and P2Y12 purinoceptors showed clear evidence of constitutive activation, especially in the electrophysiological assay. Neither point mutations at TM6 (T6.30A or N6.34A) nor substitution of the entire TM3 and TM6 regions from the MOR into the Y1 receptor increased basal WT-Y1 activation. By contrast, chimeric substitution of the third intracellular loop (ICL3) generated a constitutively active, Y1-ICL3-MOR chimera. Furthermore, the loss of stabilizing interactions from the native ICL3 enhanced the role of surrounding residues to permit basal receptor activation; because constitutive activity of the Y1-ICL3-MOR chimera was further increased by point mutation at locus 6.34, which did not alter WT-Y1 receptor activity. Our results indicate that the ICL3 stabilizes the Y1 receptor in the inactive state and confers structural properties critical for regulating Y receptor activation and signal transduction. These studies reveal the active participation of the ICL3 in the stabilization and activation of Class I GPCRs.

Constitutive activation of G protein-coupled receptors and diseases: insights into mechanisms of activation and therapeutics

The existence of constitutive activity for G protein-coupled receptors (GPCRs) was first described in 1980s. In 1991, the first naturally occurring constitutively active mutations in GPCRs that cause diseases were reported in rhodopsin. Since then, numerous constitutively active mutations that cause human diseases were reported in several additional receptors. More recently, loss of constitutive activity was postulated to also cause diseases. Animal models expressing some of these mutants confirmed the roles of these mutations in the pathogenesis of the diseases. Detailed functional studies of these naturally occurring mutations, combined with homology modeling using rhodopsin crystal structure as the template, lead to important insights into the mechanism of activation in the absence of crystal structure of GPCRs in active state. Search for inverse agonists on these receptors will be critical for correcting the diseases cause by activating mutations in GPCRs. Theoretically, these inverse agonists are better therapeutics than neutral antagonists in treating genetic diseases caused by constitutively activating mutations in GPCRs.

A novel adrenocorticotropin receptor mutation alters its structure and function, causing familial glucocorticoid deficiency

CONTEXT

Familial glucocorticoid deficiency (FGD) is an autosomal recessive disorder characterized by unresponsiveness to ACTH. In this study, two mutations of the ACTH receptor (MC2R) gene are reported in this FGD clinical case.

OBJECTIVE

The objective of the study was to characterize a novel MC2R gene mutation in a compound heterozygous patient with FGD phenotype.

DESIGN

This was a clinical case description, biochemical, molecular, and bioinformatics analysis to describe a novel MC2R gene mutation.

PATIENTS

The subject of the study was a male diagnosed with primary adrenal insufficiency. The family history showed nonconsanguineous healthy parents, three healthy siblings, and one brother affected with FGD.

MAIN OUTCOME MEASURES

The mutant MC2R-Ala126Ser showed significantly lower activity when it was stimulated with ACTH-(1-24) than did cells transfected with wild-type MC2R.

RESULTS

The molecular studies demonstrated the presence of an adenine heterozygous insertion (InsA1347) in the MC2R gene (G217fs) in the patient. This insertion was due to a frame shift mutation in one allele and a premature stop codon codifying an aberrant receptor of 247 residues (27.2 kDa). We also found a novel heterozygous mutation alanine 126 by serine. Molecular dynamic simulations showed that serine 126 side chain fluctuates forming a noncanonical intrahelical hydrogen bond in the transmembrane helix 3 of the mutated receptor. This produces a structural rearrangement of the MC2R internal cavities that may affect the ligand recognition and signal transduction throughout the G protein.

CONCLUSIONS

We propose a molecular explanation for the reduced activity exhibited by the MC2R alanine 126 by serine mutant.

Atomistic insights into rhodopsin activation from a dynamic model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.