Ligand binding and micro-switches in 7TM receptor structures

Feeding-dependent activation of enteric cells and sensory neurons by lymphatic fluid: evidence for a neurolymphocrine system

Lymphatic fluid is a plasma filtrate that can be viewed as having biological activity through the passive accumulation of molecules from the interstitial fluid. The possibility that lymphatic fluid is part of an active self-contained signaling process that parallels the endocrine system, through the activation of G-protein coupled receptors (GPCR), has remained unexplored. We show that the GPCR lysophosphatidic acid 5 (LPA5) is found in sensory nerve fibers expressing calcitonin gene-related peptide (CGRP) that innervate the lumen of lymphatic lacteals and enteric nerves. Using LPA5 as a model for nutrient-responsive GPCRs present on sensory nerves, we demonstrate that dietary protein hydrolysate (peptone) can induce c-Fos expression in enterocytes and nerves that express LPA5. Mesenteric lymphatic fluid (MLF) mobilizes intracellular calcium in cell models expressing LPA5 upon feeding in a time- and dose-dependent manner. Primary cultured neurons of the dorsal root ganglia expressing CGRP are activated by MLF, which is enhanced upon LPA5 overexpression. Activation is independent of the known LPA5 agonists, lysophosphatidic acid and farnesyl pyrophosphate. These data bring forth a pathway for the direct stimulation of sensory nerves by luminal contents and interstitial fluid. Thus, by activating LPA5 on sensory nerves, MLF provides a means for known and yet to be identified constituents of the interstitial fluid to act as signals to comprise a "neurolymphocrine" system.

Unraveling structural mechanisms of allosteric drug action

Orthosteric drugs block the active site to obstruct function; allosteric drugs modify the population of the active state, to modulate function. Available data lead us to propose that allosteric drugs can constitute anchors and drivers. The anchor docks into an allosteric pocket. The conformation with which it interacts is unchanged during the transition between the inactive and active states. The anchor provides the foundation that allows the driver to exert a 'pull' and/or 'push' action that shifts the receptor population from the inactive to the active state. The presence or absence of driver atom in an allosteric drug can exert opposite agonism. We map a strategy for driver identification and expect the allosteric trigger concept to transform agonist/antagonist drug discovery.

Retinal ligand mobility explains internal hydration and reconciles active rhodopsin structures

Rhodopsin, the mammalian dim-light receptor, is one of the best-characterized G-protein-coupled receptors, a pharmaceutically important class of membrane proteins that has garnered a great deal of attention because of the recent availability of structural information. Yet the mechanism of rhodopsin activation is not fully understood. Here, we use microsecond-scale all-atom molecular dynamics simulations, validated by solid-state (2)H nuclear magnetic resonance spectroscopy, to understand the transition between the dark and metarhodopsin I (Meta I) states. Our analysis of these simulations reveals striking differences in ligand flexibility between the two states. Retinal is much more dynamic in Meta I, adopting an elongated conformation similar to that seen in the recent activelike crystal structures. Surprisingly, this elongation corresponds to both a dramatic influx of bulk water into the hydrophobic core of the protein and a concerted transition in the highly conserved Trp265(6.48) residue. In addition, enhanced ligand flexibility upon light activation provides an explanation for the different retinal orientations observed in X-ray crystal structures of active rhodopsin.

Constitutively active rhodopsin and retinal disease

Rhodopsin is the light receptor in rod photoreceptor cells of the retina that initiates scotopic vision. In the dark, rhodopsin is bound to the chromophore 11-cis retinal, which locks the receptor in an inactive state. The maintenance of an inactive rhodopsin in the dark is critical for rod photoreceptor cells to remain highly sensitive. Perturbations by mutation or the absence of 11-cis retinal can cause rhodopsin to become constitutively active, which leads to the desensitization of photoreceptor cells and, in some instances, retinal degeneration. Constitutive activity can arise in rhodopsin by various mechanisms and can cause a variety of inherited retinal diseases including Leber congenital amaurosis, congenital night blindness, and retinitis pigmentosa. In this review, the molecular and structural properties of different constitutively active forms of rhodopsin are overviewed, and the possibility that constitutive activity can arise from different active-state conformations is discussed.

Lipid bilayer modules as determinants of K+ channel gating

The crystal structure of the sensorless pore module of a voltage-gated K(+) (Kv) channel showed that lipids occupy a crevice between subunits. We asked if individual lipid monolayers of the bilayer embody independent modules linked to channel gating modulation. Functional studies using single channel current recordings of the sensorless pore module reconstituted in symmetric and asymmetric lipid bilayers allowed us to establish the deterministic role of lipid headgroup on gating. We discovered that individual monolayers with headgroups that coat the bilayer-aqueous interface with hydroxyls stabilize the channel open conformation. The hydroxyl need not be at a terminal position and the effect is not dependent on the presence of phosphate or net charge on the lipid headgroup. Asymmetric lipid bilayers allowed us to determine that phosphoglycerides with glycerol or inositol on the extracellular facing monolayer stabilize the open conformation of the channel. This indirect effect is attributed to a change in water structure at the membrane interface. By contrast, inclusion of the positively charged lysyl-dioleoyl-phosphatidylglycerol exclusively on the cytoplasmic facing monolayer of the bilayer increases drastically the probability of finding the channel open. Such modulation is mediated by a π-cation interaction between Phe-19 of the pore module and the lysyl moiety anchored to the phosphatidylglycerol headgroup. The new findings imply that the specific chemistry of the lipid headgroup and its selective location in either monolayer of the bilayer dictate the stability of the open conformation of a Kv pore module in the absence of voltage-sensing modules.

Mapping the intramolecular signal transduction of G-protein coupled receptors

G-protein coupled receptors (GPCRs), a major gatekeeper of extracellular signals on plasma membrane, are unarguably one of the most important therapeutic targets. Given the recent discoveries of allosteric modulations, an allosteric wiring diagram of intramolecular signal transductions would be of great use to glean the mechanism of receptor regulation. Here, by evaluating betweenness centrality (CB ) of each residue, we calculate maps of information flow in GPCRs and identify key residues for signal transductions and their pathways. Compared with preexisting approaches, the allosteric hotspots that our CB -based analysis detects for A2 A adenosine receptor (A2 A AR) and bovine rhodopsin are better correlated with biochemical data. In particular, our analysis outperforms other methods in locating the rotameric microswitches, which are generally deemed critical for mediating orthosteric signaling in class A GPCRs. For A2 A AR, the inter-residue cross-correlation map, calculated using equilibrium structural ensemble from molecular dynamics simulations, reveals that strong signals of long-range transmembrane communications exist only in the agonist-bound state. A seemingly subtle variation in structure, found in different GPCR subtypes or imparted by agonist bindings or a point mutation at an allosteric site, can lead to a drastic difference in the map of signaling pathways and protein activity. The signaling map of GPCRs provides valuable insights into allosteric modulations as well as reliable identifications of orthosteric signaling pathways.

