Unifying family A GPCR theories of activation

Exploring the Activation Mechanism of a Metabotropic Glutamate Receptor Homodimer via Molecular Dynamics Simulation

Metabotropic glutamate receptors of class C GPCRs exist as constitutive dimers, which play important roles in activating excitatory synapses of the central nervous system. However, the activation mechanism induced by agonists has not been clarified in experiments. To address the problem, we used microsecond all-atom molecular dynamics (MD) simulation couple with protein structure network (PSN) to explore the glutamate-induced activation for the mGluR1 homodimer. The results indicate that glutamate binding stabilizes not only the closure of Venus flytrap domains but also the polar interaction of LB2-LB2, in turn keeping the extracelluar domain in the active state. The activation of the extracelluar domain drives transmembrane domains (TMDs) of the two protomers closer and induces asymmetric activation for the TMD domains of the two protomers. One protomer with lower binding affinity to the agonist is activated, while the other protomer with higher binding energy is still in the inactive state. The PSN analysis identifies the allosteric regulation pathway from the ligand-binding pocket in the extracellular domain to the G-protein binding site in the intracellular TMD region and further reveals that the asymmetric activation is attributed to a combination of trans-pathway and cis-pathway regulations from two glumatates, rather than a single activation pathway. These observations could provide valuable molecular information for understanding of the structure and the implications in drug efficacy for the class C GPCR dimers.

Preferential Coupling of Dopamine D2S and D2L Receptor Isoforms with Gi1 and Gi2 Proteins-In Silico Study

The dopamine D₂ receptor belongs to rhodopsin-like G protein-coupled receptors (GPCRs) and it is an important molecular target for the treatment of many disorders, including schizophrenia and Parkinson's disease. Here, computational methods were used to construct the full models of the dopamine D₂ receptor short (D_2S) and long (D_2L) isoforms (differing with 29 amino acids insertion in the third intracellular loop, ICL3) and to study their coupling with G_i1 and G_i2 proteins. It was found that the D_2L isoform preferentially couples with the G_i2 protein and D_2S isoform with the G_i1 protein, which is in accordance with experimental data. Our findings give mechanistic insight into the interplay between isoforms of dopamine D₂ receptors and G_i proteins subtypes, which is important to understand signaling by these receptors and their mediation by pharmaceuticals, in particular psychotic and antipsychotic agents.

Advances in the Treatment of Explicit Water Molecules in Docking and Binding Free Energy Calculations

Background:

The inclusion of direct effects mediated by water during the ligandreceptor recognition is a hot-topic of modern computational chemistry applied to drug discovery and development. Docking or virtual screening with explicit hydration is still debatable, despite the successful cases that have been presented in the last years. Indeed, how to select the water molecules that will be included in the docking process or how the included waters should be treated remain open questions.

Objective:

In this review, we will discuss some of the most recent methods that can be used in computational drug discovery and drug development when the effect of a single water, or of a small network of interacting waters, needs to be explicitly considered.

Results:

Here, we analyse the software to aid the selection, or to predict the position, of water molecules that are going to be explicitly considered in later docking studies. We also present software and protocols able to efficiently treat flexible water molecules during docking, including examples of applications. Finally, we discuss methods based on molecular dynamics simulations that can be used to integrate docking studies or to reliably and efficiently compute binding energies of ligands in presence of interfacial or bridging water molecules.

Conclusions:

Software applications aiding the design of new drugs that exploit water molecules, either as displaceable residues or as bridges to the receptor, are constantly being developed. Although further validation is needed, workflows that explicitly consider water will probably become a standard for computational drug discovery soon.

Common activation mechanism of class A GPCRs

Class A G-protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology. Understanding receptor activation mechanism is critical for discovering novel therapeutics since about one-third of all marketed drugs target members of this family. GPCR activation is an allosteric process that couples agonist binding to G-protein recruitment, with the hallmark outward movement of transmembrane helix 6 (TM6). However, what leads to TM6 movement and the key residue level changes of this movement remain less well understood. Here, we report a framework to quantify conformational changes. By analyzing the conformational changes in 234 structures from 45 class A GPCRs, we discovered a common GPCR activation pathway comprising of 34 residue pairs and 35 residues. The pathway unifies previous findings into a common activation mechanism and strings together the scattered key motifs such as CWxP, DRY, Na+ pocket, NPxxY and PIF, thereby directly linking the bottom of ligand-binding pocket with G-protein coupling region. Site-directed mutagenesis experiments support this proposition and reveal that rational mutations of residues in this pathway can be used to obtain receptors that are constitutively active or inactive. The common activation pathway provides the mechanistic interpretation of constitutively activating, inactivating and disease mutations. As a module responsible for activation, the common pathway allows for decoupling of the evolution of the ligand binding site and G-protein-binding region. Such an architecture might have facilitated GPCRs to emerge as a highly successful family of proteins for signal transduction in nature.

The effect of dimerization on the activation and conformational dynamics of adenosine A1 receptor

The adenosine A1 receptor (A1R) is one of four adenosine receptors in humans, which are involved in the function of the cardiovascular, respiratory and central nervous systems. Experimental results indicate that A1R can form a homodimer and that the protomer-protomer interaction in the A1R dimer is related to certain pharmacological characteristics of A1R activation. In this work, we performed docking, metadynamics simulation, conventional molecular dynamics simulations, Gaussian-accelerated molecular dynamics simulations, potential of mean force calculations, dynamic cross-correlation motions analysis and community network analysis to study the binding mode of 5'-N-ethylcarboxamidoadenosine (NECA) to A1R and the effect of dimerization on the activation of A1R. Our results show that NECA binds to A1R in a similar mode to adenosine in the A1R crystal structure and NECA in the A2AR crystal structure. The A1R homodimer can be activated by one or two agonists with NECA occupying its orthosteric pockets in one (which we call the NECA-A1R system) or both protomers (which we call the dNECA-A1R system). In the NECA-A1R system, activation is predicated in the protomer without NECA bound. In the dNECA-A1R system, only one protomer achieves the active state. These findings suggest an asymmetrical activation mechanism of the homodimer and a negative cooperativity between the two protomers. We envision that our results may further facilitate the drug development of A1R.

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.

