GPCR activation: a mutagenic spotlight on crystal structures

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.

Characterizing conformational states in GPCR structures using machine learning

G protein-coupled receptors (GPCRs) play a pivotal role in signal transduction and represent attractive targets for drug development. Recent advances in structural biology have provided insights into GPCR conformational states, which are critical for understanding their signaling pathways and facilitating structure-based drug discovery. In this study, we introduce a machine learning approach for conformational state annotation of GPCRs. We represent GPCR conformations as high-dimensional feature vectors, incorporating information about amino acid residue pairs involved in the activation pathway. Using a dataset of GPCR conformations in inactive and active states obtained through molecular dynamics simulations, we trained machine learning models to distinguish between inactive-like and active-like conformations. The developed model provides interpretable predictions and can be used for the large-scale analysis of molecular dynamics trajectories of GPCRs.

Protons taken hostage: Dynamic H-bond networks of the pH-sensing GPR68

Proton-sensing G Protein Coupled Receptors (GPCRs) sense changes in the extracellular pH to effect cell signaling for cellular homeostasis. They tend to be overexpressed in solid tumors associated with acidic extracellular pH, and are of direct interest as drug targets. How proton-sensing GPCRs sense extracellular acidification and activate upon protonation change is important to understand, because it may guide the design of therapeutics. Lack of publicly available experimental structures make it challenging to discriminate between conflicting mechanisms proposed for proton-binding, as main roles have been assigned to either an extracellular histidine cluster or to an internal carboxylic triad. Here we present a protocol to derive and evaluate structural models of the proton-sensing GPR68. This approach integrates state-of-the-art homology modeling with microsecond-timescale atomistic simulations, and with a detailed assessment of the compatibility of the structural models with known structural features of class A GPCRs. To decipher structural elements of potential interest for protonation-coupled conformational changes of GPR68, we used the best-compatible model as a starting point for independent atomistic simulations of GPR68 with different protonation states, and graph computations to characterize the response of GPR68 to changes in protonation. We found that GPR68 hosts an extended hydrogen-bond network that inter-connects the extracellular histidine cluster to the internal carboxylic triad, and which can even reach groups at the cytoplasmic G-protein binding site. Taken together, results suggest that GPR68 relies on dynamic, hydrogen-bond networks to inter-connect extracellular and internal proton-binding sites, and to elicit conformational changes at the cytoplasmic G-protein binding site.

Role of Conserved Tyrosine Lid Residues in the Activation of the M2 Muscarinic Acetylcholine Receptor

The development of subtype selective small molecule drugs for the muscarinic acetylcholine receptor (mAChR) family has been challenging. The design of more selective ligands can be improved by understanding the structure and function of key amino acid residues that line ligand binding sites. Here we study the role of three conserved key tyrosine residues [Y1043.33, Y4036.51, and Y4267.39 (Ballesteros and Weinstein numbers in superscript)] at the human M2 mAChR, located at the interface between the orthosteric and allosteric binding sites of the receptor. We specifically focused on the role of the three tyrosine hydroxyl groups in the transition between the inactive and active conformations of the receptor by making phenylalanine point mutants. Single-point mutation at either of the three positions was sufficient to reduce the affinity of agonists by ∼100-fold for the M2 mAChR, whereas the affinity of antagonists remained largely unaffected. In contrast, neither of the mutations affected the efficacy of orthosteric agonists. When mutations were combined into double and triple M2 mAChR mutants, the affinity of antagonists was reduced by more than 100-fold compared with the wild-type M2 receptor. In contrast, the affinity of allosteric modulators, either negative or positive, was retained at all single and multiple mutations, but the degree of allosteric effect exerted on the endogenous ligand acetylcholine was affected at all mutants containing Y4267.39F. These findings will provide insights to consider when designing future mAChR ligands. SIGNIFICANCE STATEMENT: Structural studies demonstrated that three tyrosine residues between the orthosteric and allosteric sites of the M2 muscarinic acetylcholine receptor (mAChR) had different hydrogen bonding networks in the inactive and active conformations. The role of hydroxyl groups of the tyrosine residues on orthosteric and allosteric ligand pharmacology was unknown. We found that hydroxyl groups of the tyrosine residues differentially affected the molecular pharmacology of orthosteric and allosteric ligands. These results provide insights to consider when designing future mAChR ligands.

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.

Impact of cancer-associated mutations in CC chemokine receptor 2 on receptor function and antagonism

CC chemokine receptor 2 (CCR2), a G protein-coupled receptor, plays a role in many cancer-related processes such as metastasis formation and immunosuppression. Since ∼ 20 % of human cancers contain mutations in G protein-coupled receptors, ten cancer-associated CCR2 mutants obtained from the Genome Data Commons were investigated for their effect on receptor functionality and antagonist binding. Mutations were selected based on either their vicinity to CCR2's orthosteric or allosteric binding sites or their presence in conserved amino acid motifs. One of the mutant receptors, namely S101P^2.63 with a mutation near the orthosteric binding site, did not express on the cell surface. All other studied mutants showed a decrease in or a lack of G protein activation in response to the main endogenous CCR2 ligand CCL2, but no change in potency was observed. Furthermore, INCB3344 and LUF7482 were chosen as representative orthosteric and allosteric antagonists, respectively. No change in potency was observed in a functional assay, but mutations located at F116^3.28 impacted orthosteric antagonist binding significantly, while allosteric antagonist binding was abolished for L134Q^3.46 and D137N^3.49 mutants. As CC chemokine receptor 2 is an attractive drug target in cancer, the negative effect of these mutations on receptor functionality and drugability should be considered in the drug discovery process.

New insights into the structure and function of class B1 GPCRs

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors. Class B1 GPCRs constitute a subfamily of 15 receptors that characteristically contain large extracellular domains (ECDs) and respond to long polypeptide hormones. Class B1 GPCRs are critical regulators of homeostasis, and as such, many are important drug targets. While most transmembrane proteins, including GPCRs, are recalcitrant to crystallization, recent advances in electron cryo-microscopy (cryo-EM) have facilitated a rapid expansion of the structural understanding of membrane proteins. As a testament to this success, structures for all the class B1 receptors bound to G proteins have been determined by cryo-EM in the past five years. Further advances in cryo-EM have uncovered dynamics of these receptors, ligands, and signalling partners. Here, we examine the recent structural underpinnings of the class B1 GPCRs with an emphasis on structure-function relationships.

Allosteric Effect of Nanobody Binding on Ligand-Specific Active States of the β2 Adrenergic Receptor

