High-throughput modeling of human G-protein coupled receptors: amino acid sequence alignment, three-dimensional model building, and receptor library screening

Characterization of a new potent and long-lasting single chain peptide agonist of RXFP1 in cells and in vivo translational models

Despite beneficial effects in acute heart failure, the full therapeutic potential of recombinant relaxin-2 has been hampered by its short half-life and the need for intravenous administration limiting its use to intensive care units. A multiparametric optimization of the relaxin B-chain led to the identification of single chain lipidated peptide agonists of RXFP1 like SA10SC-RLX with subcutaneous bioavailability and extended half-life. SA10SC-RLX has sub nanomolar activity on cells expressing human RXFP1 and molecular modeling associated with the study of different RXFP1 mutants was used to decipher the mechanism of SA10SC-RLX interaction with RXFP1. Telemetry was performed in rat where SA10SC-RLX was able to engage RXFP1 after subcutaneous administration without tachyphylaxis after repeated dosing. Renal blood flow was then used as a translational model to evaluate RXFP1 activation. SA10SC-RLX increased renal blood flow and decreased renal vascular resistance in rats as reported for relaxin in humans. In conclusion, SA10SC-RLX mimics relaxin activity in in vitro and in vivo models of acute RXFP1 engagement. SA10SC-RLX represents a new class of long-lasting RXFP1 agonist, suitable for once daily subcutaneous administration in patients and potentially paving the way to new treatments for chronic fibrotic and cardiovascular diseases.

Molecular dynamics of the histamine H3 membrane receptor reveals different mechanisms of GPCR signal transduction

In this work, we studied the mechanisms of classical activation and inactivation of signal transduction by the histamine H3 receptor, a 7-helix transmembrane bundle G-Protein Coupled Receptor through long-time-scale atomistic molecular dynamics simulations of the receptor embedded in a hydrated double layer of dipalmitoyl phosphatidyl choline, a zwitterionic polysaturated ordered lipid. Three systems were prepared: the apo receptor, representing the constitutively active receptor; and two holo-receptors-the receptor coupled to the antagonist/inverse agonist ciproxifan, representing the inactive state of the receptor, and the receptor coupled to the endogenous agonist histamine and representing the active state of the receptor. An extensive analysis of the simulation showed that the three states of H3R present significant structural and dynamical differences as well as a complex behavior given that the measured properties interact in multiple and interdependent ways. In addition, the simulations described an unexpected escape of histamine from the orthosteric binding site, in agreement with the experimental modest affinities and rapid off-rates of agonists.

Lysophosphatidic Acid Receptor Agonism: Discovery of Potent Nonlipid Benzofuran Ethanolamine Structures

Lysophosphatidic acid (LPA) is the natural ligand for two phylogenetically distinct families of receptors (LPA_1-3, LPA_4-6) whose pathways control a variety of physiologic and pathophysiological responses. Identifying the benefit of balanced activation/repression of LPA receptors has always been a challenge because of the high lability of LPA and the limited availability of selective and/or stable agonists. In this study, we document the discovery of small benzofuran ethanolamine derivatives (called CpX and CpY) behaving as LPA_1-3 agonists. Initially found as rabbit urethra contracting agents, their elusive receptors were identified from [³⁵S]GTPγS-binding and β-arrestin2 recruitment investigations and then confirmed by [³H]CpX binding studies (urethra, hLPA_1-2 membranes). Both compounds induced a calcium response in hLPA_1-3 cells within a range of 0.4-1.5-log lower potency as compared with LPA. The contractions of rabbit urethra strips induced by these compounds perfectly matched binding affinities with values reaching the two-digit nanomolar level. The antagonist, KI16425, dose-dependently antagonized CpX-induced contractions in agreement with its affinity profile (LPA₁≥LPA₃>>LPA₂). The most potent agonist, CpY, doubled intraurethral pressure in anesthetized female rats at 3 µg/kg i.v. Alternatively, CpX was shown to inhibit human preadipocyte differentiation, a process totally reversed by KI16425. Together with original molecular docking data, these findings clearly established these molecules as potent agonists of LPA_1-3 and consolidated the pivotal role of LPA₁ in urethra/prostate contraction as well as in fat cell development. The discovery of these unique and less labile LPA_1-3 agonists would offer new avenues to investigate the roles of LPA receptors. SIGNIFICANCE STATEMENT: We report the identification of benzofuran ethanolamine derivatives behaving as potent selective nonlipid LPA_1-3 agonists and shown to alter urethra muscle contraction or preadipocyte differentiation. Unique at this level of potency, selectivity, and especially stability, compared with lysophosphatidic acid, they represent more appropriate tools for investigating the physiological roles of lysophosphatidic acid receptors and starting point for optimization of drug candidates for therapeutic applications.

Unsupervised Classification of G-Protein Coupled Receptors and Their Conformational States Using IChem Intramolecular Interaction Patterns

Over the past decade, the ever-growing structural information on G-protein coupled receptors (GPCRs) has revealed the three-dimensional (3D) characteristics of a receptor structure that is competent for G-protein binding. Structural markers are now commonly used to distinguish GPCR functional states, especially when analyzing molecular dynamics simulations. In particular, the position of the sixth helix within the seven transmembrane domains (TMs) is directly related to the coupling of the G-protein. Here, we show that the structural pattern defined by transmembrane intramolecular interactions (hydrogen bonds excluding backbone/backbone interactions, ionic bonds and aromatic interactions) is suitable for comparison of GPCR 3D structures and unsupervised distinction of the receptor states. First, we analyze a microsecond long molecular dynamic simulation of the human ß2-adrenergic receptor (ADRB2). Clustering of the 3D structures by pattern similarity identifies stable states which match the conformational classes defined by structural markers. Furthermore, the method directly spots the few state-specific interactions. Transforming pattern into graph, we extend the method to the comparison of different GPCRs. Clustering all GPCR experimentally determined structures by clique relative size first separates receptors, then their conformational states, thereby suggesting that the interaction patterns are specific of the receptor sequence and that the interaction signatures of conformational states are not shared across distant homologues.

