A chemogenomic analysis of the transmembrane binding cavity of human G-protein-coupled receptors

GPR101: Modeling a constitutively active receptor linked to X-linked acrogigantism

GPR101 is a G protein-coupled receptor (GPCR) implicated in a rare form of genetic gigantism known as X-linked acrogigantism, or X-LAG. In particular, X-LAG patients harbor microduplications in the long arm of the X-chromosome that invariably include the GPR101 gene. Duplications of the GPR101 gene lead to the formation of a new chromatin domain that causes over-expression of the receptor in the pituitary tumors of the patients. Notably, GPR101 is a constitutively active receptor, which stimulates cells to produce the second messenger cyclic AMP (cAMP) in the absence of ligands. Moreover, GPR101 was recently reported to constitutively activate not only the cAMP pathway via Gs, but also other G protein subunits (Gq/11 and G12/13). Hence, chemicals that block the constitutive activity of GPR101, known as inverse agonists, have the potential to be useful for the development of pharmacological tools for the treatment of X-LAG. In this study, we provide structural insights into the putative structure of GPR101 based on in-house built homology models, as well as third party models based on the machine learning methods AlphaFold and AlphaFold-Multistate. Moreover, we report a molecular dynamics study, meant to further probe the constitutive activity of GPR101. Finally, we provide a structural comparison with the closest GPCRs, which suggests that GPR101 does not share their natural ligands. While this manuscript was under review, cryo-electron microscopy structures of GPR101 were reported. These structures are expected to enable computer-aided ligand discovery efforts targeting GPR101.

Improving the Modeling of Extracellular Ligand Binding Pockets in RosettaGPCR for Conformational Selection

G protein-coupled receptors (GPCRs) are the largest class of drug targets and undergo substantial conformational changes in response to ligand binding. Despite recent progress in GPCR structure determination, static snapshots fail to reflect the conformational space of putative binding pocket geometries to which small molecule ligands can bind. In comparative modeling of GPCRs in the absence of a ligand, often a shrinking of the orthosteric binding pocket is observed. However, the exact prediction of the flexible orthosteric binding site is crucial for adequate structure-based drug discovery. In order to improve ligand docking and guide virtual screening experiments in computer-aided drug discovery, we developed RosettaGPCRPocketSize. The algorithm creates a conformational ensemble of biophysically realistic conformations of the GPCR binding pocket between the TM bundle, which is consistent with a knowledge base of expected pocket geometries. Specifically, tetrahedral volume restraints are defined based on information about critical residues in the orthosteric binding site and their experimentally observed range of Cα-Cα-distances. The output of RosettaGPCRPocketSize is an ensemble of binding pocket geometries that are filtered by energy to ensure biophysically probable arrangements, which can be used for docking simulations. In a benchmark set, pocket shrinkage observed in the default RosettaGPCR was reduced by up to 80% and the binding pocket volume range and geometric diversity were increased. Compared to models from four different GPCR homology model databases (RosettaGPCR, GPCR-Tasser, GPCR-SSFE, and GPCRdb), the here-created models showed more accurate volumes of the orthosteric pocket when evaluated with respect to the crystallographic reference structure. Furthermore, RosettaGPCRPocketSize was able to generate an improved realistic pocket distribution. However, while being superior to other homology models, the accuracy of generated model pockets was comparable to AlphaFold2 models. Furthermore, in a docking benchmark using small-molecule ligands with a higher molecular weight between 400 and 700 Da, a higher success rate in creating native-like binding poses was observed. In summary, RosettaGPCRPocketSize can generate GPCR models with realistic orthosteric pocket volumes, which are useful for structure-based drug discovery applications.

Transfer of knowledge from model organisms to evolutionarily distant non-model organisms: The coral Pocillopora damicornis membrane signaling receptome

With the ease of gene sequencing and the technology available to study and manipulate non-model organisms, the extension of the methodological toolbox required to translate our understanding of model organisms to non-model organisms has become an urgent problem. For example, mining of large coral and their symbiont sequence data is a challenge, but also provides an opportunity for understanding functionality and evolution of these and other non-model organisms. Much more information than for any other eukaryotic species is available for humans, especially related to signal transduction and diseases. However, the coral cnidarian host and human have diverged over 700 million years ago and homologies between proteins in the two species are therefore often in the gray zone, or at least often undetectable with traditional BLAST searches. We introduce a two-stage approach to identifying putative coral homologues of human proteins. First, through remote homology detection using Hidden Markov Models, we identify candidate human homologues in the cnidarian genome. However, for many proteins, the human genome alone contains multiple family members with similar or even more divergence in sequence. In the second stage, therefore, we filter the remote homology results based on the functional and structural plausibility of each coral candidate, shortlisting the coral proteins likely to have conserved some of the functions of the human proteins. We demonstrate our approach with a pipeline for mapping membrane receptors in humans to membrane receptors in corals, with specific focus on the stony coral, P. damicornis. More than 1000 human membrane receptors mapped to 335 coral receptors, including 151 G protein coupled receptors (GPCRs). To validate specific sub-families, we chose opsin proteins, representative GPCRs that confer light sensitivity, and Toll-like receptors, representative non-GPCRs, which function in the immune response, and their ability to communicate with microorganisms. Through detailed structure-function analysis of their ligand-binding pockets and downstream signaling cascades, we selected those candidate remote homologues likely to carry out related functions in the corals. This pipeline may prove generally useful for other non-model organisms, such as to support the growing field of synthetic biology.

Mind the Gap—Deciphering GPCR Pharmacology Using 3D Pharmacophores and Artificial Intelligence

G protein-coupled receptors (GPCRs) are amongst the most pharmaceutically relevant and well-studied protein targets, yet unanswered questions in the field leave significant gaps in our understanding of their nuanced structure and function. Three-dimensional pharmacophore models are powerful computational tools in in silico drug discovery, presenting myriad opportunities for the integration of GPCR structural biology and cheminformatics. This review highlights success stories in the application of 3D pharmacophore modeling to de novo drug design, the discovery of biased and allosteric ligands, scaffold hopping, QSAR analysis, hit-to-lead optimization, GPCR de-orphanization, mechanistic understanding of GPCR pharmacology and the elucidation of ligand–receptor interactions. Furthermore, advances in the incorporation of dynamics and machine learning are highlighted. The review will analyze challenges in the field of GPCR drug discovery, detailing how 3D pharmacophore modeling can be used to address them. Finally, we will present opportunities afforded by 3D pharmacophore modeling in the advancement of our understanding and targeting of GPCRs.

Taste GPCRs and their ligands

Taste GPCRs are expressed in taste buds on the tongue and play a key role in food choice and consumption. They are also expressed extra-orally, with various physiological roles that are currently under study. Unraveling the roles of these receptors relies on the knowledge of their ligands. Combining sensory, cell-based and computational approaches enabled the discovery of numerous agonists and several antagonists. Here we provide a short overview of taste receptor families, main recent methods for ligands discovery, and current sources of information about known ligands. The future directions that are likely to impact the taste GPCR field include focus on ligand interactions with naturally occurring polymorphisms, as well as harnessing the power of CryoEM and of multiple signaling readout techniques.