Nanobody binding stabilizes G-protein-coupled receptors (GPCR) in a fully active state and modulates their affinity for bound ligands. However, the atomic-level basis for this allosteric regulation remains elusive. Here, we investigate the conformational changes induced by the binding of a nanobody (Nb80) on the active-like β2 adrenergic receptor (β2AR) via enhanced sampling molecular dynamics simulations. Dimensionality reduction analysis shows that Nb80 stabilizes structural features of the β2AR with an ∼14 Å outward movement of transmembrane helix 6 and a close proximity of transmembrane (TM) helices 5 and 7, and favors the fully active-like conformation of the receptor, independent of ligand binding, in contrast to the conditions under which no intracellular binding partner is bound, in which case the receptor is only stabilized in an intermediate-active state. This activation is supported by the residues located at hotspots located on TMs 5, 6, and 7, as shown by supervised machine learning methods. Besides, ligand-specific subtle differences in the conformations assumed by intracellular loop 2 and extracellular loop 2 are captured from the trajectories of various ligand-bound receptors in the presence of Nb80. Dynamic network analysis further reveals that Nb80 binding triggers tighter and stronger local communication networks between the Nb80 and the ligand-binding sites, primarily involving residues around ICL2 and the intracellular end of TM3, TM5, TM6, as well as ECL2, ECL3, and the extracellular ends of TM6 and TM7. In particular, we identify unique allosteric signal transmission mechanisms between the Nb80-binding site and the extracellular domains in conformations modulated by a full agonist, BI167107, and a G-protein-biased partial agonist, salmeterol, involving mainly TM1 and TM2, and TM5, respectively. Altogether, our results provide insights into the effect of intracellular binding partners on the GPCR activation mechanism, which should be taken into account in structure-based drug discovery.

Encoding mu-opioid receptor biased agonism with interaction fingerprints

Opioids are potent painkillers, however, their therapeutic use requires close medical monitoring to diminish the risk of severe adverse effects. The G-protein biased agonists of the μ-opioid receptor (MOR) have shown safer therapeutic profiles than non-biased ligands. In this work, we performed extensive all-atom molecular dynamics simulations of two markedly biased ligands and a balanced reference molecule. From those simulations, we identified a protein-ligand interaction fingerprint that characterizes biased ligands. Then, we built and virtually screened a database containing 68,740 ligands with proven or potential GPCR agonistic activity. Exemplary molecules that fulfill the interacting pattern for biased agonism are showcased, illustrating the usefulness of this work for the search of biased MOR ligands and how this contributes to the understanding of MOR biased signaling.

Functionally selective and biased agonists of muscarinic receptors

Disruption of cholinergic signalling via muscarinic receptors is associated with various pathologies, like Alzheimer's disease or schizophrenia. Selective muscarinic agonists possess therapeutic potential in the treatment of diabetes, pain or Sjögren's syndrome. The orthosteric binding site of all subtypes of the muscarinic receptor is structurally identical, making the development of affinity-based selective agonists virtually impossible. Some agonists, however, are functionally selective; they activate only a subset of receptors or signalling pathways. Others may stabilise specific conformations of the receptor leading to non-uniform modulation of individual signalling pathways (biased agonists). Functionally selective and biased agonists represent a promising approach for selective activation of individual subtypes of muscarinic receptors. In this work we review chemical structures, receptor binding and agonist-specific conformations of currently known functionally selective and biased muscarinic agonists in the context of their intricate intracellular signalling. Further, we take a perspective on the possible use of biased agonists for tissue and organ-specific activation of muscarinic receptors.

Allosteric Modulation of GPCRs of Class A by Cholesterol

G-protein coupled receptors (GPCRs) are membrane proteins that convey extracellular signals to the cellular milieu. They represent a target for more than 30% of currently marketed drugs. Here we review the effects of membrane cholesterol on the function of GPCRs of Class A. We review both the specific effects of cholesterol mediated via its direct high-affinity binding to the receptor and non-specific effects mediated by cholesterol-induced changes in the properties of the membrane. Cholesterol binds to many GPCRs at both canonical and non-canonical binding sites. It allosterically affects ligand binding to and activation of GPCRs. Additionally, it changes the oligomerization state of GPCRs. In this review, we consider a perspective of the potential for the development of new therapies that are targeted at manipulating the level of membrane cholesterol or modulating cholesterol binding sites on to GPCRs.

The Conformational Equilibrium of the Neuropeptide Y2 Receptor in Bilayer Membranes

Dynamic structural transitions within the seven-transmembrane bundle represent the mechanism by which G-protein-coupled receptors convert an extracellular chemical signal into an intracellular biological function. Here, the conformational dynamics of the neuropeptide Y receptor type 2 (Y2R) during activation was investigated. The apo, full agonist-, and arrestin-bound states of Y2R were prepared by cell-free expression, functional refolding, and reconstitution into lipid membranes. To study conformational transitions between these states, all six tryptophans of Y2R were 13 C-labeled. NMR-signal assignment was achieved by dynamic-nuclear-polarization enhancement and the individual functional states of the receptor were characterized by monitoring 13 C NMR chemical shifts. Activation of Y2R is mediated by molecular switches involving the toggle switch residue Trp2816.48 of the highly conserved SWLP motif and Trp3277.55 adjacent to the NPxxY motif. Furthermore, a conformationally preserved "cysteine lock"-Trp11623.50 was identified.

Structural and functional characterization of G protein-coupled receptors with deep mutational scanning

The >800 human G protein–coupled receptors (GPCRs) are responsible for transducing diverse chemical stimuli to alter cell state- and are the largest class of drug targets. Their myriad structural conformations and various modes of signaling make it challenging to understand their structure and function. Here, we developed a platform to characterize large libraries of GPCR variants in human cell lines with a barcoded transcriptional reporter of G protein signal transduction. We tested 7800 of 7828 possible single amino acid substitutions to the beta-2 adrenergic receptor (β₂AR) at four concentrations of the agonist isoproterenol. We identified residues specifically important for β₂AR signaling, mutations in the human population that are potentially loss of function, and residues that modulate basal activity. Using unsupervised learning, we identify residues critical for signaling, including all major structural motifs and molecular interfaces. We also find a previously uncharacterized structural latch spanning the first two extracellular loops that is highly conserved across Class A GPCRs and is conformationally rigid in both the inactive and active states of the receptor. More broadly, by linking deep mutational scanning with engineered transcriptional reporters, we establish a generalizable method for exploring pharmacogenomics, structure and function across broad classes of drug receptors.

Atomic-Level Characterization of the Methadone-Stabilized Active Conformation of µ-Opioid Receptor