Repertoires of G protein-coupled receptors for Ciona-specific neuropeptides

Neuropeptides play pivotal roles in various biological events in the nervous, neuroendocrine, and endocrine systems, and are correlated with both physiological functions and unique behavioral traits of animals. Elucidation of functional interaction between neuropeptides and receptors is a crucial step for the verification of their biological roles and evolutionary processes. However, most receptors for novel peptides remain to be identified. Here, we show the identification of multiple G protein-coupled receptors (GPCRs) for species-specific neuropeptides of the vertebrate sister group, Ciona intestinalis Type A, by combining machine learning and experimental validation. We developed an original peptide descriptor-incorporated support vector machine and used it to predict 22 neuropeptide-GPCR pairs. Of note, signaling assays of the predicted pairs identified 1 homologous and 11 Ciona-specific neuropeptide-GPCR pairs for a 41% hit rate: the respective GPCRs for Ci-GALP, Ci-NTLP-2, Ci-LF-1, Ci-LF-2, Ci-LF-5, Ci-LF-6, Ci-LF-7, Ci-LF-8, Ci-YFV-1, and Ci-YFV-3. Interestingly, molecular phylogenetic tree analysis revealed that these receptors, excluding the Ci-GALP receptor, were evolutionarily unrelated to any other known peptide GPCRs, confirming that these GPCRs constitute unprecedented neuropeptide receptor clusters. Altogether, these results verified the neuropeptide-GPCR pairs in the protochordate and evolutionary lineages of neuropeptide GPCRs, and pave the way for investigating the endogenous roles of novel neuropeptides in the closest relatives of vertebrates and the evolutionary processes of neuropeptidergic systems throughout chordates. In addition, the present study also indicates the versatility of the machine-learning-assisted strategy for the identification of novel peptide-receptor pairs in various organisms.

Engineering Salt Bridge Networks between Transmembrane Helices Confers Thermostability in G-Protein-Coupled Receptors

Introduction of specific point mutations has been an effective strategy in enhancing the thermostability of G-protein-coupled receptors (GPCRs). Our previous work showed that a specific residue position on transmembrane helix 6 (TM6) in class A GPCRs consistently yields thermostable mutants. The crystal structure of human chemokine receptor CCR5 also showed increased thermostability upon mutation of two positions, A233D^6.33 and K303E^7.59. With the goal of testing the transferability of these two thermostabilizing mutations in other chemokine receptors, we tested the mutations A237D^6.33 and R307E^7.59 in human CCR3 for thermostability and aggregation properties in detergent solution. Interestingly, the double mutant exhibited a 6-10-fold decrease in the aggregation propensity of the wild-type protein. This is in stark contrast to the two single mutants whose aggregation properties resemble the wild type (WT). Moreover, unlike in CCR5, the two single mutants separately showed no increase in thermostability compared to the wild-type CCR3, while the double-mutant A237D^6.33/R307E^7.59 confers an increase of 2.6 °C in the melting temperature compared to the WT. Extensive all-atom molecular dynamics (MD) simulations in detergent micelles show that a salt bridge network between transmembrane helices TM3, TM6, and TM7 that is absent in the two single mutants confers stability in the double mutant. The free energy surface of the double mutant shows conformational homogeneity compared to the single mutants. An annular n-dodecyl maltoside detergent layer packs tighter to the hydrophobic surface of the double-mutant CCR3 compared to the single mutants providing additional stability. The purification of other C-C chemokine receptors lacking such stabilizing residues may benefit from the incorporation of these two point mutations.

Using molecular dynamics for the refinement of atomistic models of GPCRs by homology modeling

Despite GPCRs sharing a common seven helix bundle, analysis of the diverse crystallographic structures available reveal specific features that might be relevant for ligand design. Despite the number of crystallographic structures of GPCRs steadily increasing, there are still challenges that hamper the availability of new structures. In the absence of a crystallographic structure, homology modeling remains one of the important techniques for constructing 3D models of proteins. In the present study we investigated the use of molecular dynamics simulations for the refinement of GPCRs models constructed by homology modeling. Specifically, we investigated the relevance of template selection, ligand inclusion as well as the length of the simulation on the quality of the GPCRs models constructed. For this purpose we chose the crystallographic structure of the rat muscarinic M3 receptor as reference and constructed diverse atomistic models by homology modeling, using different templates. Specifically, templates used in the present work include the human muscarinic M2; the more distant human histamine H1 and the even more distant bovine rhodopsin as shown in the GPCRs phylogenetic tree. We also investigated the use or not of a ligand in the refinement process. Hence, we conducted the refinement process of the M3 model using the M2 muscarinic as template with tiotropium or NMS docked in the orthosteric site and compared with the results obtained with a model refined without any ligand bound.

Invertebrate Gonadotropin-Releasing Hormone-Related Peptides and Their Receptors: An Update

