Three-dimensional representations of G protein-coupled receptor structures and mechanisms

CovET: A covariation-evolutionary trace method that identifies protein structure-function modules

Measuring the relative effect that any two sequence positions have on each other may improve protein design or help better interpret coding variants. Current approaches use statistics and machine learning but rarely consider phylogenetic divergences which, as shown by Evolutionary Trace studies, provide insight into the functional impact of sequence perturbations. Here, we reframe covariation analyses in the Evolutionary Trace framework to measure the relative tolerance to perturbation of each residue pair during evolution. This approach (CovET) systematically accounts for phylogenetic divergences: at each divergence event, we penalize covariation patterns that belie evolutionary coupling. We find that while CovET approximates the performance of existing methods to predict individual structural contacts, it performs significantly better at finding structural clusters of coupled residues and ligand binding sites. For example, CovET found more functionally critical residues when we examined the RNA recognition motif and WW domains. It correlates better with large-scale epistasis screen data. In the dopamine D2 receptor, top CovET residue pairs recovered accurately the allosteric activation pathway characterized for Class A G protein-coupled receptors. These data suggest that CovET ranks highest the sequence position pairs that play critical functional roles through epistatic and allosteric interactions in evolutionarily relevant structure-function motifs. CovET complements current methods and may shed light on fundamental molecular mechanisms of protein structure and function.

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.

Evolutionary analyses reveal immune cell receptor GPR84 as a conserved receptor for bacteria-derived molecules

The G protein-coupled receptor 84 (GPR84) is found in immune cells and its expression is increased under inflammatory conditions. Activation of GPR84 by medium-chain fatty acids results in pro-inflammatory responses. Here, we screened available vertebrate genome data and found that GPR84 is present in vertebrates for more than 500 million years but absent in birds and a pseudogene in bats. Cloning and functional characterization of several mammalian GPR84 orthologs in combination with evolutionary and model-based structural analyses revealed evidence for positive selection of bear GPR84 orthologs. Naturally occurring human GPR84 variants are most frequent in Asian populations causing a loss of function. Further, we identified cis- and trans-2-decenoic acid, both known to mediate bacterial communication, as evolutionary highly conserved ligands. Our integrated set of approaches contributes to a comprehensive understanding of GPR84 in terms of evolutionary and structural aspects, highlighting GPR84 as a conserved immune cell receptor for bacteria-derived molecules.

Combined docking and machine learning identify key molecular determinants of ligand pharmacological activity on β2 adrenoceptor

G protein‐coupled receptors (GPCRs) are valuable therapeutic targets for many diseases. A central question of GPCR drug discovery is to understand what determines the agonism or antagonism of ligands that bind them. Ligands exert their action via the interactions in the ligand binding pocket. We hypothesized that there is a common set of receptor interactions made by ligands of diverse structures that mediate their action and that among a large dataset of different ligands, the functionally important interactions will be over‐represented. We computationally docked ~2700 known β2AR ligands to multiple β2AR structures, generating ca 75 000 docking poses and predicted all atomic interactions between the receptor and the ligand. We used machine learning (ML) techniques to identify specific interactions that correlate with the agonist or antagonist activity of these ligands. We demonstrate with the application of ML methods that it is possible to identify the key interactions associated with agonism or antagonism of ligands. The most representative interactions for agonist ligands involve K97^2.68×67, F194^ECL2, S203^5.42×43, S204^5.43×44, S207^5.46×641, H296^6.58×58, and K305^7.32×31. Meanwhile, the antagonist ligands made interactions with W286^6.48×48 and Y316^7.43×42, both residues considered to be important in GPCR activation. The interpretation of ML analysis in human understandable form allowed us to construct an exquisitely detailed structure‐activity relationship that identifies small changes to the ligands that invert their pharmacological activity and thus helps to guide the drug discovery process. This approach can be readily applied to any drug target.

Structures of the ADGRG2-Gs complex in apo and ligand-bound forms

Adhesion G protein-coupled receptors are elusive in terms of their structural information and ligands. Here, we solved the cryogenic-electron microscopy (cryo-EM) structure of apo-ADGRG2, an essential membrane receptor for maintaining male fertility, in complex with a G_s trimer. Whereas the formations of two kinks were determinants of the active state, identification of a potential ligand-binding pocket in ADGRG2 facilitated the screening and identification of dehydroepiandrosterone (DHEA), dehydroepiandrosterone sulfate and deoxycorticosterone as potential ligands of ADGRG2. The cryo-EM structures of DHEA-ADGRG2-G_s provided interaction details for DHEA within the seven transmembrane domains of ADGRG2. Collectively, our data provide a structural basis for the activation and signaling of ADGRG2, as well as characterization of steroid hormones as ADGRG2 ligands, which might be used as useful tools for further functional studies of the orphan ADGRG2.

Angiotensin and Endothelin Receptor Structures With Implications for Signaling Regulation and Pharmacological Targeting

In conjunction with the endothelin (ET) type A (ETAR) and type B (ETBR) receptors, angiotensin (AT) type 1 (AT1R) and type 2 (AT2R) receptors, are peptide-binding class A G-protein-coupled receptors (GPCRs) acting in a physiologically overlapping context. Angiotensin receptors (ATRs) are involved in regulating cell proliferation, as well as cardiovascular, renal, neurological, and endothelial functions. They are important therapeutic targets for several diseases or pathological conditions, such as hypertrophy, vascular inflammation, atherosclerosis, angiogenesis, and cancer. Endothelin receptors (ETRs) are expressed primarily in blood vessels, but also in the central nervous system or epithelial cells. They regulate blood pressure and cardiovascular homeostasis. Pathogenic conditions associated with ETR dysfunctions include cancer and pulmonary hypertension. While both receptor groups are activated by their respective peptide agonists, pathogenic autoantibodies (auto-Abs) can also activate the AT1R and ETAR accompanied by respective clinical conditions. To date, the exact mechanisms and differences in binding and receptor-activation mediated by auto-Abs as opposed to endogenous ligands are not well understood. Further, several questions regarding signaling regulation in these receptors remain open. In the last decade, several receptor structures in the apo- and ligand-bound states were determined with protein X-ray crystallography using conventional synchrotrons or X-ray Free-Electron Lasers (XFEL). These inactive and active complexes provide detailed information on ligand binding, signal induction or inhibition, as well as signal transduction, which is fundamental for understanding properties of different activity states. They are also supportive in the development of pharmacological strategies against dysfunctions at the receptors or in the associated signaling axis. Here, we summarize current structural information for the AT1R, AT2R, and ETBR to provide an improved molecular understanding.

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.

