Coupling between protein and reaction dynamics in enzymatic processes:: Application of Grote-Hynes theory to catechol O-methyltransferase

The generalized Langevin equation (GLE)-based Grote-Hynes (GH) theory is used to calculate the transmission coefficients, kappa, for the methyl transfer from S-adenosylmethionine to catecholate both in aqueous solution and in the catechol O-methyltransferase active site. Values of kappa, which measures the deviation of the rate constants from the Transition State Theory (TST) predictions, are obtained by means of rare event molecular dynamics simulations. The results are 0.62 +/- 0.04 and 0.83 +/- 0.03 for the aqueous and enzymatic environments, respectively, while the Grote-Hynes predictions are 0.58 +/- 0.09 and 0.89 +/- 0.03, respectively. The Kramers theory estimates are much smaller, about 0.01 and 0.1, respectively. Thus, the enzymatic transmission coefficient is closer to TST predictions than the value obtained in solution. In addition, our results show that the enzymatic coefficient is also closer to its nonadiabatic (or frozen environment) limit than is the solution coefficient. These findings can be understood considering that, during the passage over the barrier top, there is a smaller coupling between the reactive system and the environment in the enzyme than in solution, as well as a smaller reorganization suffered by the enzyme. Analysis of the transition state friction kernel leads to the identification of some key vibrational modes governing the coupling between the two different environments and the reacting solute in the transition state region and insights on their relevance for the reaction dynamics' influence on the transmission coefficient.