Progress in the development of more effective and safer analgesics for pain management

Opioid analgesics have been used for thousands of years in the treatment of pain and related disorders, and have become among the most widely prescribed medications. Among opioid analgesics, mu opioid receptor (MOR) agonists are the most commonly used and are indicated for acute and chronic pain management. However, their use results in a plethora of well-described side-effects. From selective delta opioid receptor (DOR) and kappa opioid receptor (KOR) agonists to multitarget MOR/DOR and MOR/KOR ligands, medicinal chemistry provided different approaches aimed at the development of opioid analgesics with an improved pharmacological and tolerability fingerprint. The emergent medicinal chemistry strategy to develop ameliorated opioid analgesics is based upon the concept that functional selectivity for G-protein signalling is necessary for the therapeutic effect, whether β-arrestin recruitment is mainly responsible for the manifestation of side effects, including the development of tolerance after repeated administrations. This review summarises most relevant biased MOR, DOR, KOR and multitarget MOR/DOR ligands synthesised in the last decade and their pharmacological profile in "in vitro" and "in vivo" studies. Such biased ligands could have a significant impact on modern drug discovery and represent a new strategy for the development of better-tolerated drug candidates.

G protein-coupled receptors of class A harness the energy of membrane potential to increase their sensitivity and selectivity

The human genome contains about 700 genes of G protein-coupled receptors (GPCRs) of class A; these seven-helical membrane proteins are the targets of almost half of all known drugs. In the middle of the helix bundle, crystal structures reveal a highly conserved sodium-binding site, which is connected with the extracellular side by a water-filled tunnel. This binding site contains a sodium ion in those GPCRs that are crystallized in their inactive conformations but does not in those GPCRs that are trapped in agonist-bound active conformations. The escape route of the sodium ion upon the inactive-to-active transition and its very direction have until now remained obscure. Here, by modeling the available experimental data, we show that the sodium gradient over the cell membrane increases the sensitivity of GPCRs if their activation is thermodynamically coupled to the sodium ion translocation into the cytoplasm but decreases it if the sodium ion retreats into the extracellular space upon receptor activation. The model quantitatively describes the available data on both activation and suppression of distinct GPCRs by membrane voltage. The model also predicts selective amplification of the signal from (endogenous) agonists if only they, but not their (partial) analogs, induce sodium translocation. Comparative structure and sequence analyses of sodium-binding GPCRs indicate a key role for the conserved leucine residue in the second transmembrane helix (Leu2.46) in coupling sodium translocation to receptor activation. Hence, class A GPCRs appear to harness the energy of the transmembrane sodium potential to increase their sensitivity and selectivity.

Characterization, Dynamics, and Mechanism of CXCR4 Antagonists on a Constitutively Active Mutant

The G protein-coupled receptor (GPCR) CXCR4 is a co-receptor for HIV and is involved in cancers and autoimmune diseases. We characterized five purine or quinazoline core polyamine pharmacophores used for targeting CXCR4 dysregulation in diseases. All were neutral antagonists for wild-type CXCR4 and two were biased antagonists with effects on β-arrestin-2 only at high concentrations. These compounds displayed various activities for a constitutively active mutant (CAM). We use the IT1t-CXCR4 crystal structure and molecular dynamics (MD) simulations to develop two hypotheses for the activation of the N1193.35A CAM. The N1193.35A mutation facilitates increased coupling of TM helices III and VI. IT1t deactivates the CAM by disrupting the coupling between TM helices III and VI, mediated primarily by residue F872.53. Mutants of F872.53 in N1193.35A CXCR4 precluded constitutive signaling and prevented inverse agonism. This work characterizes CXCR4 ligands and provides a mechanism for N1193.35A constitutive activation.

Allosteric Modulation Mechanism of the mGluR5 Transmembrane Domain

Positive allosteric modulators (PAMs) of metabotropic glutamate receptor type 5 (mGluR₅), a prototypical class C G protein-coupled receptor (GPCR), have shown therapeutic potential for various neurological disorders. Understanding the allosteric activation mechanism is essential for the rational design of mGluR₅ PAMs. We studied the actions of positive and negative allosteric modulators within the transmembrane domain of mGluR₅, using enhance-sampling all-atom molecular dynamics simulations. We found dual binding modes of the PAM, associated with distinct shapes of the allosteric pocket. The negative allosteric modulators, in contrast, showed only one binding mode. The simulations revealed the mechanism by which the PAM activated the receptor, in the absence of the orthosteric agonist (the so-called allosteric agonism). The mechanism relied on dynamic communications between amino-acid motifs that are highly conserved across class C GPCRs. The findings may guide structure-based design and virtual screening of allosteric modulators for mGluR₅ as well as for other class C GPCRs.

Understanding Ligand Binding to G-Protein Coupled Receptors Using Multiscale Simulations

Human G-protein coupled receptors (GPCRs) convey a wide variety of extracellular signals inside the cell and they are one of the main targets for pharmaceutical intervention. Rational drug design requires structural information on these receptors; however, the number of experimental structures is scarce. This gap can be filled by computational models, based on homology modeling and docking techniques. Nonetheless, the low sequence identity across GPCRs and the chemical diversity of their ligands may limit the quality of these models and hence refinement using molecular dynamics simulations is recommended. This is the case for olfactory and bitter taste receptors, which constitute the first and third largest GPCR groups and show sequence identities with the available GPCR templates below 20%. We have developed a molecular dynamics approach, based on the combination of molecular mechanics and coarse grained (MM/CG), tailored to study ligand binding in GPCRs. This approach has been applied so far to bitter taste receptor complexes, showing significant predictive power. The protein/ligand interactions observed in the simulations were consistent with extensive mutagenesis and functional data. Moreover, the simulations predicted several binding residues not previously tested, which were subsequently verified by carrying out additional experiments. Comparison of the simulations of two bitter taste receptors with different ligand selectivity also provided some insights into the binding determinants of bitter taste receptors. Although the MM/CG approach has been applied so far to a limited number of GPCR/ligand complexes, the excellent agreement of the computational models with the mutagenesis and functional data supports the applicability of this method to other GPCRs for which experimental structures are missing. This is particularly important for the challenging case of GPCRs with low sequence identity with available templates, for which molecular docking shows limited predictive power.