A new class of serotonin 5-HT2A /5-HT2C receptor inverse agonists: Synthesis, molecular modeling, in vitro and in vivo pharmacology of novel 2-aminotetralins

BACKGROUND AND PURPOSE

The 5-HT receptor (5-HTR) subtypes 5-HT_2A and 5-HT_2C are important neurotherapeutic targets, though, obtaining selectivity over 5-HT_2B and closely related histamine H₁ Rs is challenging. Here, we delineated molecular determinants of selective binding to 5-HT_2A and 5-HT_2C Rs for novel 4-phenyl-2-dimethylaminotetralins (4-PATs).

EXPERIMENTAL APPROACH

We synthesized 42 novel 4-PATs with halogen or aryl moieties at the C(4)-phenyl meta position. Affinity, function, molecular modeling, and 5-HT_2A R mutagenesis studies were undertaken to understand structure-activity relationships at 5-HT₂ -type and H₁ Rs. Lead 4-PAT-type selective 5-HT_2A /5-HT_2C R inverse agonists were compared to pimavanserin, a selective 5-HT_2A /5-HT_2C R inverse agonist approved to treat psychoses, in the mouse head twitch response, and locomotor activity assays, as models relevant to antipsychotic drug development.

KEY RESULTS

Most 4-PAT diastereomers in the (2S,4R)-configuration bound non-selectively to 5-HT_2A , 5-HT_2C, and H₁ Rs, with >100-fold selectivity over 5-HT_2B Rs, whereas, diastereomers in the (2R,4R)-configuration bound preferentially to 5-HT_2A over 5-HT_2C Rs and had >100-fold selectivity over 5-HT_2B and H₁ Rs. Results suggest that G238^5.42 and V235^5.39 in 5-HT_2A Rs (conserved in 5-HT_2C Rs) are important for high affinity binding, whereas, interactions with T194^5.42 and W158^4.56 determine H₁ R affinity. The 4-PAT (2S,4R)-2k, a potent and selective 5-HT_2A /5-HT_2C R inverse agonist, had activity like pimavanserin in the mouse head-twitch response assay, but was distinct in not suppressing locomotor activity.

CONCLUSIONS AND IMPLICATIONS

We provide evidence that the novel 4-PAT chemotype can yield selective 5-HT_2A /5-HT_2C R inverse agonists for antipsychotic drug development by optimizing ligand-receptor interactions in transmembrane domain 5. We also show that chirality can be exploited to attain selectivity over H₁ Rs which may circumvent sedative effects.

Structures of active melanocortin-4 receptor-Gs-protein complexes with NDP-α-MSH and setmelanotide

The melanocortin-4 receptor (MC4R), a hypothalamic master regulator of energy homeostasis and appetite, is a class A G-protein-coupled receptor and a prime target for the pharmacological treatment of obesity. Here, we present cryo-electron microscopy structures of MC4R–Gs-protein complexes with two drugs recently approved by the FDA, the peptide agonists NDP-α-MSH and setmelanotide, with 2.9 Å and 2.6 Å resolution. Together with signaling data from structure-derived MC4R mutants, the complex structures reveal the agonist-induced origin of transmembrane helix (TM) 6-regulated receptor activation. The ligand-binding modes of NDP-α-MSH, a high-affinity linear variant of the endogenous agonist α-MSH, and setmelanotide, a cyclic anti-obesity drug with biased signaling toward Gq/11, underline the key role of TM3 in ligand-specific interactions and of calcium ion as a ligand-adaptable cofactor. The agonist-specific TM3 interplay subsequently impacts receptor–Gs-protein interfaces at intracellular loop 2, which also regulates the G-protein coupling profile of this promiscuous receptor. Finally, our structures reveal mechanistic details of MC4R activation/inhibition, and provide important insights into the regulation of the receptor signaling profile which will facilitate the development of tailored anti-obesity drugs.

Structures of the glucocorticoid-bound adhesion receptor GPR97-Go complex

Adhesion G-protein-coupled receptors (GPCRs) are a major family of GPCRs, but limited knowledge of their ligand regulation or structure is available^1-3. Here we report that glucocorticoid stress hormones activate adhesion G-protein-coupled receptor G3 (ADGRG3; also known as GPR97)^4-6, a prototypical adhesion GPCR. The cryo-electron microscopy structures of GPR97-G_o complexes bound to the anti-inflammatory drug beclomethasone or the steroid hormone cortisol revealed that glucocorticoids bind to a pocket within the transmembrane domain. The steroidal core of glucocorticoids is packed against the 'toggle switch' residue W^6.53, which senses the binding of a ligand and induces activation of the receptor. Active GPR97 uses a quaternary core and HLY motif to fasten the seven-transmembrane bundle and to mediate G protein coupling. The cytoplasmic side of GPR97 has an open cavity, where all three intracellular loops interact with the G_o protein, contributing to the high basal activity of GRP97. Palmitoylation at the cytosolic tail of the G_o protein was found to be essential for efficient engagement with GPR97 but is not observed in other solved GPCR complex structures. Our work provides a structural basis for ligand binding to the seven-transmembrane domain of an adhesion GPCR and subsequent G protein coupling.

Structures of the glucocorticoid-bound adhesion receptor GPR97-Go complex

Adhesion G-protein-coupled receptors (GPCRs) are a major family of GPCRs, but limited knowledge of their ligand regulation or structure is available^1-3. Here we report that glucocorticoid stress hormones activate adhesion G-protein-coupled receptor G3 (ADGRG3; also known as GPR97)^4-6, a prototypical adhesion GPCR. The cryo-electron microscopy structures of GPR97-G_o complexes bound to the anti-inflammatory drug beclomethasone or the steroid hormone cortisol revealed that glucocorticoids bind to a pocket within the transmembrane domain. The steroidal core of glucocorticoids is packed against the 'toggle switch' residue W^6.53, which senses the binding of a ligand and induces activation of the receptor. Active GPR97 uses a quaternary core and HLY motif to fasten the seven-transmembrane bundle and to mediate G protein coupling. The cytoplasmic side of GPR97 has an open cavity, where all three intracellular loops interact with the G_o protein, contributing to the high basal activity of GRP97. Palmitoylation at the cytosolic tail of the G_o protein was found to be essential for efficient engagement with GPR97 but is not observed in other solved GPCR complex structures. Our work provides a structural basis for ligand binding to the seven-transmembrane domain of an adhesion GPCR and subsequent G protein coupling.