Methadone is a synthetic opioid agonist with notoriously unique properties, such as lower abuse liability and induced relief of withdrawal symptoms and drug cravings, despite acting on the same opioid receptors triggered by classic opioids-in particular the µ-opioid receptor (MOR). Its distinct pharmacologic properties, which have recently been attributed to the preferential activation of β-arrestin over G proteins, make methadone a standard-of-care maintenance medication for opioid addiction. Although a recent biophysical study suggests that methadone stabilizes different MOR active conformations from those stabilized by classic opioid drugs or G protein-biased agonists, how this drug modulates the conformational equilibrium of MOR and what specific active conformation of the receptor it stabilizes are unknown. Here, we report the results of submillisecond adaptive sampling molecular dynamics simulations of a predicted methadone-bound MOR complex and compare them with analogous data obtained for the classic opioid morphine and the G protein-biased ligand TRV130. The model, which is supported by existing experimental data, is analyzed using Markov state models and transfer entropy analysis to provide testable hypotheses of methadone-specific conformational dynamics and activation kinetics of MOR. SIGNIFICANCE STATEMENT: Opioid addiction has reached epidemic proportions in both industrialized and developing countries. Although methadone maintenance treatment represents an effective therapeutic approach for opioid addiction, it is not as widely used as needed. In this study, we contribute an atomic-level understanding of how methadone exerts its unique function in pursuit of more accessible treatments for opioid addiction. In particular, we present details of a methadone-specific active conformation of the µ-opioid receptor that has thus far eluded experimental structural characterization.

Elucidation of Mechanism for Ligand Efficacy at Leukotriene B4 Receptor 2 (BLT2)

G protein-coupled receptors (GPCRs) have always been important drug targets in the pharmaceutical industry. One major question for the current GPCR drug discovery is how drugs have distinct efficacies at the same GPCR target. Related to this question, we studied how different ligands can have disparate efficacies at Leukotriene B₄ receptor (BLT2). By using molecular modeling studies, we predicted that Tyr271^6.51 located at TM6 of BLT2 performs as a key trigger for its activation and verified the prediction by site-directed mutagenesis, chemotactic motility studies, which included a chemical derivative of agonist CAY10583. We further identified Asn275^6.55 located at TM6 as a weak activation trigger in BLT2 and performed double mutation studies to confirm our computational results. Our results provide strong evidence for the exact mechanism of ligand efficacy at BLT2.

Agonist-Specific Conformations of the M2 Muscarinic Acetylcholine Receptor Assessed by Molecular Dynamics

Binding of muscarinic ligands, both antagonists and agonists, and their effects on the conformation of the M₂ acetylcholine receptor were modeled in silico and compared to experimental data. After docking of antagonists to the M₂ receptor in an inactive conformation (3UON, 5ZK3, 5ZKB, or 5ZKB) and agonists in an active conformation (4MQS), 100 ns of conventional molecular dynamics (MD) followed by 500 ns of accelerated MD was run. Conventional MD revealed ligand-specific interactions with the receptor. Antagonists stabilized the receptor in an inactive conformation during accelerated MD. The receptor in complex with various agonists attained different conformations specific to individual agonists. The magnitude of the TM6 movement correlated with agonist efficacy at the non-preferential G_s pathway. The shape of the intracellular opening where the receptor interacts with a G-protein was different for the classical agonist carbachol, super-agonist iperoxo, and G_i/o-biased partial agonists JR-6 and JR-7, being compatible with experimentally observed agonist bias at the G-protein level. Moreover, a wash-resistant binding of the unique agonist xanomeline associated with interactions with membrane lipids was formed during accelerated MD. Thus, accelerated MD is suitable for modeling of ligand-specific receptor binding and receptor conformations that is essential for the design of experiments aimed at identification of the secondary binding sites and understanding molecular mechanisms underlying receptor activation.

Common activation mechanism of class A GPCRs

Class A G-protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology. Understanding receptor activation mechanism is critical for discovering novel therapeutics since about one-third of all marketed drugs target members of this family. GPCR activation is an allosteric process that couples agonist binding to G-protein recruitment, with the hallmark outward movement of transmembrane helix 6 (TM6). However, what leads to TM6 movement and the key residue level changes of this movement remain less well understood. Here, we report a framework to quantify conformational changes. By analyzing the conformational changes in 234 structures from 45 class A GPCRs, we discovered a common GPCR activation pathway comprising of 34 residue pairs and 35 residues. The pathway unifies previous findings into a common activation mechanism and strings together the scattered key motifs such as CWxP, DRY, Na+ pocket, NPxxY and PIF, thereby directly linking the bottom of ligand-binding pocket with G-protein coupling region. Site-directed mutagenesis experiments support this proposition and reveal that rational mutations of residues in this pathway can be used to obtain receptors that are constitutively active or inactive. The common activation pathway provides the mechanistic interpretation of constitutively activating, inactivating and disease mutations. As a module responsible for activation, the common pathway allows for decoupling of the evolution of the ligand binding site and G-protein-binding region. Such an architecture might have facilitated GPCRs to emerge as a highly successful family of proteins for signal transduction in nature.

Crystal structure of dopamine receptor D4 bound to the subtype selective ligand, L745870

Multiple subtypes of dopamine receptors within the GPCR superfamily regulate neurological processes through various downstream signaling pathways. A crucial question about the dopamine receptor family is what structural features determine the subtype-selectivity of potential drugs. Here, we report the 3.5-angstrom crystal structure of mouse dopamine receptor D4 (DRD4) complexed with a subtype-selective antagonist, L745870. Our structure reveals a secondary binding pocket extended from the orthosteric ligand-binding pocket to a DRD4-specific crevice located between transmembrane helices 2 and 3. Additional mutagenesis studies suggest that the antagonist L745870 prevents DRD4 activation by blocking the relative movement between transmembrane helices 2 and 3. These results expand our knowledge of the molecular basis for the physiological functions of DRD4 and assist new drug design.

Basal Histamine H4 Receptor Activation: Agonist Mimicry by the Diphenylalanine Motif

Histamine H4 receptor (H4 R) orthologues are G-protein-coupled receptors (GPCRs) that exhibit species-dependent basal activity. In contrast to the basally inactive mouse H4 R (mH4 R), human H4 R (hH4 R) shows a high degree of basal activity. We have performed long-timescale molecular dynamics simulations and rigidity analyses on wild-type hH4 R, the experimentally characterized hH4 R variants S179M, F169V, F169V+S179M, F168A, and on mH4 R to investigate the molecular nature of the differential basal activity. H4 R variant-dependent differences between essential motifs of GPCR activation and structural stabilities correlate with experimentally determined basal activities and provide a molecular explanation for the differences in basal activation. Strikingly, during the MD simulations, F16945.55 dips into the orthosteric binding pocket only in the case of hH4 R, thus adopting the role of an agonist and contributing to the stabilization of the active state. The results shed new light on the molecular mechanism of basal H4 R activation that are of importance for other GPCRs.