Nociceptin Receptor Is Overexpressed in Non-small Cell Lung Cancer and Predicts Poor Prognosis

Classic opioid receptors, mu (μ), delta (δ), and kappa (κ), have been reported to be expressed in non-small cell lung cancer (NSCLC) cell lines and tumor tissues and to play a role in tumor prognosis. However, the expression and role of the non-classic opioid receptor, nociceptin receptor (NOP) in cancer are unclear. Our hypothesis was that NOP was also highly expressed in NSCLC tumor tissues and this could be correlated with patients' prognostic characters. Expression of NOP was examined in archived cancer tissues from 129 enrolled NSCLC patients by immunohistochemistry and was further analyzed with the patients' outcomes. NOP expression in NSCLC cell lines was also detected. The dataset from Kaplan-Meier Plotter was used to explore the correlation between the levels of NOP mRNA in cancerous tissue and the prognosis of NSCLC patients. Cell functional assays were performed to detect the effect of NOP activation on tumor aggressive furthers. Results showed NOP expression was highly expressed in cancer tissues and human cancer cell lines. NOP expression was not associated with patients' opioid requirement but closely with some clinicopathological indicators which reflected the malignancy. Moreover, NOP staining level was the independent poor prognostic factor for NSCLC patients receiving lobectomy, which was further verified by determining the mRNA expression levels through the online dataset. In vitro experiments revealed that NOP activation promotes the proliferation and invasion of A549 cells via PI3K/Akt signaling pathway. We conclude that NOP is overexpressed in NSCLC and is inversely correlated with patient's postoperative survival.

Structural Insight into G Protein-Coupled Receptor Signaling Efficacy and Bias between Gs and β-Arrestin

G protein-coupled receptors (GPCRs) form the largest family of membrane proteins involved in signal transduction. Because of their ability to regulate a wide range of cellular responses and their dysregulation being associated with many diseases, GPCRs remain a key therapeutic target for several clinical indications. In recent years, it has been demonstrated that ligands for a given receptor can engage distinct pathways with different relative efficacies, a concept known as biased signaling or functional selectivity. However, the structural determinants of this phenomenon remain poorly understood. Using the β2-adrenergic receptor as a model, we identified a linker residue (L124^3.43) between the known PIF and NPxxY structural motifs, that plays a central role in the differential efficacy of biased ligands toward the Gs and β-arrestin pathways. Given the high level of conservation of this linker residue, the study provides structural explanations for biased signaling that can be extrapolated to other GPCRs.

Dynamic Role of the G Protein in Stabilizing the Active State of the Adenosine A2A Receptor

Agonist binding in the extracellular region of the G protein-coupled adenosine A2A receptor increases its affinity to the G proteins in the intracellular region, and vice versa. The structural basis for this effect is not evident from the crystal structures of A_2AR in various conformational states since it stems from the receptor dynamics. Using atomistic molecular dynamics simulations on four different conformational states of the adenosine A_2A receptor, we observed that the agonists show decreased ligand mobility, lower entropy of the extracellular loops in the active-intermediate state compared with the inactive state. In contrast, the entropy of the intracellular region increases to prime the receptor for coupling the G protein. Coupling of the G protein to A_2AR shrinks the agonist binding site, making tighter receptor agonist contacts with an increase in the strength of allosteric communication compared with the active-intermediate state. These insights provide a strong basis for structure-based ligand design studies.

•

GPCR conformation dynamics reveals the forward and backward allosteric mechanism

•

Agonist binding increases the entropy in the intracellular region of the GPCR

•

G protein binding shrinks the receptor-ligand contacts in the extracellular region

•

Increased allostery between G protein and agonist in the GPCR-G protein complex

Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

The aminergic family of G protein-coupled receptors (GPCRs) plays an important role in various diseases and represents a major drug discovery target class. Structure determination of all major aminergic subfamilies has enabled structure-based ligand design for these receptors. Site-directed mutagenesis data provides an invaluable complementary source of information for elucidating the structural determinants of binding of different ligand chemotypes. The current study provides a comparative analysis of 6692 mutation data points on 34 aminergic GPCR subtypes, covering the chemical space of 540 unique ligands from mutagenesis experiments and information from experimentally determined structures of 52 distinct aminergic receptor-ligand complexes. The integrated analysis enables detailed investigation of structural receptor-ligand interactions and assessment of the transferability of combined binding mode and mutation data across ligand chemotypes and receptor subtypes. An overview is provided of the possibilities and limitations of using mutation data to guide the design of novel aminergic receptor ligands.

GPCR drug discovery: integrating solution NMR data with crystal and cryo-EM structures

The 826 G protein-coupled receptors (GPCRs) in the human proteome regulate key physiological processes and thus have long been attractive drug targets. With the crystal structures of more than 50 different human GPCRs determined over the past decade, an initial platform for structure-based rational design has been established for drugs that target GPCRs, which is currently being augmented with cryo-electron microscopy (cryo-EM) structures of higher-order GPCR complexes. Nuclear magnetic resonance (NMR) spectroscopy in solution is one of the key approaches for expanding this platform with dynamic features, which can be accessed at physiological temperature and with minimal modification of the wild-type GPCR covalent structures. Here, we review strategies for the use of advanced biochemistry and NMR techniques with GPCRs, survey projects in which crystal or cryo-EM structures have been complemented with NMR investigations and discuss the impact of this integrative approach on GPCR biology and drug discovery.

Activation Dynamics of the Neurotensin G Protein-Coupled Receptor 1

