Functional microdomains in G-protein-coupled receptors. The conserved arginine-cage motif in the gonadotropin-releasing hormone receptor

Gonadotropin-Releasing Hormone Receptor (GnRHR) and Hypogonadotropic Hypogonadism

Human sexual and reproductive development is regulated by the hypothalamic-pituitary-gonadal (HPG) axis, which is primarily controlled by the gonadotropin-releasing hormone (GnRH) acting on its receptor (GnRHR). Dysregulation of the axis leads to conditions such as congenital hypogonadotropic hypogonadism (CHH) and delayed puberty. The pathophysiology of GnRHR makes it a potential target for treatments in several reproductive diseases and in congenital adrenal hyperplasia. GnRHR belongs to the G protein-coupled receptor family and its GnRH ligand, when bound, activates several complex and tissue-specific signaling pathways. In the pituitary gonadotrope cells, it triggers the G protein subunit dissociation and initiates a cascade of events that lead to the production and secretion of the luteinizing hormone (LH) and follicle-stimulating hormone (FSH) accompanied with the phospholipase C, inositol phosphate production, and protein kinase C activation. Pharmacologically, GnRHR can be modulated by synthetic analogues. Such analogues include the agonists, antagonists, and the pharmacoperones. The agonists stimulate the gonadotropin release and lead to receptor desensitization with prolonged use while the antagonists directly block the GnRHR and rapidly reduce the sex hormone production. Pharmacoperones include the most recent GnRHR therapeutic approaches that directly correct the misfolded GnRHRs, which are caused by genetic mutations and hold serious promise for CHH treatment. Understanding of the GnRHR's genomic and protein structure is crucial for the most appropriate assessing of the mutation impact. Such mutations in the GNRHR are linked to normosmic hypogonadotropic hypogonadism and lead to various clinical symptoms, including delayed puberty, infertility, and impaired sexual development. These mutations vary regarding their mode of inheritance and can be found in the homozygous, compound heterozygous, or in the digenic state. GnRHR expression extends beyond the pituitary gland, and is found in reproductive tissues such as ovaries, uterus, and prostate and non-reproductive tissues such as heart, muscles, liver and melanoma cells. This comprehensive review explores GnRHR's multifaceted role in human reproduction and its clinical implications for reproductive disorders.

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.

Class A and C GPCR Dimers in Neurodegenerative Diseases

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

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

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

NPF activates a specific NPF receptor and regulates food intake in Pacific abalone Haliotis discus hannai

Neuropeptides function through G protein-coupled receptors (GPCRs) with high specificity, implying a significant degree of neuropeptide-GPCR coevolution. However, potential neuropeptide signaling systems in non-chordates are relatively elusive. We determined the specificity of the neuropeptide F (Hdh-NPF) signaling system with a cognate receptor (Hdh-NPFR) in the Pacific abalone, Haliotis discus hannai. Phylogenetic and exon–intron arrangement analyses of bilaterian NPF and the chordate ortholog NPY with their receptor sequences revealed a likely common ancestor, and Hdh-NPFR was similar to the NPYR2 subtype among the NPYR1, NPYR2, and NPYR5 subtypes. Among four Hdh-NPFR-related receptors, Hdh-NPFR specifically responded to Hdh-NPF peptide, supported by the dose–response luciferase reporter curve, intracellular Ca²⁺ mobilization, and phosphorylation of ERK1/2 and its inhibition with a protein kinase C inhibitor. Peptide fragmentations and shuffling of Hdh-NPF with human NPY could not activate the cellular response of Hdh-NPFR. Three-dimensional in silico modeling suggested that interaction of Hdh-NPF C-terminal amino acids with the extracellular loops of Hdh-NPFR is critical for Hdh-NPFR activation. In vivo injection of Hdh-NPF peptide increased food consumption, and knockdown of Hdh-NPF expression decreased food consumption in Pacific abalone. These findings provide evidence for co-evolution of the NPF/Y ligand-receptor system, enabling further research on mollusk orexigenic neuropeptides.

Structure reveals the activation mechanism of the MC4 receptor to initiate satiation signaling

Obesity is a global epidemic that causes morbidity and impaired quality of life. The melanocortin receptor 4 (MC4R) is at the crux of appetite, energy homeostasis, and body-weight control in the central nervous system and is a prime target for anti-obesity drugs. Here, we present the cryo-electron microscopy (cryo-EM) structure of the human MC4R-Gs signaling complex bound to the agonist setmelanotide, a cyclic peptide recently approved for the treatment of obesity. The work reveals the mechanism of MC4R activation, highlighting a molecular switch that initiates satiation signaling. In addition, our findings indicate that calcium (Ca2+) is required for agonist, but not antagonist, efficacy. These results fill a gap in the understanding of MC4R activation and could guide the design of future weight-management drugs.

Gonadotropin-Releasing Hormone and GnRH Receptor: Structure, Function and Drug Development

BACKGROUND

Gonadotropin-Releasing Hormone (GnRH) is a key element in sexual maturation and regulation of the reproductive cycle in the human organism. GnRH interacts with the pituitary cells through the activation of the Gonadotropin Releasing Hormone Receptors (GnRHR). Any impairments/dysfunctions of the GnRH-GnRHR complex lead to the development of various cancer types and disorders. Furthermore, the identification of GnRHR as a potential drug target has led to the development of agonist and antagonist molecules implemented in various treatment protocols. The development of these drugs was based on the information derived from the functional studies of GnRH and GnRHR.

OBJECTIVE

This review aims at shedding light on the versatile function of GnRH and GnRH receptor and offers an apprehensive summary regarding the development of different agonists, antagonists and non-peptide GnRH analogues.

CONCLUSION

The information derived from these studies can enhance our understanding of the GnRH-GnRHR versatile nature and offer valuable insight into the design of new more potent molecules.

Computational Investigations on the Binding Mode of Ligands for the Cannabinoid-Activated G Protein-Coupled Receptor GPR18

GPR18 is an orphan G protein-coupled receptor (GPCR) expressed in cells of the immune system. It is activated by the cannabinoid receptor (CB) agonist ∆⁹-tetrahydrocannabinol (THC). Several further lipids have been proposed to act as GPR18 agonists, but these results still require unambiguous confirmation. In the present study, we constructed a homology model of the human GPR18 based on an ensemble of three GPCR crystal structures to investigate the binding modes of the agonist THC and the recently reported antagonists which feature an imidazothiazinone core to which a (substituted) phenyl ring is connected via a lipophilic linker. Docking and molecular dynamics simulation studies were performed. As a result, a hydrophobic binding pocket is predicted to accommodate the imidazothiazinone core, while the terminal phenyl ring projects towards an aromatic pocket. Hydrophobic interaction of Cys251 with substituents on the phenyl ring could explain the high potency of the most potent derivatives. Molecular dynamics simulation studies suggest that the binding of imidazothiazinone antagonists stabilizes transmembrane regions TM1, TM6 and TM7 of the receptor through a salt bridge between Asp118 and Lys133. The agonist THC is presumed to bind differently to GPR18 than to the distantly related CB receptors. This study provides insights into the binding mode of GPR18 agonists and antagonists which will facilitate future drug design for this promising potential drug target.