Gonadotropin-releasing hormones (GnRHs) play pivotal roles in reproductive functions via the hypothalamus, pituitary, and gonad axis, namely, HPG axis in vertebrates. GnRHs and their receptors (GnRHRs) are likely to be conserved in invertebrate deuterostomes and lophotrochozoans. All vertebrate and urochordate GnRHs are composed of 10 amino acids, whereas protostome, echinoderm, and amphioxus GnRH-like peptides are 11- or 12-residue peptide containing two amino acids after an N-terminal pyro-Glu. In urochordates, Halocynthia roretzi GnRH gene encodes two GnRH peptide sequences, whereas two GnRH genes encode three different GnRH peptides in Ciona intestinalis. These findings indicate the species-specific diversification of GnRHs. Intriguingly, the major signaling pathway for GnRHRs is intracellular Ca²⁺ mobilization in chordates, echinoderms, and protostomes, whereas Ciona GnRHRs (Ci-GnRHRs) are endowed with multiple GnRHergic cAMP production pathways in a ligand-selective manner. Moreover, the ligand-specific modulation of signal transduction via heterodimerization among Ci-GnRHR paralogs suggests the species-specific development of fine-tuning of gonadal functions in ascidians. Echinoderm GnRH-like peptides show high sequence differences compared to those of protostome counterparts, leading to the difficulty in classification of peptides and receptors. These findings also show both the diversity and conservation of GnRH signaling systems in invertebrates. The lack of the HPG axis in invertebrates indicates that biological functions of GnRHs are not release of gonadotropins in current invertebrates and common ancestors of vertebrates and invertebrates. To date, authentic or putative GnRHRs have been characterized from various echinoderms and protostomes as well as chordates and the mRNAs have been found to be distributed not only reproductive organs but also other tissues. Collectively, these findings further support the notion that invertebrate GnRHs have biological roles other than the regulation of reproductive functions. Moreover, recent molecular phylogenetic analysis suggests that adipokinetic hormone (AKH), corazonin (CRZ), and AKH/CRZ-related peptide (ACP) belong to the GnRH superfamily but has led to the different classifications of these peptides and receptors using different datasets including the number of sequences and structural domains. In this review, we provide current knowledge of, and perspectives in, molecular basis and evolutionary aspects of the GnRH, AKH, CRZ, and ACP.

Structure of Class B GPCRs: new horizons for drug discovery

Class B GPCRs of the secretin family are important drug targets in many human diseases including diabetes, neurodegeneration, cardiovascular disease and psychiatric disorders. X-ray crystal structures for the glucagon receptor and corticotropin-releasing factor receptor 1 have now been published. In this review, we analyse the new structures and how they compare with each other and with Class A and F receptors. We also consider the differences in druggability and possible similarity in the activation mechanisms. Finally, we discuss the potential for the design of small-molecule modulators for these important targets in drug discovery. This new structural insight allows, for the first time, structure-based drug design methods to be applied to Class B GPCRs.

Comparative Modeling of Drug Target Proteins☆

In this perspective, we begin by describing the comparative protein structure modeling technique and the accuracy of the corresponding models. We then discuss the significant role that comparative prediction plays in drug discovery. We focus on virtual ligand screening against comparative models and illustrate the state-of-the-art by a number of specific examples.

Chemical genomics approach for GPCR-ligand interaction prediction and extraction of ligand binding determinants

Chemical genomics research has revealed that G-protein coupled receptors (GPCRs) interact with a variety of ligands and that a large number of ligands are known to bind GPCRs even with low transmembrane (TM) sequence similarity. It is crucial to extract informative binding region propensities from large quantities of bioactivity data. To address this issue, we propose a machine learning approach that enables identification of both chemical substructures and amino acid properties that are associated with ligand binding, which can be applied to virtual ligand screening on a GPCR-wide scale. We also address the question of how to select plausible negative noninteraction pairs based on a statistical approach in order to develop reliable prediction models for GPCR-ligand interactions. The key interaction sites estimated by our approach can be of great use not only for screening of active compounds but also for modification of active compounds with the aim of improving activity or selectivity.

Similarity between class A and class B G-protein-coupled receptors exemplified through calcitonin gene-related peptide receptor modelling and mutagenesis studies

Modelling class B G-protein-coupled receptors (GPCRs) using class A GPCR structural templates is difficult due to lack of homology. The plant GPCR, GCR1, has homology to both class A and class B GPCRs. We have used this to generate a class A–class B alignment, and by incorporating maximum lagged correlation of entropy and hydrophobicity into a consensus score, we have been able to align receptor transmembrane regions. We have applied this analysis to generate active and inactive homology models of the class B calcitonin gene-related peptide (CGRP) receptor, and have supported it with site-directed mutagenesis data using 122 CGRP receptor residues and 144 published mutagenesis results on other class B GPCRs. The variation of sequence variability with structure, the analysis of polarity violations, the alignment of group-conserved residues and the mutagenesis results at 27 key positions were particularly informative in distinguishing between the proposed and plausible alternative alignments. Furthermore, we have been able to associate the key molecular features of the class B GPCR signalling machinery with their class A counterparts for the first time. These include the [K/R]KLH motif in intracellular loop 1, [I/L]xxxL and KxxK at the intracellular end of TM5 and TM6, the NPXXY/VAVLY motif on TM7 and small group-conserved residues in TM1, TM2, TM3 and TM7. The equivalent of the class A DRY motif is proposed to involve Arg^2.39, His^2.43 and Glu^3.46, which makes a polar lock with T^6.37. These alignments and models provide useful tools for understanding class B GPCR function.

Comparison of fragments comprising the first two helices of the human Y4 and the yeast Ste2p G-protein-coupled receptors