Adaptive evolution of GPR39 in diverse directions in vertebrates

G protein-coupled receptors (GPCRs) play an important role in physiology and disease and represent productive drug targets. Orphan GPCRs, which have unknown endogenous ligands, are considered drug targets and consequently have attracted great interest in identifying their endogenous cognate ligands for deorphanization. However, additional studies have shown that GPCRs, including many orphan GPCRs, can constitutively activate G protein signaling in a ligand-independent manner. GPR39 is such an orphan GPCR with constitutive activity. Here, we performed a phylogenetic and selection analysis of GPR39 in vertebrates, and we found that GPR39 underwent positive selection in different branches of vertebrates. Using luciferase reporter assays, we demonstrated that human, frog and chicken GPR39 can constitutively activate Gq and G12 signaling pathways in a ligand-independent manner. Zebrafish GPR39 can constitutively activate Gs, Gq and G12 signaling pathways in a ligand-independent manner. We further found that the zebrafish-H296^7.35 site is crucial for the activity of the Gs signaling pathway. In addition, our mutagenesis studies indicated that the positive selection sites of GPR39 from different species had important effects on the constitutive activity of the receptor. Our results revealed the adaptive evolution of GPR39 in diverse directions, which led to differences in constitutive activity.

Structure and selectivity engineering of the M1 muscarinic receptor toxin complex

Muscarinic toxins (MTs) are natural toxins produced by mamba snakes that primarily bind to muscarinic acetylcholine receptors (MAChRs) and modulate their function. Despite their similar primary and tertiary structures, MTs show distinct binding selectivity toward different MAChRs. The molecular details of how MTs distinguish MAChRs are not well understood. Here, we present the crystal structure of M₁AChR in complex with MT7, a subtype-selective anti-M₁AChR snake venom toxin. The structure reveals the molecular basis of the extreme subtype specificity of MT7 for M₁AChR and the mechanism by which it regulates receptor function. Through in vitro engineering of MT7 finger regions that was guided by the structure, we have converted the selectivity from M₁AChR toward M₂AChR, suggesting that the three-finger fold is a promising scaffold for developing G protein-coupled receptor modulators.

Pseudo-Symmetric Assembly of Protodomains as a Common Denominator in the Evolution of Polytopic Helical Membrane Proteins

The polytopic helical membrane proteome is dominated by proteins containing seven transmembrane helices (7TMHs). They cannot be grouped under a monolithic fold or superfold. However, a parallel structural analysis of folds around that magic number of seven in distinct protein superfamilies (SWEET, PnuC, TRIC, FocA, Aquaporin, GPCRs) reveals a common homology, not in their structural fold, but in their systematic pseudo-symmetric construction during their evolution. Our analysis leads to guiding principles of intragenic duplication and pseudo-symmetric assembly of ancestral transmembrane helical protodomains, consisting of 3 (or 4) helices. A parallel deconstruction and reconstruction of these domains provides a structural and mechanistic framework for their evolutionary paths. It highlights the conformational plasticity inherent to fold formation itself, the role of structural as well as functional constraints in shaping that fold, and the usefulness of protodomains as a tool to probe convergent vs divergent evolution. In the case of FocA vs. Aquaporin, this protodomain analysis sheds new light on their potential divergent evolution at the protodomain level followed by duplication and parallel evolution of the two folds. GPCR domains, whose function does not seem to require symmetry, nevertheless exhibit structural pseudo-symmetry. Their construction follows the same protodomain assembly as any other pseudo-symmetric protein suggesting their potential evolutionary origins. Interestingly, all the 6/7/8TMH pseudo-symmetric folds in this study also assemble as oligomeric forms in the membrane, emphasizing the role of symmetry in evolution, revealing self-assembly and co-evolution not only at the protodomain level but also at the domain level.

Human red and green cone opsins are O-glycosylated at an N-terminal Ser/Thr-rich domain conserved in vertebrates

There are fundamental differences in the structures of outer segments between rod and cone photoreceptor cells in the vertebrate retina. Visual pigments are the only essential membrane proteins that differ between rod and cone outer segments, making it likely that they contribute to these structural differences. Human rhodopsin is N-glycosylated on Asn² and Asn¹⁵, whereas human (h) red and green cone opsins (hOPSR and hOPSG, respectively) are N-glycosylated at Asn³⁴ Here, utilizing a monoclonal antibody (7G8 mAB), we demonstrate that hOPSR and hOPSG from human retina also are O-glycosylated with full occupancy. We determined that 7G8 mAB recognizes the N-terminal sequence ²¹DSTQSSIF²⁸ of hOPSR and hOPSG from extracts of human retina, but only after their O-glycans have been removed with O-glycosidase treatment, thus revealing this post-translational modification of red and green cone opsins. In addition, we show that hOPSR and hOPSG from human retina are recognized by jacalin, a lectin that binds to O-glycans, preferentially to Gal-GalNAc. Next, we confirmed the presence of O-glycans on OPSR and OPSG from several vertebrate species, including mammals, birds, and amphibians. Finally, the analysis of bovine OPSR by MS identified an O-glycan on Ser²², a residue that is semi-conserved (Ser or Thr) among vertebrate OPSR and OPSG. These results suggest that O-glycosylation is a fundamental feature of red and green cone opsins, which may be relevant to their function or to cone cell development, and that differences in this post-translational modification also could contribute to the different morphologies of rod and cone photoreceptors.

JC Polyomavirus Entry by Clathrin-Mediated Endocytosis Is Driven by β-Arrestin