Recent structural advances of β1 and β2 adrenoceptors yield keys for ligand recognition and drug design

Because they represent attractive drug targets, adrenoceptors have been widely studied. Recent progress in structural data of β-adrenoceptors allows us to understand and predict key interactions in ligand recognition and receptor activation. Nevertheless, an important aspect of this process has only begun to be explored: the stabilization of a conformational state of these receptors upon contact with a ligand and the capacity of a ligand to influence receptor conformation through allosteric modulation, biased signaling, and selectivity. The aim of the present Perspective is to identify the well-defined orthosteric binding site and possible allosteric sites and to analyze the importance of the ligand-receptor interaction in the stabilization of certain receptor conformations. For this purpose, we have reviewed recent advances made through the use of X-ray data from ligand-β-adrenoceptor (including ADRB1 and ADRB2) crystal structures. Most importantly, implications in the medicinal chemistry field are explored in relation to drug design.

Novel insights on thyroid-stimulating hormone receptor signal transduction

The TSH receptor (TSHR) is a member of the glycoprotein hormone receptors, a subfamily of family A G protein-coupled receptors. The TSHR is of great importance for the growth and function of the thyroid gland. The TSHR and its endogenous ligand TSH are pivotal proteins with respect to a variety of physiological functions and malfunctions. The molecular events of TSHR regulation can be summarized as a process of signal transduction, including signal reception, conversion, and amplification. The steps during signal transduction from the extra- to the intracellular sites of the cell are not yet comprehensively understood. However, essential new insights have been achieved in recent years on the interrelated mechanisms at the extracellular region, the transmembrane domain, and intracellular components. This review contains a critical summary of available knowledge of the molecular mechanisms of signal transduction at the TSHR, for example, the key amino acids involved in hormone binding or in the structural conformational changes that lead to G protein activation or signaling regulation. Aspects of TSHR oligomerization, signaling promiscuity, signaling selectivity, phenotypes of genetic variations, and potential extrathyroidal receptor activity are also considered, because these are relevant to an understanding of the overall function of the TSHR, including physiological, pathophysiological, and pharmacological perspectives. Directions for future research are discussed.

GPCR activation: protonation and membrane potential

GPCR proteins represent the largest family of signaling membrane proteins in eukaryotic cells. Their importance to basic cell biology, human diseases, and pharmaceutical interventions is well established. Many crystal structures of GPCR proteins have been reported in both active and inactive conformations. These data indicate that agonist binding alone is not sufficient to trigger the conformational change of GPCRs necessary for binding of downstream G-proteins, yet other essential factors remain elusive. Based on analysis of available GPCR crystal structures, we identified a potential conformational switch around the conserved Asp2.50, which consistently shows distinct conformations between inactive and active states. Combining the structural information with the current literature, we propose an energy-coupling mechanism, in which the interaction between a charge change of the GPCR protein and the membrane potential of the living cell plays a key role for GPCR activation.

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.

Atypical chemokine receptors: from silence to sound

ACRs (atypical chemokine receptors) were initially referred to as 'silent' receptors on the basis of a lack of signalling and functional activities that are typically observed with conventional chemokine receptors. Although ACRs do not directly induce cell migration, they indirectly control leucocyte recruitment by shaping chemokine gradients in tissues through degradation, transcytosis or local concentration of their cognate ligands. Recent evidence also suggests that these biological activities are supported by G-protein-independent, β-arrestin-dependent signalling events. In the present article, we review current knowledge on structural and signalling properties of ACRs that are changing our view on this entire class of receptors from silent to endogenous β-arrestin-biased signalling receptors.

Retinal conformation governs pKa of protonated Schiff base in rhodopsin activation

We have explored the relationship between conformational energetics and the protonation state of the Schiff base in retinal, the covalently bound ligand responsible for activating the G protein-coupled receptor rhodopsin, using quantum chemical calculations. Guided by experimental structural determinations and large-scale molecular simulations on this system, we examined rotation about each bond in the retinal polyene chain, for both the protonated and deprotonated states that represent the dark and photoactivated states, respectively. Particular attention was paid to the torsional degrees of freedom that determine the shape of the molecule, and hence its interactions with the protein binding pocket. While most torsional degrees of freedom in retinal are characterized by large energetic barriers that minimize structural fluctuations under physiological temperatures, the C6-C7 dihedral defining the relative orientation of the β-ionone ring to the polyene chain has both modest barrier heights and a torsional energy surface that changes dramatically with protonation of the Schiff base. This surprising coupling between conformational degrees of freedom and protonation state is further quantified by calculations of the pKa as a function of the C6-C7 dihedral angle. Notably, pKa shifts of greater than two units arise from torsional fluctuations observed in molecular dynamics simulations of the full ligand-protein-membrane system. It follows that fluctuations in the protonation state of the Schiff base occur prior to forming the activated MII state. These new results shed light on important mechanistic aspects of retinal conformational changes that are involved in the activation of rhodopsin in the visual process.

Active-state models of ternary GPCR complexes: determinants of selective receptor-G-protein coupling

Based on the recently described crystal structure of the β₂ adrenergic receptor - G_s-protein complex, we report the first molecular-dynamics simulations of ternary GPCR complexes designed to identify the selectivity determinants for receptor-G-protein binding. Long-term molecular dynamics simulations of agonist-bound β2AR-Gα_s and D2R-Gα_i complexes embedded in a hydrated bilayer environment and computational alanine-scanning mutagenesis identified distinct residues of the N-terminal region of intracellular loop 3 to be crucial for coupling selectivity. Within the G-protein, specific amino acids of the α5-helix, the C-terminus of the Gα-subunit and the regions around αN-β1 and α4-β6 were found to determine receptor recognition. Knowledge of these determinants of receptor-G-protein binding selectivity is essential for designing drugs that target specific receptor/G-protein combinations.

Structural features for functional selectivity at serotonin receptors

