Computational methods in drug design: modeling G protein-coupled receptor monomers, dimers, and oligomers

Computational Methods to Target Protein-Protein Interactions

The chapter emphasizes the importance of understanding protein-protein interactions in cellular mechanisms and highlights the role of computational modeling in predicting these interactions. It discusses sequence-based approaches such as evolutionary trace (ET), correlated mutation analysis (CMA), and subtractive correlated mutation (SCM) for identifying crucial amino acid residues, considering interface conservation or evolutionary changes. The chapter also explores methods like differential ET, hidden-site class model, and spatial cluster detection (SCD) for interface specificity and spatial clustering. Furthermore, it examines approaches combining structural and sequential methodologies and evaluates modeled predictions through initiatives like critical assessment of prediction of interactions (CAPRI). Additionally, the chapter provides an overview of various software programs used for molecular docking, detailing their search, sampling, refinement and scoring stages, along with innovative techniques and tools like normal mode analysis (NMA) and adaptive Poisson-Boltzmann solver (APBS) for electrostatic calculations. These computational and experimental approaches are crucial for unraveling protein-protein interactions and aid in developing potential therapeutics for various diseases.

Structural insight into apelin receptor-G protein stoichiometry

The technique of cryogenic-electron microscopy (cryo-EM) has revolutionized the field of membrane protein structure and function with a focus on the dominantly observed molecular species. This report describes the structural characterization of a fully active human apelin receptor (APJR) complexed with heterotrimeric G protein observed in both 2:1 and 1:1 stoichiometric ratios. We use cryo-EM single-particle analysis to determine the structural details of both species from the same sample preparation. Protein preparations, in the presence of the endogenous peptide ligand ELA or a synthetic small molecule, both demonstrate these mixed stoichiometric states. Structural differences in G protein engagement between dimeric and monomeric APJR suggest a role for the stoichiometry of G protein-coupled receptor- (GPCR-)G protein coupling on downstream signaling and receptor pharmacology. Furthermore, a small, hydrophobic dimer interface provides a starting framework for additional class A GPCR dimerization studies. Together, these findings uncover a mechanism of versatile regulation through oligomerization by which GPCRs can modulate their signaling.

GPR6 Structural Insights: Homology Model Construction and Docking Studies

GPR6 is an orphan G protein-coupled receptor that has been associated with the cannabinoid family because of its recognition of a sub-set of cannabinoid ligands. The high abundance of GPR6 in the central nervous system, along with high constitutive activity and a link to several neurodegenerative diseases make GPR6 a promising biological target. In fact, diverse research groups have demonstrated that GPR6 represents a possible target for the treatment of neurodegenerative disorders such as Parkinson's disease, Alzheimer's disease, and Huntington's disease. Several patents have claimed the use of a wide range of pyrazine derivatives as GPR6 inverse agonists for the treatment of Parkinson's disease symptoms and other dyskinesia syndromes. However, the full pharmacological importance of GPR6 has not yet been fully explored due to the lack of high potency, readily available ligands targeting GPR6. The long-term goal of the present study is to develop such ligands. In this paper, we describe our initial steps towards this goal. A human GPR6 homology model was constructed using a suite of computational techniques. This model permitted the identification of unique GPR6 structural features and the exploration of the GPR6 binding crevice. A subset of patented pyrazine analogs were docked in the resultant GPR6 inactive state model to validate the model, rationalize the structure-activity relationships from the reported patents and identify the key residues in the binding crevice for ligand recognition. We will take this structural knowledge into the next phase of GPR6 project, in which scaffold hopping will be used to design new GPR6 ligands.

Prediction and targeting of GPCR oligomer interfaces

GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.

Towards A Molecular Understanding of The Cannabinoid Related Orphan Receptor GPR18: A Focus on Its Constitutive Activity

The orphan G-protein coupled receptor (GPCR), GPR18, has been recently proposed as a potential member of the cannabinoid family as it recognizes several endogenous, phytogenic, and synthetic cannabinoids. Potential therapeutic applications for GPR18 include intraocular pressure, metabolic disorders, and cancer. GPR18 has been reported to have high constitutive activity, i.e., activation/signaling occurs in the absence of an agonist. This activity can be reduced significantly by the A3.39N mutation. At the intracellular (IC) ends of (transmembrane helices) TMH3 and TMH6 in GPCRs, typically, a pair of oppositely charged amino acids form a salt bridge called the "ionic lock". Breaking of this salt bridge creates an IC opening for coupling with G protein. The GPR18 "ionic lock" residues (R3.50/S6.33) can form only a hydrogen bond. In this paper, we test the hypothesis that the high constitutive activity of GPR18 is due to the weakness of its "ionic lock" and that the A3.39N mutation strengthens this lock. To this end, we report molecular dynamics simulations of wild-type (WT) GPR18 and the A3.39N mutant in fully hydrated (POPC) phophatidylcholine lipid bilayers. Results suggest that in the A3.39N mutant, TMH6 rotates and brings R3.50 and S6.33 closer together, thus strengthening the GPR18 "ionic lock".

Remodeling our concept of chemokine receptor function: From monomers to oligomers