A replica-exchange protocol remarkably enhances the sampling of the activation dynamics of the neurotensin receptor type 1, a G protein-coupled receptor (GPCR) and important drug target. Our work highlights the dynamic communication between conformational changes of the agonist and the G protein-binding site, via contraction-oscillation of the orthosteric pocket. It also gives insights into the mechanism by which certain mutations diminish or stimulate activation. The replica-exchange protocol effectively enhances barrier crossing where standard brute-force molecular dynamics simulations fail. It is readily applicable to other GPCRs and represents a promising approach for virtual ligand screening, using the typical features of receptor activation as a benchmark.

Decoding Corticotropin-Releasing Factor Receptor Type 1 Crystal Structures

The structural analysis of class B G protein-coupled receptors (GPCR), cell surface proteins responding to peptide hormones, has until recently been restricted to the extracellular domain (ECD). Cor-ticotropin-releasing factor receptor type 1 (CRF1R) is a class B receptor mediating stress response and also considered a drug target for depression and anxiety. Here we report the crystal structure of the trans-membrane domain of human CRF1R in complex with the small-molecule antagonist CP-376395 in a hex-agonal setting with translational non-crystallographic symmetry. Molecular dynamics and metadynamics simulations on this novel structure and the existing TMD structure for CRF1R provides insight as to how the small molecule ligand gains access to the induced-fit allosteric binding site with implications for the observed selectivity against CRF2R. Furthermore, molecular dynamics simulations performed using a full-length receptor model point to key interactions between the ECD and extracellular loop 3 of the TMD providing insight into the full inactive state of multidomain class B GPCRs.

Structure of the complement C5a receptor bound to the extra-helical antagonist NDT9513727

The complement system is a crucial component of the host response to infection and tissue damage. Activation of the complement cascade generates anaphylatoxins including C5a and C3a. C5a exerts a pro-inflammatory effect via the G-protein-coupled receptor C5a anaphylatoxin chemotactic receptor 1 (C5aR1, also known as CD88) that is expressed on cells of myeloid origin. Inhibitors of the complement system have long been of interest as potential drugs for the treatment of diseases such as sepsis, rheumatoid arthritis, Crohn's disease and ischaemia-reperfusion injuries. More recently, a role of C5a in neurodegenerative conditions such as Alzheimer's disease has been identified. Peptide antagonists based on the C5a ligand have progressed to phase 2 trials in psoriasis and rheumatoid arthritis; however, these compounds exhibited problems with off-target activity, production costs, potential immunogenicity and poor oral bioavailability. Several small-molecule competitive antagonists for C5aR1, such as W-54011 and NDT9513727, have been identified by C5a radioligand-binding assays. NDT9513727 is a non-peptide inverse agonist of C5aR1, and is highly selective for the primate and gerbil receptors over those of other species. Here, to study the mechanism of action of C5a antagonists, we determine the structure of a thermostabilized C5aR1 (known as C5aR1 StaR) in complex with NDT9513727. We found that the small molecule bound between transmembrane helices 3, 4 and 5, outside the helical bundle. One key interaction between the small molecule and residue Trp213^5.49 seems to determine the species selectivity of the compound. The structure demonstrates that NDT9513727 exerts its inverse-agonist activity through an extra-helical mode of action.

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.

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.

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

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

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.

Approaching Pharmacological Space: Events and Components

With a view to introducing the concept of pharmacological space and its potential applications in investigating and predicting the toxic mechanisms of xenobiotics, this opening chapter describes the logical relations between conformational behavior, physicochemical properties and binding spaces, which are seen as the three key elements composing the pharmacological space. While the concept of conformational space is routinely used to encode molecular flexibility, the concepts of property spaces and, particularly, of binding spaces are more innovative. Indeed, their descriptors can find fruitful applications (a) in describing the dynamic adaptability a given ligand experiences when inserted into a specific environment, and (b) in parameterizing the flexibility a ligand retains when bound to a biological target. Overall, these descriptors can conveniently account for the often disregarded entropic factors and as such they prove successful when inserted in ligand- or structure-based predictive models. Notably, and although binding space parameters can clearly be derived from MD simulations, the chapter will illustrate how docking calculations, despite their static nature, are able to evaluate ligand's flexibility by analyzing several poses for each ligand. Such an approach, which represents the founding core of the binding space concept, can find various applications in which the related descriptors show an impressive enhancing effect on the statistical performances of the resulting predictive models.

Structural Features and Ligand Selectivity for 10 Intermediates in the Activation Process of β2-Adrenergic Receptor

It has already been suggested by researchers that there should be multiple intermediate states in the activation process for G-protein-coupled receptors (GPCRs). However, the intermediate states are very short-lived and hardly captured by the experiments, leading to very limited understanding of their structural features and drug efficacies. In this work, a novel joint strategy of targeted molecular dynamics simulation, conventional molecular dynamics simulation, and virtual screening is developed to address the problems. The results from 10 intermediate conformations obtained from the work reveal that the ligand pocket is very unstable and fluctuates between the inactive state and the active one in the case of ligand-free, in particular for ECL2 as a gate-keeper of the ligand-binding. The ligand-binding site could be stable in the active state with a small volume and a completely closed ECL2, only when the G-protein-binding region is fully activated. In addition, the activations of the ligand-binding pocket and G-protein-binding site are relatively independent and exhibit a loose allosteric coupling, which contributes to the existence of multiple intermediate conformations. Interestingly, the screening performance of the agonists does not increase on increasing the overall activity of the intermediate state, but is dependent on the activated extent of the ligand pocket. The receptor is prone to bind the agonist when closing ECL2 and reducing the ligand-binding pocket volume, whereas it is more favorable for binding the antagonist when opening ECL2 and increasing the pocket volume. These observations added to previous studies could help us better understand the activation mechanism of GPCRs and provide valuable information for drug design.

Evolutionary action and structural basis of the allosteric switch controlling β2AR functional selectivity

