Lifting the lid on GPCRs: the role of extracellular loops

Multi-target dopamine D3 receptor modulators: Actionable knowledge for drug design from molecular dynamics and machine learning

Local changes in the structure of G-protein coupled receptors (GPCR) binders largely affect their pharmacological profile. While the sought efficacy can be empirically obtained by introducing local modifications, the underlining structural explanation can remain elusive. Here, molecular dynamics (MD) simulations of the eticlopride-bound inactive state of the Dopamine D3 Receptor (D3DR) have been clustered using a machine learning-based approach in the attempt to rationalize the efficacy change in four congeneric modulators. Accumulating extended MD trajectories of receptor-ligand complexes, we observed how the increase in ligand flexibility progressively destabilized the crystal structure of the inactivated receptor. To prospectively validate this model, a partial agonist was rationally designed based on structural insights and computational modeling, and eventually synthesized and tested. Results turned out to be in line with the predictions. This case study suggests that the investigation of ligand flexibility in the framework of extended MD simulations can assist and inform drug design strategies, highlighting its potential role as a powerful in silico counterpart to functional assays.

Structural basis of ligand selectivity and disease mutations in cysteinyl leukotriene receptors

Cysteinyl leukotriene G protein-coupled receptors CysLT₁ and CysLT₂ regulate pro-inflammatory responses associated with allergic disorders. While selective inhibition of CysLT₁R has been used for treating asthma and associated diseases for over two decades, CysLT₂R has recently started to emerge as a potential drug target against atopic asthma, brain injury and central nervous system disorders, as well as several types of cancer. Here, we describe four crystal structures of CysLT₂R in complex with three dual CysLT₁R/CysLT₂R antagonists. The reported structures together with the results of comprehensive mutagenesis and computer modeling studies shed light on molecular determinants of CysLTR ligand selectivity and specific effects of disease-related single nucleotide variants.

Cysteinyl leukotriene G protein-coupled receptors CysLT₁ and CysLT₂ regulate pro-inflammatory responses associated with allergic disorders. Here, authors describe four crystal structures of CysLT₂R in complex with three dual CysLT₁R/CysLT₂R antagonists, which shed light on CysLTR ligand selectivity.

Allosteric Modulation of Cannabinoid Receptor 1-Current Challenges and Future Opportunities

The cannabinoid receptor type 1 (CB1R), a G protein-coupled receptor (GPCR), plays an essential role in the control of many physiological processes such as hunger, memory loss, gastrointestinal activity, catalepsy, fear, depression, and chronic pain. Therefore, it is an attractive target for drug discovery to manage pain, neurodegenerative disorders, obesity, and substance abuse. However, the psychoactive adverse effects, generated by CB1R activation in the brain, limit the use of the orthosteric CB1R ligands as drugs. The discovery of CB1R allosteric modulators during the last decade provided new tools to target the CB1R. Moreover, application of the site-directed mutagenesis in combination with advanced physical methods, especially X-ray crystallography and computational modeling, has opened new horizons for understanding the complexity of the structure, function, and activity of cannabinoid receptors. In this paper, we present the latest advances in research on the CB1R, its allosteric modulation and allosteric ligands, and their translational potential. We focused on structural essentials of the cannabinoid 1 receptor- ligand (drug) interactions, as well as modes of CB1R signaling regulation.

The vasotocinergic system and its role in the regulation of stress in birds

The regulation of stress in birds includes a complex interaction of neural systems affecting the hypothalamic-pituitary-adrenal (HPA) axis. In addition to the hypothalamic paraventricular nucleus, a structure called the nucleus of the hippocampal commissure likewise affects the output of pituitary stress hormones and appears to be unique to avian species. Within the anterior pituitary, the avian V1a and V1b receptors were found in corticotropes. Based on our studies with central administration of hormones in the chicken, corticotropic releasing hormone (CRH) is a more potent ACTH secretagogue than arginine vasotocin (AVT). In contrast, when applied peripherally, AVT is more efficacious. Co-administration of AVT and CRH peripherally, resulted in a synergistic stimulation of corticosterone release. Data suggest receptor oligomerization as one possible mechanism. In birds, vasotocin receptors associated with stress responses include the V1a and V1b receptors. Three-dimensional, homology-based structural models of the avian V1aR were built to test agonists and antagonists for each receptor that were screened by molecular docking to map their binding sites on each receptor. Additionally, binding affinity values for each available peptide antagonist to the V1aR and V1bR were determined. An anterior pituitary primary culture system was developed to determine how effective each antagonist blocked the function of each receptor in culture when stimulated by a combination of AVT/CRH administration. Use of an antagonist in subsequent in vivo studies identified the V1aR in regulating food intake in birds. The V1aR was likewise found in circumventricular organs of the brain, suggesting a possible function in stress.

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

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

Vaccine Targeted Alpha 1D-Adrenergic Receptor for Hypertension

The α1-AR (α1 adrenergic receptor) blockers currently on the market cannot meet clinical needs because of low-selectivity for subtypes of α1-ARs, short half-life, and uncertain role in cardiovascular end point events. The study sought to find a vaccine specifically against α1D-AR (α1D-adrenergic receptor) for treating hypertension. A short peptide ADR-004 (cgiteeagy) belonging to α1D-AR was screened, and the ADRQβ-004 vaccine was produced and injected into spontaneously hypertensive rats model (including a short-term study, 10 weeks, and a long-term observation study, 39 weeks) and NG-nitro-l-arginine methyl ester + spontaneously hypertensive rats model (15 weeks). The antihypertensive effect and target organ protection of the ADRQβ-004 vaccine were carefully evaluated. The possible immune-mediated damage was detected in normal vaccinated Sprague Dawley rats. The ADR-004 peptide has perfect immunogenicity, and the ADRQβ-004 vaccine could induce strong antibody production. In the short-term study, the ADRQβ-004 vaccine averagely decreased the systolic blood pressure of spontaneously hypertensive rats up to 15 mm Hg and that of NG-nitro-l-arginine methyl ester+spontaneously hypertensive rats up to 29 mm Hg. In the long-term observation model, the antihypertensive effect of the ADRQβ-004 vaccine was quite stable, and the average decline of systolic blood pressure was 22 mm Hg. The ADRQβ-004 vaccine effectively prevented vascular structural remodeling, cardiac hypertrophy and fibrosis, and renal injury of hypertensive animals, superior to prazosin at renal level. Moreover, the ADRQβ-004 vaccine obviously downregulated the expression of α1D-AR, but not α1A-AR. Additionally, no significant immune-mediated damage was detected in immunized animals. The present results demonstrate that the ADRQβ-004 vaccine may provide a novel and promising method for the treatment of hypertension.