Cell-autonomous light sensitivity via Opsin3 regulates fuel utilization in brown adipocytes

Opsin3 (Opn3) is a transmembrane heptahelical G protein–coupled receptor (GPCR) with the potential to produce a nonvisual photoreceptive effect. Interestingly, anatomical profiling of GPCRs reveals that Opn3 mRNA is highly expressed in adipose tissue. The photosensitive functions of Opn3 in mammals are poorly understood, and whether Opn3 has a role in fat is entirely unknown. In this study, we found that Opn3-knockout (Opn3-KO) mice were prone to diet-induced obesity and insulin resistance. At the cellular level, Opn3-KO brown adipocytes cultured in darkness had decreased glucose uptake and lower nutrient-induced mitochondrial respiration than wild-type (WT) cells. Light exposure promoted mitochondrial activity and glucose uptake in WT adipocytes but not in Opn3-KO cells. Brown adipocytes carrying a defective mutation in Opn3’s putative G protein–binding domain also exhibited a reduction in glucose uptake and mitochondrial respiration in darkness. Using RNA-sequencing, we identified several novel light-sensitive and Opn3-dependent molecular signatures in brown adipocytes. Importantly, direct exposure of brown adipose tissue (BAT) to light in living mice significantly enhanced thermogenic capacity of BAT, and this effect was diminished in Opn3-KO animals. These results uncover a previously unrecognized cell-autonomous, light-sensing mechanism in brown adipocytes via Opn3-GPCR signaling that can regulate fuel metabolism and mitochondrial respiration. Our work also provides a molecular basis for developing light-based treatments for obesity and its related metabolic disorders.

Brown adipose tissue plays a pivotal role in energy homeostasis and serves as a metabolic sink for glucose and fatty acid. This study demonstrates a novel light-sensing mechanism in mice via the photoreceptor Opsin3 that regulates fuel utilization and mitochondrial activity of brown adipocytes.

Towards A Molecular Understanding of The Cannabinoid Related Orphan Receptor GPR18: A Focus on Its Constitutive Activity

The orphan G-protein coupled receptor (GPCR), GPR18, has been recently proposed as a potential member of the cannabinoid family as it recognizes several endogenous, phytogenic, and synthetic cannabinoids. Potential therapeutic applications for GPR18 include intraocular pressure, metabolic disorders, and cancer. GPR18 has been reported to have high constitutive activity, i.e., activation/signaling occurs in the absence of an agonist. This activity can be reduced significantly by the A3.39N mutation. At the intracellular (IC) ends of (transmembrane helices) TMH3 and TMH6 in GPCRs, typically, a pair of oppositely charged amino acids form a salt bridge called the "ionic lock". Breaking of this salt bridge creates an IC opening for coupling with G protein. The GPR18 "ionic lock" residues (R3.50/S6.33) can form only a hydrogen bond. In this paper, we test the hypothesis that the high constitutive activity of GPR18 is due to the weakness of its "ionic lock" and that the A3.39N mutation strengthens this lock. To this end, we report molecular dynamics simulations of wild-type (WT) GPR18 and the A3.39N mutant in fully hydrated (POPC) phophatidylcholine lipid bilayers. Results suggest that in the A3.39N mutant, TMH6 rotates and brings R3.50 and S6.33 closer together, thus strengthening the GPR18 "ionic lock".

Glu2.53(90) of the GnRH receptor is part of the conserved G protein-coupled receptor structure and does not form a salt-bridge with Lys3.32(121)

GnRH receptor mutations, Glu^2.53(90)Lys and Glu^2.53(90)Asp, cause congenital hypogonadotropic hypogonadism. The Glu^2.53(90) side-chain has been proposed to form an intramolecular salt-bridge with Lys^3.32(121), but conserved intramolecular interaction networks in G protein-coupled receptor crystal structures predict that it interacts with Ser^3.35(124) and Trp^6.48(280). We investigated interhelical interactions of Glu^2.53(90) that stabilise GnRH receptor folding using functional analyses and computational modelling of mutant receptors. The Glu^2.53(90)Asp mutant was non-functional, but mutants with hydrophobic amino acids or Arg substituted for Glu^2.53(90) were functional, excluding a salt-bridge interaction. The Glu^2.53(90)Arg and Trp^6.48(280)Arg mutants had decreased affinity for GnRH. Models showed that congenital Glu^2.53(90)Lys and Glu^2.53(90)Asp mutations disrupt interactions with Ser^3.35(124) and Trp^6.48(280) respectively, whereas the Glu^2.53(90)Arg and Trp^6.48(280)Arg mutations preserve intramolecular contacts, but increase distance between the transmembrane helices. Our results show that disruption of interhelical contacts that are conserved in G protein-coupled receptors accounts for the effects of some disease-associated GnRH receptor mutations.

GPR62 constitutively activates cAMP signaling but is dispensable for male fertility in mice

G-protein-coupled receptors (GPCRs) participate in diverse physiological functions and are promising targets for drug discovery. However, there are still over 140 orphan GPCRs whose functions remain to be elucidated. Gpr62 is one of the orphan GPCRs that is expressed in the rat and human brain. In this study, we found that Gpr62 is also expressed in male germ cells in mice, and its expression increases along with sperm differentiation. GPR62 lacks the BBXXB and DRY motifs, which are conserved across many GPCRs, and it was able to induce cAMP signaling in the absence of a ligand. These structural and functional features are conserved among mammals, and the mutant analysis of GPR62 has revealed that lacking of these motifs is involved in the constitutive activity. We also found that GPR62 can homooligomerize, but it is not sufficient for its constitutive activity. We further investigated its physiological function by using Gpr62 knockout (Gpr62^-/-) mice. Gpr62^-/- mice were born normally and did not show any abnormality in growth and behavior. In addition, both male and female Gp62^-/- mice were fertile, and the differentiation and motility of spermatozoa were normal. We also found that Gpr61, the gene most similar to Gpr62 in the GPCR family shows a constitutive activity and an expression pattern similar to those of Gpr62 Our results suggest that GPR62 constitutively activates the cAMP pathway in male germ cells but is dispensable for male fertility, which is probably due to its functional redundancy with GPR61.

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.

Modeling and high-throughput experimental data uncover the mechanisms underlying Fshb gene sensitivity to gonadotropin-releasing hormone pulse frequency