Functional selectivity of G-protein-coupled receptors is believed to originate from ligand-specific conformations that activate only subsets of signaling effectors. In this study, to identify molecular motifs playing important roles in transducing ligand binding into distinct signaling responses, we combined in silico evolutionary lineage analysis and structure-guided site-directed mutagenesis with large-scale functional signaling characterization and non-negative matrix factorization clustering of signaling profiles. Clustering based on the signaling profiles of 28 variants of the β₂-adrenergic receptor reveals three clearly distinct phenotypical clusters, showing selective impairments of either the Gi or βarrestin/endocytosis pathways with no effect on Gs activation. Robustness of the results is confirmed using simulation-based error propagation. The structural changes resulting from functionally biasing mutations centered around the DRY, NPxxY, and PIF motifs, selectively linking these micro-switches to unique signaling profiles. Our data identify different receptor regions that are important for the stabilization of distinct conformations underlying functional selectivity.

Ligand-induced biased signaling is thought to result in part from ligand-specific receptor conformations that cause the engagement of distinct effectors. Here the authors trace and evaluate the impact of mutations of the β2–adrenergic receptor on multiple signaling outputs to provide structural-level insight into the determinants of GPCR functional selectivity.

Dynamical Binding Modes Determine Agonistic and Antagonistic Ligand Effects in the Prostate-Specific G-Protein Coupled Receptor (PSGR)

We analysed the ligand-based activation mechanism of the prostate-specific G-protein coupled receptor (PSGR), which is an olfactory receptor that mediates cellular growth in prostate cancer cells. Furthermore, it is an olfactory receptor with a known chemically near identic antagonist/agonist pair, α- and β-ionone. Using a combined theoretical and experimental approach, we propose that this receptor is activated by a ligand-induced rearrangement of a protein-internal hydrogen bond network. Surprisingly, this rearrangement is not induced by interaction of the ligand with the network, but by dynamic van der Waals contacts of the ligand with the involved amino acid side chains, altering their conformations and intraprotein connectivity. Ligand recognition in this GPCR is therefore highly stereo selective, but seemingly lacks any ligand recognition via polar contacts. A putative olfactory receptor-based drug design scheme will have to take this unique mode of protein/ligand action into account.

Extending the Structural View of Class B GPCRs

The secretin-like class B family of G protein-coupled receptors (GPCRs) are key players in hormonal homeostasis. Recent structures of various receptors in complex with a variety of orthosteric and allosteric ligands provide fundamental new insights into the function and mechanism of class B GPCRs, including: (i) ligand-induced changes in the relative orientation of the extracellular and transmembrane receptor domains; (ii) intramolecular interaction networks that stabilize conformational changes to accommodate intracellular G protein binding; and (iii) allosteric modulation of receptor activation. This review provides a comprehensive analysis of the structural, biochemical, and pharmacological data on class B GPCRs for understanding ligand-receptor interaction and modulation mechanisms and assessing the potential implications for drug discovery for the secretin-like GPCR family.

Agonist Binding to Chemosensory Receptors: A Systematic Bioinformatics Analysis

Human G-protein coupled receptors (hGPCRs) constitute a large and highly pharmaceutically relevant membrane receptor superfamily. About half of the hGPCRs' family members are chemosensory receptors, involved in bitter taste and olfaction, along with a variety of other physiological processes. Hence these receptors constitute promising targets for pharmaceutical intervention. Molecular modeling has been so far the most important tool to get insights on agonist binding and receptor activation. Here we investigate both aspects by bioinformatics-based predictions across all bitter taste and odorant receptors for which site-directed mutagenesis data are available. First, we observe that state-of-the-art homology modeling combined with previously used docking procedures turned out to reproduce only a limited fraction of ligand/receptor interactions inferred by experiments. This is most probably caused by the low sequence identity with available structural templates, which limits the accuracy of the protein model and in particular of the side-chains' orientations. Methods which transcend the limited sampling of the conformational space of docking may improve the predictions. As an example corroborating this, we review here multi-scale simulations from our lab and show that, for the three complexes studied so far, they significantly enhance the predictive power of the computational approach. Second, our bioinformatics analysis provides support to previous claims that several residues, including those at positions 1.50, 2.50, and 7.52, are involved in receptor activation.

Activation processes in ligand-activated G protein-coupled receptors: A case study of the adenosine A2A receptor

Here we review concepts related to an ensemble description of G-protein-coupled receptors (GPCRs). The ensemble is characterized by both inactive and active states, whose equilibrium populations and exchange rates depend sensitively on ligand, environment, and allosteric factors. This review focuses on the adenosine A₂ receptor (A_2A R), a prototypical class A GPCR. ¹⁹ F Nuclear Magnetic Resonance (NMR) studies show that apo A_2A R is characterized by a broad ensemble of conformers, spanning inactive to active states, and resembling states defined earlier for rhodopsin. In keeping with ideas associated with a conformational selection mechanism, addition of agonist serves to allosterically restrict the overall degrees of freedom at the G protein binding interface and bias both states and functional dynamics to facilitate G protein binding and subsequent activation. While the ligand does not necessarily "induce" activation, it does bias sampling of states, increase the cooperativity of the activation process and thus, the lifetimes of functional activation intermediates, while restricting conformational dynamics to that needed for activation.

Receptor activity-modifying protein dependent and independent activation mechanisms in the coupling of calcitonin gene-related peptide and adrenomedullin receptors to Gs

Calcitonin gene-related peptide (CGRP) or adrenomedullin (AM) receptors are heteromers of the calcitonin receptor-like receptor (CLR), a class B G protein-coupled receptor, and one of three receptor activity-modifying proteins (RAMPs). How CGRP and AM activate CLR and how this process is modulated by RAMPs is unclear. We have defined how CGRP and AM induce Gs-coupling in CLR-RAMP heteromers by measuring the effect of targeted mutagenesis in the CLR transmembrane domain on cAMP production, modeling the active state conformations of CGRP and AM receptors in complex with the Gs C-terminus and conducting molecular dynamics simulations in an explicitly hydrated lipidic bilayer. The largest effects on receptor signaling were seen with H295A^5.40b, I298A^5.43b, L302A^5.47b, N305A^5.50b, L345A^6.49b and E348A^6.52b, F349A^6.53b and H374A^7.47b (class B numbering in superscript). Many of these residues are likely to form part of a group in close proximity to the peptide binding site and link to a network of hydrophilic and hydrophobic residues, which undergo rearrangements to facilitate Gs binding. Residues closer to the extracellular loops displayed more pronounced RAMP or ligand-dependent effects. Mutation of H374^7.47b to alanine increased AM potency 100-fold in the CGRP receptor. The molecular dynamics simulation showed that TM5 and TM6 pivoted around TM3. The data suggest that hydrophobic interactions are more important for CLR activation than other class B GPCRs, providing new insights into the mechanisms of activation of this class of receptor. Furthermore the data may aid in the understanding of how RAMPs modulate the signaling of other class B GPCRs.