Drugs active at G protein-coupled receptors (GPCRs) can differentially modulate either canonical or noncanonical signaling pathways via a phenomenon known as functional selectivity or biased signaling. We report biochemical studies showing that the hallucinogen lysergic acid diethylamide, its precursor ergotamine (ERG), and related ergolines display strong functional selectivity for β-arrestin signaling at the 5-HT2B 5-hydroxytryptamine (5-HT) receptor, whereas they are relatively unbiased at the 5-HT1B receptor. To investigate the structural basis for biased signaling, we determined the crystal structure of the human 5-HT2B receptor bound to ERG and compared it with the 5-HT1B/ERG structure. Given the relatively poor understanding of GPCR structure and function to date, insight into different GPCR signaling pathways is important to better understand both adverse and favorable therapeutic activities.

The role of Cysteine 6.47 in class A GPCRs

Background

The CWxP motif of transmembrane helix 6 (x: any residue) is highly conserved in class A GPCRs. Within this motif, W6.48 is a big star in the theory of the global “toggle switch” because of its key role in the activation mechanism of GPCRs upon ligand binding. With all footlights focused on W6.48, the reason why the preceding residue, C6.47, is largely conserved is still unknown. The present study is aimed to fill up this lack of knowledge by characterizing the role of C6.47 of the CWxP motif.

Results

A complete analysis of available crystal structures has been made alongside with molecular dynamics simulations of model peptides to explore a possible structural role for C6.47.

Conclusions

We conclude that C6.47 does not modulate the conformation of the TM6 proline kink and propose that C6.47 participates in the rearrangement of the TM6 and TM7 interface accompanying activation.

Biased and constitutive signaling in the CC-chemokine receptor CCR5 by manipulating the interface between transmembrane helices 6 and 7

The equilibrium state of CCR5 is manipulated here toward either activation or inactivation by introduction of single amino acid substitutions in the transmembrane domains (TMs) 6 and 7. Insertion of a steric hindrance mutation in the center of TM7 (G286F in position VII:09/7.42) resulted in biased signaling. Thus, β-arrestin recruitment was eliminated, whereas constitutive activity was observed in Gαi-mediated signaling. Furthermore, the CCR5 antagonist aplaviroc was converted to a full agonist (a so-called efficacy switch). Computational modeling revealed that the position of the 7TM receptor-conserved Trp in TM6 (Trp-248 in position VI:13/6.48, part of the CWXP motif) was influenced by the G286F mutation, causing Trp-248 to change orientation away from TM7. The essential role of Trp-248 in CCR5 activation was supported by complete inactivity of W248A-CCR5 despite maintaining chemokine binding. Furthermore, replacing Trp-248 with a smaller aromatic amino acid (Tyr/Phe) impaired the β-arrestin recruitment, yet with maintained G protein activity (biased signaling); also, here aplaviroc switched to a full agonist. Thus, the altered positioning of Trp-248, induced by G286F, led to a constraint of G protein active, but β-arrestin inactive and thus biased, CCR5 conformation. These results provide important information on the molecular interplay and impact of TM6 and TM7 for CCR5 activity, which may be extrapolated to other chemokine receptors and possibly to other 7TM receptors.

G-protein-coupled receptor structure, ligand binding and activation as studied by solid-state NMR spectroscopy

GPCRs (G-protein-coupled receptors) are versatile signalling molecules at the cell surface and make up the largest and most diverse family of membrane receptors in the human genome. They convert a large variety of extracellular stimuli into intracellular responses through the activation of heterotrimeric G-proteins, which make them key regulatory elements in a broad range of normal and pathological processes, and are therefore one of the most important targets for pharmaceutical drug discovery. Knowledge of a GPCR structure enables us to gain a mechanistic insight into its function and dynamics, and further aid rational drug design. Despite intensive research carried out over the last three decades, resolving the structural basis of GPCR function is still a major activity. The crystal structures obtained in the last 5 years provide the first opportunity to understand how protein structure dictates the unique functional properties of these complex signalling molecules. However, owing to the intrinsic hydrophobicity, flexibility and instability of membrane proteins, it is still a challenge to crystallize GPCRs, and, when this is possible, it is no longer in its native membrane environment and no longer without modification. Furthermore, the conformational change of the transmembrane α-helices associated with the structure activation increases the difficulty of capturing the activation state of a GPCR to a higher resolution by X-ray crystallography. On the other hand, solid-state NMR may offer a unique opportunity to study membrane protein structure, ligand binding and activation at atomic resolution in the native membrane environment, as well as described functionally significant dynamics. In the present review, we discuss some recent achievements of solid-state NMR for understanding GPCRs, the largest mammalian proteome at ~1% of the total expressed proteins. Structural information, details of determination, details of ligand conformations and the consequences of ligand binding to initiate activation can all be explored with solid-state NMR.

Constitutively active CCR5 chemokine receptors differ in mediating HIV envelope-dependent fusion

The CCR5 chemokine receptor is a rhodopsin-like G protein-coupled receptor that mediates the effects of pro-inflammatory β-chemokines. CCR5 is also the major co-receptor for entry of human immunodeficiency virus (HIV) into human cells. G protein-coupled receptors exist in ensembles of active and inactive conformations. Active receptor conformations can be stabilized by mutations. Although binding of the HIV envelope protein to CCR5 stimulates cellular signaling, the CCR5 conformation that induces fusion of the viral membrane with cellular membranes is not known. We mutated conserved amino acids to generate constitutively active CCR5 receptors, which are stabilized in active conformations, and tested the ability of constitutively active CCR5 receptors to mediate HIV envelope-directed membrane fusion. Mutation of the Asp^3.49(125) and Arg^6.32(225) residues of CCR5 did not cause constitutive activity, but Lys or Pro substitutions for Thr^2.56(82), in the TxP motif, caused high basal inositol phosphate signaling. Signaling did not increase in response to MIP-1β, suggesting that the Thr^2.56(82) mutants were fully stabilized in active conformations. The Thr^2.56(82)Lys mutation severely decreased cell surface CCR5 expression. Combining the Thr^2.56(82)Lys mutation with an Arg^6.32(225)Gln mutation partially reversed the decrease in expression. Mutants with Thr^2.56(82)Lys substitutions were poor mediators of HIV envelope-directed membrane fusion, but mutants with the Thr^2.65(82)Pro substitution exhibited full co-receptor function. Our results suggest that the Thr^2.65(82)Lys and Thr^2.65(82)Pro mutations stabilize distinct constitutively active CCR5 conformations. Lys in position 2.65(82) stabilizes activated receptor conformations that appear to be constitutively internalized and do not induce envelope-dependent membrane fusion, whereas Pro stabilizes activated conformations that are not constitutively internalized and fully mediate envelope-directed membrane fusion.