G-protein-coupled receptor P2Y10 facilitates chemokine-induced CD4 T cell migration through autocrine/paracrine mediators

G-protein-coupled receptors (GPCRs), especially chemokine receptors, play a central role in the regulation of T cell migration. Various GPCRs are upregulated in activated CD4 T cells, including P2Y10, a putative lysophospholipid receptor that is officially still considered an orphan GPCR, i.e., a receptor with unknown endogenous ligand. Here we show that in mice lacking P2Y10 in the CD4 T cell compartment, the severity of experimental autoimmune encephalomyelitis and cutaneous contact hypersensitivity is reduced. P2Y10-deficient CD4 T cells show normal activation, proliferation and differentiation, but reduced chemokine-induced migration, polarization, and RhoA activation upon in vitro stimulation. Mechanistically, CD4 T cells release the putative P2Y10 ligands lysophosphatidylserine and ATP upon chemokine exposure, and these mediators induce P2Y10-dependent RhoA activation in an autocrine/paracrine fashion. ATP degradation impairs RhoA activation and migration in control CD4 T cells, but not in P2Y10-deficient CD4 T cells. Importantly, the P2Y10 pathway appears to be conserved in human T cells. Taken together, P2Y10 mediates RhoA activation in CD4 T cells in response to auto-/paracrine-acting mediators such as LysoPS and ATP, thereby facilitating chemokine-induced migration and, consecutively, T cell-mediated diseases.

Free fatty acid receptor 1: a ray of hope in the therapy of type 2 diabetes mellitus

Free fatty acid receptor 1 (FFAR1) is a G-protein coupled receptor with prominent expression on pancreatic beta cells, bones, intestinal cells as well as the nerve cells. This receptor mediates a multitude of functions in the body including release of incretins, secretion of insulin as well as sensation of pain. Since FFAR1 causes secretion of insulin and regulates glucose metabolism, efforts were made to unfold its structure followed by discovering agonists for the receptor and the utilization of these agonists in the therapy of type 2 diabetes mellitus. Development of such functional FFAR1 agonists is a necessity because the currently available therapy for type 2 diabetes mellitus has numerous drawbacks, of which, the major one is hypoglycemia. Since the most prominent effect of the FFAR1 agonists is on glucose concentration in the body, so the major research is focused on treating type 2 diabetes mellitus, though the agonists could benefit other metabolic disorders and neurological disorders as well. The agonists developed so far had one major limitation, i.e., hepatotoxicity. Although, the only agonist that could reach phase 3 clinical trials was TAK-875 developed by Takeda Pharmaceuticals but it was also withdrawn due to toxic effects on the liver. Thus, there are numerous agonists for the varied binding sites of the receptor but no drug available yet. There does seem to be a ray of hope in the drugs that target FFAR1 but a lot more efforts towards drug discovery would result in the successful management of type 2 diabetes mellitus.

Elucidation of the Physiological Functions of Membrane Proteins as Novel Drug Target Candidate Molecules

For pharmaceutical research focused on identifying novel drug target candidate molecules, it is essential to explore unknown biological phenomena, elucidate underlying molecular mechanisms and regulate biological processes based on these findings. Proteins expressed on the plasma membrane and endoplasmic reticulum (ER) membrane play important roles in linking extracellular environmental information to intracellular processes. Stimulating membranous proteins induces various kinds of changes in cells, such as alterations in gene expression levels and enzymatic activities. However, the physiological functions and endogenous ligands of many G-protein-coupled receptors (GPCRs) have not been determined, although GPCRs already constitute a large class of drug-target membrane proteins. Furthermore, the precise physiological roles played by many ER membrane proteins have not been elucidated to date. In this review article, I summarize the results of our recent studies, including the observations that the lipid sensor FFAR4/GPR120 controlled systemic energy homeostasis and that the ER membrane monovalent cation channel trimeric intracellular cation (TRIC)-B and the plasma membrane divalent cation channel transient receptor potential melastatin 7 (TRPM7) regulated bone formation. I further describe the therapeutic significance of these membranous protein-related biological processes.

Molecular basis for kinin selectivity and activation of the human bradykinin receptors

Bradykinin and kallidin are endogenous kinin peptide hormones that belong to the kallikrein-kinin system and are essential to the regulation of blood pressure, inflammation, coagulation and pain control. Des-Arg¹⁰-kallidin, the carboxy-terminal des-Arg metabolite of kallidin, and bradykinin selectively activate two G protein-coupled receptors, type 1 and type 2 bradykinin receptors (B1R and B2R), respectively. The hyperactivation of bradykinin receptors, termed 'bradykinin storm', is associated with pulmonary edema in COVID-19 patients, suggesting that bradykinin receptors are important targets for COVID-19 intervention. Here we report two G protein-coupled complex structures of human B1R and B2R bound to des-Arg¹⁰-kallidin and bradykinin, respectively. Combined with functional analysis, our structures reveal the mechanism of ligand selectivity and specific activation of the bradykinin receptor. These findings also provide a framework for guiding drug design targeting bradykinin receptors for the treatment of inflammation, cardiovascular disorders and COVID-19.

Ligands of Adrenergic Receptors: A Structural Point of View

Adrenergic receptors are G protein-coupled receptors for epinephrine and norepinephrine. They are targets of many drugs for various conditions, including treatment of hypertension, hypotension, and asthma. Adrenergic receptors are intensively studied in structural biology, displayed for binding poses of different types of ligands. Here, we summarized molecular mechanisms of ligand recognition and receptor activation exhibited by structure. We also reviewed recent advances in structure-based ligand discovery against adrenergic receptors.

Environment-Sensitive Fluorescence of 7-Nitrobenz-2-oxa-1,3-diazol-4-yl (NBD)-Labeled Ligands for Serotonin Receptors

Serotonin is a neurotransmitter that plays a crucial role in the regulation of several behavioral and cognitive functions by binding to a number of different serotonin receptors present on the cell surface. We report here the synthesis and characterization of several novel fluorescent analogs of serotonin in which the fluorescent NBD (7-nitrobenz-2-oxa-1,3-diazol-4-yl) group is covalently attached to serotonin. The fluorescent ligands compete with the serotonin_1A receptor specific radiolabeled agonist for binding to the receptor. Interestingly, these fluorescent ligands display a high environmental sensitivity of their fluorescence. Importantly, the human serotonin_1A receptor stably expressed in CHO-K1 cells could be specifically labeled with one of the fluorescent ligands with minimal nonspecific labeling. Interestingly, we show by spectral imaging that the NBD-labeled ligand exhibits a red edge excitation shift (REES) of 29 nm when bound to the receptor, implying that it is localized in a restricted microenvironment. Taken together, our results show that NBD-labeled serotonin analogs offer an attractive fluorescent approach for elucidating the molecular environment of the serotonin binding site in serotonin receptors. In view of the multiple roles played by the serotonergic systems in the central and peripheral nervous systems, these fluorescent ligands would be useful in future studies involving serotonin receptors.

