User Manual Privacy Policy Disclaimer Contact us
  Advanced SearchBrowse




Journal Article

Artificial reaction coordinate "tunneling" in free-energy calculations: the catalytic reaction of RNase H

There are no MPG-Authors available
External Ressource
No external resources are shared
Fulltext (public)
There are no public fulltexts stored in PuRe
Supplementary Material (public)
There is no public supplementary material available

Rosta, E., Woodcock, H. L., Brooks, B. R., & Hummer, G. (2009). Artificial reaction coordinate "tunneling" in free-energy calculations: the catalytic reaction of RNase H. Journal of Computational Chemistry, 30(11), 1634-1641. doi:10.1002/jcc.21312.

Cite as: http://hdl.handle.net/21.11116/0000-0007-4AEF-A
We describe a method for the systematic improvement of reaction coordinates in quantum mechanical/molecular mechanical (QM/MM) calculations of reaction free-energy profiles. In umbrella-sampling free-energy calculations, a biasing potential acting on a chosen reaction coordinate is used to sample the system in reactant, product, and transition states. Sharp, nearly discontinuous changes along the resulting reaction path are used to identify coordinates that are relevant for the reaction but not properly sampled. These degrees of freedom are then included in an extended reaction coordinate. The general formalism is illustrated for the catalytic cleavage of the RNA backbone of an RNA/DNA hybrid duplex by the RNase H enzyme of Bacillus halodurans. We find that in the initial attack of the phosphate diester by water, the oxygen-phosphorus distances alone are not sufficient as reaction coordinates, resulting in substantial hysteresis in the proton degrees of freedom and a barrier that is too low (approximately 10 kcal/mol). If the proton degrees of freedom are included in an extended reaction coordinate, we obtain a barrier of 21.6 kcal/mol consistent with the experimental rates. As the barrier is approached, the attacking water molecule transfers one of its protons to the O1P oxygen of the phosphate group. At the barrier top, the resulting hydroxide ion forms a penta-coordinated phosphate intermediate. The method used to identify important degrees of freedom, and the procedure to optimize the reaction coordinate are general and should be useful both in classical and in QM/MM free-energy calculations.