Molecular insights into the interaction of hemorphin and its targets

Hemorphins are atypical endogenous opioid peptides produced by the cleavage of hemoglobin beta chain. Several studies have reported the therapeutic potential of hemorphin in memory enhancement, blood regulation, and analgesia. However, the mode of interaction of hemorphin with its target remains largely elusive. The decapeptide LVV-hemorphin-7 is the most stable form of hemorphin. It binds with high affinity to mu-opioid receptors (MOR), angiotensin-converting enzyme (ACE) and insulin-regulated aminopeptidase (IRAP). In this study, computational methods were used extensively to elucidate the most likely binding pose of mammalian LVV-hemorphin-7 with the aforementioned proteins and to calculate the binding affinity. Additionally, alignment of mammalian hemorphin sequences showed that the hemorphin sequence of the camel harbors a variation - a Q > R substitution at position 8. This study also investigated the binding affinity and the interaction mechanism of camel LVV-hemorphin-7 with these proteins. To gain a better understanding of the dynamics of the molecular interactions between the selected targets and hemorphin peptides, 100 ns molecular dynamics simulations of the best-ranked poses were performed. Simulations highlighted major interactions between the peptides and key residues in the binding site of the proteins. Interestingly, camel hemorphin had a higher binding affinity and showed more interactions with all three proteins when compared to the canonical mammalian LVV-hemorphin-7. Thus, camel LVV-hemorphin-7 could be explored as a potent therapeutic agent for memory loss, hypertension, and analgesia.

Dual binding mode of "bitter sugars" to their human bitter taste receptor target

The 25 human bitter taste receptors (hTAS2Rs) are responsible for detecting bitter molecules present in food, and they also play several physiological and pathological roles in extraoral compartments. Therefore, understanding their ligand specificity is important both for food research and for pharmacological applications. Here we provide a molecular insight into the exquisite molecular recognition of bitter β-glycopyranosides by one of the members of this receptor subclass, hTAS2R16. Most of its agonists have in common the presence of a β-glycopyranose unit along with an extremely structurally diverse aglycon moiety. This poses the question of how hTAS2R16 can recognize such a large number of "bitter sugars". By means of hybrid molecular mechanics/coarse grained molecular dynamics simulations, here we show that the three hTAS2R16 agonists salicin, arbutin and phenyl-β-D-glucopyranoside interact with the receptor through a previously unrecognized dual binding mode. Such mechanism may offer a seamless way to fit different aglycons inside the binding cavity, while maintaining the sugar bound, similar to the strategy used by several carbohydrate-binding lectins. Our prediction is validated a posteriori by comparison with mutagenesis data and also rationalizes a wealth of structure-activity relationship data. Therefore, our findings not only provide a deeper molecular characterization of the binding determinants for the three ligands studied here, but also give insights applicable to other hTAS2R16 agonists. Together with our results for other hTAS2Rs, this study paves the way to improve our overall understanding of the structural determinants of ligand specificity in bitter taste receptors.

A Machine Learning Approach for the Discovery of Ligand-Specific Functional Mechanisms of GPCRs

G protein-coupled receptors (GPCRs) play a key role in many cellular signaling mechanisms, and must select among multiple coupling possibilities in a ligand-specific manner in order to carry out a myriad of functions in diverse cellular contexts. Much has been learned about the molecular mechanisms of ligand-GPCR complexes from Molecular Dynamics (MD) simulations. However, to explore ligand-specific differences in the response of a GPCR to diverse ligands, as is required to understand ligand bias and functional selectivity, necessitates creating very large amounts of data from the needed large-scale simulations. This becomes a Big Data problem for the high dimensionality analysis of the accumulated trajectories. Here we describe a new machine learning (ML) approach to the problem that is based on transforming the analysis of GPCR function-related, ligand-specific differences encoded in the MD simulation trajectories into a representation recognizable by state-of-the-art deep learning object recognition technology. We illustrate this method by applying it to recognize the pharmacological classification of ligands bound to the 5-HT_2A and D2 subtypes of class-A GPCRs from the serotonin and dopamine families. The ML-based approach is shown to perform the classification task with high accuracy, and we identify the molecular determinants of the classifications in the context of GPCR structure and function. This study builds a framework for the efficient computational analysis of MD Big Data collected for the purpose of understanding ligand-specific GPCR activity.

A benchmark study of loop modeling methods applied to G protein-coupled receptors

G protein-coupled receptors (GPCR) are important drug discovery targets. Despite progress, many GPCR structures have not yet been solved. For these targets, comparative modeling is used in virtual ligand screening to prioritize experimental efforts. However, the structure of extracellular loop 2 (ECL2) is often poorly predicted. This is significant due to involvement of ECL2 in ligand binding for many Class A GPCR. Here we examine the performance of loop modeling protocols available in the Rosetta (cyclic coordinate descent [CCD], KIC with fragments [KICF] and next generation KIC [NGK]) and Molecular Operating Environment (MOE) software suites (de novo search). ECL2 from GPCR crystal structures served as the structure prediction targets and were divided into four sets depending on loop length. Results suggest that KICF and NGK sampled and scored more loop models with sub-angstrom and near-atomic accuracy than CCD or de novo search for loops of 24 or fewer residues. None of the methods were able to sample loop conformations with near-atomic accuracy for the longest targets ranging from 25 to 32 residues based on 1000 models generated. For these long loop targets, increased conformational sampling is necessary. The strongly conserved disulfide bond between Cys3.25 and Cys45.50 in ECL2 proved an effective filter. Setting an upper limit of 5.1 Å on the S-S distance improved the lowest RMSD model included in the top 10 scored structures in Groups 1-4 on average between 0.33 and 1.27 Å. Disulfide bond formation and geometry optimization of ECL2 provided an additional incremental benefit in structure quality.