β-Arrestin recruitment and G protein signaling by the atypical human chemokine decoy receptor CCX-CKR

Chemokine receptors form a large subfamily of G protein-coupled receptors that predominantly activate heterotrimeric Gi proteins and are involved in immune cell migration. CCX-CKR is an atypical chemokine receptor with high affinity for CCL19, CCL21, and CCL25 chemokines, but is not known to activate intracellular signaling pathways. However, CCX-CKR acts as decoy receptor and efficiently internalizes these chemokines, thereby preventing their interaction with other chemokine receptors, like CCR7 and CCR9. Internalization of fluorescently labeled CCL19 correlated with β-arrestin2-GFP translocation. Moreover, recruitment of β-arrestins to CCX-CKR in response to CCL19, CCL21, and CCL25 was demonstrated using enzyme-fragment complementation and bioluminescence resonance energy transfer methods. To unravel why CCX-CKR is unable to activate Gi signaling, CCX-CKR chimeras were constructed by substituting its intracellular loops with the corresponding CCR7 or CCR9 domains. The signaling properties of chimeric CCX-CKR receptors were characterized using a cAMP-responsive element (CRE)-driven reporter gene assay. Unexpectedly, wild type CCX-CKR and a subset of the chimeras induced an increase in CRE activity in response to CCL19, CCL21, and CCL25 in the presence of the Gi inhibitor pertussis toxin. CCX-CKR signaling to CRE required an intact DRY motif. These data suggest that inactive Gi proteins impair CCX-CKR signaling most likely by hindering the interaction of this receptor with pertussis toxin-insensitive G proteins that transduce signaling to CRE. On the other hand, recruitment of the putative signaling scaffold β-arrestin to CCX-CKR in response to chemokines might allow activation of yet to be identified signal transduction pathways.

G-protein coupled receptor 83 (GPR83) signaling determined by constitutive and zinc(II)-induced activity

The G-protein coupled receptor 83 (GPR83) is an orphan G-protein coupled receptor for which the natural ligand(s) and signaling pathway(s) remain to be identified. Previous studies suggest a role of GPR83 in the regulation of thermogenesis and the control of circulating adiponectin. The aim of this study was to gain insights into the molecular underpinnings underlying GPR83 signaling. In particular, we aimed to assess the underlying G-protein activated signaling pathway of GPR83 and how this pathway is affected by mutational activation and zinc(II) challenge. Finally, we assessed the capacity of GPR83 for homodimerization. Our results show for the first time that mouse (m) GPR83 has high basal Gq/11 activity without affecting Gi or Gs signaling. Furthermore, we found that, under physiological conditions, zinc(II) (but not calcium(II) and magnesium(II)) potently activates mGPR83, thus identifying zinc(II) as an endogenous molecule with agonistic capability to activate mGPR83. In line with the observation that zinc(II)-ions activate mGPR83, we identified a cluster of ion-binding sensitive amino acids (e.g. His145, His204, Cys207, Glu217) in an activation sensitive receptor region of mGPR83. The occurrence of a constitutive activating mutant and a zinc(II)-binding residue at the N-terminal part corroborate the importance of this region in mGPR83 signal regulation. Finally, our results indicate that mGPR83 forms homodimers, which extend the current knowledge and molecular facets of GPR83 signaling.

Advances in methods to characterize ligand-induced ionic lock and rotamer toggle molecular switch in G protein-coupled receptors

Structural biology of GPCRs has made significant progress upon recently developed technologies for GPCRs expression/purification and elucidation of GPCRs crystal structures. The crystal structures provide a snapshot of the receptor structural disposition of GPCRs itself or with cocrystallized ligands, and the results are congruent with biophysical and computer modeling studies reported about GPCRs conformational and dynamics flexibility, regulated activation, and the various stabilizing interactions, such as "molecular switches." The molecular switches generally constitute the most conserved domains within a particular GPCR superfamily. Often agonist-induced receptor activation proceeds by the disruption of majority of these interactions, while antagonist and inverse agonist act as blockers and structural stabilizers, respectively. Several elegant studies, particularly for the β2AR, have demonstrated the relationship between ligand structure, receptor conformational changes, and corresponding pharmacological outcomes. Thus, it is of great importance to understand GPCRs activation related to cell signaling pathways. Herein, we summarize the steps to produce functional GPCRs, generate suitably fluorescent labeled GPCRs and the procedure to use that to understand if ligand-induced activation can proceed by activation of the GPCRs via ionic lock switch and/or rotamer toggle switch mechanisms. Such understanding of ligand structure and mechanism of receptor activation will provide great insight toward uncovering newer pathways of GPCR activation and aid in structure-based drug design.

Receptor clustering affects signal transduction at the membrane level in the reaction-limited regime

Many types of membrane receptors are found to be organized as clusters on the cell surface. We investigate the potential effect of such receptor clustering on the intracellular signal transduction stage. We consider a canonical pathway with a membrane receptor (R) activating a membrane-bound intracellular relay protein (G). We use Monte Carlo simulations to recreate biochemical reactions using different receptor spatial distributions and explore the dynamics of the signal transduction. Results show that activation of G by R is severely impaired by R clustering, leading to an apparent blunted biological effect compared to control. Paradoxically, this clustering decreases the half maximal effective dose (ED50) of the transduction stage, increasing the apparent affinity. We study an example of inter-receptor interaction in order to account for possible compensatory effects of clustering and observe the parameter range in which such interactions slightly counterbalance the loss of activation of G. The membrane receptors' spatial distribution affects the internal stages of signal amplification, suggesting a functional role for membrane domains and receptor clustering independently of proximity-induced receptor-receptor interactions.

Identifying ligand binding conformations of the β2-adrenergic receptor by using its agonists as computational probes

Recently available G-protein coupled receptor (GPCR) structures and biophysical studies suggest that the difference between the effects of various agonists and antagonists cannot be explained by single structures alone, but rather that the conformational ensembles of the proteins need to be considered. Here we use an elastic network model-guided molecular dynamics simulation protocol to generate an ensemble of conformers of a prototypical GPCR, β(2)-adrenergic receptor (β(2)AR). The resulting conformers are clustered into groups based on the conformations of the ligand binding site, and distinct conformers from each group are assessed for their binding to known agonists of β(2)AR. We show that the select ligands bind preferentially to different predicted conformers of β(2)AR, and identify a role of β(2)AR extracellular region as an allosteric binding site for larger drugs such as salmeterol. Thus, drugs and ligands can be used as "computational probes" to systematically identify protein conformers with likely biological significance.

