1
|
Neha Ashraf S, Henry Blackwell J, Holdgate GA, C C Lucas S, Solovyeva A, Ian Storer R, Whitehurst BC. Hit me with your best shot: Integrated hit discovery for the next generation of drug targets. Drug Discov Today 2024:104143. [PMID: 39173704 DOI: 10.1016/j.drudis.2024.104143] [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: 05/31/2024] [Revised: 08/07/2024] [Accepted: 08/16/2024] [Indexed: 08/24/2024]
Abstract
Identification of high-quality hit chemical matter is of vital importance to the success of drug discovery campaigns. However, this goal is becoming ever harder to achieve as the targets entering the portfolios of pharmaceutical and biotechnology companies are increasingly trending towards novel and traditionally challenging to drug. This demand has fuelled the development and adoption of numerous new screening approaches, whereby the contemporary hit identification toolbox comprises a growing number of orthogonal and complementary technologies including high-throughput screening, fragment-based ligand design, affinity screening (affinity-selection mass spectrometry, differential scanning fluorimetry, DNA-encoded library screening), as well as increasingly sophisticated computational predictive approaches. Herein we describe how an integrated strategy for hit discovery, whereby multiple hit identification techniques are tactically applied, selected in the context of target suitability and resource priority, represents an optimal and often essential approach to maximise the likelihood of identifying quality starting points from which to develop the next generation of medicines.
Collapse
Affiliation(s)
- S Neha Ashraf
- Hit Discovery, Discovery Science, AstraZeneca R&D, Cambridge, CB2 0AA, UK
| | - J Henry Blackwell
- Hit Discovery, Discovery Science, AstraZeneca R&D, Cambridge, CB2 0AA, UK
| | | | - Simon C C Lucas
- Hit Discovery, Discovery Science, AstraZeneca R&D, Cambridge, CB2 0AA, UK
| | - Alisa Solovyeva
- Hit Discovery, Discovery Science, AstraZeneca R&D, Gothenburg, SE-431 83, Sweden
| | - R Ian Storer
- Hit Discovery, Discovery Science, AstraZeneca R&D, Cambridge, CB2 0AA, UK.
| | | |
Collapse
|
2
|
Crivelli-Decker JE, Beckwith Z, Tom G, Le L, Khuttan S, Salomon-Ferrer R, Beall J, Gómez-Bombarelli R, Bortolato A. Machine Learning Guided AQFEP: A Fast and Efficient Absolute Free Energy Perturbation Solution for Virtual Screening. J Chem Theory Comput 2024. [PMID: 39146234 DOI: 10.1021/acs.jctc.4c00399] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/17/2024]
Abstract
Structure-based methods in drug discovery have become an integral part of the modern drug discovery process. The power of virtual screening lies in its ability to rapidly and cost-effectively explore enormous chemical spaces to select promising ligands for further experimental investigation. Relative free energy perturbation (RFEP) and similar methods are the gold standard for binding affinity prediction in drug discovery hit-to-lead and lead optimization phases, but have high computational cost and the requirement of a structural analog with a known activity. Without a reference molecule requirement, absolute FEP (AFEP) has, in theory, better accuracy for hit ID, but in practice, the slow throughput is not compatible with VS, where fast docking and unreliable scoring functions are still the standard. Here, we present an integrated workflow to virtually screen large and diverse chemical libraries efficiently, combining active learning with a physics-based scoring function based on a fast absolute free energy perturbation method. We validated the performance of the approach in the ranking of structurally related ligands, virtual screening hit rate enrichment, and active learning chemical space exploration; disclosing the largest reported collection of free energy simulations to date.
Collapse
Affiliation(s)
| | - Zane Beckwith
- SandboxAQ, Palo Alto, California 94301, United States
| | - Gary Tom
- SandboxAQ, Palo Alto, California 94301, United States
- Department of Chemistry and Department of Computer Science, University of Toronto, Toronto, ON M5S 3H6, Canada
- Vector Institute for Artificial Intelligence, Toronto, ON M5S 3H6, Canada
| | - Ly Le
- SandboxAQ, Palo Alto, California 94301, United States
| | - Sheenam Khuttan
- SandboxAQ, Palo Alto, California 94301, United States
- Department of Chemistry, Brooklyn College of the City University of New York, Brooklyn, New York 11367, United States
| | | | - Jackson Beall
- SandboxAQ, Palo Alto, California 94301, United States
| | - Rafael Gómez-Bombarelli
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | | |
Collapse
|
3
|
Sampson JM, Cannon DA, Duan J, Epstein JCK, Sergeeva AP, Katsamba PS, Mannepalli SM, Bahna FA, Adihou H, Guéret SM, Gopalakrishnan R, Geschwindner S, Rees DG, Sigurdardottir A, Wilkinson T, Dodd RB, De Maria L, Mobarec JC, Shapiro L, Honig B, Buchanan A, Friesner RA, Wang L. Robust Prediction of Relative Binding Energies for Protein-Protein Complex Mutations Using Free Energy Perturbation Calculations. J Mol Biol 2024; 436:168640. [PMID: 38844044 PMCID: PMC11339910 DOI: 10.1016/j.jmb.2024.168640] [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] [Received: 04/25/2024] [Accepted: 05/31/2024] [Indexed: 06/18/2024]
Abstract
Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how Protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.
Collapse
Affiliation(s)
- Jared M Sampson
- Schrödinger, Inc., Life Sciences Software, New York, NY, USA
| | - Daniel A Cannon
- Schrödinger, GmbH, Life Sciences Software, Mannheim, Germany
| | - Jianxin Duan
- Schrödinger, GmbH, Life Sciences Software, Mannheim, Germany
| | | | - Alina P Sergeeva
- Columbia University, Department of Systems Biology, New York, NY, USA
| | | | - Seetha M Mannepalli
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA
| | - Fabiana A Bahna
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA
| | - Hélène Adihou
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, Gothenburg, Sweden; Max Planck Institute of Molecular Physiology, AstraZeneca-MPI Satellite Unit, Dortmund, Germany
| | - Stéphanie M Guéret
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, Gothenburg, Sweden; Max Planck Institute of Molecular Physiology, AstraZeneca-MPI Satellite Unit, Dortmund, Germany
| | - Ranganath Gopalakrishnan
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, Gothenburg, Sweden; Max Planck Institute of Molecular Physiology, AstraZeneca-MPI Satellite Unit, Dortmund, Germany
| | - Stefan Geschwindner
- AstraZeneca, Mechanistic and Structural Biology, Discovery Sciences, R&D, Gothenburg, Sweden
| | - D Gareth Rees
- AstraZeneca, Biologics Engineering, R&D, Cambridge, UK
| | | | | | - Roger B Dodd
- AstraZeneca, Biologics Engineering, R&D, Cambridge, UK
| | - Leonardo De Maria
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Respiratory and Immunology, BioPharmaceuticals R&D, Gothenburg, Sweden
| | - Juan Carlos Mobarec
- AstraZeneca, Mechanistic and Structural Biology, Discovery Sciences, R&D, Cambridge, UK
| | - Lawrence Shapiro
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA; Columbia University, Department of Biochemistry and Molecular Biophysics, New York, NY, USA
| | - Barry Honig
- Columbia University, Department of Systems Biology, New York, NY, USA; Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA; Columbia University, Department of Biochemistry and Molecular Biophysics, New York, NY, USA; Columbia University, Department of Medicine, New York, NY, USA
| | | | | | - Lingle Wang
- Schrödinger, Inc., Life Sciences Software, New York, NY, USA.
| |
Collapse
|
4
|
Heinzelmann G, Huggins DJ, Gilson MK. BAT2: an Open-Source Tool for Flexible, Automated, and Low Cost Absolute Binding Free Energy Calculations. J Chem Theory Comput 2024; 20:6518-6530. [PMID: 39088306 PMCID: PMC11325538 DOI: 10.1021/acs.jctc.4c00205] [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: 08/03/2024]
Abstract
Absolute binding free energy (ABFE) calculations with all-atom molecular dynamics (MD) have the potential to greatly reduce costs in the first stages of drug discovery. Here, we introduce BAT2, the new version of the Binding Affinity Tool (BAT.py), designed to combine full automation of ABFE calculations with high-performance MD simulations, making it a potential tool for virtual screening. We describe and test several changes and new features that were incorporated into the code, such as relative restraints between the protein and the ligand instead of using fixed dummy atoms, support for the OpenMM simulation engine, a merged approach to the application/release of restraints, support for cobinders and proteins with multiple chains, and many others. We also reduced the simulation times for each ABFE calculation, assessing the effect on the expected robustness and accuracy of the calculations.
Collapse
Affiliation(s)
- Germano Heinzelmann
- Departamento de Fisica, Universidade Federal de Santa Catarina, Florianopolis 88040-970, Brasil
| | - David J Huggins
- Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York 10065, United States
- Sanders Tri-Institutional Therapeutics Discovery Institute, 1230 York Avenue, Box 122, New York, New York 10065, United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego 92093, United States
| |
Collapse
|
5
|
Wu Y, Liu F, Glenn I, Fonseca-Valencia K, Paris L, Xiong Y, Jerome SV, Brooks CL, Shoichet BK. Identifying Artifacts from Large Library Docking. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.07.17.603966. [PMID: 39071262 PMCID: PMC11275789 DOI: 10.1101/2024.07.17.603966] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
While large library docking has discovered potent ligands for multiple targets, as the libraries have grown, the very top of the hit-lists can become populated with artifacts that cheat our scoring functions. Though these cheating molecules are rare, they become ever-more dominant with library growth. Here, we investigate rescoring top-ranked molecules from docking screens with orthogonal methods to identify these artifacts, exploring implicit solvent models and absolute binding free energy perturbation (AB-FEP) as cross-filters. In retrospective studies, this approach deprioritized high-ranking non-binders for nine targets while leaving true ligands relatively unaffected. We tested the method prospectively against results from large library docking AmpC β-lactamase. From the very top of the docking hit lists, we prioritized 128 molecules for synthesis and experimental testing, a mixture of 39 molecules that rescoring flagged as likely cheaters and another 89 that were plausible true actives. None of the 39 predicted cheating compounds inhibited AmpC up to 200 μ M in enzyme assays, while 57% of the 89 plausible true actives did do so, with 19 of them inhibiting the enzyme with apparentK i values better than 50 μ M . As our libraries continue to grow, a strategy of catching docking artifacts by rescoring with orthogonal methods may find wide use in the field.
Collapse
Affiliation(s)
- Yujin Wu
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 94158, United States
| | - Fangyu Liu
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 94158, United States
| | - Isabella Glenn
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 94158, United States
| | - Karla Fonseca-Valencia
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 94158, United States
| | - Lu Paris
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 94158, United States
| | - Yuyue Xiong
- Schrödinger, Inc., 9868 Scranton Road, San Diego, California 92121, United States
| | - Steven V. Jerome
- Schrödinger, Inc., 1540 Broadway, New York, New York, 10036, United States
| | - Charles L. Brooks
- Biophysics Program, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Brian K. Shoichet
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 94158, United States
| |
Collapse
|
6
|
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] [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
|
7
|
Chen S, Xie J, Ye R, Xu DD, Yang Y. Structure-aware dual-target drug design through collaborative learning of pharmacophore combination and molecular simulation. Chem Sci 2024; 15:10366-10380. [PMID: 38994407 PMCID: PMC11234869 DOI: 10.1039/d4sc00094c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Accepted: 06/09/2024] [Indexed: 07/13/2024] Open
Abstract
Dual-target drug design has gained significant attention in the treatment of complex diseases, such as cancers and autoimmune disorders. A widely employed design strategy is combining pharmacophores to leverage the knowledge of structure-activity relationships of both targets. Unfortunately, pharmacophore combination often struggles with long and expensive trial and error, because the protein pockets of the two targets impose complex structural constraints. In this study, we propose AIxFuse, a structure-aware dual-target drug design method that learns pharmacophore fusion patterns to satisfy the dual-target structural constraints simulated by molecular docking. AIxFuse employs two self-play reinforcement learning (RL) agents to learn pharmacophore selection and fusion by comprehensive feedback including dual-target molecular docking scores. Collaboratively, the molecular docking scores are learned by active learning (AL). Through collaborative RL and AL, AIxFuse learns to generate molecules with multiple desired properties. AIxFuse is shown to outperform state-of-the-art methods in generating dual-target drugs against glycogen synthase kinase-3 beta (GSK3β) and c-Jun N-terminal kinase 3 (JNK3). When applied to another task against retinoic acid receptor-related orphan receptor γ-t (RORγt) and dihydroorotate dehydrogenase (DHODH), AIxFuse exhibits consistent performance while compared methods suffer from performance drops, leading to a 5 times higher performance in success rate. Docking studies demonstrate that AIxFuse can generate molecules concurrently satisfying the binding mode required by both targets. Further free energy perturbation calculation indicates that the generated candidates have promising binding free energies against both targets.
Collapse
Affiliation(s)
- Sheng Chen
- School of Computer Science and Engineering, Sun Yat-sen University Guangzhou 510006 China
- AixplorerBio Inc. Jiaxing 314031 China
| | - Junjie Xie
- School of Computer Science and Engineering, Sun Yat-sen University Guangzhou 510006 China
- AixplorerBio Inc. Jiaxing 314031 China
| | | | | | - Yuedong Yang
- School of Computer Science and Engineering, Sun Yat-sen University Guangzhou 510006 China
| |
Collapse
|
8
|
Mittal A, Martin MF, Levin EJ, Adams C, Yang M, Provins L, Hall A, Procter M, Ledecq M, Hillisch A, Wolff C, Gillard M, Horanyi PS, Coleman JA. Structures of synaptic vesicle protein 2A and 2B bound to anticonvulsants. Nat Struct Mol Biol 2024:10.1038/s41594-024-01335-1. [PMID: 38898101 DOI: 10.1038/s41594-024-01335-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Accepted: 05/14/2024] [Indexed: 06/21/2024]
Abstract
Epilepsy is a common neurological disorder characterized by abnormal activity of neuronal networks, leading to seizures. The racetam class of anti-seizure medications bind specifically to a membrane protein found in the synaptic vesicles of neurons called synaptic vesicle protein 2 (SV2) A (SV2A). SV2A belongs to an orphan subfamily of the solute carrier 22 organic ion transporter family that also includes SV2B and SV2C. The molecular basis for how anti-seizure medications act on SV2s remains unknown. Here we report cryo-electron microscopy structures of SV2A and SV2B captured in a luminal-occluded conformation complexed with anticonvulsant ligands. The conformation bound by anticonvulsants resembles an inhibited transporter with closed luminal and intracellular gates. Anticonvulsants bind to a highly conserved central site in SV2s. These structures provide blueprints for future drug design and will facilitate future investigations into the biological function of SV2s.
Collapse
Affiliation(s)
- Anshumali Mittal
- Department of Structural Biology, University of Pittsburgh, Pittsburgh, PA, USA
| | - Matthew F Martin
- Department of Structural Biology, University of Pittsburgh, Pittsburgh, PA, USA
| | | | | | | | | | | | | | | | | | | | | | | | - Jonathan A Coleman
- Department of Structural Biology, University of Pittsburgh, Pittsburgh, PA, USA.
| |
Collapse
|
9
|
Gutkin E, Gusev F, Gentile F, Ban F, Koby SB, Narangoda C, Isayev O, Cherkasov A, Kurnikova MG. In silico screening of LRRK2 WDR domain inhibitors using deep docking and free energy simulations. Chem Sci 2024; 15:8800-8812. [PMID: 38873063 PMCID: PMC11168082 DOI: 10.1039/d3sc06880c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Accepted: 04/10/2024] [Indexed: 06/15/2024] Open
Abstract
The Critical Assessment of Computational Hit-Finding Experiments (CACHE) Challenge series is focused on identifying small molecule inhibitors of protein targets using computational methods. Each challenge contains two phases, hit-finding and follow-up optimization, each of which is followed by experimental validation of the computational predictions. For the CACHE Challenge #1, the Leucine-Rich Repeat Kinase 2 (LRRK2) WD40 Repeat (WDR) domain was selected as the target for in silico hit-finding and optimization. Mutations in LRRK2 are the most common genetic cause of the familial form of Parkinson's disease. The LRRK2 WDR domain is an understudied drug target with no known molecular inhibitors. Herein we detail the first phase of our winning submission to the CACHE Challenge #1. We developed a framework for the high-throughput structure-based virtual screening of a chemically diverse small molecule space. Hit identification was performed using the large-scale Deep Docking (DD) protocol followed by absolute binding free energy (ABFE) simulations. ABFEs were computed using an automated molecular dynamics (MD)-based thermodynamic integration (TI) approach. 4.1 billion ligands from Enamine REAL were screened with DD followed by ABFEs computed by MD TI for 793 ligands. 76 ligands were prioritized for experimental validation, with 59 compounds successfully synthesized and 5 compounds identified as hits, yielding a 8.5% hit rate. Our results demonstrate the efficacy of the combined DD and ABFE approaches for hit identification for a target with no previously known hits. This approach is widely applicable for the efficient screening of ultra-large chemical libraries as well as rigorous protein-ligand binding affinity estimation leveraging modern computational resources.
Collapse
Affiliation(s)
- Evgeny Gutkin
- Department of Chemistry, Mellon College of Science, Carnegie Mellon University Pittsburgh PA 15213 USA
| | - Filipp Gusev
- Department of Chemistry, Mellon College of Science, Carnegie Mellon University Pittsburgh PA 15213 USA
- Computational Biology Department, School of Computer Science, Carnegie Mellon University Pittsburgh PA 15213 USA
| | - Francesco Gentile
- Department of Chemistry and Biomolecular Sciences, University of Ottawa Ottawa ON Canada
- Ottawa Institute of Systems Biology Ottawa ON Canada
| | - Fuqiang Ban
- Vancouver Prostate Centre, The University of British Columbia Vancouver BC Canada
| | - S Benjamin Koby
- Department of Chemistry, Mellon College of Science, Carnegie Mellon University Pittsburgh PA 15213 USA
| | - Chamali Narangoda
- Department of Chemistry, Mellon College of Science, Carnegie Mellon University Pittsburgh PA 15213 USA
| | - Olexandr Isayev
- Department of Chemistry, Mellon College of Science, Carnegie Mellon University Pittsburgh PA 15213 USA
- Computational Biology Department, School of Computer Science, Carnegie Mellon University Pittsburgh PA 15213 USA
| | - Artem Cherkasov
- Vancouver Prostate Centre, The University of British Columbia Vancouver BC Canada
| | - Maria G Kurnikova
- Department of Chemistry, Mellon College of Science, Carnegie Mellon University Pittsburgh PA 15213 USA
| |
Collapse
|
10
|
Adediwura VA, Koirala K, Do HN, Wang J, Miao Y. Understanding the impact of binding free energy and kinetics calculations in modern drug discovery. Expert Opin Drug Discov 2024; 19:671-682. [PMID: 38722032 PMCID: PMC11108734 DOI: 10.1080/17460441.2024.2349149] [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] [Received: 07/27/2023] [Accepted: 04/25/2024] [Indexed: 05/22/2024]
Abstract
INTRODUCTION For rational drug design, it is crucial to understand the receptor-drug binding processes and mechanisms. A new era for the use of computer simulations in predicting drug-receptor interactions at an atomic level has begun with remarkable advances in supercomputing and methodological breakthroughs. AREAS COVERED End-point free energy calculation methods such as Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) or Molecular-Mechanics/Generalized Born Surface Area (MM/GBSA), free energy perturbation (FEP), and thermodynamic integration (TI) are commonly used for binding free energy calculations in drug discovery. In addition, kinetic dissociation and association rate constants (k off and k on ) play critical roles in the function of drugs. Nowadays, Molecular Dynamics (MD) and enhanced sampling simulations are increasingly being used in drug discovery. Here, the authors provide a review of the computational techniques used in drug binding free energy and kinetics calculations. EXPERT OPINION The applications of computational methods in drug discovery and design are expanding, thanks to improved predictions of the binding free energy and kinetic rates of drug molecules. Recent microsecond-timescale enhanced sampling simulations have made it possible to accurately capture repetitive ligand binding and dissociation, facilitating more efficient and accurate calculations of ligand binding free energy and kinetics.
Collapse
Affiliation(s)
- Victor A. Adediwura
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Kushal Koirala
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Hung N. Do
- Center for Computational Biology, University of Kansas, Lawrence, KS, USA
- Present address: Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
| | - Jinan Wang
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Yinglong Miao
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| |
Collapse
|
11
|
Boresch S. On Analytical Corrections for Restraints in Absolute Binding Free Energy Calculations. J Chem Inf Model 2024; 64:3605-3609. [PMID: 38640478 PMCID: PMC11094717 DOI: 10.1021/acs.jcim.4c00442] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2024] [Revised: 04/16/2024] [Accepted: 04/17/2024] [Indexed: 04/21/2024]
Abstract
Double decoupling absolute binding free energy simulations require an intermediate state at which the ligand is held solely by restraints in a position and orientation resembling the bound state. One possible choice consists of one distance, two angle, and three dihedral angle restraints. Here, I demonstrate that in practically all cases the analytical correction derived under the rigid rotator harmonic oscillator approximation is sufficient to account for the free energy of the restraints.
Collapse
Affiliation(s)
- Stefan Boresch
- Department of Chemistry, University
of Vienna, Währinger
Straße 17, A-1090 Vienna, Austria
| |
Collapse
|
12
|
Sampson JM, Cannon DA, Duan J, Epstein JCK, Sergeeva AP, Katsamba PS, Mannepalli SM, Bahna FA, Adihou H, Guéret SM, Gopalakrishnan R, Geschwindner S, Rees DG, Sigurdardottir A, Wilkinson T, Dodd RB, De Maria L, Mobarec JC, Shapiro L, Honig B, Buchanan A, Friesner RA, Wang L. Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.04.22.590325. [PMID: 38712280 PMCID: PMC11071377 DOI: 10.1101/2024.04.22.590325] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2024]
Abstract
Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.
Collapse
Affiliation(s)
| | | | - Jianxin Duan
- Schrödinger, GmbH, Life Sciences Software, Mannheim, Germany
| | | | - Alina P. Sergeeva
- Columbia University, Department of Systems Biology, New York, NY, USA
| | | | - Seetha M. Mannepalli
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA, 10027
| | - Fabiana A. Bahna
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA, 10027
| | - Hélène Adihou
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, Gothenburg, Sweden
- Max Planck Institute of Molecular Physiology, AstraZeneca-MPI Satellite Unit, Dortmund, Germany
| | - Stéphanie M. Guéret
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, Gothenburg, Sweden
- Max Planck Institute of Molecular Physiology, AstraZeneca-MPI Satellite Unit, Dortmund, Germany
| | - Ranganath Gopalakrishnan
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, Gothenburg, Sweden
- Max Planck Institute of Molecular Physiology, AstraZeneca-MPI Satellite Unit, Dortmund, Germany
| | - Stefan Geschwindner
- AstraZeneca, Mechanistic and Structural Biology, Discovery Sciences, R&D, Cambridge, UK
| | | | | | | | - Roger B. Dodd
- AstraZeneca, Biologics Engineering, R&D, Cambridge, UK
| | - Leonardo De Maria
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Respiratory and Immunology, BioPharmaceuticals R&D, Gothenburg, Sweden
| | - Juan Carlos Mobarec
- AstraZeneca, Mechanistic and Structural Biology, Discovery Sciences, R&D, Cambridge, UK
| | - Lawrence Shapiro
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA, 10027
- Columbia University, Department of Biochemistry and Molecular Biophysics, New York, NY, USA
| | - Barry Honig
- Columbia University, Department of Systems Biology, New York, NY, USA
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA, 10027
- Columbia University, Department of Biochemistry and Molecular Biophysics, New York, NY, USA
- Columbia University, Department of Medicine, New York, NY, USA
| | | | | | - Lingle Wang
- Schrödinger, Inc., Life Sciences Software, New York, NY, USA
| |
Collapse
|
13
|
Macaya L, González D, Vöhringer-Martinez E. Nonbonded Force Field Parameters from MBIS Partitioning of the Molecular Electron Density Improve Binding Affinity Predictions of the T4-Lysozyme Double Mutant. J Chem Inf Model 2024; 64:3269-3277. [PMID: 38546407 DOI: 10.1021/acs.jcim.3c01912] [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/23/2024]
Abstract
The use of computer simulation for binding affinity prediction is growing in drug discovery. However, its wider use is constrained by the accuracy of the free energy calculations. The key sources of error are the force fields used to depict molecular interactions and insufficient sampling of the configurational space. To improve the quality of the force field, we developed a Python-based computational workflow. The workflow described here uses the minimal basis iterative stockholder (MBIS) method to determine atomic charges and Lennard-Jones parameters from the polarized molecular density. This is done by performing electronic structure calculations on various configurations of the ligand when it is both bound and unbound. In addition, we validated a simulation procedure that accounts for the protein and ligand degrees of freedom to precisely calculate binding free energies. This was achieved by comparing the self-adjusted mixture sampling and nonequilibrium thermodynamic integration methods using various protein and ligand conformations. The accuracy of predicting binding affinity is improved by using MBIS-derived force field parameters and a validated simulation procedure. This improvement surpasses the chemical precision for the eight aromatic ligands, reaching a root-mean-square error of 0.7 kcal/mol.
Collapse
Affiliation(s)
- Luis Macaya
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, 4070386 Concepción, Chile
| | - Duván González
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, 4070386 Concepción, Chile
| | - Esteban Vöhringer-Martinez
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, 4070386 Concepción, Chile
| |
Collapse
|
14
|
Atz K, Cotos L, Isert C, Håkansson M, Focht D, Hilleke M, Nippa DF, Iff M, Ledergerber J, Schiebroek CCG, Romeo V, Hiss JA, Merk D, Schneider P, Kuhn B, Grether U, Schneider G. Prospective de novo drug design with deep interactome learning. Nat Commun 2024; 15:3408. [PMID: 38649351 PMCID: PMC11035696 DOI: 10.1038/s41467-024-47613-w] [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] [Received: 09/13/2023] [Accepted: 04/02/2024] [Indexed: 04/25/2024] Open
Abstract
De novo drug design aims to generate molecules from scratch that possess specific chemical and pharmacological properties. We present a computational approach utilizing interactome-based deep learning for ligand- and structure-based generation of drug-like molecules. This method capitalizes on the unique strengths of both graph neural networks and chemical language models, offering an alternative to the need for application-specific reinforcement, transfer, or few-shot learning. It enables the "zero-shot" construction of compound libraries tailored to possess specific bioactivity, synthesizability, and structural novelty. In order to proactively evaluate the deep interactome learning framework for protein structure-based drug design, potential new ligands targeting the binding site of the human peroxisome proliferator-activated receptor (PPAR) subtype gamma are generated. The top-ranking designs are chemically synthesized and computationally, biophysically, and biochemically characterized. Potent PPAR partial agonists are identified, demonstrating favorable activity and the desired selectivity profiles for both nuclear receptors and off-target interactions. Crystal structure determination of the ligand-receptor complex confirms the anticipated binding mode. This successful outcome positively advocates interactome-based de novo design for application in bioorganic and medicinal chemistry, enabling the creation of innovative bioactive molecules.
Collapse
Affiliation(s)
- Kenneth Atz
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 4, 8093, Zurich, Switzerland
| | - Leandro Cotos
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 4, 8093, Zurich, Switzerland
| | - Clemens Isert
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 4, 8093, Zurich, Switzerland
| | - Maria Håkansson
- SARomics Biostructures AB, Medicon Village, SE-223 81, Lund, Sweden
| | - Dorota Focht
- SARomics Biostructures AB, Medicon Village, SE-223 81, Lund, Sweden
| | - Mattis Hilleke
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 4, 8093, Zurich, Switzerland
| | - David F Nippa
- Roche Pharma Research and Early Development (pRED), Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, CH-4070, Basel, Switzerland
- Department of Pharmacy, Ludwig-Maximilians-Universität München, Butenandtstrasse 5, 81377, Munich, Germany
| | - Michael Iff
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 4, 8093, Zurich, Switzerland
| | - Jann Ledergerber
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 4, 8093, Zurich, Switzerland
| | - Carl C G Schiebroek
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 4, 8093, Zurich, Switzerland
| | - Valentina Romeo
- Roche Pharma Research and Early Development (pRED), Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, CH-4070, Basel, Switzerland
| | - Jan A Hiss
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 4, 8093, Zurich, Switzerland
| | - Daniel Merk
- Department of Pharmacy, Ludwig-Maximilians-Universität München, Butenandtstrasse 5, 81377, Munich, Germany
| | - Petra Schneider
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 4, 8093, Zurich, Switzerland
| | - Bernd Kuhn
- Roche Pharma Research and Early Development (pRED), Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, CH-4070, Basel, Switzerland
| | - Uwe Grether
- Roche Pharma Research and Early Development (pRED), Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, CH-4070, Basel, Switzerland
| | - Gisbert Schneider
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 4, 8093, Zurich, Switzerland.
| |
Collapse
|
15
|
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
|
16
|
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
|
17
|
Awoonor-Williams E, Abu-Saleh AAAA. Molecular Insights into the Impact of Mutations on the Binding Affinity of Targeted Covalent Inhibitors of BTK. J Phys Chem B 2024; 128:2874-2884. [PMID: 38502552 DOI: 10.1021/acs.jpcb.4c00310] [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: 03/21/2024]
Abstract
Targeted covalent inhibitors (TCIs) have witnessed a significant resurgence in recent years, particularly in the kinase drug discovery field for treating diverse clinical indications. The inhibition of Bruton's tyrosine kinase (BTK) for treating B-cell cancers is a classic example where TCIs such as ibrutinib have had breakthroughs in targeted therapy. However, selectivity remains challenging, and the emergence of resistance mutations is a critical concern for clinical efficacy. Computational methods that can accurately predict the impact of mutations on inhibitor binding affinity could prove helpful in informing targeted approaches─providing insights into drug resistance mechanisms. In addition, such systems could help guide the systematic evaluation and impact of mutations in disease models for optimal experimental design. Here, we have employed in silico physics-based methods to understand the effects of mutations on the binding affinity and conformational dynamics of select TCIs of BTK. The TCIs studied include ibrutinib, acalabrutinib, and zanubrutinib─all of which are FDA-approved drugs for treating multiple forms of leukemia and lymphoma. Our results offer useful molecular insights into the structural determinants, thermodynamics, and conformational energies that impact ligand binding for this biological target of clinical relevance.
Collapse
Affiliation(s)
- Ernest Awoonor-Williams
- Department of Chemistry, Memorial University of Newfoundland, St. John's, NL A1B 3X7, Canada
| | - Abd Al-Aziz A Abu-Saleh
- Department of Chemistry and Biochemistry, University of Windsor, Windsor, ON N9B 3P4, Canada
| |
Collapse
|
18
|
Tu G, Fu T, Zheng G, Xu B, Gou R, Luo D, Wang P, Xue W. Computational Chemistry in Structure-Based Solute Carrier Transporter Drug Design: Recent Advances and Future Perspectives. J Chem Inf Model 2024; 64:1433-1455. [PMID: 38294194 DOI: 10.1021/acs.jcim.3c01736] [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: 02/01/2024]
Abstract
Solute carrier transporters (SLCs) are a class of important transmembrane proteins that are involved in the transportation of diverse solute ions and small molecules into cells. There are approximately 450 SLCs within the human body, and more than a quarter of them are emerging as attractive therapeutic targets for multiple complex diseases, e.g., depression, cancer, and diabetes. However, only 44 unique transporters (∼9.8% of the SLC superfamily) with 3D structures and specific binding sites have been reported. To design innovative and effective drugs targeting diverse SLCs, there are a number of obstacles that need to be overcome. However, computational chemistry, including physics-based molecular modeling and machine learning- and deep learning-based artificial intelligence (AI), provides an alternative and complementary way to the classical drug discovery approach. Here, we present a comprehensive overview on recent advances and existing challenges of the computational techniques in structure-based drug design of SLCs from three main aspects: (i) characterizing multiple conformations of the proteins during the functional process of transportation, (ii) identifying druggability sites especially the cryptic allosteric ones on the transporters for substrates and drugs binding, and (iii) discovering diverse small molecules or synthetic protein binders targeting the binding sites. This work is expected to provide guidelines for a deep understanding of the structure and function of the SLC superfamily to facilitate rational design of novel modulators of the transporters with the aid of state-of-the-art computational chemistry technologies including artificial intelligence.
Collapse
Affiliation(s)
- Gao Tu
- Chongqing Key Laboratory of Natural Product Synthesis and Drug Research, School of Pharmaceutical Sciences, Chongqing University, Chongqing 401331, China
| | - Tingting Fu
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
| | | | - Binbin Xu
- Chengdu Sintanovo Biotechnology Co., Ltd., Chengdu 610200, China
| | - Rongpei Gou
- Chongqing Key Laboratory of Natural Product Synthesis and Drug Research, School of Pharmaceutical Sciences, Chongqing University, Chongqing 401331, China
| | - Ding Luo
- Chongqing Key Laboratory of Natural Product Synthesis and Drug Research, School of Pharmaceutical Sciences, Chongqing University, Chongqing 401331, China
| | - Panpan Wang
- College of Chemistry and Pharmaceutical Engineering, Huanghuai University, Zhumadian 463000, China
| | - Weiwei Xue
- Chongqing Key Laboratory of Natural Product Synthesis and Drug Research, School of Pharmaceutical Sciences, Chongqing University, Chongqing 401331, China
| |
Collapse
|
19
|
Karrenbrock M, Rizzi V, Procacci P, Gervasio FL. Addressing Suboptimal Poses in Nonequilibrium Alchemical Calculations. J Phys Chem B 2024; 128:1595-1605. [PMID: 38323915 DOI: 10.1021/acs.jpcb.3c06516] [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: 02/08/2024]
Abstract
Alchemical transformations can be used to quantitatively estimate absolute binding free energies at a reasonable computational cost. However, most of the approaches currently in use require knowledge of the correct (crystallographic) pose. In this paper, we present a combined Hamiltonian replica exchange nonequilibrium alchemical method that allows us to reliably calculate absolute binding free energies, even when starting from suboptimal initial binding poses. Performing a preliminary Hamiltonian replica exchange enhances the sampling of slow degrees of freedom of the ligand and the target, allowing the system to populate the correct binding pose when starting from an approximate docking pose. We apply the method on 6 ligands of the first bromodomain of the BRD4 bromodomain-containing protein. For each ligand, we start nonequilibrium alchemical transformations from both the crystallographic pose and the top-scoring docked pose that are often significantly different. We show that the method produces statistically equivalent binding free energies, making it a useful tool for computational drug discovery pipelines.
Collapse
Affiliation(s)
- Maurice Karrenbrock
- School of Pharmaceutical Sciences, University of Geneva, Rue Michel-Servet 1, CH-1206 Geneva, Switzerland
| | - Valerio Rizzi
- School of Pharmaceutical Sciences, University of Geneva, Rue Michel-Servet 1, CH-1206 Geneva, Switzerland
| | - Piero Procacci
- Chemistry Department, University of Florence, Via della Lastruccia 3-13, 50019 Sesto Fiorentino, Italy
| | - Francesco Luigi Gervasio
- School of Pharmaceutical Sciences, University of Geneva, Rue Michel-Servet 1, CH-1206 Geneva, Switzerland
- Institute of Pharmaceutical Sciences of Western Switzerland, University of Geneva, CH-1206 Geneva, Switzerland
- Chemistry Department, University College London (UCL), WC1E 6BT London, U.K
- Swiss Bioinformatics Institute, University of Geneva, CH-1206 Geneva, Switzerland
| |
Collapse
|
20
|
Bansal N, Wang Y, Sciabola S. Machine Learning Methods as a Cost-Effective Alternative to Physics-Based Binding Free Energy Calculations. Molecules 2024; 29:830. [PMID: 38398581 PMCID: PMC10893267 DOI: 10.3390/molecules29040830] [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/20/2023] [Revised: 01/24/2024] [Accepted: 02/09/2024] [Indexed: 02/25/2024] Open
Abstract
The rank ordering of ligands remains one of the most attractive challenges in drug discovery. While physics-based in silico binding affinity methods dominate the field, they still have problems, which largely revolve around forcefield accuracy and sampling. Recent advances in machine learning have gained traction for protein-ligand binding affinity predictions in early drug discovery programs. In this article, we perform retrospective binding free energy evaluations for 172 compounds from our internal collection spread over four different protein targets and five congeneric ligand series. We compared multiple state-of-the-art free energy methods ranging from physics-based methods with different levels of complexity and conformational sampling to state-of-the-art machine-learning-based methods that were available to us. Overall, we found that physics-based methods behaved particularly well when the ligand perturbations were made in the solvation region, and they did not perform as well when accounting for large conformational changes in protein active sites. On the other end, machine-learning-based methods offer a good cost-effective alternative for binding free energy calculations, but the accuracy of their predictions is highly dependent on the experimental data available for training the model.
Collapse
Affiliation(s)
- Nupur Bansal
- Biotherapeutic and Medicinal Sciences, Biogen, 225 Binney Street, Cambridge, MA 02142, USA; (Y.W.); (S.S.)
| | | | | |
Collapse
|
21
|
Chaturvedi S, Pandya N, Sadhukhan S, Sonawane A. Identification of selective plant-derived natural carotenoid and flavonoids as the potential inhibitors of DHHC-mediated protein S-palmitoylation: an in silico study. J Biomol Struct Dyn 2024:1-14. [PMID: 38319030 DOI: 10.1080/07391102.2024.2306502] [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: 08/20/2023] [Accepted: 01/10/2024] [Indexed: 02/07/2024]
Abstract
Protein S-palmitoylation mediated by DHHCs is recognized as a distinct and reversible form of lipid modification connected with several health perturbations, including neurodegenerative disorders, cancer, and autoimmune conditions. However, the pharmacological characteristics of current pan-DHHC inhibitors, particularly their toxicity and off-target effects, have hindered their in-depth cellular investigations. The therapeutic properties of the natural compounds, with minimal side effects, allowed us to evaluate them as DHHC-targeting inhibitors. Here, we performed an insilico screening of 115 phytochemicals to assess their interactions with the DHHC20 binding site. Among these compounds, lutein, 5-hydroxyflavone, and 6-hydroxyflavone exhibited higher binding energy (-9.2, -8.5, and -8.5 kcal/mol) in the DHHC20 groove compared to pan-DHHC inhibitor 2-BP (-7.0 kcal/mol). Furthermore, we conducted a 100 ns MD simulation to evaluate the stability of these complexes under physiological conditions. The MDsimulation results indicated that DHHC20 formed a more stable conformation with lutein compared to 5-hydroxyflavone and 6-hyroxyflavone via hydrophobic and H-bond interactions. Conclusively, these results could serve as a promising starting point for exploring the use of these natural molecules as DHHC20 inhibitors.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Suchi Chaturvedi
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Simrol, Madhya Pradesh, India
| | - Nirali Pandya
- Department of Chemistry, National University of Singapore, Singapore, Singapore
- Department of Pharmacology and Regenerative Medicine, University of Illinois Chicago, Chicago, IL, USA
| | - Sushabhan Sadhukhan
- Department of Chemistry, Indian Institute of Technology Palakkad, Palakkad, Kerala, India
- Physical & Chemical Biology Laboratory and Department of Biological Sciences and Engineering, Indian Institute of Technology Palakkad, Palakkad, Kerala, India
| | - Avinash Sonawane
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Simrol, Madhya Pradesh, India
| |
Collapse
|
22
|
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
|
23
|
Ogawa N, Ohta M, Ikeguchi M. Conformational Selectivity of ITK Inhibitors: Insights from Molecular Dynamics Simulations. J Chem Inf Model 2023; 63:7860-7872. [PMID: 38069816 PMCID: PMC10751800 DOI: 10.1021/acs.jcim.3c01352] [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/22/2023] [Revised: 11/24/2023] [Accepted: 11/27/2023] [Indexed: 12/26/2023]
Abstract
Interleukin-2-inducible T-cell kinase (ITK) regulates the response to T-cell receptor signaling and is a drug target for inflammatory and immunological diseases. Molecules that bind preferentially to the active form of ITK have low selectivity between kinases, whereas those that bind preferentially to the inactive form have high selectivity for ITK. Therefore, computational methods to predict the conformational selectivity of compounds are required to design highly selective ITK inhibitors. In this study, we performed absolute binding free-energy perturbation (ABFEP) simulations for 11 compounds on both active and inactive forms of ITK, and the calculated binding free energies were compared with experimental data. The conformational selectivity of 10 of the 11 compounds was correctly predicted using ABFEP. To investigate the mechanism underlying the stabilization of the active and inactive structures by the compounds, we performed extensive, conventional molecular dynamics simulations, which revealed that the compound-induced stabilization of the P-loop and linkage of conformational changes in L489, V419, F501, and M410 upon compound binding were critical factors. A guideline for designing inactive-form binders is proposed based on these key structural factors. The ABFEP and the created guidelines are expected to facilitate the discovery of highly selective ITK inhibitors.
Collapse
Affiliation(s)
- Naoki Ogawa
- Graduate
School of Medicinal Life Science, Yokohama
City University, 1-7-29, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
- Central
Pharmaceutical Research Institute, Japan
Tobacco Inc., 1-1, Murasaki-cho, Takatsuki, Osaka 569-1125, Japan
| | - Masateru Ohta
- HPC-
and AI-Driven Drug Development Platform Division, Center for Computational
Science, RIKEN, 1-7-22, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Mitsunori Ikeguchi
- Graduate
School of Medicinal Life Science, Yokohama
City University, 1-7-29, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
- HPC-
and AI-Driven Drug Development Platform Division, Center for Computational
Science, RIKEN, 1-7-22, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| |
Collapse
|
24
|
Michino M, Beautrait A, Boyles NA, Nadupalli A, Dementiev A, Sun S, Ginn J, Baxt L, Suto R, Bryk R, Jerome SV, Huggins DJ, Vendome J. Shape-Based Virtual Screening of a Billion-Compound Library Identifies Mycobacterial Lipoamide Dehydrogenase Inhibitors. ACS BIO & MED CHEM AU 2023; 3:507-515. [PMID: 38144256 PMCID: PMC10739260 DOI: 10.1021/acsbiomedchemau.3c00046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 08/25/2023] [Accepted: 08/25/2023] [Indexed: 12/26/2023]
Abstract
Lpd (lipoamide dehydrogenase) in Mycobacterium tuberculosis (Mtb) is required for virulence and is a genetically validated tuberculosis (TB) target. Numerous screens have been performed over the last decade, yet only two inhibitor series have been identified. Recent advances in large-scale virtual screening methods combined with make-on-demand compound libraries have shown the potential for finding novel hits. In this study, the Enamine REAL library consisting of ∼1.12 billion compounds was efficiently screened using the GPU Shape screen method against Mtb Lpd to find additional chemical matter that would expand on the known sulfonamide inhibitor series. We identified six new inhibitors with IC50 in the range of 5-100 μM. While these compounds remained chemically close to the already known sulfonamide series inhibitors, some diversity was found in the cores of the hits. The two most potent hits were further validated by one-step potency optimization to submicromolar levels. The co-crystal structure of optimized analogue TDI-13537 provided new insights into the potency determinants of the series.
Collapse
Affiliation(s)
- Mayako Michino
- Sanders
Tri-Institutional Therapeutics Discovery Institute, 1230 York Avenue, Box 122, New York, New York 10065, United States
| | - Alexandre Beautrait
- Schrödinger,
Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Nicholas A. Boyles
- Schrödinger,
Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Aparna Nadupalli
- Schrödinger,
Inc., 12 Michigan Dr., Natick, Massachusetts 01760, United States
| | - Alexey Dementiev
- Schrödinger,
Inc., 12 Michigan Dr., Natick, Massachusetts 01760, United States
| | - Shan Sun
- Sanders
Tri-Institutional Therapeutics Discovery Institute, 1230 York Avenue, Box 122, New York, New York 10065, United States
| | - John Ginn
- Sanders
Tri-Institutional Therapeutics Discovery Institute, 1230 York Avenue, Box 122, New York, New York 10065, United States
| | - Leigh Baxt
- Sanders
Tri-Institutional Therapeutics Discovery Institute, 1230 York Avenue, Box 122, New York, New York 10065, United States
| | - Robert Suto
- Schrödinger,
Inc., 12 Michigan Dr., Natick, Massachusetts 01760, United States
| | - Ruslana Bryk
- Department
of Microbiology and Immunology, Weill Cornell
Medicine, New York, New York 10065, United States
| | - Steven V. Jerome
- Schrödinger,
Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - David J. Huggins
- Sanders
Tri-Institutional Therapeutics Discovery Institute, 1230 York Avenue, Box 122, New York, New York 10065, United States
- Department
of Physiology and Biophysics, Weill Cornell
Medicine, New York, New York 10021, United States
| | - Jeremie Vendome
- Schrödinger,
Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| |
Collapse
|
25
|
Vögele M, Zhang BW, Kaindl J, Wang L. Is the Functional Response of a Receptor Determined by the Thermodynamics of Ligand Binding? J Chem Theory Comput 2023; 19:8414-8422. [PMID: 37943175 DOI: 10.1021/acs.jctc.3c00899] [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: 11/10/2023]
Abstract
For an effective drug, strong binding to the target protein is a prerequisite, but it is not enough. To produce a particular functional response, drugs need to either block the proteins' functions or modulate their activities by changing their conformational equilibrium. The binding free energy of a compound to its target is routinely calculated, but the timescales for the protein conformational changes are prohibitively long to be efficiently modeled via physics-based simulations. Thermodynamic principles suggest that the binding free energies of the ligands with different receptor conformations may infer their efficacy. However, this hypothesis has not been thoroughly validated. We present an actionable protocol and a comprehensive study to show that binding thermodynamics provides a strong predictor of the efficacy of a ligand. We apply the absolute binding free energy perturbation method to ligands bound to active and inactive states of eight G protein-coupled receptors and a nuclear receptor and then compare the resulting binding free energies. We find that carefully designed restraints are often necessary to efficiently model the corresponding conformational ensembles for each state. Our method achieves unprecedented performance in classifying ligands as agonists or antagonists across the various investigated receptors, all of which are important drug targets.
Collapse
Affiliation(s)
- Martin Vögele
- Schrödinger, Inc., 1540 Broadway 24th Floor, New York, New York 10036, United States
| | - Bin W Zhang
- Schrödinger, Inc., 1540 Broadway 24th Floor, New York, New York 10036, United States
| | - Jonas Kaindl
- Schrödinger GmbH, Glücksteinallee 25, Mannheim 68163, Germany
| | - Lingle Wang
- Schrödinger, Inc., 1540 Broadway 24th Floor, New York, New York 10036, United States
| |
Collapse
|
26
|
Bose S, Lotz SD, Deb I, Shuck M, Lee KSS, Dickson A. How Robust Is the Ligand Binding Transition State? J Am Chem Soc 2023; 145:25318-25331. [PMID: 37943667 PMCID: PMC11059145 DOI: 10.1021/jacs.3c08940] [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] [Indexed: 11/12/2023]
Abstract
For many drug targets, it has been shown that the kinetics of drug binding (e.g., on rate and off rate) is more predictive of drug efficacy than thermodynamic quantities alone. This motivates the development of predictive computational models that can be used to optimize compounds on the basis of their kinetics. The structural details underpinning these computational models are found not only in the bound state but also in the short-lived ligand binding transition states. Although transition states cannot be directly observed experimentally due to their extremely short lifetimes, recent successes have demonstrated that modeling the ligand binding transition state is possible with the help of enhanced sampling molecular dynamics methods. Previously, we generated unbinding paths for an inhibitor of soluble epoxide hydrolase (sEH) with a residence time of 11 min. Here, we computationally modeled unbinding events with the weighted ensemble method REVO (resampling of ensembles by variation optimization) for five additional inhibitors of sEH with residence times ranging from 14.25 to 31.75 min, with average prediction accuracy within an order of magnitude. The unbinding ensembles are analyzed in detail, focusing on features of the ligand binding transition state ensembles (TSEs). We find that ligands with similar bound poses can show significant differences in their ligand binding TSEs, in terms of their spatial distribution and protein-ligand interactions. However, we also find similarities across the TSEs when examining more general features such as ligand degrees of freedom. Together these findings show significant challenges for rational, kinetics-based drug design.
Collapse
Affiliation(s)
- Samik Bose
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Samuel D Lotz
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Indrajit Deb
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Megan Shuck
- Department of Pharmacology and Toxicology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Kin Sing Stephen Lee
- Department of Pharmacology and Toxicology, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Institute of Integrative Toxicology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Alex Dickson
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
27
|
Papadourakis M, Sinenka H, Matricon P, Hénin J, Brannigan G, Pérez-Benito L, Pande V, van Vlijmen H, de Graaf C, Deflorian F, Tresadern G, Cecchini M, Cournia Z. Alchemical Free Energy Calculations on Membrane-Associated Proteins. J Chem Theory Comput 2023; 19:7437-7458. [PMID: 37902715 PMCID: PMC11017255 DOI: 10.1021/acs.jctc.3c00365] [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/31/2023] [Indexed: 10/31/2023]
Abstract
Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.
Collapse
Affiliation(s)
- Michail Papadourakis
- Biomedical
Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| | - Hryhory Sinenka
- Institut
de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg Cedex, France
| | - Pierre Matricon
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Jérôme Hénin
- Laboratoire
de Biochimie Théorique UPR 9080, CNRS and Université Paris Cité, 75005 Paris, France
| | - Grace Brannigan
- Center
for Computational and Integrative Biology, Rutgers University−Camden, Camden, New Jersey 08103, United States of America
- Department
of Physics, Rutgers University−Camden, Camden, New Jersey 08102, United States
of America
| | - Laura Pérez-Benito
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Vineet Pande
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Herman van Vlijmen
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Chris de Graaf
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Francesca Deflorian
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Gary Tresadern
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Marco Cecchini
- Institut
de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg Cedex, France
| | - Zoe Cournia
- Biomedical
Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| |
Collapse
|
28
|
Hernández González JE, de Araujo AS. Alchemical Calculation of Relative Free Energies for Charge-Changing Mutations at Protein-Protein Interfaces Considering Fixed and Variable Protonation States. J Chem Inf Model 2023; 63:6807-6822. [PMID: 37851531 DOI: 10.1021/acs.jcim.3c00972] [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/20/2023]
Abstract
The calculation of relative free energies (ΔΔG) for charge-changing mutations at protein-protein interfaces through alchemical methods remains challenging due to variations in the system's net charge during charging steps, the possibility of mutated and contacting ionizable residues occurring in various protonation states, and undersampling issues. In this study, we present a set of strategies, collectively termed TIRST/TIRST-H+, to address some of these challenges. Our approaches combine thermodynamic integration (TI) with the prediction of pKa shifts to calculate ΔΔG values. Moreover, special sets of restraints are employed to keep the alchemically transformed molecules separated. The accuracy of the devised approaches was assessed on a large and diverse data set comprising 164 point mutations of charged residues (Asp, Glu, Lys, and Arg) to Ala at the protein-protein interfaces of complexes with known three-dimensional structures. Mean absolute and root-mean-square errors ranging from 1.38 to 1.66 and 1.89 to 2.44 kcal/mol, respectively, and Pearson correlation coefficients of ∼0.6 were obtained when testing the approaches on the selected data set using the GPU-TI module of Amber18 suite and the ff14SB force field. Furthermore, the inclusion of variable protonation states for the mutated acid residues improved the accuracy of the predicted ΔΔG values. Therefore, our results validate the use of TIRST/TIRST-H+ in prospective studies aimed at evaluating the impact of charge-changing mutations to Ala on the stability of protein-protein complexes.
Collapse
|
29
|
Anichina K, Mavrova A, Vuchev D, Popova-Daskalova G, Bassi G, Rossi A, Montesi M, Panseri S, Fratev F, Naydenova E. Benzimidazoles Containing Piperazine Skeleton at C-2 Position as Promising Tubulin Modulators with Anthelmintic and Antineoplastic Activity. Pharmaceuticals (Basel) 2023; 16:1518. [PMID: 38004384 PMCID: PMC10675210 DOI: 10.3390/ph16111518] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2023] [Revised: 10/17/2023] [Accepted: 10/20/2023] [Indexed: 11/26/2023] Open
Abstract
Benzimidazole anthelmintic drugs hold promise for repurposing as cancer treatments due to their interference with tubulin polymerization and depolymerization, manifesting anticancer properties. We explored the potential of benzimidazole compounds with a piperazine fragment at C-2 as tubulin-targeting agents. In particular, we assessed their anthelmintic activity against isolated Trichinella spiralis muscle larvae and their effects on glioblastoma (U-87 MG) and breast cancer (MDA-MB-231) cell lines. Compound 7c demonstrated exceptional anthelmintic efficacy, achieving a 92.7% reduction in parasite activity at 100 μg/mL after 48 hours. In vitro cytotoxicity analysis of MDA-MB 231 and U87 MG cell lines showed that derivatives 7b, 7d, and 7c displayed lower IC50 values compared to albendazole (ABZ), the control. These piperazine benzimidazoles effectively reduced cell migration in both cell lines, with compound 7c exhibiting the most significant reduction, making it a promising candidate for further study. The binding mode of the most promising compound 7c, was determined using the induced fit docking-molecular dynamics (IFD-MD) approach. Regular docking and IFD were also employed for comparison. The IFD-MD analysis revealed that 7c binds to tubulin in a unique binding cavity near that of ABZ, but the benzimidazole ring was fitted much deeper into the binding pocket. Finally, the absolute free energy of perturbation technique was applied to evaluate the 7c binding affinity, further confirming the observed binding mode.
Collapse
Affiliation(s)
- Kameliya Anichina
- Department of Organic Synthesis, University of Chemical Technology and Metallurgy, 8 Kliment Ohridski Blvd., 1756 Sofia, Bulgaria;
| | - Anelia Mavrova
- Department of Organic Synthesis, University of Chemical Technology and Metallurgy, 8 Kliment Ohridski Blvd., 1756 Sofia, Bulgaria;
| | - Dimitar Vuchev
- Department of Infectious Diseases, Parasitology and Tropical Medicine, Medical University, 15A Vasil Aprilov Blvd., 4002 Plovdiv, Bulgaria; (D.V.); (G.P.-D.)
| | - Galya Popova-Daskalova
- Department of Infectious Diseases, Parasitology and Tropical Medicine, Medical University, 15A Vasil Aprilov Blvd., 4002 Plovdiv, Bulgaria; (D.V.); (G.P.-D.)
| | - Giada Bassi
- Institute of Science, Technology and Sustainability for Ceramics, National Research Council of Italy, Via Granarolo 64, 48018 Faenza, Italy; (G.B.); (A.R.); (M.M.); (S.P.)
- Department of Neurosciences, Imaging and Clinical Sciences, University of G. D’Annunzio, Via Luigi Polacchi, 11, 66100 Chieti, Italy
| | - Arianna Rossi
- Institute of Science, Technology and Sustainability for Ceramics, National Research Council of Italy, Via Granarolo 64, 48018 Faenza, Italy; (G.B.); (A.R.); (M.M.); (S.P.)
- Department of Chemical, Biological, Pharmaceutical and Environmental Sciences, University of Messina, Viale Ferdinando Stagno d’Alcontres, 31, 98166 Messina, Italy
| | - Monica Montesi
- Institute of Science, Technology and Sustainability for Ceramics, National Research Council of Italy, Via Granarolo 64, 48018 Faenza, Italy; (G.B.); (A.R.); (M.M.); (S.P.)
| | - Silvia Panseri
- Institute of Science, Technology and Sustainability for Ceramics, National Research Council of Italy, Via Granarolo 64, 48018 Faenza, Italy; (G.B.); (A.R.); (M.M.); (S.P.)
| | - Filip Fratev
- Micar Innovation (Micar 21) Ltd., 34B Persenk Str., 1407 Sofia, Bulgaria;
- Department of Pharmaceutical Sciences, School of Pharmacy, The University of Texas at El Paso, 1101 N Campbell St., El Paso, TX 79968, USA
| | - Emilia Naydenova
- Department of Organic Chemistry, University of Chemical Technology and Metallurgy, 8 Kliment Ohridski Blvd., 1756 Sofia, Bulgaria;
| |
Collapse
|
30
|
Salomon K, Abramyan AM, Plenge P, Wang L, Bundgaard C, Bang-Andersen B, Loland CJ, Shi L. Dynamic extracellular vestibule of human SERT: Unveiling druggable potential with high-affinity allosteric inhibitors. Proc Natl Acad Sci U S A 2023; 120:e2304089120. [PMID: 37792512 PMCID: PMC10576121 DOI: 10.1073/pnas.2304089120] [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] [Received: 03/13/2023] [Accepted: 08/15/2023] [Indexed: 10/06/2023] Open
Abstract
The serotonin transporter (SERT) tightly regulates synaptic serotonin levels and has been the primary target of antidepressants. Binding of inhibitors to the allosteric site of human SERT (hSERT) impedes the dissociation of antidepressants bound at the central site and may enhance the efficacy of such antidepressants to potentially reduce their dosage and side effects. Here, we report the identification of a series of high-affinity allosteric inhibitors of hSERT in a unique scaffold, with the lead compound, Lu AF88273 (3-(1-(2-(1H-indol-3-yl)ethyl)piperidin-4-yl)-6-chloro-1H-indole), having 2.1 nM allosteric potency in inhibiting imipramine dissociation. In addition, we find that Lu AF88273 also inhibits serotonin transport in a noncompetitive manner. The binding pose of Lu AF88273 in the allosteric site of hSERT is determined with extensive molecular dynamics simulations and rigorous absolute binding free energy perturbation (FEP) calculations, which show that a part of the compound occupies a dynamically formed small cavity. The predicted binding location and pose are validated by site-directed mutagenesis and can explain much of the structure-activity relationship of these inhibitors using the relative binding FEP calculations. Together, our findings provide a promising lead compound and the structural basis for the development of allosteric drugs targeting hSERT. Further, they demonstrate that the divergent allosteric sites of neurotransmitter transporters can be selectively targeted.
Collapse
Affiliation(s)
- Kristine Salomon
- Laboratory for Membrane Protein Dynamics, Department of Neuroscience, Faculty of Health and Medical Sciences, University of Copenhagen, DK-2200Copenhagen N, Denmark
| | - Ara M. Abramyan
- Computational Chemistry and Molecular Biophysics Section, Molecular Targets and Medications Discovery Branch, National Institute on Drug Abuse - Intramural Research Program, NIH, Baltimore, MD21224
- Schrödinger, Inc., San Diego, CA92121
| | - Per Plenge
- Laboratory for Membrane Protein Dynamics, Department of Neuroscience, Faculty of Health and Medical Sciences, University of Copenhagen, DK-2200Copenhagen N, Denmark
| | | | - Christoffer Bundgaard
- Medicinal Chemistry and Translational DMPK, H. Lundbeck A/S, DK-2500Copenhagen-Valby, Denmark
| | - Benny Bang-Andersen
- Medicinal Chemistry and Translational DMPK, H. Lundbeck A/S, DK-2500Copenhagen-Valby, Denmark
| | - Claus J. Loland
- Laboratory for Membrane Protein Dynamics, Department of Neuroscience, Faculty of Health and Medical Sciences, University of Copenhagen, DK-2200Copenhagen N, Denmark
| | - Lei Shi
- Computational Chemistry and Molecular Biophysics Section, Molecular Targets and Medications Discovery Branch, National Institute on Drug Abuse - Intramural Research Program, NIH, Baltimore, MD21224
| |
Collapse
|
31
|
Majid MF, Mohd Zaid HF, Abd Shukur MF, Ahmad A, Jumbri K. Physicochemical properties and density functional theory calculation of octahedral UiO-66 with Bis(Trifluoromethanesulfonyl)imide ionic liquids. Heliyon 2023; 9:e20743. [PMID: 37867795 PMCID: PMC10585329 DOI: 10.1016/j.heliyon.2023.e20743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Revised: 08/24/2023] [Accepted: 10/05/2023] [Indexed: 10/24/2023] Open
Abstract
In this study, the physicochemical properties and molecular interactions between zirconium-based metal-organic framework (UiO-66) and three different ionic liquids based on bis(trifluoromethanesulfonyl)imide anion (EMIM+, BMIM+ and OMIM+) was performed via a combined experimental and computational approach. The ionic liquid loaded UiO-66 or IL@UiO-66 was synthesized and characterized to understand the host-guest interaction. Density functional theory calculation was performed to analyse the electronic structure of IL@UiO-66 to provide molecular insight on the dominant interactions occurred in the hybrid material. Results showed that all ILs were successfully incorporated into the micropores of UiO-66. The 3D framework was retained even after loaded with ILs as analyzed from XRD pattern. FTIR spectrum reveals that interactions of ILs with UiO-66 influenced by the alkyl chain length of the cation. The anion has a profound affinity with the UiO-66 due to the presence of electronegative atoms. Phase transition study from DSC suggested that the incorporation of ILs has stabilized the framework of UiO-66 by shifting the endothermic peak to a higher state. These findings were further elaborated with DFT calculation. Geometrical optimizations confirmed the structural parameter changes of UiO-66 when loaded with ILs. These was mainly contributed by the non-covalent interactions which was confirmed by the reduced density gradient scattered plot. Another important findings are the strength of hydrogen bonding at the host-guest interface was influenced by the alkyl chain length. The molecular orbital analysis also shows that the size of alkyl chain influence the reactivity of the hybrid material. The present study provides fundamental insights on the molecular interaction of UiO-66 and ILs as a hybrid material, which can open new possibilities for advanced material for metal-organic framework applications in energy storage system, catalysis, gas storage and medicinal chemistry.
Collapse
Affiliation(s)
- Mohd Faridzuan Majid
- Department of Fundamental and Applied Sciences, Universiti Teknologi PETRONAS, 32610, Seri Iskandar, Perak Darul Ridzuan, Malaysia
- Centre of Innovative Nanostructures & Nanodevices (COINN), Universiti Teknologi PETRONAS, 32610, Seri Iskandar, Perak Darul Ridzuan, Malaysia
| | - Hayyiratul Fatimah Mohd Zaid
- Chemical Engineering Department, Universiti Teknologi PETRONAS, Seri Iskandar, 32610, Perak Darul Ridzuan, Malaysia
- Centre of Innovative Nanostructures & Nanodevices (COINN), Universiti Teknologi PETRONAS, 32610, Seri Iskandar, Perak Darul Ridzuan, Malaysia
| | - Muhammad Fadhlullah Abd Shukur
- Department of Fundamental and Applied Sciences, Universiti Teknologi PETRONAS, 32610, Seri Iskandar, Perak Darul Ridzuan, Malaysia
- Centre of Innovative Nanostructures & Nanodevices (COINN), Universiti Teknologi PETRONAS, 32610, Seri Iskandar, Perak Darul Ridzuan, Malaysia
| | - Azizan Ahmad
- Department of Chemical Sciences, Universiti Kebangsaan Malaysia, Bangi, 43600, Selangor, Malaysia
- Department of Physics, Faculty of Science and Technology, Airlangga University (Campus C), Mulyorejo Road, Surabaya, 60115, Indonesia
| | - Khairulazhar Jumbri
- Department of Fundamental and Applied Sciences, Universiti Teknologi PETRONAS, 32610, Seri Iskandar, Perak Darul Ridzuan, Malaysia
- Centre for Research in Ionic Liquids (CORIL), Universiti Teknologi PETRONAS, 32610, Seri Iskandar, Perak Darul Ridzuan, Malaysia
| |
Collapse
|
32
|
Ghahremanpour MM, Saar A, Tirado-Rives J, Jorgensen WL. Computation of Absolute Binding Free Energies for Noncovalent Inhibitors with SARS-CoV-2 Main Protease. J Chem Inf Model 2023; 63:5309-5318. [PMID: 37561001 DOI: 10.1021/acs.jcim.3c00874] [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: 08/11/2023]
Abstract
Accurate, routine calculation of absolute binding free energies (ABFEs) for protein-ligand complexes remains a key goal of computer-aided drug design since it can enable screening and optimization of drug candidates. For development and testing of related methods, it is important to have high-quality datasets. To this end, from our own experimental studies, we have selected a set of 16 inhibitors of the SARS-CoV-2 main protease (Mpro) with structural diversity and well-distributed BFEs covering a 5 kcal/mol range. There is also minimal structural uncertainty since X-ray crystal structures have been deposited for 12 of the compounds. For methods testing, we report ABFE results from 2 μs molecular dynamics (MD) simulations using free energy perturbation (FEP) theory. The correlation of experimental and computed results is encouraging, with a Pearson's r2 of 0.58 and a Kendall τ of 0.24. The results indicate that current FEP-based ABFE calculations can be used for identification of active compounds (hits). While their accuracy for lead optimization is not yet sufficient, this activity remains addressable in separate lead series by relative BFE calculations.
Collapse
Affiliation(s)
| | - Anastasia Saar
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - Julian Tirado-Rives
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - William L Jorgensen
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| |
Collapse
|
33
|
Clark F, Robb G, Cole DJ, Michel J. Comparison of Receptor-Ligand Restraint Schemes for Alchemical Absolute Binding Free Energy Calculations. J Chem Theory Comput 2023; 19:3686-3704. [PMID: 37285579 PMCID: PMC10308817 DOI: 10.1021/acs.jctc.3c00139] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2023] [Indexed: 06/09/2023]
Abstract
Alchemical absolute binding free energy calculations are of increasing interest in drug discovery. These calculations require restraints between the receptor and ligand to restrict their relative positions and, optionally, orientations. Boresch restraints are commonly used, but they must be carefully selected in order to sufficiently restrain the ligand and to avoid inherent instabilities. Applying multiple distance restraints between anchor points in the receptor and ligand provides an alternative framework without inherent instabilities which may provide convergence benefits by more strongly restricting the relative movements of the receptor and ligand. However, there is no simple method to calculate the free energy of releasing these restraints due to the coupling of the internal and external degrees of freedom of the receptor and ligand. Here, a method to rigorously calculate free energies of binding with multiple distance restraints by imposing intramolecular restraints on the anchor points is proposed. Absolute binding free energies for the human macrophage migration inhibitory factor/MIF180, system obtained using a variety of Boresch restraints and rigorous and nonrigorous implementations of multiple distance restraints are compared. It is shown that several multiple distance restraint schemes produce estimates in good agreement with Boresch restraints. In contrast, calculations without orientational restraints produce erroneously favorable free energies of binding by up to approximately 4 kcal mol-1. These approaches offer new options for the deployment of alchemical absolute binding free energy calculations.
Collapse
Affiliation(s)
- Finlay Clark
- EaStCHEM
School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| | - Graeme Robb
- Oncology
R&D, AstraZeneca, Cambridge CB4 0WG, United Kingdom
| | - Daniel J. Cole
- School
of Natural and Environmental Sciences, Newcastle
University, Newcastle
upon Tyne NE1 7RU, United Kingdom
| | - Julien Michel
- EaStCHEM
School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| |
Collapse
|