Conserved helix 7 tyrosine acts as a multistate conformational switch in the 5HT2C receptor. Identification of a novel "locked-on" phenotype and double revertant mutations

Identification and functional analysis of six melanocortin-4-receptor-like (MC4R-like) mutations in goldfish (Carassius auratus)

Melanocortin receptor-4 (MC4R) belongs to the G protein-coupled receptor family, characterized by a classical structure of seven transmembrane domains (7TMD). They play an important role in food intake and weight regulation. In the present study, we identified melanocortin-4-receptor-like (caMC4RL) mutants of goldfish from the Qian River in the Qin Ling region and characterized their functional properties, including the constitutive activities of the mutants, ligand-induced cAMP and ERK1/2 accumulation, and AMPK activation. The results show that six caMC4RL mutants were identified in goldfish from the Qian River in the Qin Ling region, and are located in the conserved position of the Cyprinidae MC4Rs. The mutations (E57K, P296S, and R302T/K) result in the loss of Gs signaling function. The mutations (P296 and R302T/K) exhibited biased signaling in response to ACTH stimulation in the MAPK/ERK pathway. In addition, the E57K mutant may play a role in weight regulation and could serve as molecular markers for molecular breeding. These data will provide fundamental information for functional studies of teleost GPCR mutants and MC4R isoforms.

The effect of cancer-associated mutations on ligand binding and receptor function - A case for the 5-HT2C receptor

The serotonin 5-HT2C receptor is a G protein-coupled receptor (GPCR) mainly expressed in the central nervous system. Besides regulating mood, appetite, and reproductive behavior, it has been identified as a potential target for cancer treatment. In this study, we aimed to investigate the effects of cancer patient-derived 5-HT2C receptor mutations on ligand binding and receptor functionality. By filtering the sequencing data from the Genomic Data Commons data portal (GDC), we selected 12 mutations from multiple cancer types. We found that the affinity of the endogenous agonist serotonin (5-HT) and inverse agonist mesulergine were both drastically decreased by mutations L209HECL2 and F328S6.52, which are located in the orthosteric binding pocket. In the calcium-flux assay, the potency of 5-HT was decreased at F328S6.52, while a trend of increased efficacy was observed. In contrast, 5-HT displayed higher affinity at E306K6.30 and E306A6.30, while a trend of decreased efficacy was observed. These two mutations may disrupt the conserved ionic interaction between E6.30 and R3.50, and thus increase the constitutive activity of the receptor. The inhibitory potency of mesulergine was increased at E306A6.30 but not E306K6.30. Lastly, P365H7.50 decreased the expression level of the receptor by more than ten-fold, which prevented further functional analyses. This study shows that cancer-associated mutations of 5-HT2C receptor have diverse effects on ligand binding and function. Such mutations may affect serotonin-mediated signaling in tumor cells as well as treatment strategies targeting this receptor.

Determinants of Neutral Antagonism and Inverse Agonism in the β2-Adrenergic Receptor

Free-energy profiles for the activation/deactivation of the β2-adrenergic receptor (ADRB2) with neutral antagonist and inverse agonist ligands have been determined with well-tempered multiple-walker (MW) metadynamics simulations. The inverse agonists carazolol and ICI118551 clearly favor single inactive conformational minima in both the binary and ternary ligand-receptor-G-protein complexes, in accord with the inverse-agonist activity of the ligands. The behavior of neutral antagonists is more complex, as they seem also to affect the recruitment of the G-protein. The results are analyzed in terms of the conformational states of the well-known microswitches that have been proposed as indicators of receptor activity.

All-Atom Molecular Dynamics Simulations Indicated the Involvement of a Conserved Polar Signaling Channel in the Activation Mechanism of the Type I Cannabinoid Receptor

The type I cannabinoid G protein-coupled receptor (CB1, GPCR) is an intensely investigated pharmacological target, owing to its involvement in numerous physiological functions as well as pathological processes such as cancers, neurodegenerative diseases, metabolic disorders and neuropathic pain. In order to develop modern medications that exert their effects through binding to the CB1 receptor, it is essential to understand the structural mechanism of activation of this protein. The pool of atomic resolution experimental structures of GPCRs has been expanding rapidly in the past decade, providing invaluable information about the function of these receptors. According to the current state of the art, the activity of GPCRs involves structurally distinct, dynamically interconverting functional states and the activation is controlled by a cascade of interconnecting conformational switches in the transmembrane domain. A current challenge is to uncover how different functional states are activated and what specific ligand properties are responsible for the selectivity towards those different functional states. Our recent studies of the μ-opioid and β₂-adrenergic receptors (MOP and β₂AR, respectively) revealed that the orthosteric binding pockets and the intracellular surfaces of these receptors are connected through a channel of highly conserved polar amino acids whose dynamic motions are in high correlation in the agonist- and G protein-bound active states. This and independent literature data led us to hypothesize that, in addition to consecutive conformational transitions, a shift of macroscopic polarization takes place in the transmembrane domain, which is furnished by the rearrangement of polar species through their concerted movements. Here, we examined the CB1 receptor signaling complexes utilizing microsecond scale, all-atom molecular dynamics (MD) simulations in order to see if our previous assumptions could be applied to the CB1 receptor too. Apart from the identification of the previously proposed general features of the activation mechanism, several specific properties of the CB1 have been indicated that could possibly be associated with the signaling profile of this receptor.

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.

Structural basis of sphingosine-1-phosphate receptor 1 activation and biased agonism