Design, Synthesis, and Structure-Activity Relationship Studies of (4-Alkoxyphenyl)glycinamides and Bioisosteric 1,3,4-Oxadiazoles as GPR88 Agonists

Increasing evidence implicates the orphan G protein-coupled receptor 88 (GPR88) in a number of striatal-associated disorders. In this study, we report the design and synthesis of a series of novel (4-alkoxyphenyl)glycinamides (e.g., 31) and the corresponding 1,3,4-oxadiazole bioisosteres derived from the 2-AMPP scaffold (1) as GPR88 agonists. The 5-amino-1,3,4-oxadiazole derivatives (84, 88-90) had significantly improved potency and lower lipophilicity compared to 2-AMPP. Compound 84 had an EC₅₀ of 59 nM in the GPR88 overexpressing cell-based cAMP assay. In addition, 84 had an EC₅₀ of 942 nM in the [³⁵S]GTPγS binding assay using mouse striatal membranes but was inactive in membranes from GPR88 knockout mice, even at a concentration of 100 μM. In vivo pharmacokinetic testing of 90 in rats revealed that the 5-amino-1,3,4-oxadiazole analogues may have limited brain permeability. Taken together, these results provide the basis for further optimization to develop a suitable agonist to probe GPR88 functions in the brain.

Molecular evolution of a collage of cholesterol interaction motifs in transmembrane helix V of the serotonin1A receptor

The human serotonin_1A receptor is a representative member of the superfamily of G protein-coupled receptors (GPCRs) and an important drug target for neurological disorders. Using a combination of biochemical, biophysical and molecular dynamics simulation approaches, we and others have shown that membrane cholesterol modulates the organization, dynamics and function of vertebrate serotonin_1A receptors. Previous studies have shown that the cytoplasmic portion of transmembrane helix V (TM V) and the extramembraneous intracellular loop 3 are critical for G-protein coupling, phosphorylation and desensitization of the receptor. We have recently resolved a collage of putative cholesterol interaction motifs from the amino acid sequence overlapping this region. In this paper, we explore the sequence plasticity of this fragment that may have adapted to altered membrane lipidome, after vertebrates evolved from primordial invertebrates. Since invertebrates have lower levels of membrane cholesterol relative to vertebrates, we compared TM V sequence fragments from invertebrate serotonin₁ receptors with vertebrate orthologs to infer the sequence plasticity in TM V. We report that the average number of cholesterol interaction motifs in TM V for diverse phyla represents an increasing trend that could mirror vertebrate evolution from primordial invertebrates. By statistical modeling, we propose that the collage of cholesterol interaction motifs in TM V of the human serotonin_1A receptor may have evolved from rudimentary collages, reminiscent of primordial invertebrate orthologs. Taken together, we propose that a repertoire of cholesterol-philic nonsynonymous substitutions may have enhanced collage complexity in TM V during vertebrate evolution.

Structural Insight into the Mechanism of 4-Aminoquinolines Selectivity for the alpha2A-Adrenoceptor

Background

α_2A-adrenoceptor (AR) is a potential target for the treatment of degenerative diseases of the central nervous system, and α_2A-AR agonists are effective drugs for this condition. However, the lack of high selectivity for α_2A-AR subtype of traditional drugs greatly limits their clinic usage.

Methods

A series of homobivalent 4-aminoquinolines conjugated by two 4-aminoquinoline moieties via varying alkane linker length (C2-C12) were characterized for their affinities for each α₂-AR subtype. Subsequently, docking, molecular dynamics and mutagenesis were applied to uncover the molecular mechanism.

Results

Most 4-aminoquinolines (4-aminoquinoline monomer, C2-C6, C8-C10) were selective for α_2A-AR over α_2B- and α_2C-ARs. Besides, the affinities are of similar linker length-dependence for each α₂-AR subtype. Among all the compounds tested, C10 has the highest affinity for α_2A-AR (pKi=−7.45±0.62), which is 12-fold and 60-fold selective over α_2B-AR and α_2C-AR, respectively. Docking and molecular dynamics suggest that C10 simultaneously interacts with orthosteric and “allosteric” sites of the α_2A-AR. The mutation of F205 decreases the affinity by 2-fold. The potential allosteric residues include S90, N93, E94 and W99.

Conclusion

The specificity of C10 for the α_2A-AR and the potential orthosteric and allosteric binding sites proposed in this study provide valuable guidance for the development of novel α_2A-AR subtype selective compounds.

Revealing the Positive Binding Cooperativity Mechanism between the Orthosteric and the Allosteric Antagonists of CCR2 by Metadynamics and Gaussian Accelerated Molecular Dynamics Simulations

CC chemokine receptor 2 (CCR2) and its endogenous CC chemokine ligands are associated with numerous inflammatory, neurodegenerative diseases, and cancer. CCR2 is becoming an attractive target in the treatment of autoimmune disease and neurodegenerative diseases. The orthosteric antagonist BMS-681 and allosteric antagonist CCR2-RA-[R] of CCR2 show positive binding cooperativity. We performed well-tempered metadynamics simulations and Gaussian accelerated MD simulations to reveal the influence of the orthosteric antagonist on the unbinding of allosteric antagonist of CCR2. We revealed different unbinding pathways of CCR2-RA-[R] in binary complex CCR2-VT5 and ternary complex CCR2-73R-VT5. The different unbinding pathways of CCR2-RA-[R] are due to the conformational dynamics of TM6. We obtained the significant conformational differences of the intracellular side of TM6 upon CCR2 binding to different ligands by GaMD simulation. The conformational dynamics of TM6 are consistent with the unbinding pathway analysis. GaMD simulations indicate that BMS-681 binding restricts the bend of intracellular side of TM6 by stabilizing the extracellular sides of TM6 and TM7. The charged residues Arg206^5.43 of TM5 and Glu291^7.39 of TM7 play key roles in stabling TM7 and TM6. TM6 and TM7 are crucial components in the orthosteric and allosteric binding sites. Our results illustrate the conformational details about the effect of the orthosteric antagonist on the allosteric antagonist of CCR2. The conformational dynamics of CCR2 upon binding to different ligands can provide a rational basis for development of allosteric ligands of CCR2.

The orphan receptor GPR88 blunts the signaling of opioid receptors and multiple striatal GPCRs

GPR88 is an orphan G protein-coupled receptor (GPCR) considered as a promising therapeutic target for neuropsychiatric disorders; its pharmacology, however, remains scarcely understood. Based on our previous report of increased delta opioid receptor activity in Gpr88 null mice, we investigated the impact of GPR88 co-expression on the signaling of opioid receptors in vitro and revealed that GPR88 inhibits the activation of both their G protein- and β-arrestin-dependent signaling pathways. In Gpr88 knockout mice, morphine-induced locomotor sensitization, withdrawal and supra-spinal analgesia were facilitated, consistent with a tonic inhibitory action of GPR88 on µOR signaling. We then explored GPR88 interactions with more striatal versus non-neuronal GPCRs, and revealed that GPR88 can decrease the G protein-dependent signaling of most receptors in close proximity, but impedes β-arrestin recruitment by all receptors tested. Our study unravels an unsuspected buffering role of GPR88 expression on GPCR signaling, with intriguing consequences for opioid and striatal functions.

