Mail:
Dept. of Chemistry
Ohio State University
100 W. 18th Ave.
Columbus, OH 43210

Office:
412 CBEC

Email:
herbert@
chemistry.ohio-state.edu

Computation of hydration free energies using the multiple environment single system quantum mechanical/molecular mechanical method

G. König, Y. Mei, F. Pickard IV, A. C. Simmonett, B. T. Miller, J. M. Herbert, H. L. Woodcock, B. R. Brooks, and Y. Shao
J. Chem. Theory Comput. 12, 332–344 (2016)

Abstract

A recently-developed MESS-E-QM/MM method (multiple-environment single-system quantum mechanical molecular/mechanical calculations with a Roothaan-step extrapolation) is applied to the computation of hydration free energies for the blind SAMPL4 test set and for twelve small molecules. First, free energy simulations are performed with classical molecular mechanics force field using fixed-geometry solute molecules and explicit TIP3P solvent, then the non-Boltzmann- Bennett method is employed to compute the QM/MM correction (QM/MM-NBB) to the molecular mechanical hydration free energies. For the SAMPL4 set, MESS-E-QM/MM-NBB corrections to the hydration free energy can be obtained 2 or 3 orders of magnitude faster than fully converged QM/MM-NBB corrections, and, on average, the hydration free energies predicted with MESS-E-QM/MM-NBB fall within 0.10–0.20 kcal/mol of fully-converged QM/MM-NBB results. Out of five density functionals (BLYP, B3LYP, PBE0, M06-2X, and ωB97X-D), the BLYP functional is found to be most compatible with the TIP3P solvent model and yields the most accurate hydration free energies against experimental values for solute molecules included in this study.

[DOI] [PDF]
[Supporting Info]
Last modified January 12, 2016. Proudly powered by Words. By which we mean, hand-written HTML.

This page best viewed with a browser