Sphingosine-1-phosphate receptor 1 (S1PR1) is a master regulator of lymphocyte egress from the lymph node and an established drug target for multiple sclerosis (MS). Mechanistically, therapeutic S1PR1 modulators activate the receptor yet induce sustained internalization through a potent association with β-arrestin. However, a structural basis of biased agonism remains elusive. Here, we report the cryo-electron microscopy (cryo-EM) structures of G_i-bound S1PR1 in complex with S1P, fingolimod-phosphate (FTY720-P) and siponimod (BAF312). In combination with functional assays and molecular dynamics (MD) studies, we reveal that the β-arrestin-biased ligands direct a distinct activation path in S1PR1 through the extensive interplay between the PIF and the NPxxY motifs. Specifically, the intermediate flipping of W269^6.48 and the retained interaction between F265^6.44 and N307^7.49 are the key features of the β-arrestin bias. We further identify ligand-receptor interactions accounting for the S1PR subtype specificity of BAF312. These structural insights provide a rational basis for designing novel signaling-biased S1PR modulators.

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

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

In Silico Identification of Cholesterol Binding Motifs in the Chemokine Receptor CCR3

CC motif chemokine receptor 3 (CCR3) is a Class A G protein-coupled receptor (GPCR) mainly responsible for the cellular trafficking of eosinophils. As such, it plays key roles in inflammatory conditions, such as asthma and arthritis, and the metastasis of many deadly forms of cancer. However, little is known about how CCR3 functionally interacts with its bilayer environment. Here, we investigate cholesterol binding sites in silico through Coarse-Grained Molecular Dynamics (MD) and Pylipid analysis using an extensively validated homology model based on the crystal structure of CCR5. These simulations identified several cholesterol binding sites containing Cholesterol Recognition/Interaction Amino Acid Consensus motif (CRAC) and its inversion CARC motifs in CCR3. One such site, a CARC site in TM1, in conjunction with aliphatic residues in TM7, emerged as a candidate for future investigation based on the cholesterol residency time within the binding pocket. This site forms the core of a cholesterol binding site previously observed in computational studies of CCR2 and CCR5. Most importantly, these cholesterol binding sites are conserved in other chemokine receptors and may provide clues to cholesterol regulation mechanisms in this subfamily of Class A GPCRs.

DEER Analysis of GPCR Conformational Heterogeneity

G protein-coupled receptors (GPCRs) represent a large class of transmembrane helical proteins which are involved in numerous physiological signaling pathways and therefore represent crucial pharmacological targets. GPCR function and the action of therapeutic molecules are defined by only a few parameters, including receptor basal activity, ligand affinity, intrinsic efficacy and signal bias. These parameters are encoded in characteristic receptor conformations existing in equilibrium and their populations, which are thus of paramount interest for the understanding of receptor (mal-)functions and rational design of improved therapeutics. To this end, the combination of site-directed spin labeling and EPR spectroscopy, in particular double electron-electron resonance (DEER), is exceedingly valuable as it has access to sub-Angstrom spatial resolution and provides a detailed picture of the number and populations of conformations in equilibrium. This review gives an overview of existing DEER studies on GPCRs with a focus on the delineation of structure/function frameworks, highlighting recent developments in data analysis and visualization. We introduce "conformational efficacy" as a parameter to describe ligand-specific shifts in the conformational equilibrium, taking into account the loose coupling between receptor segments observed for different GPCRs using DEER.

Correlated Motions of Conserved Polar Motifs Lay out a Plausible Mechanism of G Protein-Coupled Receptor Activation

Recent advancements in the field of experimental structural biology have provided high-resolution structures of active and inactive state G protein-coupled receptors (GPCRs), a highly important pharmaceutical target family, but the process of transition between these states is poorly understood. According to the current theory, GPCRs exist in structurally distinct, dynamically interconverting functional states of which populations are shifted upon binding of ligands and intracellular signaling proteins. However, explanation of the activation mechanism, on an entirely structural basis, gets complicated when multiple activation pathways and active receptor states are considered. Our unbiased, atomistic molecular dynamics simulations of the μ opioid receptor (MOP) revealed that transmission of external stimulus to the intracellular surface of the receptor is accompanied by subtle, concerted movements of highly conserved polar amino acid side chains along the 7th transmembrane helix. This may entail the rearrangement of polar species and the shift of macroscopic polarization in the transmembrane domain, triggered by agonist binding. Based on our observations and numerous independent indications, we suggest amending the widely accepted theory that the initiation event of GPCR activation is the shift of macroscopic polarization between the ortho- and allosteric binding pockets and the intracellular G protein-binding interface.

Characterizing Interhelical Interactions of G-Protein Coupled Receptors with the Fragment Molecular Orbital Method

G-protein coupled receptors (GPCRs) are the largest superfamily of membrane proteins, regulating almost every aspect of cellular activity and serving as key targets for drug discovery. We have identified an accurate and reliable computational method to characterize the strength and chemical nature of the interhelical interactions between the residues of transmembrane (TM) domains during different receptor activation states, something that cannot be characterized solely by visual inspection of structural information. Using the fragment molecular orbital (FMO) quantum mechanics method to analyze 35 crystal structures representing different branches of the class A GPCR family, we have identified 69 topologically equivalent TM residues that form a consensus network of 51 inter-TM interactions, providing novel results that are consistent with and help to rationalize experimental data. This discovery establishes a comprehensive picture of how defined molecular forces govern specific interhelical interactions which, in turn, support the structural stability, ligand binding, and activation of GPCRs.

Conformational plasticity of the intracellular cavity of GPCR-G-protein complexes leads to G-protein promiscuity and selectivity

