Are Mixed Explicit/Implicit Solvation Models Reliable for Studying Phosphate Hydrolysis? A Comparative Study of Continuum, Explicit and Mixed Solvation Models
- 4 May 2009
- journal article
- research article
- Published by Wiley in Chemphyschem
- Vol. 10 (7), 1125-1134
- https://doi.org/10.1002/cphc.200800753
Abstract
Are all solvation models equal? An in‐depth comparison of mixed implicit/explicit solvation models shows that while both the implicit and explicit (QM/MM–FEP) solvation models can reproduce activation free energies for phosphate diester hydrolysis in solution with high accuracy, the use of a mixed solvation model tends to be unreliable for modelling phosphate hydrolysis in solution. Phosphate hydrolysis is ubiquitous in biology. However, despite intensive research on this class of reactions, the precise nature of the reaction mechanism remains controversial. Herein, we have examined the hydrolysis of three homologous phosphate diesters. The solvation free energy was simulated by means of either an implicit solvation model (COSMO), hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM–FEP) or a mixed solvation model in which N water molecules were explicitly included in the ab initio description of the reacting system (where N=1–3), with the remainder of the solvent being implicitly modelled as a continuum. Here, both COSMO and QM/MM–FEP reproduce ΔGobs within an error of about 1 kcal mol−1. However, we demonstrate that in order to obtain any kind of reliable results from a mixed model, it is essential to carefully select the explicit water molecules from short QM/MM runs that act as a model for the true infinite system. Additionally, the mixed models tend to be increasingly unstable and miss larger entropic contributions as more explicit water molecules are placed into the system. Thus, our analysis indicates that this approach provides an unreliable way for modelling phosphate hydrolysis in solution.Keywords
This publication has 77 references indexed in Scilit:
- Progress in Ab Initio QM/MM Free-Energy Simulations of Electrostatic Energies in Proteins: Accelerated QM/MM Studies of pKa, Redox Reactions and Solvation Free EnergiesThe Journal of Physical Chemistry B, 2008
- Associative Versus Dissociative Mechanisms of Phosphate Monoester Hydrolysis: On the Interpretation of Activation EntropiesChemphyschem, 2008
- Solvation free energies of molecules. The most stable anionic tautomers of uracilPhysical Chemistry Chemical Physics, 2008
- Accelerating QM/MM Free Energy Calculations: Representing the Surroundings by an Updated Mean Charge DistributionThe Journal of Physical Chemistry B, 2008
- Theoretical Study of Specific Hydrogen-Bonding Effects on the Bridging P−OR Bond Strength of Phosphate Monoester DianionsThe Journal of Physical Chemistry B, 2008
- Transesterification Thio Effects of Phosphate Diesters: Free Energy Barriers and Kinetic and Equilibrium Isotope Effects from Density-Functional TheoryBiochemistry, 2006
- Density Functional Study of the In-Line Mechanism of Methanolysis of Cyclic Phosphate and Thiophosphate Esters in Solution: Insight into Thio Effects in RNA TransesterificationThe Journal of Physical Chemistry B, 2005
- The Substrate-Assisted General Base Catalysis Model for Phosphate Monoester Hydrolysis: Evaluation Using Reactivity ComparisonsJournal of the American Chemical Society, 2000
- The quest for free metaphosphate in solution. Racemization at phosphorus in the transfer of the phospho group from aryl phosphate monoesters to tert-butyl alcohol in acetonitrile or in tert-butyl alcoholJournal of the American Chemical Society, 1988
- Monomeric metaphosphatesChemical Reviews, 1981