Neuroendocrine control of reproduction by brain-secreted pulses of gonadotropin-releasing hormone (GnRH) represents a longstanding puzzle about extracellular signal decoding mechanisms. GnRH regulates the pituitary gonadotropin's follicle-stimulating hormone (FSH) and luteinizing hormone (LH), both of which are heterodimers specified by unique β subunits (FSHβ/LHβ). Contrary to Lhb, Fshb gene induction has a preference for low-frequency GnRH pulses. To clarify the underlying regulatory mechanisms, we developed three biologically anchored mathematical models: 1) parallel activation of Fshb inhibitory factors (e.g. inhibin α and VGF nerve growth factor-inducible), 2) activation of a signaling component with a refractory period (e.g. G protein), and 3) inactivation of a factor needed for Fshb induction (e.g. growth differentiation factor 9). Simulations with all three models recapitulated the Fshb expression levels obtained in pituitary gonadotrope cells perifused with varying GnRH pulse frequencies. Notably, simulations altering average concentration, pulse duration, and pulse frequency revealed that the apparent frequency-dependent pattern of Fshb expression in model 1 actually resulted from variations in average GnRH concentration. In contrast, models 2 and 3 showed "true" pulse frequency sensing. To resolve which components of this GnRH signal induce Fshb, we developed a high-throughput parallel experimental system. We analyzed over 4,000 samples in experiments with varying near-physiological GnRH concentrations and pulse patterns. Whereas Egr1 and Fos genes responded only to variations in average GnRH concentration, Fshb levels were sensitive to both average concentration and true pulse frequency. These results provide a foundation for understanding the role of multiple regulatory factors in modulating Fshb gene activity.

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.

Gonadotropin-Releasing Hormone (GnRH) Receptor Structure and GnRH Binding

Gonadotropin-releasing hormone (GnRH) regulates reproduction. The human GnRH receptor lacks a cytoplasmic carboxy-terminal tail but has amino acid sequence motifs characteristic of rhodopsin-like, class A, G protein-coupled receptors (GPCRs). This review will consider how recent descriptions of X-ray crystallographic structures of GPCRs in inactive and active conformations may contribute to understanding GnRH receptor structure, mechanism of activation and ligand binding. The structures confirmed that ligands bind to variable extracellular surfaces, whereas the seven membrane-spanning α-helices convey the activation signal to the cytoplasmic receptor surface, which binds and activates heterotrimeric G proteins. Forty non-covalent interactions that bridge topologically equivalent residues in different transmembrane (TM) helices are conserved in class A GPCR structures, regardless of activation state. Conformation-independent interhelical contacts account for a conserved receptor protein structure and their importance in the GnRH receptor structure is supported by decreased expression of receptors with mutations of residues in the network. Many of the GnRH receptor mutations associated with congenital hypogonadotropic hypogonadism, including the Glu^2.53(90) Lys mutation, involve amino acids that constitute the conserved network. Half of the ~250 intramolecular interactions in GPCRs differ between inactive and active structures. Conformation-specific interhelical contacts depend on amino acids changing partners during activation. Conserved inactive conformation-specific contacts prevent receptor activation by stabilizing proximity of TM helices 3 and 6 and a closed G protein-binding site. Mutations of GnRH receptor residues involved in these interactions, such as Arg^3.50(139) of the DRY/S motif or Tyr^7.53(323) of the N/DPxxY motif, increase or decrease receptor expression and efficiency of receptor coupling to G protein signaling, consistent with the native residues stabilizing the inactive GnRH receptor structure. Active conformation-specific interhelical contacts stabilize an open G protein-binding site. Progress in defining the GnRH-binding site has recently slowed, with evidence that Tyr^6.58(290) contacts Tyr⁵ of GnRH, whereas other residues affect recognition of Trp³ and Gly¹⁰NH₂. The surprisingly consistent observations that GnRH receptor mutations that disrupt GnRH binding have less effect on "conformationally constrained" GnRH peptides may now be explained by crystal structures of agonist-bound peptide receptors. Analysis of GPCR structures provides insight into GnRH receptor function.

In silico Exploration of the Conformational Universe of GPCRs

The structural plasticity of G protein coupled receptors (GPCRs) leads to a conformational universe going from inactive to active receptor states with several intermediate states. Many of them have not been captured yet and their role for GPCR activation is not well understood. The study of this conformational space and the transition dynamics between different receptor populations is a major challenge in molecular biophysics. The rational design of effector molecules that target such receptor populations allows fine-tuning receptor signalling with higher specificity to produce drugs with safer therapeutic profiles. In this minireview, we outline highly conserved receptor regions which are considered determinant for the establishment of distinct receptor states. We then discuss in-silico approaches such as dimensionality reduction methods and Markov State Models to explore the GPCR conformational universe and exploit the obtained conformations through structure-based drug design.

Intrinsic Relative Stabilities of the Neutral Tautomers of Arginine Side-Chain Models

The specific protonation state of amino acids is crucial for the physicochemical properties of proteins and their biological functions. These protonation states influence, for instance, properties related to hydrogen bonding, solubility, and folding. pKa calculations for proteins are, therefore, important and require, in principle, a specification of the most stable protonated and deprotonated forms of each titratable group. This is complicated by the existence of multiple tautomers, like the five neutral tautomers of the guanidine moiety in arginine. In this study, the compounds N-methyl-guanidine and N-ethyl-guanidine were used to model the charged and all neutral protonation states of the arginine side chain. The relative stabilities of all five neutral tautomers were investigated systematically for the first time, using quantum-mechanical calculations. These relative stabilities were obtained in vacuo, water and chloroform, by combining the quantum-mechanical calculations with a continuum solvation model. The water model was used to represent arginines exposed to an aqueous solution, whereas the chloroform model has a polarity representative of a protein core or a membrane. This allowed determining the relative pKa's associated with each neutral tautomer in these environments. A key result is that significant differences in stability are found between the neutral tautomers, in both water and chloroform. Some tautomers are consistently found to be the most stable. These findings will be helpful to refine pKa calculations in proteins.

Milestones in the Activation of a G Protein-Coupled Receptor. Insights from Molecular-Dynamics Simulations into the Human Cholecystokinin Receptor-1