While the dynamics of the intracellular surface in agonist-stimulated GPCRs is well studied, the impact of GPCR dynamics on G-protein selectivity remains unclear. Here, we combine molecular dynamics simulations with live-cell FRET and secondary messenger measurements, for 21 GPCR-G-protein combinations, to advance a dynamic model of the GPCR-G-protein interface. Our data show C terminus peptides of Gαs, Gαi, and Gαq proteins assume a small ensemble of unique orientations when coupled to their cognate GPCRs, similar to the variations observed in 3D structures of GPCR-G-protein complexes. The noncognate G proteins interface with latent intracellular GPCR cavities but dissociate due to weak and unstable interactions. Three predicted mutations in β2-adrenergic receptor stabilize binding of noncognate Gαq protein in its latent cavity, allowing promiscuous signaling through both Gαs and Gαq in a dose-dependent manner. This demonstrates that latent GPCR cavities can be evolved, by design or nature, to tune G-protein selectivity, giving insights to pluridimensional GPCR signaling.

Characterising GPCR-ligand interactions using a fragment molecular orbital-based approach

There has been fantastic progress in solving GPCR crystal structures. However, the ability of X-ray crystallography to guide the drug discovery process for GPCR targets is limited by the availability of accurate tools to explore receptor-ligand interactions. Visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum mechanical approaches (QM) are often too computationally expensive, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed and the ability to reveal key interactions that would otherwise be hard to detect. Integration of GPCR crystallography or homology modelling with FMO reveals atomistic details of the individual contributions of each residue and water molecule towards ligand binding, including an analysis of their chemical nature.

GPCRs: What Can We Learn from Molecular Dynamics Simulations?

Advances in the structural biology of G-protein Coupled Receptors have resulted in a significant step forward in our understanding of how this important class of drug targets function at the molecular level. However, it has also become apparent that they are very dynamic molecules, and moreover, that the underlying dynamics is crucial in shaping the response to different ligands. Molecular dynamics simulations can provide unique insight into the dynamic properties of GPCRs in a way that is complementary to many experimental approaches. In this chapter, we describe progress in three distinct areas that are particularly difficult to study with other techniques: atomic level investigation of the conformational changes that occur when moving between the various states that GPCRs can exist in, the pathways that ligands adopt during binding/unbinding events and finally, the influence of lipids on the conformational dynamics of GPCRs.

Exploring GPCR-Ligand Interactions with the Fragment Molecular Orbital (FMO) Method

The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity. It is essential for an efficient structure-based drug design (SBDD) process. FMO enables ab initio approaches to be applied to systems that conventional quantum-mechanical (QM) methods would find challenging. The key advantage of the Fragment Molecular Orbital Method (FMO) is that it can reveal atomistic details about the individual contributions and chemical nature of each residue and water molecule toward ligand binding which would otherwise be difficult to detect without using QM methods. In this chapter, we demonstrate the typical use of FMO to analyze 19 crystal structures of β1 and β2 adrenergic receptors with their corresponding agonists and antagonists.

Evolutionary action and structural basis of the allosteric switch controlling β2AR functional selectivity

Functional selectivity of G-protein-coupled receptors is believed to originate from ligand-specific conformations that activate only subsets of signaling effectors. In this study, to identify molecular motifs playing important roles in transducing ligand binding into distinct signaling responses, we combined in silico evolutionary lineage analysis and structure-guided site-directed mutagenesis with large-scale functional signaling characterization and non-negative matrix factorization clustering of signaling profiles. Clustering based on the signaling profiles of 28 variants of the β₂-adrenergic receptor reveals three clearly distinct phenotypical clusters, showing selective impairments of either the Gi or βarrestin/endocytosis pathways with no effect on Gs activation. Robustness of the results is confirmed using simulation-based error propagation. The structural changes resulting from functionally biasing mutations centered around the DRY, NPxxY, and PIF motifs, selectively linking these micro-switches to unique signaling profiles. Our data identify different receptor regions that are important for the stabilization of distinct conformations underlying functional selectivity.

Ligand-induced biased signaling is thought to result in part from ligand-specific receptor conformations that cause the engagement of distinct effectors. Here the authors trace and evaluate the impact of mutations of the β2–adrenergic receptor on multiple signaling outputs to provide structural-level insight into the determinants of GPCR functional selectivity.

Functional characterization of lysophosphatidic acid receptor 1 mutants identified in rat cancer tissues

Lysophosphatidic acid (LPA), an extracellular lipid mediator, exerts various cellular effects through activation of LPA receptors, LPA₁-LPA₆, in many types of cells including cancer cells. We recently found several missense mutations of Lpar1 in rat cancer tissues. One of these mutations is located at the extracellular tip of the seventh transmembrane domain of LPA₁, and another three mutations are found within the NPXXY motif in the seventh transmembrane domain. These mutants are designated F295S LPA₁ and P308S, I310T, and Y311H LPA₁, respectively. Here, we examined the functions of these LPA₁ mutants. Compared with wild-type (WT) LPA₁, F295S, P308S, and I310T LPA₁ showed decreased maximal responses in inhibition of cAMP formation, Ca²⁺ mobilization, and cytoskeletal changes. Y311H LPA₁ failed to show LPA-induced cellular responses. However, these LPA₁ mutants were internalized in response to LPA exposure. Finally, while WT and F295S LPA₁ showed a similar, broad distribution throughout the cell, P308S, I310T, and Y311H LPA₁ displayed a restricted cellular distribution and co-localized with the endoplasmic reticulum. These data suggest that the LPA₁ mutants perturb LPA signaling in cancer tissues.

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.

Application of an Integrated GPCR SAR-Modeling Platform To Explain the Activation Selectivity of Human 5-HT2C over 5-HT2B