Multiscale simulations on human Frizzled and Taste2 GPCRs

Recently, molecular dynamics simulations, from all atom and coarse grained to hybrid methods bridging the two scales, have provided exciting functional insights into class F (Frizzled and Taste2) GPCRs (about 40 members in humans). Findings include: (i) The activation of one member of the Frizzled receptors (FZD4) involves a bending of transmembrane helix TM7 far larger than that in class A GPCRs. (ii) The affinity of an anticancer drug targeting another member (Smoothened receptor) decreases in a specific drug-resistant variant, because the mutation ultimately disrupts the binding cavity and affects TM6. (iii) A novel two-state recognition mechanism explains the very large agonist diversity for at least one member of the Taste2 GPCRs, hTAS2R46.

Automated discovery of GPCR bioactive ligands

While G-protein-coupled receptors (GPCRs) constitute the largest class of membrane proteins, structures and endogenous ligands of a large portion of GPCRs remain unknown. Because of the involvement of GPCRs in various signaling pathways and physiological roles, the identification of endogenous ligands as well as designing novel drugs is of high interest to the research and medical communities. Along with highlighting the recent advances in structure-based ligand discovery, including docking and molecular dynamics, this article focuses on the latest advances for automating the discovery of bioactive ligands using machine learning. Machine learning is centered around the development and applications of algorithms that can learn from data automatically. Such an approach offers immense opportunities for bioactivity prediction as well as quantitative structure-activity relationship studies. This review describes the most recent and successful applications of machine learning for bioactive ligand discovery, concluding with an outlook on deep learning methods that are capable of automatically extracting salient information from structural data as a promising future direction for rapid and efficient bioactive ligand discovery.

Influence of Cholesterol and Its Stereoisomers on Members of the Serotonin Receptor Family

Despite the ubiquity of cholesterol within the cell membrane, the mechanism by which it influences embedded proteins remains elusive. Numerous G-protein coupled receptors exhibit dramatic responses to membrane cholesterol with regard to the ligand-binding affinity and functional properties, including the 5-HT receptor family. Here, we use over 25 μs of unbiased atomistic molecular dynamics simulations to identify cholesterol interaction sites in the 5-HT_1B and 5-HT_2B receptors and evaluate their impact on receptor structure. Susceptibility to membrane cholesterol is shown to be subtype dependent and determined by the quality of interactions between the extracellular loops. Charged residues are essential for maintaining the arrangement of the extracellular surface in 5-HT_2B; in the absence of such interactions, the extracellular surface of the 5-HT_1B is malleable, populating a number of distinct conformations. Elevated cholesterol density near transmembrane helix 4 is considered to be conducive to the conformation of extracellular loop 2. Occupation of this site is also shown to be stereospecific, illustrated by differential behavior of nat-cholesterol isomers, ent- and epi-cholesterol. In simulations containing the endogenous agonist, serotonin, cholesterol binding at transmembrane helix 4 biases bound serotonin molecules toward an unexpected binding mode in the extended binding pocket. The results highlight the capability of membrane cholesterol to influence the mobility of the extracellular surface in the 5-HT₁ receptor family and manipulate the architecture of the extracellular ligand-binding pocket.

Mammalian class I odorant receptors exhibit a conserved vestibular-binding pocket

Odorant receptors represent the largest family of mammalian G protein-coupled receptors. Phylogenetically, they are split into two classes (I and II). By analyzing the entire subclass I odorant receptors sequences, we identified two class I-specific and highly conserved motifs. These are predicted to face each other at the extra-cellular portion of the transmembrane domain, forming a vestibular site at the entrance to the orthosteric-binding cavity. Molecular dynamics simulation combined with site-directed mutagenesis and in vitro functional assays confirm the functional role of this vestibular site in ligand-driven activation. Mutations at this part of the receptor differentially affect the receptor response to four agonists. Since this vestibular site is involved in ligand recognition, it could serve ligand design that targets specifically this sub-genome of mammalian odorant receptors.

Mutagenesis of GPR139 reveals ways to create gain or loss of function receptors

GPR139 is a Gq-coupled receptor activated by the essential amino acids L-tryptophan (L-Trp) and L-phenylalanine (L-Phe). We carried out mutagenesis studies of the human GPR139 receptor to identify the critical structural motifs required for GPR139 activation. We applied site-directed and high throughput random mutagenesis approaches using a double addition normalization strategy to identify novel GPR139 sequences coding receptors that have altered sensitivity to endogenous ligands. This approach resulted in GPR139 clones with gain-of-function, reduction-of-function or loss-of-function mutations. The agonist pharmacology of these mutant receptors was characterized and compared to wild-type receptor using calcium mobilization, radioligand binding, and protein expression assays. The structure-activity data were incorporated into a homology model which highlights that many of the gain-of-function mutations are either in or immediately adjacent to the purported orthosteric ligand binding site, whereas the loss-of-function mutations were largely in the intracellular G-protein binding area or were disrupters of the helix integrity. There were also some reduction-of-function mutations in the orthosteric ligand binding site. These findings may not only facilitate the rational design of novel agonists and antagonists of GPR139, but also may guide the design of transgenic animal models to study the physiological function of GPR139.

New Binding Sites, New Opportunities for GPCR Drug Discovery

Many central biological events rely on protein-ligand interactions. The identification and characterization of protein-binding sites for ligands are crucial for the understanding of functions of both endogenous ligands and synthetic drug molecules. G protein-coupled receptors (GPCRs) typically detect extracellular signal molecules on the cell surface and transfer these chemical signals across the membrane, inducing downstream cellular responses via G proteins or β-arrestin. GPCRs mediate many central physiological processes, making them important targets for modern drug discovery. Here, we focus on the most recent breakthroughs in finding new binding sites and binding modes of GPCRs and their potentials for the development of new medicines.