Activation of G protein-coupled receptors (GPCRs) obeys an allosteric mechanism triggered by ligand binding. To understand how the signal is transduced in the cell, identification of the milestones paving the pathway between the active and the inactive states of the receptor is necessary. A model of the human cholecystokinin receptor-1 (CCK1R) has been proposed recently. The complex formed by CCK1R and an agonist ligand will serve as a paradigm of an active conformation to capture milestones in GPCR activation. To reach this goal, assuming microreversibility, the initial step toward the inactivation of CCK1R was modeled using free energy calculations, whereby the ligand is removed from the binding pocket. However accurate the reproduction of the experimental affinity constant, this simulation only represents an embryonic stage of the inactivation process. Starting from the apo receptor, an unprecedented 0.1-μs molecular dynamics trajectory was generated, bereft of experimental biases, bringing into the light key events in the inactivation of CCK1R, chief among which the hydration of its internal cavity, concomitant with the spatial rearrangement of the transmembrane helical segments. Hydration is intimately related to the isomerization of the highly conserved residue W326 of helix VI, acting as a two-state toggle switch, and of residue M121 of helix III. In the active state, the former residue obstructs the crevice, thereby preventing water leakage, which would otherwise trigger the disruption of an ionic lock between helices II and III involving the signature E/DRY motif ubiquitous to GPCRs.

Selectivity in the Use of Gi/o Proteins Is Determined by the DRF Motif in CXCR6 and Is Cell-Type Specific

CXCR6, the receptor for CXCL16, is expressed on multiple cell types and can be a coreceptor for human immunodeficiency virus 1. Except for CXCR6, all human chemokine receptors contain the D(3.49)R(3.50)Y(3.51) sequence, and all but two contain A(3.53) at the cytoplasmic terminus of the third transmembrane helix (H3C), a region within class A G protein-coupled receptors that contacts G proteins. In CXCR6, H3C contains D(3.49)R(3.50)F(3.51)I(3.52)V(3.53) at positions 126-130. We investigated the importance and interdependence of the canonical D126 and the noncanonical F128 and V130 in CXCR6 by mutating D126 to Y, F128 to Y, and V130 to A singly and in combination. For comparison, we mutated the analogous positions D142, Y144, and A146 to Y, F, and V, respectively, in CCR6, a related receptor containing the canonical sequences. Mutants were analyzed in both human embryonic kidney 293T and Jurkat E6-1 cells. Our data show that for CXCR6 and/or CCR6, mutations in H3C can affect both receptor signaling and chemokine binding; noncanonical H3C sequences are functionally linked, with dual changes mitigating the effects of single mutations; mutations in H3C that compromise receptor activity show selective defects in the use of individual Gi/o proteins; and the effects of mutations in H3C on receptor function and selectivity in Gi/o protein use can be cell-type specific. Our findings indicate that the ability of CXCR6 to make promiscuous use of the available Gi/o proteins is exquisitely dependent on sequences within the H3C and suggest that the native sequence allows for preservation of this function across different cellular environments.

Identification and characterization of cichlid TAAR genes and comparison with other teleost TAAR repertoires

BACKGROUND

TAARs (trace amine-associated receptors) are among the principal receptors expressed by the olfactory epithelium. We used the recent BROAD Institute release of the genome sequences of five representative fishes of the cichlid family to establish the complete TAAR repertoires of these species and to compare them with five other fish TAAR repertoires.

RESULTS

The genome sequences of O. niloticus, P. nyererei, H. burtoni, N. brichardi and M. zebra were analyzed by exhaustive TBLASTN searches with a set of published TAAR gene sequences used as positive bait. A second TBLASTN analysis was then performed on the candidate genes, with a set of non-TAAR class A GPCR (G protein-coupled receptors) used as negative bait. The resulting cichlid repertoire contained 44 complete TAAR genes from O. niloticus, 18 from P. nyererei, 23 from H. burtoni, 12 from N. brichardi and 20 from M. zebra, plus a number of pseudogenes, edge genes and fragments. A large proportion of these sequences (80%) consisted of two coding exons, separated in all but two cases by an intron in the interloop 1 coding sequence. We constructed phylogenetic trees. These trees indicated that TAARs constitute a distinct clade, well separated from ORs (olfactory receptors) and other class A GPCRs. Also these repertoires consist of several families and subfamilies, a number of which are common to fugu, tetraodon, stickleback and medaka. Like all other TAARs identified to date, cichlid TAARs have a characteristic two-dimensional structure and contain a number of amino-acid motifs or amino acids, such cysteine, in particular conserved positions.

CONCLUSIONS

Little is known about the functions of TAARs: in most cases their ligands have yet to be identified, partly because appropriate methods for such investigations have not been developed. Sequences analyses and comparisons of TAARs in several animal species, here fishes living in the same environment, should help reveal their roles and whether they are complementary to that of ORs.

Histidine(7.36(305)) in the conserved peptide receptor activation domain of the gonadotropin releasing hormone receptor couples peptide binding and receptor activation

Transmembrane helix seven residues of G protein-coupled receptors (GPCRs) couple agonist binding to a conserved receptor activation mechanism. Amino-terminal residues of the GnRH peptide determine agonist activity. We investigated GnRH interactions with the His(7.36(305)) residue of the GnRH receptor, using functional and computational analysis of modified GnRH receptors and peptides. Non-polar His(7.36(305)) substitutions decreased receptor affinity for GnRH four- to forty-fold, whereas GnRH signaling potency was more decreased (~150-fold). Uncharged polar His(7.36(305)) substitutions decreased GnRH potency, but not affinity. [2-Nal(3)]-GnRH retained high affinity at receptors with non-polar His(7.36(305)) substitutions, supporting a role for His(7.36(305)) in recognizing Trp(3) of GnRH. Compared with GnRH, [2-Nal(3)]-GnRH potency was lower at the wild type GnRH receptor, but unchanged or higher at mutant receptors. Results suggest that His(7.36(305)) of the GnRH receptor forms two distinct interactions that determine binding to Trp(3) and couple agonist binding to the conserved transmembrane domain network that activates GPCRs.

Mutations in the 'DRY' motif of the CB1 cannabinoid receptor result in biased receptor variants

The role of the highly conserved 'DRY' motif in the signaling of the CB1 cannabinoid receptor (CB1R) was investigated by inducing single-, double-, and triple-alanine mutations into this site of the receptor. We found that the CB1R-R3.50A mutant displays a partial decrease in its ability to activate heterotrimeric Go proteins (∼80% of WT CB1R (CB1R-WT)). Moreover, this mutant showed an enhanced basal β-arrestin2 (β-arr2) recruitment. More strikingly, the double-mutant CB1R-D3.49A/R3.50A was biased toward β-arrs, as it gained a robustly increased β-arr1 and β-arr2 recruitment ability compared with the WT receptor, while its G-protein activation was decreased. In contrast, the double-mutant CB1R-R3.50A/Y3.51A proved to be G-protein-biased, as it was practically unable to recruit β-arrs in response to agonist stimulus, while still activating G-proteins, although at a reduced level (∼70% of CB1R-WT). Agonist-induced ERK1/2 activation of the CB1R mutants showed a good correlation with their β-arr recruitment ability but not with their G-protein activation or inhibition of cAMP accumulation. Our results suggest that G-protein activation and β-arr binding of the CB1R are mediated by distinct receptor conformations, and the conserved 'DRY' motif plays different roles in the stabilization of these conformations, thus mediating both G-protein- and β-arr-mediated functions of CB1R.