Agonism of the 5-HT2C serotonin receptor has been associated with the treatment of a number of diseases including obesity, psychiatric disorders, sexual health, and urology. However, the development of effective 5-HT2C agonists has been hampered by the difficulty in obtaining selectivity over the closely related 5-HT2B receptor, agonism of which is associated with irreversible cardiac valvulopathy. Understanding how to design selective agonists requires exploration of the structural features governing the functional uniqueness of the target receptor relative to related off targets. X-ray crystallography, the major experimental source of structural information, is a slow and challenging process for integral membrane proteins, and so is currently not feasible for every GPCR or GPCR-ligand complex. Therefore, the integration of existing ligand SAR data with GPCR modeling can be a practical alternative to provide this essential structural insight. To demonstrate this, we integrated SAR data from 39 azepine series 5-HT2C agonists, comprising both selective and unselective examples, with our hierarchical GPCR modeling protocol (HGMP). Through this work we have been able to demonstrate how relatively small differences in the amino acid sequences of GPCRs can lead to significant differences in secondary structure and function, as supported by experimental data. In particular, this study suggests that conformational differences in the tilt of TM7 between 5-HT2B and 5-HT2C, which result from differences in interhelical interactions, may be the major source of selectivity in G-protein activation between these two receptors. Our approach also demonstrates how the use of GPCR models in conjunction with SAR data can be used to explain activity cliffs.

7TM Domain Structure of Adhesion GPCRs

Schematic presentation of the overall adhesion G Protein-Coupled Receptor (aGPCR) structure and functional domains, covering an extracellular N-terminal fragment (NTF), a membrane-spanning C-terminal fragment (CTF) and a GPCR proteolysis site (GPS). (Left side) aGPCR model constructed based on the seven-transmembrane (7TM) structure (blue) of secretin family glucagon receptor (GCGR) (PDB, 4L6R) [11] and the GPCR autoproteolysis inducing (GAIN) domain (magenta) structure of latrophilin 1 (PDB, 4DLQ) [9]. The β-13 strand residues are depicted in green. (Right side) The experimentally validated full-length secretin family GCGR structure combining structural and experimental information from the GCGR 7TM crystal structure (PDB, 4L6R) (blue), the GCGR extracellular domain (ECD) structure (PDB, 4ERS) (magenta) and the ECD structure of glucagon-like peptide-1 (GLP-1)-bound glucagon-like peptide-1 receptor (GLP-1R) (PDB, 3IOL) (green), complemented by site-directed mutagenesis, electron microscopy (EM), hydrogen-deuterium exchange (HDX) and cross-linking studies [11-13]) Despite the recent breakthroughs in the elucidation of the three-dimensional structures of the seven transmembrane (7TM) domain of the G protein-coupled receptor (GPCR) superfamily, a corresponding structure of a member of the adhesion GPCR (aGPCR) family has not yet been solved. In this chapter, we give an overview of the current knowledge of the 7TM domain of aGPCRs by comparative structure-based sequence similarity analyses between aGPCRs and GPCRs with known crystal structure. Of the GPCR superfamily, only the secretin family shares some sequence similarity with aGPCRs. This chapter will therefore emphasize on the comparison of these two GPCR families. Two 7TM domain structures of secretin family GPCRs are known that provide insight into the structure-function relationships of conserved sequence motifs that play important roles and are also present in most aGPCRs. This suggests that the 7TM domains of aGPCRs and secretin family GPCRs share a similar structural fold and that the conserved residues in both families may be involved in similar intermolecular interaction networks and facilitate similar conformational changes. Comparison of the residues that line the large peptide hormone binding pocket in the 7TM domain of secretin family GPCRs with corresponding residues in aGPCRs indicates that in the latter, the corresponding pocket in the 7TM domain is relatively hydrophobic and may be even larger. Improved knowledge on these conserved sequence motifs will help to understand the interactions of the aGPCR 7TM domain with ligands and gain insight into the activation mechanism of aGPCRs.

The Fragment Molecular Orbital Method Reveals New Insight into the Chemical Nature of GPCR–Ligand Interactions

Our interpretation of ligand-protein interactions is often informed by high-resolution structures, which represent the cornerstone of structure-based drug design. However, visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum Mechanics approaches are often too computationally expensive, but one method, Fragment Molecular Orbital (FMO), offers an excellent compromise and has the potential to reveal key interactions that would otherwise be hard to detect. To illustrate this, we have applied the FMO method to 18 Class A GPCR-ligand crystal structures, representing different branches of the GPCR genome. Our work reveals key interactions that are often omitted from structure-based descriptions, including hydrophobic interactions, nonclassical hydrogen bonds, and the involvement of backbone atoms. This approach provides a more comprehensive picture of receptor-ligand interactions than is currently used and should prove useful for evaluation of the chemical nature of ligand binding and to support structure-based drug design.

Transmembrane signal transduction by peptide hormones via family B G protein-coupled receptors

Although family B G protein-coupled receptors (GPCRs) contain only 15 members, they play key roles in transmembrane signal transduction of hormones. Family B GPCRs are drug targets for developing therapeutics for diseases ranging from metabolic to neurological disorders. Despite their importance, the molecular mechanism of activation of family B GPCRs remains largely unexplored due to the challenges in expression and purification of functional receptors to the quantity for biophysical characterization. Currently, there is no crystal structure available of a full-length family B GPCR. However, structures of key domains, including the extracellular ligand binding regions and seven-helical transmembrane regions, have been solved by X-ray crystallography and NMR, providing insights into the mechanisms of ligand recognition and selectivity, and helical arrangements within the cell membrane. Moreover, biophysical and biochemical methods have been used to explore functions, key residues for signaling, and the kinetics and dynamics of signaling processes. This review summarizes the current knowledge of the signal transduction mechanism of family B GPCRs at the molecular level and comments on the challenges and outlook for mechanistic studies of family B GPCRs.