The chemokines direct leukocyte recruitment in both homeostatic and inflammatory conditions, and are therefore critical for immune reactions. By binding to members of the class A G protein-coupled receptors, the chemokines play an essential role in numerous physiological and pathological processes. In the last quarter century, the field has accumulated much information regarding the implications of these molecules in different immune processes, as well as mechanistic insight into the signaling events activated through their binding to their receptors. Here, we will focus on chemokine receptors and how new methodological approaches have underscored the role of their conformations in chemokine functions. Advances in biophysical-based techniques show that chemokines and their receptors act in very complex networks and therefore should not be considered isolated entities. In this regard, the chemokine receptors can form homo- and heterodimers as well as oligomers at the cell surface. These findings are changing our view as to how chemokines influence cell biology, identify partners that regulate chemokine function, and open new avenues for therapeutic intervention.

Methods for the Development of In Silico GPCR Models

The Reggio group has constructed computer models of the inactive and G-protein-activated states of the cannabinoid CB1 and CB2 receptors, as well as, several orphan receptors that recognize a subset of cannabinoid compounds, including GPR55 and GPR18. These models have been used to design ligands, mutations, and covalent labeling studies. The resultant second-generation models have been used to design ligands with improved affinity, efficacy, and subtype selectivity. Herein, we provide a guide for the development of GPCR models using the most recent orphan receptor studied in our lab, GPR3. GPR3 is an orphan receptor that belongs to the Class A family of G-protein-coupled receptors. It shares high sequence similarity with GPR6, GPR12, the lysophospholipid receptors, and the cannabinoid receptors. GPR3 is predominantly expressed in mammalian brain and oocytes and it is known as a Gαs-coupled receptor activated constitutively in cells. GPR3 represents a possible target for the treatment of different pathological conditions such as Alzheimer's disease, oocyte maturation, or neuropathic pain. However, the lack of potent and selective GPR3 ligands is delaying the exploitation of this promising therapeutic target. In this context, we aim to develop a homology model that helps us to elucidate the structural determinants governing ligand-receptor interactions at GPR3. In this chapter, we detail the methods and rationale behind the construction of the GPR3 active-and inactive-state models. These homology models will enable the rational design of novel ligands, which may serve as research tools for further understanding of the biological role of GPR3.

Methods used to study the oligomeric structure of G-protein-coupled receptors

G-protein-coupled receptors (GPCRs), which constitute the largest family of cell surface receptors, were originally thought to function as monomers, but are now recognized as being able to act in a wide range of oligomeric states and indeed, it is known that the oligomerization state of a GPCR can modulate its pharmacology and function. A number of experimental techniques have been devised to study GPCR oligomerization including those based upon traditional biochemistry such as blue-native PAGE (BN-PAGE), co-immunoprecipitation (Co-IP) and protein-fragment complementation assays (PCAs), those based upon resonance energy transfer, FRET, time-resolved FRET (TR-FRET), FRET spectrometry and bioluminescence resonance energy transfer (BRET). Those based upon microscopy such as FRAP, total internal reflection fluorescence microscopy (TIRFM), spatial intensity distribution analysis (SpIDA) and various single molecule imaging techniques. Finally with the solution of a growing number of crystal structures, X-ray crystallography must be acknowledged as an important source of discovery in this field. A different, but in many ways complementary approach to the use of more traditional experimental techniques, are those involving computational methods that possess obvious merit in the study of the dynamics of oligomer formation and function. Here, we summarize the latest developments that have been made in the methods used to study GPCR oligomerization and give an overview of their application.

Apelin receptor homodimer-oligomers revealed by single-molecule imaging and novel G protein-dependent signaling

The apelin receptor (APJ) belongs to family A of the G protein-coupled receptors (GPCRs) and is a potential pharmacotherapeutic target for heart failure, hypertension, and other cardiovascular diseases. There is evidence APJ heterodimerizes with other GPCRs; however, the existence of APJ homodimers and oligomers remains to be investigated. Here, we measured APJ monomer-homodimer-oligomer interconversion by monitoring APJ dynamically on cells and compared their proportions, spatial arrangement, and mobility using total internal reflection fluorescence microscopy, resonance energy transfer, and proximity biotinylation. In cells with <0.3 receptor particles/μm², approximately 60% of APJ molecules were present as dimers or oligomers. APJ dimers were present on the cell surface in a dynamic equilibrium with constant formation and dissociation of receptor complexes. Furthermore, we applied interference peptides and MALDI-TOF mass spectrometry to confirm APJ homo-dimer and explore the dimer-interfaces. Peptides corresponding to transmembrane domain (TMD)1, 2, 3, and 4, but not TMD5, 6, and 7, disrupted APJ dimerization. APJ mutants in TMD1 and TMD2 also decreased bioluminescence resonance energy transfer of APJ dimer. APJ dimerization resulted in novel functional characteristics, such as a distinct G-protein binding profile and cell responses after agonist stimulation. Thus, dimerization may serve as a unique mechanism for fine-tuning APJ-mediated functions.

Chromenopyrazole, a Versatile Cannabinoid Scaffold with in Vivo Activity in a Model of Multiple Sclerosis

A combination of molecular modeling and structure-activity relationship studies has been used to fine-tune CB2 selectivity in the chromenopyrazole ring, a versatile CB1/CB2 cannabinoid scaffold. Thus, a series of 36 new derivatives covering a wide range of structural diversity has been synthesized, and docking studies have been performed for some of them. Biological evaluation of the new compounds includes, among others, cannabinoid binding assays, functional studies, and surface plasmon resonance measurements. The most promising compound [43 (PM226)], a selective and potent CB2 agonist isoxazole derivative, was tested in the acute phase of Theiler's murine encephalomyelitis virus-induced demyelinating disease (TMEV-IDD), a well-established animal model of primary progressive multiple sclerosis. Compound 43 dampened neuroinflammation by reducing microglial activation in the TMEV.