Functional characterization in vitro of twelve naturally occurring variants of the human pancreatic polypeptide receptor NPY4R

Obesity has become a global health problem and therefore understanding of the mechanisms regulating hunger and satiety is of utmost importance for the development of new treatment strategies. The Y4 receptor, encoded by the NPY4R gene, and its ligand pancreatic polypeptide (PP) have been reported to mediate a satiety signal. Multiple genetic studies have reported an association between NPY4R copy number and body weight. The gene also displays several SNP variants, many of which lead to amino acid differences, making it interesting to study. We have investigated the functional properties of 12 naturally occurring amino acid sequence variants of the Y4 and interpret the results in relation to sequence conservation and our structural model of the human Y4 receptor protein. Three receptor variants, Cys201^ECL2Tyr, Val271^6.41Leu and Asn318^7.49Asp, were found to completely lose functional response, measured as inositol phosphate turnover, while retaining membrane expression. They display high sequence conservation and have important roles in the receptor structure. For two receptor variants the potency of PP was significantly decreased, Cys34^NTSer (EC50 = 2.9 nM, p < .001) and Val135^3.46Met (EC50 = 3.0 nM, p < .01), compared to wild-type Y4 (EC50 = 0.68 nM). Cys34 forms a disulphide bond with Cys298, linking the N-terminal part to ECL3. The Val135^3.46Met variant has an amino acid replacement located in the TM3 helix, one helix turn above the highly conserved ERH motif. This position has influence on the network of residues involved in receptor activation and subsequent inactivation. Sequence conservation and the structural model are consistent with these results. The remaining seven positions had no significant effect on the receptor's functional response compared to wild-type Y4. These positions display more variation during evolution. Understanding of the interactions between the Y4 receptor and its native PP agonist and the effects of amino acid variation on its functional response will hopefully lead to future therapeutic possibilities.

Acetylcholine and cholinergic receptors

This review provides a distillate of the advances in knowledge about the neurotransmitter functions of acetylcholine over the 50-year period between 1967 and 2017, together with incremental information about the cognate nicotinic and muscarinic acetylcholine receptors, and some brief comments on possible advances in the near future. The text is supplemented by a timelines figure indicating the dates of some key advances in knowledge about acetylcholine receptors and a box-figure providing a snapshot of selected papers about acetylcholine published in the year 1967.

Fentanyl Family at the Mu-Opioid Receptor: Uniform Assessment of Binding and Computational Analysis

Interactions of 21 fentanyl derivatives with μ-opioid receptor (μOR) were studied using experimental and theoretical methods. Their binding to μOR was assessed with radioligand competitive binding assay. A uniform set of binding affinity data contains values for two novel and one previously uncharacterized derivative. The data confirms trends known so far and thanks to their uniformity, they facilitate further comparisons. In order to provide structural hypotheses explaining the experimental affinities, the complexes of the studied derivatives with μOR were modeled and subject to molecular dynamics simulations. Five common General Features (GFs) of fentanyls' binding modes stemmed from these simulations. They include: GF1) the ionic interaction between D147 and the ligands' piperidine NH⁺ moiety; GF2) the N-chain orientation towards the μOR interior; GF3) the other pole of ligands is directed towards the receptor outlet; GF4) the aromatic anilide ring penetrates the subpocket formed by TM3, TM4, ECL1 and ECL2; GF5) the 4-axial substituent (if present) is directed towards W318. Except for the ionic interaction with D147, the majority of fentanyl-μOR contacts is hydrophobic. Interestingly, it was possible to find nonlinear relationships between the binding affinity and the volume of the N-chain and/or anilide's aromatic ring. This kind of relationships is consistent with the apolar character of interactions involved in ligand⁻receptor binding. The affinity reaches the optimum for medium size while it decreases for both large and small substituents. Additionally, a linear correlation between the volumes and the average dihedral angles of W293 and W133 was revealed by the molecular dynamics study. This seems particularly important, as the W293 residue is involved in the activation processes. Further, the Y326 (OH) and D147 (Cγ) distance found in the simulations also depends on the ligands' size. In contrast, neither RMSF measures nor D114/Y336 hydrations show significant structure-based correlations. They also do not differentiate studied fentanyl derivatives. Eventually, none of 14 popular scoring functions yielded a significant correlation between the predicted and observed affinity data (R < 0.30, n = 28).

'Crystal' Clear? Lysophospholipid Receptor Structure Insights and Controversies

Lysophospholipids (LPLs), particularly sphingosine 1-phosphate (S1P) and lysophosphatidic acid (LPA), are bioactive lipid modulators of cellular homeostasis and pathology. The discovery and characterization of five S1P- and six LPA-specific G protein-coupled receptors (GPCRs), S1P1-5 and LPA1-6, have expanded their known involvement in all mammalian physiological systems. Resolution of the S1P1, LPA1, and LPA6 crystal structures has fueled the growing interest in these receptors and their ligands as targets for pharmacological manipulation. In this review, we have attempted to provide an integrated overview of the three crystallized LPL GPCRs with biochemical and physiological structure-function data. Finally, we provide a novel discussion of how chaperones for LPLs may be considered when extrapolating crystallographic and computational data toward understanding actual biological interactions and phenotypes.

Distinct but overlapping binding sites of agonist and antagonist at the relaxin family peptide 3 (RXFP3) receptor

The relaxin-3 neuropeptide activates the relaxin family peptide 3 (RXFP3) receptor to modulate stress, appetite, and cognition. RXFP3 shows promise as a target for treating neurological disorders, but realization of its clinical potential requires development of smaller RXFP3-specific drugs that can penetrate the blood-brain barrier. Designing such drugs is challenging and requires structural knowledge of agonist- and antagonist-binding modes. Here, we used structure-activity data for relaxin-3 and a peptide RXFP3 antagonist termed R3 B1-22R to guide receptor mutagenesis and develop models of their binding modes. RXFP3 residues were alanine-substituted individually and in combination and tested in cell-based binding and functional assays to refine models of agonist and antagonist binding to active- and inactive-state homology models of RXFP3, respectively. These models suggested that both agonists and antagonists interact with RXFP3 via similar residues in their B-chain central helix. The models further suggested that the B-chain Trp²⁷ inserts into the binding pocket of RXFP3 and interacts with Trp¹³⁸ and Lys²⁷¹, the latter through a salt bridge with the C-terminal carboxyl group of Trp²⁷ in relaxin-3. R3 B1-22R, which does not contain Trp²⁷, used a non-native Arg²³ residue to form cation-π and salt-bridge interactions with Trp¹³⁸ and Glu¹⁴¹ in RXFP3, explaining a key contribution of Arg²³ to affinity. Overall, relaxin-3 and R3 B1-22R appear to share similar binding residues but may differ in binding modes, leading to active and inactive RXFP3 conformational states, respectively. These mechanistic insights may assist structure-based drug design of smaller relaxin-3 mimetics to manage neurological disorders.