Functions of DPLIY motif and helix 8 of human melanocortin-3 receptor

The melanocortin-3 receptor (MC3R) is a member of the family A G protein-coupled receptors (GPCRs). The MC3R remains the most enigmatic of the melanocortin receptors with regard to its physiological functions, especially its role in energy homeostasis. The N/DPxxY motif and the eighth helix (helix 8) in the carboxyl terminus of GPCRs have been identified to be important for receptor expression, ligand binding, signal transduction and internalization. To gain a better understanding of the structure-function relationship of MC3R, we performed a systematic study of all 20 residues in this domain using alanine-scanning mutagenesis. We showed that although all mutants were expressed normally on the cell surface, eleven residues were important for ligand binding and one was indispensable for downstream cAMP generation. F347A showed constitutive activity in cAMP signaling while all the other mutants had normal basal activities. We studied the signaling capacity of nine mutants in the ERK1/2 signaling pathway. All of these mutants showed normal basal ERK1/2 phosphorylation levels. The pERK1/2 levels of six binding- or signaling-defective mutants were enhanced upon agonist stimulation. The unbalanced cAMP and pERK1/2 signaling pathways suggested the existence of biased signaling in MC3R mutants. In summary, we showed that the DPLIY motif and helix 8 was important for MC3R activation and signal transduction. Our data led to a better understanding of the structure-function relationship of MC3R.

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.

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.

Probing the metabotropic glutamate receptor 5 (mGlu₅) positive allosteric modulator (PAM) binding pocket: discovery of point mutations that engender a "molecular switch" in PAM pharmacology

Positive allosteric modulation of metabotropic glutamate receptor subtype 5 (mGlu₅) is a promising novel approach for the treatment of schizophrenia and cognitive disorders. Allosteric binding sites are topographically distinct from the endogenous ligand (orthosteric) binding site, allowing for co-occupation of a single receptor with the endogenous ligand and an allosteric modulator. Negative allosteric modulators (NAMs) inhibit and positive allosteric modulators (PAMs) enhance the affinity and/or efficacy of the orthosteric agonist. The molecular determinants that govern mGlu₅ modulator affinity versus cooperativity are not well understood. Focusing on the modulators based on the acetylene scaffold, we sought to determine the molecular interactions that contribute to PAM versus NAM pharmacology. Generation of a comparative model of the transmembrane-spanning region of mGlu₅ served as a tool to predict and interpret the impact of mutations in this region. Application of an operational model of allosterism allowed for determination of PAM and NAM affinity estimates at receptor constructs that possessed no detectable radioligand binding as well as delineation of effects on affinity versus cooperativity. Novel mutations within the transmembrane domain (TM) regions were identified that had differential effects on acetylene PAMs versus 2-methyl-6-(phenylethynyl)-pyridine, a prototypical NAM. Three conserved amino acids (Y658, T780, and S808) and two nonconserved residues (P654 and A809) were identified as key determinants of PAM activity. Interestingly, we identified two point mutations in TMs 6 and 7 that, when mutated, engender a mode switch in the pharmacology of certain PAMs.

Molecular and pharmacological characterization of serotonin 5-HT2α and 5-HT7 receptors in the salivary glands of the blowfly Calliphora vicina

Secretion in blowfly (Calliphora vicina) salivary glands is stimulated by the biogenic amine serotonin (5-hydroxytryptamine, 5-HT), which activates both inositol 1,4,5-trisphosphate (InsP(3))/Ca(2+) and cyclic adenosine 3',5'-monophosphate (cAMP) signalling pathways in the secretory cells. In order to characterize the signal-inducing 5-HT receptors, we cloned two cDNAs (Cv5-ht2α, Cv5-ht7) that share high similarity with mammalian 5-HT(2) and 5-HT(7) receptor genes, respectively. RT-PCR demonstrated that both receptors are expressed in the salivary glands and brain. Stimulation of Cv5-ht2α-transfected mammalian cells with 5-HT elevates cytosolic [Ca(2+)] in a dose-dependent manner (EC(50) = 24 nM). In Cv5-ht7-transfected cells, 5-HT produces a dose-dependent increase in [cAMP](i) (EC(50) = 4 nM). We studied the pharmacological profile for both receptors. Substances that appear to act as specific ligands of either Cv5-HT(2α) or Cv5-HT(7) in the heterologous expression system were also tested in intact blowfly salivary gland preparations. We observed that 5-methoxytryptamine (100 nM) activates only the Cv5-HT(2α) receptor, 5-carboxamidotryptamine (300 nM) activates only the Cv5-HT(7) receptor, and clozapine (1 µM) antagonizes the effects of 5-HT via Cv5-HT(7) in blowfly salivary glands, providing means for the selective activation of each of the two 5-HT receptor subtypes. This study represents the first comprehensive molecular and pharmacological characterization of two 5-HT receptors in the blowfly and permits the analysis of the physiological role of these receptors, even when co-expressed in cells, and of the modes of interaction between the Ca(2+)- and cAMP-signalling cascades.

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

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

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.

The hydrophobic amino acids in putative helix 8 in carboxy-terminus of histamine H(3) receptor are involved in receptor-G-protein coupling