A Critical Analysis of Molecular Mechanisms Underlying Membrane Cholesterol Sensitivity of GPCRs

G protein-coupled receptors (GPCRs) are the largest and a diverse family of proteins involved in signal transduction across biological membranes. GPCRs mediate a wide range of physiological processes and have emerged as major targets for the development of novel drug candidates in all clinical areas. Since GPCRs are integral membrane proteins, regulation of their organization, dynamics, and function by membrane lipids, in particular membrane cholesterol, has emerged as an exciting area of research. Cholesterol sensitivity of GPCRs could be due to direct interaction of cholesterol with the receptor (specific effect). Alternately, GPCR function could be influenced by the effect of cholesterol on membrane physical properties (general effect). In this review, we critically analyze the specific and general mechanisms of the modulation of GPCR function by membrane cholesterol, taking examples from representative GPCRs. While evidence for both the proposed mechanisms exists, there appears to be no clear-cut distinction between these two mechanisms, and a combination of these mechanisms cannot be ruled out in many cases. We conclude that classifying the mechanism underlying cholesterol sensitivity of GPCR function merely into these two mutually exclusive classes could be somewhat arbitrary. A more holistic approach could be suitable for analyzing GPCR-cholesterol interaction.

Binding of histamine to the H1 receptor-a molecular dynamics study

Binding of histamine to the G-protein coupled histamine H₁ receptor plays an important role in the context of allergic reactions; however, no crystal structure of the resulting complex is available yet. To deduce the histamine binding site, we performed unbiased molecular dynamics (MD) simulations on a microsecond time scale, which allowed to monitor one binding event, in which particularly the residues of the extracellular loop 2 were involved in the initial recognition process. The final histamine binding pose in the orthosteric pocket is characterized by interactions with Asp107^3.32, Tyr108^3.33, Thr194^5.43, Asn198^5.46, Trp428^6.48, Tyr431^6.51, Phe432^6.52, and Phe435^6.55, which is in agreement with existing mutational data. The conformational stability of the obtained complex structure was subsequently confirmed in 2 μs equilibrium MD simulations, and a metadynamics simulation proved that the detected binding site represents an energy minimum. A complementary investigation of a D107A mutant, which has experimentally been shown to abolish ligand binding, revealed that this exchange results in a significantly weaker interaction and enhanced ligand dynamics. This finding underlines the importance of the electrostatic interaction between the histamine ammonium group and the side chain of Asp107^3.32 for histamine binding.

Molecular Determinants for Ligand Selectivity of the Cell-Free Synthesized Human Endothelin B Receptor

Extracellular domains of G-protein-coupled receptors act as initial molecular selectivity filters for subtype specific ligands and drugs. Chimeras of the human endothelin-B receptor containing structural units from the extracellular domains of the endothelin-A receptor were analyzed after their co-translational insertion into preformed nanodiscs. A short β-strand and a linker region in the second extracellular loop as well as parts of the extracellular N-terminal domain were identified as molecular discrimination sites for the endothelin-B receptor-selective agonists IRL1620, sarafotoxin 6c, 4Ala-ET-1 and ET-3, but not for the non-selective agonist ET-1 recognized by both endothelin receptors. A proposed second disulfide bridge in the endothelin-B receptor tethering the N-terminal domain with the third extracellular loop was not essential for ET-1 recognition and binding, but increased the receptor thermostability. We further demonstrate an experimental approach with cell-free synthesized engineered agonists to analyze the differential discrimination of peptide ligand topologies by the two endothelin receptors. The study is based on the engineering and cell-free insertion of G-protein-coupled receptors into defined membranes and may become interesting also for other targets as an alternative platform to reveal molecular details of ligand selectivity and ligand binding mechanisms.

Chemogenetic Approach Using Ni(II) Complex-Agonist Conjugates Allows Selective Activation of Class A G-Protein-Coupled Receptors

Investigating individual G-protein-coupled receptors (GPCRs) involved in various signaling cascades can unlock a myriad of invaluable physiological findings. One of the promising strategies for addressing the activity of each subtype of receptor is to design chemical turn-on switches on the target receptors. However, valid methods to selectively control class A GPCRs, the largest receptor family encoded in the human genome, remain limited. Here, we describe a novel approach to chemogenetically manipulate activity of engineered class A GPCRs carrying a His₄ tag, using metal complex-agonist conjugates (MACs). This manipulation is termed coordination tethering. With the assistance of coordination bonds, MACs showed 10-100-fold lower EC₅₀ values in the engineered receptors, compared with wild-type receptors. Such coordination tethering enabled selective activation of β₂-adrenoceptors and muscarinic acetylcholine receptors, without loss of natural receptor responses, in living mammalian cells, including primary cultured astrocytes. Our generalized, modular chemogenetic approach should facilitate more precise control and deeper understanding of individual GPCR signaling pathways in living systems.

Parathyroid Hormone Senses Extracellular Calcium To Modulate Endocrine Signaling upon Binding to the Family B GPCR Parathyroid Hormone 1 Receptor

Parathyroid hormone (PTH) binds to a family B G protein coupled receptor, parathyroid hormone 1 receptor (PTH1R). One of its functions is to regulate Ca2+ homeostasis in bone remodeling, during which Ca2+ can reach up to 40 mM. A truncated version of PTH, PTH(1-34), can fully activate PTH1R and has been used for osteoporosis treatments. Here, we used fluorescence anisotropy to examine the binding of PTH(1-34) to PTH1R purified in nanodiscs (PTH1R-ND) and found that the affinity increases 5-fold in the presence of 15 mM Ca2+. However, PTHrP(1-36), another truncated endogenous agonist for PTH1R, does not show this Ca2+ effect. Mutations of Glu19 and Glu22 in PTH(1-34) that are not conserved in PTHrP(1-36) largely abolished the Ca2+ effect. The results support that PTH(1-34) not only activates PTH1R but also uniquely senses Ca2+. This dual function of a peptide hormone is a novel observation that couples changes in extracellular environment with endocrine signaling. Understanding this can potentially reveal the complex role of PTH signaling in bone remodeling and improve the PTH(1-34) treatment for osteoporosis.