JC polyomavirus (JCPyV) establishes a persistent, lifelong, asymptomatic infection within the kidney of the majority of the human population. Under conditions of severe immunosuppression or immune modulation, JCPyV can reactivate in the central nervous system (CNS) and cause progressive multifocal leukoencephalopathy (PML), a fatal demyelinating disease. Initiation of infection is mediated through viral attachment to α2,6-sialic acid-containing lactoseries tetrasaccharide c (LSTc) on the surface of host cells. JCPyV internalization is dependent on serotonin 5-hydroxytryptamine subfamily 2 receptors (5-HT₂Rs), and entry is thought to occur by clathrin-mediated endocytosis (CME). However, the JCPyV entry process and the cellular factors involved in viral internalization remain poorly understood. Treatment of cells with small-molecule chemical inhibitors and RNA interference of 5-HT₂R endocytic machinery, including β-arrestin, clathrin, AP2, and dynamin, significantly reduced JCPyV infection. However, infectivity of the polyomavirus simian virus 40 (SV40) was not affected by CME-specific treatments. Inhibition of clathrin or β-arrestin specifically reduced JCPyV internalization but did not affect viral attachment. Furthermore, mutagenesis of a β-arrestin binding domain (Ala-Ser-Lys) within the intracellular C terminus of 5-HT_2AR severely diminished internalization and infection, suggesting that β-arrestin interactions with 5-HT_2AR are critical for JCPyV infection and entry. These conclusions illuminate key host factors that regulate clathrin-mediated endocytosis of JCPyV, which is necessary for viral internalization and productive infection.IMPORTANCE Viruses usurp cellular factors to invade host cells. Activation and utilization of these proteins upon initiation of viral infection are therefore required for productive infection and resultant viral disease. The majority of healthy individuals are asymptomatically infected by JC polyomavirus (JCPyV), but if the host immune system is compromised, JCPyV can cause progressive multifocal leukoencephalopathy (PML), a rare, fatal, demyelinating disease. Individuals infected with HIV or taking prolonged immunomodulatory therapies have a heightened risk for developing PML. The cellular proteins and pathways utilized by JCPyV to mediate viral entry are poorly understood. Our findings further characterize how JCPyV utilizes the clathrin-mediated endocytosis pathway to invade host cells. We have identified specific components of this pathway that are necessary for the viral entry process and infection. Collectively, the conclusions increase our understanding of JCPyV infection and pathogenesis and may contribute to the future development of novel therapeutic strategies for PML.

Classics in Chemical Neuroscience: Chlorpromazine

The discovery of chlorpromazine in the early 1950s revolutionized the clinical treatment of schizophrenia, galvanized the development of psychopharmacology, and standardized protocols used for testing the clinical efficacy of antipsychotics. Furthermore, chlorpromazine expanded our understanding of the role of chemical messaging in neurotransmission and reduced the stigma associated with mental illness, facilitating deinstitutionalization in the 1960s and 1970s. In this review, we will discuss the synthesis, manufacturing, metabolism and pharmacokinetics, pharmacology, structure-activity relationship, and adverse effects of chlorpromazine. In conclusion, we summarize the history and significant contributions of chlorpromazine that have resulted in this potent first-generation antipsychotic maintaining its clinical relevance for nearly 70 years.

The Model Structures of the Complement Component 5a Receptor (C5aR) Bound to the Native and Engineered hC5a

The interaction of ^hC5a with C5aR, previously hypothesized to involve a “two-site” binding, (i) recognition of the bulk of ^hC5a by the N-terminus (NT) of C5aR (“site1”), and (ii) recognition of C-terminus (CT) of ^hC5a by the extra cellular surface (ECS) of the C5aR (“site2”). However, the pharmacological landscapes of such recognition sites are yet to be illuminated at atomistic resolution. In the context, unique model complexes of C5aR, harboring pharmacophores of diverse functionality at the “site2” has recently been described. The current study provides a rational illustration of the “two-site” binding paradigm in C5aR, by recruiting the native agonist ^hC5a and engineered antagonist ^hC5a(A8). The ^hC5a-C5aR and ^hC5a(A8)-C5aR complexes studied over 250 ns of molecular dynamics (MD) each in POPC bilayer illuminate the hallmark of activation mechanism in C5aR. The intermolecular interactions in the model complexes are well supported by the molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) based binding free energy calculation, strongly correlating with the reported mutational studies. Exemplified in two unique and contrasting molecular complexes, the study provides an exceptional understanding of the pharmacological divergence observed in C5aR, which will certainly be useful for search and optimization of new generation “neutraligands” targeting the ^hC5a-C5aR interaction.

The G Protein-Coupled Receptor UT of the Neuropeptide Urotensin II Displays Structural and Functional Chemokine Features

The urotensinergic system was previously considered as being linked to numerous physiopathological states, including atherosclerosis, heart failure, hypertension, pre-eclampsia, diabetes, renal disease, as well as brain vascular lesions. Thus, it turns out that the actions of the urotensin II (UII)/G protein-coupled receptor UT system in animal models are currently not predictive enough in regard to their effects in human clinical trials and that UII analogs, established to target UT, were not as beneficial as expected in pathological situations. Thus, many questions remain regarding the overall signaling profiles of UT leading to complex involvement in cardiovascular and inflammatory responses as well as cancer. We address the potential UT chemotactic structural and functional definition under an evolutionary angle, by the existence of a common conserved structural feature among chemokine receptorsopioïdergic receptors and UT, i.e., a specific proline position in the transmembrane domain-2 TM2 (P2.58) likely responsible for a kink helical structure that would play a key role in chemokine functions. Even if the last decade was devoted to the elucidation of the cardiovascular control by the urotensinergic system, we also attempt here to discuss the role of UII on inflammation and migration, likely providing a peptide chemokine status for UII. Indeed, our recent work established that activation of UT by a gradient concentration of UII recruits Gαi/o and Gα13 couplings in a spatiotemporal way, controlling key signaling events leading to chemotaxis. We think that this new vision of the urotensinergic system should help considering UT as a chemotactic therapeutic target in pathological situations involving cell chemoattraction.

Structure-Based Sequence Alignment of the Transmembrane Domains of All Human GPCRs: Phylogenetic, Structural and Functional Implications