Analysis of positive and negative allosteric modulation in metabotropic glutamate receptors 4 and 5 with a dual ligand

As class C GPCRs and regulators of synaptic activity, human metabotropic glutamate receptors (mGluRs) 4 and 5 are prime targets for allosteric modulation, with mGlu5 inhibition or mGlu4 stimulation potentially treating conditions like chronic pain and Parkinson’s disease. As an allosteric modulator that can bind both receptors, 2-Methyl-6-(phenylethynyl)pyridine (MPEP) is able to negatively modulate mGlu5 or positively modulate mGlu4. At a structural level, how it elicits these responses and how mGluRs undergo activation is unclear. Here, we employ homology modelling and 30 µs of atomistic molecular dynamics (MD) simulations to probe allosteric conformational change in mGlu4 and mGlu5, with and without docked MPEP. Our results identify several structural differences between mGlu4 and mGlu5, as well as key differences responsible for MPEP-mediated positive and negative allosteric modulation, respectively. A novel mechanism of mGlu4 activation is revealed, which may apply to all mGluRs in general. This involves conformational changes in TM3, TM4 and TM5, separation of intracellular loop 2 (ICL2) from ICL1/ICL3, and destabilization of the ionic-lock. On the other hand, mGlu5 experiences little disturbance when MPEP binds, maintaining its inactive state with reduced conformational fluctuation. In addition, when MPEP is absent, a lipid molecule can enter the mGlu5 allosteric pocket.

Identification of Histamine H3 Receptor Ligands Using a New Crystal Structure Fragment-based Method

Virtual screening offers an efficient alternative to high-throughput screening in the identification of pharmacological tools and lead compounds. Virtual screening is typically based on the matching of target structures or ligand pharmacophores to commercial or in-house compound catalogues. This study provides the first proof-of-concept for our recently reported method where pharmacophores are instead constructed based on the inference of residue-ligand fragments from crystal structures. We demonstrate its unique utility for G protein-coupled receptors, which represent the largest families of human membrane proteins and drug targets. We identified five neutral antagonists and one inverse agonist for the histamine H₃ receptor with potencies of 0.7-8.5 μM in a recombinant receptor cell-based inositol phosphate accumulation assay and validated their activity using a radioligand competition binding assay. H₃ receptor antagonism is of large therapeutic value and our ligands could serve as starting points for further lead optimisation. The six ligands exhibit four chemical scaffolds, whereof three have high novelty in comparison to the known H₃ receptor ligands in the ChEMBL database. The complete pharmacophore fragment library is freely available through the GPCR database, GPCRdb, allowing the successful application herein to be repeated for most of the 285 class A GPCR targets. The method could also easily be adapted to other protein families.

Drug Binding Poses Relate Structure with Efficacy in the μ Opioid Receptor

The μ-opioid receptor (MOPr) is a clinically important G protein-coupled receptor that couples to Gi/o proteins and arrestins. At present, the receptor conformational changes that occur following agonist binding and activation are poorly understood. This study has employed molecular dynamics simulations to investigate the binding mode and receptor conformational changes induced by structurally similar opioid ligands of widely differing intrinsic agonist efficacy, norbuprenorphine, buprenorphine, and diprenorphine. Bioluminescence resonance energy transfer assays for Gi activation and arrestin-3 recruitment in human embryonic kidney 293 cells confirmed that norbuprenorphine is a high efficacy agonist, buprenorphine a low efficacy agonist, and diprenorphine an antagonist at the MOPr. Molecular dynamics simulations revealed that these ligands adopt distinct binding poses and engage different subsets of residues, despite sharing a common morphinan scaffold. Notably, norbuprenorphine interacted with sodium ion-coordinating residues W2936.48 and N1503.35, whilst buprenorphine and diprenorphine did not. Principal component analysis of the movements of the receptor transmembrane domains showed that the buprenorphine-bound receptor occupied a distinct set of conformations to the norbuprenorphine-bound receptor. Addition of an allosteric sodium ion caused the receptor and ligand to adopt an inactive conformation. The differences in ligand-residue interactions and receptor conformations observed here may underlie the differing efficacies for cellular signalling outputs for these ligands.

Crystal structure of the GLP-1 receptor bound to a peptide agonist

Glucagon-like peptide 1 (GLP-1) regulates glucose homeostasis through the control of insulin release from the pancreas. GLP-1 peptide agonists are efficacious drugs for the treatment of diabetes. To gain insight into the molecular mechanism of action of GLP-1 peptides, here we report the crystal structure of the full-length GLP-1 receptor bound to a truncated peptide agonist. The peptide agonist retains an α-helical conformation as it sits deep within the receptor-binding pocket. The arrangement of the transmembrane helices reveals hallmarks of an active conformation similar to that observed in class A receptors. Guided by this structural information, we design peptide agonists with potent in vivo activity in a mouse model of diabetes.

Activation mechanisms of the first sphingosine-1-phosphate receptor