A functional selectivity mechanism at the serotonin-2A GPCR involves ligand-dependent conformations of intracellular loop 2

With recent progress in determination of G protein-coupled receptor (GPCR) structure with crystallography, a variety of other experimental approaches (e.g., NMR spectroscopy, fluorescent-based assays, mass spectrometry techniques) are also being used to characterize state-specific and ligand-specific conformational states. MD simulations offer a powerful complementary approach to elucidate the dynamic features associated with ligand-specific GPCR conformations. To shed light on the conformational elements and dynamics of the important aspect of GPCR functional selectivity, we carried out unbiased microsecond-length MD simulations of the human serotonin 2A receptor (5-HT(2A)R) in the absence of ligand and bound to four distinct serotonergic agonists. The 5-HT(2A)R is a suitable system to study the structural features involved in the ligand-dependent conformational heterogeneity of GPCRs because it is well-characterized experimentally and exhibits a strong agonist-specific phenotype in that some 5-HT(2A)R agonists induce LSD-like hallucinations, while others lack this psychoactive property entirely. Here we report evidence for structural and dynamic differences in 5-HT(2A)R interacting with such pharmacologically distinct ligands, hallucinogens, and nonhallucinogens obtained from all-atom MD simulations. Differential ligand binding contacts were identified for structurally similar hallucinogens and nonhallucinogens and found to correspond to different conformations in the intracellular loop 2 (ICL2). From the different ICL2 conformations, functional selective phenotypes are suggested through effects on dimerization and/or distinct direct interaction with effector proteins. The findings are presented in the context of currently proposed hallucinogenesis mechanisms, and ICL2 is proposed as a fine-tuning selective switch that can differentiates modes of 5-HT(2A)R activation.

Understanding the effects on constitutive activation and drug binding of a D130N mutation in the β2 adrenergic receptor via molecular dynamics simulation

G-protein-coupled receptors (GPCRs) are currently one of the largest families of drug targets. The constitutive activation induced by mutation of key GPCR residues is associated closely with various diseases. However, the structural basis underlying such activation and its role in drug binding has remained unclear. Herein, we used all-atom molecular dynamics simulations and free energy calculations to study the effects of a D130N mutation on the structure of β2 adrenergic receptor (β2AR) and its binding of the agonist salbutamol. The results indicate that the mutation caused significant changes in some key helices. In particular, the mutation leads to the departure of transmembrane 3 (TM3) from transmembrane 6 (TM6) and marked changes in the NPxxY region as well as the complete disruption of a key ionic lock, all of which contribute to the observed constitutive activation. In addition, the D130N mutation weakens some important H-bonds, leading to structural changes in these regions. Binding free energy calculations indicate that van der Waals and electrostatic interactions are the main driving forces in binding salbutamol; however, binding strength in the mutant β2AR is significantly enhanced mainly through modifying electrostatic interactions. Further analysis revealed that the increase in binding energy upon mutation stems mainly from the H-bonds formed between the hydroxyl group of salbutamol and the serine residues of TM5. This observation suggests that modifications of the H-bond groups of this drug could significantly influence drug efficacy in the treatment of diseases associated with this mutation.

A yeast screening method to decipher the interaction between the adenosine A2B receptor and the C-terminus of different G protein α-subunits

The expression of human G protein-coupled receptors (GPCRs) in Saccharomyces cerevisiae containing chimeric yeast/mammalian Gα subunits provides a useful tool for the study of GPCR activation. In this study, we used a one-GPCR-one-G protein yeast screening method in combination with molecular modeling and mutagenesis studies to decipher the interaction between GPCRs and the C-terminus of different α-subunits of G proteins. We chose the human adenosine A2B receptor (hA2BR) as a paradigm, a typical class A GPCR that shows promiscuous behavior in G protein coupling in this yeast system. The wild-type hA2BR and five mutant receptors were expressed in 8 yeast strains with different humanized G proteins, covering the four major classes: Gαi, Gαs, Gαq, and Gα12. Our experiments showed that a tyrosine residue (Y) at the C-terminus of the Gα subunit plays an important role in controlling the activation of GPCRs. Receptor residues R103(3.50) and I107(3.54) are vital too in G protein-coupling and the activation of the hA2BR, whereas L213(IL3) is more important in G protein inactivation. Substitution of S235(6.36) to alanine provided the most divergent G protein-coupling profile. Finally, L236(6.37) substitution decreased receptor activation in all G protein pathways, although to a different extent. In conclusion, our findings shed light on the selectivity of receptor/G protein coupling, which may help in further understanding GPCR signaling.

How the dynamic properties and functional mechanisms of GPCRs are modulated by their coupling to the membrane environment

Experimental observations of the dependence of function and organization of G protein-coupled receptors (GPCRs) on their lipid environment have stimulated new quantitative studies of the coupling between the proteins and the membrane. It is important to develop such a quantitative understanding at the molecular level because the effects of the coupling are seen to be physiologically and clinically significant. Here we review findings that offer insight into how membrane-GPCR coupling is connected to the structural characteristics of the GPCR, from sequence to 3D structural detail, and how this coupling is involved in the actions of ligands on the receptor. The application of a recently developed computational approach designed for quantitative evaluation of membrane remodeling and the energetics of membrane-protein interactions brings to light the importance of the radial asymmetry of the membrane-facing surface of GPCRs in their interaction with the surrounding membrane. As the radial asymmetry creates adjacencies of hydrophobic and polar residues at specific sites of the GPCR, the ability of membrane remodeling to achieve complete hydrophobic matching is limited, and the residual mismatch carries a significant energy cost. The adjacencies are shown to be affected by ligand-induced conformational changes. Thus, functionally important organization of GPCRs in the cell membrane can depend both on ligand-determined properties and on the lipid composition of various membrane regions with different remodeling capacities. That this functionally important reorganization can be driven by oligomerization patterns that reduce the energy cost of the residual mismatch, suggests a new perspective on GPCR dimerization and ligand-GPCR interactions. The relation between the modulatory effects on GPCRs from the binding of specific cell-membrane components, e.g., cholesterol, and those produced by the non-local energetics of hydrophobic mismatch are discussed in this context.

Amino acid conservation and interactions in rhodopsin: probing receptor activation by NMR spectroscopy