Difference and Influence of Inactive and Active States of Cannabinoid Receptor Subtype CB2: From Conformation to Drug Discovery

Cannabinoid receptor 2 (CB2), a G protein-coupled receptor (GPCR), is a promising target for the treatment of neuropathic pain, osteoporosis, immune system, cancer, and drug abuse. The lack of an experimental three-dimensional CB2 structure has hindered not only the development of studies of conformational differences between the inactive and active CB2 but also the rational discovery of novel functional compounds targeting CB2. In this work, we constructed models of both inactive and active CB2 by homology modeling. Then we conducted two comparative 100 ns molecular dynamics (MD) simulations on the two systems-the active CB2 bound with both the agonist and G protein and the inactive CB2 bound with inverse agonist-to analyze the conformational difference of CB2 proteins and the key residues involved in molecular recognition. Our results showed that the inactive CB2 and the inverse agonist remained stable during the MD simulation. However, during the MD simulations, we observed dynamical details about the breakdown of the "ionic lock" between R131(3.50) and D240(6.30) as well as the outward/inward movements of transmembrane domains of the active CB2 that bind with G proteins and agonist (TM5, TM6, and TM7). All of these results are congruent with the experimental data and recent reports. Moreover, our results indicate that W258(6.48) in TM6 and residues in TM4 (V164(4.56)-L169(4.61)) contribute greatly to the binding of the agonist on the basis of the binding energy decomposition, while residues S180-F183 in extracellular loop 2 (ECL2) may be of importance in recognition of the inverse agonist. Furthermore, pharmacophore modeling and virtual screening were carried out for the inactive and active CB2 models in parallel. Among all 10 hits, two compounds exhibited novel scaffolds and can be used as novel chemical probes for future studies of CB2. Importantly, our studies show that the hits obtained from the inactive CB2 model mainly act as inverse agonist(s) or neutral antagonist(s) at low concentration. Moreover, the hit from the active CB2 model also behaves as a neutral antagonist at low concentration. Our studies provide new insight leading to a better understanding of the structural and conformational differences between two states of CB2 and illuminate the effects of structure on virtual screening and drug design.

Computational approaches for modeling GPCR dimerization

Growing experimental evidences suggest that dimerization and oligomerization are important for G Protein- Coupled Receptors (GPCRs) function. The detailed structural information of dimeric/oligomeric GPCRs would be very important to understand their function. Although it is encouraging that recently several experimental GPCR structures in oligomeric forms have appeared, experimental determination of GPCR structures in oligomeric forms is still a big challenge, especially in mimicking the membrane environment. Therefore, development of computational approaches to predict dimerization of GPCRs will be highly valuable. In this review, we summarize computational approaches that have been developed and used for modeling of GPCR dimerization. In addition, we introduce a novel two-dimensional Brownian Dynamics based protein docking approach, which we have recently adapted, for GPCR dimer prediction.

Description of a bivalent cannabinoid ligand with hypophagic properties

A series of bivalent cannabinoid ligands is proposed. The synthesis of double amides based on the rimonabant structure separated by an alkyl chain and the evaluation of their affinities for cannabinoid receptors are reported. The data of 4d confirmed that a bivalent structure is a suitable scaffold for CB1 cannabinoid receptor binding. The compound 4d was selected for in vitro and in vivo pharmacological evaluations. Moreover, intraperitoneal administration of 4d to food-deprived rats resulted in a dose-dependent inhibition of feeding that was maintained up to 240 min.

Prediction of the Human EP1 Receptor Binding Site by Homology Modeling and Molecular Dynamics Simulation

The prostanoid receptor EP1 is a G-protein-coupled receptor (GPCR) known to be involved in a variety of pathological disorders such as pain, fever and inflammation. These receptors are important drug targets, but design of subtype specific agonists and antagonists has been partially hampered by the absence of three-dimensional structures for these receptors. To understand the molecular interactions of the PGE2, an endogen ligand, with the EP1 receptor, a homology model of the human EP1 receptor (hEP1R) with all connecting loops was constructed from the 2.6 Å resolution crystal structure (PDB code: 1L9H) of bovine rhodopsin. The initial model generated by MODELLER was subjected to molecular dynamics simulation to assess quality of the model. Also, a step by step ligand-supported model refinement was performed, including initial docking of PGE2 and iloprost in the putative binding site, followed by several rounds of energy minimizations and molecular dynamics simulations. Docking studies were performed for PGE2 and some other related compounds in the active site of the final hEP1 receptor model. The docking enabled us to identify key molecular interactions supported by the mutagenesis data. Also, the correlation of r²=0.81 was observed between the Ki values and the docking scores of 15 prostanoid compounds. The results obtained in this study may provide new insights toward understanding the active site conformation of the hEP1 receptor and can be used for the structure-based design of novel specific ligands.

Reflections on the theory of "silver bullet" octreotide tracers: implications for ligand-receptor interactions in the age of peptides, heterodimers, receptor mosaics, truncated receptors, and multifractal analysis

The classical attitude of Nuclear Medicine practitioners on matters of peptide-receptor interactions has maintained an intrinsic monogamic character since many years. New advances in the field of biochemistry and even in clinical Nuclear Medicine have challenged this type of thinking, which prompted me to work on this review. The central issue of this paper will be the use of somatostatin analogs, i.e., octreotide, in clinical imaging procedures as well as in relation to neuroendocirne tumors. Newly described characteristics of G-protein coupled receptors such as the formation of receptor mosaics will be discussed. A small section will enumerate the regulatory processes found in the cell membrane. Possible new interpretations, other than tumor detection, based on imaging procedures with somatostatin analogs will be presented. The readers will be taken to situations such as inflammation, nociception, mechanosensing, chemosensing, fibrosis, taste, and vascularity where somatostatin is involved. Thyroid-associated orbitopathy will be used as a model for the development of multi-agent therapeutics. The final graphical summary depicts the multifactorial properties of ligand binding.