Pharmacological Properties and Biological Functions of the GPR17 Receptor, a Potential Target for Neuro-Regenerative Medicine

In 2006, cells heterologously expressing the "orphan" receptor GPR17 were shown to acquire responses to both uracil nucleotides and cysteinyl-leukotrienes, two families of signaling molecules accumulating in brain or heart as a result of hypoxic/traumatic injuries. In subsequent years, evidence of GPR17 key role in oligodendrogenesis and myelination has highlighted it as a "model receptor" for new therapies in demyelinating and neurodegenerative diseases. The apparently contrasting evidence in the literature about the role of GPR17 in promoting or inhibiting myelination can be due to its transient expression in the intermediate stages of differentiation, exerting a pro-differentiating function in early oligodendrocyte precursor cells (OPCs), and an inhibitory role in late stage maturing cells. Meanwhile, several papers extended the initial data on GPR17 pharmacology, highlighting a "promiscuous" behavior of this receptor; indeed, GPR17 is able to respond to other emergency signals like oxysterols or the pro-inflammatory cytokine SDF-1, underlying GPR17 ability to adapt its responses to changes of the surrounding extracellular milieu, including damage conditions. Here, we analyze the available literature on GPR17, in an attempt to summarize its emerging biological roles and pharmacological properties.

An acid-compatible co-polymer for the solubilization of membranes and proteins into lipid bilayer-containing nanoparticles

The fundamental importance of membrane proteins in drug discovery has meant that membrane mimetic systems for studying membrane proteins are of increasing interest. One such system has been the amphipathic, negatively charged poly(styrene-co-maleic acid) (SMA) polymer to form "SMA Lipid Particles" (SMALPs) which have been widely adopted to solubilize membrane proteins directly from the cell membrane. However, SMALPs are only soluble under basic conditions and precipitate in the presence of divalent cations required for many downstream applications. Here, we show that the positively charged poly(styrene-co-maleimide) (SMI) forms similar nanoparticles with comparable efficiency to SMA, whilst remaining functional at acidic pH and compatible with high concentrations of divalent cations. We have performed a detailed characterization of the performance of SMI that enables a direct comparison with similar data published for SMA. We also demonstrate that SMI is capable of extracting proteins directly from the cell membrane and can solubilize functional human G-protein coupled receptors (GPCRs) expressed in cultured HEK 293T cells. "SMILPs" thus provide an alternative membrane solubilization method that successfully overcomes some of the limitations of the SMALP method.

Identification of avian vasotocin receptor subtype-specific antagonists involved in the stress response of the chicken, Gallus gallus

Vasotocin 1a and 1b receptors (V1aR and V1bR) have been shown to play important roles in the neuroendocrine regulation of stress responses via the anterior pituitary (AP) of birds. To identify effective subtype-specific antagonists for the chicken V1aR (cV1aR) and cV1bR, potential antagonists to the mammalian V1R were screened against the cV1aR and cV1bR 3D structural models by molecular docking analysis with determination of binding pocket/amino acid residues involved in the interaction. The antagonistic effects of the selected ligands were examined by measuring pro-opiomelanocortin (POMC) heteronuclear RNA (hnPOMC) levels following the in vitro stress administration to primary chicken AP cells. Results of in silico analysis showed that the Manning compound and several other antagonists were bound to cV1bR with higher affinity than the natural agonist, arginine vasotocin (AVT). Similarities and differences in the antagonist-receptor binding interface with receptors were characterized for each ligand. Non-peptide mammalian V1bR antagonists, SSR-149415 and L-368899, were shown to be effective and had an additive effect in blocking POMC hnRNA expression in pituitary cell culture studies. SR-49059 antagonized the effect(s) of AVT/CRH on the downregulation of the cV1aR and the upregulation of the cCRH-R2 expression but not the cV1bR and cCRH-R1. The Manning compound antagonized the downregulation of cV1aR, cV1bR and cCRH-R1 and the upregulation of cCRH-R2 expression. The specificity of antagonists apparently resulted from unique differences in the interacting residues and their binding affinities. Collectively, these results provide valuable leads for future development of novel compounds capable of blocking or attenuating the AP stress response of avian species and perhaps other non-mammalian vertebrates as well.

A Novel G-Protein-Coupled Receptors Gene from Upland Cotton Enhances Salt Stress Tolerance in Transgenic Arabidopsis

Plants have developed a number of survival strategies which are significant for enhancing their adaptation to various biotic and abiotic stress factors. At the transcriptome level, G-protein-coupled receptors (GPCRs) are of great significance, enabling the plants to detect a wide range of endogenous and exogenous signals which are employed by the plants in regulating various responses in development and adaptation. In this research work, we carried out genome-wide analysis of target of Myb1 (TOM1), a member of the GPCR gene family. The functional role of TOM1 in salt stress tolerance was studied using a transgenic Arabidopsis plants over-expressing the gene. By the use of the functional domain PF06454, we obtained 16 TOM genes members in Gossypium hirsutum, 9 in Gossypium arboreum, and 11 in Gossypium raimondii. The genes had varying physiochemical properties, and it is significant to note that all the grand average of hydropathy (GRAVY) values were less than one, indicating that all are hydrophobic in nature. In all the genes analysed here, both the exonic and intronic regions were found. The expression level of Gh_A07G0747 (GhTOM) was significantly high in the transgenic lines as compared to the wild type; a similar trend in expression was observed in all the salt-related genes tested in this study. The study in epidermal cells confirmed the localization of the protein coded by the gene TOM1 in the plasma membrane. Analysis of anti-oxidant enzymes showed higher concentrations of antioxidants in transgenic lines and relatively lower levels of oxidant substances such as H₂O₂. The low malondialdehyde (MDA) level in transgenic lines indicated that the transgenic lines had relatively low level of oxidative damage compared to the wild types. The results obtained indicate that Gh_A07G0747 (GhTOM) can be a putative target gene for enhancing salt stress tolerance in plants and could be exploited in the future for the development of salt stress-tolerant cotton cultivars.