Solution NMR techniques are used to determine the structure and the topology of micelle integration of a large fragment of the Y4 receptor, a human G-protein-coupled receptor, that contains the entire N-terminal domain plus the first two transmembrane (TM) segments. The structure calculations reveal that the putative TM helices are indeed helical to a large extent, but that interruptions of secondary structure occur close to internal polar or charged residues. This view is supported by (15)N relaxation data, amide-water exchange rates, and attenuations from micelle-integrating spin labels. No contacts between different helices are observed. This is in contrast to a similar TM1-TM2 fragment from the yeast Ste2p receptor for which locations of the secondary and the tertiary structure agreed well with the predictions from a homology model. The difference in structure is discussed in terms of principal biophysical properties of residues within central regions of the putative TM helices. Overall, using the biophysical scale of Wimley and White the TM regions of Ste2p display much more favorable free energies for membrane integration. Accordingly, the full secondary structure and the tertiary structure in TM1-TM2 of the Y4 receptor is likely to be formed only when tertiary contacts with other TM segments are created during folding of the receptor.

Calcitonin and calcitonin receptor-like receptors: common themes with family B GPCRs?

The calcitonin receptor (CTR) and calcitonin receptor-like receptor (CLR) are two of the 15 human family B (or Secretin-like) GPCRs. CTR and CLR are of considerable biological interest as their pharmacology is moulded by interactions with receptor activity-modifying proteins. They also have therapeutic relevance for many conditions, such as osteoporosis, diabetes, obesity, lymphatic insufficiency, migraine and cardiovascular disease. In light of recent advances in understanding ligand docking and receptor activation in both the family as a whole and in CLR and CTR specifically, this review reflects how applicable general family B GPCR themes are to these two idiosyncratic receptors. We review the main functional domains of the receptors; the N-terminal extracellular domain, the juxtamembrane domain and ligand interface, the transmembrane domain and the intracellular C-terminal domain. Structural and functional findings from the CLR and CTR along with other family B GPCRs are critically appraised to gain insight into how these domains may function. The ability for CTR and CLR to interact with receptor activity-modifying proteins adds another level of sophistication to these receptor systems but means careful consideration is needed when trying to apply generic GPCR principles. This review encapsulates current thinking in the realm of family B GPCR research by highlighting both conflicting and recurring themes and how such findings relate to two unusual but important receptors, CTR and CLR.

Topogenesis of heterologously expressed fragments of the human Y4 GPCR

Fragments of large membrane proteins have the potential to facilitate structural analysis by NMR, but their folding state remains a concern. Here we determined the quality of folding upon heterologous expression for a series of N- or C-terminally truncated fragments of the human Y4 G-protein coupled receptor, amounting to six different complementation pairs. As the individual fragments lack a specific function that could be used to ascertain proper folding, we instead assessed folding on a basic level by studying their membrane topology and by comparing it to well-established structural models of GPCRs. The topology of the fragments was determined using a reporter assay based on C-terminal green fluorescent protein- or alkaline phosphatase-fusions. N-terminal fusions to Lep or Mistic were used if a periplasmic orientation of the N-terminus of the fragments was expected based on predictions. Fragments fused to Mistic expressed at comparably high levels, whereas Lep fusions were produced to a much lower extent. Though none of the fragments exclusively adopted one orientation, often the correct topology predominated. In addition, systematic analysis of the fragment series suggested that the C-terminal half of the Y4 receptor is more important for adopting the correct topology than the N-terminal part. Using the detergent dodecylphosphocholine, selected fragments were solubilized from the membrane and proved sufficiently stable to allow purification. Finally, as a first step toward reconstituting a functional receptor from two fragments, we observed a physical interaction between complementing fragments pairs upon co-expression.

Structural determinants of allosteric antagonism at metabotropic glutamate receptor 2: mechanistic studies with new potent negative allosteric modulators

BACKGROUND AND PURPOSE

Altered glutamatergic neurotransmission is linked to several neurological and psychiatric disorders. Metabotropic glutamate receptor 2 (mGlu₂) plays an important role on the presynaptic control of glutamate release and negative allosteric modulators (NAMs) acting on mGlu₂/₃ receptors are under assessment for their potential as antidepressants, neurogenics and cognitive enhancers. Two new potent mGlu₂/₃ NAMs, RO4988546 and RO5488608, are described in this study and the allosteric binding site in the transmembrane (TM) domain of mGlu₂ is characterized.

EXPERIMENTAL APPROACH

Site directed mutagenesis, functional measurements and β₂-adrenoceptor-based modelling of mGlu₂ were employed to identify important molecular determinants of two new potent mGlu₂/₃ NAMs.

KEY RESULTS

RO4988546 and RO5488608 affected both [³H]-LY354740 agonist binding at the orthosteric site and the binding of a tritiated positive allosteric modulator (³H-PAM), indicating that NAMs and PAMs could have overlapping binding sites in the mGlu₂ TM domain. We identified eight residues in the allosteric binding pocket that are crucial for non-competitive antagonism of agonist-dependent activation of mGlu₂ and directly interact with the NAMs: Arg³·²⁸, Arg³·²⁹, Phe³·³⁶, His(E2.52) , Leu⁵·⁴³, Trp⁶·⁴⁸, Phe⁶·⁵⁵ and Val⁷·⁴³. The mGlu₂ specific residue His(E2.52) is likely to be involved in selectivity and residues located in the outer part of the binding pocket are more important for [³H]-LY354740 agonist binding inhibition, which is independent of the highly conserved Trp⁶·⁴⁸ residue.

CONCLUSIONS AND IMPLICATIONS

This is the first complete molecular investigation of the allosteric binding pocket of mGlu₂ and Group II mGluRs and provides new information on what determines mGlu₂ NAMs selective interactions and effects.

Structure-based discovery of allosteric modulators of two related class B G-protein-coupled receptors