Full rescue of an inactive olfactory receptor mutant by elimination of an allosteric ligand-gating site

Ligand-gating has recently been proposed as a novel mechanism to regulate olfactory receptor sensitivity. TAAR13c, the zebrafish olfactory receptor activated by the death-associated odor cadaverine, appears to possess an allosteric binding site for cadaverine, which was assumed to block progress of the ligand towards the internal orthosteric binding-and-activation site. Here we have challenged the suggested gating mechanism by modeling the entry tunnel for the ligand as well as the ligand path inside the receptor. We report an entry tunnel, whose opening is blocked by occupation of the external binding site by cadaverine, confirming the hypothesized gating mechanism. A multistep docking algorithm suggested a plausible path for cadaverine from the allosteric to the orthosteric binding-and-activation site. Furthermore we have combined a gain-of-function gating site mutation and a loss-of-function internal binding site mutation in one recombinant receptor. This receptor had almost wildtype ligand affinities, consistent with modeling results that showed localized effects for each mutation. A novel mutation of the suggested gating site resulted in increased receptor ligand affinity. In summary both the experimental and the modeling results provide further evidence for the proposed gating mechanism, which surprisingly exhibits pronounced similarity to processes described for some metabotropic neurotransmitter receptors.

Exploring the mechanism of F282L mutation-caused constitutive activity of GPCR by a computational study

G-protein-coupled receptors (GPCRs) are important drug targets and generally activated by ligands. However, some experiments found that GPCRs also give rise to constitutive activity through some mutations (viz., CAM), which are usually associated with different kinds of diseases. However, the mechanisms of CAMs and their roles in interactions with drug-ligands are unclear in experiments. Herein, we used microsecond molecular dynamics simulations to study the effect of one important F282L mutation on β2AR in order to address the questions above. With the aid of principle component and correlation analysis, our results revealed that the F282L mutation could increase the instability of the overall structure, increase the dramatic fluctuations of NPxxY and extracellular loops, and decrease restraint of the helices through weakening interhelical H-bonding and correlations between residues, which could partly contribute to the constitutive activity reported by the experiments. The observations from the protein structure network (PSN) analysis indicate that the mutant exhibits less information flow than the wild β2AR and weakens the role of TM5 and TM6 in the signal transmission, but it enhances the impact of TM3 on the orthosteric pathway and TM4 on the allosteric one. In addition, the results from the virtual screening reveal that the mutant prefers to select agonists rather than antagonists, similar to the active state but opposite of the inactive state, further confirming that the F282L mutation advances the activation of β2AR. Our observations provide valuable information for understanding the mechanism of the mutation-caused constitutive activity of GPCR and related drug-design.

Structural Connection between Activation Microswitch and Allosteric Sodium Site in GPCR Signaling

Sodium ions are endogenous allosteric modulators of many G-protein-coupled receptors (GPCRs). Mutation of key residues in the sodium binding motif causes a striking effect on G-protein signaling. We report the crystal structures of agonist complexes for two variants in the first sodium coordination shell of the human A_2A adenosine receptor, D52^2.50N and S91^3.39A. Both structures present an overall active-like conformation; however, the variants show key changes in the activation motif NPxxY. Changes in the hydrogen bonding network in this microswitch suggest a possible mechanism for modified G-protein signaling and enhanced thermal stability. These structures, signaling data, and thermal stability analysis with a panel of pharmacological ligands provide a basis for understanding the role of the sodium-coordinating residues on stability and G-protein signaling. Utilizing the D^2.50N variant is a promising method for stabilizing class A GPCRs to accelerate structural efforts and drug discovery.

Apelinergic System Structure and Function

Apelin and apela (ELABELA/ELA/Toddler) are two peptide ligands for a class A G-protein-coupled receptor named the apelin receptor (AR/APJ/APLNR). Ligand-AR interactions have been implicated in regulation of the adipoinsular axis, cardiovascular system, and central nervous system alongside pathological processes. Each ligand may be processed into a variety of bioactive isoforms endogenously, with apelin ranging from 13 to 55 amino acids and apela from 11 to 32, typically being cleaved C-terminal to dibasic proprotein convertase cleavage sites. The C-terminal region of the respective precursor protein is retained and is responsible for receptor binding and subsequent activation. Interestingly, both apelin and apela exhibit isoform-dependent variability in potency and efficacy under various physiological and pathological conditions, but most studies focus on a single isoform. Biophysical behavior and structural properties of apelin and apela isoforms show strong correlations with functional studies, with key motifs now well determined for apelin. Unlike its ligands, the AR has been relatively difficult to characterize by biophysical techniques, with most characterization to date being focused on effects of mutagenesis. This situation may improve following a recently reported AR crystal structure, but there are still barriers to overcome in terms of comprehensive biophysical study. In this review, we summarize the three components of the apelinergic system in terms of structure-function correlation, with a particular focus on isoform-dependent properties, underlining the potential for regulation of the system through multiple endogenous ligands and isoforms, isoform-dependent pharmacological properties, and biological membrane-mediated receptor interaction. © 2018 American Physiological Society. Compr Physiol 8:407-450, 2018.

Pathophysiological consequences of receptor mistraffic: Tales from the platelet P2Y12 receptor

Genetic variations in G protein-coupled receptor (GPCR) genes can disrupt receptor function in a wide variety of human genetic diseases, including platelet bleeding disorders. Platelets are critical for haemostasis with inappropriate platelet activation leading to the development of arterial thrombosis, which can result in heart attack and stroke whilst decreased platelet activity is associated with an increased risk of bleeding. GPCRs expressed on the surface of platelets play key roles in regulating platelet activity and therefore function. Receptors include purinergic receptors (P2Y₁ and P2Y₁₂), proteinase-activated receptor (PAR1 and PAR4) and thromboxane receptors (TPα), among others. Pharmacological blockade of these receptors forms a powerful therapeutic tool in the treatment and prevention of arterial thrombosis. With the advance of genomic technologies, there has been a substantial increase in the identification of naturally occurring rare and common GPCR variants. These variants include single-nucleotide polymorphisms (SNPs) and insertion or deletions that have the potential to alter GPCR expression or function. A number of defects in platelet GPCRs that disrupt receptor function have now been characterized in patients with mild bleeding disorders. This review will focus on rare, function-disrupting variants of platelet GPCRs with particular emphasis upon mutations in the P2Y₁₂ receptor gene that affect receptor traffic to modulate platelet function. Further this review will outline how the identification and characterization of function-disrupting GPCR mutations provides an essential link in translating our detailed understanding of receptor traffic and function in cell line studies into relevant human biological systems.

