QM/MM Reweighting Free Energy SCF for Geometry Optimization on Extensive Free Energy Surface of Enzymatic Reaction

• Published in: Journal of chemical theory and computation Vol. 8; no. 1; pp. 322 - 334
• Main Authors: Kosugi, Takahiro, Hayashi, Shigehiko
• Format: Journal Article
• Language: English
• Published: United States American Chemical Society 10.01.2012
• Subjects: Biomolecular Systems
• ISSN: 1549-9618; 1549-9626
• DOI: 10.1021/ct2005837

We developed a quantum mechanical/molecular mechanical (QM/MM) free energy geometry optimization method by which the geometry of a quantum chemically treated (QM) molecule is optimized on a free energy surface defined with thermal distribution of the surrounding molecular environment obtained by molecular dynamics simulation with a molecular mechanics (MM) force field. The method called QM/MM reweighting free energy self-consistent field combines a mean field theory of QM/MM free energy geometry optimization developed by Yamamoto (Yamamoto, T. J. Chem. Phys. 2008, 129, 244104) with a reweighting scheme for updating the MM distribution introduced by Hu et al. (Hu, H., et al. J. Chem. Phys. 2008, 128, 034105) and features high computational efficiency suitable for exploring the reaction free energy surface of extensive protein conformational space. The computational efficiency with improved treatment of a long-range electrostatic (ES) interaction using the Ewald summation technique permits one to take into account global conformational relaxation of an entire protein of an enzyme in the free energy geometry optimization of its reaction center. We applied the method to an enzymatic reaction of a substrate complex of psychrophilic α-amylase from Antarctic bacterium Pseudoalteromonas haloplanktis and succeeded in geometry optimizations of the reactant and the product of the catalytic reaction that involve large conformational changes of protein loops adjacent to the reaction center on time scales reaching sub-microseconds. We found that the adjacent loops in the reactant and the product form in different conformations and produce catalytic ES potentials on the reaction center.