Molecular dynamics simulations and docking studies on 3D models of the heterodimeric and homodimeric 5-HT2A receptor subtype

Background: G-protein coupled receptors may exist as functional homodimers, heterodimers and even as higher aggregates. In this work, we investigate the 5-HT2A receptor, which is a known target for antipsychotic drugs. Recently, 5-HT2A has been shown to form functional homodimers and heterodimers with the mGluR2 receptor. The objective of this study is to build up 3D models of the 5-HT2A/mGluR2 heterodimer and of the 5-HT2A-5-HT2A homodimer, and to evaluate the impact of the dimerization interface on the shape of the 5-HT2A binding pocket by using molecular dynamics simulations and docking studies. Results and discussion: The heterodimer, homodimer and monomeric 5-HT2A receptors were simulated by molecular dynamics for 40 ns each. The trajectories were clustered and representative structures of six clusters for each system were generated. Inspection of the these representative structures clearly indicate an effect of the dimerization interface on the topology of the binding pocket. Docking studies allowed to generate receiver operating characteristic curves for a set of 5-HT2A ligands, indicating that different complexes prefer different classes of 5-HT2A ligands. Conclusion: This study clearly indicates that the presence of a dimerization interface must explicitly be considered when studying G-protein coupled receptors known to exist as dimers. Molecular dynamics simulation and cluster analysis are appropriate tools to study the phenomenon.

GPCR Conformations: Implications for Rational Drug Design

G protein-coupled receptors (GPCRs) comprise a large class of transmembrane proteins that play critical roles in both normal physiology and pathophysiology. These critical roles offer targets for therapeutic intervention, as exemplified by the substantial fraction of current pharmaceutical agents that target members of this family. Tremendous contributions to our understanding of GPCR structure and dynamics have come from both indirect and direct structural characterization techniques. Key features of GPCR conformations derived from both types of characterization techniques are reviewed.

Structures of the CXCR4 chemokine GPCR with small-molecule and cyclic peptide antagonists

Chemokine receptors are critical regulators of cell migration in the context of immune surveillance, inflammation, and development. The G protein-coupled chemokine receptor CXCR4 is specifically implicated in cancer metastasis and HIV-1 infection. Here we report five independent crystal structures of CXCR4 bound to an antagonist small molecule IT1t and a cyclic peptide CVX15 at 2.5 to 3.2 angstrom resolution. All structures reveal a consistent homodimer with an interface including helices V and VI that may be involved in regulating signaling. The location and shape of the ligand-binding sites differ from other G protein-coupled receptors and are closer to the extracellular surface. These structures provide new clues about the interactions between CXCR4 and its natural ligand CXCL12, and with the HIV-1 glycoprotein gp120.

A lipid pathway for ligand binding is necessary for a cannabinoid G protein-coupled receptor

Recent isothiocyanate covalent labeling studies have suggested that a classical cannabinoid, (-)-7'-isothiocyanato-11-hydroxy-1',1'dimethylheptyl-hexahydrocannabinol (AM841), enters the cannabinoid CB2 receptor via the lipid bilayer (Pei, Y., Mercier, R. W., Anday, J. K., Thakur, G. A., Zvonok, A. M., Hurst, D., Reggio, P. H., Janero, D. R., and Makriyannis, A. (2008) Chem. Biol. 15, 1207-1219). However, the sequence of steps involved in such a lipid pathway entry has not yet been elucidated. Here, we test the hypothesis that the endogenous cannabinoid sn-2-arachidonoylglycerol (2-AG) attains access to the CB2 receptor via the lipid bilayer. To this end, we have employed microsecond time scale all-atom molecular dynamics (MD) simulations of the interaction of 2-AG with CB2 via a palmitoyl-oleoyl-phosphatidylcholine lipid bilayer. Results suggest the following: 1) 2-AG first partitions out of bulk lipid at the transmembrane alpha-helix (TMH) 6/7 interface; 2) 2-AG then enters the CB2 receptor binding pocket by passing between TMH6 and TMH7; 3) the entrance of the 2-AG headgroup into the CB2 binding pocket is sufficient to trigger breaking of the intracellular TMH3/6 ionic lock and the movement of the TMH6 intracellular end away from TMH3; and 4) subsequent to protonation at D3.49/D6.30, further 2-AG entry into the ligand binding pocket results in both a W6.48 toggle switch change and a large influx of water. To our knowledge, this is the first demonstration via unbiased molecular dynamics that a ligand can access the binding pocket of a class A G protein-coupled receptor via the lipid bilayer and the first demonstration via molecular dynamics of G protein-coupled receptor activation triggered by a ligand binding event.

Critical review of general guidelines for membrane proteins model building and analysis

Membrane proteins play major roles in many biological processes such as signalling, transport, etc. They have been shown to be involved in the development of many diseases and have become important drug targets per se. The understanding of their functional properties may be facilitated if a 3D structure is available. However, in the case of membrane proteins, only a few 3D structures have been solved to date. Bioinformatics and molecular modelling approaches are thus powerful alternatives to fill the gap between the sequence and the structure. Here, a review of the most recent approaches is proposed together with guidelines on how to use them. In addition, insofar as important biological processes require conformational changes, we discuss some interesting methods aimed at exploring the dynamic behaviour of proteins in their membrane environment. The paper ends with a brief description of useful approaches for determining oligomerisation or ligand binding sites.