Targeting Lysophosphatidic Acid in Cancer: The Issues in Moving from Bench to Bedside

Since the clear demonstration of lysophosphatidic acid (LPA)'s pathological roles in cancer in the mid-1990s, more than 1000 papers relating LPA to various types of cancer were published. Through these studies, LPA was established as a target for cancer. Although LPA-related inhibitors entered clinical trials for fibrosis, the concept of targeting LPA is yet to be moved to clinical cancer treatment. The major challenges that we are facing in moving LPA application from bench to bedside include the intrinsic and complicated metabolic, functional, and signaling properties of LPA, as well as technical issues, which are discussed in this review. Potential strategies and perspectives to improve the translational progress are suggested. Despite these challenges, we are optimistic that LPA blockage, particularly in combination with other agents, is on the horizon to be incorporated into clinical applications.

A collage of cholesterol interaction motifs in the serotonin1A receptor: An evolutionary implication for differential cholesterol interaction

The serotonin_1A receptor is a representative member of the G protein-coupled receptor (GPCR) superfamily and acts as an important drug target. In our previous work, we comprehensively demonstrated that membrane cholesterol is necessary in the organization, dynamics and function of the serotonin_1A receptor. In this context, analysis of high-resolution GPCR crystal structures in general and in silico studies of the serotonin_1A receptor in particular, have suggested the presence of cholesterol interaction sites (hotspots) in various regions of the receptor. In this work, we have identified an evolutionarily conserved collage of four categories of cholesterol interaction motifs associated with transmembrane helix V and the adjacent intracellular loop 3 fragment of the vertebrate serotonin_1A receptor. This collage of motifs represents a total of twenty diverse context-dependent cholesterol interaction configurations. We envision that the gamut of cholesterol interaction sites, characterized by sequence plasticity in cholesterol interaction, could be relevant in receptor-cholesterol interaction in membranes of varying cholesterol content and organization, as found in diverse cell types. We conclude that an evolutionarily conserved mechanism of GPCR-cholesterol interaction allows the serotonin_1A receptor to adapt to diverse membrane cholesterol levels during natural evolution.

Molecular features characterizing non-peptide selectivity to the human B1 and B2 bradykinin receptors

Bradykinin is produced in response to inflammation, trauma, burns, shock, allergy and some cardiovascular diseases. Actions of this peptide are mediated through two different G-protein coupled receptors, named B1 and B2 that have different pharmacological characteristics. The former is up-regulated during inflammation episodes or tissue trauma whereas, the latter is constitutively expressed in a variety of cell types. In a previous work we have characterized the molecular features that explain the observed structure-activity results for both receptors by means of molecular modeling studies, using diverse ligands for both receptors. These results were summarized in the form of two different pharmacophores that provided new insights to be used for the design of novel molecules with antagonistic profile. In the present work, we compare these pharmacophores to understand the features that characterize ligand selectivity to the two bradykinin receptors. The study shows that most of the residues involved in the binding pocket are similar in both receptors and consequently are the pharmacophores obtained. The main difference between the two pharmacophores remains on point #5 that involves a polar moiety for the B1 receptor and an aromatic ring for the B2 receptor. Accordingly, analysis of the prospective bound conformation of several non-selective small molecule ligands of the bradykinin receptors permits to conclude that fulfilment of point#5 is a requirement to produce selective ligands. However, the study also shows that this is a necessary condition only, since ligands need also to be bulky enough to avoid binding to these receptors in diverse poses. These results provide new insights for a better understanding of the molecular features that ligands are required to exhibit to be selective bradykinin ligands.

Salt Bridge in Ligand-Protein Complexes-Systematic Theoretical and Statistical Investigations

Although the salt bridge is the strongest among all known noncovalent molecular interactions, no comprehensive studies have been conducted to date to examine its role and significance in drug design. Thus, a systematic study of the salt bridge in biological systems is reported herein, with a broad analysis of publicly available data from Protein Data Bank, DrugBank, ChEMBL, and GPCRdb. The results revealed the distance and angular preferences as well as privileged molecular motifs of salt bridges in ligand-receptor complexes, which could be used to design the strongest interactions. Moreover, using quantum chemical calculations at the MP2 level, the energetic, directionality, and spatial variabilities of salt bridges were investigated using simple model systems mimicking salt bridges in a biological environment. Additionally, natural orbitals for chemical valence (NOCV) combined with the extended-transition-state (ETS) bond-energy decomposition method (ETS-NOCV) were analyzed and indicated a strong covalent contribution to the salt bridge interaction. The present results could be useful for implementation in rational drug design protocols.

GPCR-SAS: A web application for statistical analyses on G protein-coupled receptors sequences

G protein-coupled receptors (GPCRs) are one of the largest protein families in mammals. They mediate signal transduction across cell membranes and are important targets for the pharmaceutical industry. The G Protein-Coupled Receptors-Sequence Analysis and Statistics (GPCR-SAS) web application provides a set of tools to perform comparative analysis of sequence positions between receptors, based on a curated structural-informed multiple sequence alignment. The analysis tools include: (i) percentage of occurrence of an amino acid or motif and entropy at a position or range of positions, (ii) covariance of two positions, (iii) correlation between two amino acids in two positions (or two sequence motifs in two ranges of positions), and (iv) snake-plot representation for a specific receptor or for the consensus sequence of a group of selected receptors. The analysis of conservation of residues and motifs across transmembrane (TM) segments may guide the design of more selective ligands or help to rationalize activation mechanisms, among others. As an example, here we analyze the amino acids of the "transmission switch", that initiates receptor activation following ligand binding. The tool is freely accessible at http://lmc.uab.cat/gpcrsas/.

Chemical Diversity in the G Protein-Coupled Receptor Superfamily

G protein-coupled receptors (GPCRs) are the largest family of cell signaling transmembrane proteins that can be modulated by a plethora of chemical compounds. Systematic cheminformatics analysis of structurally and pharmacologically characterized GPCR ligands shows that cocrystallized GPCR ligands cover a significant part of chemical ligand space, despite their limited number. Many GPCR ligands and substructures interact with multiple receptors, providing a basis for polypharmacological ligand design. Experimentally determined GPCR structures represent a variety of binding sites and receptor-ligand interactions that can be translated to chemically similar ligands for which structural data are lacking. This integration of structural, pharmacological, and chemical information on GPCR-ligand interactions enables the extension of the structural GPCR-ligand interactome and the structure-based design of novel modulators of GPCR function.