Despite the availability of X-ray crystal structure data for several members of the G-protein-coupled receptor (GPCR) superfamily, structure-based discovery of GPCR ligands has been exclusively restricted to class A (rhodopsin-like) receptors. Herein we report the identification, by a docking-based virtual screening approach, of noncompetitive ligands for two related class B (secretin-like) GPCRs: the glucagon receptor (GLR) and the glucagon-like peptide 1 receptor (GLP-1R). Starting from a knowledge-based three-dimensional model of the GLR, a database of 1.9 million commercially available drug-like compounds was screened for chemical similarity to existing GLR noncompetitive antagonists and docked to the transmembrane cavity of the GLR; 23 compounds were then selected based on protein-ligand interaction fingerprints, and were then purchased and evaluated for in vitro binding to GLR and modulation of glucagon-induced cAMP release. Two of the 23 compounds inhibited the effect of glucagon in a dose-dependent manner, with one inhibitor exhibiting the same potency as L-168 049, a reference noncompetitive GLR antagonist, in a whole-cell-based functional assay. Interestingly, one virtual hit that was inactive at the GLR was shown to bind to GLP-1R and potentiate the response to the endogenous GLP-1 ligand.

Snooker: a structure-based pharmacophore generation tool applied to class A GPCRs

G-protein coupled receptors (GPCRs) are important drug targets for various diseases and of major interest to pharmaceutical companies. The function of individual members of this protein family can be modulated by the binding of small molecules at the extracellular side of the structurally conserved transmembrane (TM) domain. Here, we present Snooker, a structure-based approach to generate pharmacophore hypotheses for compounds binding to this extracellular side of the TM domain. Snooker does not require knowledge of ligands, is therefore suitable for apo-proteins, and can be applied to all receptors of the GPCR protein family. The method comprises the construction of a homology model of the TM domains and prioritization of residues on the probability of being ligand binding. Subsequently, protein properties are converted to ligand space, and pharmacophore features are generated at positions where protein ligand interactions are likely. Using this semiautomated knowledge-driven bioinformatics approach we have created pharmacophore hypotheses for 15 different GPCRs from several different subfamilies. For the beta-2-adrenergic receptor we show that ligand poses predicted by Snooker pharmacophore hypotheses reproduce literature supported binding modes for ∼75% of compounds fulfilling pharmacophore constraints. All 15 pharmacophore hypotheses represent interactions with essential residues for ligand binding as observed in mutagenesis experiments and compound selections based on these hypotheses are shown to be target specific. For 8 out of 15 targets enrichment factors above 10-fold are observed in the top 0.5% ranked compounds in a virtual screen. Additionally, prospectively predicted ligand binding poses in the human dopamine D3 receptor based on Snooker pharmacophores were ranked among the best models in the community wide GPCR dock 2010.

Residues within the transmembrane domain of the glucagon-like peptide-1 receptor involved in ligand binding and receptor activation: modelling the ligand-bound receptor

The C-terminal regions of glucagon-like peptide-1 (GLP-1) bind to the N terminus of the GLP-1 receptor (GLP-1R), facilitating interaction of the ligand N terminus with the receptor transmembrane domain. In contrast, the agonist exendin-4 relies less on the transmembrane domain, and truncated antagonist analogs (e.g. exendin 9-39) may interact solely with the receptor N terminus. Here we used mutagenesis to explore the role of residues highly conserved in the predicted transmembrane helices of mammalian GLP-1Rs and conserved in family B G protein coupled receptors in ligand binding and GLP-1R activation. By iteration using information from the mutagenesis, along with the available crystal structure of the receptor N terminus and a model of the active opsin transmembrane domain, we developed a structural receptor model with GLP-1 bound and used this to better understand consequences of mutations. Mutation at Y152 [transmembrane helix (TM) 1], R190 (TM2), Y235 (TM3), H363 (TM6), and E364 (TM6) produced similar reductions in affinity for GLP-1 and exendin 9-39. In contrast, other mutations either preferentially [K197 (TM2), Q234 (TM3), and W284 (extracellular loop 2)] or solely [D198 (TM2) and R310 (TM5)] reduced GLP-1 affinity. Reduced agonist affinity was always associated with reduced potency. However, reductions in potency exceeded reductions in agonist affinity for K197A, W284A, and R310A, while H363A was uncoupled from cAMP generation, highlighting critical roles of these residues in translating binding to activation. Data show important roles in ligand binding and receptor activation of conserved residues within the transmembrane domain of the GLP-1R. The receptor structural model provides insight into the roles of these residues.

Allosteric model of maraviroc binding to CC chemokine receptor 5 (CCR5)

Maraviroc is a nonpeptidic small molecule human immunodeficiency virus type 1 (HIV-1) entry inhibitor that has just entered the therapeutic arsenal for the treatment of patients. We recently demonstrated that maraviroc binding to the HIV-1 coreceptor, CC chemokine receptor 5 (CCR5), prevents it from binding the chemokine CCL3 and the viral envelope glycoprotein gp120 by an allosteric mechanism. However, incomplete knowledge of ligand-binding sites and the lack of CCR5 crystal structures have hampered an in-depth molecular understanding of how the inhibitor works. Here, we addressed these issues by combining site-directed mutagenesis (SDM) with homology modeling and docking. Six crystal structures of G-protein-coupled receptors were compared for their suitability for CCR5 modeling. All CCR5 models had equally good geometry, but that built from the recently reported dimeric structure of the other HIV-1 coreceptor CXCR4 bound to the peptide CVX15 (Protein Data Bank code 3OE0) best agreed with the SDM data and discriminated CCR5 from non-CCR5 binders in a virtual screening approach. SDM and automated docking predicted that maraviroc inserts deeply in CCR5 transmembrane cavity where it can occupy three different binding sites, whereas CCL3 and gp120 lie on distinct yet overlapped regions of the CCR5 extracellular loop 2. Data suggesting that the transmembrane cavity remains accessible for maraviroc in CCL3-bound and gp120-bound CCR5 help explain our previous observation that the inhibitor enhances dissociation of preformed ligand-CCR5 complexes. Finally, we identified residues in the predicted CCR5 dimer interface that are mandatory for gp120 binding, suggesting that receptor dimerization might represent a target for new CCR5 entry inhibitors.