Hetero-oligomerization of CCR2, CCR5, and CXCR4 and the protean effects of "selective" antagonists

Chemokine receptors constitute an attractive family of drug targets in the frame of inflammatory diseases. However, targeting specific chemokine receptors may be complicated by their ability to form dimers or higher order oligomers. Using a combination of luminescence complementation and bioluminescence resonance energy transfer assays, we demonstrate for the first time the existence of hetero-oligomeric complexes composed of at least three chemokine receptors (CCR2, CCR5, and CXCR4). We show in T cells and monocytes that negative binding cooperativity takes place between the binding pockets of these receptors, demonstrating their functional interaction in leukocytes. We also show that specific antagonists of one receptor (TAK-779 or AMD3100) lead to functional cross-inhibition of the others. Finally, using the air pouch model in mice, we show that the CCR2 and CCR5 antagonist TAK-779 inhibits cell recruitment promoted by the CXCR4 agonist SDF-1 alpha, demonstrating that cross-inhibition by antagonists also occurs in vivo. Thus, antagonists of the therapeutically important chemokine receptors regulate the functional properties of other receptors to which they do not bind directly with important implications for the use of these agents in vivo.

Molecular dynamics simulation of the heterodimeric mGluR2/5HT(2A) complex. An atomistic resolution study of a potential new target in psychiatric conditions

Homo- and heterodimerization is becoming an assessed concept in G-protein coupled receptor (GPCR) pharmacology, and the notion that GPCRs may dimerize or oligomerize is allowing for a reinterpretation of some inconsistencies or anomalies and is providing medicinal chemists with potentially relevant novel molecular targets for a variety of therapeutic conditions. Recently, it has been reported that two unrelated GPCRs, namely class C metabotropic glutamate receptor type-2 (mGluR2) and class A 5HT(2A) serotoninergic receptor, can heterodimerize at the transmembrane domain level. We performed a 40 ns molecular dynamics simulation of the mGluR2/5HT(2A) heterocomplex constructed around a TM4/TM5 interface and embedded in an explicit phospholipidic bilayer surrounded by water molecules. In a separate experiment, the monomeric 5HT(2A) receptor was simulated for additional 40 ns under the same conditions. The analysis and the comparison of the two simulations allowed us to clearly identify a cross-talk between the two protomers and to put forward an effect of the heterodimerization on the shape of the binding pocket of 5HT(2A). This result provides the first molecular explanation for the reported allosteric effect of mGluR2 on 5HT(2A)-mediated response and suggests that the heterocomplex can be a more suitable target for in silico screening than the monomeric protomers.

Increasingly accurate dynamic molecular models of G-protein coupled receptor oligomers: Panacea or Pandora's box for novel drug discovery?

For years, conventional drug design at G-protein coupled receptors (GPCRs) has mainly focused on the inhibition of a single receptor at a usually well-defined ligand-binding site. The recent discovery of more and more physiologically relevant GPCR dimers/oligomers suggests that selectively targeting these complexes or designing small molecules that inhibit receptor-receptor interactions might provide new opportunities for novel drug discovery. To uncover the fundamental mechanisms and dynamics governing GPCR dimerization/oligomerization, it is crucial to understand the dynamic process of receptor-receptor association, and to identify regions that are suitable for selective drug binding. This minireview highlights current progress in the development of increasingly accurate dynamic molecular models of GPCR oligomers based on structural, biochemical, and biophysical information that has recently appeared in the literature. In view of this new information, there has never been a more exciting time for computational research into GPCRs than at present. Information-driven modern molecular models of GPCR complexes are expected to efficiently guide the rational design of GPCR oligomer-specific drugs, possibly allowing researchers to reach for the high-hanging fruits in GPCR drug discovery, i.e. more potent and selective drugs for efficient therapeutic interventions.

Structural biology of human cannabinoid receptor-2 helix 6 in membrane-mimetic environments

We detail the structure and dynamics of a synthetic peptide corresponding to transmembrane helix 6 (TMH6) of human cannabinoid receptor-2 (hCB2) in biomembrane-mimetic environments. The peptide's NMR structural biology is characterized by two alpha-helical domains bridged by a flexible, nonhelical hinge region containing a highly-conserved CWFP motif with an environmentally sensitive, Pro-based conformational switch. Buried within the peptide's flexible region, W(258) may hydrogen-bond with L(255) to help stabilize the Pro-kinked hCB2 TMH6 structure and position C(257) advantageously for interaction with agonist ligands. These characteristics of hCB2 TMH6 are potential structural features of ligand-induced hCB2 activation in vivo.

Solid-state NMR and molecular dynamics characterization of cannabinoid receptor-1 (CB1) helix 7 conformational plasticity in model membranes