Critical hydrogen bond formation for activation of the angiotensin II type 1 receptor

G protein-coupled receptors contain selectively important residues that play central roles in the conformational changes that occur during receptor activation. Asparagine 111 (N111(3.35)) is such a residue within the angiotensin II type 1 (AT(1)) receptor. Substitution of N111(3.35) for glycine leads to a constitutively active receptor, whereas substitution for tryptophan leads to an inactivable receptor. Here, we analyzed the AT(1) receptor and two mutants (N111G and N111W) by molecular dynamics simulations, which revealed a novel molecular switch involving the strictly conserved residue D74(2.50). Indeed, D74(2.50) forms a stable hydrogen bond (H-bond) with the residue in position 111(3.35) in the wild-type and the inactivable receptor. However, in the constitutively active mutant N111G-AT(1) receptor, residue D74 is reoriented to form a new H-bond with another strictly conserved residue, N46(1.50). When expressed in HEK293 cells, the mutant N46G-AT(1) receptor was poorly activable, although it retained a high binding affinity. Interestingly, the mutant N46G/N111G-AT(1) receptor was also inactivable. Molecular dynamics simulations also revealed the presence of a cluster of hydrophobic residues from transmembrane domains 2, 3, and 7 that appears to stabilize the inactive form of the receptor. Whereas this hydrophobic cluster and the H-bond between D74(2.50) and W111(3.35) are more stable in the inactivable N111W-AT(1) receptor, the mutant N111W/F77A-AT(1) receptor, designed to weaken the hydrophobic core, showed significant agonist-induced signaling. These results support the potential for the formation of an H-bond between residues D74(2.50) and N46(1.50) in the activation of the AT(1) receptor.

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

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

Structure-function of the G protein-coupled receptor superfamily

During the past few years, crystallography of G protein-coupled receptors (GPCRs) has experienced exponential growth, resulting in the determination of the structures of 16 distinct receptors-9 of them in 2012 alone. Including closely related subtype homology models, this coverage amounts to approximately 12% of the human GPCR superfamily. The adrenergic, rhodopsin, and adenosine receptor systems are also described by agonist-bound active-state structures, including a structure of the receptor-G protein complex for the β(2)-adrenergic receptor. Biochemical and biophysical techniques, such as nuclear magnetic resonance and hydrogen-deuterium exchange coupled with mass spectrometry, are providing complementary insights into ligand-dependent dynamic equilibrium between different functional states. Additional details revealed by high-resolution structures illustrate the receptors as allosteric machines that are controlled not only by ligands but also by ions, lipids, cholesterol, and water. This wealth of data is helping redefine our knowledge of how GPCRs recognize such a diverse array of ligands and how they transmit signals 30 angstroms across the cell membrane; it also is shedding light on a structural basis of GPCR allosteric modulation and biased signaling.

PheVI:09 (Phe6.44) as a sliding microswitch in seven-transmembrane (7TM) G protein-coupled receptor activation

In seven-transmembrane (7TM), G protein-coupled receptors, highly conserved residues function as microswitches, which alternate between different conformations and interaction partners in an extended allosteric interface between the transmembrane segments performing the large scale conformational changes upon receptor activation. Computational analysis using x-ray structures of the β(2)-adrenergic receptor demonstrated that PheVI:09 (6.44), which in the inactive state is locked between the backbone and two hydrophobic residues in transmembrane (TM)-III, upon activation slides ∼2 Å toward TM-V into a tight pocket generated by five hydrophobic residues protruding from TM-III and TM-V. Of these, the residue in position III:16 (3.40) (often an Ile or Val) appears to function as a barrier or gate for the transition between inactive and active conformation. Mutational analysis showed that PheVI:09 is essential for the constitutive and/or agonist-induced signaling of the ghrelin receptor, GPR119, the β(2)-adrenergic receptor, and the neurokinin-1 receptor. Substitution of the residues constituting the hydrophobic pocket between TM-III and TM-V in the ghrelin receptor in four of five positions impaired receptor signaling. In GPR39, representing the 12% of 7TM receptors lacking an aromatic residue at position VI:09, unchanged agonist-induced signaling was observed upon Ala substitution of LeuVI:09 despite reduced cell surface expression of the mutant receptor. It is concluded that PheVI:09 constitutes an aromatic microswitch that stabilizes the active, outward tilted conformation of TM-VI relative to TM-III by sliding into a tight hydrophobic pocket between TM-III and TM-V and that the hydrophobic residue in position III:16 constitutes a gate for this transition.

A homology modeling study toward the understanding of three-dimensional structure and putative pharmacological profile of the G-protein coupled receptor GPR55

The orphan G-protein coupled receptor GPR55 was shown to bind to certain cannabinoid compounds which led to its initial classification as the third type of cannabinoid receptor. Later studies showed that lysophosphatidylinositol (LPI) also activated GPR55, in particular 2-arachidonoyl-LPI was proposed to be its endogenous ligand. However, the results of pharmacological studies regarding GPR55 have been quite inconsistent. Despite its contradictory pharmacological profile, GPR55 has been implicated in various disease states including inflammatory and neuropathic pain, metabolic bone diseases, and cancer. Herein, we report the ligand binding properties of GPR55 by applying homology modeling and automated docking algorithms in order to understand its pharmacological profile. The 3D homology model of GPR55 was built based on the adenosine A(2A) receptor crystal structure. Docking studies of several types of reported ligands were carried out afterwards. The results indicated that both hydrogen bonding and hydrophobic interactions contributed significantly for its ligand binding and the amino acid residue Lys80 seemed to be the anchor residue for receptor recognition. In addition, its putative agonist and antagonist appeared to recognize different domains of the receptor corresponding to their reported pharmacological activities.

Emerging opportunities for allosteric modulation of G-protein coupled receptors