Functional roles of putative helix 8 in the carboxy-terminal tail of the human histamine H(3) receptor were investigated using deleted and alanine-substituted mutant receptors. While the deletion of the carboxy-terminal tail did not decrease the total expression level, surface expression, or ligand binding affinity, the agonist-stimulated cAMP response, [((35))S] GTPγS binding, and MAPK activation were totally abolished. The receptor lacking the carboxy-terminal tail also failed to respond to an inverse agonist, thioperamide, suggesting that the carboxy-terminal tail is involved in the regulation of receptor activity by changing G-protein coupling with the receptor. Site-directed mutagenesis revealed that hydrophobic amino acids in the putative helix 8 such as phenylalanines at position 419 (F7.60) and 423 (F7.64) or leucines at 426 (L7.67) and 427 (L7.68) were important for the agonist-induced activation of H(3) receptor. Substitution of F7.60 also resulted in a receptor that was less responsive to inactivation by the inverse agonist, implying the existence of an intermediate conformation that can be either activated or inactivated. Our results suggest that hydrophobic interface of putative helix 8 is important for the regulation of H(3) receptor activity, presumably by stabilizing the helix to the plasma membrane.

Structural determinants in the second intracellular loop of the human cannabinoid CB1 receptor mediate selective coupling to G(s) and G(i)

BACKGROUND AND PURPOSE

The cannabinoid CB(1) receptor is primarily thought to be functionally coupled to the G(i) form of G proteins, through which it negatively regulates cAMP accumulation. Here, we investigated the dual coupling properties of CB(1) receptors and characterized the structural determinants that mediate selective coupling to G(s) and G(i).

EXPERIMENTAL APPROACH

A cAMP-response element reporter gene system was employed to quantitatively analyze cAMP change. CB(1)/CB(2) receptor chimeras and site-directed mutagenesis combined with functional assays and computer modelling were used to determine the structural determinants mediating selective coupling to G(s) and G(i).

KEY RESULTS

CB(1) receptors could couple to both G(s)-mediated cAMP accumulation and G(i)-induced activation of ERK1/2 and Ca(2+) mobilization, whereas CB(2) receptors selectively coupled to G(i) and inhibited cAMP production. Using CB(1)/CB(2) chimeric receptors, the second intracellular loop (ICL2) of the CB(1) receptor was identified as primarily responsible for mediating G(s) and G(i) coupling specificity. Furthermore, mutation of Leu-222 in ICL2 to either Ala or Pro switched G protein coupling from G(s) to G(i), while to Ile or Val led to balanced coupling of the mutant receptor with G(s) and G(i) .

CONCLUSIONS AND IMPLICATIONS

The ICL2 of CB(1) receptors and in particular Leu-222, which resides within a highly conserved DRY(X)(5) PL motif, played a critical role in G(s) and G(i) protein coupling and specificity. Our studies provide new insight into the mechanisms governing the coupling of CB(1) receptors to G proteins and cannabinoid-induced tolerance.

Modeling GPCR active state conformations: the β(2)-adrenergic receptor

The recent publication of several G protein-coupled receptor (GPCR) structures has increased the information available for homology modeling inactive class A GPCRs. Moreover, the opsin crystal structure shows some active features. We have therefore combined information from these two sources to generate an extensively validated model of the active conformation of the β(2)-adrenergic receptor. Experimental information on fully active GPCRs from zinc binding studies, site-directed spin labeling, and other spectroscopic techniques has been used in molecular dynamics simulations. The observed conformational changes reside mainly in transmembrane helix 6 (TM6), with additional small but significant changes in TM5 and TM7. The active model has been validated by manual docking and is in agreement with a large amount of experimental work, including site-directed mutagenesis information. Virtual screening experiments show that the models are selective for β-adrenergic agonists over other GPCR ligands, for (R)- over (S)-β-hydroxy agonists and for β(2)-selective agonists over β(1)-selective agonists. The virtual screens reproduce interactions similar to those generated by manual docking. The C-terminal peptide from a model of the stimulatory G protein, readily docks into the active model in a similar manner to which the C-terminal peptide from transducin, docks into opsin, as shown in a recent opsin crystal structure. This GPCR-G protein model has been used to explain site-directed mutagenesis data on activation. The agreement with experiment suggests a robust model of an active state of the β(2)-adrenergic receptor has been produced. The methodology used here should be transferable to modeling the active state of other GPCRs.

The cannabinoid type-1 receptor carboxyl-terminus, more than just a tail

The cannabinoid type-1 (CB(1)) receptor is a G protein-coupled receptor that binds the main active ingredient of marijuana, Δ(9)-tetrahydrocannabinol, and has been implicated in several disease states, including drug addiction, anxiety, depression, obesity, and chronic pain. In the two decades since the discovery of CB(1), studies at the molecular level have centered on the transmembrane core. This interest has now expanded as we discover that other regions of CB(1), including the CB(1) carboxyl-terminus, have critical structures that are important for CB(1) activity and regulation. Following the recent description of the three dimensional structure of the full-length CB(1) carboxyl-terminal tail [Biopolymers (2009) vol. 91, pp. 565-573], several residues and structural motifs including two α-helices (termed H8 and H9) have been postulated to interact with common G protein-coupled receptor accessory proteins, such as G-proteins and β-arrestins. This discourse will focus on the CB(1) carboxyl-terminus; our current understanding of the structural features of this region, evidence for its interaction with proteins, and the impact of structure on the binding and regulatory function of CB(1) accessory proteins. The involvement of the carboxyl-terminus in the receptor life cycle including activation, desensitization, and internalization will be highlighted.

Helix 8 of the M1 muscarinic acetylcholine receptor: scanning mutagenesis delineates a G protein recognition site

