1
|
Stampelou M, Ladds G, Kolocouris A. Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A 3 Receptor. J Phys Chem B 2024; 128:914-936. [PMID: 38236582 DOI: 10.1021/acs.jpcb.3c05986] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.
Collapse
Affiliation(s)
- Margarita Stampelou
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Graham Ladds
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K
| | - Antonios Kolocouris
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| |
Collapse
|
2
|
Jabeen A, Vijayram R, Ranganathan S. BIO-GATS: A Tool for Automated GPCR Template Selection Through a Biophysical Approach for Homology Modeling. Front Mol Biosci 2021; 8:617176. [PMID: 33898512 PMCID: PMC8059640 DOI: 10.3389/fmolb.2021.617176] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Accepted: 02/24/2021] [Indexed: 11/13/2022] Open
Abstract
G protein-coupled receptors (GPCRs) are the largest family of membrane proteins with more than 800 members. GPCRs are involved in numerous physiological functions within the human body and are the target of more than 30% of the United States Food and Drug Administration (FDA) approved drugs. At present, over 400 experimental GPCR structures are available in the Protein Data Bank (PDB) representing 76 unique receptors. The absence of an experimental structure for the majority of GPCRs demand homology models for structure-based drug discovery workflows. The generation of good homology models requires appropriate templates. The commonly used methods for template selection are based on sequence identity. However, there exists low sequence identity among the GPCRs. Sequences with similar patterns of hydrophobic residues are often structural homologs, even with low sequence identity. Extending this, we propose a biophysical approach for template selection based principally on hydrophobicity correspondence between the target and the template. Our approach takes into consideration other relevant parameters, including resolution, similarity within the orthosteric binding pocket of GPCRs, and structure completeness, for template selection. The proposed method was implemented in the form of a free tool called Bio-GATS, to provide the user with easy selection of the appropriate template for a query GPCR sequence. Bio-GATS was successfully validated with recent published benchmarking datasets. An application to an olfactory receptor to select an appropriate template has also been provided as a case study.
Collapse
Affiliation(s)
- Amara Jabeen
- Department of Molecular Sciences, Macquarie University, Sydney, NSW, Australia
| | - Ramya Vijayram
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, India
| | - Shoba Ranganathan
- Department of Molecular Sciences, Macquarie University, Sydney, NSW, Australia
| |
Collapse
|
3
|
Jacobson KA, Gao Z, Matricon P, Eddy MT, Carlsson J. Adenosine A2Areceptor antagonists: from caffeine to selective non‐xanthines. Br J Pharmacol 2020; 179:3496-3511. [DOI: 10.1111/bph.15103] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2020] [Revised: 04/30/2020] [Accepted: 05/03/2020] [Indexed: 12/12/2022] Open
Affiliation(s)
- Kenneth A. Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic ChemistryNational Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health Bethesda MD USA
| | - Zhan‐Guo Gao
- Molecular Recognition Section, Laboratory of Bioorganic ChemistryNational Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health Bethesda MD USA
| | - Pierre Matricon
- Department of Cell and Molecular Biology, Science for Life LaboratoryUppsala University Uppsala Sweden
| | - Matthew T. Eddy
- Department of ChemistryUniversity of Florida Gainesville FL USA
| | - Jens Carlsson
- Department of Cell and Molecular Biology, Science for Life LaboratoryUppsala University Uppsala Sweden
| |
Collapse
|
4
|
Salmaso V, Jacobson KA. In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors. Biomolecules 2020; 10:E812. [PMID: 32466404 PMCID: PMC7356333 DOI: 10.3390/biom10060812] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 05/21/2020] [Accepted: 05/22/2020] [Indexed: 12/14/2022] Open
Abstract
Molecular modeling has contributed to drug discovery for purinergic GPCRs, including adenosine receptors (ARs) and P2Y receptors (P2YRs). Experimental structures and homology modeling have proven to be useful in understanding and predicting structure activity relationships (SAR) of agonists and antagonists. This review provides an excursus on molecular dynamics (MD) simulations applied to ARs and P2YRs. The binding modes of newly synthesized A1AR- and A3AR-selective nucleoside derivatives, potentially of use against depression and inflammation, respectively, have been predicted to recapitulate their SAR and the species dependence of A3AR affinity. P2Y12R and P2Y1R crystallographic structures, respectively, have provided a detailed understanding of the recognition of anti-inflammatory P2Y14R antagonists and a large group of allosteric and orthosteric antagonists of P2Y1R, an antithrombotic and neuroprotective target. MD of A2AAR (an anticancer and neuroprotective target), A3AR, and P2Y1R has identified microswitches that are putatively involved in receptor activation. The approach pathways of different ligands toward A2AAR and P2Y1R binding sites have also been explored. A1AR, A2AAR, and A3AR were utilizes to study allosteric phenomena, but locating the binding site of structurally diverse allosteric modulators, such as an A3AR enhancer LUF6000, is challenging. Ligand residence time, a predictor of in vivo efficacy, and the structural role of water were investigated through A2AAR MD simulations. Thus, new MD and other modeling algorithms have contributed to purinergic GPCR drug discovery.
Collapse
Affiliation(s)
| | - Kenneth A. Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA;
| |
Collapse
|
5
|
Orgován Z, Ferenczy GG, Keserű GM. The role of water and protein flexibility in the structure-based virtual screening of allosteric GPCR modulators: an mGlu 5 receptor case study. J Comput Aided Mol Des 2019; 33:787-797. [PMID: 31542869 PMCID: PMC6825653 DOI: 10.1007/s10822-019-00224-w] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2019] [Accepted: 09/15/2019] [Indexed: 12/22/2022]
Abstract
Stabilizing unique receptor conformations, allosteric modulators of G-protein coupled receptors (GPCRs) might open novel treatment options due to their new pharmacological action, their enhanced specificity and selectivity in both binding and signaling. Ligand binding occurs at intrahelical allosteric sites and involves significant induced fit effects that include conformational changes in the local protein environment and water networks. Based on the analysis of available crystal structures of metabotropic glutamate receptor 5 (mGlu5) we investigated these effects in the binding of mGlu5 receptor negative allosteric modulators. A large set of retrospective virtual screens revealed that the use of multiple protein structures and the inclusion of selected water molecules improves virtual screening performance compared to conventional docking strategies. The role of water molecules and protein flexibility in ligand binding can be taken into account efficiently by the proposed docking protocol that provided reasonable enrichment of true positives. This protocol is expected to be useful also for identifying intrahelical allosteric modulators for other GPCR targets.
Collapse
Affiliation(s)
- Zoltán Orgován
- Medicinal Chemistry Research Group, Research Centre for Natural Sciences, Hungarian Academy of Sciences, Magyar tudósok krt. 2, Budapest, 1117, Hungary
| | - György G Ferenczy
- Medicinal Chemistry Research Group, Research Centre for Natural Sciences, Hungarian Academy of Sciences, Magyar tudósok krt. 2, Budapest, 1117, Hungary
| | - György M Keserű
- Medicinal Chemistry Research Group, Research Centre for Natural Sciences, Hungarian Academy of Sciences, Magyar tudósok krt. 2, Budapest, 1117, Hungary.
| |
Collapse
|
6
|
Deb PK, Chandrasekaran B, Mailavaram R, Tekade RK, Jaber AMY. Molecular modeling approaches for the discovery of adenosine A2B receptor antagonists: current status and future perspectives. Drug Discov Today 2019; 24:1854-1864. [DOI: 10.1016/j.drudis.2019.05.011] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2019] [Revised: 03/26/2019] [Accepted: 05/10/2019] [Indexed: 12/13/2022]
|
7
|
Deb PK. Recent Updates in the Computer Aided Drug Design Strategies for the Discovery of Agonists and Antagonists of Adenosine Receptors. Curr Pharm Des 2019; 25:747-749. [DOI: 10.2174/1381612825999190515120510] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Pran Kishore Deb
- Faculty of Pharmacy, Philadelphia University, PO Box 1, 19392, Amman, Jordan
| |
Collapse
|
8
|
Congreve M, Brown GA, Borodovsky A, Lamb ML. Targeting adenosine A2A receptor antagonism for treatment of cancer. Expert Opin Drug Discov 2018; 13:997-1003. [DOI: 10.1080/17460441.2018.1534825] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Miles Congreve
- Heptares Therapeutics Limited, Steinmetz Building, Cambridge, Granta Park, UK
| | - Giles A. Brown
- Heptares Therapeutics Limited, Steinmetz Building, Cambridge, Granta Park, UK
| | | | - Michelle L. Lamb
- Medicinal Chemistry, Oncology, IMED Biotech Unit, AstraZeneca, Boston, MA, USA
| |
Collapse
|
9
|
Breakthrough in GPCR Crystallography and Its Impact on Computer-Aided Drug Design. Methods Mol Biol 2018; 1705:45-72. [PMID: 29188558 DOI: 10.1007/978-1-4939-7465-8_3] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Recent crystallographic structures of G protein-coupled receptors (GPCRs) have greatly advanced our understanding of the recognition of their diverse agonist and antagonist ligands. We illustrate here how this applies to A2A adenosine receptors (ARs) and to P2Y1 and P2Y12 receptors (P2YRs) for ADP. These X-ray structures have impacted the medicinal chemistry aimed at discovering new ligands for these two receptor families, including receptors that have not yet been crystallized but are closely related to the known structures. In this Chapter, we discuss recent structure-based drug design projects that led to the discovery of: (a) novel A3AR agonists based on a highly rigidified (N)-methanocarba scaffold for the treatment of chronic neuropathic pain and other conditions, (b) fluorescent probes of the ARs and P2Y14R, as chemical tools for structural probing of these GPCRs and for improving assay capabilities, and (c) new more drug-like antagonists of the inflammation-related P2Y14R. We also describe the computationally enabled molecular recognition of positive (for A3AR) and negative (P2Y1R) allosteric modulators that in some cases are shown to be consistent with structure-activity relationship (SAR) data. Thus, computational modeling has become an essential tool for the design of purine receptor ligands.
Collapse
|
10
|
Lagarias P, Vrontaki E, Lambrinidis G, Stamatis D, Convertino M, Ortore G, Mavromoustakos T, Klotz KN, Kolocouris A. Discovery of Novel Adenosine Receptor Antagonists through a Combined Structure- and Ligand-Based Approach Followed by Molecular Dynamics Investigation of Ligand Binding Mode. J Chem Inf Model 2018; 58:794-815. [PMID: 29485875 DOI: 10.1021/acs.jcim.7b00455] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
An intense effort is made by pharmaceutical and academic research laboratories to identify and develop selective antagonists for each adenosine receptor (AR) subtype as potential clinical candidates for "soft" treatment of various diseases. Crystal structures of subtypes A2A and A1ARs offer exciting opportunities for structure-based drug design. In the first part of the present work, Maybridge HitFinder library of 14400 compounds was utilized to apply a combination of structure-based against the crystal structure of A2AAR and ligand-based methodologies. The docking poses were rescored by CHARMM energy minimization and calculation of the desolvation energy using Poisson-Boltzmann equation electrostatics. Out of the eight selected and tested compounds, five were found positive hits (63% success). Although the project was initially focused on targeting A2AAR, the identified antagonists exhibited low micromolar or micromolar affinity against A2A/A3, ARs, or A3AR, respectively. Based on these results, 19 compounds characterized by novel chemotypes were purchased and tested. Sixteen of them were identified as AR antagonists with affinity toward combinations of the AR family isoforms (A2A/A3, A1/A3, A1/A2A/A3, and A3). The second part of this work involves the performance of hundreds of molecular dynamics (MD) simulations of complexes between the ARs and a total of 27 ligands to resolve the binding interactions of the active compounds, which were not achieved by docking calculations alone. This computational work allowed the prediction of stable and unstable complexes which agree with the experimental results of potent and inactive compounds, respectively. Of particular interest is that the 2-amino-thiophene-3-carboxamides, 3-acylamino-5-aryl-thiophene-2-carboxamides, and carbonyloxycarboximidamide derivatives were found to be selective and possess a micromolar to low micromolar affinity for the A3 receptor.
Collapse
Affiliation(s)
- Panagiotis Lagarias
- Division of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences , National and Kapodistrian University of Athens , Panepistimiopolis-Zografou , 15771 Athens , Greece
| | - Eleni Vrontaki
- Division of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences , National and Kapodistrian University of Athens , Panepistimiopolis-Zografou , 15771 Athens , Greece
| | - George Lambrinidis
- Division of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences , National and Kapodistrian University of Athens , Panepistimiopolis-Zografou , 15771 Athens , Greece
| | - Dimitrios Stamatis
- Division of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences , National and Kapodistrian University of Athens , Panepistimiopolis-Zografou , 15771 Athens , Greece
| | - Marino Convertino
- Department of Biochemistry & Biophysics , University of North Carolina at Chapel Hill , 120 Mason Farm Road , Chapel Hill , North Carolina 27599 , United States
| | - Gabriella Ortore
- Department of Pharmacy , University of Pisa , 56126 Pisa , Italy
| | - Thomas Mavromoustakos
- Division of Organic Chemistry, Department of Chemistry, School of Science , National and Kapodistrian University of Athens , Panepistimiopolis-Zografou , 15771 Athens , Greece
| | - Karl-Norbert Klotz
- Institute of Pharmacology and Toxicology , University of Würzburg Versbacher Str. 9 , 97078 Würzburg , Germany
| | - Antonios Kolocouris
- Division of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences , National and Kapodistrian University of Athens , Panepistimiopolis-Zografou , 15771 Athens , Greece
| |
Collapse
|
11
|
Badshah SL, Ullah A, Al-Showiman SS, Mabkhot YN. Better agonist for the opioid receptors. Chem Cent J 2018; 12:13. [PMID: 29417377 PMCID: PMC5803155 DOI: 10.1186/s13065-018-0383-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2016] [Accepted: 01/30/2018] [Indexed: 12/29/2022] Open
Abstract
This commentary highlights the recent work published in journal Nature on the structural based discovery of novel analgesic compounds for opioid receptors with minimal effects. Manglik et al. selectively targeted the Gi based μOR pathway instead of the β-arrestin pathway of the opioids. The computational screening of millions of compounds showed a list of several competent ligands. From these ligands they synthesized the compounds with the best docking score, which were further optimized by adding side residues for better interaction with the μOR. A promising compound, PZM21, was a selective agonist of μOR. It has better analgesic properties with minimal side effects of respiratory depression and constipation. This work is a step towards better drug designing and synthesizing in terms of efficacy, specificity with least side effects of targeted GPCR proteins present in the human proteome.
Collapse
Affiliation(s)
- Syed Lal Badshah
- Department of Chemistry, Islamia College University, Peshawar, 25120, Pakistan.
| | - Asad Ullah
- Department of Chemistry, Islamia College University, Peshawar, 25120, Pakistan
| | - Salim S Al-Showiman
- Department of Chemistry, College of Sciences, King Saud University, Riyadh, Saudi Arabia
| | - Yahia Nasser Mabkhot
- Department of Chemistry, College of Sciences, King Saud University, Riyadh, Saudi Arabia.
| |
Collapse
|
12
|
Lemos A, Melo R, Preto AJ, Almeida JG, Moreira IS, Cordeiro MNDS. In Silico Studies Targeting G-protein Coupled Receptors for Drug Research Against Parkinson's Disease. Curr Neuropharmacol 2018; 16:786-848. [PMID: 29521236 PMCID: PMC6080095 DOI: 10.2174/1570159x16666180308161642] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2017] [Revised: 02/16/2018] [Accepted: 02/02/2018] [Indexed: 11/22/2022] Open
Abstract
Parkinson's Disease (PD) is a long-term neurodegenerative brain disorder that mainly affects the motor system. The causes are still unknown, and even though currently there is no cure, several therapeutic options are available to manage its symptoms. The development of novel antiparkinsonian agents and an understanding of their proper and optimal use are, indeed, highly demanding. For the last decades, L-3,4-DihydrOxyPhenylAlanine or levodopa (L-DOPA) has been the gold-standard therapy for the symptomatic treatment of motor dysfunctions associated to PD. However, the development of dyskinesias and motor fluctuations (wearing-off and on-off phenomena) associated with long-term L-DOPA replacement therapy have limited its antiparkinsonian efficacy. The investigation for non-dopaminergic therapies has been largely explored as an attempt to counteract the motor side effects associated with dopamine replacement therapy. Being one of the largest cell membrane protein families, G-Protein-Coupled Receptors (GPCRs) have become a relevant target for drug discovery focused on a wide range of therapeutic areas, including Central Nervous System (CNS) diseases. The modulation of specific GPCRs potentially implicated in PD, excluding dopamine receptors, may provide promising non-dopaminergic therapeutic alternatives for symptomatic treatment of PD. In this review, we focused on the impact of specific GPCR subclasses, including dopamine receptors, adenosine receptors, muscarinic acetylcholine receptors, metabotropic glutamate receptors, and 5-hydroxytryptamine receptors, on the pathophysiology of PD and the importance of structure- and ligand-based in silico approaches for the development of small molecules to target these receptors.
Collapse
Affiliation(s)
- Agostinho Lemos
- LAQV/REQUIMTE, Department of Chemistry and Biochemistry, Faculty of Sciences, University of Porto, Rua do Campo Alegre s/n, 4169-007Porto, Portugal
- GIGA Cyclotron Research Centre In Vivo Imaging, University of Liège, 4000Liège, Belgium
| | - Rita Melo
- CNC - Center for Neuroscience and Cell Biology, Faculty of Medicine, University of Coimbra, Rua Larga, 3004-517Coimbra, Portugal
- Centro de Ciências e Tecnologias Nucleares, Instituto Superior Técnico, Universidade de Lisboa, Estrada Nacional 10 (ao km 139,7), 2695-066 Bobadela LRS, Portugal
| | - Antonio Jose Preto
- CNC - Center for Neuroscience and Cell Biology, Faculty of Medicine, University of Coimbra, Rua Larga, 3004-517Coimbra, Portugal
| | - Jose Guilherme Almeida
- CNC - Center for Neuroscience and Cell Biology, Faculty of Medicine, University of Coimbra, Rua Larga, 3004-517Coimbra, Portugal
| | - Irina Sousa Moreira
- CNC - Center for Neuroscience and Cell Biology, Faculty of Medicine, University of Coimbra, Rua Larga, 3004-517Coimbra, Portugal
- Bijvoet Center for Biomolecular Research, Faculty of Science - Chemistry, Utrecht University, Utrecht, 3584CH, The Netherlands
| | - Maria Natalia Dias Soeiro Cordeiro
- LAQV/REQUIMTE, Department of Chemistry and Biochemistry, Faculty of Sciences, University of Porto, Rua do Campo Alegre s/n, 4169-007Porto, Portugal
| |
Collapse
|
13
|
Loo JSE, Emtage AL, Ng KW, Yong ASJ, Doughty SW. Assessing GPCR homology models constructed from templates of various transmembrane sequence identities: Binding mode prediction and docking enrichment. J Mol Graph Model 2017; 80:38-47. [PMID: 29306746 DOI: 10.1016/j.jmgm.2017.12.017] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2017] [Revised: 11/27/2017] [Accepted: 12/26/2017] [Indexed: 11/15/2022]
Abstract
GPCR crystal structures have become more readily accessible in recent years. However, homology models of GPCRs continue to play an important role as many GPCR structures remain unsolved. The new crystal structures now available provide not only additional templates for homology modelling but also the opportunity to assess the performance of homology models against their respective crystal structures and gain insight into the performance of such models. In this study we have constructed homology models from templates of various transmembrane sequence identities for eight GPCR targets to better understand the relationship between transmembrane sequence identity and model quality. Model quality was assessed relative to the crystal structure in terms of structural accuracy as well as performance in two typical structure-based drug design applications: ligand binding pose prediction and docking enrichment in virtual screening. Crystal structures significantly outperformed homology models in both assessments. Accurate ligand binding pose prediction was possible but difficult to achieve using homology models, even with the use of induced fit docking. In virtual screening using homology models still conferred significant enrichment compared to random selection, with a clear benefit also observed in using models optimized through induced fit docking. Our results indicate that while homology models that are reasonably accurate structurally can be constructed, without significant refinement homology models will be outperformed by crystal structures in ligand binding pose prediction and docking enrichment regardless of the template used, primarily due to the extremely high level of structural accuracy needed for such applications.
Collapse
Affiliation(s)
- Jason S E Loo
- School of Pharmacy, Taylor's University, No.1 Jalan Taylor's, 47500 Subang Jaya, Selangor, Malaysia.
| | - Abigail L Emtage
- School of Pharmacy, The University of Nottingham Malaysia Campus, Jalan Broga, 43500 Semenyih, Selangor, Malaysia
| | - Kar Weng Ng
- School of Pharmacy, Taylor's University, No.1 Jalan Taylor's, 47500 Subang Jaya, Selangor, Malaysia
| | - Alene S J Yong
- School of Pharmacy, Taylor's University, No.1 Jalan Taylor's, 47500 Subang Jaya, Selangor, Malaysia
| | - Stephen W Doughty
- Penang Medical College, No. 4 Jalan Sepoy Lines, 10450 George Town, Penang, Malaysia
| |
Collapse
|
14
|
Computer-Aided Drug Design Approaches to Study Key Therapeutic Targets in Alzheimer’s Disease. ACTA ACUST UNITED AC 2017. [DOI: 10.1007/978-1-4939-7404-7_3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
15
|
Yadav R, Deepshikha D, Srivastava P. Homology Modeling and Protein Interaction Map of CHRNA7 Neurogenesis Protein. Ann Neurosci 2017; 24:173-179. [PMID: 28867899 DOI: 10.1159/000477155] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Accepted: 04/10/2017] [Indexed: 11/19/2022] Open
Abstract
CHRNA7 is a neurodevelopmental protein involved in differentiation and neurogenesis, which is also named as nicotinic acetylcholine receptors, cholinergic receptor, nicotinic, alpha 7 (neuronal). The protein encoded by this gene forms a homo-oligomeric channel. It is a major component of brain nicotinic receptors displays that are blocked by and sensitive to alpha-bungarotoxin. Studies reports involvement of CHRNA7 protein in different neurological diseases. Non-availability of 3-dimensional (3D) structure leads the study toward structure 3D prediction along with its interaction analysis. The current paper is focused on the structure prediction through homology modeling of CHRNA7 along with binding site prediction using Schrödinger software suite. In continuation of the study, protein-protein interaction analysis is carried out by using string database. Tertiary structure along with binding sites was obtained, and visualized CHRAN7 protein have interaction with CHRNA protein family along with JAK2, AKT1, PICK1 protein that are involved in neurological disease. Structure formation analysis is an important aspect of proteomics studies. Hence, this predicted structure can be used for further advance studies and drug designing. Protein interaction analysis shows that CHRNA7 protein also interact with AKT1 protein which regulate neuronal differentiation and development, that signifies the role of CHRNA7 protein in neurological diseases.
Collapse
Affiliation(s)
- Ruchi Yadav
- AMITY Institute of Biotechnology, AMITY University, Lucknow, India
| | | | | |
Collapse
|
16
|
In search of novel ligands using a structure-based approach: a case study on the adenosine A 2A receptor. J Comput Aided Mol Des 2016; 30:863-874. [PMID: 27629350 PMCID: PMC5083784 DOI: 10.1007/s10822-016-9963-7] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2016] [Accepted: 09/02/2016] [Indexed: 12/03/2022]
Abstract
In this work, we present a case study to explore the challenges associated with finding novel molecules for a receptor that has been studied in depth and has a wealth of chemical information available. Specifically, we apply a previously described protocol that incorporates explicit water molecules in the ligand binding site to prospectively screen over 2.5 million drug-like and lead-like compounds from the commercially available eMolecules database in search of novel binders to the adenosine A2A receptor (A2AAR). A total of seventy-one compounds were selected for purchase and biochemical assaying based on high ligand efficiency and high novelty (Tanimoto coefficient ≤0.25 to any A2AAR tested compound). These molecules were then tested for their affinity to the adenosine A2A receptor in a radioligand binding assay. We identified two hits that fulfilled the criterion of ~50 % radioligand displacement at a concentration of 10 μM. Next we selected an additional eight novel molecules that were predicted to make a bidentate interaction with Asn2536.55, a key interacting residue in the binding pocket of the A2AAR. None of these eight molecules were found to be active. Based on these results we discuss the advantages of structure-based methods and the challenges associated with finding chemically novel molecules for well-explored targets.
Collapse
|
17
|
Guo D, Heitman LH, IJzerman AP. Kinetic Aspects of the Interaction between Ligand and G Protein-Coupled Receptor: The Case of the Adenosine Receptors. Chem Rev 2016; 117:38-66. [DOI: 10.1021/acs.chemrev.6b00025] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Dong Guo
- Division of Medicinal Chemistry,
Leiden Academic Centre for Drug Research (LACDR), Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Laura H. Heitman
- Division of Medicinal Chemistry,
Leiden Academic Centre for Drug Research (LACDR), Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Adriaan P. IJzerman
- Division of Medicinal Chemistry,
Leiden Academic Centre for Drug Research (LACDR), Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
18
|
Chen JF, Lee CF, Chern Y. Adenosine receptor neurobiology: overview. INTERNATIONAL REVIEW OF NEUROBIOLOGY 2015; 119:1-49. [PMID: 25175959 DOI: 10.1016/b978-0-12-801022-8.00001-5] [Citation(s) in RCA: 90] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Adenosine is a naturally occurring nucleoside that is distributed ubiquitously throughout the body as a metabolic intermediary. In the brain, adenosine functions as an important upstream neuromodulator of a broad spectrum of neurotransmitters, receptors, and signaling pathways. By acting through four G-protein-coupled receptors, adenosine contributes critically to homeostasis and neuromodulatory control of a variety of normal and abnormal brain functions, ranging from synaptic plasticity, to cognition, to sleep, to motor activity to neuroinflammation, and cell death. This review begun with an overview of the gene and genome structure and the expression pattern of adenosine receptors (ARs). We feature several new developments over the past decade in our understanding of AR functions in the brain, with special focus on the identification and characterization of canonical and noncanonical signaling pathways of ARs. We provide an update on functional insights from complementary genetic-knockout and pharmacological studies on the AR control of various brain functions. We also highlight several novel and recent developments of AR neurobiology, including (i) recent breakthrough in high resolution of three-dimension structure of adenosine A2A receptors (A2ARs) in several functional status, (ii) receptor-receptor heterodimerization, (iii) AR function in glial cells, and (iv) the druggability of AR. We concluded the review with the contention that these new developments extend and strengthen the support for A1 and A2ARs in brain as therapeutic targets for neurologic and psychiatric diseases.
Collapse
Affiliation(s)
- Jiang-Fan Chen
- Department of Neurology, Boston University School of Medicine, Boston, Massachusetts, USA.
| | - Chien-fei Lee
- Division of Neuroscience, Institute of Biomedical Sciences, Academia Sinica, Taipei, Taiwan
| | - Yijuang Chern
- Division of Neuroscience, Institute of Biomedical Sciences, Academia Sinica, Taipei, Taiwan
| |
Collapse
|
19
|
Duan X, Zhang M, Zhang X, Wang F, Lei M. Molecular modeling and docking study on dopamine D2-like and serotonin 5-HT2A receptors. J Mol Graph Model 2015; 57:143-55. [DOI: 10.1016/j.jmgm.2015.01.014] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2014] [Revised: 01/23/2015] [Accepted: 01/29/2015] [Indexed: 01/22/2023]
|
20
|
Beuming T, Lenselink B, Pala D, McRobb F, Repasky M, Sherman W. Docking and Virtual Screening Strategies for GPCR Drug Discovery. Methods Mol Biol 2015; 1335:251-76. [PMID: 26260606 DOI: 10.1007/978-1-4939-2914-6_17] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
Progress in structure determination of G protein-coupled receptors (GPCRs) has made it possible to apply structure-based drug design (SBDD) methods to this pharmaceutically important target class. The quality of GPCR structures available for SBDD projects fall on a spectrum ranging from high resolution crystal structures (<2 Å), where all water molecules in the binding pocket are resolved, to lower resolution (>3 Å) where some protein residues are not resolved, and finally to homology models that are built using distantly related templates. Each GPCR project involves a distinct set of opportunities and challenges, and requires different approaches to model the interaction between the receptor and the ligands. In this review we will discuss docking and virtual screening to GPCRs, and highlight several refinement and post-processing steps that can be used to improve the accuracy of these calculations. Several examples are discussed that illustrate specific steps that can be taken to improve upon the docking and virtual screening accuracy. While GPCRs are a unique target class, many of the methods and strategies outlined in this review are general and therefore applicable to other protein families.
Collapse
Affiliation(s)
- Thijs Beuming
- Schrödinger, Inc., 120 West 45th Street, New York, NY, 10036, USA,
| | | | | | | | | | | |
Collapse
|
21
|
Yuan G, Gedeon NG, Jankins TC, Jones GB. Novel approaches for targeting the adenosine A2Areceptor. Expert Opin Drug Discov 2014; 10:63-80. [DOI: 10.1517/17460441.2015.971006] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
22
|
Congreve M, Dias JM, Marshall FH. Structure-based drug design for G protein-coupled receptors. PROGRESS IN MEDICINAL CHEMISTRY 2014; 53:1-63. [PMID: 24418607 DOI: 10.1016/b978-0-444-63380-4.00001-9] [Citation(s) in RCA: 57] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Our understanding of the structural biology of G protein-coupled receptors has undergone a transformation over the past 5 years. New protein-ligand complexes are described almost monthly in high profile journals. Appreciation of how small molecules and natural ligands bind to their receptors has the potential to impact enormously how medicinal chemists approach this major class of receptor targets. An outline of the key topics in this field and some recent examples of structure- and fragment-based drug design are described. A table is presented with example views of each G protein-coupled receptor for which there is a published X-ray structure, including interactions with small molecule antagonists, partial and full agonists. The possible implications of these new data for drug design are discussed.
Collapse
Affiliation(s)
- Miles Congreve
- Heptares Therapeutics Ltd, BioPark, Welwyn Garden City, Hertfordshire, United Kingdom
| | - João M Dias
- Heptares Therapeutics Ltd, BioPark, Welwyn Garden City, Hertfordshire, United Kingdom
| | - Fiona H Marshall
- Heptares Therapeutics Ltd, BioPark, Welwyn Garden City, Hertfordshire, United Kingdom
| |
Collapse
|
23
|
Lenselink EB, Beuming T, Sherman W, van Vlijmen HWT, IJzerman AP. Selecting an optimal number of binding site waters to improve virtual screening enrichments against the adenosine A2A receptor. J Chem Inf Model 2014; 54:1737-46. [PMID: 24835542 DOI: 10.1021/ci5000455] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
A major challenge in structure-based virtual screening (VS) involves the treatment of explicit water molecules during docking in order to improve the enrichment of active compounds over decoys. Here we have investigated this in the context of the adenosine A2A receptor, where water molecules have previously been shown to be important for achieving high enrichment rates with docking, and where the positions of some binding site waters are known from a high-resolution crystal structure. The effect of these waters (both their presence and orientations) on VS enrichment was assessed using a carefully curated set of 299 high affinity A2A antagonists and 17,337 decoys. We show that including certain crystal waters greatly improves VS enrichment and that optimization of water hydrogen positions is needed in order to achieve the best results. We also show that waters derived from a molecular dynamics simulation - without any knowledge of crystallographic waters - can improve enrichments to a similar degree as the crystallographic waters, which makes this strategy applicable to structures without experimental knowledge of water positions. Finally, we used decision trees to select an ensemble of structures with different water molecule positions and orientations that outperforms any single structure with water molecules. The approach presented here is validated against independent test sets of A2A receptor antagonists and decoys from the literature. In general, this water optimization strategy could be applied to any target with waters-mediated protein-ligand interactions.
Collapse
Affiliation(s)
- Eelke B Lenselink
- Division of Medicinal Chemistry, Leiden Academic Centre for Drug Research, Leiden University , Einsteinweg 55, 2333 CC Leiden, The Netherlands
| | | | | | | | | |
Collapse
|
24
|
Jacobson KA. Crystal structures of the A 2A adenosine receptor and their use in medicinal chemistry. In Silico Pharmacol 2013; 1:22. [PMID: 24660138 PMCID: PMC3956660 DOI: 10.1186/2193-9616-1-22] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
New insights into drug design are derived from the X-ray crystallographic structures of G protein-coupled receptors (GPCRs), and the adenosine receptors (ARs) are at the forefront of this effort. The 3D knowledge of receptor binding and activation promises to enable drug discovery for GPCRs in general, and specifically for the ARs. The predictability of modeling based on the X-ray structures of the A2AAR has been well demonstrated in the identification, design and modification of both known and novel AR agonists and antagonists. It is expected that structure-based design of drugs acting through ARs will provide new avenues to clinically useful agents.
Collapse
|
25
|
Andrews SP, Brown GA, Christopher JA. Structure-Based and Fragment-Based GPCR Drug Discovery. ChemMedChem 2013; 9:256-75. [DOI: 10.1002/cmdc.201300382] [Citation(s) in RCA: 56] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2013] [Revised: 11/15/2013] [Indexed: 01/05/2023]
|
26
|
Furtmann N, Bajorath J. Evaluation of molecular model-based discovery of ecto-5′-nucleotidase inhibitors on the basis of X-ray structures. Bioorg Med Chem 2013; 21:6616-22. [DOI: 10.1016/j.bmc.2013.08.021] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2013] [Revised: 06/26/2013] [Accepted: 08/08/2013] [Indexed: 11/25/2022]
|
27
|
Drabczyńska A, Karcz T, Szymańska E, Köse M, Müller CE, Paskaleva M, Karolak-Wojciechowska J, Handzlik J, Yuzlenko O, Kieć-Kononowicz K. Synthesis, biological activity and molecular modelling studies of tricyclic alkylimidazo-, pyrimido- and diazepinopurinediones. Purinergic Signal 2013; 9:395-414. [PMID: 23543220 PMCID: PMC3757144 DOI: 10.1007/s11302-013-9358-3] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2012] [Accepted: 02/28/2013] [Indexed: 01/30/2023] Open
Abstract
Syntheses and biological activities of imidazo-, pyrimido- and diazepino[2,1-f]purinediones containing N-alkyl substituents (with straight, branched or unsaturated chains) are described. Tricyclic derivatives were synthesized by the cyclization of 8-bromo-substituted 7-(2-bromoethyl)-, 7-(3-chloropropyl)- or 7-(4-bromobutyl)-theophylline with primary amines under various conditions. Compound 22 with an ethenyl substituent was synthesized by dehydrohalogenation of 9-(2-bromoethyl)-1,3-dimethyltetrahydropyrimido[2,1-f]purinedione. The obtained derivatives (5-35) were initially evaluated for their affinity at rat A1 and A2A adenosine receptors (AR), showing moderate affinity for both adenosine receptor subtypes. The best ligands were diazepinopurinedione 28 (K i = 0.28 μM) with fivefold A2A selectivity and the non-selective A1/A2A AR ligand pyrimidopurinedione 35 (K i A1 = 0.28 μM and K i A2A = 0.30 μM). The compounds were also evaluated for their affinity at human A1, A2A, A2B and A3 ARs. All of the obtained compounds were docked to the A2A AR X-ray structure in complex with the xanthine-based, potent adenosine receptor antagonist-XAC. The likely interactions of imidazo-, pyrimido- and diazepino[2,1-f]purinediones with the residues forming the A2A binding pocket were discussed. Furthermore, the new compounds were tested in vivo as anticonvulsants in maximal electroshock, subcutaneous pentylenetetrazole (ScMet) and TOX tests in mice (i.p.). Pyrimidopurinediones showed anticonvulsant activity mainly in the ScMet test. The best derivative was compound 11, showing 100 % protection at a dose of 100 mg/kg without symptoms of neurotoxicity. Compounds 6, 7, 8 and 14 with short substituents showed neurotoxicity and caused death. In rat tests (p.o.), 9 was characterized by a high protection index (>13.3). AR affinity did not apparently correlate with the antiepileptic potency of the compounds.
Collapse
Affiliation(s)
- Anna Drabczyńska
- Department of Technology and Biotechnology of Drugs, Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland
| | - Tadeusz Karcz
- Department of Technology and Biotechnology of Drugs, Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland
| | - Ewa Szymańska
- Department of Technology and Biotechnology of Drugs, Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland
| | - Meryem Köse
- PharmaCenter Bonn, Pharmaceutical Institute, University of Bonn, An der Immenburg 4, 53121 Bonn, Germany
| | - Christa E. Müller
- PharmaCenter Bonn, Pharmaceutical Institute, University of Bonn, An der Immenburg 4, 53121 Bonn, Germany
| | - Minka Paskaleva
- PharmaCenter Bonn, Pharmaceutical Institute, University of Bonn, An der Immenburg 4, 53121 Bonn, Germany
| | | | - Jadwiga Handzlik
- Department of Technology and Biotechnology of Drugs, Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland
| | - Olga Yuzlenko
- Department of Technology and Biotechnology of Drugs, Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland
- Present Address: Department of Chemistry, City College of New York, City University of New York, 160 Convent Avenue, New York, NY 10031 USA
| | - Katarzyna Kieć-Kononowicz
- Department of Technology and Biotechnology of Drugs, Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9, 30-688 Kraków, Poland
| |
Collapse
|
28
|
Piccoli S, Suku E, Garonzi M, Giorgetti A. Genome-wide Membrane Protein Structure Prediction. Curr Genomics 2013; 14:324-9. [PMID: 24403851 PMCID: PMC3763683 DOI: 10.2174/13892029113149990009] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2013] [Revised: 07/19/2013] [Accepted: 07/22/2013] [Indexed: 01/25/2023] Open
Abstract
Transmembrane proteins allow cells to extensively communicate with the external world in a very accurate and specific way. They form principal nodes in several signaling pathways and attract large interest in therapeutic intervention, as the majority pharmaceutical compounds target membrane proteins. Thus, according to the current genome annotation methods, a detailed structural/functional characterization at the protein level of each of the elements codified in the genome is also required. The extreme difficulty in obtaining high-resolution three-dimensional structures, calls for computational approaches. Here we review to which extent the efforts made in the last few years, combining the structural characterization of membrane proteins with protein bioinformatics techniques, could help describing membrane proteins at a genome-wide scale. In particular we analyze the use of comparative modeling techniques as a way of overcoming the lack of high-resolution three-dimensional structures in the human membrane proteome.
Collapse
Affiliation(s)
- Stefano Piccoli
- Applied Bioinformatics Group, Dept. of Biotechnology, University of Verona, strada Le grazie 15, 37134, Verona,
Italy
| | - Eda Suku
- Applied Bioinformatics Group, Dept. of Biotechnology, University of Verona, strada Le grazie 15, 37134, Verona,
Italy
| | - Marianna Garonzi
- Applied Bioinformatics Group, Dept. of Biotechnology, University of Verona, strada Le grazie 15, 37134, Verona,
Italy
| | - Alejandro Giorgetti
- Applied Bioinformatics Group, Dept. of Biotechnology, University of Verona, strada Le grazie 15, 37134, Verona,
Italy
- German Research School for Simulation Sciences, Juelich, Germany
- Center for Biomedical Computing (CBMC), University of Verona, strada Le grazie 8, 37134, Verona, Italy
| |
Collapse
|
29
|
Jacobson KA. Structure-based approaches to ligands for G-protein-coupled adenosine and P2Y receptors, from small molecules to nanoconjugates. J Med Chem 2013; 56:3749-67. [PMID: 23597047 PMCID: PMC3701956 DOI: 10.1021/jm400422s] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Adenosine receptor (ARs) and P2Y receptors (P2YRs) that respond to extracellular nucleosides/nucleotides are associated with new directions for therapeutics. The X-ray structures of the A2AAR complexes with agonists and antagonists are examined in relationship to the G-protein-coupled receptor (GPCR) superfamily and applied to drug discovery. Much of the data on AR ligand structure from early SAR studies now are explainable from the A2AAR X-ray crystallography. The ligand-receptor interactions in related GPCR complexes can be identified by means of modeling approaches, e.g., molecular docking. Thus, molecular recognition in binding and activation processes has been studied effectively using homology modeling and applied to ligand design. Virtual screening has yielded new nonnucleoside AR antagonists, and existing ligands have been improved with knowledge of the receptor interactions. New agonists are being explored for central nervous system and peripheral therapeutics based on in vivo activity, such as chronic neuropathic pain. Ligands for receptors more distantly related to the X-ray template, i.e., P2YRs, have been introduced and are mainly used as pharmacological tools for elucidating the physiological role of extracellular nucleotides. Other ligand tools for drug discovery include fluorescent probes, radioactive probes, multivalent probes, and functionalized nanoparticles.
Collapse
Affiliation(s)
- Kenneth A Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health , Bethesda, Maryland 20892, USA.
| |
Collapse
|
30
|
Hagiwara Y, Ohno K, Kamohara M, Takasaki J, Watanabe T, Fukunishi Y, Nakamura H, Orita M. Molecular modeling of vasopressin receptor and in silico screening of V1b receptor antagonists. Expert Opin Drug Discov 2013; 8:951-64. [PMID: 23682717 DOI: 10.1517/17460441.2013.799134] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
INTRODUCTION G protein-coupled receptors (GPCRs) are integral membrane proteins which contain seven-transmembrane-spanning alpha-helices. GPCR-mediated signaling has been associated with various human diseases, positioning GPCRs as attractive targets in the drug discovery field. Recently, through advances in protein engineering and crystallography, the number of resolved GPCR structures has increased dramatically. This growing availability of GPCR structures has greatly accelerated structure-based drug design (SBDD) and in silico screening for GPCR-targeted drug discovery. AREAS COVERED The authors introduce the current status of X-ray crystallography of GPCRs and what has been revealed from the resolved crystal structures. They also review the recent advances in SBDD and in silico screening for GPCR-targeted drug discovery and discuss a docking study, using homology modeling, with the discovery of potent antagonists of the vasopressin 1b receptor. EXPERT OPINION Several innovative protein engineering techniques and crystallographic methods have greatly accelerated SBDD, not only for already-resolved GPCRs but also for those structures which remain unclear. These technological advances are expected to enable the determination of GPCR-fragment complexes, making it practical to perform fragment-based drug discovery. This paves the way for a new era of GPCR-targeted drug discovery.
Collapse
Affiliation(s)
- Yohsuke Hagiwara
- Drug Discovery Research, Astellas Pharma, Inc., 21, Miyukigaoka, Ibaraki 305-8585, Japan.
| | | | | | | | | | | | | | | |
Collapse
|
31
|
Insight into the binding mode and the structural features of the pyrimidine derivatives as human A2A adenosine receptor antagonists. Biosystems 2013; 115:13-22. [PMID: 23665268 DOI: 10.1016/j.biosystems.2013.04.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2013] [Revised: 04/07/2013] [Accepted: 04/24/2013] [Indexed: 12/12/2022]
Abstract
The interaction of 278 monocyclic and bicyclic pyrimidine derivatives with human A2A adenosine receptor (AR) was investigated by employing molecular dynamics, thermodynamic analysis and three-dimensional quantitative structure-activity relationship (3D-QSAR) approaches. The binding analysis reveals that the pyrimidine derivatives are anchored in TM2, 3, 5, 6 and 7 of A2A AR by the aromatic stacking and hydrogen bonding interactions. The key residues involving Phe168, Glu169, and Asn253 stabilize the monocyclic and bicyclic cores of inhibitors. The thermodynamic analysis by molecular mechanics/Poisson Boltzmann surface area (MM-PBSA) approach also confirms the reasonableness of the binding modes. In addition, the ligand-/receptor-based comparative molecular similarity indices analysis (CoMSIA) models of high statistical significance were generated and the resulting contour maps correlate well with the structural features of the antagonists essential for high A2A AR affinity. A minor/bulky group with negative charge at C2/C6 of pyrimidine ring respectively enhances the activity for all these pyrimidine derivatives. Particularly, the higher electron density of the ring in the bicyclic derivatives, the more potent the antagonists. The obatined results might be helpful in rational design of novel candidate of A2A adenosine receptor antagonist for treatment of Parkinson's disease.
Collapse
|
32
|
Homology models of melatonin receptors: challenges and recent advances. Int J Mol Sci 2013; 14:8093-121. [PMID: 23584026 PMCID: PMC3645733 DOI: 10.3390/ijms14048093] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2013] [Revised: 03/28/2013] [Accepted: 03/28/2013] [Indexed: 12/15/2022] Open
Abstract
Melatonin exerts many of its actions through the activation of two G protein-coupled receptors (GPCRs), named MT1 and MT2. So far, a number of different MT1 and MT2 receptor homology models, built either from the prototypic structure of rhodopsin or from recently solved X-ray structures of druggable GPCRs, have been proposed. These receptor models differ in the binding modes hypothesized for melatonin and melatonergic ligands, with distinct patterns of ligand-receptor interactions and putative bioactive conformations of ligands. The receptor models will be described, and they will be discussed in light of the available information from mutagenesis experiments and ligand-based pharmacophore models. The ability of these ligand-receptor complexes to rationalize structure-activity relationships of known series of melatonergic compounds will be commented upon.
Collapse
|
33
|
Pala D, Beuming T, Sherman W, Lodola A, Rivara S, Mor M. Structure-based virtual screening of MT2 melatonin receptor: influence of template choice and structural refinement. J Chem Inf Model 2013; 53:821-35. [PMID: 23541165 DOI: 10.1021/ci4000147] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Developing GPCR homology models for structure-based virtual screening requires the choice of a suitable template and refinement of binding site residues. We explored this systematically for the MT2 melatonin receptor, with the aim to build a receptor homology model that is optimized for the enrichment of active melatoninergic ligands. A set of 12 MT2 melatonin receptor models was built using different GPCR X-ray structural templates and submitted to a virtual screening campaign on a set of compounds composed of 29 known melatonin receptor ligands and 2560 drug-like decoys. To evaluate the effect of including a priori information in receptor models, 12 representative melatonin receptor ligands were placed into the MT2 receptor models in poses consistent with known mutagenesis data and with assessed pharmacophore models. The receptor structures were then adapted to the ligands by induced-fit docking. Most of the 144 ligand-adapted MT2 receptor models showed significant improvements in screening enrichments compared to the unrefined homology models, with some template/refinement combinations giving excellent enrichment factors. The discriminating ability of the models was further tested on the 29 active ligands plus a set of 21 inactive or low-affinity compounds from the same chemical classes. Rotameric states of side chains for some residues, presumed to be involved in the binding process, were correlated with screening effectiveness, suggesting the existence of specific receptor conformations able to recognize active compounds. The top MT2 receptor model was able to identify 24 of 29 active ligands among the first 2% of the screened database. This work provides insights into the use of refined GPCR homology models for virtual screening.
Collapse
Affiliation(s)
- Daniele Pala
- Dipartimento di Farmacia, Università degli Studi di Parma, Parco Area delle Scienze 27/A, I-43124 Parma, Italy
| | | | | | | | | | | |
Collapse
|
34
|
Topiol S. X-ray structural information of GPCRs in drug design: what are the limitations and where do we go? Expert Opin Drug Discov 2013; 8:607-20. [PMID: 23537065 DOI: 10.1517/17460441.2013.783815] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
INTRODUCTION In 2007, the X-ray structural determination of non-rhodopsin G-Protein coupled receptors (GPCRs), considered the most extensively targeted protein class for marketed drugs, commenced. With the relatively rapid availability of additional structures, an assessment of the progression made is needed in addition to the assessment of the understandings gleaned, deployment successes and forthcoming prospects. AREAS COVERED The author reviews the approaches and tools that have made it possible to determine the three dimensional structures of GPCRs using X-ray crystallography. Furthermore, the author describes the methods suited for crystallization of membrane bound GPCR proteins including the lipidic cubic phase and various protein modification approaches. The author also provides highlights, from the literature, of the structures determined to date including targets solved, the nature of the content provided (such as selectivity, activating vs. inactivating determinants) and how these structural features relate to drug design strategies. EXPERT OPINION The GPCR X-ray structures that have been so far determined have yielded significant information. This has presented dramatic evidence concerning their ability to impact the discovery of compounds through their action as traditional, orthosteric modulators. It is, however, noted that more challenging design strategies, such as identifying biased agonists and the use of sites remote from the orthosteric site for allosteric modulation, are still in their infancy.
Collapse
Affiliation(s)
- Sid Topiol
- 3D-2Drug LLC, PO Box 184, Fair Lawn, NJ 07410, USA.
| |
Collapse
|
35
|
Roy A, Taddese B, Vohra S, Thimmaraju PK, Illingworth CJR, Simpson LM, Mukherjee K, Reynolds CA, Chintapalli SV. Identifying subset errors in multiple sequence alignments. J Biomol Struct Dyn 2013; 32:364-71. [PMID: 23527867 DOI: 10.1080/07391102.2013.770371] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
Multiple sequence alignment (MSA) accuracy is important, but there is no widely accepted method of judging the accuracy that different alignment algorithms give. We present a simple approach to detecting two types of error, namely block shifts and the misplacement of residues within a gap. Given a MSA, subsets of very similar sequences are generated through the use of a redundancy filter, typically using a 70-90% sequence identity cut-off. Subsets thus produced are typically small and degenerate, and errors can be easily detected even by manual examination. The errors, albeit minor, are inevitably associated with gaps in the alignment, and so the procedure is particularly relevant to homology modelling of protein loop regions. The usefulness of the approach is illustrated in the context of the universal but little known [K/R]KLH motif that occurs in intracellular loop 1 of G protein coupled receptors (GPCR); other issues relevant to GPCR modelling are also discussed.
Collapse
Affiliation(s)
- Aparna Roy
- a School of Biological Sciences, University of Essex , Wivenhoe Park, Colchester , CO4 3SQ , UK
| | | | | | | | | | | | | | | | | |
Collapse
|
36
|
Levit A, Barak D, Behrens M, Meyerhof W, Niv MY. Homology model-assisted elucidation of binding sites in GPCRs. Methods Mol Biol 2013; 914:179-205. [PMID: 22976029 DOI: 10.1007/978-1-62703-023-6_11] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
G protein-coupled receptors (GPCRs) are important mediators of cell signaling and a major family of drug targets. Despite recent breakthroughs, experimental elucidation of GPCR structures remains a formidable challenge. Homology modeling of 3D structures of GPCRs provides a practical tool for elucidating the structural determinants governing the interactions of these important receptors with their ligands. The working model of the binding site can then be used for virtual screening of additional ligands that may fit this site, for determining and comparing specificity profiles of related receptors, and for structure-based design of agonists and antagonists. The current review presents the protocol and enumerates the steps for modeling and validating the residues involved in ligand binding. The main stages include (a) modeling the receptor structure using an automated fragment-based approach, (b) predicting potential binding pockets, (c) docking known binders, (d) analyzing predicted interactions and comparing with positions that have been shown to bind ligands in other receptors, (e) validating the structural model by mutagenesis.
Collapse
Affiliation(s)
- Anat Levit
- Institute of Biochemistry, Food Science, and Nutrition, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University, Rehovot, Israel
| | | | | | | | | |
Collapse
|
37
|
Kooistra AJ, Roumen L, Leurs R, de Esch IJ, de Graaf C. From Heptahelical Bundle to Hits from the Haystack. Methods Enzymol 2013; 522:279-336. [DOI: 10.1016/b978-0-12-407865-9.00015-7] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
38
|
Taddese B, Simpson LM, Wall ID, Blaney FE, Reynolds CA. Modeling Active GPCR Conformations. Methods Enzymol 2013; 522:21-35. [DOI: 10.1016/b978-0-12-407865-9.00002-9] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
39
|
Andrews SP, Tehan B. Stabilised G protein-coupled receptors in structure-based drug design: a case study with adenosine A2A receptor. MEDCHEMCOMM 2013. [DOI: 10.1039/c2md20164j] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The first example of structure-based drug design with stabilised GPCRs has enabled the identification of a preclinical candidate for the treatment of Parkinson's disease.
Collapse
Affiliation(s)
| | - Benjamin Tehan
- Heptares Therapeutics Limited
- BioPark
- Welwyn Garden City
- UK
| |
Collapse
|
40
|
Beuming T, Sherman W. Current Assessment of Docking into GPCR Crystal Structures and Homology Models: Successes, Challenges, and Guidelines. J Chem Inf Model 2012; 52:3263-77. [PMID: 23121495 DOI: 10.1021/ci300411b] [Citation(s) in RCA: 69] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Thijs Beuming
- Schrödinger, Inc., 120
West 45th Street, New York, New York, United States
| | - Woody Sherman
- Schrödinger, Inc., 120
West 45th Street, New York, New York, United States
| |
Collapse
|
41
|
Sirci F, Goracci L, Rodríguez D, van Muijlwijk-Koezen J, Gutiérrez-de-Terán H, Mannhold R. Ligand-, structure- and pharmacophore-based molecular fingerprints: a case study on adenosine A1, A2A, A2B, and A3 receptor antagonists. J Comput Aided Mol Des 2012; 26:1247-66. [DOI: 10.1007/s10822-012-9612-8] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2012] [Accepted: 09/20/2012] [Indexed: 10/27/2022]
|
42
|
Jacobson KA, Balasubramanian R, Deflorian F, Gao ZG. G protein-coupled adenosine (P1) and P2Y receptors: ligand design and receptor interactions. Purinergic Signal 2012; 8:419-36. [PMID: 22371149 PMCID: PMC3360101 DOI: 10.1007/s11302-012-9294-7] [Citation(s) in RCA: 67] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2011] [Accepted: 01/30/2012] [Indexed: 12/17/2022] Open
Abstract
The medicinal chemistry and pharmacology of the four subtypes of adenosine receptors (ARs) and the eight subtypes of P2Y receptors (P2YRs, activated by a range of purine and pyrimidine mono- and dinucleotides) has recently advanced significantly leading to selective ligands. X-ray crystallographic structures of both agonist- and antagonist-bound forms of the A(2A)AR have provided unprecedented three-dimensional detail concerning molecular recognition in the binding site and the conformational changes in receptor activation. It is apparent that this ubiquitous cell signaling system has implications for understanding and treating many diseases. ATP and other nucleotides are readily released from intracellular sources under conditions of injury and organ stress, such as hypoxia, ischemia, or mechanical stress, and through channels and vesicular release. Adenosine may be generated extracellularly or by cellular release. Therefore, depending on pathophysiological factors, in a given tissue, there is often a tonic activation of one or more of the ARs or P2YRs that can be modulated by exogenous agents for a beneficial effect. Thus, this field has provided fertile ground for pharmaceutical development, leading to clinical trials of selective receptor ligands as imaging agents or for conditions including cardiac arrhythmias, ischemia/reperfusion injury, diabetes, pain, thrombosis, Parkinson's disease, rheumatoid arthritis, psoriasis, dry eye disease, pulmonary diseases such as cystic fibrosis, glaucoma, cancer, chronic hepatitis C, and other diseases.
Collapse
Affiliation(s)
- Kenneth A Jacobson
- Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA.
| | | | | | | |
Collapse
|
43
|
Jacobson KA, Costanzi S. New insights for drug design from the X-ray crystallographic structures of G-protein-coupled receptors. Mol Pharmacol 2012; 82:361-71. [PMID: 22695719 DOI: 10.1124/mol.112.079335] [Citation(s) in RCA: 72] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
Methodological advances in X-ray crystallography have made possible the recent solution of X-ray structures of pharmaceutically important G protein-coupled receptors (GPCRs), including receptors for biogenic amines, peptides, a nucleoside, and a sphingolipid. These high-resolution structures have greatly increased our understanding of ligand recognition and receptor activation. Conformational changes associated with activation common to several receptors entail outward movements of the intracellular side of transmembrane helix 6 (TM6) and movements of TM5 toward TM6. Movements associated with specific agonists or receptors have also been described [e.g., extracellular loop (EL) 3 in the A(2A) adenosine receptor]. The binding sites of different receptors partly overlap but differ significantly in ligand orientation, depth, and breadth of contact areas in TM regions and the involvement of the ELs. A current challenge is how to use this structural information for the rational design of novel potent and selective ligands. For example, new chemotypes were discovered as antagonists of various GPCRs by subjecting chemical libraries to in silico docking in the X-ray structures. The vast majority of GPCR structures and their ligand complexes are still unsolved, and no structures are known outside of family A GPCRs. Molecular modeling, informed by supporting information from site-directed mutagenesis and structure-activity relationships, has been validated as a useful tool to extend structural insights to related GPCRs and to analyze docking of other ligands in already crystallized GPCRs.
Collapse
Affiliation(s)
- Kenneth A Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institutes of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892-0810, USA.
| | | |
Collapse
|
44
|
Tosh DK, Deflorian F, Phan K, Gao ZG, Wan TC, Gizewski E, Auchampach JA, Jacobson KA. Structure-guided design of A(3) adenosine receptor-selective nucleosides: combination of 2-arylethynyl and bicyclo[3.1.0]hexane substitutions. J Med Chem 2012; 55:4847-60. [PMID: 22559880 PMCID: PMC3371665 DOI: 10.1021/jm300396n] [Citation(s) in RCA: 69] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
(N)-Methanocarba adenosine 5'-methyluronamides containing known A(3) AR (adenosine receptor)-enhancing modifications, i.e., 2-(arylethynyl)adenine and N(6)-methyl or N(6)-(3-substituted-benzyl), were nanomolar full agonists of human (h) A(3)AR and highly selective (K(i) ∼0.6 nM, N(6)-methyl 2-(halophenylethynyl) analogues 13 and 14). Combined 2-arylethynyl-N(6)-3-chlorobenzyl substitutions preserved A(3)AR affinity/selectivity in the (N)-methanocarba series (e.g., 3,4-difluoro full agonist MRS5698 31, K(i) 3 nM, human and mouse A(3)) better than that for ribosides. Polyaromatic 2-ethynyl N(6)-3-chlorobenzyl analogues, such as potent linearly extended 2-p-biphenylethynyl MRS5679 34 (K(i) hA(3) 3.1 nM; A(1), A(2A), inactive) and fluorescent 1-pyrene adduct MRS5704 35 (K(i) hA(3) 68.3 nM), were conformationally rigid; receptor docking identified a large, mainly hydrophobic binding region. The vicinity of receptor-bound C2 groups was probed by homology modeling based on recent X-ray structure of an agonist-bound A(2A)AR, with a predicted helical rearrangement requiring an agonist-specific outward displacement of TM2 resembling opsin. Thus, the X-ray structure of related A(2A)AR is useful in guiding the design of new A(3)AR agonists.
Collapse
Affiliation(s)
- Dilip K. Tosh
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892
| | - Francesca Deflorian
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892
| | - Khai Phan
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892
| | - Zhan-Guo Gao
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892
| | - Tina C. Wan
- Department of Pharmacology, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226
| | - Elizabeth Gizewski
- Department of Pharmacology, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226
| | - John A. Auchampach
- Department of Pharmacology, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226
| | - Kenneth A. Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892
| |
Collapse
|
45
|
Zhu M, Li M. Revisiting the homology modeling of G-protein coupled receptors: β1-adrenoceptor as an example. MOLECULAR BIOSYSTEMS 2012; 8:1686-93. [PMID: 22454032 DOI: 10.1039/c2mb05491d] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
G-protein coupled receptors (GPCRs) are recognized to constitute the largest family of membrane proteins. Due to the disproportion in the quantity of crystal structures and their amino acid sequences, homology modeling contributes a reasonable and feasible approach to GPCR theoretical coordinates. With the brand new crystal structures resolved recently, herein we deliberated how to designate them as templates to carry out homology modeling in four aspects: (1) various sequence alignment methods; (2) protein weight matrix; (3) different sets of multiple templates; (4) active and inactive state of templates. The accuracy of models was evaluated by comparing the similarity of stereo conformation and molecular docking results between models and the experimental structure of Meleagris gallopavo β(1)-adrenergic receptor (Mg_Adrb1) that we desired to develop as an example. Our results proposed that: (1) Cobalt and MAFFT, two algorithms of sequence alignment, were suitable for single- and multiple-template modeling, respectively; (2) Blosum30 is applicable to align sequences in the case of low sequence identity; (3) multiple-template modeling is not always better than single-template one; (4) the state of template is an influential factor in simulating the GPCR structures as well.
Collapse
Affiliation(s)
- Mengyuan Zhu
- Department of Medicinal Chemistry, School of Pharmacy, Shandong University, Jinan, Shandong, China
| | | |
Collapse
|
46
|
Tang H, Wang XS, Hsieh JH, Tropsha A. Do crystal structures obviate the need for theoretical models of GPCRs for structure-based virtual screening? Proteins 2012; 80:1503-21. [PMID: 22275072 DOI: 10.1002/prot.24035] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2011] [Revised: 12/06/2011] [Accepted: 12/29/2011] [Indexed: 12/28/2022]
Abstract
Recent highly expected structural characterizations of agonist-bound and antagonist-bound beta-2 adrenoreceptor (β2AR) by X-ray crystallography have been widely regarded as critical advances to enable more effective structure-based discovery of GPCRs ligands. It appears that this very important development may have undermined many previous efforts to develop 3D theoretical models of GPCRs. To address this question directly, we have compared several historical β2AR models versus the inactive state and nanobody-stabilized active state of β2AR crystal structures in terms of their structural similarity and effectiveness of use in virtual screening for β2AR specific agonists and antagonists. Theoretical models, incluing both homology and de novo types, were collected from five different groups who have published extensively in the field of GPCRs modeling. All models were built before X-ray structures became available. In general, β2AR theoretical models differ significantly from the crystal structure in terms of TMH definition and the global packing. Nevertheless, surprisingly, several models afforded hit rates resulting from virtual screening of large chemical library enriched by known β2AR ligands that exceeded those using X-ray structures. The hit rates were particularly higher for agonists. Furthemore, the screening performance of models is associated with local structural quality, such as the RMSDs for binding pocket residues and the ability to capture accurately, most if not all critical protein/ligand interactions. These results suggest that carefully built models of GPCRs could capture critical chemical and structural features of the binding pocket, and thus may be even more useful for practical structure-based drug discovery than X-ray structures.
Collapse
Affiliation(s)
- Hao Tang
- Laboratory for Molecular Modeling, Division of Medicinal Chemistry and Natural Products and Carolina Exploratory Center for Cheminformatics Research, School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599-7360, USA
| | | | | | | |
Collapse
|
47
|
Prakash A, Luthra PM. Insilico study of the A2AR–D2R kinetics and interfacial contact surface for heteromerization. Amino Acids 2012; 43:1451-64. [DOI: 10.1007/s00726-012-1218-x] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2011] [Accepted: 01/04/2012] [Indexed: 12/28/2022]
|
48
|
Wooding S. Signatures of natural selection in a primate bitter taste receptor. J Mol Evol 2012; 73:257-65. [PMID: 22218679 DOI: 10.1007/s00239-011-9481-0] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2011] [Accepted: 12/14/2011] [Indexed: 01/15/2023]
Abstract
Bitter taste receptors (TAS2Rs) enable animals to detect and avoid toxins in the environment, including noxious defense compounds produced by plants. This suggests that TAS2Rs are under complex pressures from natural selection. To investigate these pressures, we examined signatures of selection in the primate TAS2R38 gene. Whole-gene (1,002 bp) sequences from 40 species representing all major primate taxa uncovered extensive variation. Nucleotide substitutions occurred at 448 positions, resulting in 201 amino acid changes. Two single-nucleotide deletions, one three-nucleotide in-frame deletion, and one premature stop codon were also observed. The rate of non-synonymous substitution (ω = dN/dS), was high in TAS2R38 (ω = 0.60) compared to other genes, but significantly lower than expected under neutrality (P = 4.0 × 10(-9)), indicating that purifying selection has maintained the basic structure of the receptor. However, differences were present among receptor subregions. Non-synonymous rates were significantly lower than expected in transmembrane domains (ω = 0.55, P = 1.18 × 10(-12)) and internal loops (ω = 0.51, P = 7.04 × 10(-5)), but not external loops (ω = 1.16, P = 0.53), and evidence of positive selection was found in external loop 2, which exhibited a high rate (ω = 2.53) consistent with rapid shifts in ligand targeting. These patterns point to a history of rapid yet constrained change in bitter taste responses in the course of primate evolution.
Collapse
Affiliation(s)
- Stephen Wooding
- McDermott Center for Human Growth and Development, University of Texas Southwestern Medical Center, Dallas, TX 75390-8591, USA.
| |
Collapse
|
49
|
Doré AS, Robertson N, Errey JC, Ng I, Hollenstein K, Tehan B, Hurrell E, Bennett K, Congreve M, Magnani F, Tate CG, Weir M, Marshall FH. Structure of the adenosine A(2A) receptor in complex with ZM241385 and the xanthines XAC and caffeine. Structure 2011; 19:1283-93. [PMID: 21885291 DOI: 10.1016/j.str.2011.06.014] [Citation(s) in RCA: 443] [Impact Index Per Article: 34.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2011] [Revised: 06/08/2011] [Accepted: 06/13/2011] [Indexed: 02/01/2023]
Abstract
Methylxanthines, including caffeine and theophylline, are among the most widely consumed stimulant drugs in the world. These effects are mediated primarily via blockade of adenosine receptors. Xanthine analogs with improved properties have been developed as potential treatments for diseases such as Parkinson's disease. Here we report the structures of a thermostabilized adenosine A(2A) receptor in complex with the xanthines xanthine amine congener and caffeine, as well as the A(2A) selective inverse agonist ZM241385. The receptor is crystallized in the inactive state conformation as defined by the presence of a salt bridge known as the ionic lock. The complete third intracellular loop, responsible for G protein coupling, is visible consisting of extended helices 5 and 6. The structures provide new insight into the features that define the ligand binding pocket of the adenosine receptor for ligands of diverse chemotypes as well as the cytoplasmic regions that interact with signal transduction proteins.
Collapse
Affiliation(s)
- Andrew S Doré
- Heptares Therapeutics Ltd, BioPark, Welwyn Garden City, Herts, AL7 3AX, UK
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
50
|
Lane JR, Klein Herenbrink C, van Westen GJP, Spoorendonk JA, Hoffmann C, IJzerman AP. A novel nonribose agonist, LUF5834, engages residues that are distinct from those of adenosine-like ligands to activate the adenosine A(2a) receptor. Mol Pharmacol 2011; 81:475-87. [PMID: 22188926 DOI: 10.1124/mol.111.075937] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
The recent publication of both the antagonist- and agonist-bound structures of the adenosine A(2A) receptor have revealed much about how a ligand may bind to a receptor and cause the conformational changes associated with agonist-mediated activation. In particular, the agonist-bound structure revealed key interactions between the ribose group of adenosine-derived agonists and amino acids in the receptor binding pocket that lead to receptor activation. However, agonists without a ribose group also exist, and we wondered whether such compounds occupy the same agonist binding site. Therefore we used a mutagenesis approach in this study to investigate the mode of binding of 2-amino-4-(4-hydroxyphenyl)- 6-(1H-imidazol-2-ylmethylsulfanyl)pyridine-3,5-dicarbonitrile (LUF5834), a potent partial agonist without a ribose moiety, compared with the adenosine-derived reference agonist 2-[p-(2-carboxyethyl)phenyl-ethylamino]-5'-N-ethylcarboxamidoadenosine (CGS21680). Mutation of the orthosteric residue Phe168 to alanine abrogated the function of both agonists. However, mutation to alanine of residues Thr88 and Ser277 shown by the crystal structures to interact with the ribose group of adenosine-like ligands had no effect on the potency of LUF5834. Furthermore, alanine mutation of Asn253, which makes a hydrogen-bonding interaction with the exocyclic nitrogen of the adenine ring, had minimal effect on LUF5834 affinity but removed agonist activity of this ligand. Mutation of other residues, such as the highly conserved Trp246 or Glu13, had significant deleterious effects on the function of CGS21680 but little effect on LUF5834. In summary, our findings suggest that this class of agonist interacts with distinct residues to activate the receptor compared with classic adenosine derived agonists.
Collapse
Affiliation(s)
- J Robert Lane
- Division of Medicinal Chemistry, Leiden/Amsterdam Center for Drug Research, Leiden University, Leiden, the Netherlands.
| | | | | | | | | | | |
Collapse
|