Peptide recognition, signaling and modulation of class B G protein-coupled receptors

Class B G protein-coupled receptors (GPCRs) are important drug targets in many human diseases, including type 2 diabetes, obesity, cardiovascular disease and neurodegeneration. Peptide hormones bind to these receptors through interactions with both the extracellular domain and transmembrane domain. Despite remarkable advances in structural studies of GPCRs, structural characterization of the full-length class B receptors remains extremely challenging due to their conformational complexity. The recently solved structures of class B GPCRs reveal the structural basis of peptide ligand recognition and modulation mechanisms of small molecule allosteric modulators. Furthermore, these structures provide essential insights into molecular mechanisms of class B GPCR signal transduction and modulation.

Independent Evolution of Strychnine Recognition by Bitter Taste Receptor Subtypes

The 25 human bitter taste receptors (hT2Rs) recognize thousands of structurally and chemically diverse bitter substances. The binding modes of human bitter taste receptors hT2R10 and hT2R46, which are responsible for strychnine recognition, were previously established using site-directed mutagenesis, functional assays, and molecular modeling. Here we construct a phylogenetic tree and reconstruct ancestral sequences of the T2R10 and T2R46 clades. We next analyze the binding sites in view of experimental data to predict their ability to recognize strychnine. This analysis suggests that the common ancestor of hT2R10 and hT2R46 is unlikely to bind strychnine in the same mode as either of its two descendants. Estimation of relative divergence times shows that hT2R10 evolved earlier than hT2R46. Strychnine recognition was likely acquired first by the earliest common ancestor of the T2R10 clade before the separation of primates from other mammals, and was highly conserved within the clade. It was probably independently acquired by the common ancestor of T2R43-47 before the homo-ape speciation, lost in most T2Rs within this clade, but enhanced in the hT2R46 after humans diverged from the rest of primates. Our findings suggest hypothetical strychnine T2R receptors in several species, and serve as an experimental guide for further study. Improved understanding of how bitter taste receptors acquire the ability to be activated by particular ligands is valuable for the development of sensors for bitterness and for potential toxicity.

Plasticity of the ligand binding pocket in the bitter taste receptor T2R7

Bitter taste receptors (T2Rs) are a group of 25 G protein-coupled receptors (GPCRs) in humans. The cognate agonists and the mechanism of ligand binding to the majority of the T2Rs remain unknown. Here we report the first structure-function analysis of T2R7 and study the ability of this receptor to bind to different agonists by site-directed mutagenesis. Screening of ligands for T2R7 in calcium based assays lead to the identification of novel compounds that activate this receptor. Quinine, diphenidol, dextromethorphan and diphenhydramine showed substantial activation of T2R7. Interestingly, these bitter compounds showed different pharmacological characteristics. To investigate the structural features in T2R7 that might contribute to the observed differences in agonist specificities, molecular model guided ligand docking and site-directed mutagenesis was pursued. Amino acids D65, D86, W89, N167, T169, W170, S181, T255 and E271 in the ligand-binding pocket were replaced and the mutants characterized pharmacologically. Our results suggest D86, S181 and W170 present on the extracellular side of transmembrane 3 (TM3), TM5 and in extracellular loop 2 (ECL2) are essential for agonist binding in T2R7. Mutations of these amino acids lead to loss-of-function. We also identified gain-of-function residues that are agonist specific. These results suggest that agonists bind at an extracellular site rather than deep within the TM core involving residues present in both ECL2 and TM helices in T2R7. Similar to majority of the Class A GPCRs, ECL2 in T2R7 plays a significant role in agonist binding and activation.

Predicting ligand binding poses for low-resolution membrane protein models: Perspectives from multiscale simulations

Membrane receptors constitute major targets for pharmaceutical intervention. Drug design efforts rely on the identification of ligand binding poses. However, the limited experimental structural information available may make this extremely challenging, especially when only low-resolution homology models are accessible. In these cases, the predictions may be improved by molecular dynamics simulation approaches. Here we review recent developments of multiscale, hybrid molecular mechanics/coarse-grained (MM/CG) methods applied to membrane proteins. In particular, we focus on our in-house MM/CG approach. It is especially tailored for G-protein coupled receptors, the largest membrane receptor family in humans. We show that our MM/CG approach is able to capture the atomistic details of the receptor/ligand binding interactions, while keeping the computational cost low by representing the protein frame and the membrane environment in a highly simplified manner. We close this review by discussing ongoing improvements and challenges of the current implementation of our MM/CG code.

Probing the cooperative mechanism of the μ–δ opioid receptor heterodimer by multiscale simulation

The activation-cooperativity of the μ–δ opioid receptor heterodimer was probed by multiscale simulation coupled with a protein structure network.

Outlining migrainous through dihydroergotamine–serotonin receptor interactions using quantum biochemistry

We present the electronic structure of the complex dihydroergotamine–serotonin receptor to unveil new medications to treat migraine and related diseases.

Leukotriene B4 receptor 2 gene polymorphism (rs1950504, Asp196Gly) leads to enhanced cell motility under low-dose ligand stimulation

Recently, single-nucleotide polymorphisms (SNPs) in G-protein-coupled receptors (GPCRs) have been suggested to contribute to physiopathology and therapeutic effects. Leukotriene B₄ receptor 2 (BLT2), a member of the GPCR family, plays a critical role in the pathogenesis of several inflammatory diseases, including cancer and asthma. However, no studies on BLT2 SNP effects have been reported to date. In this study, we demonstrate that the BLT2 SNP (rs1950504, Asp196Gly), a Gly-196 variant of BLT2 (BLT2 D196G), causes enhanced cell motility under low-dose stimulation of its ligands. In addition, we demonstrated that Akt activation and subsequent production of reactive oxygen species (ROS), both of which act downstream of BLT2, are also increased by BLT2 D196G in response to low-dose ligand stimulation. Furthermore, we observed that the ligand binding affinity of BLT2 D196G was enhanced compared with that of BLT2. Through homology modeling analysis, it was predicted that BLT2 D196G loses ionic interaction with R197, potentially resulting in increased agonist-receptor interaction. To the best of our knowledge, this report is the first to describe a SNP study on BLT2 and shows that BLT2 D196G enhances ligand sensitivity, thereby increasing cell motility in response to low-dose ligand stimulation.