Their ubiquitous nature, wide cellular distribution and versatile molecular recognition and signalling help make G-protein binding receptors (GPCRs) the most important class of membrane proteins in clinical medicine, accounting for ∼40% of all current therapeutics. A large percentage of current drugs target the endogenous ligand binding (orthosteric) site, which are structurally and evolutionarily conserved, particularly among members of the same GPCR subfamily. With the recent advances in GPCR X-ray crystallography, new opportunities for developing novel subtype selective drugs have emerged. Given the increasing recognition that the extracellular surface conformation changes in response to ligand binding, it is likely that all GPCRs possess an allosteric site(s) capable of regulating GPCR signalling. Allosteric sites are less structurally conserved than their corresponding orthosteric site and thus provide new opportunities for the development of more selective drugs. Constitutive oligomerisation (dimerisation) identified in many of the GPCRs investigated, adds another dimension to the structural and functional complexity of GPCRs. In this review, we compare 60 crystal structures of nine GPCR subtypes (rhodopsin, ß₂-AR, ß₁-AR, A(2a)-AR, CXCR4, D₃R, H₁R, M₂R, M₃R) across four subfamilies of Class A GPCRs, and discuss mechanisms involved in receptor activation and potential allosteric binding sites across the highly variable extracellular surface of these GPCRs. This analysis has identified a new extracellular salt bridge (ESB-2) that might be exploited in the design of allosteric modulators.

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.

Modeling of mammalian olfactory receptors and docking of odorants

Olfactory receptors (ORs) belong to the superfamily of G protein-coupled receptors (GPCRs), the second largest class of genes after those related to immunity, and account for about 3 % of mammalian genomes. ORs are present in all multicellular organisms and represent more than half the GPCRs in mammalian species (e.g., the mouse OR repertoire contains >1,000 functional genes). ORs are mainly expressed in the olfactory epithelium where they detect odorant molecules, but they are also expressed in a number of other cells, such as sperm cells, although their functions in these cells remain mostly unknown. It has recently been reported that ORs are present in tumoral tissues where they are expressed at different levels than in healthy tissues. A specific OR is over-expressed in prostate cancer cells, and activation of this OR has been shown to inhibit the proliferation of these cells. Odorant stimulation of some of these receptors results in inhibition of cell proliferation. Even though their biological role has not yet been elucidated, these receptors might constitute new targets for diagnosis and therapeutics. It is important to understand the activation mechanism of these receptors at the molecular level, in particular to be able to predict which ligands are likely to activate a particular receptor ('deorphanization') or to design antagonists for a given receptor. In this review, we describe the in silico methodologies used to model the three-dimensional (3D) structure of ORs (in the more general framework of GPCR modeling) and to dock ligands into these 3D structures.

Reversed binding of a small molecule ligand in homologous chemokine receptors - differential role of extracellular loop 2

BACKGROUND AND PURPOSE

The majority of small molecule compounds targeting chemokine receptors share a similar pharmacophore with a centrally located aliphatic positive charge and flanking aromatic moieties. Here we describe a novel piperidine-based compound with structural similarity to previously described CCR8-specific agonists, but containing a unique phenyl-tetrazol moiety which, in addition to activity at CCR8 was also active at CCR1.

EXPERIMENTAL APPROACH

Single point mutations were introduced in CCR1 and CCR8, and their effect on small molecule ligand-induced receptor activation was examined through inositol trisphosphate (IP(3) ) accumulation. The molecular interaction profile of the agonist was verified by molecular modeling.

KEY RESULTS

The chemokine receptor conserved glutamic acid in TM-VII served as a common anchor for the positively charged amine in the piperidine ring. However, whereas the phenyl-tetrazol group interacted with TyrIV:24 (Tyr(172) ) and TyrIII:09 (Tyr(114) ) in the major binding pocket (delimited by TM-III to VII) of CCR8, it also interacted with TrpII:20 (Trp(90) ) and LysII:24 (Lys(94) ) in the minor counterpart (delimited TM-I to III, plus TM-VII) in CCR1. A straightening of TM-II by Ala-substitution of ProII:18 confirmed its unique role in CCR1. The extracellular loop 2 (ECL-2) contributed directly to the small molecule binding site in CCR1, whereas it contributed to efficacy, but not potency in CCR8.

CONCLUSION AND IMPLICATIONS

Despite high ligand potency and efficacy and receptor similarity, this dual-active and bitopic compound binds oppositely in CCR1 and CCR8 with different roles of ECL-2, thereby expanding and diversifying the influence of extracellular receptor regions in drug action.

Comparison of the activation kinetics of the M3 acetylcholine receptor and a constitutively active mutant receptor in living cells

Activation of G-protein-coupled receptors is the first step of the signaling cascade triggered by binding of an agonist. Here we compare the activation kinetics of the G(q)-coupled M(3) acetylcholine receptor (M(3)-AChR) with that of a constitutively active mutant receptor (M(3)-AChR-N514Y) using M(3)-AChR constructs that report receptor activation by changes in the fluorescence resonance energy transfer (FRET) signal. We observed a leftward shift in the concentration-dependent FRET response for acetylcholine and carbachol with M(3)-AChR-N514Y. Consistent with this result, at submaximal agonist concentrations, the activation kinetics of M(3)-AChR-N514Y were significantly faster, whereas at maximal agonist concentrations the kinetics of receptor activation were identical. Receptor deactivation was significantly faster with carbachol than with acetylcholine and was significantly delayed by the N514Y mutation. Receptor-G-protein interaction was measured by FRET between M(3)-AChR-yellow fluorescent protein (YFP) and cyan fluorescent protein (CFP)-Gγ(2). Agonist-induced receptor-G-protein coupling was of a time scale similar to that of receptor activation. As observed for receptor deactivation, receptor-G-protein dissociation was slower for acetylcholine than that for carbachol. Acetylcholine-stimulated increases in receptor-G-protein coupling of M(3)-AChR-N514Y reached only 12% of that of M(3)-AChR and thus cannot be kinetically analyzed. G-protein activation was measured using YFP-tagged Gα(q) and CFP-tagged Gγ(2). Activation of G(q) was significantly slower than receptor activation and indistinguishable for the two agonists. However, G(q) deactivation was significantly prolonged for acetylcholine compared with that for carbachol. Consistent with decreased agonist-stimulated coupling to G(q), agonist-stimulated G(q) activation by M(3)-AChR-N514Y was not detected. Taken together, these results indicate that the N514Y mutation produces constitutive activation of M(3)-AChR by decreasing the rate of receptor deactivation, while having minimal effect on receptor activation.

Modulation of constitutive activity and signaling bias of the ghrelin receptor by conformational constraint in the second extracellular loop

