Combining a Dissociative Water Model with a Hybrid QM/MM Approach—A Simulation Strategy for the Study of Proton Transfer Reactions in Solution
Authors:
Open Access
- 30 March 2012
- journal article
- research article
- Published by American Chemical Society (ACS) in Journal of Chemical Theory and Computation
- Vol. 8 (10) , 3586-3595
Abstract
An implementation strategy for a dissociative water potential in hybrid QM/MM simulations is outlined. As the knowledge of the time-dependent topology is crucial for the assignment of solvent molecules to the QM or MM subregion, proton transfer events and the associated change of the molecular composition have to be monitored as the simulation is progressing. A simple and effective update criterion is proposed, which was found to be an efficient tool to identify sustained proton transfer reactions. The resulting topology data enable the application of the dissociative solvent model in QM/MM simulations and serve as a reference for the analysis of time-dependent properties such as the proton hopping rate and the diffusion coefficient. For the latter, an interpolation scheme is proposed linking subsequent proton transfer events into a single diffusive entity. Suitable settings of key parameters for the topology update and the interpolation scheme have been determined by analyzing MD trajectories of an excess proton in water and 1 M HCl. The resulting values for the proton hopping rate and the diffusion coefficient are well within the range estimated by EVB models and CPMD approaches. An investigation of the hydrolytic conversion of As(III) to [As(OH)2]+ serves as an exemplary application of the dissociative model in a QM/MM simulation study.Keywords
This publication has 57 references indexed in Scilit:
- Isobaric−Isothermal Molecular Dynamics Simulations Utilizing Density Functional Theory: An Assessment of the Structure and Density of Water at Near-Ambient ConditionsThe Journal of Physical Chemistry B, 2009
- On the phase diagram of water with density functional theory potentials: The melting temperature of ice Ih with the Perdew–Burke–Ernzerhof and Becke–Lee–Yang–Parr functionalsThe Journal of Chemical Physics, 2009
- QM/MM studies of enzymesCurrent Opinion in Chemical Biology, 2007
- QM/MM: what have we learned, where are we, and where do we go from here?Theoretical Chemistry Accounts, 2006
- Molecular Electronic‐Structure Theory2000
- Potential of mean force for the isomerization of DMF in aqueous solution: a Monte Carlo QM/MM simulation studyJournal of the American Chemical Society, 1993
- A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulationsJournal of Computational Chemistry, 1990
- Unified Approach for Molecular Dynamics and Density-Functional TheoryPhysical Review Letters, 1985
- Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozymeJournal of Molecular Biology, 1976
- Quantisierung als EigenwertproblemAnnalen der Physik, 1926