Structural insight to mutation effects uncover a common allosteric site in class C GPCRs

Motivation

Class C G protein-coupled receptors (GPCRs) regulate important physiological functions and allosteric modulators binding to the transmembrane domain constitute an attractive and, due to a lack of structural insight, a virtually unexplored potential for therapeutics and the food industry. Combining pharmacological site-directed mutagenesis data with the recent class C GPCR experimental structures will provide a foundation for rational design of new therapeutics.

Results

We uncover one common site for both positive and negative modulators with different amino acid layouts that can be utilized to obtain selectivity. Additionally, we show a large potential for structure-based modulator design, especially for four orphan receptors with high similarity to the crystal structures.

Availability and Implementation

All collated mutagenesis data is available in the GPCRdb mutation browser at http://gpcrdb.org/mutations/ and can be analyzed online or downloaded in excel format.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Mutagenesis Analysis Reveals Distinct Amino Acids of the Human Serotonin 5-HT2C Receptor Underlying the Pharmacology of Distinct Ligands

While exploring the structure-activity relationship of 4-phenyl-2-dimethylaminotetralins (PATs) at serotonin 5-HT_2C receptors, we discovered that relatively minor modification of PAT chemistry impacts function at 5-HT_2C receptors. In HEK293 cells expressing human 5-HT_2C-INI receptors, for example, (-)-trans-3'-Br-PAT and (-)-trans-3'-Cl-PAT are agonists regarding Gα_q-inositol phosphate signaling, whereas (-)-trans-3'-CF₃-PAT is an inverse agonist. To investigate the ligand-receptor interactions that govern this change in function, we performed site-directed mutagenesis of 14 amino acids of the 5-HT_2C receptor based on molecular modeling and reported G protein-coupled receptor crystal structures, followed by molecular pharmacology studies. We found that S3.36, T3.37, and F5.47 in the orthosteric binding pocket are critical for affinity (K_i) of all PATs tested, we also found that F6.44, M6.47, C7.45, and S7.46 are primarily involved in regulating EC/IC₅₀ functional potencies of PATs. We discovered that when residue S5.43, N6.55, or both are mutated to alanine, (-)-trans-3'-CF₃-PAT switches from inverse agonist to agonist function, and when N6.55 is mutated to leucine, (-)-trans-3'-Br-PAT switches from agonist to inverse agonist function. Notably, most point-mutations that affected PAT pharmacology did not significantly alter affinity (K_D) of the antagonist radioligand [³H]mesulergine, but every mutation tested negatively impacted serotonin binding. Also, amino acid mutations differentially affected the pharmacology of other commercially available 5-HT_2C ligands tested. Collectively, the data show that functional outcomes shared by different ligands are mediated by different amino acids and that some 5-HT_2C receptor residues important for pharmacology of one ligand are not necessarily important for another ligand.

GNAI3: Another Candidate Gene to Screen in Persons with Ocular Albinism

Ocular albinism type 1 (OA), caused by mutations in the OA1 gene, encodes a G-protein coupled receptor, OA1, localized in melanosomal membranes of the retinal pigment epithelium (RPE). This disorder is characterized by both RPE macro-melanosomes and abnormal decussation of ganglion cell axons at the brain's optic chiasm. We demonstrated previously that Oa1 specifically activates Gαi3, which also signals in the Oa1 transduction pathway that regulates melanosomal biogenesis. In this study, we screened the human Gαi3 gene, GNAI3, in DNA samples from 26 patients who had all clinical characteristics of OA but in whom a specific mutation in the OA1 gene had not been found, and in 6 normal control individuals. Using the Agilent HaloPlex Target Enrichment System and next-generation sequencing (NGS) on the Illumina MiSeq platform, we identified 518 variants after rigorous filtering. Many of these variants were corroborated by Sanger sequencing. Overall, 98.8% coverage of the GNAI3 gene was obtained by the HaloPlex amplicons. Of all variants, 6 non-synonymous and 3 synonymous were in exons, 41 in a non-coding exon embedded in the 3' untranslated region (UTR), 6 in the 5' UTR, and 462 in introns. These variants included novel SNVs, insertions, deletions, and a frameshift mutation. All were found in at least one patient but none in control samples. Using computational methods, we modeled the GNAI3 protein and its non-synonymous exonic mutations and determined that several of these may be the cause of disease in the patients studied. Thus, we have identified GNAI3 as a second gene possibly responsible for X-linked ocular albinism.

Use of Fluorescence Indicators in Receptor Ligands

Fluorescence techniques can provide insights into the environment of fluorescence indicators incorporated within a ligand as it is bound to its receptor. Fluorescence indicators of different sizes and chemical characteristics can provide insights into the nature of the binding environment, the surrounding structures, and even into conformational changes associated with receptor activation. Methods for determining fluorescence spectral analysis, fluorescence quenching, fluorescence anisotropy, fluorescence lifetimes, and red edge excitation shifts of the ligand probes are described. The applications of these techniques to the CCK1 receptor occupied by alexa(488)-CCK and aladan-CCK, as examples of probes developed (1) by derivatization of an existing peptide and (2) by incorporation during peptide synthesis, are utilized as examples. These methods represent powerful tools to expand our understanding of the structure and molecular basis of ligand activation of G protein-coupled receptors.

Activation and conformational dynamics of a class B G-protein-coupled glucagon receptor

The human glucagon receptor (GCGR) is a class B G-protein-coupled receptor (GPCR). The GCGR can be activated by glucagon and regulates the release of glucose. The GCGR has been proposed to be an important drug target for type 2 diabetes. Based on the structural model of a full-length glucagon-bound GCGR (glu-GCGR), we performed accelerated molecular dynamics (aMD) simulations, potential of mean force (PMF) calculations, cross-correlation analysis and community network analysis to study the activation mechanism and the conformational dynamics during the activation process. The PMF map depicts three different conformational states of the GCGR: the inactive, intermediate and active states. The activation of the GCGR is characterized by the outward movement of the intracellular side of helix VI. In the active state of the GCGR, the Arg173(2.46)-Ser350(6.41) and Glu245(3.50)-Thr351(6.42) hydrogen bonds break, and the χ1 rotamer of Phe322(5.54) changes from perpendicular to parallel to helix VI. The binding of the agonist glucagon decreases the correlated motions of the extracellular loops (ELCs) and the helices around the glucagon-binding site. During the activation of the GCGR, the connections between the intracellular sides of helices become weaker, and the connections between glucagon and ECLs and the extracellular sides of helices become stronger. These facilitate G-protein coupling on the intracellular side and glucagon binding on the extracellular side, and stabilize the GCGR in the active state. We expect that this study can provide useful information on the activation mechanism of the GCGR and facilitate the future design of GCGR inhibitors.