The 3D Structure of Human DP Prostaglandin G-Protein-Coupled Receptor Bound to Cyclopentanoindole Antagonist, Predicted Using the DuplexBiHelix Modification of the GEnSeMBLE Method

Prostaglandins play a critical physiological role in both cardiovascular and immune systems, acting through their interactions with 9 prostanoid G protein-coupled receptors (GPCRs). These receptors are important therapeutic targets for a variety of diseases including arthritis, allergies, type 2 diabetes, and cancer. The DP prostaglandin receptor is of interest because it has unique structural and physiological properties. Most notably, DP does not have the 3-6 ionic lock common to Class A GPCRs. However, the lack of X-ray structures for any of the 9 prostaglandin GPCRs hampers the application of structure-based drug design methods to develop more selective and active medications to specific receptors. We predict here 3D structures for the DP prostaglandin GPCR, based on the GEnSeMBLE complete sampling with hierarchical scoring (CS-HS) methodology. This involves evaluating the energy of 13 trillion packings to finally select the best 20 that are stable enough to be relevant for binding to antagonists, agonists, and modulators. To validate the predicted structures, we predict the binding site for the Merck cyclopentanoindole (CPI) selective antagonist docked to DP. We find that the CPI binds vertically in the 1-2-7 binding pocket, interacting favorably with residues R310^7.40 and K76^2.54 with additional interactions with S313^7.43, S316^7.46, S19^1.35, etc. This binding site differs significantly from that of antagonists to known Class A GPCRs where the ligand binds in the 3-4-5-6 region. We find that the predicted binding site leads to reasonable agreement with experimental Structure-Activity Relationship (SAR). We suggest additional mutation experiments including K76^2.54, E129^3.49, L123^3.43, M270^6.40, F274^6.44 to further validate the structure, function, and activation mechanism of receptors in the prostaglandin family. Our structures and binding sites are largely consistent and improve upon the predictions by Li et al. ( J. Am. Chem. Soc. 2007 , 129 ( 35 ), 10720 ) that used our earlier MembStruk prediction methodology.

Comparative molecular field analysis and molecular dynamics studies of the dopamine D2 receptor antagonists without a protonatable nitrogen atom

The dopaminergic hypothesis of schizophrenia is the main concept explaining the direct reasons of schizophrenia and the effectiveness of current antipsychotics. All antipsychotics present on the market are potent dopamine D₂ receptor antagonists or partial agonists. In this work we investigate a series of dopamine D₂ receptor antagonists which do not fulfill the criteria of the classical pharmacophore model as they do not possess a protonatable nitrogen atom necessary to interact with the conserved Asp(3.32). Such compounds are interesting, inter alia, due to possible better pharmacokinetic profile when compared to basic, ionizable molecules. By means of homology modeling, molecular docking and molecular dynamics we determined that the compounds investigated interact with Asp(3.32) via their amide nitrogen atom. It was found that the studied compounds stabilize the receptor inactive conformation through the effect on the ionic lock, which is typical for GPCR antagonists. We constructed a CoMFA model for the studied compounds with the following statistics: R² = 0.95, Q² = 0.63. The quality of the CoMFA model was confirmed by high value of R² of the test set, equal 0.96. The CoMFA model indicated two regions where bulky substituents are favored and two regions where bulky substituents are not beneficial. Two red contour regions near carbonyl groups were identified meaning that negative charge would be favored here. Furthermore, the S-oxide group is connected with blue contour region meaning that positive charge is favored in this position. These findings may be applied for further optimization of the studied compound series.

The molecular basis of subtype selectivity of human kinin G-protein-coupled receptors

G-protein-coupled receptors (GPCRs) are the most important signal transducers in higher eukaryotes. Despite considerable progress, the molecular basis of subtype-specific ligand selectivity, especially for peptide receptors, remains unknown. Here, by integrating DNP-enhanced solid-state NMR spectroscopy with advanced molecular modeling and docking, the mechanism of the subtype selectivity of human bradykinin receptors for their peptide agonists has been resolved. The conserved middle segments of the bound peptides show distinct conformations that result in different presentations of their N and C termini toward their receptors. Analysis of the peptide-receptor interfaces reveals that the charged N-terminal residues of the peptides are mainly selected through electrostatic interactions, whereas the C-terminal segments are recognized via both conformations and interactions. The detailed molecular picture obtained by this approach opens a new gateway for exploring the complex conformational and chemical space of peptides and peptide analogs for designing GPCR subtype-selective biochemical tools and drugs.

A Structural Framework for GPCR Chemogenomics: What's In a Residue Number?

The recent surge of crystal structures of G protein-coupled receptors (GPCRs), as well as comprehensive collections of sequence, structural, ligand bioactivity, and mutation data, has enabled the development of integrated chemogenomics workflows for this important target family. This chapter will focus on cross-family and cross-class studies of GPCRs that have pinpointed the need for, and the implementation of, a generic numbering scheme for referring to specific structural elements of GPCRs. Sequence- and structure-based numbering schemes for different receptor classes will be introduced and the remaining caveats will be discussed. The use of these numbering schemes has facilitated many chemogenomics studies such as consensus binding site definition, binding site comparison, ligand repurposing (e.g. for orphan receptors), sequence-based pharmacophore generation for homology modeling or virtual screening, and class-wide chemogenomics studies of GPCRs.

Characterization of Ligand Binding to GPCRs Through Computational Methods

The recent increase in available G protein-coupled receptor structures now contributes decisively to the structure-based ligand design. In this context, computational approaches in combination with medicinal chemistry and pharmacology are extremely helpful. Here, we provide an update on our structure-based computational protocols, used to answer key questions related to GPCR-ligand binding. All combined, these techniques can shed light on ligand binding modes, determine the molecular basis of conformational selection, for agonists and antagonists, as well as of subtype selectivity. To illustrate each of these questions, we will consider examples from existing projects on three families of class A (rhodopsin-like) GPCRs: one small-molecule (nucleotide-like) family, i.e., the adenosine receptors, and two peptide-binding receptors: neuropeptide-Y and angiotensin II receptors. The successful application of the same computational protocols to investigate this diverse group of receptor families gives an idea of the general applicability of our methodology in the characterization of GPCR-ligand binding.

Modeling and Deorphanization of Orphan GPCRs

Despite tremendous efforts, approximately 120 GPCRs remain orphan. Their physiological functions and their potential roles in diseases are poorly understood. Orphan GPCRs are extremely important because they may provide novel therapeutic targets for unmet medical needs. As a complement to experimental approaches, molecular modeling and virtual screening are efficient techniques to discover synthetic surrogate ligands which can help to elucidate the role of oGPCRs. Constitutively activated mutants and recently published active structures of GPCRs provide stimulating opportunities for building active molecular models for oGPCRs and identifying activators using virtual screening of compound libraries. We describe the molecular modeling and virtual screening process we have applied in the discovery of surrogate ligands, and provide examples for CCKA, a simulated oGPCR, and for two oGPCRs, GPR52 and GPR34.

