A comparison of QM/MM simulations with and without the Drude oscillator model based on hydration free energies of simple solutes
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. Copyright 2018 by the authors.Identifier to cite or link to this item
https://www.scopus.com/inward/record.uri?eid=2-s2.0-85055202116&doi=10.3390%2fmolecules23102695&partnerID=40&md5=4e360b184dd740c1f05722ee918ba433; http://hdl.handle.net/10713/9405ae974a485f413a2113503eed53cd6c53
10.3390/molecules23102695
Scopus Count
Collections
Related articles
- Computation of Hydration Free Energies Using the Multiple Environment Single System Quantum Mechanical/Molecular Mechanical Method.
- Authors: König G, Mei Y, Pickard FC 4th, Simmonett AC, Miller BT, Herbert JM, Woodcock HL, Brooks BR, Shao Y
- Issue date: 2016 Jan 12
- Combining the polarizable Drude force field with a continuum electrostatic Poisson-Boltzmann implicit solvation model.
- Authors: Aleksandrov A, Lin FY, Roux B, MacKerell AD Jr
- Issue date: 2018 Aug 15
- 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.
- Authors: König G, Pickard FC 4th, Huang J, Simmonett AC, Tofoleanu F, Lee J, Dral PO, Prasad S, Jones M, Shao Y, Thiel W, Brooks BR
- Issue date: 2016 Nov
- Predicting hydration free energies with a hybrid QM/MM approach: an evaluation of implicit and explicit solvation models in SAMPL4.
- Authors: König G, Pickard FC 4th, Mei Y, Brooks BR
- Issue date: 2014 Mar
- Polarizable Force Field for Molecular Ions Based on the Classical Drude Oscillator.
- Authors: Lin FY, Lopes PEM, Harder E, Roux B, MacKerell AD Jr
- Issue date: 2018 May 29