Help Privacy Policy Disclaimer
  Advanced SearchBrowse


  A Comparison of QM/MM Simulations with and without the Drude Oscillator Model Based on Hydration Free Energies of Simple Solutes

König, G., Pickard, F. C., Huang, J., Thiel, W., MacKerell, A. D., Brooks, B. R., et al. (2018). A Comparison of QM/MM Simulations with and without the Drude Oscillator Model Based on Hydration Free Energies of Simple Solutes. Molecules, 23(10): 2695. doi:10.3390/molecules23102695.

Item is


show hide
Genre: Journal Article


show Files




König, Gerhard1, 2, 3, Author              
Pickard, Frank C.1, Author
Huang, Jing1, 4, 5, Author
Thiel, Walter3, Author              
MacKerell, Alexander D.4, Author
Brooks, Bernard R.1, Author
York, Darrin M.2, Author
1Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA, ou_persistent22              
2Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA, ou_persistent22              
3Research Department Thiel, Max-Planck-Institut für Kohlenforschung, Max Planck Society, ou_1445590              
4Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA, ou_persistent22              
5School of Life Sciences, Westlake University, 18 Shilongshan Street, Hangzhou 310024, China, ou_persistent22              


Free keywords: hydration free energy; QM/MM; polarization
 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.


Language(s): eng - English
 Dates: 2018-10-152018-09-082018-10-162018-10-19
 Publication Status: Published online
 Pages: 25
 Publishing info: -
 Table of Contents: -
 Rev. Type: Peer
 Identifiers: DOI: 10.3390/molecules23102695
 Degree: -



Legal Case


Project information


Source 1

Title: Molecules
Source Genre: Journal
Publ. Info: Basel : MDPI
Pages: - Volume / Issue: 23 (10) Sequence Number: 2695 Start / End Page: - Identifier: ISSN: 1420-3049
CoNE: https://pure.mpg.de/cone/journals/resource/954925623244