1
|
Pederson JP, McDaniel JG. PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics. J Chem Phys 2024; 161:034103. [PMID: 39007371 DOI: 10.1063/5.0219851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2024] [Accepted: 06/24/2024] [Indexed: 07/16/2024] Open
Abstract
PyDFT-QMMM is a Python-based package for performing hybrid quantum mechanics/molecular mechanics (QM/MM) simulations at the density functional level of theory. The program is designed to treat short-range and long-range interactions through user-specified combinations of electrostatic and mechanical embedding procedures within periodic simulation domains, providing necessary interfaces to external quantum chemistry and molecular dynamics software. To enable direct embedding of long-range electrostatics in periodic systems, we have derived and implemented force terms for our previously described QM/MM/PME approach [Pederson and McDaniel, J. Chem. Phys. 156, 174105 (2022)]. Communication with external software packages Psi4 and OpenMM is facilitated through Python application programming interfaces (APIs). The core library contains basic utilities for running QM/MM molecular dynamics simulations, and plug-in entry-points are provided for users to implement custom energy/force calculation and integration routines, within an extensible architecture. The user interacts with PyDFT-QMMM primarily through its Python API, allowing for complex workflow development with Python scripting, for example, interfacing with PLUMED for free energy simulations. We provide benchmarks of forces and energy conservation for the QM/MM/PME and alternative QM/MM electrostatic embedding approaches. We further demonstrate a simple example use case for water solute in a water solvent system, for which radial distribution functions are computed from 100 ps QM/MM simulations; in this example, we highlight how the solvation structure is sensitive to different basis-set choices due to under- or over-polarization of the QM water molecule's electron density.
Collapse
Affiliation(s)
- John P Pederson
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Jesse G McDaniel
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
2
|
Zubatyuk R, Biczysko M, Ranasinghe K, Moriarty NW, Gokcan H, Kruse H, Poon BK, Adams PD, Waller MP, Roitberg AE, Isayev O, Afonine PV. AQuaRef: Machine learning accelerated quantum refinement of protein structures. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.07.21.604493. [PMID: 39071315 PMCID: PMC11275739 DOI: 10.1101/2024.07.21.604493] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
Cryo-EM and X-ray crystallography provide crucial experimental data for obtaining atomic-detail models of biomacromolecules. Refining these models relies on library- based stereochemical restraints, which, in addition to being limited to known chemical entities, do not include meaningful noncovalent interactions relying solely on nonbonded repulsions. Quantum mechanical (QM) calculations could alleviate these issues but are too expensive for large molecules. We present a novel AI-enabled Quantum Refinement (AQuaRef) based on AIMNet2 neural network potential mimicking QM at substantially lower computational costs. By refining 41 cryo-EM and 30 X-ray structures, we show that this approach yields atomic models with superior geometric quality compared to standard techniques, while maintaining an equal or better fit to experimental data.
Collapse
|
3
|
Yan S, Wang B, Lin H. Reshaping the QM Region On-the-Fly: Adaptive-Shape QM/MM Dynamic Simulations of a Hydrated Proton in Bulk Water. J Chem Theory Comput 2024; 20:3462-3472. [PMID: 38671391 DOI: 10.1021/acs.jctc.4c00164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/28/2024]
Abstract
Adaptive quantum mechanics/molecular mechanics (QM/MM) reclassifies on-the-fly a molecule or molecular fragment as QM or MM during dynamics simulations without abrupt changes in the energy or forces. Notably, the permuted adaptive-partitioning (PAP) algorithms have been applied to simulate a hydrated proton, with a mobile QM zone anchored at a pseudoatom called a proton indicator. The position of the proton indicator approximates the location of the delocalized excess proton, yielding a smooth trajectory of the proton diffusing via the Grotthuss mechanism in aqueous solutions. The mobile QM zone, which has been taken to be a sphere with a preset radius, follows the proton wherever it goes. Although the simulations are successful, the use of a spherical QM zone has one disadvantage: A large preset radius must be utilized to minimize the chance of missing water molecules that are important to proton translocation. A large radius leads to a large QM zone, which is computationally expensive. In this work, we report a new way to set up the QM zone, where one includes only the water molecules important to proton transfer. The importance of a given water molecule is quantified by its "weight" that depends on its relation to the reaction path of proton transfer. The weight varies smoothly, ensuring that a water molecule gradually appears in or disappears from the QM zone without abrupt changes, as required by the PAP method. Consequently, the shape of the QM zone evolves on-the-fly, keeping the QM zone as small as possible and as large as necessary. Test simulations demonstrate that the new algorithm significantly improves the computation efficiency while maintaining the proper descriptions of proton transfer in bulk water.
Collapse
Affiliation(s)
- Shengheng Yan
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, P. R. China
| | - Binju Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, P. R. China
| | - Hai Lin
- Department of Chemistry, CB 194, University of Colorado Denver, Denver, P.O. Box 173364, Colorado 80217, United States
| |
Collapse
|
4
|
Tran AL, Guidez EB, Lin H. Adaptive-Partitioning Multilayer Dynamics Simulations: 2. Implementations of the Permuted and Interpolated Adaptive-Partitioning Gradients. J Phys Chem A 2023; 127:10320-10333. [PMID: 38058156 PMCID: PMC10712430 DOI: 10.1021/acs.jpca.3c05600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2023] [Revised: 10/06/2023] [Accepted: 10/17/2023] [Indexed: 12/08/2023]
Abstract
Recently, an adaptive-partitioning multilayer Q1/Q2/MM method was proposed, where Q1 and Q2 denote, respectively, two distinct quantum-mechanical levels of theory and MM, the molecular-mechanical force fields. Such a multilayer model resembles the ONIOM (our own N-layered integrated molecular orbital and molecular mechanics) model by Morokuma and co-workers, but it is distinguished by on-the-fly reclassifying atoms to be Q1, Q2, or MM in dynamics simulations. To smoothly blend the levels of descriptions of the atoms, buffer zones are introduced between adjacent layers, and the energy is smoothly interpolated. In particular, the Q1/Q2 interaction energy was expressed in two different formalisms: permuted and interpolated adaptive-partitioning (PAP and IAP), respectively. While the PAP energy is based on a weighted many-body expansion, the IAP energy is derived via alchemical quantum calculations with interpolated Fock and overlap matrices. In this article, we examine in-depth the irregularities in the IAP energy near the boundary between the buffer and Q2 zones, which were found prominent in some calculations. These irregularities are due to basis-set linear dependencies, which can be effectively suppressed using a cutoff for the weighted atomic orbital coefficients. Furthermore, we derived and implemented the gradients for both PAP and IAP. Test calculations on a series of water cluster models show perfectly smooth gradients in PAP, while a minor discontinuity occurs in IAP gradients at the buffer/Q2 boundary. The energy and gradient discontinuities in IAP become smaller when moving the buffer/Q2 boundary further away from the Q1 center and when increasing the size of the basis sets used. Overall, those discontinuities are controllable, and possible ways to further diminish them are discussed.
Collapse
Affiliation(s)
- Anh L. Tran
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217, United States
| | - Emilie B. Guidez
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217, United States
| | - Hai Lin
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217, United States
| |
Collapse
|
5
|
Csizi K, Reiher M. Universal
QM
/
MM
approaches for general nanoscale applications. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
| | - Markus Reiher
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| |
Collapse
|
6
|
Yan S, Wang B, Lin H. Tracking the Delocalized Proton in Concerted Proton Transfer in Bulk Water. J Chem Theory Comput 2023; 19:448-459. [PMID: 36630655 DOI: 10.1021/acs.jctc.2c01097] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
A solvated proton in water is often characterized as a charge or structural defect, and it is important to track its evolution on-the-fly in certain dynamics simulations. Previously, we introduced the proton indicator, a pseudo-atom, whose position approximates the location of the excess proton modeled as a structural defect. The proton indicator generally yields a smooth trajectory of a hydrated proton diffusing in aqueous solutions, including in the events of stepwise proton transfer via the Grotthuss mechanism; however, the proton indicator did not perform well in the events of concerted proton transfer, for which it occasionally yielded large position displacements between two successive time steps. To overcome this hurdle, we develop a new algorithm of a proton indicator with greatly enhanced performance for concerted proton transfer in bulk water. A protocol is proposed to exhaustively explore the hydrogen-bonding network of the water wires over which the excess proton is delocalized and to properly account for the contributions of the water molecules in this network as the geometry evolves. The new proton indicator (called Indicator 2.0) is assessed in dynamics simulations of an excess proton in bulk water and in specially constructed model systems of more complex architectures. The results demonstrate that the new indicator yields a smooth trajectory in both stepwise and concerted proton transfers.
Collapse
Affiliation(s)
- Shengheng Yan
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen360015P. R. China
| | - Binju Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen360015P. R. China
| | - Hai Lin
- Department of Chemistry, CB 194, University of Colorado Denver, P.O. Box 173364, Denver, Colorado80217, United States
| |
Collapse
|
7
|
Nazemi A, Steeves AH, Kastner DW, Kulik HJ. Influence of the Greater Protein Environment on the Electrostatic Potential in Metalloenzyme Active Sites: The Case of Formate Dehydrogenase. J Phys Chem B 2022; 126:4069-4079. [PMID: 35609244 DOI: 10.1021/acs.jpcb.2c02260] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The Mo/W-containing metalloenzyme formate dehydrogenase (FDH) is an efficient and selective natural catalyst that reversibly converts CO2 to formate under ambient conditions. In this study, we investigate the impact of the greater protein environment on the electrostatic potential (ESP) of the active site. To model the enzyme environment, we used a combination of classical molecular dynamics and multiscale quantum-mechanical (QM)/molecular-mechanical (MM) simulations. We leverage charge shift analysis to systematically construct QM regions and analyze the electronic environment of the active site by evaluating the degree of charge transfer between the core active site and the protein environment. The contribution of the terminal chalcogen ligand to the ESP of the metal center is substantial and dependent on the chalcogen identity, with similar, less negative ESPs for Se and S terminal chalcogens in comparison to O regardless of whether the metal is Mo or W. The orientation of the side chains and conformations of the cofactor also affect the ESP, highlighting the importance of sampling dynamic fluctuations in the protein. Overall, our observations suggest that the terminal chalcogen ligand identity plays an important role in the enzymatic activity of FDH, suggesting opportunities for a rational bioinspired catalyst design.
Collapse
Affiliation(s)
- Azadeh Nazemi
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Adam H Steeves
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - David W Kastner
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
8
|
Saputri WD, Pranowo HD, Hofer TS. Can’t we negotiate the importance of electron correlation? HF vs RIMP2 in ab initio quantum mechanical charge field molecular dynamics simulations of Cu+ in pure liquid ammonia. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118286] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
9
|
Mato J, Duster AW, Guidez EB, Lin H. Adaptive-Partitioning Multilayer Dynamics Simulations: 1. On-the-Fly Switch between Two Quantum Levels of Theory. J Chem Theory Comput 2021; 17:5456-5465. [PMID: 34448578 DOI: 10.1021/acs.jctc.1c00556] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We propose to generalize the previously developed two-layer permuted adaptive-partitioning quantum-mechanics/molecular-mechanics (QM/MM), which reclassifies atoms as QM or MM on-the-fly in dynamics simulations, to multilayer adaptive-partitioning algorithms that enable multiple levels of theory. In this work, we formulate two new algorithms that smoothly interpolate the energy between two QM (Q1 and Q2) levels of theory. The first "permuted adaptive-partitioning" scheme is based on the weighted many-body expansion of the potential, as in the adaptive-partitioning QM/MM. Unconventional and potentially more efficient, the second "interpolated adaptive-partitioning" method employs alchemical QM calculations with Q1/Q2-mixed basis sets, Fock matrices, and overlap matrices. To our knowledge, this is the first time that such alchemical calculations are performed in QM, although they are routinely done in MM. Test calculations on water-cluster models show that both new algorithms indeed yield smooth energy curves when water molecules shift between Q1 and Q2.
Collapse
Affiliation(s)
- Joani Mato
- Department of Chemistry, University of Colorado, Denver, Denver, Colorado 80217, United States
| | - Adam W Duster
- Department of Chemistry, University of Colorado, Denver, Denver, Colorado 80217, United States
| | - Emilie B Guidez
- Department of Chemistry, University of Colorado, Denver, Denver, Colorado 80217, United States
| | - Hai Lin
- Department of Chemistry, University of Colorado, Denver, Denver, Colorado 80217, United States
| |
Collapse
|
10
|
Lambros E, Dasgupta S, Palos E, Swee S, Hu J, Paesani F. General Many-Body Framework for Data-Driven Potentials with Arbitrary Quantum Mechanical Accuracy: Water as a Case Study. J Chem Theory Comput 2021; 17:5635-5650. [PMID: 34370954 DOI: 10.1021/acs.jctc.1c00541] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a general framework for the development of data-driven many-body (MB) potential energy functions (MB-QM PEFs) that represent the interactions between small molecules at an arbitrary quantum-mechanical (QM) level of theory. As a demonstration, a family of MB-QM PEFs for water is rigorously derived from density functionals belonging to different rungs across Jacob's ladder of approximations within density functional theory (MB-DFT) and from Møller-Plesset perturbation theory (MB-MP2). Through a systematic analysis of individual MB contributions to the interaction energies of water clusters, we demonstrate that all MB-QM PEFs preserve the same accuracy as the corresponding ab initio calculations, with the exception of those derived from density functionals within the generalized gradient approximation (GGA). The differences between the DFT and MB-DFT results are traced back to density-driven errors that prevent GGA functionals from accurately representing the underlying molecular interactions for different cluster sizes and hydrogen-bonding arrangements. We show that this shortcoming may be overcome, within the MB formalism, by using density-corrected functionals (DC-DFT) that provide a more consistent representation of each individual MB contribution. This is demonstrated through the development of a MB-DFT PEF derived from DC-PBE-D3 data, which more accurately reproduce the corresponding ab initio results.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Etienne Palos
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Steven Swee
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Jie Hu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States.,Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States.,San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
11
|
Wappett DA, Goerigk L. A guide to benchmarking enzymatically catalysed reactions: the importance of accurate reference energies and the chemical environment. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02770-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|
12
|
Brunken C, Reiher M. Automated Construction of Quantum–Classical Hybrid Models. J Chem Theory Comput 2021; 17:3797-3813. [DOI: 10.1021/acs.jctc.1c00178] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Christoph Brunken
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
13
|
Schlick T, Portillo-Ledesma S. Biomolecular modeling thrives in the age of technology. NATURE COMPUTATIONAL SCIENCE 2021; 1:321-331. [PMID: 34423314 PMCID: PMC8378674 DOI: 10.1038/s43588-021-00060-9] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 03/22/2021] [Indexed: 12/12/2022]
Abstract
The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, New York, NY, USA
- Courant Institute of Mathematical Sciences, New York University, New York, NY, USA
- New York University–East China Normal University Center for Computational Chemistry at New York University Shanghai, Shanghai, China
| | | |
Collapse
|
14
|
Cruzeiro VWD, Manathunga M, Merz KM, Götz AW. Open-Source Multi-GPU-Accelerated QM/MM Simulations with AMBER and QUICK. J Chem Inf Model 2021; 61:2109-2115. [PMID: 33913331 DOI: 10.1021/acs.jcim.1c00169] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The quantum mechanics/molecular mechanics (QM/MM) approach is an essential and well-established tool in computational chemistry that has been widely applied in a myriad of biomolecular problems in the literature. In this publication, we report the integration of the QUantum Interaction Computational Kernel (QUICK) program as an engine to perform electronic structure calculations in QM/MM simulations with AMBER. This integration is available through either a file-based interface (FBI) or an application programming interface (API). Since QUICK is an open-source GPU-accelerated code with multi-GPU parallelization, users can take advantage of "free of charge" GPU-acceleration in their QM/MM simulations. In this work, we discuss implementation details and give usage examples. We also investigate energy conservation in typical QM/MM simulations performed at the microcanonical ensemble. Finally, benchmark results for two representative systems in bulk water, the N-methylacetamide (NMA) molecule and the photoactive yellow protein (PYP), show the performance of QM/MM simulations with QUICK and AMBER using a varying number of CPU cores and GPUs. Our results highlight the acceleration obtained from a single or multiple GPUs; we observed speedups of up to 53× between a single GPU vs a single CPU core and of up to 2.6× when comparing four GPUs to a single GPU. Results also reveal speedups of up to 3.5× when the API is used instead of FBI.
Collapse
Affiliation(s)
- Vinícius Wilian D Cruzeiro
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States.,Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Madushanka Manathunga
- Department of Chemistry, Department of Biochemistry and Molecular Biology, Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| | - Kenneth M Merz
- Department of Chemistry, Department of Biochemistry and Molecular Biology, Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
15
|
Wang JN, Liu W, Li P, Mo Y, Hu W, Zheng J, Pan X, Shao Y, Mei Y. Accelerated Computation of Free Energy Profile at Ab Initio Quantum Mechanical/Molecular Mechanics Accuracy via a Semiempirical Reference Potential. 4. Adaptive QM/MM. J Chem Theory Comput 2021; 17:1318-1325. [PMID: 33593057 PMCID: PMC8335528 DOI: 10.1021/acs.jctc.0c01149] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Although quantum mechanical/molecular mechanics (QM/MM) methods are now routinely applied to the studies of chemical reactions in condensed phases and enzymatic reactions, they may experience technical difficulties when the reactive region is varying over time. For instance, when the solvent molecules are directly participating in the reaction, the exchange of water molecules between the QM and MM regions may occur on a time scale comparable to the reaction time. To cope with this situation, several adaptive QM/MM schemes have been proposed. However, these methods either add significantly to the computational cost or introduce artificial restraints to the system. In this work, we developed a novel adaptive QM/MM scheme and applied it to the study of a nucleophilic addition reaction. In this scheme, the configuration sampling was performed with a small QM region (without solvent molecules), and the thermodynamic properties under another potential energy function with a larger QM region (with a certain number of solvent molecules and/or different levels of QM theory) are computed via extrapolation using the reference-potential method. Our simulation results show that this adaptive QM/MM scheme is numerically stable, at least for the case studied in this work. Furthermore, this method also offers an inexpensive way to examine the convergence of the QM/MM calculation with respect to the size of the QM region.
Collapse
Affiliation(s)
- Jia-Ning Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Wei Liu
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Pengfei Li
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Yan Mo
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic 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
| | - Wenxin Hu
- The Computer Center, School of Data Science & Engineering, East China Normal University, Shanghai 200062, China
| | - Jun Zheng
- The Computer Center, School of Data Science & Engineering, East China Normal University, Shanghai 200062, China
| | - Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic 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
| |
Collapse
|
16
|
Lambros E, Lipparini F, Cisneros GA, Paesani F. A Many-Body, Fully Polarizable Approach to QM/MM Simulations. J Chem Theory Comput 2020; 16:7462-7472. [PMID: 33213149 PMCID: PMC8131112 DOI: 10.1021/acs.jctc.0c00932] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
We present a new development in quantum mechanics/molecular mechanics (QM/MM) methods by replacing conventional MM models with data-driven many-body (MB) representations rigorously derived from high-level QM calculations. The new QM/MM approach builds on top of mutually polarizable QM/MM schemes developed for polarizable force fields with inducible dipoles and uses permutationally invariant polynomials to effectively account for quantum-mechanical contributions (e.g., exchange-repulsion and charge transfer and penetration) that are difficult to describe by classical expressions adopted by conventional MM models. Using the many-body MB-pol and MB-DFT potential energy functions for water, which include explicit two-body and three-body terms fitted to reproduce the corresponding CCSD(T) and PBE0 two-body and three-body energies for water, we demonstrate a smooth energetic transition as molecules are transferred between QM and MM regions, without the need of a transition layer. By effectively elevating the accuracy of both the MM region and the QM/MM interface to that of the QM region, the new QM/MB-MM approach achieves an accuracy comparable to that obtained with a fully QM treatment of the entire system.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | | | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
17
|
Xin X, Niu X, Liu W, Wang D. Hybrid Solvation Model with First Solvation Shell for Calculation of Solvation Free Energy. Chemphyschem 2020; 21:762-769. [PMID: 32154979 DOI: 10.1002/cphc.202000039] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2020] [Revised: 02/15/2020] [Indexed: 02/03/2023]
Abstract
We present a hybrid solvation model with first solvation shell to calculate solvation free energies. This hybrid model combines the quantum mechanics and molecular mechanics methods with the analytical expression based on the Born solvation model to calculate solvation free energies. Based on calculated free energies of solvation and reaction profiles in gas phase, we set up a unified scheme to predict reaction profiles in solution. The predicted solvation free energies and reaction barriers are compared with experimental results for twenty bimolecular nucleophilic substitution reactions. These comparisons show that our hybrid solvation model can predict reliable solvation free energies and reaction barriers for chemical reactions of small molecules in aqueous solution.
Collapse
Affiliation(s)
- Xin Xin
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| | - Xiao Niu
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| | - Wanqi Liu
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| | - Dunyou Wang
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| |
Collapse
|
18
|
Brunken C, Reiher M. Self-Parametrizing System-Focused Atomistic Models. J Chem Theory Comput 2020; 16:1646-1665. [DOI: 10.1021/acs.jctc.9b00855] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Christoph Brunken
- Laboratory for Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Laboratory for Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
19
|
Dohn AO, Jónsson EÖ, Jónsson H. Polarizable Embedding with a Transferable H 2O Potential Function II: Application to (H 2O) n Clusters and Liquid Water. J Chem Theory Comput 2019; 15:6578-6587. [PMID: 31692344 DOI: 10.1021/acs.jctc.9b00778] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The incorporation of polarization in multiscale quantum-mechanics/molecular-mechanics (QM/MM) simulations is important for a variety of applications, for example, charge-transfer reactions. A recently developed formalism based on a density functional theory description of the QM region and a potential energy function for H2O molecules that includes quadrupole as well as dipole polarizability of the MM region is used to simulate liquid water and water clusters. Analysis of the energy, atomic forces, MM polarization, and structure is presented. A quantitative assessment of the QM/MM-MM/MM interaction energy differences of all possible QM/MM configurations of (H2O)n clusters shows that the interquartile range of the distributions of the QM/MM binding energies is never more than 20 meV/molecule higher or lower than the binding energies produced with either of the single-model results. Comparing these interaction energy differences with the QM/MM induction differences show that they are not systematically caused by the induced MM moments of our polarizable embedding scheme. Optimized hexamer geometries as well as the liquid water structure are shown to be improved in comparison with results obtained using point-charge based embedding models neglecting polarization.
Collapse
Affiliation(s)
- Asmus Ougaard Dohn
- Science Institute and Faculty of Physical Sciences , University of Iceland , Reykjavík 107 , Iceland
| | - Elvar Örn Jónsson
- Science Institute and Faculty of Physical Sciences , University of Iceland , Reykjavík 107 , Iceland
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences , University of Iceland , Reykjavík 107 , Iceland
| |
Collapse
|
20
|
Jónsson EÖ, Dohn AO, Jónsson H. Polarizable Embedding with a Transferable H 2O Potential Function I: Formulation and Tests on Dimer. J Chem Theory Comput 2019; 15:6562-6577. [PMID: 31689104 DOI: 10.1021/acs.jctc.9b00777] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The incorporation of mutual polarization in multiscale simulations where different regions of the system are treated at different level of theory is important in studies of, for example, electronic excitations and charge transfer processes. We present here an energy functional for describing a quantum mechanics/molecular mechanics (QM/MM) scheme that includes reciprocal polarization between the two subsystems. The inclusion of polarization alleviates shortcomings inherent in electrostatic embedding QM/MM models based on point-charge force fields. A density functional theory (DFT) description of the QM subsystem is coupled to a single center multipole expansion (SCME) description of H2O molecules in the MM subsystem that includes anisotropic dipole and quadrupole polarizability as well as static multipoles up to and including the hexadecapole. The energy functional and the coupling scheme is general and can be extended to arbitrary order in terms of both the static and induced moments. Tests of the energy surface for the H2O dimer show that the QM/MM results lie in between the pure DFT and pure SCME values. The consistency of the many-body contributions to the energy and analytical forces is demonstrated for an H2O pentamer.
Collapse
Affiliation(s)
- Elvar Örn Jónsson
- Science Institute and Faculty of Physical Sciences, VR-III , University of Iceland , Reykjavík 107 , Iceland
| | - Asmus Ougaard Dohn
- Science Institute and Faculty of Physical Sciences, VR-III , University of Iceland , Reykjavík 107 , Iceland
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, VR-III , University of Iceland , Reykjavík 107 , Iceland
| |
Collapse
|
21
|
Advances and challenges in modeling solvated reaction mechanisms for renewable fuels and chemicals. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1446] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
22
|
Duster AW, Lin H. Tracking Proton Transfer through Titratable Amino Acid Side Chains in Adaptive QM/MM Simulations. J Chem Theory Comput 2019; 15:5794-5809. [DOI: 10.1021/acs.jctc.9b00649] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Adam W. Duster
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Hai Lin
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| |
Collapse
|
23
|
Zhang B, Altarawy D, Barnes T, Turney JM, Schaefer HF. Janus: An Extensible Open-Source Software Package for Adaptive QM/MM Methods. J Chem Theory Comput 2019; 15:4362-4373. [DOI: 10.1021/acs.jctc.9b00182] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Boyi Zhang
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Doaa Altarawy
- The Molecular Sciences Software Institute, Virginia Tech, Blacksburg, Virginia 24060, United States
- Department of Computer and Systems Engineering, Alexandria University, Alexandria 21544, Egypt
| | - Taylor Barnes
- The Molecular Sciences Software Institute, Virginia Tech, Blacksburg, Virginia 24060, United States
| | - Justin M. Turney
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Henry F. Schaefer
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
24
|
Gupta AK, Thapa B, Raghavachari K. Exploring Reaction Energy Profiles Using the Molecules-in-Molecules Fragmentation-Based Approach. J Chem Theory Comput 2019; 15:3991-4002. [PMID: 31181886 DOI: 10.1021/acs.jctc.9b00152] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Molecules-in-Molecules (MIM) fragmentation-based approach has been successfully used in previous studies to obtain the energies, optimized geometries, and spectroscopic properties of large molecular systems. The present work delineates a protocol to study the potential energy profiles for multistep chemical reactions using the MIM methodology. In a complex multistep chemical reaction, the fragmentation scheme needs to be changed as the reacting species transition into a new reaction step, resulting in a discontinuity in the potential energy curve of the reaction. In our approach, the fragmentation scheme for a particular step in a reaction is chosen on the basis of the nature of the bonding changes associated with that step. Thus, the reactant, transition state, and product are treated consistently throughout the reaction step, leading to an accurate energy barrier for that step. The discontinuity now occurs in describing the energies of reaction intermediates at the transition point between two reaction steps that are treated by two different fragmentation schemes. To address this issue, we propose a systematic procedure for obtaining continuous potential energy curves that are least shifted from their initial positions. The corrected MIM potential energy curves are continuous with activation energies preserved. Following this approach, energy profiles of complex reactions involving large molecular species can be obtained at high levels of theory with a reasonable computational cost.
Collapse
Affiliation(s)
- Ankur Kumar Gupta
- Department of Chemistry , Indiana University , Bloomington , Indiana 47405 , United States
| | - Bishnu Thapa
- Department of Chemistry , Indiana University , Bloomington , Indiana 47405 , United States
| | - Krishnan Raghavachari
- Department of Chemistry , Indiana University , Bloomington , Indiana 47405 , United States
| |
Collapse
|
25
|
Duster AW, Garza CM, Aydintug BO, Negussie MB, Lin H. Adaptive Partitioning QM/MM for Molecular Dynamics Simulations: 6. Proton Transport through a Biological Channel. J Chem Theory Comput 2019; 15:892-905. [DOI: 10.1021/acs.jctc.8b01128] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Adam W. Duster
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Christina M. Garza
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Baris O. Aydintug
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Mikias B. Negussie
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Hai Lin
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| |
Collapse
|
26
|
Hofer TS, de Visser SP. Editorial: Quantum Mechanical/Molecular Mechanical Approaches for the Investigation of Chemical Systems - Recent Developments and Advanced Applications. Front Chem 2018; 6:357. [PMID: 30271768 PMCID: PMC6146044 DOI: 10.3389/fchem.2018.00357] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2018] [Accepted: 07/30/2018] [Indexed: 12/28/2022] Open
Affiliation(s)
- Thomas S Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Sam P de Visser
- School of Chemical Engineering and Analytical Science, Manchester Institute of Biotechnology, The University of Manchester, Manchester, United Kingdom
| |
Collapse
|
27
|
Duster AW, Wang CH, Lin H. Adaptive QM/MM for Molecular Dynamics Simulations: 5. On the Energy-Conserved Permuted Adaptive-Partitioning Schemes. Molecules 2018; 23:E2170. [PMID: 30154373 PMCID: PMC6225285 DOI: 10.3390/molecules23092170] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Revised: 08/23/2018] [Accepted: 08/24/2018] [Indexed: 11/16/2022] Open
Abstract
In combined quantum-mechanical/molecular-mechanical (QM/MM) dynamics simulations, the adaptive-partitioning (AP) schemes reclassify atoms on-the-fly as QM or MM in a smooth manner. This yields a mobile QM subsystem with contents that are continuously updated as needed. Here, we tailor the Hamiltonian adaptive many-body correction (HAMBC) proposed by Boreboom et al. [J. Chem. Theory Comput.2016, 12, 3441] to the permuted AP (PAP) scheme. The treatments lead to the HAMBC-PAP method (HPAP), which both conserves energy and produces accurate solvation structures in the test of "water-in-water" model system.
Collapse
Affiliation(s)
- Adam W Duster
- Department of Chemistry, University of Colorado Denver, Denver, CO 80217, USA.
| | - Chun-Hung Wang
- Department of Chemistry, University of Colorado Denver, Denver, CO 80217, USA.
| | - Hai Lin
- Department of Chemistry, University of Colorado Denver, Denver, CO 80217, USA.
| |
Collapse
|
28
|
Delle Site L. Simulation of Many-Electron Systems That Exchange Matter with the Environment. ADVANCED THEORY AND SIMULATIONS 2018. [DOI: 10.1002/adts.201800056] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Luigi Delle Site
- Institute for Mathematics; Freie Universität Berlin; D-14195 Berlin Germany
| |
Collapse
|
29
|
Zhang YJ, Khorshidi A, Kastlunger G, Peterson AA. The potential for machine learning in hybrid QM/MM calculations. J Chem Phys 2018; 148:241740. [DOI: 10.1063/1.5029879] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Yin-Jia Zhang
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Alireza Khorshidi
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| | - Georg Kastlunger
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| | - Andrew A. Peterson
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
30
|
Hofer TS, Hünenberger PH. Absolute proton hydration free energy, surface potential of water, and redox potential of the hydrogen electrode from first principles: QM/MM MD free-energy simulations of sodium and potassium hydration. J Chem Phys 2018; 148:222814. [DOI: 10.1063/1.5000799] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Thomas S. Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, Centre for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | | |
Collapse
|
31
|
Morzan UN, Alonso de Armiño DJ, Foglia NO, Ramírez F, González Lebrero MC, Scherlis DA, Estrin DA. Spectroscopy in Complex Environments from QM–MM Simulations. Chem Rev 2018; 118:4071-4113. [DOI: 10.1021/acs.chemrev.8b00026] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Diego J. Alonso de Armiño
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Nicolás O. Foglia
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Francisco Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Darío A. Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| |
Collapse
|
32
|
Boereboom JM, Fleurat-Lessard P, Bulo RE. Explicit Solvation Matters: Performance of QM/MM Solvation Models in Nucleophilic Addition. J Chem Theory Comput 2018; 14:1841-1852. [PMID: 29438621 PMCID: PMC6023263 DOI: 10.1021/acs.jctc.7b01206] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
![]()
Nucleophilic
addition onto a carbonyl moiety is strongly affected
by solvent, and correctly simulating this solvent effect is often
beyond the capability of single-scale quantum mechanical (QM) models.
This work explores multiscale approaches for the description of the
reversible and highly solvent-sensitive nucleophilic N|···C=O
bond formation in an Me2N–(CH2)3–CH=O molecule. In the first stage of this work, we
rigorously compare and test four recent quantum mechanical/molecular
mechanical (QM/MM) explicit solvation models, employing a QM description
of water molecules in spherical regions around both the oxygen and
the nitrogen atom of the solute. The accuracy of the models is benchmarked
against a reference QM simulation, focusing on properties of the solvated
Me2N–(CH2)3–CH=O
molecule in its ring-closed form. In the second stage, we select one
of the models (continuous adaptive QM/MM) and use it to obtain a reliable
free energy profile for the N|···C bond formation reaction.
We find that the dual-sphere approach allows the model to accurately
account for solvent reorganization along the entire reaction path.
In contrast, a simple microsolvation model cannot adapt to the changing
conditions and provides an incorrect description of the reaction process.
Collapse
Affiliation(s)
- Jelle M Boereboom
- Inorganic Chemistry and Catalysis Group, Debye Institute for Nanomaterials Science , Utrecht University , Universiteitsweg 99 , 3584 CG Utrecht , The Netherlands
| | - Paul Fleurat-Lessard
- Institut de Chimie Moléculaire de l'Université de Bourgogne (ICMUB, UMR-CNRS 6302) , Université de Bourgogne Franche-Comté , 9 Avenue Alain Savary , 21078 Dijon Cedex , France
| | - Rosa E Bulo
- Inorganic Chemistry and Catalysis Group, Debye Institute for Nanomaterials Science , Utrecht University , Universiteitsweg 99 , 3584 CG Utrecht , The Netherlands
| |
Collapse
|
33
|
Zheng M, Waller MP. Yoink:An interaction-based partitioning API. J Comput Chem 2018; 39:799-806. [DOI: 10.1002/jcc.25146] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2017] [Revised: 11/30/2017] [Accepted: 12/05/2017] [Indexed: 01/08/2023]
Affiliation(s)
- Min Zheng
- Department of Physics and International Centre for Quantum and Molecular Structures; Shanghai University; Shanghai 200444 China
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40; Münster 48149 Germany
| | - Mark P. Waller
- Department of Physics and International Centre for Quantum and Molecular Structures; Shanghai University; Shanghai 200444 China
| |
Collapse
|
34
|
Trnka T, Tvaroška I, Koča J. Automated Training of ReaxFF Reactive Force Fields for Energetics of Enzymatic Reactions. J Chem Theory Comput 2017; 14:291-302. [DOI: 10.1021/acs.jctc.7b00870] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Igor Tvaroška
- Institute
of Chemistry, Slovak Academy of Sciences, 845 38 Bratislava, Slovak Republic
| | | |
Collapse
|
35
|
Zheng M, Moriarty NW, Xu Y, Reimers JR, Afonine PV, Waller MP. Solving the scalability issue in quantum-based refinement: Q|R#1. Acta Crystallogr D Struct Biol 2017; 73:1020-1028. [PMID: 29199981 PMCID: PMC5713877 DOI: 10.1107/s2059798317016746] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2017] [Accepted: 11/20/2017] [Indexed: 12/27/2022] Open
Abstract
Accurately refining biomacromolecules using a quantum-chemical method is challenging because the cost of a quantum-chemical calculation scales approximately as nm, where n is the number of atoms and m (≥3) is based on the quantum method of choice. This fundamental problem means that quantum-chemical calculations become intractable when the size of the system requires more computational resources than are available. In the development of the software package called Q|R, this issue is referred to as Q|R#1. A divide-and-conquer approach has been developed that fragments the atomic model into small manageable pieces in order to solve Q|R#1. Firstly, the atomic model of a crystal structure is analyzed to detect noncovalent interactions between residues, and the results of the analysis are represented as an interaction graph. Secondly, a graph-clustering algorithm is used to partition the interaction graph into a set of clusters in such a way as to minimize disruption to the noncovalent interaction network. Thirdly, the environment surrounding each individual cluster is analyzed and any residue that is interacting with a particular cluster is assigned to the buffer region of that particular cluster. A fragment is defined as a cluster plus its buffer region. The gradients for all atoms from each of the fragments are computed, and only the gradients from each cluster are combined to create the total gradients. A quantum-based refinement is carried out using the total gradients as chemical restraints. In order to validate this interaction graph-based fragmentation approach in Q|R, the entire atomic model of an amyloid cross-β spine crystal structure (PDB entry 2oNA) was refined.
Collapse
Affiliation(s)
- Min Zheng
- International Center for Quantum and Molecular Structures, Shanghai University, Shanghai 200444, People’s Republic of China
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, 48149 Münster, Germany
| | - Nigel W. Moriarty
- Molecular Biosciences and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Yanting Xu
- International Center for Quantum and Molecular Structures, Shanghai University, Shanghai 200444, People’s Republic of China
| | - Jeffrey R. Reimers
- International Center for Quantum and Molecular Structures, Shanghai University, Shanghai 200444, People’s Republic of China
- School of Mathematical and Physical Sciences, University of Technology Sydney, NSW 2007 Australia
| | - Pavel V. Afonine
- International Center for Quantum and Molecular Structures, Shanghai University, Shanghai 200444, People’s Republic of China
- Molecular Biosciences and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Mark P. Waller
- International Center for Quantum and Molecular Structures, Shanghai University, Shanghai 200444, People’s Republic of China
| |
Collapse
|
36
|
Dohn AO, Jónsson EÖ, Levi G, Mortensen JJ, Lopez-Acevedo O, Thygesen KS, Jacobsen KW, Ulstrup J, Henriksen NE, Møller KB, Jónsson H. Grid-Based Projector Augmented Wave (GPAW) Implementation of Quantum Mechanics/Molecular Mechanics (QM/MM) Electrostatic Embedding and Application to a Solvated Diplatinum Complex. J Chem Theory Comput 2017; 13:6010-6022. [PMID: 29083921 DOI: 10.1021/acs.jctc.7b00621] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
A multiscale density functional theory-quantum mechanics/molecular mechanics (DFT-QM/MM) scheme is presented, based on an efficient electrostatic coupling between the electronic density obtained from a grid-based projector augmented wave (GPAW) implementation of density functional theory and a classical potential energy function. The scheme is implemented in a general fashion and can be used with various choices for the descriptions of the QM or MM regions. Tests on H2O clusters, ranging from dimer to decamer show that no systematic energy errors are introduced by the coupling that exceeds the differences in the QM and MM descriptions. Over 1 ns of liquid water, Born-Oppenheimer QM/MM molecular dynamics (MD) are sampled combining 10 parallel simulations, showing consistent liquid water structure over the QM/MM border. The method is applied in extensive parallel MD simulations of an aqueous solution of the diplatinum [Pt2(P2O5H2)4]4- complex (PtPOP), spanning a total time period of roughly half a nanosecond. An average Pt-Pt distance deviating only 0.01 Å from experimental results, and a ground-state Pt-Pt oscillation frequency deviating by <2% from experimental results were obtained. The simulations highlight a remarkable harmonicity of the Pt-Pt oscillation, while also showing clear signs of Pt-H hydrogen bonding and directional coordination of water molecules along the Pt-Pt axis of the complex.
Collapse
Affiliation(s)
- A O Dohn
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland
| | - E Ö Jónsson
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland
| | - G Levi
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - J J Mortensen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - O Lopez-Acevedo
- Department of Applied Physics, Aalto University , 02150 Espoo, Finland
| | - K S Thygesen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - K W Jacobsen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - J Ulstrup
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - N E Henriksen
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - K B Møller
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - H Jónsson
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland.,Department of Applied Physics, Aalto University , 02150 Espoo, Finland
| |
Collapse
|
37
|
Duster AW, Lin H. Restrained Proton Indicator in Combined Quantum-Mechanics/Molecular-Mechanics Dynamics Simulations of Proton Transfer through a Carbon Nanotube. J Phys Chem B 2017; 121:8585-8592. [PMID: 28820594 DOI: 10.1021/acs.jpcb.7b06657] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Recently, a collective variable "proton indicator" was purposed for tracking an excess proton solvated in bulk water in molecular dynamics simulations. In this work, we demonstrate the feasibility of utilizing the position of this proton indicator as a reaction coordinate to model an excess proton migrating through a hydrophobic carbon nanotube in combined quantum-mechanics/molecular-mechanics simulations. Our results indicate that applying a harmonic restraint to the proton indicator in the bulk solvent near the nanotube pore entrance leads to the recruitment of water molecules into the pore. This is consistent with an earlier study that employed a multistate empirical valence bond potential and a different representation (center of excess charge) of the proton. We attribute this water recruitment to the delocalized nature of the solvated proton, which prefers to be in high-dielectric bulk solvent. While water recruitment into the pore is considered an artifact in the present simulations (because of the artificially imposed restraint on the proton), if the proton were naturally restrained, it could assist in building water wires prior to proton transfer through the pore. The potential of mean force for a proton translocation through the water-filled pore was computed by umbrella sampling, where the bias potentials were applied to the proton indicator. The free energy curve and barrier heights agree reasonably with those in the literature. The results suggest that the proton indicator can be used as a reaction coordinate in simulations of proton transport in confined environments.
Collapse
Affiliation(s)
- Adam W Duster
- Chemistry Department, CB 194, University of Colorado Denver , Denver, Colorado 80217, United States
| | - Hai Lin
- Chemistry Department, CB 194, University of Colorado Denver , Denver, Colorado 80217, United States
| |
Collapse
|
38
|
Duster AW, Wang C, Garza CM, Miller DE, Lin H. Adaptive quantum/molecular mechanics: what have we learned, where are we, and where do we go from here? WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1310] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Adam W. Duster
- Chemistry Department University of Colorado Denver Denver CO USA
| | - Chun‐Hung Wang
- Chemistry Department University of Colorado Denver Denver CO USA
| | | | | | - Hai Lin
- Chemistry Department University of Colorado Denver Denver CO USA
| |
Collapse
|
39
|
Affiliation(s)
- Martin J. Field
- Dynamo Team/DYNAMOP Group,
UMR5075, Université Grenoble I, CEA, CNRS, Institut de Biologie Structurale, 71 Avenue des Martyrs, CS 10090, 38044 Grenoble Cedex 9, France
| |
Collapse
|
40
|
de Castro AA, Assis LC, Silva DR, Corrêa S, Assis TM, Gajo GC, Soares FV, Ramalho TC. Computational enzymology for degradation of chemical warfare agents: promising technologies for remediation processes. AIMS Microbiol 2017; 3:108-135. [PMID: 31294152 PMCID: PMC6604975 DOI: 10.3934/microbiol.2017.1.108] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2016] [Accepted: 02/14/2017] [Indexed: 11/18/2022] Open
Abstract
Chemical weapons are a major worldwide problem, since they are inexpensive, easy to produce on a large scale and difficult to detect and control. Among the chemical warfare agents, we can highlight the organophosphorus compounds (OP), which contain the phosphorus element and that have a large number of applications. They affect the central nervous system and can lead to death, so there are a lot of works in order to design new effective antidotes for the intoxication caused by them. The standard treatment includes the use of an anticholinergic combined to a central nervous system depressor and an oxime. Oximes are compounds that reactivate Acetylcholinesterase (AChE), a regulatory enzyme responsible for the transmission of nerve impulses, which is one of the molecular targets most vulnerable to neurotoxic agents. Increasingly, enzymatic treatment becomes a promising alternative; therefore, other enzymes have been studied for the OP degradation function, such as phosphotriesterase (PTE) from bacteria, human serum paraoxonase 1 (HssPON1) and diisopropyl fluorophosphatase (DFPase) that showed significant performances in OP detoxification. The understanding of mechanisms by which enzymes act is of extreme importance for the projection of antidotes for warfare agents, and computational chemistry comes to aid and reduce the time and costs of the process. Molecular Docking, Molecular Dynamics and QM/MM (quantum-mechanics/molecular-mechanics) are techniques used to investigate the molecular interactions between ligands and proteins.
Collapse
Affiliation(s)
| | - Letícia C. Assis
- Department of Chemistry, Federal University of Lavras, 37200-000, Lavras, Brazil
| | - Daniela R. Silva
- Department of Chemistry, Federal University of Lavras, 37200-000, Lavras, Brazil
| | - Silviana Corrêa
- Department of Chemistry, Federal University of Lavras, 37200-000, Lavras, Brazil
| | - Tamiris M. Assis
- Department of Chemistry, Federal University of Lavras, 37200-000, Lavras, Brazil
| | - Giovanna C. Gajo
- Department of Chemistry, Federal University of Lavras, 37200-000, Lavras, Brazil
| | - Flávia V. Soares
- Department of Chemistry, Federal University of Lavras, 37200-000, Lavras, Brazil
| | - Teodorico C. Ramalho
- Department of Chemistry, Federal University of Lavras, 37200-000, Lavras, Brazil
- Center for Basic and Applied Research, Faculty of Informatics and Management, University of Hradec Kralove, Rokitanskeho 62, 50003, Czech Republic
| |
Collapse
|
41
|
A. de Castro A, C. Assis L, R. Silva D, Corrêa S, M. Assis T, C. Gajo G, V. Soares F, C. Ramalho T. Computational enzymology for degradation of chemical warfare agents: promising technologies for remediation processes. AIMS Microbiol 2017. [DOI: 10.3934/microbiol.2017.2.108] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|