Little direct information is available regarding the influence of membrane environment on transmembrane (TM) G-protein-coupled receptor (GPCR) conformation and dynamics. The human CB1 cannabinoid receptor (hCB1) is a prominent GPCR pharmacotherapeutic target in which helix 7 appears critical to ligand recognition. We have chemically synthesized a hCB1 peptide corresponding to a segment of TM helix 7 and the entire contiguous helix 8 domain (fourth cytoplasmic loop) and reconstituted it in defined phospholipid-bilayer model membranes. Using an NMR-based strategy combined with molecular dynamics simulations, we provide the first direct experimental description of the orientation of hCB1 helix 7 in phospholipid membranes of varying thickness and the mechanism by which helix-7 conformation adjusts to avoid hydrophobic mismatch. Solid-state (15)N NMR data show that hCB1 helices 7 and 8 reconstituted into phospholipid bilayers are oriented in a TM and in-plane (i.e., parallel to the phospholipid membrane surface) fashion, respectively. TM helix orientation is influenced by the thickness of the hydrophobic membrane bilayer as well as the interaction of helix 8 with phospholipid polar headgroups. Molecular dynamics simulations show that a decrease in phospholipid chain-length induces a kink at P394 in TM helix 7 to avoid hydrophobic mismatch. Thus, the NP(X)nY motif found in hCB1 and highly conserved throughout the GPCR superfamily is important for flexing helix 7 to accommodate bilayer thickness. Dynamic modulation of hCB1-receptor TM helix conformation by its membrane environment may have general relevance to GPCR structure and function.

Structure of the third intracellular loop of the vasopressin V2 receptor and conformational changes upon binding to gC1qR

The V2 vasopressin receptor is a G-protein-coupled receptor that regulates the renal antidiuretic response. Its third intracellular loop is involved in the coupling not only with the GalphaS protein but also with gC1qR, a potential chaperone of G-protein-coupled receptors. In this report, we describe the NMR solution structure of the V2 i3 loop under a cyclized form (i3_cyc) and characterize its interaction with gC1qR. i3_cyc formed a left-twisted alpha-helical hairpin structure. The building of a model of the entire V2 receptor including the i3_cyc NMR structure clarified the side-chain orientation of charged residues, in agreement with literature mutagenesis reports. In the model, the i3 loop formed a rigid helical column, protruding deep inside the cytoplasm, as does the i3 loop in the recently elucidated structure of squid rhodopsin. However, its higher packing angle resulted in a different structural motif at the intracellular interface, which may be important for the specific recognition of GalphaS. Moreover, we could estimate the apparent K(d) of the i3_cyc/gC1qR complex by anisotropy fluorescence. Using a shorter and more soluble version of i3_cyc, which encompassed the putative site of gC1qR binding, we showed by NMR saturation transfer difference spectroscopy that the binding surface corresponded to the central arginine cluster. Binding to gC1qR induced the folding of the otherwise disordered short peptide into a spiral-like path formed by a succession of I and IV turns. Our simulations suggested that this folding would rigidify the arginine cluster in the entire i3 loop and would alter the conformation of the cytosolic extensions of TM V and TM VI helices. In agreement with this conformational rearrangement, we observed that binding of gC1qR to the full-length receptor modifies the intrinsic tryptophan fluorescence binding curves of V2 to an antagonist.

Computational molecular biology approaches to ligand-target interactions

Binding of small molecules to their targets triggers complex pathways. Computational approaches are keys for predictions of the molecular events involved in such cascades. Here we review current efforts at characterizing the molecular determinants in the largest membrane-bound receptor family, the G-protein-coupled receptors (GPCRs). We focus on odorant receptors, which constitute more than half GPCRs. The work presented in this review uncovers structural and energetic aspects of components of the cellular cascade. Finally, a computational approach in the context of radioactive boron-based antitumoral therapies is briefly described.

Ligand-binding architecture of human CB2 cannabinoid receptor: evidence for receptor subtype-specific binding motif and modeling GPCR activation

The extensive physiological influence of transmission through the CB2 cannabinoid receptor makes this G protein-coupled receptor (GPCR) a promising therapeutic target for treating neuropathic pain, inflammation, and immune disorders. However, there is little direct structural information pertaining to either GPCR or CB2-receptor ligand recognition and activation. The present work helps characterize experimentally the ligand-binding interactions of the human CB2 (hCB2) receptor. This study illustrates how our overall experimental approach, "ligand-assisted protein structure" (LAPS), affords direct determination of the requirements for ligand binding to the hCB2 receptor and discrimination among the binding motifs for ligands that activate therapeutically relevant GPCRs.

A novel adrenocorticotropin receptor mutation alters its structure and function, causing familial glucocorticoid deficiency

CONTEXT

Familial glucocorticoid deficiency (FGD) is an autosomal recessive disorder characterized by unresponsiveness to ACTH. In this study, two mutations of the ACTH receptor (MC2R) gene are reported in this FGD clinical case.

OBJECTIVE

The objective of the study was to characterize a novel MC2R gene mutation in a compound heterozygous patient with FGD phenotype.

DESIGN

This was a clinical case description, biochemical, molecular, and bioinformatics analysis to describe a novel MC2R gene mutation.

PATIENTS

The subject of the study was a male diagnosed with primary adrenal insufficiency. The family history showed nonconsanguineous healthy parents, three healthy siblings, and one brother affected with FGD.

MAIN OUTCOME MEASURES

The mutant MC2R-Ala126Ser showed significantly lower activity when it was stimulated with ACTH-(1-24) than did cells transfected with wild-type MC2R.

RESULTS

The molecular studies demonstrated the presence of an adenine heterozygous insertion (InsA1347) in the MC2R gene (G217fs) in the patient. This insertion was due to a frame shift mutation in one allele and a premature stop codon codifying an aberrant receptor of 247 residues (27.2 kDa). We also found a novel heterozygous mutation alanine 126 by serine. Molecular dynamic simulations showed that serine 126 side chain fluctuates forming a noncanonical intrahelical hydrogen bond in the transmembrane helix 3 of the mutated receptor. This produces a structural rearrangement of the MC2R internal cavities that may affect the ligand recognition and signal transduction throughout the G protein.