Inhibitors of Helicobacter pylori protease HtrA found by 'virtual ligand' screening combat bacterial invasion of epithelia

Background

The human pathogen Helicobacter pylori (H. pylori) is a main cause for gastric inflammation and cancer. Increasing bacterial resistance against antibiotics demands for innovative strategies for therapeutic intervention.

Methodology/Principal Findings

We present a method for structure-based virtual screening that is based on the comprehensive prediction of ligand binding sites on a protein model and automated construction of a ligand-receptor interaction map. Pharmacophoric features of the map are clustered and transformed in a correlation vector (‘virtual ligand’) for rapid virtual screening of compound databases. This computer-based technique was validated for 18 different targets of pharmaceutical interest in a retrospective screening experiment. Prospective screening for inhibitory agents was performed for the protease HtrA from the human pathogen H. pylori using a homology model of the target protein. Among 22 tested compounds six block E-cadherin cleavage by HtrA in vitro and result in reduced scattering and wound healing of gastric epithelial cells, thereby preventing bacterial infiltration of the epithelium.

Conclusions/Significance

This study demonstrates that receptor-based virtual screening with a permissive (‘fuzzy’) pharmacophore model can help identify small bioactive agents for combating bacterial infection.

An inverse docking approach for identifying new potential anti-cancer targets

Inverse docking is a relatively new technique that has been used to identify potential receptor targets of small molecules. Our docking software package MDock is well suited for such an application as it is both computationally efficient, yet simultaneously shows adequate results in binding affinity predictions and enrichment tests. As a validation study, we present the first stage results of an inverse-docking study which seeks to identify potential direct targets of PRIMA-1. PRIMA-1 is well known for its ability to restore mutant p53's tumor suppressor function, leading to apoptosis in several types of cancer cells. For this reason, we believe that potential direct targets of PRIMA-1 identified in silico should be experimentally screened for their ability to inhibit cancer cell growth. The highest-ranked human protein of our PRIMA-1 docking results is oxidosqualene cyclase (OSC), which is part of the cholesterol synthetic pathway. The results of two followup experiments which treat OSC as a possible anti-cancer target are promising. We show that both PRIMA-1 and Ro 48-8071, a known potent OSC inhibitor, significantly reduce the viability of BT-474 and T47-D breast cancer cells relative to normal mammary cells. In addition, like PRIMA-1, we find that Ro 48-8071 results in increased binding of p53 to DNA in BT-474 cells (which express mutant p53). For the first time, Ro 48-8071 is shown as a potent agent in killing human breast cancer cells. The potential of OSC as a new target for developing anticancer therapies is worth further investigation.

The mouse eugenol odorant receptor: structural and functional plasticity of a broadly tuned odorant binding pocket

Molecular interactions of odorants with their olfactory receptors (ORs) are of central importance for the ability of the mammalian olfactory system to detect and discriminate a vast variety of odors with a limited set of receptors. How a particular OR binds and distinguishes different odorant molecules remains largely unknown on a structural basis. Here we investigated this question for the mouse eugenol receptor (mOR-EG). By screening a large odorant library, we discovered a wide range of chemical structures activating the receptor in heterologous mammalian cells. Potent agonists comprise (i) benzene, (ii) cyclohexane, or (iii) polycyclic structures substituted with alcohol, aldehyde, keto, ether, or esterified carboxylic groups. To detect those amino acids within the receptor that are in contact with a particular bound odorant molecule, we investigated how distinct mOR-EG point mutants were activated by the different odorant agonists found for the wild-type receptor. We identified 11 amino acids as a part of the receptor's ligand binding pocket. Molecular modeling predicted 10 of these residues in transmembrane helices TM3-TM6 and one in the extracellular loop between TM2 and TM3. These amino acids participate in odorant binding with variable importance depending on the type of odorant, revealing functional "fingerprints" of ligand-receptor interactions.

Evidence that interaction between conserved residues in transmembrane helices 2, 3, and 7 are crucial for human VPAC1 receptor activation

The VPAC(1) receptor belongs to family B of G protein-coupled receptors (GPCR-B) and is activated upon binding of the vasoactive intestinal peptide (VIP). Despite the recent determination of the structure of the N terminus of several members of this receptor family, little is known about the structure of the transmembrane (TM) region and about the molecular mechanisms leading to activation. In the present study, we designed a new structural model of the TM domain and combined it with experimental mutagenesis experiments to investigate the interaction network that governs ligand binding and receptor activation. Our results suggest that this network involves the cluster of residues Arg(188) in TM2, Gln(380) in TM7, and Asn(229) in TM3. This cluster is expected to be altered upon VIP binding, because Arg(188) has been shown previously to interact with Asp(3) of VIP. Several point mutations at positions 188, 229, and 380 were experimentally characterized and were shown to severely affect VIP binding and/or VIP-mediated cAMP production. Double mutants built from reciprocal residue exchanges exhibit strong cooperative or anticooperative effects, thereby indicating the spatial proximity of residues Arg(188), Gln(380), and Asn(229). Because these residues are highly conserved in the GPCR-B family, they can moreover be expected to have a general role in mediating function.

GPCR 3D homology models for ligand screening: lessons learned from blind predictions of adenosine A2a receptor complex