Based on a rare, natural Glu for Ala-204(C+6) variant located six residues after the conserved Cys residue in extracellular loop 2b (ECL2b) associated with selective elimination of the high constitutive signaling of the ghrelin receptor, this loop was subjected to a detailed structure functional analysis. Introduction of Glu in different positions demonstrated that although the constitutive signaling was partly reduced when introduced in position 205(C+7) it was only totally eliminated in position 204(C+6). No charge-charge interaction partner could be identified for the Glu(C+6) variant despite mutational analysis of a number of potential partners in the extracellular loops and outer parts of the transmembrane segments. Systematic probing of position 204(C+6) with amino acid residues of different physicochemical properties indicated that a positively charged Lys surprisingly provided phenotypes similar to those of the negatively charged Glu residue. Computational chemistry analysis indicated that the propensity for the C-terminal segment of extracellular loop 2b to form an extended α-helix was increased from 15% in the wild type to 89 and 82% by introduction in position 204(C+6) of a Glu or a Lys residue, respectively. Moreover, the constitutive activity of the receptor was inhibited by Zn(2+) binding in an engineered metal ion site, stabilizing an α-helical conformation of this loop segment. It is concluded that the high constitutive activity of the ghrelin receptor is dependent upon flexibility in the C-terminal segment of extracellular loop 2 and that mutations or ligand binding that constrains this segment and thereby conceivably the movements of transmembrane domain V relative to transmembrane domain III inhibits the high constitutive signaling.

Structural basis for allosteric regulation of GPCRs by sodium ions

Pharmacological responses of G protein-coupled receptors (GPCRs) can be fine-tuned by allosteric modulators. Structural studies of such effects have been limited due to the medium resolution of GPCR structures. We reengineered the human A(2A) adenosine receptor by replacing its third intracellular loop with apocytochrome b(562)RIL and solved the structure at 1.8 angstrom resolution. The high-resolution structure allowed us to identify 57 ordered water molecules inside the receptor comprising three major clusters. The central cluster harbors a putative sodium ion bound to the highly conserved aspartate residue Asp(2.50). Additionally, two cholesterols stabilize the conformation of helix VI, and one of 23 ordered lipids intercalates inside the ligand-binding pocket. These high-resolution details shed light on the potential role of structured water molecules, sodium ions, and lipids/cholesterol in GPCR stabilization and function.

Structural Diversity in Conserved Regions Like the DRY-Motif among Viral 7TM Receptors-A Consequence of Evolutionary Pressure?

Several herpes- and poxviruses have captured chemokine receptors from their hosts and modified these to their own benefit. The human and viral chemokine receptors belong to class A 7 transmembrane (TM) receptors which are characterized by several structural motifs like the DRY-motif in TM3 and the C-terminal tail. In the DRY-motif, the arginine residue serves important purposes by being directly involved in G protein coupling. Interestingly, among the viral receptors there is a greater diversity in the DRY-motif compared to their endogenous receptor homologous. The C-terminal receptor tail constitutes another regulatory region that through a number of phosphorylation sites is involved in signaling, desensitization, and internalization. Also this region is more variable among virus-encoded 7TM receptors compared to human class A receptors. In this review we will focus on these two structural motifs and discuss their role in viral 7TM receptor signaling compared to their endogenous counterparts.

The arginine of the DRY motif in transmembrane segment III functions as a balancing micro-switch in the activation of the β2-adrenergic receptor

Recent high resolution x-ray structures of the β2-adrenergic receptor confirmed a close salt-bridge interaction between the suspected micro-switch residue ArgIII:26 (Arg3.50) and the neighboring AspIII:25 (Asp3.49). However, neither the expected "ionic lock" interactions between ArgIII:26 and GluVI:-06 (Glu6.30) in the inactive conformation nor the interaction with TyrV:24 (Tyr5.58) in the active conformation were observed in the x-ray structures. Here we find through molecular dynamics simulations, after removal of the stabilizing T4 lysozyme, that the expected salt bridge between ArgIII:26 and GluVI:-06 does form relatively easily in the inactive receptor conformation. Moreover, mutational analysis of GluVI:-06 in TM-VI and the neighboring AspIII:25 in TM-III demonstrated that these two residues do function as locks for the inactive receptor conformation as we observed increased G(s) signaling, arrestin mobilization, and internalization upon alanine substitutions. Conversely, TyrV:24 appears to play a role in stabilizing the active receptor conformation as loss of function of G(s) signaling, arrestin mobilization, and receptor internalization was observed upon alanine substitution of TyrV:24. The loss of function of the TyrV:24 mutant could partly be rescued by alanine substitution of either AspIII:25 or GluVI:-06 in the double mutants. Surprisingly, removal of the side chain of the ArgIII:26 micro-switch itself had no effect on G(s) signaling and internalization and only reduced arrestin mobilization slightly. It is suggested that ArgIII:26 is equally important for stabilizing the inactive and the active conformation through interaction with key residues in TM-III, -V, and -VI, but that the ArgIII:26 micro-switch residue itself apparently is not essential for the actual G protein activation.

Evolutionary Analysis of Functional Divergence among Chemokine Receptors, Decoy Receptors, and Viral Receptors

Chemokine receptors (CKRs) function in the inflammatory response and in vertebrate homeostasis. Decoy and viral receptors are two types of CKR homologs with modified functions from those of the typical CKRs. The decoy receptors are able to bind ligands without signaling. On the other hand, the viral receptors show constitutive signaling without ligands. We examined the sites related to the functional difference. At first, the decoy and viral receptors were each classified into five groups, based on the molecular phylogenetic analysis. A multiple amino acid sequence alignment between each group and the CKRs was then constructed. The difference in the amino acid composition between the group and the CKRs was evaluated as the Kullback-Leibler (KL) information value at each alignment site. The KL information value is considered to reflect the difference in the functional constraints at the site. The sites with the top 5% of KL information values were selected and mapped on the structure of a CKR. The comparisons with decoy receptor groups revealed that the detected sites were biased on the intracellular side. In contrast, the sites detected from the comparisons with viral receptor groups were found on both the extracellular and intracellular sides. More sites were found in the ligand binding pocket in the analyses of the viral receptor groups, as compared to the decoy receptor groups. Some of the detected sites were located in the GPCR motifs. For example, the DRY motif of the decoy receptors was often degraded, although the motif of the viral receptors was basically conserved. The observations for the viral receptor groups suggested that the constraints in the pocket region are loose and that the sites on the intracellular side are different from those for the decoy receptors, which may be related to the constitutive signaling activity of the viral receptors.

Ligand-mimicking receptor variant discloses binding and activation mode of prolactin-releasing peptide