The understanding of G-protein coupled receptors (GPCRs) is undergoing a revolution due to increased information about their signaling and the experimental determination of structures for more than 25 receptors. The availability of at least one receptor structure for each of the GPCR classes, well separated in sequence space, enables an integrated superfamily-wide analysis to identify signatures involving the role of conserved residues, conserved contacts, and downstream signaling in the context of receptor structures. In this study, we align the transmembrane (TM) domains of all experimental GPCR structures to maximize the conserved inter-helical contacts. The resulting superfamily-wide GpcR Sequence-Structure (GRoSS) alignment of the TM domains for all human GPCR sequences is sufficient to generate a phylogenetic tree that correctly distinguishes all different GPCR classes, suggesting that the class-level differences in the GPCR superfamily are encoded at least partly in the TM domains. The inter-helical contacts conserved across all GPCR classes describe the evolutionarily conserved GPCR structural fold. The corresponding structural alignment of the inactive and active conformations, available for a few GPCRs, identifies activation hot-spot residues in the TM domains that get rewired upon activation. Many GPCR mutations, known to alter receptor signaling and cause disease, are located at these conserved contact and activation hot-spot residue positions. The GRoSS alignment places the chemosensory receptor subfamilies for bitter taste (TAS2R) and pheromones (Vomeronasal, VN1R) in the rhodopsin family, known to contain the chemosensory olfactory receptor subfamily. The GRoSS alignment also enables the quantification of the structural variability in the TM regions of experimental structures, useful for homology modeling and structure prediction of receptors. Furthermore, this alignment identifies structurally and functionally important residues in all human GPCRs. These residues can be used to make testable hypotheses about the structural basis of receptor function and about the molecular basis of disease-associated single nucleotide polymorphisms.

G-protein coupled receptors (GPCRs) are a large superfamily of integral membrane proteins that share a characteristic 7 transmembrane helix fold. They detect various molecules outside of the cell and signal their presence to the inside of the cell. At least half of the 800 human GPCRs are potential drug targets, so understanding their structure and function is critical. Experimental structures are now available for at least one receptor from each GPCR class. The structure of the 7 helix fold is highly conserved even for receptors with very low sequence similarity. We analyze the available experimental structures and compare the common inter-helical contacts. Our analysis leads to a unified sequence-structure alignment of the GPCR superfamily that can then be used as the starting point for structure prediction of all other GPCRs. A key result of our analysis is a list of conserved contact residues and activation “hot-spots” residues that are critical for GPCR folding and function. We propose that mutations and natural variants of amino acids at these locations in the GPCRs can dramatically influence their activation state and alter intracellular signaling. This provides hypotheses for the molecular mechanisms underlying disease causing mutants for any GPCR.

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

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

Functional elements of the gastric inhibitory polypeptide receptor: Comparison between secretin- and rhodopsin-like G protein-coupled receptors

Innovative crystallographic techniques have resulted in an exponential growth in the number of solved G-protein coupled receptor (GPCR) structures and a better understanding of the mechanisms of class A receptor activation and G protein binding. The recent release of the type 1 receptor for the corticotropin-releasing factor and the glucagon receptor structures, two members of the secretin-like family, gives the opportunity to understand these mechanisms of activation in this family of GPCRs. Here, we addressed the comparison of the functional elements of class A and secretin-like GPCRs, using the glucose-dependent insulinotropic polypeptide receptor (GIPR) as a model receptor. Inactive and active models of GIPR permitted to select, by structural homology with class A GPCRs, several residues that may form key interactions presumably involved in receptor activation and Gs coupling, for pharmacological evaluation. Mutants on these amino acids were expressed in HEKT 293 cells and characterized in terms of GIP-induced cAMP production. We identified various functional domains spanning from the peptide-binding to the G protein pockets: including: a network linking the extracellular part of transmembrane (TM) 6 with TMs 2 and 7; a polar lock that resembles the ionic-lock in class A GPCRs; an interaction between TMs 3 and 7 that favors activation; and two clusters of polar/charged and of hydrophobic residues that interact with the C-terminus of the Gα. The results show that despite the low degree of sequence similarity between rhodopsin- and secretin-like GPCRs, the two families share conserved elements in their mechanisms of activation and G protein binding.

Model structures of inactive and peptide agonist bound C5aR: Insights into agonist binding, selectivity and activation

C5a receptor (C5aR) is one of the major chemoattractant receptors of the druggable proteome that binds C5a, the proinflammatory polypeptide of complement cascade, triggering inflammation and SEPSIS. Here, we report the model structures of C5aR in both inactive and peptide agonist (YSFKPMPLaR; a=D-Ala) bound meta-active state. Assembled in CYANA and evolved over molecular dynamics (MD) in POPC bilayer, the inactive C5aR demonstrates a topologically unique compact heptahelical bundle topology harboring a β-hairpin in extracellular loop 2 (ECL2), derived from the atomistic folding simulations. The peptide agonist bound meta-active C5aR deciphers the “site2” at an atomistic resolution in the extracellular surface (ECS), in contrast to the previously hypothesized inter-helical crevice. With estimated Ki≈2.75 μM, the meta-active C5aR excellently rationalizes the IC₅₀ (0.1–13 μM) and EC₅₀ (0.01–6 μM) values, displayed by the peptide agonist in several signaling studies. Moreover, with Ki≈5.3×10⁵ μM, the “site2” also illustrates selectivity, by discriminating the stereochemical mutant peptide (YSFkPMPLaR; k=D-Lys), known to be inert toward C5aR, up to 1 mM concentration. Topologically juxtaposed between the structures of rhodopsin and CXCR1, the C5aR models also display excellent structural correlations with the other G-protein coupled receptors (GPCRs). The models elaborated in the current study unravel many important structural insights previously not known for regulating the agonist binding and activation mechanism of C5aR.

•

Topologically unique inactive and meta-active atomistic models of C5aR.

•

Model demonstrates excellent structural correlation with the other reported GPCRs.

•

Model deciphers the “site2” in the ECS and also demonstrates agonist selectivity.

•

Agonist binding and activation requires “cation–π” interaction with F275 of C5aR.

•

Inactive to meta-active transition involves TM3–TM6 movements (ΔΘ≈+11.1°) in C5aR.

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

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

Cholesterol as a co-solvent and a ligand for membrane proteins

