Free energy calculation on enzyme reactions with an efficient iterative procedure to determine minimum energy paths on a combined ab initio QM/MM potential energy surface

Journal Article (Journal Article)

A new practical approach to studying enzyme reactions by combining ab initio QM/MM calculations with free energy perturbation is presented. An efficient iterative optimization procedure has been developed to determine optimized structures and minimum energy paths for a system with thousands of atoms on the ab initio QM/MM potential: the small QM sub-system is optimized using a quasi-Newton minimizer in redundant internal coordinates with ab initio QM/MM calculations, while the large MM sub-system is minimized by the truncated Newton method in Cartesian coordinates with only molecular mechanical calculations. The above two optimization procedures are performed iteratively until they converge. With the determined minimum energy paths, free energy perturbation calculations are carried out to determine the change in free energy along the reaction coordinate. Critical to the success of the iterative optimization procedure and the free energy calculations is the smooth connection between the QM and MM regions provided by a recently proposed pseudobond QM/MM approach [J. Chem. Phys. 110, 46 (1999)], The methods have been demonstrated by studying the initial proton transfer step in the reaction catalyzed by the enzyme triosephosphate isomerase (TIM). © 2000 American Institute of Physics.

Full Text

Duke Authors

Cited Authors

  • Zhang, Y; Liu, H; Yang, W

Published Date

  • February 22, 2000

Published In

Volume / Issue

  • 112 / 8

Start / End Page

  • 3483 - 3492

International Standard Serial Number (ISSN)

  • 0021-9606

Digital Object Identifier (DOI)

  • 10.1063/1.480503

Citation Source

  • Scopus