We have used alanine-scanning mutagenesis followed by functional expression and molecular modeling to analyze the roles of the 14 residues, Asn422 to Cys435, C-terminal to transmembrane (TM) helix 7 of the M(1) muscarinic acetylcholine receptor. The results suggest that they form an eighth (H8) helix, associated with the cytoplasmic surface of the cell membrane in the active state of the receptor. We suggest that the amide side chain of Asn422 may act as a cap to the C terminus of TM7, stabilizing its junction with H8, whereas the side chain of Phe429 may restrict the relative movements of H8 and the C terminus of TM7 in the inactive ground state of the receptor. We have identified four residues, Phe425, Arg426, Thr428, and Leu432, which are important for G protein binding and signaling. These may form a docking site for the C-terminal helix of the G protein α subunit, and collaborate with G protein recognition residues elsewhere in the cytoplasmic domain of the receptor to form a coherent surface for G protein binding in the activated state of the receptor.

Human cannabinoid 1 GPCR C-terminal domain interacts with bilayer phospholipids to modulate the structure of its membrane environment

G protein-coupled receptors (GPCRs) play critical physiological and therapeutic roles. The human cannabinoid 1 GPCR (hCB1) is a prime pharmacotherapeutic target for addiction and cardiometabolic disease. Our prior biophysical studies on the structural biology of a synthetic peptide representing the functionally significant hCB1 transmembrane helix 7 (TMH7) and its cytoplasmic extension, helix 8 (H8), [hCB1(TMH7/H8)] demonstrated that the helices are oriented virtually perpendicular to each other in membrane-mimetic environments. We identified several hCB1(TMH7/H8) structure-function determinants, including multiple electrostatic amino-acid interactions and a proline kink involving the highly conserved NPXXY motif. In phospholipid bicelles, TMH7 structure, orientation, and topology relative to H8 are dynamically modulated by the surrounding membrane phospholipid bilayer. These data provide a contextual basis for the present solid-state NMR study to investigate whether intermolecular interactions between hCB1(TMH7/H8) and its phospholipid environment may affect membrane-bilayer structure. For this purpose, we measured (1)H-(13)C heteronuclear dipolar couplings for the choline, glycerol, and acyl-chain regions of dimyristoylphosphocholine in a magnetically aligned hCB1(TMH7/H8) bicelle sample. The results identify discrete regional interactions between hCB1(TMH7/H8) and membrane lipid molecules that increase phospholipid motion and decrease phospholipid order, indicating that the peptide's partial traversal of the bilayer alters membrane structure. These data offer new insight into hCB1(TMH7/H8) properties and support the concept that the membrane bilayer itself may serve as a mechanochemical mediator of hCB1/GPCR signal transduction. Since interaction with its membrane environment has been implicated in hCB1 function and its modulation by small-molecule therapeutics, our work should help inform hCB1 pharmacology and the design of hCB1-targeted drugs.

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

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

Identification of specific transmembrane residues and ligand-induced interface changes involved in homo-dimer formation of a yeast G protein-coupled receptor

The Saccharomyces cerevisiae alpha-factor pheromone receptor, Ste2p, has been studied as a model for G protein-coupled receptor (GPCR) structure and function. Dimerization has been demonstrated for many GPCRs, although the role(s) of dimerization in receptor function is disputed. Transmembrane domains one (TM1) and four (TM4) of Ste2p were shown previously to play a role in dimerization. In this study, single cysteine substitutions were introduced into a Cys-less Ste2p, and disulfide-mediated dimerization was assessed. Six residues in TM1 (L64 to M69) that had not been previously investigated and 19 residues in TM7 (T278 to A296) of which 15 were not previously investigated were mutated to create 25 single Cys-containing Ste2p molecules. Ste2p mutants V68C in TM1 and nine mutants in TM7 (cysteine substituted into residues 278, 285, 289, and 291 to 296) showed increased dimerization upon addition of an oxidizing agent in comparison to the background dimers formed by the Cys-less receptor. The formation of dimers was decreased for TM7 mutant receptors in the presence of alpha-factor indicating that ligand binding resulted in a conformational change that influenced dimerization. The effect of ligand on dimer formation suggests that dimers are formed in the resting state and the activated state of the receptor by different TM interactions.

Light activation of rhodopsin: insights from molecular dynamics simulations guided by solid-state NMR distance restraints

Structural restraints provided by solid-state NMR measurements of the metarhodopsin II intermediate are combined with molecular dynamics simulations to help visualize structural changes in the light activation of rhodopsin. Since the timescale for the formation of the metarhodopsin II intermediate (>1 ms) is beyond that readily accessible by molecular dynamics, we use NMR distance restraints derived from 13C dipolar recoupling measurements to guide the simulations. The simulations yield a working model for how photoisomerization of the 11-cis retinylidene chromophore bound within the interior of rhodopsin is coupled to transmembrane helix motion and receptor activation. The mechanism of activation that emerges is that multiple switches on the extracellular (or intradiscal) side of rhodopsin trigger structural changes that converge to disrupt the ionic lock between helices H3 and H6 on the intracellular side of the receptor.

A G protein-coupled receptor at work: the rhodopsin model