As of mid 2013 a Medline search on "cholesterol" yielded over 200,000 hits, reflecting the prominence of this lipid in numerous aspects of animal cell biology and physiology under conditions of health and disease. Aberrations in cholesterol homeostasis underlie both a number of rare genetic disorders and contribute to common sporadic and complex disorders including heart disease, stroke, type II diabetes, and Alzheimer's disease. The corresponding author of this review and his lab stumbled only recently into the sprawling area of cholesterol research when they discovered that the amyloid precursor protein (APP) binds cholesterol, a topic covered by the Hans Neurath Award lecture at the 2013 Protein Society Meeting. Here, we first provide a brief overview of cholesterol-protein interactions and then offer our perspective on how and why binding of cholesterol to APP and its C99 domain (β-CTF) promotes the amyloidogenic pathway, which is closely related to the etiology of Alzheimer's disease.

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.

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.

Interactions of vasopressin and oxytocin receptors with vasopressin analogues substituted in position 2 with 3,3'-diphenylalanine--a molecular docking study

Vasopressin and oxytocin receptors belong to the superfamily of G protein-coupled receptors and play an important role in many physiological functions. They are also involved in a number of pathological conditions being important drug targets. In this work, four vasopressin analogues substituted at position 2 with 3,3'-diphenylalanine have been docked into partially flexible vasopressin and oxytocin receptors. The bulky residue at position 2 acts as a structural restraint much stronger in the oxytocin receptor (OTR) than in the vasopressin V2 receptor (V2R), resulting in a different location of the analogues in these receptors. This explains the different, either agonistic or antagonistic, activities of the analogues in V2R and OTR, respectively. In all complexes, the conserved polar residues serve as anchor points for the ligand both in OTR and V2R. Strong interactions of the C-terminus of analogue II ([Mpa(1) ,d-Dpa(2) ,Val(4) ,d-Arg(8) ]VP) with extracellular loop 3 may be responsible for its highest activity at V2R. It also appears that V2R adapts more readily to the docking analogues by conformational changes in the aromatic side chains triggering receptor activation. A weak activity at V1a vasopressin receptor appears to be caused by weak receptor-ligand interactions. Results of this study may facilitate a rational design of new analogues with the highest activity/selectivity at vasopressin and OTRs.

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.

The essential role for aromatic cluster in the β3 adrenergic receptor

AIM

To explore the function of the conserved aromatic cluster F213(5.47), F308(6.51), and F309(6.52) in human β3 adrenergic receptor (hβ3AR).

METHODS

Point mutation technology was used to produce plasmid mutations of hβ3AR. HEK-293 cells were transiently co-transfected with the hβ3AR (wild-type or mutant) plasmids and luciferase reporter vector pCRE-luc. The expression levels of hβ3AR in the cells were determined by Western blot analysis. The constitutive signalling and the signalling induced by the β3AR selective agonist, BRL (BRL37344), were then evaluated. To further explore the interaction mechanism between BRL and β3AR, a three-dimensional complex model of β3AR and BRL was constructed by homology modelling and molecular docking.

RESULTS

For F308(6.51), Ala and Leu substitution significantly decreased the constitutive activities of β3AR to approximately 10% of that for the wild-type receptor. However, both the potency and maximal efficacy were unchanged by Ala substitution. In the F308(6.51)L construct, the EC(50) value manifested as a "right shift" of approximately two orders of magnitude with an increased E(max). Impressively, the molecular pharmacological phenotype was similar to the wild-type receptor for the introduction of Tyr at position 308(6.51), though the EC(50) value increased by approximately five-fold for the mutant. For F309(6.52), the constitutive signalling for both F309(6.52)A and F309(6.52)L constructs were strongly impaired. In the F309(6.52)A construct, BRL-stimulated signalling showed a normal E(max) but reduced potency. Leu substitution of F309(6.52) reduced both the E(max) and potency. When F309(6.52) was mutated to Tyr, the constitutive activity was decreased approximately three-fold, and BRL-stimulated signalling was significantly impaired. Furthermore, the double mutant (F308(6.51)A_F309(6.52)A) caused the total loss of β3AR function. The predicted binding mode between β3AR and BRL revealed that both F308(6.51) and F309(6.52) were in the BRL binding pocket of β3AR, while F213(5.47) and W305(6.48) were distant from the binding site.

CONCLUSION

These results revealed that aromatic residues, especially F308(6.51) and F309(6.52), play essential roles in the function of β3AR. Aromatic residues maintained the receptor in a partially activated state and significantly contributed to ligand binding. The results supported the common hypothesis that the aromatic cluster F[Y]5.47/F[Y]6.52/F[Y]6.51 conserved in class A G protein-coupled receptor (GPCR) plays an important role in the structural stability and activation of GPCRs.

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.

Do crystal structures obviate the need for theoretical models of GPCRs for structure-based virtual screening?

Recent highly expected structural characterizations of agonist-bound and antagonist-bound beta-2 adrenoreceptor (β2AR) by X-ray crystallography have been widely regarded as critical advances to enable more effective structure-based discovery of GPCRs ligands. It appears that this very important development may have undermined many previous efforts to develop 3D theoretical models of GPCRs. To address this question directly, we have compared several historical β2AR models versus the inactive state and nanobody-stabilized active state of β2AR crystal structures in terms of their structural similarity and effectiveness of use in virtual screening for β2AR specific agonists and antagonists. Theoretical models, incluing both homology and de novo types, were collected from five different groups who have published extensively in the field of GPCRs modeling. All models were built before X-ray structures became available. In general, β2AR theoretical models differ significantly from the crystal structure in terms of TMH definition and the global packing. Nevertheless, surprisingly, several models afforded hit rates resulting from virtual screening of large chemical library enriched by known β2AR ligands that exceeded those using X-ray structures. The hit rates were particularly higher for agonists. Furthemore, the screening performance of models is associated with local structural quality, such as the RMSDs for binding pocket residues and the ability to capture accurately, most if not all critical protein/ligand interactions. These results suggest that carefully built models of GPCRs could capture critical chemical and structural features of the binding pocket, and thus may be even more useful for practical structure-based drug discovery than X-ray structures.

Developing a high-quality scoring function for membrane protein structures based on specific inter-residue interactions