Rhodopsin is a classical two-state G protein-coupled receptor (GPCR). In the dark, its 11-cis retinal chromophore serves as an inverse agonist to lock the receptor in an inactive state. Retinal-protein and protein-protein interactions have evolved to reduce the basal activity of the receptor in order to achieve low dark noise in the visual system. In contrast, absorption of light triggers rapid isomerization of the retinal, which drives the conversion of the receptor to a fully active conformation. Several specific protein-protein interactions have evolved that maintain the lifetime of the active state in order to increase the sensitivity of this receptor for dim-light vision in vertebrates. In this article, we review the molecular interactions that stabilize rhodopsin in the dark-state and describe the use of solid-state NMR spectroscopy for probing the structural changes that occur upon light-activation. Amino acid conservation provides a guide for those interactions that are common in the class A GPCRs as well as those that are unique to the visual system. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

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.

Structural features of the G-protein/GPCR interactions

BACKGROUND

The details of the functional interaction between G proteins and the G protein coupled receptors (GPCRs) have long been subjected to extensive investigations with structural and functional assays and a large number of computational studies.

SCOPE OF REVIEW

The nature and sites of interaction in the G-protein/GPCR complexes, and the specificities of these interactions selecting coupling partners among the large number of families of GPCRs and G protein forms, are still poorly defined.

MAJOR CONCLUSIONS

Many of the contact sites between the two proteins in specific complexes have been identified, but the three dimensional molecular architecture of a receptor-Gα interface is only known for one pair. Consequently, many fundamental questions regarding this macromolecular assembly and its mechanism remain unanswered.

GENERAL SIGNIFICANCE

In the context of current structural data we review the structural details of the interfaces and recognition sites in complexes of sub-family A GPCRs with cognate G-proteins, with special emphasis on the consequences of activation on GPCR structure, the prevalence of preassembled GPCR/G-protein complexes, the key structural determinants for selective coupling and the possible involvement of GPCR oligomerization in this process.

International Union of Basic and Clinical Pharmacology. [corrected]. LXXXVII. Complement peptide C5a, C4a, and C3a receptors

The activation of the complement cascade, a cornerstone of the innate immune response, produces a number of small (74-77 amino acid) fragments, originally termed anaphylatoxins, that are potent chemoattractants and secretagogues that act on a wide variety of cell types. These fragments, C5a, C4a, and C3a, participate at all levels of the immune response and are also involved in other processes such as neural development and organ regeneration. Their primary function, however, is in inflammation, so they are important targets for the development of antiinflammatory therapies. Only three receptors for complement peptides have been found, but there are no satisfactory antagonists as yet, despite intensive investigation. In humans, there is a single receptor for C3a (C3a receptor), no known receptor for C4a, and two receptors for C5a (C5a₁ receptor and C5a₂ receptor). The most recently characterized receptor, the C5a₂ receptor (previously known as C5L2 or GPR77), has been regarded as a passive binding protein, but signaling activities are now ascribed to it, so we propose that it be formally identified as a receptor and be given a name to reflect this. Here, we describe the complex biology of the complement peptides, introduce a new suggested nomenclature, and review our current knowledge of receptor pharmacology.

Type-1 cannabinoid receptor signaling in neuronal development

The type-1 cannabinoid receptor (CB1R) was initially identified as the neuronal target of Δ(9)-tetrahydrocannabinol (THC), the major psychoactive substance of marijuana. This receptor is one of the most abundant G-protein-coupled receptors in the adult brain, the target of endocannabinoid ligands and a well-characterized retrograde synaptic regulator. However, CB1Rs are also highly and often transiently expressed in neuronal populations in the embryonic and early postnatal brain, even before the formation of synapses. This suggests important physiological roles for CB1Rs during neuronal development. Several recent reviews have summarized our knowledge about the role of the endocannabinoid (eCB) system in neurodevelopment and neurotransmission by focusing on the metabolism of endocannabinoid molecules. Here, we review current knowledge about the effects of the modulation of CB1R signaling during the different phases of brain development. More precisely, we focus on reports that directly implicate CB1Rs during progenitor cell migration and differentiation, neurite outgrowth, axonal pathfinding and synaptogenesis. Based on theoretical considerations and on the reviewed experimental data, we propose a new model to explain the diversity of experimental findings on eCB signaling on neurite growth and axonal pathfinding. In our model, cell-autonomus and paracrine eCBs acting on CB1Rs are part of a global inhibitory network of cytoskeletal effectors, which act in concert with positive-feedback local-excitation loops, to ultimately yield highly polarized neurons.

Molecular cloning and pharmacological characterization of two novel GnRH receptors in the lamprey (Petromyzon marinus)

This paper reports the identification, expression, binding kinetics, and functional studies of two novel type III lamprey GnRH receptors (lGnRH-R-2 and lGnRH-R-3) in the sea lamprey, a basal vertebrate. These novel GnRH receptors share the structural features and amino acid motifs common to other known gnathostome GnRH receptors. The ligand specificity and activation of intracellular signaling studies showed ligands lGnRH-II and -III induced an inositol phosphate (IP) response at lGnRH-R-2 and lGnRH-R-3, whereas the ligand lGnRH-I did not stimulate an IP response. lGnRH-II was a more potent activator of lGnRH-R-3 than lGnRH-III. Stimulation of lGnRH-R-2 and lGnRH-R-3 testing all three lGnRH ligands did not elicit a cAMP response. lGnRH-R-2 has a higher binding affinity in response to lGnRH-III than lGnRH-II, whereas lGnRH-R-3 has a higher binding affinity in response to lGnRH-II than IGnRH-III. lGnRH-R-2 precursor transcript was detected in a wide variety of tissues including the pituitary whereas lGnRH-R-3 precursor transcript was not as widely expressed and primarily expressed in the brain and eye of male and female lampreys. From our phylogenetic analysis, we propose that lGnRH-R-1 evolved from a common ancestor of all vertebrate GnRH receptors and lGnRH-R-2 and lGnRH-R-3 likely occurred due to a gene duplication within the lamprey lineage. In summary, we propose from our findings of receptor subtypes in the sea lamprey that the evolutionary recruitment of specific pituitary GnRH receptor subtypes for particular physiological functions seen in later evolved vertebrates was an ancestral character that first arose in a basal vertebrate.

Congenital hypogonadotropic hypogonadism due to GnRH receptor mutations in three brothers reveal sites affecting conformation and coupling

