1
|
Lagardère L, Maurin L, Adjoua O, El Hage K, Monmarché P, Piquemal JP, Hénin J. Lambda-ABF: Simplified, Portable, Accurate, and Cost-Effective Alchemical Free-Energy Computation. J Chem Theory Comput 2024; 20:4481-4498. [PMID: 38805379 DOI: 10.1021/acs.jctc.3c01249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
We introduce the lambda-Adaptive Biasing Force (lambda-ABF) method for the computation of alchemical free-energy differences. We propose a software implementation and showcase it on biomolecular systems. The method arises from coupling multiple-walker adaptive biasing force with λ-dynamics. The sampling of the alchemical variable is continuous and converges toward a uniform distribution, making manual optimization of the λ schedule unnecessary. Contrary to most other approaches, alchemical free-energy estimates are obtained immediately without any postprocessing. Free diffusion of λ improves orthogonal relaxation compared to fixed-λ thermodynamic integration or free-energy perturbation. Furthermore, multiple walkers provide generic orthogonal space coverage with minimal user input and negligible computational overhead. We show that our high-performance implementations coupling the Colvars library with NAMD and Tinker-HP can address real-world cases including ligand-receptor binding with both fixed-charge and polarizable models, with a demonstrably richer sampling than fixed-λ methods. The implementation is fully open-source, publicly available, and readily usable by practitioners of current alchemical methods. Thanks to the portable Colvars library, lambda-ABF presents a unified user interface regardless of the back-end (NAMD, Tinker-HP, or any software to be interfaced in the future), sparing users the effort of learning multiple interfaces. Finally, the Colvars Dashboard extension of the visual molecular dynamics (VMD) software provides an interactive monitoring and diagnostic tool for lambda-ABF simulations.
Collapse
Affiliation(s)
- Louis Lagardère
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Sorbonne Université, Institut Parisien de Chimie Physique et Théorique, FR2622 CNRS, 75005 Paris, France
- Qubit Pharmaceuticals, 29 rue du Faubourg Saint Jacques, 75014 Paris, France
| | - Lise Maurin
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Sorbonne Université, Laboratoire Jacques-Louis Lions, UMR 7589 CNRS, 75005 Paris, France
| | - Olivier Adjoua
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
| | - Krystel El Hage
- Qubit Pharmaceuticals, 29 rue du Faubourg Saint Jacques, 75014 Paris, France
| | - Pierre Monmarché
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Sorbonne Université, Laboratoire Jacques-Louis Lions, UMR 7589 CNRS, 75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Qubit Pharmaceuticals, 29 rue du Faubourg Saint Jacques, 75014 Paris, France
| | - Jérôme Hénin
- Laboratoire de Biochimie Théorique, Université Paris Cité, CNRS, UPR 9080, 75005 Paris, France
| |
Collapse
|
2
|
Ma L, Li M, Lv J, Yuan Q, Yin X, Lu W, Lin W, Wang P, Cui J, Lv Q, Liu J, Hu L. Design, synthesis, and biological evaluation of novel sesquiterpene lactone derivatives as PKM2 activators with potent anti-ulcerative colitis activities. Eur J Med Chem 2024; 272:116426. [PMID: 38718622 DOI: 10.1016/j.ejmech.2024.116426] [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: 01/10/2024] [Revised: 04/12/2024] [Accepted: 04/12/2024] [Indexed: 05/27/2024]
Abstract
Pyruvate kinase isoform 2 (PKM2) is closely related to the regulation of Th17/Treg balance, which is considered to be an effective strategy for UC therapy. Parthenolide (PTL), a natural product, only possesses moderate PKM2-activating activity. Thus, five series of PTL derivatives are designed and synthesized to improve PKM2-activated activities and anti-UC abilities. Through detailed structure optimization, B4 demonstrates potent T-cell anti-proliferation activity (IC50 = 0.43 μM) and excellent PKM2-activated ability (AC50 = 0.144 μM). Subsequently, through mass spectrometry analysis, B4 is identified to interact with Cys423 of PKM2 via covalent-bond. Molecular docking and molecular dynamic simulation results reveal that the trifluoromethoxy of B4 forms a stronger hydrophobic interaction with Ala401, Pro402, and Ile403. In addition, B4 has a significant effect only on Th17 cell differentiation, thereby regulating the Th17/Treg balance. The effect of B4 on Th17/Treg imbalance can be attributed to inhibition of PKM2 dimer translocation and suppression of glucose metabolism. Finally, B4 can notably ameliorate the symptoms of dextran sulfate sodium (DSS)-induced colitis in mouse model in vivo. Thus, B4 is confirmed as a potent PKM2 activator, and has the potential to develop as a novel anti-UC agent.
Collapse
Affiliation(s)
- Lingyu Ma
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Mengting Li
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Jiahao Lv
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Qingxin Yuan
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Xunkai Yin
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Wenyu Lu
- School of Artificial Intelligence and Information Technology, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Weijiang Lin
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Ping Wang
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Jian Cui
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Qi Lv
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China.
| | - Jian Liu
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China; School of Artificial Intelligence and Information Technology, Nanjing University of Chinese Medicine, Nanjing, 210023, China.
| | - Lihong Hu
- Jiangsu Key Laboratory for Functional Substance of Chinese Medicine, School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, 210023, China.
| |
Collapse
|
3
|
Schreiber S, Gercke D, Lenz F, Jose J. Application of an alchemical free energy method for the prediction of thermostable DuraPETase variants. Appl Microbiol Biotechnol 2024; 108:305. [PMID: 38643427 PMCID: PMC11033240 DOI: 10.1007/s00253-024-13144-z] [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: 01/16/2024] [Revised: 03/25/2024] [Accepted: 04/09/2024] [Indexed: 04/22/2024]
Abstract
Non-equilibrium (NEQ) alchemical free energy calculations are an emerging tool for accurately predicting changes in protein folding free energy resulting from amino acid mutations. In this study, this method in combination with the Rosetta ddg monomer tool was applied to predict more thermostable variants of the polyethylene terephthalate (PET) degrading enzyme DuraPETase. The Rosetta ddg monomer tool efficiently enriched promising mutations prior to more accurate prediction by NEQ alchemical free energy calculations. The relative change in folding free energy of 96 single amino acid mutations was calculated by NEQ alchemical free energy calculation. Experimental validation of ten of the highest scoring variants identified two mutations (DuraPETaseS61M and DuraPETaseS223Y) that increased the melting temperature (Tm) of the enzyme by up to 1 °C. The calculated relative change in folding free energy showed an excellent correlation with experimentally determined Tm resulting in a Pearson's correlation coefficient of r = - 0.84. Limitations in the prediction of strongly stabilizing mutations were, however, encountered and are discussed. Despite these challenges, this study demonstrates the practical applicability of NEQ alchemical free energy calculations in prospective enzyme engineering projects. KEY POINTS: • Rosetta ddg monomer enriches stabilizing mutations in a library of DuraPETase variants • NEQ free energy calculations accurately predict changes in Tm of DuraPETase • The DuraPETase variants S223Y, S42M, and S61M have increased Tm.
Collapse
Affiliation(s)
- Sebastian Schreiber
- University of Münster, Institute of Pharmaceutical and Medicinal Chemistry, PharmaCampus, Corrensstr. 48, 48149, Münster, Germany
| | - David Gercke
- University of Münster, Institute of Pharmaceutical and Medicinal Chemistry, PharmaCampus, Corrensstr. 48, 48149, Münster, Germany
| | - Florian Lenz
- University of Münster, Institute of Pharmaceutical and Medicinal Chemistry, PharmaCampus, Corrensstr. 48, 48149, Münster, Germany
| | - Joachim Jose
- University of Münster, Institute of Pharmaceutical and Medicinal Chemistry, PharmaCampus, Corrensstr. 48, 48149, Münster, Germany.
| |
Collapse
|
4
|
Wan S, Bhati AP, Wade AD, Coveney PV. Ensemble-Based Approaches Ensure Reliability and Reproducibility. J Chem Inf Model 2023; 63:6959-6963. [PMID: 37965695 PMCID: PMC10685440 DOI: 10.1021/acs.jcim.3c01654] [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: 10/13/2023] [Indexed: 11/16/2023]
Abstract
It is increasingly widely recognized that ensemble-based approaches are required to achieve reliability, accuracy, and precision in molecular dynamics calculations. The purpose of the present article is to address a frequently raised question: what is the optimal way to perform ensemble simulation to calculate quantities of interest?
Collapse
Affiliation(s)
- Shunzhou Wan
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U. K
| | - Agastya P. Bhati
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U. K
| | - Alexander D. Wade
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U. K
| | - Peter V. Coveney
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U. K
- Advanced
Research Computing Centre, University College
London, London WC1H 0AJ, U.K.
- Institute
for Informatics, Faculty of Science, University
of Amsterdam, 1098XH Amsterdam, The Netherlands
| |
Collapse
|
5
|
Wan S, Bhati AP, Coveney PV. Comparison of Equilibrium and Nonequilibrium Approaches for Relative Binding Free Energy Predictions. J Chem Theory Comput 2023; 19:7846-7860. [PMID: 37862058 PMCID: PMC10653111 DOI: 10.1021/acs.jctc.3c00842] [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/02/2023] [Indexed: 10/21/2023]
Abstract
Alchemical relative binding free energy calculations have recently found important applications in drug optimization. A series of congeneric compounds are generated from a preidentified lead compound, and their relative binding affinities to a protein are assessed in order to optimize candidate drugs. While methods based on equilibrium thermodynamics have been extensively studied, an approach based on nonequilibrium methods has recently been reported together with claims of its superiority. However, these claims pay insufficient attention to the basis and reliability of both methods. Here we report a comparative study of the two approaches across a large data set, comprising more than 500 ligand transformations spanning in excess of 300 ligands binding to a set of 14 diverse protein targets. Ensemble methods are essential to quantify the uncertainty in these calculations, not only for the reasons already established in the equilibrium approach but also to ensure that the nonequilibrium calculations reside within their domain of validity. If and only if ensemble methods are applied, we find that the nonequilibrium method can achieve accuracy and precision comparable to those of the equilibrium approach. Compared to the equilibrium method, the nonequilibrium approach can reduce computational costs but introduces higher computational complexity and longer wall clock times. There are, however, cases where the standard length of a nonequilibrium transition is not sufficient, necessitating a complete rerun of the entire set of transitions. This significantly increases the computational cost and proves to be highly inconvenient during large-scale applications. Our findings provide a key set of recommendations that should be adopted for the reliable implementation of nonequilibrium approaches to relative binding free energy calculations in ligand-protein systems.
Collapse
Affiliation(s)
- Shunzhou Wan
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U.K.
| | - Agastya P. Bhati
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U.K.
| | - Peter V. Coveney
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U.K.
- Advanced
Research Computing Centre, University College
London, London WC1H 0AJ, U.K.
- Computational
Science Laboratory, Institute for Informatics, Faculty of Science, University of Amsterdam, Amsterdam 1012 WP, Netherlands
| |
Collapse
|
6
|
Çınaroğlu S, Biggin PC. Computed Protein-Protein Enthalpy Signatures as a Tool for Identifying Conformation Sampling Problems. J Chem Inf Model 2023; 63:6095-6108. [PMID: 37759363 PMCID: PMC10565830 DOI: 10.1021/acs.jcim.3c01041] [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: 07/10/2023] [Indexed: 09/29/2023]
Abstract
Understanding the thermodynamic signature of protein-peptide binding events is a major challenge in computational chemistry. The complexity generated by both components possessing many degrees of freedom poses a significant issue for methods that attempt to directly compute the enthalpic contribution to binding. Indeed, the prevailing assumption has been that the errors associated with such approaches would be too large for them to be meaningful. Nevertheless, we currently have no indication of how well the present methods would perform in terms of predicting the enthalpy of binding for protein-peptide complexes. To that end, we carefully assembled and curated a set of 11 protein-peptide complexes where there is structural and isothermal titration calorimetry data available and then computed the absolute enthalpy of binding. The initial "out of the box" calculations were, as expected, very modest in terms of agreement with the experiment. However, careful inspection of the outliers allows for the identification of key sampling problems such as distinct conformations of peptide termini not being sampled or suboptimal cofactor parameters. Additional simulations guided by these aspects can lead to a respectable correlation with isothermal titration calorimetry (ITC) experiments (R2 of 0.88 and an RMSE of 1.48 kcal/mol overall). Although one cannot know prospectively whether computed ITC values will be correct or not, this work shows that if experimental ITC data are available, then this in conjunction with computed ITC, can be used as a tool to know if the ensemble being simulated is representative of the true ensemble or not. That is important for allowing the correct interpretation of the detailed dynamics of the system with respect to the measured enthalpy. The results also suggest that computational calorimetry is becoming increasingly feasible. We provide the data set as a resource for the community, which could be used as a benchmark to help further progress in this area.
Collapse
Affiliation(s)
| | - Philip C. Biggin
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| |
Collapse
|
7
|
Jiang W. Enhanced Configurational Sampling Approaches to Alchemical Ligand Binding Free Energy Simulations: Current Status and Challenges. J Phys Chem B 2023; 127:6835-6841. [PMID: 37499215 DOI: 10.1021/acs.jpcb.3c02020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
Ligand binding free energy simulations (LB-FES) have been routine tasks in modern drug discovery campaign. A long-standing challenge for LB-FES is the difficulty in adequately sampling nontrivial environmental reorganizations in response to ligand binding. Therefore, various enhanced configurational sampling (ECS) approaches were devised to speed up fluctuations of relevant slow degrees of freedom (SDOF) and ensure simulation convergence. However, in contrast to the achievements in parametrization, software performance, and workflow automation, efficient ECS methodology suitable for high throughput screening remains in an early stage of development. Here, a review of ECS developments with LB-FES is presented, revisiting current approaches and underlining the major technical pitfalls and challenges. This Perspective focuses on alchemical LB-FES on account of their predominant role in high throughput drug screening as well as the established partnership with ECS. The critical aspects of designing ECS approaches, from both theoretical and applied perspectives, are described. This work is intended to provide a contemporary review of the scientific, technical, and practical issues associated with the accelerating convergence of alchemical LB-FES.
Collapse
Affiliation(s)
- Wei Jiang
- Computational Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Building 240, Argonne, Illinois 60439, United States
| |
Collapse
|
8
|
Lima Silva WJ, Freitas de Freitas R. Assessing the performance of docking, FEP, and MM/GBSA methods on a series of KLK6 inhibitors. J Comput Aided Mol Des 2023:10.1007/s10822-023-00515-3. [PMID: 37378817 DOI: 10.1007/s10822-023-00515-3] [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: 03/30/2023] [Accepted: 06/21/2023] [Indexed: 06/29/2023]
Abstract
Kallikrein 6 (KLK6) is an attractive drug target for the treatment of neurological diseases and for various cancers. Herein, we explore the accuracy and efficiency of different computational methods and protocols to predict the free energy of binding (ΔGbind) for a series of 49 inhibitors of KLK6. We found that the performance of the methods varied strongly with the tested system. For only one of the three KLK6 datasets, the docking scores obtained with rDock were in good agreement (R2 ≥ 0.5) with experimental values of ΔGbind. A similar result was obtained with MM/GBSA (using the ff14SB force field) calculations based on single minimized structures. Improved binding affinity predictions were obtained with the free energy perturbation (FEP) method, with an overall MUE and RMSE of 0.53 and 0.68 kcal/mol, respectively. Furthermore, in a simulation of a real-world drug discovery project, FEP was able to rank the most potent compounds at the top of the list. These results indicate that FEP can be a promising tool for the structure-based optimization of KLK6 inhibitors.
Collapse
Affiliation(s)
- Wemenes José Lima Silva
- Centro de Ciências Naturais e Humanas, Universidade Federal do ABC, São Bernardo do Campo, SP, Brazil
| | - Renato Freitas de Freitas
- Centro de Ciências Naturais e Humanas, Universidade Federal do ABC, São Bernardo do Campo, SP, Brazil.
| |
Collapse
|
9
|
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
|
10
|
Bieniek M, Wade AD, Bhati AP, Wan S, Coveney PV. TIES 2.0: A Dual-Topology Open Source Relative Binding Free Energy Builder with Web Portal. J Chem Inf Model 2023; 63:718-724. [PMID: 36719676 PMCID: PMC9930115 DOI: 10.1021/acs.jcim.2c01596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Relative binding free energy (RBFE) calculations are widely used to aid the process of drug discovery. TIES, Thermodynamic Integration with Enhanced Sampling, is a dual-topology approach to RBFE calculations with support for NAMD and OpenMM molecular dynamics engines. The software has been thoroughly validated on publicly available datasets. Here we describe the open source software along with a web portal (https://ccs-ties.org) that enables users to perform such calculations correctly and rapidly.
Collapse
Affiliation(s)
- Mateusz
K. Bieniek
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom,School
of Natural and Environmental Sciences, Newcastle
University, Newcastle upon Tyne NE1 7RU, United
Kingdom
| | - Alexander D. Wade
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Agastya P. Bhati
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Shunzhou Wan
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Peter V. Coveney
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom,Advanced
Research Computing Centre, University College
London, London WC1H 0AJ, United
Kingdom,Institute
for Informatics, Faculty of Science, University
of Amsterdam, 1098XH Amsterdam, The Netherlands,E-mail:
| |
Collapse
|
11
|
Lee TS, Tsai HC, Ganguly A, York DM. ACES: Optimized Alchemically Enhanced Sampling. J Chem Theory Comput 2023; 19:10.1021/acs.jctc.2c00697. [PMID: 36630672 PMCID: PMC10333454 DOI: 10.1021/acs.jctc.2c00697] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
We present an alchemical enhanced sampling (ACES) method implemented in the GPU-accelerated AMBER free energy MD engine. The methods hinges on the creation of an "enhanced sampling state" by reducing or eliminating selected potential energy terms and interactions that lead to kinetic traps and conformational barriers while maintaining those terms that curtail the need to otherwise sample large volumes of phase space. For example, the enhanced sampling state might involve transforming regions of a ligand and/or protein side chain into a noninteracting "dummy state" with internal electrostatics and torsion angle terms turned off. The enhanced sampling state is connected to a real-state end point through a Hamiltonian replica exchange (HREMD) framework that is facilitated by newly developed alchemical transformation pathways and smoothstep softcore potentials. This creates a counterdiffusion of real and enhanced-sampling states along the HREMD network. The effect of a differential response of the environment to the real and enhanced-sampling states is minimized by leveraging the dual topology framework in AMBER to construct a counterbalancing HREMD network in the opposite alchemical direction with the same (or similar) real and enhanced sampling states at inverted end points. The method has been demonstrated in a series of test cases of increasing complexity where traditional MD, and in several cases alternative REST2-like enhanced sampling methods, are shown to fail. The hydration free energy for acetic acid was shown to be independent of the starting conformation, and the values for four additional edge case molecules from the FreeSolv database were shown to have a significantly closer agreement with experiment using ACES. The method was further able to handle different rotamer states in a Cdk2 ligand identified as fractionally occupied in crystal structures. Finally, ACES was applied to T4-lysozyme and demonstrated that the side chain distribution of V111χ1 could be reliably reproduced for the apo state, bound to p-xylene, and in p-xylene→ benzene transformations. In these cases, the ACES method is shown to be highly robust and superior to a REST2-like enhanced sampling implementation alone.
Collapse
Affiliation(s)
- Tai-Sung Lee
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Hsu-Chun Tsai
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Abir Ganguly
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| |
Collapse
|
12
|
Horton J, Boothroyd S, Wagner J, Mitchell JA, Gokey T, Dotson DL, Behara PK, Ramaswamy VK, Mackey M, Chodera JD, Anwar J, Mobley DL, Cole DJ. Open Force Field BespokeFit: Automating Bespoke Torsion Parametrization at Scale. J Chem Inf Model 2022; 62:5622-5633. [PMID: 36351167 PMCID: PMC9709916 DOI: 10.1021/acs.jcim.2c01153] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
The development of accurate transferable force fields is key to realizing the full potential of atomistic modeling in the study of biological processes such as protein-ligand binding for drug discovery. State-of-the-art transferable force fields, such as those produced by the Open Force Field Initiative, use modern software engineering and automation techniques to yield accuracy improvements. However, force field torsion parameters, which must account for many stereoelectronic and steric effects, are considered to be less transferable than other force field parameters and are therefore often targets for bespoke parametrization. Here, we present the Open Force Field QCSubmit and BespokeFit software packages that, when combined, facilitate the fitting of torsion parameters to quantum mechanical reference data at scale. We demonstrate the use of QCSubmit for simplifying the process of creating and archiving large numbers of quantum chemical calculations, by generating a dataset of 671 torsion scans for druglike fragments. We use BespokeFit to derive individual torsion parameters for each of these molecules, thereby reducing the root-mean-square error in the potential energy surface from 1.1 kcal/mol, using the original transferable force field, to 0.4 kcal/mol using the bespoke version. Furthermore, we employ the bespoke force fields to compute the relative binding free energies of a congeneric series of inhibitors of the TYK2 protein, and demonstrate further improvements in accuracy, compared to the base force field (MUE reduced from 0.560.390.77 to 0.420.280.59 kcal/mol and R2 correlation improved from 0.720.350.87 to 0.930.840.97).
Collapse
Affiliation(s)
- Joshua
T. Horton
- School
of Natural and Environmental Sciences, Newcastle
University, Newcastle
upon TyneNE1 7RU, United
Kingdom
| | - Simon Boothroyd
- Boothroyd
Scientific Consulting Ltd., 71-75 Shelton Street, LondonWC2H 9JQ, Greater London, United Kingdom
| | - Jeffrey Wagner
- The
Open Force Field Initiative, Open Molecular
Software Foundation, Davis, California95616, United States
| | - Joshua A. Mitchell
- The
Open Force Field Initiative, Open Molecular
Software Foundation, Davis, California95616, United States
| | - Trevor Gokey
- Department
of Chemistry, University of California, Irvine, California92697, United States
| | - David L. Dotson
- The
Open Force Field Initiative, Open Molecular
Software Foundation, Davis, California95616, United States
| | - Pavan Kumar Behara
- Department
of Pharmaceutical Sciences, University of
California, Irvine, California92697, United States
| | | | - Mark Mackey
- Cresset, New Cambridge House, Bassingbourn
Road, LitlingtonSG8 0SS, Cambridgeshire, United Kingdom
| | - John D. Chodera
- Computational
& Systems Biology Program, Sloan Kettering
Institute, Memorial Sloan Kettering Cancer Center, New
York, New York10065, United States
| | - Jamshed Anwar
- Department
of Chemistry, Lancaster University, LancasterLA1 4YW, United Kingdom
| | - David L. Mobley
- Department
of Chemistry, University of California, Irvine, California92697, United States,Department
of Pharmaceutical Sciences, University of
California, Irvine, California92697, United States
| | - Daniel J. Cole
- School
of Natural and Environmental Sciences, Newcastle
University, Newcastle
upon TyneNE1 7RU, United
Kingdom,
| |
Collapse
|
13
|
The performance of ensemble-based free energy protocols in computing binding affinities to ROS1 kinase. Sci Rep 2022; 12:10433. [PMID: 35729177 PMCID: PMC9211793 DOI: 10.1038/s41598-022-13319-6] [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: 03/16/2022] [Accepted: 05/23/2022] [Indexed: 11/08/2022] Open
Abstract
Optimization of binding affinities for compounds to their target protein is a primary objective in drug discovery. Herein we report on a collaborative study that evaluates a set of compounds binding to ROS1 kinase. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to rank the binding free energies. The predicted binding free energies from ESMACS simulations show good correlations with experimental data for subsets of the compounds. Consistent binding free energy differences are generated for TIES and ESMACS. Although an unexplained overestimation exists, we obtain excellent statistical rankings across the set of compounds from the TIES protocol, with a Pearson correlation coefficient of 0.90 between calculated and experimental activities.
Collapse
|
14
|
Wade A, Bhati AP, Wan S, Coveney PV. Alchemical Free Energy Estimators and Molecular Dynamics Engines: Accuracy, Precision, and Reproducibility. J Chem Theory Comput 2022; 18:3972-3987. [PMID: 35609233 PMCID: PMC9202356 DOI: 10.1021/acs.jctc.2c00114] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2022] [Indexed: 11/28/2022]
Abstract
The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat a myriad of diseases. In this work, we examine the computation of alchemical relative binding free energies with an eye for assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2, and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2, and NAMD3 alpha. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between them of 0.50 kcal/mol. The correlation between packages is very good, with the lowest Spearman's, Pearson's and Kendall's tau correlation coefficients being 0.92, 0.91, and 0.76, respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.
Collapse
Affiliation(s)
- Alexander
D. Wade
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, UK
| | - Agastya P. Bhati
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, UK
| | - Shunzhou Wan
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, UK
| | - Peter V. Coveney
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, UK
- Informatics
Institute, University of Amsterdam, Amsterdam 1098XH, The Netherlands
- Advanced
Research Computing Centre, University College
London, London WC1H 0AJ, UK
| |
Collapse
|
15
|
Wan S, Bhati AP, Wright DW, Wall ID, Graves AP, Green D, Coveney PV. Ensemble Simulations and Experimental Free Energy Distributions: Evaluation and Characterization of Isoxazole Amides as SMYD3 Inhibitors. J Chem Inf Model 2022; 62:2561-2570. [PMID: 35508076 PMCID: PMC9131449 DOI: 10.1021/acs.jcim.2c00255] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Optimization of binding affinities for ligands to their target protein is a primary objective in rational drug discovery. Herein, we report on a collaborative study that evaluates various compounds designed to bind to the SET and MYND domain-containing protein 3 (SMYD3). SMYD3 is a histone methyltransferase and plays an important role in transcriptional regulation in cell proliferation, cell cycle, and human carcinogenesis. Experimental measurements using the scintillation proximity assay show that the distributions of binding free energies from a large number of independent measurements exhibit non-normal properties. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to predict the binding free energies and to provide a detailed chemical insight into the nature of ligand-protein binding. Our results show that the 1-trajectory ESMACS protocol works well for the set of ligands studied here. Although one unexplained outlier exists, we obtain excellent statistical ranking across the set of compounds from the ESMACS protocol and good agreement between calculations and experiments for the relative binding free energies from the TIES protocol. ESMACS and TIES are again found to be powerful protocols for the accurate comparison of the binding free energies.
Collapse
Affiliation(s)
- Shunzhou Wan
- Centre for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U.K
| | - Agastya P Bhati
- Centre for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U.K
| | - David W Wright
- Centre for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U.K
| | - Ian D Wall
- GlaxoSmithKline, Gunnels Wood Road, Stevenage, Hertfordshire SG1 2NY, U.K
| | - Alan P Graves
- GlaxoSmithKline, 1250 South Collegeville Road, Collegeville, Pennsylvania 19426, United States
| | - Darren Green
- GlaxoSmithKline, Gunnels Wood Road, Stevenage, Hertfordshire SG1 2NY, U.K
| | - Peter V Coveney
- Centre for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, U.K.,Advanced Research Computing Centre, University College London, London WC1H 0AJ U.K.,Institute for Informatics, Faculty of Science, University of Amsterdam, 1098XH Amsterdam, The Netherlands
| |
Collapse
|