1
|
Scardino V, Di Filippo JI, Cavasotto CN. How good are AlphaFold models for docking-based virtual screening? iScience 2023; 26:105920. [PMID: 36686396 PMCID: PMC9852548 DOI: 10.1016/j.isci.2022.105920] [Citation(s) in RCA: 43] [Impact Index Per Article: 43.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Revised: 11/12/2022] [Accepted: 12/28/2022] [Indexed: 12/31/2022] Open
Abstract
A crucial component in structure-based drug discovery is the availability of high-quality three-dimensional structures of the protein target. Whenever experimental structures were not available, homology modeling has been, so far, the method of choice. Recently, AlphaFold (AF), an artificial-intelligence-based protein structure prediction method, has shown impressive results in terms of model accuracy. This outstanding success prompted us to evaluate how accurate AF models are from the perspective of docking-based drug discovery. We compared the high-throughput docking (HTD) performance of AF models with their corresponding experimental PDB structures using a benchmark set of 22 targets. The AF models showed consistently worse performance using four docking programs and two consensus techniques. Although AlphaFold shows a remarkable ability to predict protein architecture, this might not be enough to guarantee that AF models can be reliably used for HTD, and post-modeling refinement strategies might be key to increase the chances of success.
Collapse
Affiliation(s)
- Valeria Scardino
- Meton AI, Inc, Wilmington, DE 19801, USA
- Austral Institute for Applied Artificial Intelligence, Universidad Austral, Pilar, Buenos Aires, Argentina
| | - Juan I. Di Filippo
- Austral Institute for Applied Artificial Intelligence, Universidad Austral, Pilar, Buenos Aires, Argentina
- Computational Drug Design and Biomedical Informatics Laboratory, Instituto de Investigaciones en Medicina Traslacional (IIMT), Universidad Austral-CONICET, Pilar, Buenos Aires, Argentina
| | - Claudio N. Cavasotto
- Austral Institute for Applied Artificial Intelligence, Universidad Austral, Pilar, Buenos Aires, Argentina
- Computational Drug Design and Biomedical Informatics Laboratory, Instituto de Investigaciones en Medicina Traslacional (IIMT), Universidad Austral-CONICET, Pilar, Buenos Aires, Argentina
- Facultad de Ciencias Biomédicas, and Facultad de Ingeniería, Universidad Austral, Pilar, Buenos Aires, Argentina
| |
Collapse
|
2
|
Xu Q, Zhao Z, Liang P, Wang S, Li F, Jin S, Zhang J. Identification of novel nematode succinate dehydrogenase inhibitors: Virtual screening based on ligand-pocket interactions. Chem Biol Drug Des 2023; 101:9-23. [PMID: 34981652 DOI: 10.1111/cbdd.14019] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Revised: 12/10/2021] [Accepted: 12/29/2021] [Indexed: 12/15/2022]
Abstract
To discover new nematicidal succinate dehydrogenase (SDH) inhibitors with novel structures, we conducted a virtual screening of the ChemBridge library with 1.7 million compounds based on ligand-pocket interactions. The homology model of Caenorhabditis elegans SDH was established, along with a pharmacophore model based on ligand-pocket interactions. After the pharmacophore-based and docking-based screening, 19 compounds were selected for the subsequent enzymatic assays. The results showed that compound 1 (ID: 7607321) exhibited inhibitory activity against SDH with a determined IC50 value of 19.6 μM. Structural modifications and nematicidal activity studies were then carried out, which provided further evidence that compound 1 exhibited excellent nematicidal activity. Molecular dynamics simulations were then conducted to investigate the underlying molecular basis for the potency of these inhibitors against SDH. This work provides a reliable strategy and useful information for the future design of nematode SDH inhibitors.
Collapse
Affiliation(s)
- Qingbo Xu
- Innovation Center of Pesticide Research, Department of Applied Chemistry, College of Science, China Agricultural University, Beijing, China
| | - Zhixiang Zhao
- Innovation Center of Pesticide Research, Department of Applied Chemistry, College of Science, China Agricultural University, Beijing, China
| | - Peibo Liang
- Innovation Center of Pesticide Research, Department of Applied Chemistry, College of Science, China Agricultural University, Beijing, China
| | - Simin Wang
- Innovation Center of Pesticide Research, Department of Applied Chemistry, College of Science, China Agricultural University, Beijing, China
| | - Fang Li
- Innovation Center of Pesticide Research, Department of Applied Chemistry, College of Science, China Agricultural University, Beijing, China
| | - Shuhui Jin
- Innovation Center of Pesticide Research, Department of Applied Chemistry, College of Science, China Agricultural University, Beijing, China
| | - Jianjun Zhang
- Innovation Center of Pesticide Research, Department of Applied Chemistry, College of Science, China Agricultural University, Beijing, China
| |
Collapse
|
3
|
Chen L, Hu B, Wang H, Li W, Wang S, Luan J, Liu H, Wang J, Cheng M. Selectivity mechanism of muscarinic acetylcholine receptor antagonism through in silico investigation. Phys Chem Chem Phys 2022; 24:26269-26287. [PMID: 36281693 DOI: 10.1039/d2cp02972c] [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: 06/16/2023]
Abstract
Structures of muscarinic acetylcholine receptors illustrate the strikingly high degree of homology of the residues among isoforms, thus leading to difficulty in achieving subtype selectivity when targeting these receptors and causing undesired side effects when treating the corresponding diseases. Considering the urgent need for more selective and potency therapies, this study is aimed at revealing the selectivity mechanism against M4/5 via in silico strategies, revealing crucial molecular interactions such as hydrogen bonds and pi-cation interaction formed between the key residues TYR416, ASN417, and TRP435 of M4, respectively, hydrophobic pocket formed by the key residues, especially CYS484 of M5. Besides, the water around TYR416M4 and ASN459M5, which can be replaced by substituent groups which can form the hydrogen bond interaction network by simulated bridging water and the water around ASP112M4, whose replacement maybe not contribute to the increase in binding affinity of the compound, may affect the inhibitory selectivity among M4/5 in the aspect of the solvent. Moreover, from the point of inhibitors, compounds with a positively ionizable group could selectively bind to M4 receptors, while hydrophobic molecules may bind preferably to M5. We believe that the current study would provide a basis for the design of subsequent M4/5 selective antagonists.
Collapse
Affiliation(s)
- Lu Chen
- Key Laboratory of Structure-Based Drug Design & Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China.
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang 110016, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China
| | - Baichun Hu
- Key Laboratory of Structure-Based Drug Design & Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China.
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang 110016, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China
| | - Hanxun Wang
- Key Laboratory of Structure-Based Drug Design & Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China.
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang 110016, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China
| | - Weixia Li
- Key Laboratory of Structure-Based Drug Design & Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China.
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang 110016, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China
| | - Shizun Wang
- Key Laboratory of Structure-Based Drug Design & Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China.
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang 110016, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China
| | - Jiasi Luan
- Key Laboratory of Structure-Based Drug Design & Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China.
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang 110016, China
- School of Medical Devices, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China
| | - Haihan Liu
- Key Laboratory of Structure-Based Drug Design & Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China.
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang 110016, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China
| | - Jian Wang
- Key Laboratory of Structure-Based Drug Design & Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China.
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang 110016, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China
| | - Maosheng Cheng
- Key Laboratory of Structure-Based Drug Design & Discovery of Ministry of Education, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China.
- Key Laboratory of Intelligent Drug Design and New Drug Discovery of Liaoning Province, Shenyang Pharmaceutical University, Shenyang 110016, China
- School of Pharmaceutical Engineering, Shenyang Pharmaceutical University, Shenyang 110016, People's Republic of China
| |
Collapse
|
4
|
Dunkel A, Hofmann T, Di Pizio A. In Silico Investigation of Bitter Hop-Derived Compounds and Their Cognate Bitter Taste Receptors. JOURNAL OF AGRICULTURAL AND FOOD CHEMISTRY 2020; 68:10414-10423. [PMID: 32027492 DOI: 10.1021/acs.jafc.9b07863] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
The typical bitter taste of beer is caused by adding hops (Humulus lupulus L.) during the wort boiling process. The bitter taste of hop-derived compounds was found to be mediated by three bitter taste receptors: TAS2R1, TAS2R14, and TAS2R40. In this work, structural bioinformatics analyses were used to characterize the binding modes of trans-isocohumulone, trans-isohumulone, trans-isoadhumulone, cis-isocohumulone, cis-isohumulone, cis-isoadhumulone, cohumulone, humulone, adhumulone, and 8-prenylnaringenin into the orthosteric binding site of their cognate receptors. A conserved asparagine in transmembrane 3 was found to be essential for the recognition of hop-derived compounds, whereas the surrounding residues in the binding site of the three receptors encode the ligand specificity. Hop-derived compounds are renowned bioactive molecules and are considered as potential hit molecules for drug discovery to treat metabolic diseases. A chemoinformatics analysis revealed that hop-derived compounds cluster in a different region of the chemical space compared to known bitter food-derived compounds, pinpointing hop-derived compounds as a very peculiar class of bitter compounds.
Collapse
Affiliation(s)
- Andreas Dunkel
- Leibniz-Institute for Food Systems Biology at the Technical University of Munich, Lise-Meitner Straße 34, D-85354 Freising, Germany
| | - Thomas Hofmann
- Chair of Food Chemistry and Molecular Sensory Science, Technical University of Munich, Lise-Meitner-Straße 34, D-85354 Freising, Germany
| | - Antonella Di Pizio
- Leibniz-Institute for Food Systems Biology at the Technical University of Munich, Lise-Meitner Straße 34, D-85354 Freising, Germany
| |
Collapse
|
5
|
Cross JB. Methods for Virtual Screening of GPCR Targets: Approaches and Challenges. METHODS IN MOLECULAR BIOLOGY (CLIFTON, N.J.) 2017; 1705:233-264. [PMID: 29188566 DOI: 10.1007/978-1-4939-7465-8_11] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Virtual screening (VS) has become an integral part of the drug discovery process and is a valuable tool for finding novel chemical starting points for GPCR targets. Ligand-based VS makes use of biochemical data for known, active compounds and has been applied successfully to many diverse GPCRs. Recent progress in GPCR X-ray crystallography has made it possible to incorporate detailed structural information into the VS process. This chapter outlines the latest VS techniques along with examples that highlight successful applications of these methods. Best practices for increasing the likelihood of VS success, as well as ongoing challenges, are also discussed.
Collapse
Affiliation(s)
- Jason B Cross
- University of Texas MD Anderson Cancer Center, Houston, TX, 77054, USA.
| |
Collapse
|
6
|
Coudrat T, Simms J, Christopoulos A, Wootten D, Sexton PM. Improving virtual screening of G protein-coupled receptors via ligand-directed modeling. PLoS Comput Biol 2017; 13:e1005819. [PMID: 29131821 PMCID: PMC5708846 DOI: 10.1371/journal.pcbi.1005819] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2017] [Revised: 11/30/2017] [Accepted: 10/12/2017] [Indexed: 11/22/2022] Open
Abstract
G protein-coupled receptors (GPCRs) play crucial roles in cell physiology and pathophysiology. There is increasing interest in using structural information for virtual screening (VS) of libraries and for structure-based drug design to identify novel agonist or antagonist leads. However, the sparse availability of experimentally determined GPCR/ligand complex structures with diverse ligands impedes the application of structure-based drug design (SBDD) programs directed to identifying new molecules with a select pharmacology. In this study, we apply ligand-directed modeling (LDM) to available GPCR X-ray structures to improve VS performance and selectivity towards molecules of specific pharmacological profile. The described method refines a GPCR binding pocket conformation using a single known ligand for that GPCR. The LDM method is a computationally efficient, iterative workflow consisting of protein sampling and ligand docking. We developed an extensive benchmark comparing LDM-refined binding pockets to GPCR X-ray crystal structures across seven different GPCRs bound to a range of ligands of different chemotypes and pharmacological profiles. LDM-refined models showed improvement in VS performance over origin X-ray crystal structures in 21 out of 24 cases. In all cases, the LDM-refined models had superior performance in enriching for the chemotype of the refinement ligand. This likely contributes to the LDM success in all cases of inhibitor-bound to agonist-bound binding pocket refinement, a key task for GPCR SBDD programs. Indeed, agonist ligands are required for a plethora of GPCRs for therapeutic intervention, however GPCR X-ray structures are mostly restricted to their inactive inhibitor-bound state. G protein-coupled receptors (GPCRs) are a major target for drug discovery. These receptors are highly dynamic membrane proteins, and have had limited tractability using with biophysical screens that are widely adopted for globular protein targets. Thus, structure-based virtual screening (SBVS) holds great promise as a complement to physical screening for rational design of novel drugs. Indeed, the increasing number of atomic-detail GPCR X-ray crystal structures has coincided with an increase in prospective SBVS studies that have identified novel compounds. However, experimentally solved GPCR structures do not meet the full demand for SBVS, as the GPCR structural landscape is incomplete, lacking both in coverage of available GPCRs, and diversity in both receptor conformations and the chemistry of co-crystalised ligands. Here we present a novel computational GPCR binding pocket refinement method that can generate predictive GPCR/ligand complexes with improved SBVS performance. This ligand-directed modeling workflow uses parallel processing and efficient algorithms to search the GPCR/ligand conformational space faster and more efficiently than the widely used protein refinement method molecular dynamics. In this study, the resulting models are evaluated both structurally, and in retrospective SBVS. We demonstrate improved performance of refined models over their starting structures in the majority of our test cases.
Collapse
Affiliation(s)
- Thomas Coudrat
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - John Simms
- School of Life and Health Sciences, Aston University, Birmingham, United Kingdom
| | - Arthur Christopoulos
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
| | - Denise Wootten
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
- * E-mail: (DW); (PMS)
| | - Patrick M. Sexton
- Drug Discovery Biology and Department of Pharmacology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, Australia
- * E-mail: (DW); (PMS)
| |
Collapse
|
7
|
Fierro F, Suku E, Alfonso-Prieto M, Giorgetti A, Cichon S, Carloni P. Agonist Binding to Chemosensory Receptors: A Systematic Bioinformatics Analysis. Front Mol Biosci 2017; 4:63. [PMID: 28932739 PMCID: PMC5592726 DOI: 10.3389/fmolb.2017.00063] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2017] [Accepted: 08/22/2017] [Indexed: 12/17/2022] Open
Abstract
Human G-protein coupled receptors (hGPCRs) constitute a large and highly pharmaceutically relevant membrane receptor superfamily. About half of the hGPCRs' family members are chemosensory receptors, involved in bitter taste and olfaction, along with a variety of other physiological processes. Hence these receptors constitute promising targets for pharmaceutical intervention. Molecular modeling has been so far the most important tool to get insights on agonist binding and receptor activation. Here we investigate both aspects by bioinformatics-based predictions across all bitter taste and odorant receptors for which site-directed mutagenesis data are available. First, we observe that state-of-the-art homology modeling combined with previously used docking procedures turned out to reproduce only a limited fraction of ligand/receptor interactions inferred by experiments. This is most probably caused by the low sequence identity with available structural templates, which limits the accuracy of the protein model and in particular of the side-chains' orientations. Methods which transcend the limited sampling of the conformational space of docking may improve the predictions. As an example corroborating this, we review here multi-scale simulations from our lab and show that, for the three complexes studied so far, they significantly enhance the predictive power of the computational approach. Second, our bioinformatics analysis provides support to previous claims that several residues, including those at positions 1.50, 2.50, and 7.52, are involved in receptor activation.
Collapse
Affiliation(s)
- Fabrizio Fierro
- Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum JülichJülich, Germany
| | - Eda Suku
- Department of Biotechnology, University of VeronaVerona, Italy
| | - Mercedes Alfonso-Prieto
- Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum JülichJülich, Germany.,Cécile and Oskar Vogt Institute for Brain Research, Medical Faculty, Heinrich Heine University DüsseldorfDüsseldorf, Germany
| | - Alejandro Giorgetti
- Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum JülichJülich, Germany.,Department of Biotechnology, University of VeronaVerona, Italy
| | - Sven Cichon
- Institute of Neuroscience and Medicine INM-1, Forschungszentrum JülichJülich, Germany.,Institute for Human Genetics, Department of Genomics, Life&Brain Center, University of BonnBonn, Germany.,Division of Medical Genetics, Department of Biomedicine, University of BaselBasel, Switzerland
| | - Paolo Carloni
- Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum JülichJülich, Germany.,Department of Physics, Rheinisch-Westfälische Technische Hochschule AachenAachen, Germany.,VNU Key Laboratory "Multiscale Simulation of Complex Systems", VNU University of Science, Vietnam National UniversityHanoi, Vietnam
| |
Collapse
|
8
|
Ligand binding modes from low resolution GPCR models and mutagenesis: chicken bitter taste receptor as a test-case. Sci Rep 2017; 7:8223. [PMID: 28811548 PMCID: PMC5557796 DOI: 10.1038/s41598-017-08344-9] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2017] [Accepted: 07/10/2017] [Indexed: 11/08/2022] Open
Abstract
Bitter taste is one of the basic taste modalities, warning against consuming potential poisons. Bitter compounds activate members of the bitter taste receptor (Tas2r) subfamily of G protein-coupled receptors (GPCRs). The number of functional Tas2rs is species-dependent. Chickens represent an intriguing minimalistic model, because they detect the bitter taste of structurally different molecules with merely three bitter taste receptor subtypes. We investigated the binding modes of several known agonists of a representative chicken bitter taste receptor, ggTas2r1. Because of low sequence similarity between ggTas2r1 and crystallized GPCRs (~10% identity, ~30% similarity at most), the combination of computational approaches with site-directed mutagenesis was used to characterize the agonist-bound conformation of ggTas2r1 binding site between TMs 3, 5, 6 and 7. We found that the ligand interactions with N93 in TM3 and/or N247 in TM5, combined with hydrophobic contacts, are typically involved in agonist recognition. Next, the ggTas2r1 structural model was successfully used to identify three quinine analogues (epiquinidine, ethylhydrocupreine, quinidine) as new ggTas2r1 agonists. The integrated approach validated here may be applicable to additional cases where the sequence identity of the GPCR of interest and the existing experimental structures is low.
Collapse
|
9
|
Lupala CS, Rasaeifar B, Gomez-Gutierrez P, Perez JJ. Using molecular dynamics for the refinement of atomistic models of GPCRs by homology modeling. J Biomol Struct Dyn 2017; 36:2436-2448. [PMID: 28728517 DOI: 10.1080/07391102.2017.1357503] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
Despite GPCRs sharing a common seven helix bundle, analysis of the diverse crystallographic structures available reveal specific features that might be relevant for ligand design. Despite the number of crystallographic structures of GPCRs steadily increasing, there are still challenges that hamper the availability of new structures. In the absence of a crystallographic structure, homology modeling remains one of the important techniques for constructing 3D models of proteins. In the present study we investigated the use of molecular dynamics simulations for the refinement of GPCRs models constructed by homology modeling. Specifically, we investigated the relevance of template selection, ligand inclusion as well as the length of the simulation on the quality of the GPCRs models constructed. For this purpose we chose the crystallographic structure of the rat muscarinic M3 receptor as reference and constructed diverse atomistic models by homology modeling, using different templates. Specifically, templates used in the present work include the human muscarinic M2; the more distant human histamine H1 and the even more distant bovine rhodopsin as shown in the GPCRs phylogenetic tree. We also investigated the use or not of a ligand in the refinement process. Hence, we conducted the refinement process of the M3 model using the M2 muscarinic as template with tiotropium or NMS docked in the orthosteric site and compared with the results obtained with a model refined without any ligand bound.
Collapse
Affiliation(s)
- Cecylia S Lupala
- a Department of Chemical Engineering (ETSEIB) , Universitat Politecnica de Catalunya , Av. Diagonal, 647. 08028 Barcelona , Spain
| | - Bahareh Rasaeifar
- a Department of Chemical Engineering (ETSEIB) , Universitat Politecnica de Catalunya , Av. Diagonal, 647. 08028 Barcelona , Spain
| | - Patricia Gomez-Gutierrez
- a Department of Chemical Engineering (ETSEIB) , Universitat Politecnica de Catalunya , Av. Diagonal, 647. 08028 Barcelona , Spain
| | - Juan J Perez
- a Department of Chemical Engineering (ETSEIB) , Universitat Politecnica de Catalunya , Av. Diagonal, 647. 08028 Barcelona , Spain
| |
Collapse
|
10
|
Thomas T, Fang Y, Yuriev E, Chalmers DK. Ligand Binding Pathways of Clozapine and Haloperidol in the Dopamine D2 and D3 Receptors. J Chem Inf Model 2016; 56:308-21. [PMID: 26690887 DOI: 10.1021/acs.jcim.5b00457] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The binding of a small molecule ligand to its protein target is most often characterized by binding affinity and is typically viewed as an on/off switch. The more complex reality is that binding involves the ligand passing through a series of intermediate states between the solution phase and the fully bound pose. We have performed a set of 29 unbiased molecular dynamics simulations to model the binding pathways of the dopamine receptor antagonists clozapine and haloperidol binding to the D2 and D3 dopamine receptors. Through these simulations we have captured the binding pathways of clozapine and haloperidol from the extracellular vestibule to the orthosteric binding site and thereby, we also predict the bound pose of each ligand. These are the first long time scale simulations of haloperidol or clozapine binding to dopamine receptors. From these simulations, we have identified several important stages in the binding pathway, including the involvement of Tyr7.35 in a "handover" mechanism that transfers the ligand between the extracellular vestibule and Asp3.32. We have also performed interaction and cluster analyses to determine differences in binding pathways between the D2 and D3 receptors and identified metastable states that may be of use in drug design.
Collapse
Affiliation(s)
- Trayder Thomas
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University , 381 Royal Pde, Parkville, Victoria 3052, Australia
| | - Yu Fang
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University , 381 Royal Pde, Parkville, Victoria 3052, Australia
| | - Elizabeth Yuriev
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University , 381 Royal Pde, Parkville, Victoria 3052, Australia
| | - David K Chalmers
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University , 381 Royal Pde, Parkville, Victoria 3052, Australia
| |
Collapse
|
11
|
Thomas T, Chalmers DK, Yuriev E. Homology Modeling and Docking Evaluation of Human Muscarinic Acetylcholine Receptors. NEUROMETHODS 2016. [DOI: 10.1007/978-1-4939-2858-3_2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
12
|
Xu W, Lim J, Goh CY, Suen JY, Jiang Y, Yau MK, Wu KC, Liu L, Fairlie DP. Repurposing Registered Drugs as Antagonists for Protease-Activated Receptor 2. J Chem Inf Model 2015; 55:2079-84. [DOI: 10.1021/acs.jcim.5b00500] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Weijun Xu
- Division
of Chemistry and
Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Junxian Lim
- Division
of Chemistry and
Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Chai-Yeen Goh
- Division
of Chemistry and
Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Jacky Y. Suen
- Division
of Chemistry and
Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Yuhong Jiang
- Division
of Chemistry and
Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Mei-Kwan Yau
- Division
of Chemistry and
Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Kai-Chen Wu
- Division
of Chemistry and
Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Ligong Liu
- Division
of Chemistry and
Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - David P. Fairlie
- Division
of Chemistry and
Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Queensland 4072, Australia
| |
Collapse
|
13
|
Towards predictive docking at aminergic G-protein coupled receptors. J Mol Model 2015; 21:284. [PMID: 26453085 DOI: 10.1007/s00894-015-2824-9] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2015] [Accepted: 09/15/2015] [Indexed: 12/23/2022]
Abstract
G protein-coupled receptors (GPCRs) are hard to crystallize. However, attempts to predict their structure have boomed as a result of advancements in crystallographic techniques. This trend has allowed computer-aided molecular modeling of GPCRs. We analyzed the performance of four molecular modeling programs in pose evaluation of re-docked antagonists / inverse agonists to 11 original crystal structures of aminergic GPCRs using an induced fit-docking procedure. AutoDock and Glide were used for docking. AutoDock binding energy function, GlideXP, Prime MM-GB/SA, and YASARA binding function were used for pose scoring. Root mean square deviation (RMSD) of the best pose ranged from 0.09 to 1.58 Å, and median RMSD of the top 60 poses ranged from 1.47 to 3.83 Å. However, RMSD of the top pose ranged from 0.13 to 7.33 Å and ranking of the best pose ranged from the 1st to 60th out of 60 poses. Moreover, analysis of ligand-receptor interactions of top poses revealed substantial differences from interactions found in crystallographic structures. Bad ranking of top poses and discrepancies between top docked poses and crystal structures render current simple docking methods unsuitable for predictive modeling of receptor-ligand interactions. Prime MM-GB/SA optimized for 3NY9 by multiple linear regression did not work well at 3NY8 and 3NYA, structures of the same receptor with different ligands. However, 9 of 11 trajectories of molecular dynamics simulations by Desmond of top poses converged with trajectories of crystal structures. Key interactions were properly detected for all structures. This procedure also worked well for cross-docking of tested β2-adrenergic antagonists. Thus, this procedure represents a possible way to predict interactions of antagonists with aminergic GPCRs.
Collapse
|
14
|
Spyrakis F, Cavasotto CN. Open challenges in structure-based virtual screening: Receptor modeling, target flexibility consideration and active site water molecules description. Arch Biochem Biophys 2015; 583:105-19. [DOI: 10.1016/j.abb.2015.08.002] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2015] [Revised: 08/03/2015] [Accepted: 08/03/2015] [Indexed: 01/05/2023]
|
15
|
Yuriev E, Holien J, Ramsland PA. Improvements, trends, and new ideas in molecular docking: 2012-2013 in review. J Mol Recognit 2015; 28:581-604. [PMID: 25808539 DOI: 10.1002/jmr.2471] [Citation(s) in RCA: 159] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2014] [Revised: 01/16/2015] [Accepted: 02/05/2015] [Indexed: 12/11/2022]
Abstract
Molecular docking is a computational method for predicting the placement of ligands in the binding sites of their receptor(s). In this review, we discuss the methodological developments that occurred in the docking field in 2012 and 2013, with a particular focus on the more difficult aspects of this computational discipline. The main challenges and therefore focal points for developments in docking, covered in this review, are receptor flexibility, solvation, scoring, and virtual screening. We specifically deal with such aspects of molecular docking and its applications as selection criteria for constructing receptor ensembles, target dependence of scoring functions, integration of higher-level theory into scoring, implicit and explicit handling of solvation in the binding process, and comparison and evaluation of docking and scoring methods.
Collapse
Affiliation(s)
- Elizabeth Yuriev
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria, 3052, Australia
| | - Jessica Holien
- ACRF Rational Drug Discovery Centre and Structural Biology Laboratory, St. Vincent's Institute of Medical Research, Fitzroy, Victoria, 3065, Australia
| | - Paul A Ramsland
- Centre for Biomedical Research, Burnet Institute, Melbourne, Victoria, 3004, Australia.,Department of Surgery Austin Health, University of Melbourne, Melbourne, Victoria, 3084, Australia.,Department of Immunology, Monash University, Alfred Medical Research and Education Precinct, Melbourne, Victoria, 3004, Australia.,School of Biomedical Sciences, CHIRI Biosciences, Curtin University, Perth, Western Australia, 6845, Australia
| |
Collapse
|
16
|
Xu J, Yuan H, Ran T, Zhang Y, Liu H, Lu S, Xiong X, Xu A, Jiang Y, Lu T, Chen Y. A selectivity study of sodium-dependent glucose cotransporter 2/sodium-dependent glucose cotransporter 1 inhibitors by molecular modeling. J Mol Recognit 2015; 28:467-79. [DOI: 10.1002/jmr.2464] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2014] [Revised: 01/15/2015] [Accepted: 01/15/2015] [Indexed: 01/01/2023]
Affiliation(s)
- Jinxing Xu
- Laboratory of Molecular Design and Drug Discovery, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
| | - Haoliang Yuan
- Key Laboratory of Nuclear Medicine, Ministry of Health, Jiangsu Key Laboratory of Molecular Nuclear Medicine; Jiangsu Institute of Nuclear Medicine; Wuxi 214063 China
| | - Ting Ran
- Laboratory of Molecular Design and Drug Discovery, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
| | - Yanmin Zhang
- Laboratory of Molecular Design and Drug Discovery, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
| | - Haichun Liu
- Laboratory of Molecular Design and Drug Discovery, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
| | - Shuai Lu
- Laboratory of Molecular Design and Drug Discovery, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
| | - Xiao Xiong
- Laboratory of Molecular Design and Drug Discovery, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
| | - Anyang Xu
- Laboratory of Molecular Design and Drug Discovery, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
| | - Yulei Jiang
- Laboratory of Molecular Design and Drug Discovery, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
| | - Tao Lu
- Laboratory of Molecular Design and Drug Discovery, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
- State Key Laboratory of Natural Medicines, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
| | - Yadong Chen
- Laboratory of Molecular Design and Drug Discovery, School of Science; China Pharmaceutical University; 639 Longmian Avenue Nanjing 211198 China
| |
Collapse
|
17
|
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
|
18
|
Cavasotto CN, Palomba D. Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models. Chem Commun (Camb) 2015; 51:13576-94. [DOI: 10.1039/c5cc05050b] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.
Collapse
Affiliation(s)
- Claudio N. Cavasotto
- Instituto de Investigación en Biomedicina de Buenos Aires (IBioBA) - CONICET - Partner Institute of the Max Planck Society
- Buenos Aires
- Argentina
| | - Damián Palomba
- Instituto de Investigación en Biomedicina de Buenos Aires (IBioBA) - CONICET - Partner Institute of the Max Planck Society
- Buenos Aires
- Argentina
| |
Collapse
|
19
|
Chin SP, Buckle MJC, Chalmers DK, Yuriev E, Doughty SW. Toward activated homology models of the human M1 muscarinic acetylcholine receptor. J Mol Graph Model 2014; 49:91-8. [PMID: 24631873 DOI: 10.1016/j.jmgm.2014.02.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2013] [Revised: 02/14/2014] [Accepted: 02/17/2014] [Indexed: 10/25/2022]
Abstract
Structure-based virtual screening offers a good opportunity for the discovery of selective M1 muscarinic acetylcholine receptor (mAChR) agonists for the treatment of Alzheimer's disease. However, no 3-D structure of an M1 mAChR is yet available and the homology models that have been previously reported are only able to identify antagonists in virtual screening experiments. In this study, we generated a homology model of the human M1 mAChR, based on the crystal structure of an M3 mAChR as the template. This initial model was modified, using the agonist-bound crystal structure of a β2-adrenergic receptor as a guide, to give two possible activated structures. The T192 side chain was adjusted in both structures and one of the structures also had the whole of transmembrane (TM) 5 rotated and tilted toward the inner channel of the transmembrane region. The binding sites of all three structures were then refined by induced-fit docking (IFD) with acetylcholine. Virtual screening experiments showed that all three refined models could efficiently differentiate agonists from decoy molecules, with the TM5-modified models also giving good agonist/antagonist selectivity. The whole range of agonists and antagonists was observed to bind within the orthosteric site of the structure obtained by IFD refinement alone, implying that it has inactive state character. In contrast, the two TM5-modified structures were unable to accommodate the antagonists, supporting the proposition that they possess activated state character.
Collapse
Affiliation(s)
- Sek Peng Chin
- Department of Pharmacy, Faculty of Medicine, University of Malaya, 50603 Kuala Lumpur, Malaysia
| | - Michael J C Buckle
- Department of Pharmacy, Faculty of Medicine, University of Malaya, 50603 Kuala Lumpur, Malaysia.
| | - David K Chalmers
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University (Parkville Campus), 381 Royal Parade, Parkville, Victoria 3052, Australia
| | - Elizabeth Yuriev
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University (Parkville Campus), 381 Royal Parade, Parkville, Victoria 3052, Australia
| | - Stephen W Doughty
- The School of Pharmacy, University of Nottingham Malaysia Campus, Jalan Broga, 43500 Semenyih, Selangor, Malaysia.
| |
Collapse
|