1
|
Bonfrate S, Ferré N, Huix-Rotllant M. Analytic Gradients for the Electrostatic Embedding QM/MM Model in Periodic Boundary Conditions Using Particle-Mesh Ewald Sums and Electrostatic Potential Fitted Charge Operators. J Chem Theory Comput 2024; 20:4338-4349. [PMID: 38712506 DOI: 10.1021/acs.jctc.4c00201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2024]
Abstract
Long-range electrostatic effects are fundamental for describing chemical reactivity in the condensed phase. Here, we present the methodology of an efficient quantum mechanical/molecular mechanical (QM/MM) model in periodic boundary conditions (PBC) compatible with QM/MM boundaries at chemical bonds. The method combines electrostatic potential fitted charge operators and electrostatic potentials derived from the smooth particle-mesh Ewald (PME) sum approach. The total energy and its analytic first derivatives with respect to QM, MM, and lattice vectors allow QM/MM molecular dynamics (MD) in the most common thermodynamic ensembles. We demonstrate the robustness of the method by performing a QM/MM MD equilibration of methanol in water. We simulate the cis/trans isomerization free-energy profiles in water of proline amino acid and a proline-containing oligopeptide, showing a correct description of the reaction barrier. Our PBC-compatible QM/MM model can efficiently be used to study the chemical reactivity in the condensed phase and enzymatic catalysis.
Collapse
Affiliation(s)
| | - Nicolas Ferré
- Aix-Marseille Univ, CNRS, ICR, Marseille 13013, France
| | | |
Collapse
|
2
|
Gallmetzer HG, Hofer TS. Probing the range of applicability of structure- and energy-adjusted QM/MM link bonds II: Optimized link bond parameters for density functional tight binding approaches. J Comput Chem 2022; 43:746-756. [PMID: 35239208 PMCID: PMC9314059 DOI: 10.1002/jcc.26830] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2021] [Revised: 02/08/2022] [Accepted: 02/10/2022] [Indexed: 11/21/2022]
Abstract
Optimized link bond parameters for the Cα Cβ bond of 22 different capped amino acid model systems have been determined at SCC DFTB/mio (self-consistent charge density functional tight-binding), SCC DFTB/3ob and GFNn-xTB (n = 0, 1, and 2) level in conjunction with the AMBER 99SB, 14SB, and 19B force fields. The resulting parameter sets have been compared to newly calculated reference data obtained via resolution-of-identity 2nd order Møller-Plesset perturbation theory. The data collected in this work suggests that the optimized values in this study provide a more suitable setup of the QM/MM link bonds compared to the use of a single global setting applied to every amino acid fragmented by the QM/MM interface. The results also imply that a transfer of the ideal link bond settings between different levels of theory is not advised. In contrast, virtually identical parameters were obtained in calculations employing different variants of the AMBER force field. Considering the increasing success of tight binding based approaches being inter alia a results of their exceptional accuracy/effort ratio the provided collection of link atoms parameters provides a valuable resource for QM/MM studies of biomacromolecular systems as demonstrated in an exemplary QM/MM MD simulation of the β-amyloid/Zn2+ complex.
Collapse
Affiliation(s)
- Hans Georg Gallmetzer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical ChemistryCenter for Chemistry and Biomedicine, University of InnsbruckInnsbruckAustria
| | - Thomas S. Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical ChemistryCenter for Chemistry and Biomedicine, University of InnsbruckInnsbruckAustria
| |
Collapse
|
3
|
Gheorghiu A, Coveney PV, Arabi AA. The influence of base pair tautomerism on single point mutations in aqueous DNA. Interface Focus 2020; 10:20190120. [PMID: 33178413 PMCID: PMC7653342 DOI: 10.1098/rsfs.2019.0120] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/24/2020] [Indexed: 12/11/2022] Open
Abstract
The relationship between base pair hydrogen bond proton transfer and the rate of spontaneous single point mutations at ambient temperatures and pressures in aqueous DNA is investigated. By using an ensemble-based multiscale computational modelling method, statistically robust rates of proton transfer for the A:T and G:C base pairs within a solvated DNA dodecamer are calculated. Several different proton transfer pathways are observed within the same base pair. It is shown that, in G:C, the double proton transfer tautomer is preferred, while the single proton transfer process is favoured in A:T. The reported range of rate coefficients for double proton transfer is consistent with recent experimental data. Notwithstanding the approximately 1000 times more common presence of single proton transfer products from A:T, observationally there is bias towards G:C to A:T mutations in a wide range of living organisms. We infer that the double proton transfer reactions between G:C base pairs have a negligible contribution towards this bias for the following reasons: (i) the maximum half-life of the G*:C* tautomer is in the range of picoseconds, which is significantly smaller than the milliseconds it takes for DNA to unwind during replication, (ii) statistically, the majority of G*:C* tautomers revert back to their canonical forms through a barrierless process, and (iii) the thermodynamic instability of the tautomers with respect to the canonical base pairs. Through similar reasoning, we also deduce that proton transfer in the A:T base pair does not contribute to single point mutations in DNA.
Collapse
Affiliation(s)
- A Gheorghiu
- Centre for Computational Science, University College London, London, UK
| | - P V Coveney
- Centre for Computational Science, University College London, London, UK.,Informatics Institute, University of Amsterdam, Amsterdam, The Netherlands
| | - A A Arabi
- Centre for Computational Science, University College London, London, UK.,College of Medicine and Health Sciences, Biochemistry Department, United Arab Emirates University, PO Box 17666, Al Ain, United Arab Emirates
| |
Collapse
|
4
|
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
|