Membrane proteins are of particular biological and pharmaceutical importance, and computational modeling and structure prediction approaches play an important role in studies of membrane proteins. Developing an accurate model quality assessment program is of significance to the structure prediction of membrane proteins. Few such programs are proposed that can be applied to a broad range of membrane protein classes and perform with high accuracy. We developed a new model scoring function Interaction-based Quality assessment (IQ), based on the analysis of four types of inter-residue interactions within the transmembrane domains of helical membrane proteins. This function was tested using three high-quality model sets: all 206 models of GPCR Dock 2008, all 284 models of GPCR Dock 2010, and all 92 helical membrane protein models of the HOMEP set. For all three sets, the scoring function can select the native structures among all of the models with the success rates of 93, 85, and 100% respectively. For comparison, these three model sets were also adopted for a recently published model assessment program for membrane protein structures, ProQM, which gave the success rates of 85, 79, and 92% separately. These results suggested that IQ outperforms ProQM when only the transmembrane regions of the models are considered. This scoring function should be useful for the computational modeling of membrane proteins.

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 substrate-driven transition to an inward-facing conformation in the functional mechanism of the dopamine transporter

Background

The dopamine transporter (DAT), a member of the neurotransmitter:Na⁺ symporter (NSS) family, terminates dopaminergic neurotransmission and is a major molecular target for psychostimulants such as cocaine and amphetamine, and for the treatment of attention deficit disorder and depression. The crystal structures of the prokaryotic NSS homolog of DAT, the leucine transporter LeuT, have provided critical structural insights about the occluded and outward-facing conformations visited during the substrate transport, but only limited clues regarding mechanism. To understand the transport mechanism in DAT we have used a homology model based on the LeuT structure in a computational protocol validated previously for LeuT, in which steered molecular dynamics (SMD) simulations guide the substrate along a pathway leading from the extracellular end to the intracellular (cytoplasmic) end.

Methodology/Principal Findings

Key findings are (1) a second substrate binding site in the extracellular vestibule, and (2) models of the conformational states identified as occluded, doubly occupied, and inward-facing. The transition between these states involve a spatially ordered sequence of interactions between the two substrate-binding sites, followed by rearrangements in structural elements located between the primary binding site and the cytoplasmic end. These rearrangements are facilitated by identified conserved hinge regions and a reorganization of interaction networks that had been identified as gates.

Conclusions/Significance

Computational simulations supported by information available from experiments in DAT and other NSS transporters have produced a detailed mechanistic proposal for the dynamic changes associated with substrate transport in DAT. This allosteric mechanism is triggered by the binding of substrate in the S2 site in the presence of the substrate in the S1 site. Specific structural elements involved in this mechanism, and their roles in the conformational transitions illuminated here describe, a specific substrate-driven allosteric mechanism that is directly amenable to experiment as shown previously for LeuT.

Probing the structural determinants for the function of intracellular loop 2 in structurally cognate G-protein-coupled receptors

Intracellular loop 2 (IL2) in G-protein-coupled receptors (GPCRs) is functionally important, e.g., in binding to G-protein and β-arrestin. Differences in secondary structure of IL2 in the crystal structures of the very similar β(1)- and β(2)-adrenergic receptors (β(1)AR and β(2)AR, respectively), i.e., an α-helix and an L-shaped strand, respectively, emphasize the need to understand the structural basis for IL2 functionality. We studied the properties of IL2 in the context of experimental data using a Monte Carlo-based ab initio method. The procedure was validated first by verifying that the IL2 structures in β(1)AR and β(2)AR crystals were correctly reproduced, even after conformational ensemble searches at >1200 K where most secondary structure had been lost. We found that IL2 in β(1)AR and β(2)AR sampled each other's conformation but adopted different energetically preferred conformations, consistent with the crystal structures. The results indicate a persistent contextual preference for the structure of IL2, which was conserved when the IL2 sequences were interchanged between the receptors. We conclude that the protein environment, more than the IL2 sequence, regulates the IL2 structures. We extended the approach to the molecular model of 5-HT(2A)R for which no crystal structure is available and found that IL2 is predominantly helical, similar to IL2 in β(1)AR. Because the P3.57A mutation in IL2 had been shown to decrease β-arrestin binding and internalization, we predicted the effects of the mutation and found that it decreased the propensity of IL2 to form helix, identifying the helical IL2 as a component of the GPCR active form.

Cholesterol modulates the membrane effects and spatial organization of membrane-penetrating ligands for G-protein coupled receptors

The ligands of certain G-protein coupled receptors (GPCRs) are membrane soluble and reach their target from the lipid bilayer. Lipid composition and dynamics will therefore modulate the activity of these receptors, but specific roles of lipid components, including the ubiquitous cholesterol (Chol), are not clear. We have probed the organization and dynamics of such a lipid-bilayer-penetrating ligand, the endogenous ligand for the κ-opioid receptor (KOR) dynorphin A (1-17) (DynA), using molecular dynamics (MD) simulations of DynA in cholesterol-depleted and cholesterol-enriched model membranes. DynA is found to penetrate deep inside fluid dimyristoylphosphatidylcholine (DMPC) bilayers, and resides with its N-terminal helix at ∼6 Å away from the bilayer midplane, in a tilted orientation, at an ∼50° angle with respect to the membrane normal. In contrast, DynA inside DMPC/Chol membranes with 20% cholesterol (DMPC/Chol) is situated with its helical segment ∼5 A higher, i.e., closer to the lipid/water interface and in a relatively vertical orientation. The DMPC membrane shows greater thinning around the insertion and permits a stronger influx of water inside the hydrocarbon core than the DMPC/Chol membranes. Relating these results to data about key GPCR residues that have been implicated in interactions with membrane-inserting GPCR ligands, we conclude that the position of DynA in DMPC/Chol, but not in pure DMPC, correlates with generally proposed GPCR ligand entry pathways. Our predictions provide a possible mechanistic explanation as to why DynA binding to KOR, and the subsequent activation of the receptor, is facilitated in cholesterol-enriched environments. A quantitative description of DynA-induced membrane deformations is obtained with a continuum theory of membrane deformations (CTMD) that is based on hydrophobic matching. Comparison with the MD data reveals the significance of the lipid tail packing energy contribution in the DMPC/Chol mixtures in predicting equilibrium membrane shape around DynA. On this basis, specific corrections are introduced to this energy term within the CTMD framework, thereby extending the applicability of the CTMD framework to lipid raft mixtures and their interactions with GPCR proteins and their ligands.