CONCLUSIONS

We propose a molecular explanation for the reduced activity exhibited by the MC2R alanine 126 by serine mutant.

Conformational memories with variable bond angles

Conformational Memories (CM) is a simulated annealing/Monte Carlo method that explores peptide and protein dihedral conformational space completely and efficiently, independent of the original conformation. Here we extend the CM method to include the variation of a randomly chosen bond angle, in addition to the standard variation of two or three randomly chosen dihedral angles, in each Monte Carlo trial of the CM exploratory and biased phases. We test the hypothesis that the inclusion of variable bond angles in CM leads to an improved sampling of conformational space. We compare the results with variable bond angles to CM with no bond angle variation for the following systems: (1) the pentapeptide Met-enkephalin, which is a standard test case for conformational search methods; (2) the proline ring pucker in a 17mer model peptide, (Ala)(8)Pro(Ala)(8); and (3) the conformations of the Ser 7.39 chi(1) in transmembrane helix 7 (TMH7) of the cannabinoid CB1 receptor, a 25-residue system. In each case, analysis of the CM results shows that the inclusion of variable bond angles results in sampling of regions of conformational space that are inaccessible to CM calculations with only variable dihedral angles, and/or a shift in conformational populations from those calculated when variable bond angles are not included. The incorporation of variable bond angles leads to an improved sampling of conformational space without loss of efficiency. Our examples show that this improved sampling leads to better exploration of biologically relevant conformations that have been experimentally validated.

Computational chemistry approaches to drug discovery in signal transduction

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

Solution structure of urotensin-II receptor extracellular loop III and characterization of its interaction with urotensin-II

Urotensin-II (U-II) is a vasoactive hormone that acts through a G-protein-coupled receptor named UT. Recently, we have shown, using the surface plasmon resonance technology that human U-II (hU-II) interacts with the hUT(281-300) fragment, a segment containing the extracellular loop III (EC-III) and short extensions of the transmembrane domains VI and VII (TM-VI and TM-VII). To further investigate the interaction of UT receptor with U-II, we have determined the solution structure of hUT(281-300) by high-resolution NMR and molecular modeling and we have examined, also using NMR, the binding with hU-II at residue level. In the presence of dodecylphosphocholine micelles, hUT(281-300) exhibited a type III beta-turn (Q285-L288), followed by an -helical structure (A289-L299), the latter including a stretch of transmembrane helix VII. Upon addition of hU-II, significant chemical shift perturbations were observed for residues located just on the N-terminal side of the beta-turn (end of TM-VI/beginning of EC-III) and on one face of the -helix (end of EC-III/beginning of TM-VII). These data, in conjunction with intermolecular NOEs, suggest that the initiation site of EC-III, as well as the upstream portion of helix VII, would be involved in agonist binding and allow to propose points of interaction in the ligand-receptor complex.

G protein-coupled receptor dimers: functional consequences, disease states and drug targets

With an ever-expanding need for reliable therapeutic agents that are highly effective and exhibit minimal deleterious side effects, a greater understanding of the mechanisms underlying G protein-coupled receptor (GPCR) regulation is fundamental. GPCRs comprise more than 30% of all therapeutic drug targets and it is likely that this will only increase as more orphan GPCRs are identified. The past decade has seen a dramatic shift in the prevailing concept of how GPCRs function, in particular the growing acceptance that GPCRs are capable of interacting with one another at a molecular level to form complexes, with significantly different pharmacological properties to their monomeric selves. While the ability of like-receptors to associate and form homodimers raises some interesting mechanistic issues, the possibility that unlike-receptors could heterodimerise in certain tissue types, producing a functionally unique signalling complex that binds specific ligands, provides an invaluable opportunity to refine and redefine pharmacological interventions with greater specificity and efficacy.

Advances in the Development and Application of Computational Methodologies for Structural Modeling of G-Protein Coupled Receptors

BACKGROUND: Despite the large amount of experimental data accumulated in the past decade on G-protein coupled receptor (GPCR) structure and function, understanding of the molecular mechanisms underlying GPCR signaling is still far from being complete, thus impairing the design of effective and selective pharmaceuticals. OBJECTIVE: Understanding of GPCR function has been challenged even further by more recent experimental evidence that several of these receptors are organized in the cell membrane as homo- or hetero-oligomers, and that they may exhibit unique pharmacological properties. Given the complexity of these new signaling systems, researcher's efforts are turning increasingly to molecular modeling, bioinformatics and computational simulations for mechanistic insights of GPCR functional plasticity. METHODS: We review here current advances in the development and application of computational approaches to improve prediction of GPCR structure and dynamics, thus enhancing current understanding of GPCR signaling. RESULTS/CONCLUSIONS: Models resulting from use of these computational approaches further supported by experiments are expected to help elucidate the complex allosterism that propagates through GPCR complexes, ultimately aiming at successful structure-based rational drug design.

Computational modeling study of functional microdomains in cannabinoid receptor type 1