Open Boundary Simulations of Proteins and Their Hydration Shells by Hamiltonian Adaptive Resolution Scheme

The recently proposed Hamiltonian adaptive resolution scheme (H-AdResS) allows the performance of molecular simulations in an open boundary framework. It allows changing, on the fly, the resolution of specific subsets of molecules (usually the solvent), which are free to diffuse between the atomistic region and the coarse-grained reservoir. So far, the method has been successfully applied to pure liquids. Coupling the H-AdResS methodology to hybrid models of proteins, such as the molecular mechanics/coarse-grained (MM/CG) scheme, is a promising approach for rigorous calculations of ligand binding free energies in low-resolution protein models. Toward this goal, here we apply for the first time H-AdResS to two atomistic proteins in dual-resolution solvent, proving its ability to reproduce structural and dynamic properties of both the proteins and the solvent, as obtained from atomistic simulations.

A monoclonal antibody raised against a thermo-stabilised β1-adrenoceptor interacts with extracellular loop 2 and acts as a negative allosteric modulator of a sub-set of β1-adrenoceptors expressed in stable cell lines

Recent interest has focused on antibodies that can discriminate between different receptor conformations. Here we have characterised the effect of a monoclonal antibody (mAb3), raised against a purified thermo-stabilised turkey β₁-adrenoceptor (β₁AR-m23 StaR), on β₁-ARs expressed in CHO-K1 or HEK 293 cells. Immunohistochemical and radioligand-binding studies demonstrated that mAb3 was able to bind to ECL2 of the tβ₁-AR, but not its human homologue. Specific binding of mAb3 to tβ₁-AR was inhibited by a peptide based on the turkey, but not the human, ECL2 sequence. Studies with [³H]-CGP 12177 demonstrated that mAb3 prevented the binding of orthosteric ligands to a subset (circa 40%) of turkey β₁-receptors expressed in both CHO K1 and HEK 293 cells. MAb3 significantly reduced the maximum specific binding capacity of [³H]-CGP-12177 without influencing its binding affinity. Substitution of ECL2 of tβ₁-AR with its human equivalent, or mutation of residues D186S, P187D, Q188E prevented the inhibition of [³H]-CGP 12177 binding by mAb3. MAb3 also elicited a negative allosteric effect on agonist-stimulated cAMP responses. The identity of the subset of turkey β₁-adrenoceptors influenced by mAb3 remains to be established but mAb3 should become an important tool to investigate the nature of β₁-AR conformational states and oligomeric complexes.

Understanding the molecular functions of the second extracellular loop (ECL2) of the calcitonin gene-related peptide (CGRP) receptor using a comprehensive mutagenesis approach

The extracellular loop 2 (ECL2) region is the most conserved of the three ECL domains in family B G protein-coupled receptors (GPCRs) and has a fundamental role in ligand binding and activation across the receptor super-family. ECL2 is fundamental for ligand-induced activation of the calcitonin gene related peptide (CGRP) receptor, a family B GPCR implicated in migraine and heart disease. In this study we apply a comprehensive targeted non-alanine substitution analysis method and molecular modelling to the functionally important residues of ECL2 to reveal key molecular interactions. We identified an interaction network between R274/Y278/D280/W283. These amino acids had the biggest reduction in signalling following alanine substitution analysis and comprise a group of basic, acidic and aromatic residues conserved in the wider calcitonin family of class B GPCRs. This study identifies key and varied constraints at each locus, including diverse biochemical requirements for neighbouring tyrosine residues and a W283H substitution that recovered wild-type (WT) signalling, despite the strictly conserved nature of the central ECL2 tryptophan and the catastrophic effects on signalling of W283A substitution. In contrast, while the distal end of ECL2 requires strict conservation of hydrophobicity or polarity in each position, mutation of these residues never has a large effect. This approach has revealed linked networks of amino acids, consistent with structural models of ECL2 and likely to represent a shared structural framework at an important ligand-receptor interface that is present across the family B GPCRs.

Membrane-induced organization and dynamics of the N-terminal domain of chemokine receptor CXCR1: insights from atomistic simulations

The CXC chemokine receptor 1 (CXCR1) is an important member of the G protein-coupled receptor (GPCR) family in which the extracellular N-terminal domain has been implicated in ligand binding and selectivity. The structure of this domain has not yet been elucidated due to its inherent dynamics, but experimental evidence points toward membrane-dependent organization and dynamics. To gain molecular insight into the interaction of the N-terminal domain with the membrane bilayer, we performed a series of microsecond time scale atomistic simulations of the N-terminal domain of CXCR1 in the presence and absence of POPC bilayers. Our results show that the peptide displays a high propensity to adopt a β-sheet conformation in the presence of the membrane bilayer. The interaction of the peptide with the membrane bilayer was found to be transient in our simulations. Interestingly, a scrambled peptide, containing the same residues in a randomly varying sequence, did not exhibit membrane-modulated structural dynamics. These results suggest that sequence-dependent electrostatics, modulated by the membrane, could play an important role in folding of the N-terminal domain. We believe that our results reinforce the emerging paradigm that cellular membranes could be important modulators of function of G protein-coupled receptors such as CXCR1.

Serotonin 2A receptor disulfide bridge integrity is crucial for ligand binding to different signalling states but not for its homodimerization

