pubmed.ncbi.nlm.nih.gov

Enzymatic Kinetic Isotope Effects from Path-Integral Free Energy Perturbation Theory - PubMed

Enzymatic Kinetic Isotope Effects from Path-Integral Free Energy Perturbation Theory

J Gao. Methods Enzymol. 2016.

Abstract

Path-integral free energy perturbation (PI-FEP) theory is presented to directly determine the ratio of quantum mechanical partition functions of different isotopologs in a single simulation. Furthermore, a double averaging strategy is used to carry out the practical simulation, separating the quantum mechanical path integral exactly into two separate calculations, one corresponding to a classical molecular dynamics simulation of the centroid coordinates, and another involving free-particle path-integral sampling over the classical, centroid positions. An integrated centroid path-integral free energy perturbation and umbrella sampling (PI-FEP/UM, or simply, PI-FEP) method along with bisection sampling was summarized, which provides an accurate and fast convergent method for computing kinetic isotope effects for chemical reactions in solution and in enzymes. The PI-FEP method is illustrated by a number of applications, to highlight the computational precision and accuracy, the rule of geometrical mean in kinetic isotope effects, enhanced nuclear quantum effects in enzyme catalysis, and protein dynamics on temperature dependence of kinetic isotope effects.

Keywords: Combined QM/MM; Dual-level potential; Enzyme kinetics; Kinetic isotope effects; Nuclear quantum effects; PI-FEP; Path-integral simulation.

© 2016 Elsevier Inc. All rights reserved.

PubMed Disclaimer

Figures

Fig. 1
Fig. 1

Computed nuclear quantum effects along the potential of mean force for the decarboxylation reaction of N-methyl picolinate ion in water. 13C denotes isotope substitution at the carboxyl carbon position.

Fig. 2
Fig. 2

Computed ratio of quantum partition function between 12C and 13C isotope substitutions at the carboxyl carbon position along the centroid reaction coordinate for the decarboxylation of N-methyl picolinate in water.

Fig. 3
Fig. 3

Classical (blue) and quantum mechanical potential of mean force for the proton (red) and deuteron (green) transfer from nitroethane to acetate ion in water (dashed curves), and to Asp402 (solid curves) in the active site of nitroalkane oxidase.

Fig. 4
Fig. 4

Superposition comparison of snapshots of the Michaelis complex (blue) and transition state (brown) structures for the wild-type enzyme dihydrofolate reductase (left) and the M42W/G121V double mutant (right).

Fig. 5
Fig. 5

Computed and experimental H/D primary kinetic isotope effects for the hydride transfer process catalyzed by wild-type and M42W/G121V mutant dihydrofolate reductase at different temperatures.

Similar articles

Cited by

References

    1. Alhambra C, Corchado J, Sanchez ML, Garcia-Viloca M, Gao J, Truhlar DG. Canonical variational theory for enzyme kinetics with the protein mean force and multidimensional quantum mechanical tunneling dynamics. Theory and application to liver alcohol dehydrogenase. The Journal of Physical Chemistry. B. 2001;105:11326–11340.
    1. Bernasconi CF. The principle of nonperfect synchronization: More than a qualitative concept? Accounts of Chemical Research. 1992;25:9–16.
    1. Bhabha G, Lee J, Ekiert DC, Gam J, Wilson IA, et al. A dynamic knockout reveals that conformational fluctuations influence the chemical step of enzyme catalysis. Science. 2011;332:234–238. - PMC - PubMed
    1. Bigeleisen J. Statistical mechanics of isotopic systems with small quantum corrections. I. General considerations and the rule of the geometric mean. The Journal of Chemical Physics. 1955;23:2264.
    1. Boehr DD, McElheny D, Dyson HJ, Wright PE. The dynamic energy landscape of dihydrofolate reductase catalysis. Science. 2006;313:1638–1642. - PubMed

Publication types

MeSH terms

Substances