Activation of the first sphingosine-1-phosphate receptor (S1PR₁ ) promotes permeability of the blood brain barrier, astrocyte and neuronal protection, and lymphocyte egress from secondary lymphoid tissues. Although an agonist often activates the S1PR₁ , the receptor exhibits high levels of basal activity. In this study, we performed long-timescale molecular dynamics and accelerated molecular dynamics (aMD) simulations to investigate activation mechanisms of the ligand-free (apo) S1PR₁ . In the aMD enhanced sampling simulations, we observed four independent events of activation, which is characterized by close interaction between Y311^7.53 and Y221^5.58 and increased distance between the intracellular ends of transmembrane (TM) helices 3 and 6. Although TM helices TM3, TM6, TM5 and, TM7 are associated with GPCR activation, we discovered that their movements are not necessarily correlated during activation. Instead, TM5 showed a decreased correlation with each of these regions during activation. During activation of the apo receptor, Y221^5.58 and Y311^7.53 became more solvated, because a water channel formed in the intracellular pocket. Additionally, a lipid molecule repeatedly entered the receptor between the extracellular ends of TM1 and TM7, providing important insights into the pathway of ligand entry into the S1PR₁ .

Conformational equilibria of light-activated rhodopsin in nanodiscs

Conformational equilibria of G-protein-coupled receptors (GPCRs) are intimately involved in intracellular signaling. Here conformational substates of the GPCR rhodopsin are investigated in micelles of dodecyl maltoside (DDM) and in phospholipid nanodiscs by monitoring the spatial positions of transmembrane helices 6 and 7 at the cytoplasmic surface using site-directed spin labeling and double electron-electron resonance spectroscopy. The photoactivated receptor in DDM is dominated by one conformation with weak pH dependence. In nanodiscs, however, an ensemble of pH-dependent conformational substates is observed, even at pH 6.0 where the MIIbH⁺ form defined by proton uptake and optical spectroscopic methods is reported to be the sole species present in native disk membranes. In nanodiscs, the ensemble of substates in the photoactivated receptor spontaneously decays to that characteristic of the inactive state with a lifetime of ∼16 min at 20 °C. Importantly, transducin binding to the activated receptor selects a subset of the ensemble in which multiple substates are apparently retained. The results indicate that in a native-like lipid environment rhodopsin activation is not analogous to a simple binary switch between two defined conformations, but the activated receptor is in equilibrium between multiple conformers that in principle could recognize different binding partners.

Structural heterogeneity of the μ-opioid receptor's conformational ensemble in the apo state

G-protein coupled receptors (GPCRs) are the largest and most pharmaceutically relevant family of membrane proteins. Here, fully unbiased, enhanced sampling simulations of a constitutively active mutant (CAM) of a class A GPCR, the μ-opioid receptor (μOR), demonstrates repeated transitions between the inactive (IS) and active-like (AS-L) states. The interconversion features typical activation/inactivation patterns involving established conformational rearrangements of conserved residues. By contrast, wild-type μOR remains in IS during the same course of simulation, consistent with the low basal activity of the protein. The simulations point to an important role of residue W293^6.48 at the "toggle switch" in the mutation-induced constitutive activation. Such role has been already observed for other CAMs of class A GPCRs. We also find a significantly populated intermediate state, rather similar to IS. Based on the remarkable accord between simulations and experiments, we suggest here that this state, which has escaped so far experimental characterization, might constitute an early step in the activation process of the apo μOR CAM.

Functional roles of T3.37 and S5.46 in the activation mechanism of the dopamine D1 receptor

The activation mechanism of dopamine receptors is unknown. The amino acids S5.42, S5.43, and S5.46 located in helix 5 appear to be crucial, but their specific roles in receptor activation have not been studied. We modeled the D1 dopamine receptor using the crystal structures of the D3 dopamine and β2 adrenergic receptors. Molecular dynamics simulations show that the interaction of dopamine with the D1 receptor leads to the formation of a hydrogen-bond network with its catechol group and helices 3, 5, and 6, including water molecules. The para hydroxyl group of dopamine binds directly to S5.42 and N6.55, the latter also interacting with S5.43. Unexpectedly, S5.46 does not interact directly with the catechol; instead, it interacts through a water molecule with S5.42 and directly with T3.37. The formation of this hydrogen-bond network, part of which was previously observed in docking studies with dopamine agonists, triggers the opening of the E6.30-R3.60 ionic lock associated with the activation of GPCRs. These changes do not occur in the unbonded (apo) receptor or when it is in a complex with the antagonist 3-methoxy-5,6,7,8,9,14-hexahydrodibenz[d,g]azecine. Our results provide valuable insight into the T3.37-S5.46-water-S5.43-ligand interaction, which may be crucial to the activation of the D1 dopamine receptor and should be considered during the design of novel agonists. Graphical Abstract General representation of the relationship between the formation of the HBN and the opening of the R3.50-E6.30 ionic lock.

7×7 RMSD matrix: A new method for quantitative comparison of the transmembrane domain structures in the G-protein coupled receptors

The G-protein coupled receptors (GPCRs) share a conserved heptahelical fold in the transmembrane (TM) region, but the exact arrangements of the seven TM helices vary with receptors and their activation states. The differences or the changes have been observed in the experimentally solved structures, but have not been systematically and quantitatively investigated due to lack of suitable methods. In this work, we describe a novel method, called 7×7 RMSD matrix that is proposed specifically for comparing the characteristic 7TM bundle structures of GPCRs. Compared to the commonly used overall TM bundle RMSD as a single parameter, a 7×7 RMSD matrix contains 49 parameters, which reveal changes of the relative orientations of the seven TMs. We demonstrate the novelty and advantages of this method by tackling two problems that are challenging for the existing methods. With this method, we are able to identify and quantify the helix movements in the activated receptor structures and reveal structural conservation and divergence as well as the structural relationships of different GPCRs in terms of the relative orientations of the seven TMs.

Unifying view of mechanical and functional hotspots across class A GPCRs

G protein-coupled receptors (GPCRs) are the largest superfamily of signaling proteins. Their activation process is accompanied by conformational changes that have not yet been fully uncovered. Here, we carry out a novel comparative analysis of internal structural fluctuations across a variety of receptors from class A GPCRs, which currently has the richest structural coverage. We infer the local mechanical couplings underpinning the receptors’ functional dynamics and finally identify those amino acids whose virtual deletion causes a significant softening of the mechanical network. The relevance of these amino acids is demonstrated by their overlap with those known to be crucial for GPCR function, based on static structural criteria. The differences with the latter set allow us to identify those sites whose functional role is more clearly detected by considering dynamical and mechanical properties. Of these sites with a genuine mechanical/dynamical character, the top ranking is amino acid 7x52, a previously unexplored, and experimentally verifiable key site for GPCR conformational response to ligand binding.