The serotonin 2_A (5-HT_2A) receptor is a G-protein coupled receptor (GPCR) with a conserved disulfide bridge formed by Cys¹⁴⁸ (transmembrane helix 3, TM3) and Cys²²⁷ (extracellular loop 2, ECL-2). We hypothesized that disulfide bridges may determine serotonin 5-HT_2A receptor functions such as receptor activation, functional selectivity and ligand recognition. We used the reducing agent dithiothreitol (DTT) to determine how the reduction of disulfide bridges affects radioligand binding, second messenger mobilization and receptor dimerization. A DTT-induced decrease in the number of binding sites (1190 ± 63.55 fmol/mg protein for control cells compared with 921.2 ± 60.84 fmol/mg protein for DTT-treated cells) as well as in the efficacy of both signalling pathways characterized was observed, although the affinity and potency were unchanged. Bioluminiscence resonance energy transfer (BRET) assays revealed the DTT treatment did not modify the homodimeric nature of serotonin 5-HT_2A receptors. In molecular dynamic simulations, the ECL-2 of the receptor with a broken cysteine bond adopts a wider variety of conformations, some of which protrude deeper into the receptor orthosteric binding pocket leading to collapse of the pocket. A shrunken binding pocket would be incapable of accommodating lysergic acid diethylamide (LSD). Our findings suggest that the decrease of efficacy may be due to disruption of disulfide bridge between TM3 and ECL-2. This reveals the integrity of the ECL-2 epitope, which should be explored in the development of novel ligands acting as allosteric modulators of serotonin 5-HT_2A receptors.

Structural insight into the activation of a class B G-protein-coupled receptor by peptide hormones in live human cells

The activation mechanism of class B G-protein-coupled receptors (GPCRs) remains largely unknown. To characterize conformational changes induced by peptide hormones, we investigated interactions of the class B corticotropin-releasing factor receptor type 1 (CRF1R) with two peptide agonists and three peptide antagonists obtained by N-truncation of the agonists. Surface mapping with genetically encoded photo-crosslinkers and pair-wise crosslinking revealed distinct footprints of agonists and antagonists on the transmembrane domain (TMD) of CRF1R and identified numerous ligand-receptor contact sites, directly from the intact receptor in live human cells. The data enabled generating atomistic models of CRF- and CRF(12-41)-bound CRF1R, further explored by molecular dynamics simulations. We show that bound agonist and antagonist adopt different folds and stabilize distinct TMD conformations, which involves bending of helices VI and VII around flexible glycine hinges. Conservation of these glycine hinges among all class B GPCRs suggests their general role in activation of these receptors.

DOI:http://dx.doi.org/10.7554/eLife.27711.001

Understanding the common themes and diverse roles of the second extracellular loop (ECL2) of the GPCR super-family

The extracellular loops (ECLs) of G protein-coupled receptors (GPCRs) can bind directly to docked orthosteric or allosteric ligands, they can contain transient contact points for ligand entry into the transmembrane (TM) bundle and they can regulate the activation of the receptor signalling pathways. Of the three ECLs, ECL2 is the largest and most structurally diverse reflecting its functional importance. This has been shown through biochemical techniques and has been supported by the many subsequent crystal structures of GPCRs bound to both agonists and antagonists. ECL2 shares common structural features between (and sometimes across) receptor sub-families and can facilitate ligand entry to the TM core or act directly as a surface of the ligand-binding pocket. Structural similarities seem to underpin common binding mechanisms; however, where these exist, variations in primary sequence ensure ligand-binding specificity. This review will compare current understanding of the structural themes and main functional roles of ECL2 in ligand binding, activation and regulation of the major families of GPCRs.

Distinct activation modes of the Relaxin Family Peptide Receptor 2 in response to insulin-like peptide 3 and relaxin

Relaxin family peptide receptor 2 (RXFP2) is a GPCR known for its role in reproductive function. It is structurally related to the human relaxin receptor RXFP1 and can be activated by human gene-2 (H2) relaxin as well as its cognate ligand insulin-like peptide 3 (INSL3). Both receptors possess an N-terminal low-density lipoprotein type a (LDLa) module that is necessary for activation and is joined to a leucine-rich repeat domain by a linker. This linker has been shown to be important for H2 relaxin binding and activation of RXFP1 and herein we investigate the role of the equivalent region of RXFP2. We demonstrate that the linker’s highly-conserved N-terminal region is essential for activation of RXFP2 in response to both ligands. In contrast, the linker is necessary for H2 relaxin, but not INSL3, binding. Our results highlight the distinct mechanism by which INSL3 activates RXFP2 whereby ligand binding mediates reorientation of the LDLa module by the linker region to activate the RXFP2 transmembrane domains in conjunction with the INSL3 A-chain. In contrast, relaxin activation of RXFP2 involves a more RXFP1-like mechanism involving binding to the LDLa-linker, reorientation of the LDLa module and activation of the transmembrane domains by the LDLa alone.

Outer membrane phospholipase A's roles in Helicobacter pylori acid adaptation

Background

The pH of the human gastric mucosa varies around 2.5 so that only bacteria with strong acidic stress tolerance can colonize it. The ulcer causing Helicobacter pylori thrives in the gastric mucosa. We analyse the roles of the key outer membrane protein OMPLA in its roles in acid tolerance.

Results

The homology model of Helicobacter pylori outer membrane phospholipase A (OMPLA) reveals a twelve stranded β-barrel with a pore that allows molecules to pass with a diameter up to 4 Å. Structure based multiple sequence alignments revealed the functional roles of many amino acids, and led to the suggestion that OMPLA has multiple functions. Besides its role as phospholipase it lets urea enter and ammonium exit the periplasm. Combined with an extensive literature study, our work leads to a comprehensive model for H. pylori’s acid tolerance. This model is based on the conversion of urea into ammonium, and it includes multiple roles for OMPLA and involves two hitherto little studied membrane channels in the OMPLA operon.

Conclusion

The three-dimensional model of OMPLA predicts a transmembrane pore that can aid H. pylori’s acid tolerance through urea influx and ammonium efflux. After urea passes through OMPLA into the periplasm, it passes through the pH-gated inner membrane channel UreI into the cytoplasm where urease hydrolyses it into NH₃ and CO₂. Most of the NH₃ becomes NH₄⁺ that is likely to need an inner membrane channel to reach the periplasm. Two genes that are co-regulated with OMPLA in gastric Helicobacter operons could aid this transport. The NH₄⁺ that might leave the cell through the OMPLA pore has been implicated in H. pylor’s pathogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13099-017-0184-y) contains supplementary material, which is available to authorized users.

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

Motivation

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

Results

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

Availability and Implementation

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

Supplementary information

Supplementary data are available at Bioinformatics online.