Pharmacological relationships and ligand discovery of G protein-coupled receptors revealed by simultaneous ligand and receptor clustering

Conventional ligand and receptor similarity methods have been extensively used for exposing pharmacological relationships and drug lead discovery. They may in some cases neglect minor relationships useful for target hopping particularly against the remote family members. To complement the conventional methods for capturing these minor relationships, we developed a new method that uses a SLARC (Simultaneous Ligand And Receptor Clustering) 2D map to simultaneously characterize the ligand structural and receptor binding-site sequence relationships of a receptor family. The SLARC maps of the rhodopsin-like GPCR family comprehensively revealed scaffold hopping, target hopping, and multi-target relationships for the ligands of both homologous and remote family members. Their usefulness in new ligand discovery was validated by guiding the prospective discovery of novel indole piperazinylpyrimidine dual-targeting adenosine A2A receptor antagonist and dopamine D2 agonist compounds. The SLARC approach is useful for revealing pharmacological relationships and discovering new ligands at target family levels.

Structure and Function of Peptide-Binding G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors and are important human drug targets. Of the 826 human GPCRs, 118 of them recognize endogenous peptide or protein ligands, and 30 of the 118 are targeted by approved drug molecules, including the very high-profile class B glucagon-like peptide 1 receptor. In this review, we analyze the 21 experimentally determined three-dimensional structures of the known peptide-binding GPCRs in relation to the endogenous peptides and drug molecules that modulate their cell signaling processes. Our integrated analyses reveal that half of the marketed drugs and most of the drugs in clinical trials that interact with peptide GPCRs are small molecules with a wide range of binding modes distinct from those of large peptide ligands. As we continue to collect additional data on these receptors from orthogonal approaches, including nuclear magnetic resonance and electron microscopy, we are beginning to understand how these receptors interact with their ligands at the molecular level and how improving the pharmacology of GPCR signal transduction requires us to study these receptors using multiple biophysical techniques.

Mechanism Exploration of Arylpiperazine Derivatives Targeting the 5-HT2A Receptor by In Silico Methods

As a G-protein coupled receptor, the 5-hydroxytryptamine 2A (5-HT_2A) receptor is known for its critical role in the cognitive, behavioural and physiological functions, and thus is a primary molecular target to treat psychiatric diseases, including especially depression. With purpose to explore the structural traits affecting the inhibitory activity, currently a dataset of 109 arylpiperazine derivatives as promising 5-HT_2A antagonists was built, based on which the ligand-based three-dimensional quantitative structure-activity relationship (3D-QSAR) study by using both comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) approaches was carried out. The resultant optimal CoMSIA model displays proper validity and predictability with cross-validated correlation coefficient Q² = 0.587, non-cross-validated correlation coefficient R²_ncv = 0.900 and predicted correlation coefficient for the test set of compounds R²_pre = 0.897, respectively. Besides, molecular docking was also conducted to investigate the binding mode between these ligands and the active site of the 5-HT_2A receptor. Meanwhile, as a docking supplementary tool to study the antagonists' conformation in the binding cavity, molecular dynamics (MD) simulation was also performed, providing further elucidation about the changes in the ligand-receptor complex. Lastly, some new molecules were also newly-designed based on the above results that are potential arylpiperazine antagonists of 5-HT_2A receptor. We hope that the present models and derived information may be of help for facilitating the optimization and design of novel potent antagonists as antidepressant drugs as well as exploring the interaction mechanism of 5-HT_2A antagonists.

Development and validation of a high-throughput calcium mobilization assay for the orphan receptor GPR88

Background

GPR88 is an orphan G protein-coupled receptor highly expressed in the striatum and is implicated in basal ganglia-associated disorders. However, the receptor functions of GPR88 are still largely unknown due to the lack of potent and selective ligands appropriate for central nervous system investigation. Development of a high-throughput screening assay for GPR88 should facilitate the discovery of novel ligands to probe GPR88 functions.

Methods

In this paper, we describe the development of a CHO-Gα_qi5-GPR88 cell-based calcium mobilization assay. The assay takes advantage of functional coupling of GPR88 with the promiscuous Gα_qi5 protein and consequent mobilization of intracellular calcium, which can be measured in a 384-well format with a Fluorescent Imaging Plate Reader.

Results

The CHO-Gα_qi5-GPR88 cell-based calcium mobilization assay was validated by the structure-activity relationship study of known GPR88 agonist (1R,2R)-2-PCCA analogues. The assay was automated and miniaturized to a 384-well format, and was deemed robust and reproducible with a Z’-factor of 0.72 and tolerated dimethyl sulfoxide to a final concentration of 2%. Screening a pilot neurotransmitter library consisting of 228 compounds yielded 10 hits, but none of the hits were confirmed as GPR88 agonists in follow-up assays.

Conclusions

We have developed a high-throughput calcium mobilization assay for the orphan receptor GPR88. This calcium mobilization assay can be used to identify several different types of GPR88 ligands including agonists, competitive and noncompetitive antagonists, inverse agonists, and allosteric modulators. These ligands will serve as valuable tools to probe signaling mechanisms and in vivo functions of GPR88, and could expedite development of novel therapies for diseases potentially mediated by GPR88.

Electronic supplementary material

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

The impact of in silico screening in the discovery of novel and safer drug candidates

Drug discovery is a multidisciplinary and multivariate optimization endeavor. As such, in silico screening tools have gained considerable importance to archive, analyze and exploit the vast and ever-increasing amount of experimental data generated throughout the process. The current review will focus on the computer-aided prediction of the numerous properties that need to be controlled during the discovery of a preliminary hit and its promotion to a viable clinical candidate. It does not pretend to the almost impossible task of an exhaustive report but will highlight a few key points that need to be collectively addressed both by chemists and biologists to fuel the drug discovery pipeline with innovative and safe drug candidates.

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.

Orphan receptor ligand discovery by pickpocketing pharmacological neighbors

Understanding the pharmacological similarity of G protein-coupled receptors (GPCRs) is paramount for predicting ligand off-target effects, drug repurposing, and ligand discovery for orphan receptors. Phylogenetic relationships do not always correctly capture pharmacological similarity. Previous family-wide attempts to define pharmacological relationships were based on three-dimensional structures and/or known receptor-ligand pairings, both unavailable for orphan GPCRs. Here, we present GPCR-CoINPocket, a novel contact-informed neighboring pocket metric of GPCR binding-site similarity that is informed by patterns of ligand-residue interactions observed in crystallographically characterized GPCRs. GPCR-CoINPocket is applicable to receptors with unknown structure or ligands and accurately captures known pharmacological relationships between GPCRs, even those undetected by phylogeny. When applied to orphan receptor GPR37L1, GPCR-CoINPocket identified its pharmacological neighbors, and transfer of their pharmacology aided in discovery of the first surrogate ligands for this orphan with a 30% success rate. Although primarily designed for GPCRs, the method is easily transferable to other protein families.

Oligomerization of GPCRs involved in endocrine regulation