G protein-coupled receptors (GPCRs) are ubiquitous signal transducers in cell membranes, as well as important drug targets. Interaction with extracellular agonists turns the seven transmembrane helix (7TM) scaffold of a GPCR into a catalyst for GDP and GTP exchange in heterotrimeric Galphabetagamma proteins. Activation of the model GPCR, rhodopsin, is triggered by photoisomerization of its retinal ligand. From the augmentation of biochemical and biophysical studies by recent high-resolution 3D structures, its activation intermediates can now be interpreted as the stepwise engagement of protein domains. Rearrangement of TM5-TM6 opens a crevice at the cytoplasmic side of the receptor into which the C terminus of the Galpha subunit can bind. The Galpha C-terminal helix is used as a transmission rod to the nucleotide binding site. The mechanism relies on dynamic interactions between conserved residues and could therefore be common to other GPCRs.

NMR solution structure of human cannabinoid receptor-1 helix 7/8 peptide: candidate electrostatic interactions and microdomain formation

We report the NMR solution structure of a synthetic 40-mer (T(377)-E(416)) that encompasses human cannabinoid receptor-1 (hCB1) transmembrane helix 7 (TMH7) and helix 8 (H8) [hCB1(TMH7/H8)] in 30% trifluoroethanol/H(2)O. Structural features include, from the peptide's amino terminus, a hydrophobic alpha-helix (TMH7); a loop-like, 11 residue segment featuring a pronounced Pro-kink within the conserved NPxxY motif; a short amphipathic alpha-helix (H8) orthogonal to TMH7 with cationic and hydrophobic amino-acid clusters; and an unstructured C-terminal end. The hCB1(TMH7/H8) NMR solution structure suggests multiple electrostatic amino-acid interactions, including an intrahelical H8 salt bridge and a hydrogen-bond network involving the peptide's loop-like region. Potential cation-pi and cation-phenolic OH interactions between Y(397) in the TMH7 NPxxY motif and R(405) in H8 are identified as candidate structural forces promoting interhelical microdomain formation. This microdomain may function as a flexible molecular hinge during ligand-induced hCB1 conformer transitions.

Structural and dynamic effects of cholesterol at preferred sites of interaction with rhodopsin identified from microsecond length molecular dynamics simulations

An unresolved question about GPCR function is the role of membrane components in receptor stability and activation. In particular, cholesterol is known to affect the function of membrane proteins, but the details of its effect on GPCRs are still elusive. Here, we describe how cholesterol modulates the behavior of the TM1-TM2-TM7-helix 8(H8) functional network that comprises the highly conserved NPxxY(x)(5,6)F motif, through specific interactions with the receptor. The inferences are based on the analysis of microsecond length molecular dynamics (MD) simulations of rhodopsin in an explicit membrane environment. Three regions on the rhodopsin exhibit the highest cholesterol density throughout the trajectory: the extracellular end of TM7, a location resembling the high-density sterol area from the electron microscopy data; the intracellular parts of TM1, TM2, and TM4, a region suggested as the cholesterol binding site in the recent X-ray crystallography data on beta(2)-adrenergic GPCR; and the intracellular ends of TM2-TM3, a location that was categorized as the high cholesterol density area in multiple independent 100 ns MD simulations of the same system. We found that cholesterol primarily affects specific local perturbations of the helical TM domains such as the kinks in TM1, TM2, and TM7. These local distortions, in turn, relate to rigid-body motions of the TMs in the TM1-TM2-TM7-H8 bundle. The specificity of the effects stems from the nonuniform distribution of cholesterol around the protein. Through correlation analysis we connect local effects of cholesterol on structural perturbations with a regulatory role of cholesterol in the structural rearrangements involved in GPCR function.

A novel melanocortin-4 receptor gene mutation in a female patient with severe childhood obesity

This study targeted the identification of mutations of melanocortin-4 receptor gene (MC4R) in obese children. Fifty-one unrelated probands with early onset severe obesity (body mass index (BMI) > 99th percentile; 21 girls, mean age 10.6 +/- 3.6 years) were analyzed for nucleotide variations in the MC4R coding region, by the polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) method followed by direct DNA sequencing. MC4R variants were detected in three patients: the known I169S variant was found in heterozygote state in two patients and a novel heterozygous Y302F mutation was detected in one 12-year-old girl (BMI = 34 kg/m(2), BMI z-score 2.7) who has been overweight since the second year of life and suffered from hyperinsulinemia (at the age of 12: fasting insulin 45 mU/ml, after oral glucose load max. 300 mU/ml). The mutation also appears in the father, although both parents are obese (BMI father: 30.2 kg/m(2); mother: 31.9 kg/m(2)). This novel mutation is located in the functionally important NPXXY motif of the seventh transmembrane domain of the receptor. Functional characterization revealed reduction in cell surface expression and an alteration in signal transduction properties. These results add to the growing list of loss-of-function MC4R mutations in early onset obese patients and suggest an orexigenic effect of novel Y302F mutation.

Role of helix 8 in G protein-coupled receptors based on structure-function studies on the type 1 angiotensin receptor

G protein-coupled receptors (GPCRs) are transmembrane receptors that convert extracellular stimuli to intracellular signals. The type 1 angiotensin II receptor is a widely studied GPCR with roles in blood pressure regulation,water and salt balance and cell growth. The complex molecular and structural changes that underpin receptor activation and signaling are the focus of intense research. Increasingly, there is an appreciation that the plasma membrane participates in receptor function via direct, physical interactions that reciprocally modulate both lipid and receptor and provide microdomains for specialized activities. Reversible protein:lipid interactions are commonly mediated by amphipathic -helices in proteins and one such motif - a short helix, referred to as helix VIII/8 (H8), located at the start of the carboxyl (C)-terminus of GPCRs - is gaining recognition for its importance to GPCR function. Here, we review the identification of H8 in GPCRs and examine its capacity to sense and interact with diverse proteins and lipid environment, most notably with acidic lipids that include phosphatidylinositol phosphates.