Allosteric effects of sodium ion binding on activation of the m3 muscarinic g-protein-coupled receptor

G-protein-coupled receptors (GPCRs) are important membrane proteins that mediate cellular signaling and represent primary targets for about one-third of currently marketed drugs. Recent x-ray crystallographic studies identified distinct conformations of GPCRs in the active and inactive states. An allosteric sodium ion was found bound to a highly conserved D2.50 residue in inactive GPCRs, whereas the D2.50 allosteric pocket became collapsed in active GPCR structures. However, the dynamic mechanisms underlying these observations remain elusive. In this study, we aimed to understand the mechanistic effects of sodium ion binding on dynamic activation of the M3 muscarinic GPCR through long-timescale accelerated molecular dynamics (aMD) simulations. Results showed that with the D2.50 residue deprotonated, the M3 receptor is bound by an allosteric sodium ion and confined mostly in the inactive state with remarkably reduced flexibility. In contrast, the D2.50-protonated receptor does not exhibit sodium ion binding to the D2.50 allosteric site and samples a significantly larger conformational space. The receptor activation is captured and characterized by large-scale structural rearrangements of the transmembrane helices via dynamic hydrogen bond and salt bridge interactions. The residue motions are highly correlated during receptor activation. Further network analysis revealed that the allosteric signaling between residue D2.50 and key residues in the intracellular, extracellular, and orthosteric pockets is significantly weakened upon sodium ion binding.

How Can Mutations Thermostabilize G-Protein-Coupled Receptors?

Structures of over 30 different G-protein-coupled receptors (GPCRs) have advanced our understanding of cell signaling and have provided a foundation for structure-guided drug design. This exciting progress has required the development of three complementary methods to facilitate GPCR crystallization, one of which is the thermostabilization of receptors by systematic mutagenesis. However, the reason why a particular mutation, or combination of mutations, stabilizes the receptor is not always evident from a static crystal structure. Molecular dynamics (MD) simulations have been used to identify and estimate the energetic factors that affect thermostability through comparing the dynamics of the thermostabilized receptors with structure-based models of the wild-type receptor. The data indicate that receptors are stabilized through a combination of factors, including an increase in receptor rigidity, a decrease in collective motion, reduced stress at specific residues, and the presence of ordered water molecules. Predicting thermostabilizing mutations computationally represents a major challenge for the field.

The mode of agonist binding to a G protein-coupled receptor switches the effect that voltage changes have on signaling

Signaling by many heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs) is either enhanced or attenuated by changes in plasma membrane potential. To identify structural correlates of the voltage sensitivity of GPCR signaling, we chose muscarinic acetylcholine receptors (the M1, M3, and M5 isoforms) as a model system. We combined molecular docking analysis with Förster resonance energy transfer (FRET)-based assays that monitored receptor activity under voltage clamp conditions. When human embryonic kidney (HEK) 293 cells expressing the individual receptors were stimulated with the agonist carbachol, membrane depolarization enhanced signaling by the M1 receptor but attenuated signaling by the M3 and M5 receptors. Furthermore, whether membrane depolarization enhanced or inhibited receptor signaling depended on the type of agonist. Membrane depolarization attenuated M3 receptor signaling when the receptor was bound to carbachol or acetylcholine, whereas depolarization enhanced signaling when the receptor was bound to either choline or pilocarpine. Docking calculations predicted that there were two distinct binding modes for these ligands, which were associated with the effect of depolarization on receptor function. From these calculations, we identified a residue in the M3 receptor that, when mutated, would alter the binding mode of carbachol to resemble that of pilocarpine in silico. Introduction of this mutated M3 receptor into cells confirmed that the membrane depolarization enhanced, rather than attenuated, signaling by the carbachol-bound receptor. Together, these data suggest that the directionality of the voltage sensitivity of GPCR signaling is defined by the specific binding mode of each ligand to the receptor.

Non-equivalence of Key Positively Charged Residues of the Free Fatty Acid 2 Receptor in the Recognition and Function of Agonist Versus Antagonist Ligands

Short chain fatty acids (SCFAs) are produced in the gut by bacterial fermentation of poorly digested carbohydrates. A key mediator of their actions is the G protein-coupled free fatty acid 2 (FFA2) receptor, and this has been suggested as a therapeutic target for the treatment of both metabolic and inflammatory diseases. However, a lack of understanding of the molecular determinants dictating how ligands bind to this receptor has hindered development. We have developed a novel radiolabeled FFA2 antagonist to probe ligand binding to FFA2, and in combination with mutagenesis and molecular modeling studies, we define how agonist and antagonist ligands interact with the receptor. Although both agonist and antagonist ligands contain negatively charged carboxylates that interact with two key positively charged arginine residues in transmembrane domains V and VII of FFA2, there are clear differences in how these interactions occur. Specifically, although agonists require interaction with both arginine residues to bind the receptor, antagonists require an interaction with only one of the two. Moreover, different chemical series of antagonist interact preferentially with different arginine residues. A homology model capable of rationalizing these observations was developed and provides a tool that will be invaluable for identifying improved FFA2 agonists and antagonists to further define function and therapeutic opportunities of this receptor.

Modulation of Glucagon Receptor Pharmacology by Receptor Activity-modifying Protein-2 (RAMP2)

The glucagon and glucagon-like peptide-1 (GLP-1) receptors play important, opposing roles in regulating blood glucose levels. Consequently, these receptors have been identified as targets for novel diabetes treatments. However, drugs acting at the GLP-1 receptor, although having clinical efficacy, have been associated with severe adverse side-effects, and targeting of the glucagon receptor has yet to be successful. Here we use a combination of yeast reporter assays and mammalian systems to provide a more complete understanding of glucagon receptor signaling, considering the effect of multiple ligands, association with the receptor-interacting protein receptor activity-modifying protein-2 (RAMP2), and the role of individual G protein α-subunits. We demonstrate that RAMP2 alters both ligand selectivity and G protein preference of the glucagon receptor. Importantly, we also uncover novel cross-reactivity of therapeutically used GLP-1 receptor ligands at the glucagon receptor that is abolished by RAMP2 interaction. This study reveals the glucagon receptor as a previously unidentified target for GLP-1 receptor agonists and highlights a role for RAMP2 in regulating its pharmacology. Such previously unrecognized functions of RAMPs highlight the need to consider all receptor-interacting proteins in future drug development.