Accurate prediction of the burial status of transmembrane residues of α-helix membrane protein by incorporating the structural and physicochemical features

Predicting the burial status (the residue exposure to the lipid bilayer or buried within the protein core) of transmembrane (TM) residues of α-helix membrane protein (αHMP) is of great importance for genome-wide annotation and for experimental researchers to elucidate diverse physiological processes. In this work, we developed a new computational model that can be used for predicting the burial status of TM residues of αHMP. By incorporating physicochemical scales and conservation index, an efficient prediction model using least squares support vector machine (LS-SVM) was developed. The model was developed from 43 protein chains and its prediction ability was evaluated by an independent test set of other non-redundant ten protein chains. The prediction accuracy of our method was much better than the results of the reported works. Our results demonstrate that the LS-SVM prediction model incorporating structural and physicochemical features derived from sequence information could greatly improve the prediction accuracy.

Third extracellular loop (EC3)-N terminus interaction is important for seven-transmembrane domain receptor function: implications for an activation microswitch region

The canonical heptahelical bundle architecture of seven-transmembrane domain (7TM) receptors is intertwined by three intra- and three extracellular loops, whose local conformations are important in receptor signaling. Many 7TM receptors contain a cysteine residue in the third extracellular loop (EC3) and a complementary cysteine residue on the N terminus. The functional role of such EC3-N terminus conserved cysteine pairs remains unclear. This study explores the role of the EC3-N terminus cysteine pairs on receptor conformation and G protein activation by disrupting them in the chemokine receptor CXCR4, while engineering a novel EC3-N terminus cysteine pair into the complement factor 5a receptor (C5aR), a chemo attractant receptor that lacks it. Mutated CXCR4 and C5aRs were expressed in engineered yeast. Mutation of the cysteine pair with the serine pair (C28S/C274S) in constitutively active mutant CXCR4 abrogated the receptor activation, whereas mutation with the aromatic pair (C28F-C274F) or the salt bridge pair (C28R/C274E), respectively, rescued or retained the receptor activation in response to CXCL12. In this context, the cysteine pair (Cys(30) and Cys(272)) engineered into the EC3-N terminus (Ser(30) and Ser(272)) of a novel constitutively active mutant of C5aR restrained the constitutive signaling without affecting the C5a-induced activation. Further mutational studies demonstrated a previously unappreciated role for Ser(272) on EC3 of C5aR and its interaction with the N terminus, thus defining a new microswitch region within the C5aR. Similar results were obtained with mutated CXCR4 and C5aRs expressed in COS-7 cells. These studies demonstrate a novel role of the EC3-N terminus cysteine pairs in G protein-coupled receptor activation and signaling.

Structure and activation of the visual pigment rhodopsin

Rhodopsin is a specialized G protein-coupled receptor (GPCR) found in vertebrate rod cells. Absorption of light by its 11-cis retinal chromophore leads to rapid photochemical isomerization and receptor activation. Recent results from protein crystallography and NMR spectroscopy show how structural changes on the extracellular side of rhodopsin induced by retinal isomerization are coupled to the motion of membrane-spanning helices to create a G protein binding pocket on the intracellular side of the receptor. The signaling pathway provides a comprehensive explanation for the conservation of specific amino acids and structural motifs across the class A family of GPCRs, as well as for the conservation of selected residues within the visual receptor subfamily. The emerging model of activation indicates that, rather than being unique, the visual receptors provide a basis for understanding the common structural and dynamic elements in the class A GPCRs.

Engineering a GPCR-ligand pair that simulates the activation of D(2L) by Dopamine

In the past decade, engineered G-protein-coupled receptors activated solely by synthetic ligands (RASSLs) have been implemented as a new means to study neurotransmission, which is controlled by G-protein-coupled receptors in vitro and in vivo. In this study, we report an engineered dopamine receptor D(2L) F390(6.52)W, which is the first identified RASSL for the dopamine receptor family. The mutant receptor is characterized by a disrupted ligand binding and complete loss of efficacy for the endogenous ligand, dopamine, which is putatively due to a sterically induced perturbation of H-bonding with conserved serine residues in TM5. Based on this model, we rationally developed an aminoindane-derived set of agonists. Because these agonists forgo analogous H-bonding functionalities, their binding energy does not depend on the respective interactions. Binding affinity and potency were optimized by ligand modifications bearing molecular appendages that obviously interact with a secondary recognition site provided by four hydrophobic residues in TM2 and TM3. Thus, the ferrocenyl carboxamide 5b (FAUC 185) was identified as a synthetic agonist that is able to stimulate the mutant receptor in a manner similar to that by which endogenous dopamine activates the D(2L) wild-type receptor. The engineered dopamine receptor D(2L) F390(6.52)W in combination with FAUC 185 (5b) provides a new tool to probe GPCR functions selectively in specific cell populations in vitro and in vivo.

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.

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.

Comparing four different approaches for the determination of inter-residue interactions provides insight for the structure prediction of helical membrane proteins

Studying inter-residue interactions provides insight into the folding and stability of both soluble and membrane proteins and is essential for developing computational tools for protein structure prediction. As the first step, various approaches for elucidating such interactions within protein structures have been proposed and proven useful. Since different approaches may grasp different aspects of protein structural folds, it is of interest to systematically compare them. In this work, we applied four approaches for determining inter-residue interactions to the analysis of three distinct structure datasets of helical membrane proteins and compared their correlation to the three individual quality measures of structures in these datasets. These datasets included one of 35 structures of rhodopsin receptors and bacterial rhodopsins determined at various resolutions, one derived from the HOMEP benchmark dataset previously reported, and one comprising of 139 homology models. It was found that the correlation between the average number of inter-residue interactions obtained by applying the four approaches and the available structure quality measures varied quite significantly among them. The best correlation was achieved by the approach focusing exclusively on favorable inter-residue interactions. These results provide interesting insight for the development of objective quality measure for the structure prediction of helical membrane proteins.