The biological functionality of several receptors and enzymes depends on their capability to sustain large-scale structural fluctuations and adopt different conformational states in response to ligand binding. This is the case for G protein-coupled receptors (GPCRs), the largest superfamily of signaling proteins in mammals and a primary pharmaceutical target. To better understand the functional dynamics of GPCRs, we have analysed the inter-residue distance variations across the available structures for several receptors of the rhodopsin-like family (class A). We first reconstructed the network of mechanical, rigid-like couplings between nearby amino acids and then identified those acting as dynamical/mechanical hubs. These were the sites whose virtual removal led to a significant softening of the overall mechanical network. After validating the biological relevance of these sites by comparison against known key functional sites, we singled out those regions which emerge as prominent mechanical hubs and yet have an otherwise still unknown functional role. The most relevant of such novel putative functional sites, which could be probed by mutagenesis experiments, is at interface of two transmembrane helices and we expect it to be crucial for assisting GPCRs conformational response to agonist binding.

Gonadotropin-Releasing Hormone (GnRH) Receptor Structure and GnRH Binding

Gonadotropin-releasing hormone (GnRH) regulates reproduction. The human GnRH receptor lacks a cytoplasmic carboxy-terminal tail but has amino acid sequence motifs characteristic of rhodopsin-like, class A, G protein-coupled receptors (GPCRs). This review will consider how recent descriptions of X-ray crystallographic structures of GPCRs in inactive and active conformations may contribute to understanding GnRH receptor structure, mechanism of activation and ligand binding. The structures confirmed that ligands bind to variable extracellular surfaces, whereas the seven membrane-spanning α-helices convey the activation signal to the cytoplasmic receptor surface, which binds and activates heterotrimeric G proteins. Forty non-covalent interactions that bridge topologically equivalent residues in different transmembrane (TM) helices are conserved in class A GPCR structures, regardless of activation state. Conformation-independent interhelical contacts account for a conserved receptor protein structure and their importance in the GnRH receptor structure is supported by decreased expression of receptors with mutations of residues in the network. Many of the GnRH receptor mutations associated with congenital hypogonadotropic hypogonadism, including the Glu^2.53(90) Lys mutation, involve amino acids that constitute the conserved network. Half of the ~250 intramolecular interactions in GPCRs differ between inactive and active structures. Conformation-specific interhelical contacts depend on amino acids changing partners during activation. Conserved inactive conformation-specific contacts prevent receptor activation by stabilizing proximity of TM helices 3 and 6 and a closed G protein-binding site. Mutations of GnRH receptor residues involved in these interactions, such as Arg^3.50(139) of the DRY/S motif or Tyr^7.53(323) of the N/DPxxY motif, increase or decrease receptor expression and efficiency of receptor coupling to G protein signaling, consistent with the native residues stabilizing the inactive GnRH receptor structure. Active conformation-specific interhelical contacts stabilize an open G protein-binding site. Progress in defining the GnRH-binding site has recently slowed, with evidence that Tyr^6.58(290) contacts Tyr⁵ of GnRH, whereas other residues affect recognition of Trp³ and Gly¹⁰NH₂. The surprisingly consistent observations that GnRH receptor mutations that disrupt GnRH binding have less effect on "conformationally constrained" GnRH peptides may now be explained by crystal structures of agonist-bound peptide receptors. Analysis of GPCR structures provides insight into GnRH receptor function.

Structural Analysis of the Histamine H1 Receptor

The crystal structure of the human histamine H₁ receptor (H₁R) has been determined in complex with its inverse agonist doxepin, a first-generation antihistamine. The crystal structure showed that doxepin sits deeply inside the ligand-binding pocket and predominantly interacts with residues highly conserved among other aminergic receptors. This binding mode is considered to result in the low selectivity of the first-generation antihistamines for H₁R. The crystal structure also revealed the mechanism of receptor inactivation by the inverse agonist doxepin. On the other hand, the crystal structure elucidated the anion-binding site near the extracellular portion of the receptor. This site consists of residues not conserved among other aminergic receptors, which are specific for H₁R. Docking simulation and biochemical experimentation demonstrated that a carboxyl group on the second-generation antihistamines interacts with the anion-binding site. These results imply that the anion-binding site is a key site for the development of highly selective antihistamine drugs.

Intracellular allosteric antagonism of the CCR9 receptor

Chemokines and their G-protein-coupled receptors play a diverse role in immune defence by controlling the migration, activation and survival of immune cells. They are also involved in viral entry, tumour growth and metastasis and hence are important drug targets in a wide range of diseases. Despite very significant efforts by the pharmaceutical industry to develop drugs, with over 50 small-molecule drugs directed at the family entering clinical development, only two compounds have reached the market: maraviroc (CCR5) for HIV infection and plerixafor (CXCR4) for stem-cell mobilization. The high failure rate may in part be due to limited understanding of the mechanism of action of chemokine antagonists and an inability to optimize compounds in the absence of structural information. CC chemokine receptor type 9 (CCR9) activation by CCL25 plays a key role in leukocyte recruitment to the gut and represents a therapeutic target in inflammatory bowel disease. The selective CCR9 antagonist vercirnon progressed to phase 3 clinical trials in Crohn's disease but efficacy was limited, with the need for very high doses to block receptor activation. Here we report the crystal structure of the CCR9 receptor in complex with vercirnon at 2.8 Å resolution. Remarkably, vercirnon binds to the intracellular side of the receptor, exerting allosteric antagonism and preventing G-protein coupling. This binding site explains the need for relatively lipophilic ligands and describes another example of an allosteric site on G-protein-coupled receptors that can be targeted for drug design, not only at CCR9, but potentially extending to other chemokine receptors.