More than 800 different human membrane-spanning G-protein-coupled receptors (GPCRs) serve as signal transducers at biological barriers. These receptors are activated by a wide variety of ligands such as peptides, ions and hormones, and are able to activate a diverse set of intracellular signaling pathways. GPCRs are of central importance in endocrine regulation, which underpins the significance of comprehensively studying these receptors and interrelated systems. During the last decade, the capacity for multimerization of GPCRs was found to be a common and functionally relevant property. The interaction between GPCR monomers results in higher order complexes such as homomers (identical receptor subtype) or heteromers (different receptor subtypes), which may be present in a specific and dynamic monomer/oligomer equilibrium. It is widely accepted that the oligomerization of GPCRs is a mechanism for determining the fine-tuning and expansion of cellular processes by modification of ligand action, expression levels, and related signaling outcome. Accordingly, oligomerization provides exciting opportunities to optimize pharmacological treatment with respect to receptor target and tissue selectivity or for the development of diagnostic tools. On the other hand, GPCR heteromerization may be a potential reason for the undesired side effects of pharmacological interventions, faced with numerous and common mutual signaling modifications in heteromeric constellations. Finally, detailed deciphering of the physiological occurrence and relevance of specific GPCR/GPCR-ligand interactions poses a future challenge. This review will tackle the aspects of GPCR oligomerization with specific emphasis on family A GPCRs involved in endocrine regulation, whereby only a subset of these receptors will be discussed in detail.

Elimination of a ligand gating site generates a supersensitive olfactory receptor

Olfaction poses one of the most complex ligand-receptor matching problems in biology due to the unparalleled multitude of odor molecules facing a large number of cognate olfactory receptors. We have recently deorphanized an olfactory receptor, TAAR13c, as a specific receptor for the death-associated odor cadaverine. Here we have modeled the cadaverine/TAAR13c interaction, exchanged predicted binding residues by site-directed mutagenesis, and measured the activity of the mutant receptors. Unexpectedly we observed a binding site for cadaverine at the external surface of the receptor, in addition to an internal binding site, whose mutation resulted in complete loss of activity. In stark contrast, elimination of the external binding site generated supersensitive receptors. Modeling suggests this site to act as a gate, limiting access of the ligand to the internal binding site and thereby downregulating the affinity of the native receptor. This constitutes a novel mechanism to fine-tune physiological sensitivity to socially relevant odors.

Discovery of a novel small molecule agonist scaffold for the APJ receptor

The apelinergic system includes a series of endogenous peptides apelin, ELABELA/TODDLER and their 7-transmembrane G-protein coupled apelin receptor (APJ, AGTRL-1, APLNR). The APJ receptor is an attractive therapeutic target because of its involvement in cardiovascular diseases and potentially other disorders including liver fibrosis, obesity, diabetes, and neuroprotection. To date, pharmacological characterization of the APJ receptor has been limited due to the lack of small molecule functional agonists or antagonists. Through focused screening we identified a drug-like small molecule agonist hit 1 with a functional EC50 value of 21.5±5μM and binding affinity (Ki) of 5.2±0.5μM. Initial structure-activity studies afforded compound 22 having a 27-fold enhancement in potency and the first sub-micromolar full agonist with an EC50 value of 800±0.1nM and Ki of 1.3±0.3μM. Preliminary SAR, synthetic methodology, and in vitro pharmacological characterization indicate this scaffold will serve as a favorable starting point for further refinement of APJ potency and selectivity.

New paradigms in chemokine receptor signal transduction: Moving beyond the two-site model

Chemokine receptor (CKR) signaling forms the basis of essential immune cellular functions, and dysregulated CKR signaling underpins numerous disease processes of the immune system and beyond. CKRs, which belong to the seven transmembrane domain receptor (7TMR) superfamily, initiate signaling upon binding of endogenous, secreted chemokine ligands. Chemokine-CKR interactions are traditionally described by a two-step/two-site mechanism, in which the CKR N-terminus recognizes the chemokine globular core (i.e. site 1 interaction), followed by activation when the unstructured chemokine N-terminus is inserted into the receptor TM bundle (i.e. site 2 interaction). Several recent studies challenge the structural independence of sites 1 and 2 by demonstrating physical and allosteric links between these supposedly separate sites. Others contest the functional independence of these sites, identifying nuanced roles for site 1 and other interactions in CKR activation. These developments emerge within a rapidly changing landscape in which CKR signaling is influenced by receptor PTMs, chemokine and CKR dimerization, and endogenous non-chemokine ligands. Simultaneous advances in the structural and functional characterization of 7TMR biased signaling have altered how we understand promiscuous chemokine-CKR interactions. In this review, we explore new paradigms in CKR signal transduction by considering studies that depict a more intricate architecture governing the consequences of chemokine-CKR interactions.

Developmental and adult expression patterns of the G-protein-coupled receptor GPR88 in the rat: Establishment of a dual nuclear-cytoplasmic localization

GPR88 is a neuronal cerebral orphan G-protein-coupled receptor (GPCR) that has been linked to various psychiatric disorders. However, no extensive description of its localization has been provided so far. Here, we investigate the spatiotemporal expression of the GPR88 in prenatal and postnatal rat tissues by using in situ hybridization and immunohistochemistry. GPR88 protein was initially detected at embryonic day 16 (E16) in the striatal primordium. From E16-E20 to adulthood, the highest expression levels of both protein and mRNA were observed in striatum, olfactory tubercle, nucleus accumbens, amygdala, and neocortex, whereas in spinal cord, pons, and medulla GPR88 expression remains discrete. We observed an intracellular redistribution of GPR88 during cortical lamination. In the cortical plate of the developing cortex, GPR88 presents a classical GPCR plasma membrane/cytoplasmic localization that shifts, on the day of birth, to nuclei of neurons progressively settling in layers V to II. This intranuclear localization remains throughout adulthood and was also detected in monkey and human cortex as well as in the amygdala and hypothalamus of rats. Apart from the central nervous system, GPR88 was transiently expressed at high levels in peripheral tissues, including adrenal cortex (E16-E21) and cochlear ganglia (E19-P3), and also at moderate levels in retina (E18-E19) and spleen (E21-P7). The description of the GPR88 anatomical expression pattern may provide precious functional insights into this novel receptor. Furthermore, the GRP88 nuclear localization suggests nonclassical GPCR modes of action of the protein that could be relevant for cortical development and psychiatric disorders. J. Comp. Neurol. 524:2776-2802, 2016. © 2016 Wiley Periodicals, Inc.

ECL1i, d(LGTFLKC), a novel, small peptide that specifically inhibits CCL2-dependent migration