Quantifying conformational changes in GPCRs: glimpse of a common functional mechanism

Background

G-protein-coupled receptors (GPCRs) are important drug targets and a better understanding of their molecular mechanisms would be desirable. The crystallization rate of GPCRs has accelerated in recent years as techniques have become more sophisticated, particularly with respect to Class A GPCRs interacting with G-proteins. These developments have made it possible for a quantitative analysis of GPCR geometrical features and binding-site conformations, including a statistical comparison between Class A GPCRs in active (agonist-bound) and inactive (antagonist-bound) states.

Results

Here we implement algorithms for the analysis of interhelical angles, distances, interactions and binding-site volumes in the transmembrane domains of 25 Class A GPCRs (7 active and 18 inactive). Two interhelical angles change in a statistically significant way between average inactive and active states: TM3-TM6 (by -9°) and TM6-TM7 (by +12°). A third interhelical angle: TM5-TM6 shows a trend, changing by -9°. In the transition from inactive to active states, average van der Waals interactions between TM3 and TM7 significantly increase as the average distance between them decreases by >2 Å. Average H-bonding between TM3 and TM6 decreases but is seemingly compensated by an increase in H-bonding between TM5 and TM6. In five Class A GPCRs, crystallized in both active and inactive states, increased H-bonding of agonists to TM6 and TM7, relative to antagonists, is observed. These protein-agonist interactions likely favour a change in the TM6-TM7 angle, which creates a narrowing in the binding pocket of activated receptors and an average ~200 ų reduction in volume.

Conclusions

In terms of similar conformational changes and agonist binding pattern, Class A GPCRs appear to share a common mechanism of activation, which can be exploited in future drug development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0567-3) contains supplementary material, which is available to authorized users.

Solution structure of the transmembrane 2 domain of the human melanocortin-4 receptor in sodium dodecyl sulfate (SDS) micelles and the functional implication of the D90N mutant

The melanocortin receptors (MCRs) are members of the G protein-coupled receptor (GPCR) 1 superfamily with seven transmembrane (TM) domains. Among them, the melanocortin-4 receptor (MC4R) subtype has been highlighted recently by genetic studies in obese humans. In particular, in a patient with severe early-onset obesity, a novel heterozygous mutation in the MC4R gene was found in an exchange of Asp to Asn in the 90th amino acid residue located in the TM 2 domain (MC4RD90N). Mutations in the MC4R gene are the most frequent monogenic causes of severe obesity and are described as heterozygous with loss of function. We determine solution structures of the TM 2 domain of MC4R (MC4RTM2) and compared secondary structure of Asp90 mutant (MC4RTM2-D90N) in a micelle environment by nuclear magnetic resonance (NMR) spectroscopy. NMR structure shows that MC4RTM2 forms a long α-helix with a kink at Gly98. Interestingly, the structure of MC4RTM2-D90N is similar to that of MC4RTM2 based on data from CD and NMR spectrum. However, the thermal stability and homogeneity of MC4RD90N is quite different from those of MC4R. The structure from molecular modeling suggests that Asp90(2.50) plays a key role in allosteric sodium ion binding. Our data suggest that the sodium ion interaction of Asp90(2.50) in the allosteric pocket of MC4R is essential to its function, explaining the loss of function of the MC4RD90N mutant.

Proton transfer-mediated GPCR activation

G-protein coupled receptors (GPCRs) play essential roles in signal transduction from the environment into the cell. While many structural features have been elucidated in great detail, a common functional mechanism on how the ligand-binding signal is converted into a conformational change on the cytoplasmic face resulting in subsequent activation of downstream effectors remain to be established. Based on available structural and functional data of the activation process in class-A GPCRs, we propose here that a change in protonation status, together with proton transfer via conserved structural elements located in the transmembrane region, are the key elements essential for signal transduction across the membrane.

Development of a highly selective allosteric antagonist radioligand for the type 1 cholecystokinin receptor and elucidation of its molecular basis of binding

Understanding the molecular basis of ligand binding to receptors provides insights useful for rational drug design. This work describes development of a new antagonist radioligand of the type 1 cholecystokinin receptor (CCK1R), (2-fluorophenyl)-2,3-dihydro-3-[(3-isoquinolinylcarbonyl)amino]-6-methoxy-2-oxo-l-H-indole-3-propanoate (T-0632), and exploration of the molecular basis of its binding. This radioligand bound specifically with high affinity within an allosteric pocket of CCK1R. T-0632 fully inhibited binding and action of CCK at this receptor, while exhibiting no saturable binding to the closely related type 2 cholecystokinin receptor (CCK2R). Chimeric CCK1R/CCK2R constructs were used to explore the molecular basis of T-0632 binding. Exchanging exonic regions revealed the functional importance of CCK1R exon 3, extending from the bottom of transmembrane segment (TM) 3 to the top of TM5, including portions of the intramembranous pocket as well as the second extracellular loop region (ECL2). However, CCK1R mutants in which each residue facing the pocket was changed to that present in CCK2R had no negative impact on T-0632 binding. Extending the chimeric approach to ECL2 established the importance of its C-terminal region, and site-directed mutagenesis of each nonconserved residue in this region revealed the importance of Ser(208) at the top of TM5. A molecular model of T-0632-occupied CCK1R was consistent with these experimental determinants, also identifying Met(121) in TM3 and Arg(336) in TM6 as important. Although these residues are conserved in CCK2R, mutating them had a distinct impact on the two closely related receptors, suggesting differential orientation. This establishes the molecular basis of binding of a highly selective nonpeptidyl allosteric antagonist of CCK1R, illustrating differences in docking that extend beyond determinants attributable to distinct residues lining the intramembranous pocket in the two receptor subtypes.

Structure-based drug design for G protein-coupled receptors

Our understanding of the structural biology of G protein-coupled receptors has undergone a transformation over the past 5 years. New protein-ligand complexes are described almost monthly in high profile journals. Appreciation of how small molecules and natural ligands bind to their receptors has the potential to impact enormously how medicinal chemists approach this major class of receptor targets. An outline of the key topics in this field and some recent examples of structure- and fragment-based drug design are described. A table is presented with example views of each G protein-coupled receptor for which there is a published X-ray structure, including interactions with small molecule antagonists, partial and full agonists. The possible implications of these new data for drug design are discussed.