1
|
Qian R, Xue J, Xu Y, Huang J. Alchemical Transformations and Beyond: Recent Advances and Real-World Applications of Free Energy Calculations in Drug Discovery. J Chem Inf Model 2024; 64:7214-7237. [PMID: 39360948 DOI: 10.1021/acs.jcim.4c01024] [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: 10/15/2024]
Abstract
Computational methods constitute efficient strategies for screening and optimizing potential drug molecules. A critical factor in this process is the binding affinity between candidate molecules and targets, quantified as binding free energy. Among various estimation methods, alchemical transformation methods stand out for their theoretical rigor. Despite challenges in force field accuracy and sampling efficiency, advancements in algorithms, software, and hardware have increased the application of free energy perturbation (FEP) calculations in the pharmaceutical industry. Here, we review the practical applications of FEP in drug discovery projects since 2018, covering both ligand-centric and residue-centric transformations. We show that relative binding free energy calculations have steadily achieved chemical accuracy in real-world applications. In addition, we discuss alternative physics-based simulation methods and the incorporation of deep learning into free energy calculations.
Collapse
Affiliation(s)
- Runtong Qian
- Westlake AI Therapeutics Lab, Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, Zhejiang 310024, China
| | - Jing Xue
- Westlake AI Therapeutics Lab, Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, Zhejiang 310024, China
| | - You Xu
- Westlake AI Therapeutics Lab, Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, Zhejiang 310024, China
| | - Jing Huang
- Westlake AI Therapeutics Lab, Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, Zhejiang 310024, China
| |
Collapse
|
2
|
Gu X, Aranganathan A, Tiwary P. Empowering AlphaFold2 for protein conformation selective drug discovery with AlphaFold2-RAVE. eLife 2024; 13:RP99702. [PMID: 39240197 PMCID: PMC11379456 DOI: 10.7554/elife.99702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2024] Open
Abstract
Small-molecule drug design hinges on obtaining co-crystallized ligand-protein structures. Despite AlphaFold2's strides in protein native structure prediction, its focus on apo structures overlooks ligands and associated holo structures. Moreover, designing selective drugs often benefits from the targeting of diverse metastable conformations. Therefore, direct application of AlphaFold2 models in virtual screening and drug discovery remains tentative. Here, we demonstrate an AlphaFold2-based framework combined with all-atom enhanced sampling molecular dynamics and Induced Fit docking, named AF2RAVE-Glide, to conduct computational model-based small-molecule binding of metastable protein kinase conformations, initiated from protein sequences. We demonstrate the AF2RAVE-Glide workflow on three different mammalian protein kinases and their type I and II inhibitors, with special emphasis on binding of known type II kinase inhibitors which target the metastable classical DFG-out state. These states are not easy to sample from AlphaFold2. Here, we demonstrate how with AF2RAVE these metastable conformations can be sampled for different kinases with high enough accuracy to enable subsequent docking of known type II kinase inhibitors with more than 50% success rates across docking calculations. We believe the protocol should be deployable for other kinases and more proteins generally.
Collapse
Affiliation(s)
- Xinyu Gu
- Institute for Physical Science and Technology, University of Maryland, College Park, United States
- University of Maryland Institute for Health Computing, Bethesda, United States
| | - Akashnathan Aranganathan
- Institute for Physical Science and Technology, University of Maryland, College Park, United States
- Biophysics Program, University of Maryland, College Park, United States
| | - Pratyush Tiwary
- Institute for Physical Science and Technology, University of Maryland, College Park, United States
- University of Maryland Institute for Health Computing, Bethesda, United States
- Department of Chemistry and Biochemistry, University of Maryland, College Park, United States
| |
Collapse
|
3
|
King JE, Koes DR. Interpreting forces as deep learning gradients improves quality of predicted protein structures. Biophys J 2024; 123:2730-2739. [PMID: 38104241 PMCID: PMC11393680 DOI: 10.1016/j.bpj.2023.12.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Revised: 10/30/2023] [Accepted: 12/12/2023] [Indexed: 12/19/2023] Open
Abstract
Protein structure predictions from deep learning models like AlphaFold2, despite their remarkable accuracy, are likely insufficient for direct use in downstream tasks like molecular docking. The functionality of such models could be improved with a combination of increased accuracy and physical intuition. We propose a new method to train deep learning protein structure prediction models using molecular dynamics force fields to work toward these goals. Our custom PyTorch loss function, OpenMM-Loss, represents the potential energy of a predicted structure. OpenMM-Loss can be applied to any all-atom representation of a protein structure capable of mapping into our software package, SidechainNet. We demonstrate our method's efficacy by finetuning OpenFold. We show that subsequently predicted protein structures, both before and after a relaxation procedure, exhibit comparable accuracy while displaying lower potential energy and improved structural quality as assessed by MolProbity metrics.
Collapse
Affiliation(s)
- Jonathan Edward King
- Joint PhD Program in Computational Biology, Carnegie Mellon University-University of Pittsburgh, Pittsburgh, Pennsylvania
| | - David Ryan Koes
- Computational & Systems Biology, University of Pittsburgh, Pittsburgh, Pennsylvania.
| |
Collapse
|
4
|
Pala D, Clark DE. Caught between a ROCK and a hard place: current challenges in structure-based drug design. Drug Discov Today 2024; 29:104106. [PMID: 39029868 DOI: 10.1016/j.drudis.2024.104106] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2024] [Revised: 06/27/2024] [Accepted: 07/13/2024] [Indexed: 07/21/2024]
Abstract
The discipline of structure-based drug design (SBDD) is several decades old and it is tempting to think that the proliferation of experimental structures for many drug targets might make computer-aided drug design (CADD) straightforward. However, this is far from true. In this review, we illustrate some of the challenges that CADD scientists face every day in their work, even now. We use Rho-associated protein kinase (ROCK), and public domain structures and data, as an example to illustrate some of the challenges we have experienced during our project targeting this protein. We hope that this will help to prevent unrealistic expectations of what CADD can accomplish and to educate non-CADD scientists regarding the challenges still facing their CADD colleagues.
Collapse
Affiliation(s)
- Daniele Pala
- Medicinal Chemistry and Drug Design Technologies Department, Chiesi Farmaceutici S.p.A, Research Center, Largo Belloli 11/a, 43122 Parma, Italy
| | - David E Clark
- Charles River, 6-9 Spire Green Centre, Flex Meadow, Harlow CM19 5TR, UK.
| |
Collapse
|
5
|
Lyu J, Kapolka N, Gumpper R, Alon A, Wang L, Jain MK, Barros-Álvarez X, Sakamoto K, Kim Y, DiBerto J, Kim K, Glenn IS, Tummino TA, Huang S, Irwin JJ, Tarkhanova OO, Moroz Y, Skiniotis G, Kruse AC, Shoichet BK, Roth BL. AlphaFold2 structures guide prospective ligand discovery. Science 2024; 384:eadn6354. [PMID: 38753765 PMCID: PMC11253030 DOI: 10.1126/science.adn6354] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2023] [Accepted: 04/24/2024] [Indexed: 05/18/2024]
Abstract
AlphaFold2 (AF2) models have had wide impact but mixed success in retrospective ligand recognition. We prospectively docked large libraries against unrefined AF2 models of the σ2 and serotonin 2A (5-HT2A) receptors, testing hundreds of new molecules and comparing results with those obtained from docking against the experimental structures. Hit rates were high and similar for the experimental and AF2 structures, as were affinities. Success in docking against the AF2 models was achieved despite differences between orthosteric residue conformations in the AF2 models and the experimental structures. Determination of the cryo-electron microscopy structure for one of the more potent 5-HT2A ligands from the AF2 docking revealed residue accommodations that resembled the AF2 prediction. AF2 models may sample conformations that differ from experimental structures but remain low energy and relevant for ligand discovery, extending the domain of structure-based drug design.
Collapse
Affiliation(s)
- Jiankun Lyu
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
- The Evnin Family Laboratory of Computational Molecular Discovery, The Rockefeller University, New York, NY 10065, USA
| | - Nicholas Kapolka
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599, USA
| | - Ryan Gumpper
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599, USA
| | - Assaf Alon
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA
| | - Liang Wang
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94035, USA
| | - Manish K. Jain
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599, USA
| | - Ximena Barros-Álvarez
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94035, USA
| | - Kensuke Sakamoto
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599, USA
- National Institute of Mental Health Psychoactive Drug Screening Program (NIMH PDSP), School of Medicine, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599, USA
| | - Yoojoong Kim
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599, USA
| | - Jeffrey DiBerto
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599, USA
| | - Kuglae Kim
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599, USA
| | - Isabella S. Glenn
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
| | - Tia A. Tummino
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
| | - Sijie Huang
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
| | - John J. Irwin
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
| | | | - Yurii Moroz
- Chemspace LLC, Kyiv 02094, Ukraine
- Taras Shevchenko National University of Kyiv, Kyiv 01601, Ukraine
- Enamine Ltd., Kyiv 02094, Ukraine
| | - Georgios Skiniotis
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94035, USA
- Department of Structural Biology, Stanford University School of Medicine, Stanford, CA 94304, USA
| | - Andrew C. Kruse
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA
| | - Brian K. Shoichet
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
| | - Bryan L. Roth
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599, USA
- National Institute of Mental Health Psychoactive Drug Screening Program (NIMH PDSP), School of Medicine, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599, USA
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| |
Collapse
|
6
|
Gu S, Yang Y, Zhao Y, Qiu J, Wang X, Tong HHY, Liu L, Wan X, Liu H, Hou T, Kang Y. Evaluation of AlphaFold2 Structures for Hit Identification across Multiple Scenarios. J Chem Inf Model 2024; 64:3630-3639. [PMID: 38630855 DOI: 10.1021/acs.jcim.3c01976] [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: 04/19/2024]
Abstract
The introduction of AlphaFold2 (AF2) has sparked significant enthusiasm and generated extensive discussion within the scientific community, particularly among drug discovery researchers. Although previous studies have addressed the performance of AF2 structures in virtual screening (VS), a more comprehensive investigation is still necessary considering the paramount importance of structural accuracy in drug design. In this study, we evaluate the performance of AF2 structures in VS across three common drug discovery scenarios: targets with holo, apo, and AF2 structures; targets with only apo and AF2 structures; and targets exclusively with AF2 structures. We utilized both the traditional physics-based Glide and the deep-learning-based scoring function RTMscore to rank the compounds in the DUD-E, DEKOIS 2.0, and DECOY data sets. The results demonstrate that, overall, the performance of VS on AF2 structures is comparable to that on apo structures but notably inferior to that on holo structures across diverse scenarios. Moreover, when a target has solely AF2 structure, selecting the holo structure of the target from different subtypes within the same protein family produces comparable results with the AF2 structure for VS on the data set of the AF2 structures, and significantly better results than the AF2 structures on its own data set. This indicates that utilizing AF2 structures for docking-based VS may not yield most satisfactory outcomes, even when solely AF2 structures are available. Moreover, we rule out the possibility that the variations in VS performance between the binding pockets of AF2 and holo structures arise from the differences in their biological assembly composition.
Collapse
Affiliation(s)
- Shukai Gu
- Faculty of Applied Science, Macao Polytechnic University, Macao 999078, SAR, China
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Yuwei Yang
- Faculty of Applied Science, Macao Polytechnic University, Macao 999078, SAR, China
| | - Yihao Zhao
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Jiayue Qiu
- Faculty of Applied Science, Macao Polytechnic University, Macao 999078, SAR, China
| | - Xiaorui Wang
- State Key Laboratory of Quality Re-search in Chinese Medicine, Macau Institute for Applied Research in Medicine and Health, Macau University of Science and Technology, Macao 999078, China
| | - Henry Hoi Yee Tong
- Faculty of Applied Science, Macao Polytechnic University, Macao 999078, SAR, China
| | - Liwei Liu
- Advanced Computing and Storage Laboratory, Central Research Institute, 2012 Laboratories, Huawei Technologies Co., Ltd., Nanjing 210000, Jiangsu, China
| | - Xiaozhe Wan
- Advanced Computing and Storage Laboratory, Central Research Institute, 2012 Laboratories, Huawei Technologies Co., Ltd., Nanjing 210000, Jiangsu, China
| | - Huanxiang Liu
- Faculty of Applied Science, Macao Polytechnic University, Macao 999078, SAR, China
| | - Tingjun Hou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Yu Kang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| |
Collapse
|
7
|
Cournia Z, Chipot C. Applications of Free-Energy Calculations to Biomolecular Processes. A Collection. J Phys Chem B 2024; 128:3299-3301. [PMID: 38600851 DOI: 10.1021/acs.jpcb.4c01283] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2024]
Affiliation(s)
- Zoe Cournia
- Biomedical Research Foundation, Academy of Athens, 11527 Athens, Greece
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n◦7019, Université de Lorraine, B.P. 70239, 54506 Vandœuvre-lès-Nancy Cedex, France
- Department of Biochemistry and Molecular Biology, The University of Chicago, 929 East 57th Street W225, Chicago, Illinois 60637, United States
- Theoretical and Computational Biophysics Group, Beckman Institute, and Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
8
|
Cournia Z, Chipot C. Applications of Free-Energy Calculations to Biomolecular Processes. A Collection. J Chem Inf Model 2024; 64:2129-2131. [PMID: 38587007 DOI: 10.1021/acs.jcim.4c00349] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/09/2024]
Affiliation(s)
- Zoe Cournia
- Biomedical Research Foundation, Academy of Athens, 11527 Athens, Greece
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n◦7019, Université de Lorraine, B.P. 70239, 54506 Vandœuvre-lès-Nancy Cedex, France
- Department of Biochemistry and Molecular Biology, The University of Chicago, 929 East 57th Street W225, Chicago, Illinois 60637, United States
- Theoretical and Computational Biophysics Group, Beckman Institute, and Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
9
|
Lyu J, Kapolka N, Gumpper R, Alon A, Wang L, Jain MK, Barros-Álvarez X, Sakamoto K, Kim Y, DiBerto J, Kim K, Tummino TA, Huang S, Irwin JJ, Tarkhanova OO, Moroz Y, Skiniotis G, Kruse AC, Shoichet BK, Roth BL. AlphaFold2 structures template ligand discovery. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.12.20.572662. [PMID: 38187536 PMCID: PMC10769324 DOI: 10.1101/2023.12.20.572662] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/09/2024]
Abstract
AlphaFold2 (AF2) and RosettaFold have greatly expanded the number of structures available for structure-based ligand discovery, even though retrospective studies have cast doubt on their direct usefulness for that goal. Here, we tested unrefined AF2 models prospectively, comparing experimental hit-rates and affinities from large library docking against AF2 models vs the same screens targeting experimental structures of the same receptors. In retrospective docking screens against the σ2 and the 5-HT2A receptors, the AF2 structures struggled to recapitulate ligands that we had previously found docking against the receptors' experimental structures, consistent with published results. Prospective large library docking against the AF2 models, however, yielded similar hit rates for both receptors versus docking against experimentally-derived structures; hundreds of molecules were prioritized and tested against each model and each structure of each receptor. The success of the AF2 models was achieved despite differences in orthosteric pocket residue conformations for both targets versus the experimental structures. Intriguingly, against the 5-HT2A receptor the most potent, subtype-selective agonists were discovered via docking against the AF2 model, not the experimental structure. To understand this from a molecular perspective, a cryoEM structure was determined for one of the more potent and selective ligands to emerge from docking against the AF2 model of the 5-HT2A receptor. Our findings suggest that AF2 models may sample conformations that are relevant for ligand discovery, much extending the domain of applicability of structure-based ligand discovery.
Collapse
Affiliation(s)
- Jiankun Lyu
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
- The Evnin Family Laboratory of Computational Molecular Discovery, The Rockefeller University, New York, NY 10065, USA (present address)
| | - Nicholas Kapolka
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599-7365, USA
| | - Ryan Gumpper
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599-7365, USA
| | - Assaf Alon
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, MA, USA
- Pharmacology Department, Yale School of Medicine, New Haven, CT 06510, USA (present address)
| | - Liang Wang
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA, USA
| | - Manish K Jain
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599-7365, USA
| | - Ximena Barros-Álvarez
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA, USA
| | - Kensuke Sakamoto
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599-7365, USA
- National Institute of Mental Health Psychoactive Drug Screening Program (NIMH PDSP), School of Medicine, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599-7365, USA
| | - Yoojoong Kim
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599-7365, USA
| | - Jeffrey DiBerto
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599-7365, USA
| | - Kuglae Kim
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599-7365, USA
- Department of Pharmacy, College of Pharmacy, Yonsei University, Incheon 21983, Korea (present address)
| | - Tia A Tummino
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
| | - Sijie Huang
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
| | - John J Irwin
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
| | | | - Yurii Moroz
- Chemspace LLC, Kyiv, 02094, Ukraine
- Taras Shevchenko National University of Kyiv, Kyiv, 01601, Ukraine
- Enamine Ltd., Kyiv, 02094, Ukraine
| | - Georgios Skiniotis
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA, USA
- Department of Structural Biology, Stanford University School of Medicine, Stanford, CA, US
| | - Andrew C Kruse
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, MA, USA
| | - Brian K Shoichet
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158, USA
| | - Bryan L Roth
- Department of Pharmacology, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599-7365, USA
- National Institute of Mental Health Psychoactive Drug Screening Program (NIMH PDSP), School of Medicine, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27599-7365, USA
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599-7360, USA
| |
Collapse
|
10
|
Schapira M, Halabelian L, Arrowsmith CH, Harding RJ. Big data and benchmarking initiatives to bridge the gap from AlphaFold to drug design. Nat Chem Biol 2024:10.1038/s41589-024-01570-z. [PMID: 38459278 DOI: 10.1038/s41589-024-01570-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/10/2024]
Affiliation(s)
- Matthieu Schapira
- Structural Genomics Consortium, University of Toronto, Toronto, Ontario, Canada
- Department of Pharmacology & Toxicology, University of Toronto, Toronto, Ontario, Canada
| | - Levon Halabelian
- Structural Genomics Consortium, University of Toronto, Toronto, Ontario, Canada
- Department of Pharmacology & Toxicology, University of Toronto, Toronto, Ontario, Canada
| | - Cheryl H Arrowsmith
- Structural Genomics Consortium, University of Toronto, Toronto, Ontario, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Princess Margaret Cancer Centre, Toronto, Ontario, Canada
| | - Rachel J Harding
- Structural Genomics Consortium, University of Toronto, Toronto, Ontario, Canada.
- Department of Pharmacology & Toxicology, University of Toronto, Toronto, Ontario, Canada.
| |
Collapse
|
11
|
Yang Y, Cheng Y, Bai T, Liu S, Du Q, Xia W, Liu Y, Wang X, Chen X. Optimizing Trilobatin Production via Screening and Modification of Glycosyltransferases. Molecules 2024; 29:643. [PMID: 38338387 PMCID: PMC10856287 DOI: 10.3390/molecules29030643] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2023] [Revised: 01/26/2024] [Accepted: 01/26/2024] [Indexed: 02/12/2024] Open
Abstract
Trilobatin (TBL) is a key sweet compound from the traditional Chinese sweet tea plant (Rubus suavissimus S. Lee). Because of its intense sweetness, superior taste profile, and minimal caloric value, it serves as an exemplary natural dihydrochalcone sweetener. It also has various health benefits, including anti-inflammatory and glucose-lowering effects. It is primarily produced through botanical extraction, which impedes its scalability and cost-effectiveness. In a novel biotechnological approach, phloretin is used as a precursor that is transformed into TBL by the glycosyltransferase enzyme ph-4'-OGT. However, this enzyme's low catalytic efficiency and by-product formation limit the large-scale synthesis of TBL. In our study, the enzyme Mdph-4'-OGT was used to screen 17 sequences across species for TBL synthesis, of which seven exhibited catalytic activity. Notably, PT577 exhibited an unparalleled 97.3% conversion yield within 3 h. We then optimized the reaction conditions of PT577, attaining a peak TBL bioproduction of 163.3 mg/L. By employing virtual screening, we identified 25 mutation sites for PT577, thereby creating mutant strains that reduced by-products by up to 50%. This research enhances the enzymatic precision for TBL biosynthesis and offers a robust foundation for its industrial-scale production, with broader implications for the engineering and in silico analysis of glycosyltransferases.
Collapse
Affiliation(s)
- Yue Yang
- Jiaxing Synbiolab Biotechnology Co., Ltd., Jiaxing 314006, China; (Y.Y.); (T.B.); (S.L.); (Q.D.)
- School of Ecology and Environment, Northwestern Polytechnical University, Xi’an 710072, China; (Y.C.); (W.X.); (Y.L.)
| | - Yuhan Cheng
- School of Ecology and Environment, Northwestern Polytechnical University, Xi’an 710072, China; (Y.C.); (W.X.); (Y.L.)
| | - Tao Bai
- Jiaxing Synbiolab Biotechnology Co., Ltd., Jiaxing 314006, China; (Y.Y.); (T.B.); (S.L.); (Q.D.)
| | - Shimeng Liu
- Jiaxing Synbiolab Biotechnology Co., Ltd., Jiaxing 314006, China; (Y.Y.); (T.B.); (S.L.); (Q.D.)
| | - Qiuhui Du
- Jiaxing Synbiolab Biotechnology Co., Ltd., Jiaxing 314006, China; (Y.Y.); (T.B.); (S.L.); (Q.D.)
| | - Wenhao Xia
- School of Ecology and Environment, Northwestern Polytechnical University, Xi’an 710072, China; (Y.C.); (W.X.); (Y.L.)
| | - Yi Liu
- School of Ecology and Environment, Northwestern Polytechnical University, Xi’an 710072, China; (Y.C.); (W.X.); (Y.L.)
| | - Xiao Wang
- Jiaxing Synbiolab Biotechnology Co., Ltd., Jiaxing 314006, China; (Y.Y.); (T.B.); (S.L.); (Q.D.)
| | - Xianqing Chen
- Jiaxing Synbiolab Biotechnology Co., Ltd., Jiaxing 314006, China; (Y.Y.); (T.B.); (S.L.); (Q.D.)
| |
Collapse
|
12
|
Coskun D, Lihan M, Rodrigues JPGLM, Vass M, Robinson D, Friesner RA, Miller EB. Using AlphaFold and Experimental Structures for the Prediction of the Structure and Binding Affinities of GPCR Complexes via Induced Fit Docking and Free Energy Perturbation. J Chem Theory Comput 2024; 20:477-489. [PMID: 38100422 DOI: 10.1021/acs.jctc.3c00839] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2023]
Abstract
Free energy perturbation (FEP) remains an indispensable method for computationally assaying prospective compounds in advance of synthesis. However, before FEP can be deployed prospectively, it must demonstrate retrospective recapitulation of known experimental data where the subtle details of the atomic ligand-receptor model are consequential. An open question is whether AlphaFold models can serve as useful initial models for FEP in the regime where there exists a congeneric series of known chemical matter but where no experimental structures are available either of the target or of close homologues. As AlphaFold structures are provided without a bound ligand, we employ induced fit docking to refine the AlphaFold models in the presence of one or more congeneric ligands. In this work, we first validate the performance of our latest induced fit docking technology, IFD-MD, on a retrospective set of public experimental GPCR structures with 95% of cross-docks producing a pose with a ligand RMSD ≤ 2.5 Å in the top two predictions. We then apply IFD-MD and FEP on AlphaFold models of the somatostatin receptor family of GPCRs. We use AlphaFold models produced prior to the availability of any experimental structure from this family. We arrive at FEP-validated models for SSTR2, SSTR4, and SSTR5, with RMSE around 1 kcal/mol, and explore the challenges of model validation under scenarios of limited ligand data, ample ligand data, and categorical data.
Collapse
Affiliation(s)
- Dilek Coskun
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Muyun Lihan
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | | | - Márton Vass
- Schrödinger Technologies Limited, Davidson House, First Floor, Reading RG1 3 EU, U.K
| | - Daniel Robinson
- Schrödinger Technologies Limited, Davidson House, First Floor, Reading RG1 3 EU, U.K
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3110, New York, New York 10036, United States
| | - Edward B Miller
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| |
Collapse
|
13
|
Ahmed M, Maldonado AM, Durrant JD. From byte to bench to bedside: molecular dynamics simulations and drug discovery. BMC Biol 2023; 21:299. [PMID: 38155355 PMCID: PMC10755930 DOI: 10.1186/s12915-023-01791-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2023] [Accepted: 12/01/2023] [Indexed: 12/30/2023] Open
Affiliation(s)
- Mayar Ahmed
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, PA, USA
| | - Alex M Maldonado
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, PA, USA
| | - Jacob D Durrant
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, PA, USA.
| |
Collapse
|
14
|
Haufe Y, Kuruva V, Samanani Z, Lokaj G, Kamnesky G, Shadamarshan P, Shahoei R, Katz D, Sampson JM, Pusch M, Brik A, Nicke A, Leffler AE. Basic Residues at Position 11 of α-Conotoxin LvIA Influence Subtype Selectivity between α3β2 and α3β4 Nicotinic Receptors via an Electrostatic Mechanism. ACS Chem Neurosci 2023; 14:4311-4322. [PMID: 38051211 DOI: 10.1021/acschemneuro.3c00506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2023] Open
Abstract
Understanding the determinants of α-conotoxin (α-CTX) selectivity for different nicotinic acetylcholine receptor (nAChR) subtypes is a prerequisite for the design of tool compounds to study nAChRs. However, selectivity optimization of these small, disulfide-rich peptides is difficult not only because of an absence of α-CTX/nAChR co-structures but also because it is challenging to predict how a mutation to an α-CTX will alter its potency and selectivity. As a prototypical system to investigate selectivity, we employed the α-CTX LvIA that is 25-fold selective for the α3β2 nAChR over the related α3β4 nAChR subtype, which is a target for nicotine addiction. Using two-electrode voltage clamp electrophysiology, we identified LvIA[D11R] that is 2-fold selective for the α3β4 nAChR, reversing the subtype preference. This effect is specifically due to the change in charge and not shape of LvIA[D11R], as substitution of D11 with citrulline retains selectivity for the α3β2 nAChR. Furthermore, LvIA[D11K] shows a stronger reversal, with 4-fold selectivity for the α3β4 nAChR. Motivated by these findings, using site-directed mutagenesis, we found that β2[K79A] (I79 on β4), but not β2[K78A] (N78 on β4), largely restores the potency of basic mutants at position 11. Finally, to understand the structural basis of this effect, we used AlphaFold2 to generate models of LvIA in complex with both nAChR subtypes. Both models confirm the plausibility of an electrostatic mechanism to explain the data and also reproduce a broad range of potency and selectivity structure-activity relationships for LvIA mutants, as measured using free energy perturbation simulations. Our work highlights how electrostatic interactions can drive α-CTX selectivity and may serve as a strategy for optimizing the selectivity of LvIA and other α-CTXs.
Collapse
Affiliation(s)
- Yves Haufe
- Faculty of Medicine, Walther Straub Institute of Pharmacology and Toxicology, Ludwig-Maximilians-Universität München, Munich 80539, Germany
| | - Veeresh Kuruva
- Schulich Faculty of Chemistry, Technion-Israel Institute of Technology, Haifa 3200008, Israel
| | - Ziyana Samanani
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| | - Gonxhe Lokaj
- Faculty of Medicine, Walther Straub Institute of Pharmacology and Toxicology, Ludwig-Maximilians-Universität München, Munich 80539, Germany
| | - Guy Kamnesky
- Schulich Faculty of Chemistry, Technion-Israel Institute of Technology, Haifa 3200008, Israel
| | - PranavKumar Shadamarshan
- Faculty of Medicine, Walther Straub Institute of Pharmacology and Toxicology, Ludwig-Maximilians-Universität München, Munich 80539, Germany
| | - Rezvan Shahoei
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| | - Dana Katz
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| | - Jared M Sampson
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| | - Michael Pusch
- Istituto di Biofisica, CNR, Via De Marini 6, Genova 16149, Italy
| | - Ashraf Brik
- Schulich Faculty of Chemistry, Technion-Israel Institute of Technology, Haifa 3200008, Israel
| | - Annette Nicke
- Faculty of Medicine, Walther Straub Institute of Pharmacology and Toxicology, Ludwig-Maximilians-Universität München, Munich 80539, Germany
| | - Abba E Leffler
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| |
Collapse
|
15
|
Binbay FA, Rathod DC, George AAP, Imhof D. Quality Assessment of Selected Protein Structures Derived from Homology Modeling and AlphaFold. Pharmaceuticals (Basel) 2023; 16:1662. [PMID: 38139789 PMCID: PMC10747200 DOI: 10.3390/ph16121662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Revised: 11/21/2023] [Accepted: 11/24/2023] [Indexed: 12/24/2023] Open
Abstract
With technology advancing, many prediction algorithms have been developed to facilitate the modeling of inherently dynamic and flexible macromolecules such as proteins. Improvements in the prediction of protein structures have attracted a great deal of attention due to the advantages they offer, e.g., in drug design. While trusted experimental methods, such as X-ray crystallography, NMR spectroscopy, and electron microscopy, are preferred structure analysis techniques, in silico approaches are also being widely used. Two computational methods, which are on opposite ends of the spectrum with respect to their modus operandi, i.e., homology modeling and AlphaFold, have been established to provide high-quality structures. Here, a comparative study of the quality of structures either predicted by homology modeling or by AlphaFold is presented based on the characteristics determined by experimental studies using structure validation servers to fulfill the purpose. Although AlphaFold is able to predict high-quality structures, high-confidence parts are sometimes observed to be in disagreement with experimental data. On the other hand, while the structures obtained from homology modeling are successful in incorporating all aspects of the experimental structure used as a template, this method may struggle to accurately model a structure in the absence of a suitable template. In general, although both methods produce high-quality models, the criteria by which they are superior to each other are different and thus discussed in detail.
Collapse
Affiliation(s)
- Furkan Ayberk Binbay
- Pharmaceutical Biochemistry and Bioanalytics, Pharmaceutical Institute, University of Bonn, An der Immenburg 4, 53121 Bonn, Germany
| | - Dhruv Chetanbhai Rathod
- Pharmaceutical Biochemistry and Bioanalytics, Pharmaceutical Institute, University of Bonn, An der Immenburg 4, 53121 Bonn, Germany
| | | | - Diana Imhof
- Pharmaceutical Biochemistry and Bioanalytics, Pharmaceutical Institute, University of Bonn, An der Immenburg 4, 53121 Bonn, Germany
| |
Collapse
|
16
|
Nunes-Alves A, Merz K. AlphaFold2 in Molecular Discovery. J Chem Inf Model 2023; 63:5947-5949. [PMID: 37807755 DOI: 10.1021/acs.jcim.3c01459] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/10/2023]
Affiliation(s)
- Ariane Nunes-Alves
- Institute of Chemistry, Technische Universität Berlin, Berlin 10623, Germany
| | - Kenneth Merz
- Department of Chemistry, Michigan State University, East Lansing 48824, Michigan, United States
| |
Collapse
|
17
|
Stuart DD, Guzman-Perez A, Brooijmans N, Jackson EL, Kryukov GV, Friedman AA, Hoos A. Precision Oncology Comes of Age: Designing Best-in-Class Small Molecules by Integrating Two Decades of Advances in Chemistry, Target Biology, and Data Science. Cancer Discov 2023; 13:2131-2149. [PMID: 37712571 PMCID: PMC10551669 DOI: 10.1158/2159-8290.cd-23-0280] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Revised: 04/27/2023] [Accepted: 07/28/2023] [Indexed: 09/16/2023]
Abstract
Small-molecule drugs have enabled the practice of precision oncology for genetically defined patient populations since the first approval of imatinib in 2001. Scientific and technology advances over this 20-year period have driven the evolution of cancer biology, medicinal chemistry, and data science. Collectively, these advances provide tools to more consistently design best-in-class small-molecule drugs against known, previously undruggable, and novel cancer targets. The integration of these tools and their customization in the hands of skilled drug hunters will be necessary to enable the discovery of transformational therapies for patients across a wider spectrum of cancers. SIGNIFICANCE Target-centric small-molecule drug discovery necessitates the consideration of multiple approaches to identify chemical matter that can be optimized into drug candidates. To do this successfully and consistently, drug hunters require a comprehensive toolbox to avoid following the "law of instrument" or Maslow's hammer concept where only one tool is applied regardless of the requirements of the task. Combining our ever-increasing understanding of cancer and cancer targets with the technological advances in drug discovery described below will accelerate the next generation of small-molecule drugs in oncology.
Collapse
Affiliation(s)
| | | | | | | | | | | | - Axel Hoos
- Scorpion Therapeutics, Boston, Massachusetts
| |
Collapse
|
18
|
Croft D, Lodhia P, Lourenco S, MacKay C. Effectively utilizing publicly available databases for cancer target evaluation. NAR Cancer 2023; 5:zcad035. [PMID: 37457379 PMCID: PMC10346432 DOI: 10.1093/narcan/zcad035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Revised: 05/12/2023] [Accepted: 06/20/2023] [Indexed: 07/18/2023] Open
Abstract
The majority of compounds designed against cancer drug targets do not progress to become approved drugs, mainly due to lack of efficacy and/or unmanageable toxicity. Robust target evaluation is therefore required before progressing through the drug discovery process to reduce the high attrition rate. There are a wealth of publicly available databases that can be mined to generate data as part of a target evaluation. It can, however, be challenging to learn what databases are available, how and when they should be used, and to understand the associated limitations. Here, we have compiled and present key, freely accessible and easy-to-use databases that house informative datasets from in vitro, in vivo and clinical studies. We also highlight comprehensive target review databases that aim to bring together information from multiple sources into one-stop portals. In the post-genomics era, a key objective is to exploit the extensive cell, animal and patient characterization datasets in order to deliver precision medicine on a patient-specific basis. Effective utilization of the highlighted databases will go some way towards supporting the cancer research community achieve these aims.
Collapse
Affiliation(s)
- Daniel Croft
- Cancer Research Horizons, The Cancer Research UK Beatson Institute, Glasgow, G61 1BD, UK
| | - Puja Lodhia
- Cancer Research Horizons, The Francis Crick Institute, London, NW1 1AT, UK
| | - Sofia Lourenco
- Cancer Research Horizons, The Francis Crick Institute, London, NW1 1AT, UK
| | - Craig MacKay
- Cancer Research Horizons, The Cancer Research UK Beatson Institute, Glasgow, G61 1BD, UK
| |
Collapse
|
19
|
Aithani L, Alcaide E, Bartunov S, Cooper CDO, Doré AS, Lane TJ, Maclean F, Rucktooa P, Shaw RA, Skerratt SE. Advancing structural biology through breakthroughs in AI. Curr Opin Struct Biol 2023; 80:102601. [PMID: 37182397 DOI: 10.1016/j.sbi.2023.102601] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2022] [Revised: 03/06/2023] [Accepted: 04/03/2023] [Indexed: 05/16/2023]
Abstract
The past century has witnessed an exponential increase in our atomic-level understanding of molecular and cellular mechanisms from a structural perspective, with multiple landmark achievements contributing to the field. This, coupled with recent and continuing breakthroughs in artificial intelligence methods such as AlphaFold2, and enhanced computational power, is enabling our understanding of protein structure and function at unprecedented levels of accuracy and predictivity. Here, we describe some of the major recent advances across these fields, and describe, as these technologies coalesce, the potential to utilise our enhanced knowledge of intricate cellular and molecular systems to discover novel therapeutics to alleviate human suffering.
Collapse
Affiliation(s)
- Laksh Aithani
- CHARM Therapeutics Ltd., The Stanley Building, 7 St. Pancras Square, London, N1C 4AG, UK.
| | - Eric Alcaide
- CHARM Therapeutics Ltd., The Stanley Building, 7 St. Pancras Square, London, N1C 4AG, UK
| | - Sergey Bartunov
- CHARM Therapeutics Ltd., The Stanley Building, 7 St. Pancras Square, London, N1C 4AG, UK
| | - Christopher D O Cooper
- CHARM Therapeutics Ltd., B900, Babraham Research Campus, Babraham, Cambridge, CB22 3AT, UK
| | - Andrew S Doré
- CHARM Therapeutics Ltd., B900, Babraham Research Campus, Babraham, Cambridge, CB22 3AT, UK
| | - Thomas J Lane
- CHARM Therapeutics Ltd., B900, Babraham Research Campus, Babraham, Cambridge, CB22 3AT, UK
| | - Finlay Maclean
- CHARM Therapeutics Ltd., The Stanley Building, 7 St. Pancras Square, London, N1C 4AG, UK
| | - Prakash Rucktooa
- CHARM Therapeutics Ltd., B900, Babraham Research Campus, Babraham, Cambridge, CB22 3AT, UK
| | - Robert A Shaw
- CHARM Therapeutics Ltd., The Stanley Building, 7 St. Pancras Square, London, N1C 4AG, UK
| | - Sarah E Skerratt
- CHARM Therapeutics Ltd., B900, Babraham Research Campus, Babraham, Cambridge, CB22 3AT, UK.
| |
Collapse
|
20
|
Masrati G, Kessel A, Ben-Tal N. Cation/proton antiporters: novel structure-driven pharmaceutical opportunities. Trends Pharmacol Sci 2023; 44:258-262. [PMID: 36934025 DOI: 10.1016/j.tips.2023.02.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2023] [Revised: 02/23/2023] [Accepted: 02/23/2023] [Indexed: 03/18/2023]
Abstract
Cation/proton antiporters (CPAs) regulate cells' salt concentration and pH. Their malfunction is associated with a range of human pathologies, yet only a handful of CPA-targeting therapeutics are presently in clinical development. Here, we discuss how recently published mammalian protein structures and emerging computational technologies may help to bridge this gap.
Collapse
Affiliation(s)
- Gal Masrati
- Tel Aviv University, George S. Wise Faculty of Life Sciences, Department of Biochemistry and Molecular Biology, Tel Aviv, Israel
| | - Amit Kessel
- Tel Aviv University, George S. Wise Faculty of Life Sciences, Department of Biochemistry and Molecular Biology, Tel Aviv, Israel
| | - Nir Ben-Tal
- Tel Aviv University, George S. Wise Faculty of Life Sciences, Department of Biochemistry and Molecular Biology, Tel Aviv, Israel.
| |
Collapse
|
21
|
Díaz-Rovira AM, Martín H, Beuming T, Díaz L, Guallar V, Ray SS. Are Deep Learning Structural Models Sufficiently Accurate for Virtual Screening? Application of Docking Algorithms to AlphaFold2 Predicted Structures. J Chem Inf Model 2023; 63:1668-1674. [PMID: 36892986 DOI: 10.1021/acs.jcim.2c01270] [Citation(s) in RCA: 17] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/10/2023]
Abstract
Machine learning-based protein structure prediction algorithms, such as RosettaFold and AlphaFold2, have greatly impacted the structural biology field, arousing a fair amount of discussion around their potential role in drug discovery. While there are few preliminary studies addressing the usage of these models in virtual screening, none of them focus on the prospect of hit-finding in a real-world virtual screen with a model based on low prior structural information. In order to address this, we have developed an AlphaFold2 version where we exclude all structural templates with more than 30% sequence identity from the model-building process. In a previous study, we used those models in conjunction with state-of-the-art free energy perturbation methods and demonstrated that it is possible to obtain quantitatively accurate results. In this work, we focus on using these structures in rigid receptor-ligand docking studies. Our results indicate that using out-of-the-box Alphafold2 models is not an ideal scenario for virtual screening campaigns; in fact, we strongly recommend to include some post-processing modeling to drive the binding site into a more realistic holo model.
Collapse
Affiliation(s)
- Anna M Díaz-Rovira
- Barcelona Supercomputing Center, Jordi Girona 29, E-08034 Barcelona, Spain
| | | | - Thijs Beuming
- Latham Biopharm Group, 101 Main Street, Suite 1400, Cambridge, Massachusetts 02142, United States
| | - Lucía Díaz
- Nostrum Biodiscovery S.L., E-08029 Barcelona, Spain
| | - Victor Guallar
- Barcelona Supercomputing Center, Jordi Girona 29, E-08034 Barcelona, Spain.,Nostrum Biodiscovery S.L., E-08029 Barcelona, Spain.,ICREA, Passeig Lluís Companys 23, E-08010 Barcelona, Spain
| | - Soumya S Ray
- RA Capital, 200 Berkeley Street, Boston, Massachusetts 02116, United States.,3-Dimensional Consulting, 134 Franklin Avenue, Quincy, Massachusetts 02170, United States
| |
Collapse
|
22
|
Muegge I, Hu Y. Recent Advances in Alchemical Binding Free Energy Calculations for Drug Discovery. ACS Med Chem Lett 2023; 14:244-250. [PMID: 36923913 PMCID: PMC10009785 DOI: 10.1021/acsmedchemlett.2c00541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2023] [Accepted: 02/07/2023] [Indexed: 02/18/2023] Open
Abstract
Rigorous physics-based methods to calculate binding free energies of protein-ligand complexes have become a valued component of structure-based drug design. Relative and absolute binding free energy calculations have been deployed prospectively in support of solving diverse drug discovery challenges. Here we review recent applications of binding free energy calculations to fragment growing and linking, scaffold hopping, binding pose validation, virtual screening, covalent enzyme inhibition, and positional analogue scanning. Furthermore, we discuss the merits of using protein models and highlight recent efforts to replace costly binding free energy calculations with predictions from machine learning models trained on a limited number of free energy perturbation or thermodynamic integration calculations thereby allowing for extended chemical space exploration.
Collapse
Affiliation(s)
- Ingo Muegge
- Alkermes,
Inc, 852 Winter Street, Waltham, Massachusetts 02451-1420, United States
| | - Yuan Hu
- Frontier
Medicines Corp, 451 D
Street, Suite 207, Boston, Massachusetts 02210, United States
| |
Collapse
|
23
|
Kersten C, Clower S, Barthels F. Hic Sunt Dracones: Molecular Docking in Uncharted Territories with Structures from AlphaFold2 and RoseTTAfold. J Chem Inf Model 2023; 63:2218-2225. [PMID: 36884022 DOI: 10.1021/acs.jcim.2c01400] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/09/2023]
Abstract
AlphaFold2 and RoseTTAfold impress with their high accuracy in protein structure prediction. However, for structure-based virtual screenings, not only the overall structure but especially the binding sites need to be accurately predicted. In this work, the docking performance for 66 targets with known ligands but without experimental structures available in the protein data bank was elucidated. The results suggest that using an experimental surrogate-ligand complex is often superior over homology models, and only at low sequence identity to the closest homologue AlphaFold2 structures show an equal performance. The generally high fluctuation of receiver operating characteristic area under the curve values obtained for different homology models suggests that multiple combinations of docking programs and homology models should be tested prior to prospective virtual screenings, and in some cases post-processing of crude models might be necessary.
Collapse
Affiliation(s)
- Christian Kersten
- Institute of Pharmaceutical and Biomedical Sciences, Johannes Gutenberg-University Mainz, Staudingerweg 5, 55128 Mainz, Germany
| | - Steven Clower
- Institute of Pharmaceutical and Biomedical Sciences, Johannes Gutenberg-University Mainz, Staudingerweg 5, 55128 Mainz, Germany
| | - Fabian Barthels
- Institute of Pharmaceutical and Biomedical Sciences, Johannes Gutenberg-University Mainz, Staudingerweg 5, 55128 Mainz, Germany
| |
Collapse
|
24
|
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: 56] [Impact Index Per Article: 56.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
|