Congenital hypogonadotropic hypogonadism (CHH) is characterized by low gonadotropins and failure to progress normally through puberty. Mutations in the gene encoding the GnRH receptor (GNRHR1) result in CHH when present as compound heterozygous or homozygous inactivating mutations. This study identifies and characterizes the properties of two novel GNRHR1 mutations in a family in which three brothers display normosmic CHH while their sister was unaffected. Molecular analysis in the proband and the affected brothers revealed two novel non-synonymous missense GNRHR1 mutations, present in a compound heterozygous state, whereas their unaffected parents possessed only one inactivating mutation, demonstrating the autosomal recessive transmission in this kindred and excluding X-linked inheritance equivocally suggested by the initial pedigree analysis. The first mutation at c.845 C>G introduces an Arg substitution for the conserved Pro 282 in transmembrane domain (TMD) 6. The Pro282Arg mutant is unable to bind radiolabeled GnRH analogue. As this conserved residue is important in receptor conformation, it is likely that the mutation perturbs the binding pocket and affects trafficking to the cell surface. The second mutation at c.968 A>G introduces a Cys substitution for Tyr 323 in the functionally crucial N/DPxxY motif in TMD 7. The Tyr323Cys mutant has an increased GnRH binding affinity but reduced receptor expression at the plasma membrane and impaired G protein-coupling. Inositol phosphate accumulation assays demonstrated absent and impaired Gα_q/11 signal transduction by Pro282Arg and Tyr323Cys mutants, respectively. Pretreatment with the membrane permeant GnRHR antagonist NBI-42902, which rescues cell surface expression of many GNRHR1 mutants, significantly increased the levels of radioligand binding and intracellular signaling of the Tyr323Cys mutant but not Pro282Arg. Immunocytochemistry confirmed that both mutants are present on the cell membrane albeit at low levels. Together these molecular deficiencies of the two novel GNRHR1 mutations lead to the CHH phenotype when present as a compound heterozygote.

Ligand-dependent conformations and dynamics of the serotonin 5-HT(2A) receptor determine its activation and membrane-driven oligomerization properties

From computational simulations of a serotonin 2A receptor (5-HT_2AR) model complexed with pharmacologically and structurally diverse ligands we identify different conformational states and dynamics adopted by the receptor bound to the full agonist 5-HT, the partial agonist LSD, and the inverse agonist Ketanserin. The results from the unbiased all-atom molecular dynamics (MD) simulations show that the three ligands affect differently the known GPCR activation elements including the toggle switch at W6.48, the changes in the ionic lock between E6.30 and R3.50 of the DRY motif in TM3, and the dynamics of the NPxxY motif in TM7. The computational results uncover a sequence of steps connecting these experimentally-identified elements of GPCR activation. The differences among the properties of the receptor molecule interacting with the ligands correlate with their distinct pharmacological properties. Combining these results with quantitative analysis of membrane deformation obtained with our new method (Mondal et al, Biophysical Journal 2011), we show that distinct conformational rearrangements produced by the three ligands also elicit different responses in the surrounding membrane. The differential reorganization of the receptor environment is reflected in (i)-the involvement of cholesterol in the activation of the 5-HT_2AR, and (ii)-different extents and patterns of membrane deformations. These findings are discussed in the context of their likely functional consequences and a predicted mechanism of ligand-specific GPCR oligomerization.

The 5-HT_2A receptor for the neurotransmitter serotonin (5-HT) belongs to family A (rhodopsin-like) G-protein coupled receptors (GPCRs), one of the most important classes of membrane proteins that are targeted by an extensive and diverse collection of external stimuli. Recently we learned that different ligands targeting the same GPCR can elicit different biological responses, but the mechanisms remain unknown. We address this fundamental question for the serotonin 5-HT_2A receptor, because it is known to respond to the binding of structurally diverse ligands by producing similar stimuli in the cell, and to the binding of quite similar ligands with dramatically different responses. Molecular dynamics simulations of molecular models of the serotonin 5-HT_2A receptor in complex with pharmacologically distinct ligands show the dynamic rearrangements of the receptor molecule to be different for these ligands, and the nature and extents of the rearrangements reflect the known pharmacological properties of the ligands as full, partial or inverse activators of the receptor. The different rearrangements of the receptor molecule are shown to produce different rearrangements of the surrounding membrane, a remodeling of the environment that can have differential ligand-determined effects on receptor function and association in the cell's membrane.

Allosteric modulator ORG27569 induces CB1 cannabinoid receptor high affinity agonist binding state, receptor internalization, and Gi protein-independent ERK1/2 kinase activation