The seven transmembrane helices (TMH) G-protein-coupled receptors (GPCRs) constitute one of the largest superfamily of signaling proteins found in mammals. Some of its members, in which the cannabinoid (CB) receptors are included, stand out because their functional states can be modulated by a broad spectrum of effector molecules. The relative ligand promiscuity exhibited by these receptors could be related with particular attributes conferred by their molecular architecture and represents a motivating issue to be explored. In this regard, this study represents an effort to investigate the cannabinoid receptor type 1 (CB1) ligand recognition plasticity, using comparative modeling, molecular dynamics (MD) simulations and docking. Our results suggest that a cooperative set of subtle structural rearrangements within the TMHs provide to the CB1 protein the plasticity to reach alternate configurations. These changes include the relaxation of intramolecular constraints, the rotations, translations and kinks of the majority of TMHs and the reorganization of the ligand binding cavities.

Solid-state NMR spectroscopy as a tool for drug design: from membrane-embedded targets to amyloid fibrils

Structure-based design has gained credibility as a valuable component of the modern drug discovery process. The technique of SSNMR (solid-state NMR) promises to be a useful counterpart to the conventional experimental techniques of X-ray crystallography and solution-state NMR for providing structural features of drug targets that can guide medicinal chemistry towards drug candidates. This article highlights some recent SSNMR approaches from our group for identifying active compounds, such as enzyme inhibitors, receptor antagonists and peptide agents, that prevent the aggregation of amyloid proteins involved in neurodegenerative diseases. It is anticipated that the use of SSNMR in drug discovery will become more widespread in the wake of advances in hardware and methodological developments.

Analysis of the activation mechanism of the guinea-pig Histamine H1-receptor

The Histamine H(1)-receptor (H1R), belonging to the amine receptor-class of family A of the G-protein coupled receptors (GPCRs) gets activated by agonists. The consequence is a conformational change of the receptor, which may involve the binding-pocket. So, for a good prediction of the binding-mode of an agonist, it is necessary to have knowledge about these conformational changes. Meanwhile some experimental data about the structural changes of GPCRs during activation exist. Based on homology modeling of the guinea-pig H1R (gpH1R), using the crystal structure of bovine rhodopsin as template, we performed several MD simulations with distance restraints in order to get an inactive and an active structure of the gpH1R. The calculations led to a Phe6.44/Trp6.48/Phe6.52-switch and linearization of the proline kinked transmembrane helix VI during receptor activation. Our calculations showed that the Trp6.48/Phe6.52-switch induces a conformational change in Phe6.44, which slides between transmembrane helices III and VI. Additionally we observed a hydrogen bond interaction of Ser3.39 with Asn7.45 in the inactive gpH1R, but because of a counterclockwise rotation of transmembrane helix III Ser3.39 establishes a water-mediated hydrogen bond to Asp2.50 in the active gpH1R. Additionally we simulated a possible mechanism for receptor activation with a modified LigPath-algorithm.

Allosteric properties of G protein-coupled receptor oligomers

Allosteric regulation of ligand binding is a well-established mechanism regulating the function of G protein-coupled receptors (GPCR). Allosteric modulators have been considered so far as molecules binding to an allosteric site, distinct from that of the reference ligand (orthosteric site), and able to modulate the binding affinity at the orthosteric site and/or the signaling properties resulting from orthosteric site occupancy. Given that most GPCR are known to form dimers or higher order oligomers, we explored whether allosteric interactions could also occur between protomers within oligomeric arrays, thereby influencing binding and signaling receptor properties. Two main conclusions emerged from such studies. First, allosteric modulators can affect one receptor by binding to another receptor within a dimeric or oligomeric complex. Second, allosteric modulators might act on a given receptor by targeting the "orthosteric site" in another receptor of the complex. Allosteric regulation within di(oligo)mers thus implies that the pharmacological properties of a given receptor subtype can be influenced by the array of dimerization partners coexpressed in each particular cell type. Ligands could thus act as agonists or antagonists on 1 receptor, while modulating allosterically the function of a variety of other receptors to which they do not bind directly. Allosteric regulation across GPCR oligomeric interfaces is expected to greatly influence the practice of pharmacology. It will likely affect the design of drug discovery programs, which rely mostly on the overexpression of the receptor of interest in a cell line, thereby focusing on homo-oligomers and ignoring the potential effects of other partners.

Requirements and ontology for a G protein-coupled receptor oligomerization knowledge base

BACKGROUND

G Protein-Coupled Receptors (GPCRs) are a large and diverse family of membrane proteins whose members participate in the regulation of most cellular and physiological processes and therefore represent key pharmacological targets. Although several bioinformatics resources support research on GPCRs, most of these have been designed based on the traditional assumption that monomeric GPCRs constitute the functional receptor unit. The increase in the frequency and number of reports about GPCR dimerization/oligomerization and the implication of oligomerization in receptor function makes necessary the ability to store and access information about GPCR dimers/oligomers electronically.

RESULTS

We present here the requirements and ontology (the information scheme to describe oligomers and associated concepts and their relationships) for an information system that can manage the elements of information needed to describe comprehensively the phenomena of both homo- and hetero-oligomerization of GPCRs. The comprehensive information management scheme that we plan to use for the development of an intuitive and user-friendly GPCR-Oligomerization Knowledge Base (GPCR-OKB) is the result of a community dialog involving experimental and computational colleagues working on GPCRs.

CONCLUSION

Our long term goal is to disseminate to the scientific community organized, curated, and detailed information about GPCR dimerization/oligomerization and its related structural context. This information will be reported as close to the data as possible so the user can make his own judgment on the conclusions drawn for a particular study. The requirements and ontology described here will facilitate the development of future information systems for GPCR oligomers that contain both computational and experimental information about GPCR oligomerization. This information is freely accessible at http://www.gpcr-okb.org.