CC chemokine receptor type 2 (CCR2) is a key molecule in inflammatory diseases and is an obvious drug target for the treatment of inflammation. A number of nonpeptidic, competitive CCR2 antagonists have been developed, but none has yet been approved for clinical use. Our aim was to identify a short peptide that showed allosteric antagonism against human and mouse CCR2. On the basis of sequence analysis and 3-dimensional modeling, we identified an original 7-d-amino acid peptidic CCR2 inhibitor that we have called extracellular loop 1 inverso (ECL1i), d(LGTFLKC). In vitro, ECL1i selectively and potently inhibits CC chemokine ligand type 2 (CCL2)-triggered chemotaxis (IC50, 2 µM) but no other conventional CCL2-associated events. We used the classic competitive CCR2 antagonist, BMS22 {2-[(isopropylaminocarbonyl)amino]-N-[2-[[cis-2-[[4-(methylthio)benzoyl]amino]cyclohexyl]amino]-2-oxoethyl]-5-(trifluoromethyl)benzamide}, as positive control and inhibited CCL2-dependent chemotaxis with an IC50 of 18 nM. As negative control, we used a peptide with the same composition as ECL1i, but in a different sequence, d(FKLTLCG). In vivo, ECL1i (4 mg/kg) interfered with CCR2-positive cell recruitment and attenuated disease progression in experimental autoimmune encephalomyelitis, a mouse model of multiple sclerosis. This study establishes ECL1i as the first allosteric inhibitor of CCR2 with functional selectivity. ECL1i is a promising new agent in therapeutic development, and it may, by its selective effect, increase our understanding of CCR2 signaling pathways and functions.-Auvynet, C., Baudesson de Chanville, C., Hermand, P., Dorgham, K., Piesse, C., Pouchy, C., Carlier, L., Poupel, L., Barthélémy, S., Felouzis, V., Lacombe, C., Sagan, S., Salomon, B., Deterre, P., Sennlaub, F., Combadière, C. ECL1i, d(LGTFLKC), a novel, small peptide that specifically inhibits CCL2-dependent migration.

Non-classical amine recognition evolved in a large clade of olfactory receptors

Biogenic amines are important signaling molecules, and the structural basis for their recognition by G Protein-Coupled Receptors (GPCRs) is well understood. Amines are also potent odors, with some activating olfactory trace amine-associated receptors (TAARs). Here, we report that teleost TAARs evolved a new way to recognize amines in a non-classical orientation. Chemical screens de-orphaned eleven zebrafish TAARs, with agonists including serotonin, histamine, tryptamine, 2-phenylethylamine, putrescine, and agmatine. Receptors from different clades contact ligands through aspartates on transmembrane α-helices III (canonical Asp^3.32) or V (non-canonical Asp^5.42), and diamine receptors contain both aspartates. Non-classical monoamine recognition evolved in two steps: an ancestral TAAR acquired Asp^5.42, gaining diamine sensitivity, and subsequently lost Asp^3.32. Through this transformation, the fish olfactory system dramatically expanded its capacity to detect amines, ecologically significant aquatic odors. The evolution of a second, alternative solution for amine detection by olfactory receptors highlights the tremendous structural versatility intrinsic to GPCRs.

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

Many organisms make molecules called biogenic amines. These molecules, which include the human hormones adrenaline and histamine, have important roles in regulating the biology and behaviour of many animals. Some biogenic amines bind to receptor proteins called GPCRs on the surface of cells. Many drugs can affect the activity of GPCRs, so understanding how different GPCRs work is an important goal of the pharmaceutical industry. Like all proteins, GPCRs are made of chains of molecules called amino acids. The GPCRs that can detect biogenic amines use a particular amino acid named Asp^3.32, and when this amino acid is mutated, these GPCRs become unable to bind to their target amine.

Trace amine-associated receptors (TAARs) are a type of GPCR that are found in many animals to detect odors. Most TAARs in mammals contain the Asp^3.32 residue, and recognize amine odors. However, many fish TAARs do not contain Asp^3.32, and it was not clear what molecules these fish receptors detect.

Here Li et al. find that these fish TAARs also recognize amines, and use a different amino acid called Asp^5.42. Also, some TAARs contain both Asp^3.32 and Asp^5.42, and recognize chemicals with two amines named diamines. Some diamines that bind to TAARs are foul smelling odors; for example, cadaverine and putrescine are repulsive smells emitted by decomposing flesh. In total, the experiments identified amines that can bind to eleven zebrafish TAARs that previously had no odor partner.

Li et al. propose that some fish TAARs lost the Asp^3.32 during the course of evolution to leave the Asp^5.42 as the main interaction site for amines. This change dramatically altered how these TAARs interact with amines, which probably expanded the number of different amines that fish can detect. These findings open up new ways to study how the fish brain processes information about its surroundings.

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

Covalent molecular probes for class A G protein-coupled receptors: advances and applications

Covalent modification of G protein-coupled receptors (GPCRs) by employing specific molecular probes has for decades provided a successful strategy to facilitate the elucidation of the structure and function of this pharmacologically important class of membrane proteins. The ligands typically comprise a pharmacophore that generates affinity for a given GPCR and contain a reactive functionality that may form a covalent bond with a suitably positioned amino acid residue. Covalent ligands have been successfully applied to circumvent poor affinity of compounds when stable labeling of receptor populations was required, and they have been used in the isolation, purification, and pharmacological characterization of specific subtypes of GPCRs. Recently, structural studies have demonstrated that covalent molecular probes are effective at stabilizing GPCRs to obtain X-ray crystal structures, thus providing valuable insights for the development of novel therapeutics. Herein, we review covalently binding molecular probes for class A GPCRs with a focus on ligands comprising cross-linking groups that do not require photoactivation and further highlight their significant and diverse applications.

Studying the binding interactions of allosteric agonists and antagonists of the CXCR4 receptor

Several examples of allosteric modulators of GPCRs have been reported recently in the literature, but understanding their molecular mechanism presents a new challenge for medicinal chemistry. For the specific case of the cellular receptor CXCR4, it is known that pepducins (lipidated fragments of intracellular GPCR loops) such as ATI-2341 modulate CXCR4 activity agonistically via an allosteric mechanism. Moreover, there are also examples of small organic molecules such as AMD11070 and GSK812397 which may also act as allosteric antagonists. However, incomplete knowledge of the ligand-binding sites has hampered a detailed molecular understanding of how these inhibitors work. Here, we attempt to answer this question by analysing the binding interactions between the CXCR4 receptor and the above-mentioned allosteric modulators. We propose two different allosteric binding sites, one located in the intracellular loops 1, 2 and 3 (ICL1, ICL2 and ICL3) which binds the pepducin agonist ATI-2341, and the other at a subsite of the main extracellular orthosteric binding pocket between extracellular loops 1 and 2 and the N-terminus, which binds the antagonists AMD11070 and GSK812397. Allosteric interactions between the CXCR4 and ATI-2341 were predicted by combining different modeling approaches. First, a rotational blind docking search was applied and the best poses were subsequently refined using flexible docking methods and molecular dynamic simulations. For the AMD11070 and GSK812397 antagonists, the entire CXCR4 protein surface was explored by blind docking in order to define the binding region. A second docking analysis by subsites was then performed to refine the allosteric interactions. Finally, we identified the binding residues that appear to be essential for CXCR4 allosteric modulators.