Proteins of the G-protein coupled receptor (GPCR) family present numerous attractive targets for rational drug design, but also a formidable challenge for identification and conformational modeling of their 3D structure. A recently performed assessment of blind predictions of adenosine A2a receptor (AA2AR) structure in complex with ZM241385 (ZMA) antagonist provided a first example of unbiased evaluation of the current modeling algorithms on a GPCR target with approximately 30% sequence identity to the closest structural template. Several of the 29 groups participating in this assessment exercise (Michino et al., doi: 10.1038/nrd2877) successfully predicted the overall position of the ligand ZMA in the AA2AR ligand binding pocket, however models from only three groups captured more than 40% the ligand-receptor contacts. Here we describe two of these top performing approaches, in which all-atom models of the AA2AR were generated by homology modeling followed by ligand guided backbone ensemble receptor optimization (LiBERO). The resulting AA2AR-ZMA models, along with the best models from other groups are assessed here for their vitual ligand screening (VLS) performance on a large set of GPCR ligands. We show that ligand guided optimization was critical for improvement of both ligand-receptor contacts and VLS performance as compared to the initial raw homology models. The best blindly predicted models performed on par with the crystal structure of AA2AR in selecting known antagonists from decoys, as well as from antagonists for other adenosine subtypes and AA2AR agonists. These results suggest that despite certain inaccuracies, the optimized homology models can be useful in the drug discovery process.

Phe369(7.38) at human 5-HT(7) receptors confers interspecies selectivity to antagonists and partial agonists

BACKGROUND AND PURPOSE

Human and rat 5-HT(7) receptors were studied with a particular emphasis on the molecular interactions involved in ligand binding, searching for an explanation to the interspecies selectivity observed for a set of compounds. We performed affinity studies, molecular modelling and site-directed mutagenesis, with special focus on residue Phe(7.38) of the human 5-HT(7) receptor [Cys(7.38) in rat].

EXPERIMENTAL APPROACH

Competition binding studies were performed for seven 5-HT(7) receptor ligands at three different 5-HT(7) receptors. The functional behaviour was evaluated by measuring 5-carboxytryptamine-stimulated cAMP production. Computational simulations were carried out to explore the structural bases in ligand binding observed for these compounds.

KEY RESULTS

Competition experiments showed a remarkable selectivity for the human receptor when compared with the rat receptor. These results indicate that mutating Cys to Phe at position 7.38 profoundly affects the binding affinities at the 5-HT(7) receptor. Computational simulations provide a structural interpretation for this key finding. Pharmacological characterization of compounds mr25020, mr25040 and mr25053 revealed a competitive antagonistic behaviour. Compounds mr22423, mr22433, mr23284 and mr25052 behaved as partial agonists.

CONCLUSIONS AND IMPLICATIONS

We propose that the interspecies difference in binding affinities observed for the compounds at human and rat 5-HT(7) receptors is due to the nature of the residue at position 7.38. Our molecular modelling simulations suggest that Phe(7.38) in the human receptor is integrated in the hydrophobic pocket in the central part of the binding site [Phe(6.51)-Phe(6.52)] and allows a tighter binding of the ligands when compared with the rat receptor.

Molecular determinants of ligand binding to H4R species variants

The histamine H(4) receptor (H(4)R) is the latest identified histamine receptor to emerge as a potential drug target for inflammatory diseases. Animal models are employed to validate this potential drug target. Concomitantly, various H(4)R orthologs have been cloned, including the human, mouse, rat, guinea pig, monkey, pig, and dog H(4)Rs. In this article, we expressed all these H(4)R orthologs in human embryonic kidney 293T cells and compared their interactions with currently used standard H(4)R ligands, including the H(4)R agonists histamine, 4-methylhistamine, guanidinylethyl isothiourea (VUF 8430), the H(4)R antagonists 1-[(5-chloro-1H-indol-2-yl)carbonyl]-4-methylpiperazine (JNJ 7777120) and [(5-chloro-1H-benzimidazol-2-yl)carbonyl]-4-methylpiperazine (VUF 6002), and the inverse H(4)R agonist thioperamide. Most of the evaluated ligands display significantly different affinities at the different H(4)R orthologs. These "natural mutants" of H(4)R were used to study ligand-receptor interactions by using chimeric human-pig-human and pig-human-pig H(4)R proteins and site-directed mutagenesis. Our results are a useful reference for ligand selection for studies in animal models of diseases and offer new insights in the understanding of H(4)R-ligand receptor interactions.

Increasing selectivity of CC chemokine receptor 8 antagonists by engineering nondesolvation related interactions with the intended and off-target binding sites

The metabolic stability and selectivity of a series of CCR8 antagonists against binding to the hERG ion channel and cytochrome Cyp2D6 are studied by principal component analysis. It is demonstrated that an efficient way of increasing metabolic stability and selectivity of this series is to decrease compound lipophilicity by engineering nondesolvation related attractive interactions with CCR8, as rationalized by three-dimensional receptor models. Although such polar interactions led to increased compound selectivity, such a strategy could also jeopardize the DMPK profile of compounds. However, once increased potency is found, the lipophilicity can be readjusted by engineering hydrophobic substituents that fit to CCR8 but do not fit to hERG. Several such lipophilic fragments are identified by two-dimensional fragment-based QSAR analysis. Electrophysiological measurements and site-directed mutagenesis studies indicated that the repulsive interactions of these fragments with hERG are caused by steric hindrances with residue F656.

Toward the discovery of vaccine adjuvants: coupling in silico screening and in vitro analysis of antagonist binding to human and mouse CCR4 receptors

Background

Adjuvants enhance or modify an immune response that is made to an antigen. An antagonist of the chemokine CCR4 receptor can display adjuvant-like properties by diminishing the ability of CD4+CD25+ regulatory T cells (Tregs) to down-regulate immune responses.

Methodology