The cannabinoid receptor 1 (CB1), a member of the class A G protein-coupled receptor family, is expressed in brain tissue where agonist stimulation primarily activates the pertussis toxin-sensitive inhibitory G protein (G(i)). Ligands such as CP55940 ((1R,3R,4R)-3-[2-hydroxy-4-(1,1-dimethylheptyl)phenyl]-4-(3- hydroxypropyl)cyclohexan-1-ol) and Δ(9)-tetrahydrocannabinol are orthosteric agonists for the receptor, bind the conventional binding pocket, and trigger G(i)-mediated effects including inhibition of adenylate cyclase. ORG27569 (5-chloro-3-ethyl-1H-indole-2-carboxylic acid [2-(4-piperidin-1-yl-phenyl)ethyl]amide) has been identified as an allosteric modulator that displays positive cooperativity for CP55940 binding to CB1 yet acts as an antagonist of G protein coupling. To examine this apparent conundrum, we used the wild-type CB1 and two mutants, T210A and T210I (D'Antona, A. M., Ahn, K. H., and Kendall, D. A. (2006) Biochemistry 45, 5606-5617), which collectively cover a spectrum of receptor states from inactive to partially active to more fully constitutively active. Using these receptors, we demonstrated that ORG27569 induces a CB1 receptor state that is characterized by enhanced agonist affinity and decreased inverse agonist affinity consistent with an active conformation. Also consistent with this conformation, the impact of ORG27569 binding was most dramatic on the inactive T210A receptor and less pronounced on the already active T210I receptor. Although ORG27569 antagonized CP55940-induced guanosine 5'-3-O-(thio)triphosphate binding, which is indicative of G protein coupling inhibition in a concentration-dependent manner, the ORG27569-induced conformational change of the CB1 receptor led to cellular internalization and downstream activation of ERK signaling, providing the first case of allosteric ligand-biased signaling via CB1. ORG27569-induced ERK phosphorylation persisted even after pertussis toxin treatment to abrogate G(i) and occurs in HEK293 and neuronal cells.

Evolution of vertebrate GnRH receptors from the perspective of a Basal vertebrate

This minireview provides the current status on gonadotropin-releasing hormone receptors (GnRH-R) in vertebrates, from the perspective of a basal vertebrate, the sea lamprey, and provides an evolutionary scheme based on the recent advance of whole genome sequencing. In addition, we provide a perspective on the functional divergence and evolution of the receptors. In this review we use the phylogenetic classification of vertebrate GnRH receptors that groups them into three clusters: type I (mammalian and non-mammalian), type II, and type III GnRH receptors. New findings show that the sea lamprey has two type III-like GnRH receptors and an ancestral type GnRH receptor that is more closely related to the type II-like receptors. These two novel GnRH receptors along with lGnRH-R-1 share similar structural features and amino acid motifs common to other known gnathostome type II/III receptors. Recent data analyses of the lamprey genome provide strong evidence that two whole rounds of genome duplication (2R) occurred prior to the gnathostome-agnathan split. Based on our current knowledge, it is proposed that lGnRH-R-1 evolved from an ancestor of the type II receptor following a vertebrate-shared genome duplication and that the two type III receptors resulted from a duplication within lamprey of a gene derived from a lineage shared by many vertebrates.

The structure of active opsin as a basis for identification of GPCR agonists by dynamic homology modelling and virtual screening assays

Most of the currently available G protein-coupled receptor (GPCR) crystal structures represent an inactive receptor state, which has been considered to be suitable only for the discovery of antagonists and inverse agonists in structure-based computational ligand screening. Using the β(2)-adrenergic receptor (B2AR) as a model system, we show that a dynamic homology model based on an "active" opsin structure without further incorporation of experimental data performs better than the crystal structure of the inactive B2AR in finding agonists over antagonists/inverse agonists. Such "active-like state" dynamic homology models can therefore be used to selectively identify GPCR agonists in in silico ligand libraries.

Conserved amino acids participate in the structure networks deputed to intramolecular communication in the lutropin receptor

The luteinizing hormone receptor (LHR) is a G protein-coupled receptor (GPCR) particularly susceptible to spontaneous pathogenic gain-of-function mutations. Protein structure network (PSN) analysis on wild-type LHR and two constitutively active mutants, combined with in vitro mutational analysis, served to identify key amino acids that are part of the regulatory network responsible for propagating communication between the extracellular and intracellular poles of the receptor. Highly conserved amino acids in the rhodopsin family GPCRs participate in the protein structural stability as network hubs in both the inactive and active states. Moreover, they behave as the most recurrent nodes in the communication paths between the extracellular and intracellular sides in both functional states with emphasis on the active one. In this respect, non-conservative loss-of-function mutations of these amino acids is expected to impair the most relevant way of communication between activating mutation sites or hormone-binding domain and G protein recognition regions.

Structural and functional roles of small group-conserved amino acids present on helix-H7 in the β(2)-adrenergic receptor

Sequence analysis of the class A G protein-coupled receptors (GPCRs) reveals that most of the highly conserved sites are located in the transmembrane helices. A second level of conservation exists involving those residues that are conserved as a group characterized by small and/or weakly polar side chains (Ala, Gly, Ser, Cys, Thr). These positions can have group conservation levels of up to 99% across the class A GPCRs and have been implicated in mediating helix-helix interactions in membrane proteins. We have previously shown that mutation of group-conserved residues present on transmembrane helices H2-H4 in the β(2)-adrenergic receptor (β(2)-AR) can influence both receptor expression and function. We now target the group-conserved sites, Gly315(7.42) and Ser319(7.46), on H7 for structure-function analysis. Replacing Ser319(7.46) with smaller amino acids (Ala or Gly) did not influence the ability of the mutant receptors to bind to the antagonist dihydroalprenolol (DHA) but resulted in ~15-20% agonist-independent activity. Replacement of Ser319(7.46) with the larger amino acid leucine lowered the expression of the S319L mutant and its ability to bind DHA. Both the G315A and G315S mutants also exhibited agonist-independent signaling, while the G315L mutant did not show specific binding to DHA. These data indicate that Gly315(7.42) and Ser319(7.46) are stabilizing β(2)-AR in an inactive conformation. We discuss our results in the context of van der Waals interactions of Gly315(7.42) with Trp286(6.48) and hydrogen bonding interactions of Ser319(7.46) with amino acids on H1-H2-H7 and with structural water.

Highly conserved tyrosine stabilizes the active state of rhodopsin

Light-induced isomerization of the 11-cis-retinal chromophore in the visual pigment rhodopsin triggers displacement of the second extracellular loop (EL2) and motion of transmembrane helices H5, H6, and H7 leading to the active intermediate metarhodopsin II (Meta II). We describe solid-state NMR measurements of rhodopsin and Meta II that target the molecular contacts in the region of the ionic lock involving these three helices. We show that a contact between Arg135(3.50) and Met257(6.40) forms in Meta II, consistent with the outward rotation of H6 and breaking of the dark-state Glu134(3.49)-Arg135(3.50)-Glu247(6.30) ionic lock. We also show that Tyr223(5.58) and Tyr306(7.53) form molecular contacts with Met257(6.40). Together these results reveal that the crystal structure of opsin in the region of the ionic lock reflects the active state of the receptor. We further demonstrate that Tyr223(5.58) and Ala132(3.47) in Meta II stabilize helix H5 in an active orientation. Mutation of Tyr223(5.58) to phenylalanine or mutation of Ala132(3.47) to leucine decreases the lifetime of the Meta II intermediate. Furthermore, the Y223F mutation is coupled to structural changes in EL2. In contrast, mutation of Tyr306(7.53) to phenylalanine shows only a moderate influence on the Meta II lifetime and is not coupled to EL2.

Functional differences of invariant and highly conserved residues in the extracellular domain of the glycoprotein hormone receptors

Multiple interactions exist between human follicle-stimulating hormone (FSH) and the N-terminal hormone-binding fragment of the human FSH receptor (FSHR) extracellular domain (ECD). Binding of the other human glycoprotein hormones to their cognate human receptors (luteinizing hormone receptor (LHR) and thyroid-stimulating hormone receptor (TSHR)) was expected to be similar. This study focuses on amino acid residues in β-strands 2 (Lys(74)), 4 (Tyr(124), Asn(129), and Thr(130)), and 5 (Asp(150) and Asp(153)) of the FSHR ECD identified in the human FSH·FSHR ECD crystal structure as contact sites with the common glycoprotein hormone α-subunit, and on noncontact residues in β-strands 2 (Ser(78)) and 8 (Asp(224) and Ser(226)) as controls. These nine residues are either invariant or highly conserved in LHR and TSHR. Mutagenesis and functional characterization of these residues in all three human receptors allowed an assessment of their contribution to binding and receptor activation. Surprisingly, the six reported α-subunit contact residues of the FSHR ECD could be replaced without significant loss of FSH binding, while cAMP signaling potency was diminished significantly with several replacements. Comparative studies of the homologous residues in LHR and TSHR revealed both similarities and differences. The results for FSH/FSHR were analyzed on the basis of the crystal structure of the FSH·FSHR ECD complex, and comparative modeling was used to generate structures for domains, proteins, and complexes for which no structures were available. Although structural information of hormone-receptor interaction allowed the identification of hormone-receptor contact sites, functional analysis of each contact site was necessary to assess its contribution to hormone binding and receptor activation.

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.