The prolactin-releasing peptide receptor and its bioactive RF-amide peptide (PrRP20) have been investigated to explore the ligand binding mode of peptide G-protein-coupled receptors (GPCRs). By receptor mutagenesis, we identified the conserved aspartate in the upper transmembrane helix 6 (Asp(6.59)) of the receptor as the first position that directly interacts with arginine 19 of the ligand (Arg(19)). Replacement of Asp(6.59) with Arg(19) of PrRP20 led to D6.59R, which turned out to be a constitutively active receptor mutant (CAM). This suggests that the mutated residue at the top of transmembrane helix 6 mimics Arg(19) by interacting with additional binding partners in the receptor. Next, we generated an initial comparative model of this CAM because no ligand docking was required, and we selected the next set of receptor mutants to find the engaged partners of the binding pocket. In an iterative process, we identified two acidic residues and two hydrophobic residues that form the peptide ligand binding pocket. As all residues are localized on top or in the upper part of the transmembrane domains, we clearly can show that the extracellular surface of the receptor is sufficient for full signal transduction for prolactin-releasing peptide, rather than a deep, membrane-embedded binding pocket. This contributes to the knowledge of the binding of peptide ligands to GPCRs and might facilitate the development of GPCR ligands, but it also provides new targeting of CAMs involved in hereditary diseases.

Automatic modeling of mammalian olfactory receptors and docking of odorants

We present a procedure that (i) automates the homology modeling of mammalian olfactory receptors (ORs) based on the six three-dimensional (3D) structures of G protein-coupled receptors (GPCRs) available so far and (ii) performs the docking of odorants on these models, using the concept of colony energy to score the complexes. ORs exhibit low-sequence similarities with other GPCR and current alignment methods often fail to provide a reliable alignment. Here, we use a fold recognition technique to obtain a robust initial alignment. We then apply our procedure to a human OR that we have previously functionally characterized. The analysis of the resulting in silico complexes, supported by receptor mutagenesis and functional assays in a heterologous expression system, suggests that antagonists dock in the upper part of the binding pocket whereas agonists dock in the narrow lower part. We propose that the potency of agonists in activating receptors depends on their ability to establish tight interactions with the floor of the binding pocket. We developed a web site that allows the user to upload a GPCR sequence, choose a ligand in a library and obtain the 3D structure of the free receptor and ligand-receptor complex (http://genome.jouy.inra.fr/GPCRautomodel).

Homology modeling and molecular dynamics simulations of the active state of the nociceptin receptor reveal new insights into agonist binding and activation

The opioid receptor-like receptor, also known as the nociceptin receptor (NOP), is a class A G protein-coupled receptor (GPCR) in the opioid receptor family. Although NOP shares a significant homology with the other opioid receptors, it does not bind known opioid ligands and has been shown to have a distinct mechanism of activation compared to the closely related opioid receptors mu, delta, and kappa. Previously reported homology models of the NOP receptor, based on the inactive-state GPCR crystal structures, give limited information on the activation and selectivity features of this fourth member of the opioid receptor family. We report here the first active-state homology model of the NOP receptor based on the opsin GPCR crystal structure. An inactive-state homology model of NOP was also built using a multiple template approach. Molecular dynamics simulation of the active-state NOP model and comparison to the inactive-state model suggest that NOP activation involves movements of transmembrane (TM)3 and TM6 and several activation microswitches, consistent with GPCR activation. Docking of the selective nonpeptidic NOP agonist ligand Ro 64-6198 into the active-state model reveals active-site residues in NOP that play a role in the high selectivity of this ligand for NOP over the other opioid receptors. Docking the shortest active fragment of endogenous agonist nociceptin/orphaninFQ (residues 1-13) shows that the NOP extracellular loop 2 (EL2) loop interacts with the positively charged residues (8-13) of N/OFQ. Both agonists show extensive polar interactions with residues at the extracellular end of the TM domain and EL2 loop, suggesting agonist-induced reorganization of polar networks, during receptor activation.

Molecular modeling of the M3 acetylcholine muscarinic receptor and its binding site

The present study reports the results of a combined computational and site mutagenesis study designed to provide new insights into the orthosteric binding site of the human M3 muscarinic acetylcholine receptor. For this purpose a three-dimensional structure of the receptor at atomic resolution was built by homology modeling, using the crystallographic structure of bovine rhodopsin as a template. Then, the antagonist N-methylscopolamine was docked in the model and subsequently embedded in a lipid bilayer for its refinement using molecular dynamics simulations. Two different lipid bilayer compositions were studied: one component palmitoyl-oleyl phosphatidylcholine (POPC) and two-component palmitoyl-oleyl phosphatidylcholine/palmitoyl-oleyl phosphatidylserine (POPC-POPS). Analysis of the results suggested that residues F222 and T235 may contribute to the ligand-receptor recognition. Accordingly, alanine mutants at positions 222 and 235 were constructed, expressed, and their binding properties determined. The results confirmed the role of these residues in modulating the binding affinity of the ligand.

Rare MTNR1B variants impairing melatonin receptor 1B function contribute to type 2 diabetes

Genome-wide association studies have revealed that common noncoding variants in MTNR1B (encoding melatonin receptor 1B, also known as MT(2)) increase type 2 diabetes (T2D) risk(1,2). Although the strongest association signal was highly significant (P < 1 × 10(-20)), its contribution to T2D risk was modest (odds ratio (OR) of ∼1.10-1.15)(1-3). We performed large-scale exon resequencing in 7,632 Europeans, including 2,186 individuals with T2D, and identified 40 nonsynonymous variants, including 36 very rare variants (minor allele frequency (MAF) <0.1%), associated with T2D (OR = 3.31, 95% confidence interval (CI) = 1.78-6.18; P = 1.64 × 10(-4)). A four-tiered functional investigation of all 40 mutants revealed that 14 were non-functional and rare (MAF < 1%), and 4 were very rare with complete loss of melatonin binding and signaling capabilities. Among the very rare variants, the partial- or total-loss-of-function variants but not the neutral ones contributed to T2D (OR = 5.67, CI = 2.17-14.82; P = 4.09 × 10(-4)). Genotyping the four complete loss-of-function variants in 11,854 additional individuals revealed their association with T2D risk (8,153 individuals with T2D and 10,100 controls; OR = 3.88, CI = 1.49-10.07; P = 5.37 × 10(-3)). This study establishes a firm functional link between MTNR1B and T2D risk.