Here, we have used protein modelling to create a plausible chemokine receptor model with the aim of using virtual screening to identify potential small molecule chemokine antagonists. A combination of homology modelling and molecular docking was used to create a model of the CCR4 receptor in order to investigate potential lead compounds that display antagonistic properties. Three-dimensional structure-based virtual screening of the CCR4 receptor identified 116 small molecules that were calculated to have a high affinity for the receptor; these were tested experimentally for CCR4 antagonism. Fifteen of these small molecules were shown to inhibit specifically CCR4-mediated cell migration, including that of CCR4⁺ Tregs.

Significance

Our CCR4 antagonists act as adjuvants augmenting human T cell proliferation in an in vitro immune response model and compound SP50 increases T cell and antibody responses in vivo when combined with vaccine antigens of Mycobacterium tuberculosis and Plasmodium yoelii in mice.

Dual activities of odorants on olfactory and nuclear hormone receptors

We have screened an odorant compound library and discovered molecules acting as chemical signals that specifically activate both G-protein-coupled olfactory receptors (ORs) on the cell surface of olfactory sensory neurons and the human nuclear estrogen receptor alpha (ER) involved in transcriptional regulation of cellular differentiation and proliferation in a wide variety of tissues. Hence, these apparent dual active odorants induce distinct signal transduction pathways at different subcellular localizations, which affect both neuronal signaling, resulting in odor perception, and the ER-dependent transcriptional control of specific genes. We demonstrate these effects using fluorescence-based in vitro and cellular assays. Among these odorants, we have identified synthetic sandalwood compounds, an important class of molecules used in the fragrance industry. For one estrogenic odorant we have also identified the cognate OR. This prompted us to compare basic molecular recognition principles of odorants on the two structurally and apparent functionally non-related receptors using computational modeling in combination with functional assays. Faced with the increasing evidence that ORs may perform chemosensory functions in a number of tissues outside of the nasal olfactory epithelium, the unraveling of these molecular ligand-receptor interaction principles is of critical importance. In addition the evidence that certain olfactory sensory neurons naturally co-express ORs and ERs may provide a direct functional link between the olfactory and hormonal systems in humans. Our results are therefore useful for defining the structural and functional characteristics of ER-specific odorants and the role of odorant molecules in cellular processes other than olfaction.

On the applicability of GPCR homology models to computer-aided drug discovery: a comparison between in silico and crystal structures of the beta2-adrenergic receptor

The publication of the crystal structure of the beta2-adrenergic receptor (beta2-AR) proved that G protein-coupled receptors (GPCRs) share a structurally conserved rhodopsin-like 7TM core. Here, to probe to which extent realistic GPCR structures can be recreated through modeling, carazolol was docked at two rhodopsin-based homology models of the human beta 2-AR. The first featured a rhodopsin-like second extracellular loop, which interfered with ligand docking and with the orientation of several residues in the binding pocket. The second featured a second extracellular loop built completely de novo, which afforded a more accurate model of the binding pocket and a better docking of the ligand. Furthermore, incorporating available biochemical and computational data to the model by correcting the conformation of a single residue lining the binding pocket --Phe290(6.52)--, resulted in significantly improved docking poses. These results support the applicability of GPCR modeling to the design of site-directed mutagenesis experiments and to drug discovery.

Computational chemistry approaches to drug discovery in signal transduction

The advent of therapeutic strategies aimed at targeting specific macromolecular components of deregulated signaling pathways associated with particular disease states has given rise to the idea that it should be possible to design ligands as drug candidates to these targets from first principles. This concept has been beckoning for a long time but structure-based ligand design only became feasible once it was possible to determine the 3-D structures of molecular targets at atomic resolution. However, structure-based design turned out to be difficult, chiefly because under physiological conditions both receptors and ligands are not static but they behave dynamically. While it is possible to design ligands with high steric and electronic complementarity to a receptor site, it is always uncertain how biologically relevant the assumed conformations of both ligand and receptor actually are. The fact that it remains beyond our current abilities to predict with sufficient accuracy the affinity between hypothetical ligand and receptor poses is in part connected with this problem and continues to confound the reliable prediction of drug-like ligands for therapeutic targets. Nevertheless, significant progress has been made and so-called virtual screening methods that use computational methods to dock candidate ligands into receptor sites and to score the resulting complexes are now used routinely as one of the components in drug discovery screening campaigns. Here an overview is given of the underlying principles, implementations, and applications of structure-guided computational design technologies. Although the emphasis is on receptor-based strategies, mention will also be made of some of the more established ligand-based approaches, such as similarity analyses and quantitative structure-activity relationship methods.

Molecular modeling of the second extracellular loop of G-protein coupled receptors and its implication on structure-based virtual screening

The current study describes the validation of high-throughput modeling procedures for the construction of the second extracellular loop (ecl2) of all nonolfactory human G Protein-coupled receptors. Our modeling flowchart is based on the alignment of essential residues determining the particular ecl2 fold observed in the bovine rhodopsin (bRho) crystal structure. For a set of GPCR targets, the dopamine D2 receptor (DRD2), adenosine A3 receptor (AA3R), and the thromboxane A2 receptor (TA2R), the implications of including ecl2 atomic coordinates is evaluated in terms of structure-based virtual screening accuracy: the suitability of the 3D models to distinguish between known antagonists and randomly chosen decoys using automated docking approaches. The virtual screening results of different models describing increasingly exhaustive receptor representations (seven helices only, seven helices and ecl2 loop, full model) have been compared. Explicit modeling of the ecl2 loop was found to be important in only one of three test cases whereas a loopless model was shown to be accurate enough in the two other receptors. An exhaustive comparison of ecl2 loops of 365 receptors to that of bRho suggests that explicit ecl2 loop modeling should be reserved to receptors where loop building can be guided by experimental restraints.