Show simple item record

dc.contributor.authorKönig, G.
dc.contributor.authorPickard, F.C.
dc.contributor.authorHuang, J.
dc.date.accessioned2019-06-05T18:28:13Z
dc.date.available2019-06-05T18:28:13Z
dc.date.issued2018
dc.identifier.urihttps://www.scopus.com/inward/record.uri?eid=2-s2.0-85055202116&doi=10.3390%2fmolecules23102695&partnerID=40&md5=4e360b184dd740c1f05722ee918ba433
dc.identifier.urihttp://hdl.handle.net/10713/9405
dc.description.abstractMaintaining 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.en_US
dc.description.urihttps://dx.doi.org/10.3390/molecules23102695en_US
dc.language.isoen-USen_US
dc.publisherMDPI AGen_US
dc.relation.ispartofMolecules
dc.subjectHydration free energyen_US
dc.subjectPolarizationen_US
dc.subjectQM/MMen_US
dc.titleA comparison of QM/MM simulations with and without the Drude oscillator model based on hydration free energies of simple solutesen_US
dc.typeArticleen_US
dc.identifier.doi10.3390/molecules23102695
dc.identifier.pmid30347691


This item appears in the following Collection(s)

Show simple item record