1
|
El Yaagoubi OM, Ezzemani W, Oularbi L, Samaki H, Aboudkhil S. In silico identification of 20S proteasome-β5 subunit inhibitors using structure-based virtual screening. J Biomol Struct Dyn 2024; 42:6165-6173. [PMID: 37403265 DOI: 10.1080/07391102.2023.2232041] [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: 10/16/2021] [Accepted: 06/26/2023] [Indexed: 07/06/2023]
Abstract
Proteasome inhibitors have effective anti-tumor activity in cell culture and can induce apoptosis by interfering with the degradation of cell cycle proteins. 20S Proteasome is acknowledged to be a satisfactory target that has persistent properties against the human immune defense and is obligatory for the degradation of some vital proteins. This study aimed to identify potential inhibitors against 20S proteasome, specifically the β5 subunit, using structure-based virtual screening and molecular docking to reduce the number of ligands that should be eligible for experimental assays. A total of 4961 molecules with anticancer activity were screened from the ASINEX database. The filtered compounds that showed higher docking affinity were then used in more sophisticated molecular docking simulations with AutoDock Vina for validation. Finally, six drug molecules (BDE 28974746, BDE 25657353, BDE 29746159, BDD 27844484, BDE 29746109, and BDE 29746162) exhibited highly significant interactions compared to the positive controls were retained. Among these six molecules, three molecules (BDE 28974746, BDE 25657353, and BDD 27844484) showed high binding affinity and binding energy compared with Carfilzomib and Bortezomib. Molecular simulation and dynamics studies of the top three drug molecules in each case allowed us to draw further conclusions about their stability with the β5 subunit. Computed absorption, distribution, metabolism, excretion and toxicity studies on these derivatives showed encouraging results with very low toxicity, distribution, and absorption. These compounds may serve as potential hits for further biological evaluation in the development of new proteasome inhibitors.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Ouadie Mohamed El Yaagoubi
- Laboratory of Biochemistry, Environment and Agri-Food (URAC 36), Faculty of Sciences and Techniques-Mohammedia, Hassan II University of Casablanca, Morocco
| | - Wahiba Ezzemani
- Virology Unit, Viral Hepatitis Laboratory, Institut Pasteur du Maroc, Casablanca, Morocco
- Laboratoire de Biologie et Santé (URAC34), Département de Biologie, Faculté des Sciences Ben Msik, Hassan II University of Casablanca, Morocco
| | - Larbi Oularbi
- Laboratory of Materials Membranes and Environment, Faculty of Sciences and Techniques-Mohammedia, Hassan II University of Casablanca, Morocco
- Supramolecular Nanomaterials Group (SNG), Mohammed VI Polytechnic University (UM6P), Lot 660, HayMoulayRachid, BenGuerir, Morocco
| | - Hamid Samaki
- National Institute of Social Action (INAS), Tangier, Morocco
| | - Souad Aboudkhil
- Laboratory of Biochemistry, Environment and Agri-Food (URAC 36), Faculty of Sciences and Techniques-Mohammedia, Hassan II University of Casablanca, Morocco
| |
Collapse
|
2
|
Khadka D, Jayasinghe-Arachchige VM, Prabhakar R, Ramamurthy V. Application of molecular dynamic simulations in modeling the excited state behavior of confined molecules. Photochem Photobiol Sci 2023:10.1007/s43630-023-00486-2. [PMID: 37843722 DOI: 10.1007/s43630-023-00486-2] [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: 06/27/2023] [Accepted: 09/20/2023] [Indexed: 10/17/2023]
Abstract
Relative to isotropic organic solvent medium, the structure and conformation of a reactant molecule in an organized and confining medium are often different. In addition, because of the rigidity of the immediate environment, the reacting molecule have a little freedom to undergo large changes even upon gaining energy or modifications in the electronic structure. These alterations give rise to differences in the photochemistry of a molecular and supramolecular species. In this study, one such example is presented. α-Alkyl dibenzylketones upon excitation in isotropic solvents give products via Norrish type I and type II reactions that are independent of the chain length of the alkyl substituent. On the other hand, when these molecules are enclosed within an organic capsule of volume ~ 550 Å3, they give products that are strikingly dependent on the length of the α-alkyl substitution. These previously reported experimental observations are rationalized based on the structures generated by molecular modeling (docking and molecular dynamics (MD) simulations). It is shown that MD simulations that are utilized extensively in biologically important macromolecules can also be useful to understand the excited state behavior of reactive molecules that are part of supramolecular assemblies. These simulations can provide structural information of the reactant molecule and the surroundings complementing that with the one obtained from 1 and 2D NMR experiments. MD simulated structures of seven α-alkyl dibenzylketones encapsulated within the octa acid capsule provide a clear understanding of their unique behavior in this restricted medium. Because of the rigidity of the medium, these structures although generated in the ground state can rationalize the photochemical behavior of the molecules in the excited state.
Collapse
Affiliation(s)
- Dipendra Khadka
- Department of Chemistry, University of Miami, Coral Gables, FL, 33124, USA
| | | | - Rajeev Prabhakar
- Department of Chemistry, University of Miami, Coral Gables, FL, 33124, USA.
| | | |
Collapse
|
3
|
Sun Z, He Q, Gong Z, Kalhor P, Huai Z, Liu Z. A General Picture of Cucurbit[8]uril Host–Guest Binding: Recalibrating Bonded Interactions. Molecules 2023; 28:molecules28073124. [PMID: 37049887 PMCID: PMC10095826 DOI: 10.3390/molecules28073124] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2023] [Revised: 03/15/2023] [Accepted: 03/27/2023] [Indexed: 04/03/2023] Open
Abstract
Atomic-level understanding of the dynamic feature of host–guest interactions remains a central challenge in supramolecular chemistry. The remarkable guest binding behavior of the Cucurbiturils family of supramolecular containers makes them promising drug carriers. Among Cucurbit[n]urils, Cucurbit[8]uril (CB8) has an intermediate portal size and cavity volume. It can exploit almost all host–guest recognition motifs formed by this host family. In our previous work, an extensive computational investigation of the binding of seven commonly abused and structurally diverse drugs to the CB8 host was performed, and a general dynamic binding picture of CB8-guest interactions was obtained. Further, two widely used fixed-charge models for drug-like molecules were investigated and compared in great detail, aiming at providing guidelines in choosing an appropriate charge scheme in host-guest modelling. Iterative refitting of atomic charges leads to improved binding thermodynamics and the best root-mean-squared deviation from the experimental reference is 2.6 kcal/mol. In this work, we focus on a thorough evaluation of the remaining parts of classical force fields, i.e., the bonded interactions. The widely used general Amber force fields are assessed and refitted with generalized force-matching to improve the intra-molecular conformational preference, and thus the description of inter-molecular host–guest interactions. The interaction pattern and binding thermodynamics show a significant dependence on the modelling parameters. The refitted system-specific parameter set improves the consistency of the modelling results and the experimental reference significantly. Finally, combining the previous charge-scheme comparison and the current force-field refitting, we provide general guidelines for the theoretical modelling of host–guest binding.
Collapse
|
4
|
Temel M, Tayfuroglu O, Kocak A. The performance of ANI-ML potentials for ligand-n(H 2 O) interaction energies and estimation of hydration free energies from end-point MD simulations. J Comput Chem 2023; 44:559-569. [PMID: 36324248 DOI: 10.1002/jcc.27022] [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: 07/15/2022] [Revised: 09/07/2022] [Accepted: 10/05/2022] [Indexed: 11/06/2022]
Abstract
Here, we investigate the performance of "Accurate NeurAl networK engINe for Molecular Energies" (ANI), trained on small organic compounds, on bulk systems including non-covalent interactions and applicability to estimate solvation (hydration) free energies using the interaction between the ligand and explicit solvent (water) from single-step MD simulations. The method is adopted from ANI using the Atomic Simulation Environment (ASE) and predicts the non-covalent interaction energies at the accuracy of wb97x/6-31G(d) level by a simple linear scaling for the conformations sampled by molecular dynamics (MD) simulations of ligand-n(H2 O) systems. For the first time, we test ANI potentials' abilities to reproduce solvation free energies using linear interaction energy (LIE) formulism by modifying the original LIE equation. Our results on ~250 different complexes show that the method can be accurate and have a correlation of R2 = 0.88-0.89 (MAE <1.0 kcal/mol) to the experimental solvation free energies, outperforming current end-state methods. Moreover, it is competitive to other conventional free energy methods such as FEP and BAR with 15-20 × fold reduced computational cost.
Collapse
Affiliation(s)
- Mütesir Temel
- Department of Chemistry, Gebze Technical University, Kocaeli, Turkey
| | - Omer Tayfuroglu
- Department of Chemistry, Gebze Technical University, Kocaeli, Turkey
| | - Abdulkadir Kocak
- Department of Chemistry, Gebze Technical University, Kocaeli, Turkey
| |
Collapse
|
5
|
Dutkiewicz Z. Computational methods for calculation of protein-ligand binding affinities in structure-based drug design. PHYSICAL SCIENCES REVIEWS 2022. [DOI: 10.1515/psr-2020-0034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Abstract
Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.
Collapse
Affiliation(s)
- Zbigniew Dutkiewicz
- Department of Chemical Technology of Drugs , Poznan University of Medical Sciences , ul. Grunwaldzka 6 , 60-780 Poznań , Poznan , 60-780, Poland
| |
Collapse
|
6
|
Sun Z, He Q. Seeding the multi-dimensional nonequilibrium pulling for Hamiltonian variation: indirect nonequilibrium free energy simulations at QM levels. Phys Chem Chem Phys 2022; 24:8800-8819. [PMID: 35352744 DOI: 10.1039/d2cp00355d] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The combination of free energy simulations in the alchemical and configurational spaces provides a feasible route to access the thermodynamic profiles under a computationally demanding target Hamiltonian. Normally, due to the significant differences between the computational cost of ab initio quantum mechanics (QM) calculations and those of semi-empirical quantum mechanics (SQM) and molecular mechanics (MM), this indirect method could be used to obtain the QM thermodynamics by combining the SQM or MM results and the SQM-to-QM or MM-to-QM corrections. In our previous work, a multi-dimensional nonequilibrium pulling framework for Hamiltonian variations was introduced based on bidirectional pulling and bidirectional reweighting. The method performs nonequilibrium free energy simulations in the configurational space to obtain the thermodynamic profile along the conformational change pathway under a selected computationally efficient Hamiltonian, and uses the nonequilibrium alchemical method to correct or perturb the thermodynamic profile to that under the target Hamiltonian. The BAR-based method is designed to achieve the best generality and transferability and thus leads to modest (∼20 fold) speedup. In this work, we explore the possibility of further accelerating the nonequilibrium free energy simulation by employing unidirectional pulling and using the selection criterion to obtain the initial configurations used to initiate nonequilibrium trajectories following the idea of adaptive steered molecular dynamics (ASMD). A single initial condition is used to seed the whole multi-dimensional nonequilibrium free energy simulation and the sampling is performed fully in the nonequilibrium ensemble. Introducing very short ps-length equilibrium sampling to grab more initial seeds could also be helpful. The ASMD scheme estimates the free energy difference with the unidirectional exponential average (EXP), but it does not follow exactly the requirements of the EXP estimator. Another deficiency of the seeding simulation is the inherently sequential or serial pulling due to the inter-segment dependency, which triggers some problems in the parallelizability of the simulation. Numerical tests are performed to grasp some insights and guidelines for using this selection-criterion-based ASMD scheme. The presented selection-criterion-based multi-dimensional ASMD scheme follows the same perturbation network of the BAR-based method, and thus could be used in various Hamiltonian-variation cases.
Collapse
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China.
| | - Qiaole He
- AI Department of Enzymaster (Ningbo) Bio-Engineering Co., Ltd, North Century Avenue 333, 315100 Ningbo, China
| |
Collapse
|
7
|
Sun Z, Huai Z, He Q, Liu Z. A General Picture of Cucurbit[8]uril Host-Guest Binding. J Chem Inf Model 2021; 61:6107-6134. [PMID: 34818004 DOI: 10.1021/acs.jcim.1c01208] [Citation(s) in RCA: 32] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Describing, understanding, and designing complex interaction networks within macromolecular systems remain challenging in modern chemical research. Host-guest systems, despite their relative simplicity in both the structural feature and interaction patterns, still pose problems in theoretical modeling. The barrel-shaped supramolecular container cucurbit[8]uril (CB8) shows promising functionalities in various areas, e.g., catalysis and molecular recognition. It can stably coordinate a series of structurally diverse guests with high affinities. In this work, we examine the binding of seven commonly abused drugs to the CB8 host, aiming at providing a general picture of CB8-guest binding. Extensive sampling of the configurational space of these host-guest systems is performed, and the binding pathway and interaction patterns of CB8-guest complexes are investigated. A thorough comparison of widely used fixed-charge models for drug-like molecules is presented. Iterative refitting of the atomic charges suggests significant conformation dependence of charge generation. The initial model generated at the original conformation could be inaccurate for new conformations explored during conformational search, and the newly fitted charge set improves the prediction-experiment correlation significantly. Our investigations of the configurational space of CB8-drug complexes suggest that the host-guest interactions are more complex than expected. Despite the structural simplicities of these molecules, the conformational fluctuations of the host and the guest molecules and orientations of functional groups lead to the existence of an ensemble of binding modes. The insights of the binding thermodynamics, performance of fixed-charge models, and binding patterns of the CB8-guest systems are useful for studying and elucidating the binding mechanism of other host-guest complexes.
Collapse
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Zhe Huai
- XtalPi-AI Research Center (XARC), 9F, Tower A, Dongsheng Building, No. 8, Zhongguancun East Road, Haidian District, Beijing 100083, P.R. China
| | - Qiaole He
- AI Department of Enzymaster (Ningbo) Bio-Engineering Co., Ltd., North Century Avenue 333, Ningbo 315100, China
| | - Zhirong Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
8
|
Rizzi A, Carloni P, Parrinello M. Targeted Free Energy Perturbation Revisited: Accurate Free Energies from Mapped Reference Potentials. J Phys Chem Lett 2021; 12:9449-9454. [PMID: 34555284 DOI: 10.1021/acs.jpclett.1c02135] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We present an approach that extends the theory of targeted free energy perturbation (TFEP) to calculate free energy differences and free energy surfaces at an accurate quantum mechanical level of theory from a cheaper reference potential. The convergence is accelerated by a mapping function that increases the overlap between the target and the reference distributions. Building on recent work, we show that this map can be learned with a normalizing flow neural network, without requiring simulations with the expensive target potential but only a small number of single-point calculations, and, crucially, avoiding the systematic error that was found previously. We validate the method by numerically evaluating the free energy difference in a system with a double-well potential and by describing the free energy landscape of a simple chemical reaction in the gas phase.
Collapse
Affiliation(s)
- Andrea Rizzi
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
- Atomistic Simulations, Italian Institute of Technology, Via Morego 30, Genova 16163, Italy
| | - Paolo Carloni
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
- Molecular Neuroscience and Neuroimaging (INM-11), Forschungszentrum Jülich GmbH, Jülich 52428, Germany
- Department of Physics and Universitätsklinikum, RWTH Aachen University, Aachen 52074, Germany
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Via Morego 30, Genova 16163, Italy
| |
Collapse
|
9
|
Sun Z, Liu Z. BAR‐Based Multi‐Dimensional Nonequilibrium Pulling for Indirect Construction of QM/MM Free Energy Landscapes: Varying the QM Region. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202100185] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| | - Zhirong Liu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| |
Collapse
|
10
|
Rocha-Santos A, Chaves EJ, Grillo IB, de Freitas AS, Araújo DAM, Rocha GB. Thermochemical and Quantum Descriptor Calculations for Gaining Insight into Ricin Toxin A (RTA) Inhibitors. ACS OMEGA 2021; 6:8764-8777. [PMID: 33842748 PMCID: PMC8027999 DOI: 10.1021/acsomega.0c02588] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/01/2020] [Accepted: 12/30/2020] [Indexed: 05/03/2023]
Abstract
In this work, we performed a study to assess the interactions between the ricin toxin A (RTA) subunit of ricin and some of its inhibitors using modern semiempirical quantum chemistry and ONIOM quantum mechanics/molecular mechanics (QM/MM) methods. Two approaches were followed (calculation of binding enthalpies, ΔH bind, and reactivity quantum chemical descriptors) and compared with the respective half-maximal inhibitory concentration (IC50) experimental data, to gain insight into RTA inhibitors and verify which quantum chemical method would better describe RTA-ligand interactions. The geometries for all RTA-ligand complexes were obtained after running classical molecular dynamics simulations in aqueous media. We found that single-point energy calculations of ΔH bind with the PM6-DH+, PM6-D3H4, and PM7 semiempirical methods and ONIOM QM/MM presented a good correlation with the IC50 data. We also observed, however, that the correlation decreased significantly when we calculated ΔH bind after full-atom geometry optimization with all semiempirical methods. Based on the results from reactivity descriptors calculations for the cases studied, we noted that both types of interactions, molecular overlap and electrostatic interactions, play significant roles in the overall affinity of these ligands for the RTA binding pocket.
Collapse
Affiliation(s)
- Acassio Rocha-Santos
- Department
of Chemistry, Federal University of Paraíba, Cidade Universitária, João Pessoa, PB 58051-900, Brazil
| | - Elton José
Ferreira Chaves
- Department
of Biotechnology, Federal University of
Paraíba, Cidade Universitária, João Pessoa, PB 58051-900, Brazil
| | - Igor Barden Grillo
- Department
of Chemistry, Federal University of Paraíba, Cidade Universitária, João Pessoa, PB 58051-900, Brazil
| | - Amanara Souza de Freitas
- Department
of Chemical Engineering, Federal University
of Paraíba, Cidade Universitária, João Pessoa, PB 58051-900, Brazil
| | | | - Gerd Bruno Rocha
- Department
of Chemistry, Federal University of Paraíba, Cidade Universitária, João Pessoa, PB 58051-900, Brazil
- . Phone/Fax: +55-83-3216-7437
| |
Collapse
|
11
|
Abstract
QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.
Collapse
Affiliation(s)
- Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Tanmoy Pal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Luke Xie
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
12
|
Lee TS, Allen BK, Giese TJ, Guo Z, Li P, Lin C, McGee TD, Pearlman DA, Radak BK, Tao Y, Tsai HC, Xu H, Sherman W, York DM. Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery. J Chem Inf Model 2020; 60:5595-5623. [PMID: 32936637 PMCID: PMC7686026 DOI: 10.1021/acs.jcim.0c00613] [Citation(s) in RCA: 175] [Impact Index Per Article: 43.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Predicting protein-ligand binding affinities and the associated thermodynamics of biomolecular recognition is a primary objective of structure-based drug design. Alchemical free energy simulations offer a highly accurate and computationally efficient route to achieving this goal. While the AMBER molecular dynamics package has successfully been used for alchemical free energy simulations in academic research groups for decades, widespread impact in industrial drug discovery settings has been minimal because of the previous limitations within the AMBER alchemical code, coupled with challenges in system setup and postprocessing workflows. Through a close academia-industry collaboration we have addressed many of the previous limitations with an aim to improve accuracy, efficiency, and robustness of alchemical binding free energy simulations in industrial drug discovery applications. Here, we highlight some of the recent advances in AMBER20 with a focus on alchemical binding free energy (BFE) calculations, which are less computationally intensive than alternative binding free energy methods where full binding/unbinding paths are explored. In addition to scientific and technical advances in AMBER20, we also describe the essential practical aspects associated with running relative alchemical BFE calculations, along with recommendations for best practices, highlighting the importance not only of the alchemical simulation code but also the auxiliary functionalities and expertise required to obtain accurate and reliable results. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative BFE simulations in AMBER20, with a focus on real-world drug discovery applications.
Collapse
Affiliation(s)
- Tai-Sung Lee
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Bryce K. Allen
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Timothy J. Giese
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Zhenyu Guo
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Pengfei Li
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Charles Lin
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - T. Dwight McGee
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - David A. Pearlman
- QSimulate Incorporated, Cambridge, Massachusetts 02139, United States
| | - Brian K. Radak
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Yujun Tao
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Hsu-Chun Tsai
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Huafeng Xu
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Woody Sherman
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Darrin M. York
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| |
Collapse
|
13
|
Ito S, Cui Q. Multi-level free energy simulation with a staged transformation approach. J Chem Phys 2020; 153:044115. [PMID: 32752685 DOI: 10.1063/5.0012494] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Combining multiple levels of theory in free energy simulations to balance computational accuracy and efficiency is a promising approach for studying processes in the condensed phase. While the basic idea has been proposed and explored for quite some time, it remains challenging to achieve convergence for such multi-level free energy simulations as it requires a favorable distribution overlap between different levels of theory. Previous efforts focused on improving the distribution overlap by either altering the low-level of theory for the specific system of interest or ignoring certain degrees of freedom. Here, we propose an alternative strategy that first identifies the degrees of freedom that lead to gaps in the distributions of different levels of theory and then treats them separately with either constraints or restraints or by introducing an intermediate model that better connects the low and high levels of theory. As a result, the conversion from the low level to the high level model is done in a staged fashion that ensures a favorable distribution overlap along the way. Free energy components associated with different steps are mostly evaluated explicitly, and thus, the final result can be meaningfully compared to the rigorous free energy difference between the two levels of theory with limited and well-defined approximations. The additional free energy component calculations involve simulations at the low level of theory and therefore do not incur high computational costs. The approach is illustrated with two simple but non-trivial solution examples, and factors that dictate the reliability of the result are discussed.
Collapse
Affiliation(s)
- Shingo Ito
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | - Qiang Cui
- Departments of Chemistry, Physics and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| |
Collapse
|
14
|
Pecina A, Eyrilmez SM, Köprülüoğlu C, Miriyala VM, Lepšík M, Fanfrlík J, Řezáč J, Hobza P. SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design. Chempluschem 2020; 85:2362-2371. [PMID: 32609421 DOI: 10.1002/cplu.202000120] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2020] [Revised: 05/27/2020] [Indexed: 12/17/2022]
Abstract
Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.
Collapse
Affiliation(s)
- Adam Pecina
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Saltuk M Eyrilmez
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Cemal Köprülüoğlu
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Vijay Madhav Miriyala
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Martin Lepšík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| |
Collapse
|
15
|
Sylvetsky N. Toward Simple, Predictive Understanding of Protein-Ligand Interactions: Electronic Structure Calculations on Torpedo Californica Acetylcholinesterase Join Forces with the Chemist's Intuition. Sci Rep 2020; 10:9218. [PMID: 32513975 PMCID: PMC7280257 DOI: 10.1038/s41598-020-65984-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2020] [Accepted: 05/13/2020] [Indexed: 11/09/2022] Open
Abstract
Contemporary efforts for empirically-unbiased modeling of protein-ligand interactions entail a painful tradeoff - as reliable information on both noncovalent binding factors and the dynamic behavior of a protein-ligand complex is often beyond practical limits. We demonstrate that information drawn exclusively from static molecular structures can be used for reproducing and predicting experimentally-measured binding affinities for protein-ligand complexes. In particular, inhibition constants (Ki) were calculated for seven different competitive inhibitors of Torpedo californica acetylcholinesterase using a multiple-linear-regression-based model. The latter, incorporating five independent variables - drawn from QM cluster, DLPNO-CCSD(T) calculations and LED analyses on the seven complexes, each containing active amino-acid residues found within interacting distance (3.5 Å) from the corresponding ligand - is shown to recover 99.9% of the sum of squares for measured Ki values, while having no statistically-significant residual errors. Despite being fitted to a small number of data points, leave-one-out cross-validation statistics suggest that it possesses surprising predictive value (Q2LOO=0.78, or 0.91 upon removal of a single outlier). This thus challenges ligand-invariant definitions of active sites, such as implied in the lock-key binding theory, as well as in alternatives highlighting shape-complementarity without taking electronic effects into account. Broader implications of the current work are discussed in dedicated appendices.
Collapse
Affiliation(s)
- Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001, Rehovot, Israel.
| |
Collapse
|
16
|
Sun Z, Wang X, Zhang JZH. Theoretical understanding of the thermodynamics and interactions in transcriptional regulator TtgR-ligand binding. Phys Chem Chem Phys 2019; 22:1511-1524. [PMID: 31872826 DOI: 10.1039/c9cp05980f] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The transcriptional regulator TtgR belongs to the TetR family of transcriptional repressors. It depresses the transcription of the TtgABC operon and itself and thus regulates the extrusion of noxious chemicals with efflux pumps in bacterial cells. As the ligand-binding domain of TtgR is rather flexible, it can bind with a number of structurally diverse ligands, such as antibiotics, flavonoids and aromatic solvents. In the current work, we perform equilibrium and nonequilibrium alchemical free energy simulation to predict the binding affinities of a series of ligands targeting the TtgR protein and an agreement between the theoretical prediction and the experimental result is observed. End-point methods MM/PBSA and MM/GBSA are also employed for comparison. We further study the interaction maps and contacts between the protein and the ligand and identify important interactions in the protein-ligand binding cases. The dynamics fluctuation and secondary structures are also investigated. The current work sheds light on atomic and thermodynamic understanding of the TtgR-ligand interactions.
Collapse
Affiliation(s)
- Zhaoxi Sun
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich 52425, Germany. and State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xiaohui Wang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China and Institute of Computational Science, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
| | - John Z H Zhang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China and NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China and Department of Chemistry, New York University, NY, NY 10003, USA.
| |
Collapse
|
17
|
Sun Z. BAR-based multi-dimensional nonequilibrium pulling for indirect construction of QM/MM free energy landscapes: from semi-empirical to ab initio. Phys Chem Chem Phys 2019; 21:21942-21959. [PMID: 31552953 DOI: 10.1039/c9cp04113c] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
The indirect method for the construction of quantum mechanics (QM)/molecular mechanics (MM) free energy landscapes provides a cheaper alternative for free energy simulations at the QM level. The indirect method features a direct calculation of the free energy profile with a computationally efficient but less accurate Hamiltonian (i.e. low-level Hamiltonian) and a low-level-to-high-level correction. In the thermodynamic cycle, the direct low-level calculation along the physically meaningful reaction coordinate is corrected via the alchemical method, which is often achieved with perturbation-based techniques. In our previous work, a multi-dimensional nonequilibrium pulling framework is proposed for the indirect construction of QM/MM free energy landscapes. Previously, we focused on obtaining semi-empirical QM (SQM) results indirectly from direct MM simulations and MM to SQM corrections. In this work, we apply this method to obtain results under ab initio QM Hamiltonians by combining direct SQM results and SQM to QM corrections. A series of SQM and QM Hamiltonians are benchmarked. It is observed that PM6 achieves the best performance among the low-level Hamiltonians. Therefore, we recommend using PM6 as the low-level theory in the indirect free energy simulation. Considering its higher similarity to the high-level Hamiltonians, PM6 corrected with the bond charge correction could be more accurate than the existing AM1-BCC model. Another central result in the current work is a basic protocol of choosing the strength of restraints and an appropriate time step in nonequilibrium free energy simulation at the stiff spring limit. We provide theoretical derivations to emphasize the importance of using a sufficiently large force constant and choosing an appropriate time step. It is worth noting that a general rule of thumb for choosing the time step, according to our derivation, is that a time step of 1 fs or smaller should be used, as long as the stiff spring approximation is employed, even in simulations with constraints on bonds involving hydrogen atoms.
Collapse
Affiliation(s)
- Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| |
Collapse
|
18
|
Obeidat A. Free energy of formation of SPC/E-water and TIP4P-water using BAR and TI in MD and MC. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.111274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
19
|
Giese TJ, York DM. Development of a Robust Indirect Approach for MM → QM Free Energy Calculations That Combines Force-Matched Reference Potential and Bennett's Acceptance Ratio Methods. J Chem Theory Comput 2019; 15:5543-5562. [PMID: 31507179 DOI: 10.1021/acs.jctc.9b00401] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We use the PBE0/6-31G* density functional method to perform ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations under periodic boundary conditions with rigorous electrostatics using the ambient potential composite Ewald method in order to test the convergence of MM → QM/MM free energy corrections for the prediction of 17 small-molecule solvation free energies and eight ligand binding free energies to T4 lysozyme. The "indirect" thermodynamic cycle for calculating free energies is used to explore whether a series of reference potentials improve the statistical quality of the predictions. Specifically, we construct a series of reference potentials that optimize a molecular mechanical (MM) force field's parameters to reproduce the ab initio QM/MM forces from a QM/MM simulation. The optimizations form a systematic progression of successively expanded parameters that include bond, angle, dihedral, and charge parameters. For each reference potential, we calculate benchmark quality reference values for the MM → QM/MM correction by performing the mixed MM and QM/MM Hamiltonians at 11 intermediate states, each for 200 ps. We then compare forward and reverse application of Zwanzig's relation, thermodynamic integration (TI), and Bennett's acceptance ratio (BAR) methods as a function of reference potential, simulation time, and the number of simulated intermediate states. We find that Zwanzig's equation is inadequate unless a large number of intermediate states are explicitly simulated. The TI and BAR mean signed errors are very small even when only the end-state simulations are considered, and the standard deviations of the TI and BAR errors are decreased by choosing a reference potential that optimizes the bond and angle parameters. We find a robust approach for the data sets of fairly rigid molecules considered here is to use bond + angle reference potential together with the end-state-only BAR analysis. This requires QM/MM simulations to be performed in order to generate reference data to parametrize the bond + angle reference potential, and then this same simulation serves a dual purpose as the full QM/MM end state. The convergence of the results with respect to time suggests that computational resources may be used more efficiently by running multiple simulations for no more than 50 ps, rather than running one long simulation.
Collapse
Affiliation(s)
- Timothy J Giese
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854-8087 , United States
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854-8087 , United States
| |
Collapse
|
20
|
Hudson PS, Woodcock HL, Boresch S. Use of Interaction Energies in QM/MM Free Energy Simulations. J Chem Theory Comput 2019; 15:4632-4645. [PMID: 31142113 DOI: 10.1021/acs.jctc.9b00084] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The use of the most accurate (i.e., QM or QM/MM) levels of theory for free energy simulations (FES) is typically not possible. Primarily, this is because the computational cost associated with the extensive configurational sampling needed for converging FES is prohibitive. To ensure the feasibility of QM-based FES, the "indirect" approach is generally taken, necessitating a free energy calculation between the MM and QM/MM potential energy surfaces. Ideally, this step is performed with standard free energy perturbation (Zwanzig's equation) as it only requires simulations be carried out at the low level of theory; however, work from several groups over the past few years has conclusively shown that Zwanzig's equation is ill-suited to this task. As such, many approximations have arisen to mitigate difficulties with Zwanzig's equation. One particularly popular notion is that the convergence of Zwanzig's equation can be improved by using interaction energy differences instead of total energy differences. Although problematic numerical fluctuations (a major problem when using Zwanzig's equation) are indeed reduced, our results and analysis demonstrate that this "interaction energy approximation" (IEA) is theoretically incorrect, and the implicit approximation invoked is spurious at best. Herein, we demonstrate this via solvation free energy calculations using IEA from two different low levels of theory to the same target high level. Results from this proof-of-concept consistently yield the wrong results, deviating by ∼1.5 kcal/mol from the rigorously obtained value.
Collapse
Affiliation(s)
- Phillip S Hudson
- Department of Chemistry , University of South Florida , 4202 East Fowler Avenue, CHE205 , Tampa , Florida 33620-5250 , United States.,Laboratory of Computational Biology , National Institutes of Health, National Heart, Lung and Blood Institute , 12 South Drive, Rm 3053 , Bethesda , Maryland 20892-5690 , United States
| | - H Lee Woodcock
- Department of Chemistry , University of South Florida , 4202 East Fowler Avenue, CHE205 , Tampa , Florida 33620-5250 , United States
| | - Stefan Boresch
- Faculty of Chemistry, Department of Computational Biological Chemistry , University of Vienna , Währingerstraße 17 , Vienna A-1090 , Austria
| |
Collapse
|
21
|
Wang X, He Q, Sun Z. BAR-based multi-dimensional nonequilibrium pulling for indirect construction of a QM/MM free energy landscape. Phys Chem Chem Phys 2019; 21:6672-6688. [PMID: 30855611 DOI: 10.1039/c8cp07012a] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Construction of free energy landscapes at the quantum mechanics (QM) level is computationally demanding. As shown in previous studies, by employing an indirect scheme (i.e. constructing a thermodynamic cycle connecting QM states via an alchemical pathway), simulations are converged with much less computational burden. The indirect scheme makes QM/molecular mechanics (MM) free energy simulation orders of magnitude faster than the direct QM/MM schemes. However, the indirect QM/MM simulations were mostly equilibrium sampling based and the nonequilibrium methods were merely exploited in one-dimensional alchemical QM/MM end-state correction at two end states. In this work, we represent a multi-dimensional nonequilibrium pulling scheme for indirect QM/MM free energy simulations, where the whole free energy simulation is performed only with nonequilibrium methods. The collective variable (CV) space we explore is a combination of one alchemical CV and one physically meaningful CV. The current nonequilibrium indirect QM/MM simulation method can be seen as the generalization of equilibrium perturbation based indirect QM/MM methods. The test systems include one backbone dihedral case and one distance case. The two cases are significantly different in size, enabling us to investigate the dependence of the speedup of the indirect scheme on the size of the system. It is shown that the speedup becomes larger when the size of the system becomes larger, which is consistent with the scaling behavior of QM Hamiltonians.
Collapse
Affiliation(s)
- Xiaohui Wang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | | | | |
Collapse
|
22
|
Wang M, Mei Y, Ryde U. Host-Guest Relative Binding Affinities at Density-Functional Theory Level from Semiempirical Molecular Dynamics Simulations. J Chem Theory Comput 2019; 15:2659-2671. [PMID: 30811192 DOI: 10.1021/acs.jctc.8b01280] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Relative free energies for the binding of nine cyclic carboxylate ligands to the octa-acid deep-cavity host were calculated at the combined density-functional theory and molecular mechanics (DFT/MM) level of theory. The DFT calculations employed the BLYP functional and the 6-31G* basis set for the ligand. We employed free-energy perturbations (FEP) with the reference-potential approach and used molecular dynamics (MD) simulations with the semiempirical quantum mechanical (SQM) PM6-DH+ method for the ligand as an intermediate level between MM and DFT/MM to improve the convergence. Thus, the relative binding free energy of two ligands was first calculated at the MM level by an alchemical transformation from one ligand to another in both the bound and unbound states. Then, for each ligand the free-energy correction for going from the MM to the SQM/MM potentials was calculated using explicit SQM/MM MD simulations. Finally, the free-energy correction for going from the SQM/MM to the DFT/MM potentials was estimated with FEP without running any DFT/MM simulations. Instead, the free energy was calculated by single-step exponential averaging (ssEA) or employing the cumulant approximation to the second order (CA). The results show that CA converges much better than ssEA, and with 500-4500 DFT/MM single-point energy calculations, converged free energies with a precision of 0.3 kJ/mol can be obtained. These free energies reproduce the experimental binding free energy differences with a mean absolute deviation of 3.4 kJ/mol, a correlation ( R2) of 0.97, and correct signs for all of the eight free-energy differences. This is appreciably better than the results obtained at the SQM/MM level of theory and also slightly better than those obtained with MM. We show that the convergence of the SQM/MM → DFT/MM perturbations can be monitored by the use of Wu and Kofke's bias metric Π and by the standard deviation of the difference between the SQM/MM and DFT/MM energies. Finally, we show that the use of the intermediate SQM/MM MD simulations improves the convergence of the free energies by a factor of at least two, compared to doing direct MM → DFT/MM perturbations.
Collapse
Affiliation(s)
- Meiting Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science , East China Normal University , Shanghai 200062 , China.,Department of Theoretical Chemistry , Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science , East China Normal University , Shanghai 200062 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062 , China.,Collaborative Innovation Center of Extreme Optics , Shanxi University , Taiyuan , Shanxi 030006 , China
| | - Ulf Ryde
- Department of Theoretical Chemistry , Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| |
Collapse
|
23
|
Wang M, Mei Y, Ryde U. Predicting Relative Binding Affinity Using Nonequilibrium QM/MM Simulations. J Chem Theory Comput 2018; 14:6613-6622. [DOI: 10.1021/acs.jctc.8b00685] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Meiting Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU−ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
24
|
König G, Pickard FC, Huang J, Thiel W, MacKerell AD, Brooks BR, York DM. A Comparison of QM/MM Simulations with and without the Drude Oscillator Model Based on Hydration Free Energies of Simple Solutes. Molecules 2018; 23:E2695. [PMID: 30347691 PMCID: PMC6222909 DOI: 10.3390/molecules23102695] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2018] [Revised: 10/15/2018] [Accepted: 10/16/2018] [Indexed: 12/01/2022] Open
Abstract
Maintaining a proper balance between specific intermolecular interactions and non-specific solvent interactions is of critical importance in molecular simulations, especially when predicting binding affinities or reaction rates in the condensed phase. The most rigorous metric for characterizing solvent affinity are solvation free energies, which correspond to a transfer from the gas phase into solution. Due to the drastic change of the electrostatic environment during this process, it is also a stringent test of polarization response in the model. Here, we employ both the CHARMM fixed charge and polarizable force fields to predict hydration free energies of twelve simple solutes. The resulting classical ensembles are then reweighted to obtain QM/MM hydration free energies using a variety of QM methods, including MP2, Hartree⁻Fock, density functional methods (BLYP, B3LYP, M06-2X) and semi-empirical methods (OM2 and AM1 ). Our simulations test the compatibility of quantum-mechanical methods with molecular-mechanical water models and solute Lennard⁻Jones parameters. In all cases, the resulting QM/MM hydration free energies were inferior to purely classical results, with the QM/MM Drude force field predictions being only marginally better than the QM/MM fixed charge results. In addition, the QM/MM results for different quantum methods are highly divergent, with almost inverted trends for polarizable and fixed charge water models. While this does not necessarily imply deficiencies in the QM models themselves, it underscores the need to develop consistent and balanced QM/MM interactions. Both the QM and the MM component of a QM/MM simulation have to match, in order to avoid artifacts due to biased solute⁻solvent interactions. Finally, we discuss strategies to improve the convergence and efficiency of multi-scale free energy simulations by automatically adapting the molecular-mechanics force field to the target quantum method.
Collapse
Affiliation(s)
- Gerhard König
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA.
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Jing Huang
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
- Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA.
- School of Life Sciences, Westlake University, 18 Shilongshan Street, Hangzhou 310024, China.
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
| | - Alexander D MacKerell
- Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA.
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA.
| |
Collapse
|
25
|
Hudson PS, Han K, Woodcock HL, Brooks BR. Force matching as a stepping stone to QM/MM CB[8] host/guest binding free energies: a SAMPL6 cautionary tale. J Comput Aided Mol Des 2018; 32:983-999. [PMID: 30276502 PMCID: PMC6867086 DOI: 10.1007/s10822-018-0165-3] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2018] [Accepted: 09/14/2018] [Indexed: 10/28/2022]
Abstract
Use of quantum mechanical/molecular mechanical (QM/MM) methods in binding free energy calculations, particularly in the SAMPL challenge, often fail to achieve improvement over standard additive (MM) force fields. Frequently, the implementation is through use of reference potentials, or the so-called "indirect approach", and inherently relies on sufficient overlap existing between MM and QM/MM configurational spaces. This overlap is generally poor, particularly for the use of free energy perturbation to perform the MM to QM/MM free energy correction at the end states of interest (e.g., bound and unbound states). However, by utilizing MM parameters that best reproduce forces obtained at the desired QM level of theory, it is possible to lessen the configurational disparity between MM and QM/MM. To this end, we sought to use force matching to generate MM parameters for the SAMPL6 CB[8] host-guest binding challenge, classically compute binding free energies, and apply energetic end state corrections to obtain QM/MM binding free energy differences. For the standard set of 11 molecules and the bonus set (including three additional challenge molecules), error statistics, such as the root mean square deviation (RMSE) were moderately poor (5.5 and 5.4 kcal/mol). Correlation statistics, however, were in the top two for both standard and bonus set submissions ([Formula: see text] of 0.42 and 0.26, [Formula: see text] of 0.64 and 0.47 respectively). High RMSE and moderate correlation strongly indicated the presence of systematic error. Identifiable issues were ameliorated for two of the guest molecules, resulting in a reduction of error and pointing to strong prospects for the future use of this methodology.
Collapse
Affiliation(s)
- Phillip S Hudson
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20852, USA.
- Department of Chemistry, University of South Florida, Tampa, Florida, 33620, USA.
| | - Kyungreem Han
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20852, USA
| | - H Lee Woodcock
- Department of Chemistry, University of South Florida, Tampa, Florida, 33620, USA
| | - Bernard R Brooks
- Department of Chemistry, University of South Florida, Tampa, Florida, 33620, USA
| |
Collapse
|
26
|
Caldararu O, Olsson MA, Misini Ignjatović M, Wang M, Ryde U. Binding free energies in the SAMPL6 octa-acid host-guest challenge calculated with MM and QM methods. J Comput Aided Mol Des 2018; 32:1027-1046. [PMID: 30203229 DOI: 10.1007/s10822-018-0158-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Accepted: 08/31/2018] [Indexed: 01/15/2023]
Abstract
We have estimated free energies for the binding of eight carboxylate ligands to two variants of the octa-acid deep-cavity host in the SAMPL6 blind-test challenge (with or without endo methyl groups on the four upper-rim benzoate groups, OAM and OAH, respectively). We employed free-energy perturbation (FEP) for relative binding energies at the molecular mechanics (MM) and the combined quantum mechanical (QM) and MM (QM/MM) levels, the latter obtained with the reference-potential approach with QM/MM sampling for the MM → QM/MM FEP. The semiempirical QM method PM6-DH+ was employed for the ligand in the latter calculations. Moreover, binding free energies were also estimated from QM/MM optimised structures, combined with COSMO-RS estimates of the solvation energy and thermostatistical corrections from MM frequencies. They were performed at the PM6-DH+ level of theory with the full host and guest molecule in the QM system (and also four water molecules in the geometry optimisations) for 10-20 snapshots from molecular dynamics simulations of the complex. Finally, the structure with the lowest free energy was recalculated using the dispersion-corrected density-functional theory method TPSS-D3, for both the structure and the energy. The two FEP approaches gave similar results (PM6-DH+/MM slightly better for OAM), which were among the five submissions with the best performance in the challenge and gave the best results without any fit to data from the SAMPL5 challenge, with mean absolute deviations (MAD) of 2.4-5.2 kJ/mol and a correlation coefficient (R2) of 0.77-0.93. This is the first time QM/MM approaches give binding free energies that are competitive to those obtained with MM for the octa-acid host. The QM/MM-optimised structures gave somewhat worse performance (MAD = 3-8 kJ/mol and R2 = 0.1-0.9), but the results were improved compared to previous studies of this system with similar methods.
Collapse
Affiliation(s)
- Octav Caldararu
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Martin A Olsson
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Majda Misini Ignjatović
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Meiting Wang
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden.
| |
Collapse
|
27
|
A new approach for the acceleration of large-scale serial quantum chemical calculations of docking complexes. Russ Chem Bull 2018. [DOI: 10.1007/s11172-018-2186-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
28
|
König G, Brooks BR, Thiel W, York DM. On the convergence of multi-scale free energy simulations. MOLECULAR SIMULATION 2018; 44:1062-1081. [PMID: 30581251 DOI: 10.1080/08927022.2018.1475741] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Abstract
In this work we employ simple model systems to evaluate the relative performance of two of the most important free energy methods: The Zwanzig equation (also known as "Free energy perturbation") and Bennett's acceptance ratio method (BAR). Although our examples should be transferable to other kinds of free energy simulations, we focus on applications of multi-scale free energy simulations. Such calculations are especially complex, since they connect two different levels of theory with very different requirements in terms of speed, accuracy, sampling and parallelizability. We try to reconcile all those different factors by developing some simple criteria to guide the early stages of the development of a free energy protocol. This is accomplished by quantifying how many λ intermediate steps and how many potential energy evaluations are necessary in order to reach a certain level of convergence.
Collapse
Affiliation(s)
- Gerhard König
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany, EU.,Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA.,Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany, EU
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| |
Collapse
|
29
|
Steinmann C, Olsson MA, Ryde U. Relative Ligand-Binding Free Energies Calculated from Multiple Short QM/MM MD Simulations. J Chem Theory Comput 2018; 14:3228-3237. [PMID: 29768915 DOI: 10.1021/acs.jctc.8b00081] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
We have devised a new efficient approach to compute combined quantum mechanical (QM) and molecular mechanical (MM, i.e. QM/MM) ligand-binding relative free energies. Our method employs the reference-potential approach with free-energy perturbation both at the MM level (between the two ligands) and from MM to QM/MM (for each ligand). To ensure that converged results are obtained for the MM → QM/MM perturbations, explicit QM/MM molecular dynamics (MD) simulations are performed with two intermediate mixed states. To speed up the calculations, we utilize the fact that the phase space can be extensively sampled at the MM level. Therefore, we run many short QM/MM MD simulations started from snapshots of the MM simulations, instead of a single long simulation. As a test case, we study the binding of nine cyclic carboxylate ligands to the octa-acid deep cavitand. Only the ligand is in the QM system, treated with the semiempirical PM6-DH+ method. We show that for eight of the ligands, we obtain well converged results with short MD simulations (1-15 ps). However, in one case, the convergence is slower (∼50 ps) owing to a mismatch between the conformational preferences of the MM and QM/MM potentials. We test the effect of initial minimization, the need of equilibration, and how many independent simulations are needed to reach a certain precision. The results show that the present approach is about four times faster than using standard MM → QM/MM free-energy perturbations with the same accuracy and precision.
Collapse
Affiliation(s)
- Casper Steinmann
- Department of Chemistry and Bioscience , Aalborg University , Frederik Bajers Vej 7H , DK-9220 Aalborg , Denmark.,Department of Theoretical Chemistry , Lund University , Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| | - Martin A Olsson
- Department of Theoretical Chemistry , Lund University , Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry , Lund University , Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| |
Collapse
|
30
|
Manzoni F, Ryde U. Assessing the stability of free-energy perturbation calculations by performing variations in the method. J Comput Aided Mol Des 2018. [PMID: 29536221 PMCID: PMC5889414 DOI: 10.1007/s10822-018-0110-5] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
We have calculated relative binding affinities for eight tetrafluorophenyl-triazole-thiogalactoside inhibitors of galectin-3 with the alchemical free-energy perturbation approach. We obtain a mean absolute deviation from experimental estimates of only 2-3 kJ/mol and a correlation coefficient (R2) of 0.5-0.8 for seven relative affinities spanning a range of up to 11 kJ/mol. We also studied the effect of using different methods to calculate the charges of the inhibitor and different sizes of the perturbed group (the atoms that are described by soft-core potentials and are allowed to have differing coordinates). However, the various approaches gave rather similar results and it is not possible to point out one approach as consistently and significantly better than the others. Instead, we suggest that such small and reasonable variations in the computational method can be used to check how stable the calculated results are and to obtain a more accurate estimate of the uncertainty than if performing only one calculation with a single computational setup.
Collapse
Affiliation(s)
- Francesco Manzoni
- Theoretical Chemistry, Department of Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden
| | - Ulf Ryde
- Theoretical Chemistry, Department of Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden.
| |
Collapse
|
31
|
Ryde U. How Many Conformations Need To Be Sampled To Obtain Converged QM/MM Energies? The Curse of Exponential Averaging. J Chem Theory Comput 2017; 13:5745-5752. [DOI: 10.1021/acs.jctc.7b00826] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Ulf Ryde
- Department of Theoretical
Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
32
|
Wang M, Li P, Jia X, Liu W, Shao Y, Hu W, Zheng J, Brooks BR, Mei Y. Efficient Strategy for the Calculation of Solvation Free Energies in Water and Chloroform at the Quantum Mechanical/Molecular Mechanical Level. J Chem Inf Model 2017; 57:2476-2489. [DOI: 10.1021/acs.jcim.7b00001] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Meiting Wang
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Pengfei Li
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Xiangyu Jia
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Wei Liu
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Yihan Shao
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman Oklahoma 73019, United States
| | - Wenxin Hu
- The Computer
Center, School of Computer Science and Software Engineering, East China Normal University, Shanghai 200062, China
| | - Jun Zheng
- The Computer
Center, School of Computer Science and Software Engineering, East China Normal University, Shanghai 200062, China
| | - Bernard R. Brooks
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Ye Mei
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
33
|
Olsson MA, Ryde U. Comparison of QM/MM Methods To Obtain Ligand-Binding Free Energies. J Chem Theory Comput 2017; 13:2245-2253. [DOI: 10.1021/acs.jctc.6b01217] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Martin A. Olsson
- Department of Theoretical
Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221
00 Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical
Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221
00 Lund, Sweden
| |
Collapse
|
34
|
Cave-Ayland C, Skylaris CK, Essex JW. A Monte Carlo Resampling Approach for the Calculation of Hybrid Classical and Quantum Free Energies. J Chem Theory Comput 2017; 13:415-424. [DOI: 10.1021/acs.jctc.6b00506] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
| | | | - Jonathan W. Essex
- School of Chemistry, University of Southampton, Hampshire, SO17 1BJ, United Kingdom
| |
Collapse
|
35
|
Caldararu O, Olsson MA, Riplinger C, Neese F, Ryde U. Binding free energies in the SAMPL5 octa-acid host-guest challenge calculated with DFT-D3 and CCSD(T). J Comput Aided Mol Des 2017; 31:87-106. [PMID: 27600554 PMCID: PMC5239813 DOI: 10.1007/s10822-016-9957-5] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2016] [Accepted: 08/25/2016] [Indexed: 10/27/2022]
Abstract
We have tried to calculate the free energy for the binding of six small ligands to two variants of the octa-acid deep cavitand host in the SAMPL5 blind challenge. We employed structures minimised with dispersion-corrected density-functional theory with small basis sets and energies were calculated using large basis sets. Solvation energies were calculated with continuum methods and thermostatistical corrections were obtained from frequencies calculated at the HF-3c level. Care was taken to minimise the effects of the flexibility of the host by keeping the complexes as symmetric and similar as possible. In some calculations, the large net charge of the host was reduced by removing the propionate and benzoate groups. In addition, the effect of a restricted molecular dynamics sampling of structures was tested. Finally, we tried to improve the energies by using the DLPNO-CCSD(T) approach. Unfortunately, results of quite poor quality were obtained, with no correlation to the experimental data, systematically too positive affinities (by ~50 kJ/mol) and a mean absolute error (after removal of the systematic error) of 11-16 kJ/mol. DLPNO-CCSD(T) did not improve the results, so the accuracy is not limited by the energy function. Instead, four likely sources of errors were identified: first, the minimised structures were often incorrect, owing to the omission of explicit solvent. They could be partly improved by performing the minimisations in a continuum solvent with four water molecules around the charged groups of the ligands. Second, some ligands could bind in several different conformations, requiring sampling of reasonable structures. Third, there is an indication the continuum-solvation model has problems to accurately describe the binding of both the negatively and positively charged guest molecules. Fourth, different methods to calculate the thermostatistical corrections gave results that differed by up to 30 kJ/mol and there is an indication that HF-3c overestimates the entropy term. In conclusion, it is a challenge to calculate binding affinities for this octa-acid system with quantum-mechanical methods.
Collapse
Affiliation(s)
- Octav Caldararu
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden
| | - Martin A Olsson
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden
| | - Christoph Riplinger
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstraße 34-36, 45470, Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstraße 34-36, 45470, Mülheim an der Ruhr, Germany
| | - Ulf Ryde
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden.
| |
Collapse
|
36
|
Gresh N, Perahia D, de Courcy B, Foret J, Roux C, El-Khoury L, Piquemal JP, Salmon L. Complexes of a Zn-metalloenzyme binding site with hydroxamate-containing ligands. A case for detailed benchmarkings of polarizable molecular mechanics/dynamics potentials when the experimental binding structure is unknown. J Comput Chem 2016; 37:2770-2782. [DOI: 10.1002/jcc.24503] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2016] [Revised: 08/31/2016] [Accepted: 09/04/2016] [Indexed: 12/15/2022]
Affiliation(s)
- Nohad Gresh
- Laboratoire de Chimie Théorique; Sorbonne Universités; UPMC, UMR 7616 CNRS Paris France
- Chemistry and Biology, Nucleo(s)tides and Immunology for Therapy (CBNIT); UMR 8601 CNRS, UFR Biomédicale; Paris France
| | - David Perahia
- Laboratoire de Biologie et Pharmacologie Appliquées (LBPA), UMR 8113; Ecole Normale Supérieure Cachan France
| | - Benoit de Courcy
- Laboratoire de Chimie Théorique; Sorbonne Universités; UPMC, UMR 7616 CNRS Paris France
- Chemistry and Biology, Nucleo(s)tides and Immunology for Therapy (CBNIT); UMR 8601 CNRS, UFR Biomédicale; Paris France
| | - Johanna Foret
- Laboratoire de Chimie Bioorganique et Bioinorganique; Institut de Chimie Moléculaire et des Matériaux d'Orsay (ICMMO), Univ Paris-Saclay, Univ Paris-Sud, UMR 8182 CNRS; rue du Doyen Georges Poitou Orsay F-91405 France
| | - Céline Roux
- Laboratoire de Chimie Bioorganique et Bioinorganique; Institut de Chimie Moléculaire et des Matériaux d'Orsay (ICMMO), Univ Paris-Saclay, Univ Paris-Sud, UMR 8182 CNRS; rue du Doyen Georges Poitou Orsay F-91405 France
| | - Lea El-Khoury
- Laboratoire de Chimie Théorique; Sorbonne Universités; UPMC, UMR 7616 CNRS Paris France
- Centre d'Analyses et de Recherche; UR EGFEM, LSIM, Faculté de Sciences, Saint Joseph University of Beirut; BP 11-514, Riad El Solh Beirut 1116-2050 Lebanon
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique; Sorbonne Universités; UPMC, UMR 7616 CNRS Paris France
- Department of Biomolecular Engineering; The University of Texas at Austin; Texas 78712
| | - Laurent Salmon
- Laboratoire de Chimie Bioorganique et Bioinorganique; Institut de Chimie Moléculaire et des Matériaux d'Orsay (ICMMO), Univ Paris-Saclay, Univ Paris-Sud, UMR 8182 CNRS; rue du Doyen Georges Poitou Orsay F-91405 France
| |
Collapse
|
37
|
Pickard FC, König G, Tofoleanu F, Lee J, Simmonett AC, Shao Y, Ponder JW, Brooks BR. Blind prediction of distribution in the SAMPL5 challenge with QM based protomer and pK a corrections. J Comput Aided Mol Des 2016; 30:1087-1100. [PMID: 27646286 DOI: 10.1007/s10822-016-9955-7] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2016] [Accepted: 08/25/2016] [Indexed: 12/14/2022]
Abstract
The computation of distribution coefficients between polar and apolar phases requires both an accurate characterization of transfer free energies between phases and proper accounting of ionization and protomerization. We present a protocol for accurately predicting partition coefficients between two immiscible phases, and then apply it to 53 drug-like molecules in the SAMPL5 blind prediction challenge. Our results combine implicit solvent QM calculations with classical MD simulations using the non-Boltzmann Bennett free energy estimator. The OLYP/DZP/SMD method yields predictions that have a small deviation from experiment (RMSD = 2.3 [Formula: see text] D units), relative to other participants in the challenge. Our free energy corrections based on QM protomer and [Formula: see text] calculations increase the correlation between predicted and experimental distribution coefficients, for all methods used. Unfortunately, these corrections are overly hydrophilic, and fail to account for additional effects such as aggregation, water dragging and the presence of polar impurities in the apolar phase. We show that, although expensive, QM-NBB free energy calculations offer an accurate and robust method that is superior to standard MM and QM techniques alone.
Collapse
Affiliation(s)
- Frank C Pickard
- Laboratory of Computational Biology, National Institutes of Health - National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, MD, 20852, USA.
| | - Gerhard König
- Laboratory of Computational Biology, National Institutes of Health - National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, MD, 20852, USA.,Max Planck Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, NRW, Germany
| | - Florentina Tofoleanu
- Laboratory of Computational Biology, National Institutes of Health - National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, MD, 20852, USA
| | - Juyong Lee
- Laboratory of Computational Biology, National Institutes of Health - National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, MD, 20852, USA
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Institutes of Health - National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, MD, 20852, USA
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, OK, 73019, USA
| | - Jay W Ponder
- Department of Chemistry, Washington University, St. Louis, MO, 63130, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Institutes of Health - National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, MD, 20852, USA
| |
Collapse
|
38
|
König G, Pickard FC, Huang J, Simmonett AC, Tofoleanu F, Lee J, Dral PO, Prasad S, Jones M, Shao Y, Thiel W, Brooks BR. Calculating distribution coefficients based on multi-scale free energy simulations: an evaluation of MM and QM/MM explicit solvent simulations of water-cyclohexane transfer in the SAMPL5 challenge. J Comput Aided Mol Des 2016; 30:989-1006. [PMID: 27577746 DOI: 10.1007/s10822-016-9936-x] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2016] [Accepted: 08/09/2016] [Indexed: 11/29/2022]
Abstract
One of the central aspects of biomolecular recognition is the hydrophobic effect, which is experimentally evaluated by measuring the distribution coefficients of compounds between polar and apolar phases. We use our predictions of the distribution coefficients between water and cyclohexane from the SAMPL5 challenge to estimate the hydrophobicity of different explicit solvent simulation techniques. Based on molecular dynamics trajectories with the CHARMM General Force Field, we compare pure molecular mechanics (MM) with quantum-mechanical (QM) calculations based on QM/MM schemes that treat the solvent at the MM level. We perform QM/MM with both density functional theory (BLYP) and semi-empirical methods (OM1, OM2, OM3, PM3). The calculations also serve to test the sensitivity of partition coefficients to solute polarizability as well as the interplay of the quantum-mechanical region with the fixed-charge molecular mechanics environment. Our results indicate that QM/MM with both BLYP and OM2 outperforms pure MM. However, this observation is limited to a subset of cases where convergence of the free energy can be achieved.
Collapse
Affiliation(s)
- Gerhard König
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA. .,Max-Planck-Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, Germany.
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Jing Huang
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Florentina Tofoleanu
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Juyong Lee
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Pavlo O Dral
- Max-Planck-Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, Germany
| | - Samarjeet Prasad
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Michael Jones
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Yihan Shao
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, Germany
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| |
Collapse
|
39
|
Pickard FC, König G, Simmonett AC, Shao Y, Brooks BR. An efficient protocol for obtaining accurate hydration free energies using quantum chemistry and reweighting from molecular dynamics simulations. Bioorg Med Chem 2016; 24:4988-4997. [PMID: 27667551 DOI: 10.1016/j.bmc.2016.08.031] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2016] [Revised: 08/18/2016] [Accepted: 08/19/2016] [Indexed: 11/15/2022]
Abstract
The non-Boltzmann Bennett (NBB) free energy estimator method is applied to 21 molecules from the blind subset of the SAMPL4 challenge. When NBB is applied with the SMD implicit solvent model, and the OLYP/DZP level of quantum chemistry, highly accurate hydration free energy calculations are obtained with respect to experiment (RMSD=0.89kcal·mol-1). Other quantum chemical methods are also tested, and the effects of solvent model, density functional, basis set are explored in this benchmarking study, providing a framework for improvements in calculating hydration free energies. We provide a practical guide for using the best QM-NBB protocols that are consistently more accurate than either pure QM or pure MM alone. In situations where high accuracy hydration free energy predictions are needed, the QM-NBB method with SMD implicit solvent should be the first choice of quantum chemists.
Collapse
Affiliation(s)
- Frank C Pickard
- National Institutes of Health - National Heart, Lung and Blood Institute, Laboratory of Computational Biology, 5635 Fishers Lane, T-900 Suite, Rockville, MD 20852, USA
| | - Gerhard König
- National Institutes of Health - National Heart, Lung and Blood Institute, Laboratory of Computational Biology, 5635 Fishers Lane, T-900 Suite, Rockville, MD 20852, USA; Max Planck Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, NRW, Germany
| | - Andrew C Simmonett
- National Institutes of Health - National Heart, Lung and Blood Institute, Laboratory of Computational Biology, 5635 Fishers Lane, T-900 Suite, Rockville, MD 20852, USA
| | - Yihan Shao
- Department of Chemistry and Biochemistry University of Oklahoma Norman, OK 73019, USA
| | - Bernard R Brooks
- National Institutes of Health - National Heart, Lung and Blood Institute, Laboratory of Computational Biology, 5635 Fishers Lane, T-900 Suite, Rockville, MD 20852, USA
| |
Collapse
|
40
|
Olsson MA, Söderhjelm P, Ryde U. Converging ligand-binding free energies obtained with free-energy perturbations at the quantum mechanical level. J Comput Chem 2016; 37:1589-600. [PMID: 27117350 PMCID: PMC5074236 DOI: 10.1002/jcc.24375] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Revised: 03/08/2016] [Accepted: 03/09/2016] [Indexed: 12/15/2022]
Abstract
In this article, the convergence of quantum mechanical (QM) free-energy simulations based on molecular dynamics simulations at the molecular mechanics (MM) level has been investigated. We have estimated relative free energies for the binding of nine cyclic carboxylate ligands to the octa-acid deep-cavity host, including the host, the ligand, and all water molecules within 4.5 Å of the ligand in the QM calculations (158-224 atoms). We use single-step exponential averaging (ssEA) and the non-Boltzmann Bennett acceptance ratio (NBB) methods to estimate QM/MM free energy with the semi-empirical PM6-DH2X method, both based on interaction energies. We show that ssEA with cumulant expansion gives a better convergence and uses half as many QM calculations as NBB, although the two methods give consistent results. With 720,000 QM calculations per transformation, QM/MM free-energy estimates with a precision of 1 kJ/mol can be obtained for all eight relative energies with ssEA, showing that this approach can be used to calculate converged QM/MM binding free energies for realistic systems and large QM partitions. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Martin A Olsson
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| | - Pär Söderhjelm
- Department of Biophysical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| |
Collapse
|
41
|
Ryde U, Söderhjelm P. Ligand-Binding Affinity Estimates Supported by Quantum-Mechanical Methods. Chem Rev 2016; 116:5520-66. [DOI: 10.1021/acs.chemrev.5b00630] [Citation(s) in RCA: 175] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ulf Ryde
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Pär Söderhjelm
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
42
|
|