Solvent D2O and (18)O kinetic isotope effects on RNA 2'-O-transphosphorylation catalyzed by Zn(2+) demonstrate an altered transition state relative to specific base catalysis. A recent model from DFT calculations involving inner sphere coordination to the non-bridging and leaving group oxygens is consistent with the data.

Biological catalysis involves interactions distant from the site of chemistry that can position the substrate for reaction. Catalysis of RNA 2′-O-transphosphorylation by the hepatitis delta virus (HDV) ribozyme is sensitive to the identity of the N(–1) nucleotide flanking the reactive phosphoryl group. However, the interactions that affect the conformation of this position, and in turn the 2′O nucleophile, are unclear. Here, we describe the application of multiple substrate internal competition kinetic analyses to understand how the N(–1) nucleobase contributes to HDV catalysis and test the utility of this approach for RNA structure–function studies. Internal competition reactions containing all four substrate sequence variants at the N(–1) position in reactions using ribozyme active site mutations at A77 and A78 were used to test a proposed base-pairing interaction. Mutants A78U, A78G, and A79G retain significant catalytic activity but do not alter the specificity for the N(–1) nucleobase. Effects of nucleobase analog substitutions at N(–1) indicate that U is preferred due to the ability to donate an H-bond in the Watson–Crick face and avoid minor groove steric clash. The results provide information essential for evaluating models of the HDV active site and illustrate multiple substrate kinetic analyses as a practical approach for characterizing structure–function relationships in RNA reactions.

Although there have been great strides in defining the mechanisms of RNA strand cleavage by 2'-O-transphosphorylation, long-standing questions remain. How do different catalytic modes such as acid/base and metal ion catalysis influence transition state charge distribution? Does the large rate enhancement characteristic of biological catalysis result in different transition states relative to solution reactions? Answering these questions is important for understanding biological catalysis in general, and revealing principles for designing small molecule inhibitors. Recent application of linear free energy relationships and kinetic isotope effects together with multi-scale computational simulations are providing tentative answers to these questions for this fundamentally important class of phosphoryl transfer reactions.

Ribonuclease A (RNase A) is a widely used model system illustrating basic principles of enzyme catalysis. Yet, ambiguity regarding the transition state it stabilizes precludes a complete understanding of its mechanism. Here, we present kinetic isotope effect analyses that provide direct evidence for a late, anionic transition state for RNase A catalyzed RNA 2′‐O‐transphosphorylation similar to catalysis by hydroxide in solution, but with dramatically less charge accumulation on the 5′O. Molecular simulations indicate stabilization of negative charge on the NPOs by H‐bonding with a protonated His12 and Lys41. Quantum chemical calculations consistent with the experimental KIEs provide a unified mechanistic interpretation whereby P‐O5′ bond cleavage is concerted with proton transfer from a protonated His119. Thus, theory and experiment provide strong support that RNase A exploits both direct electrostatic stabilization of the anionic transition state and general acid catalysis, resulting in an altered transition state relative to that observed in solution. This work was supported by grants NIH GM096000 to M.E.H., NSF MCB0717850 to V.E.A., NIH AI081987 to J.A.P. and NIH GM064288 to D.M.Y. D.L.K was supported by NIH training grant T32‐GM008056.

Enzymes function by stabilizing reaction transition states; therefore, comparison of the transition states of enzymatic and nonenzymatic model reactions can provide insight into biological catalysis. Catalysis of RNA 2′- O -transphosphorylation by ribonuclease A is proposed to involve electrostatic stabilization and acid/base catalysis, although the structure of the rate-limiting transition state is uncertain. Here, we describe coordinated kinetic isotope effect (KIE) analyses, molecular dynamics simulations, and quantum mechanical calculations to model the transition state and mechanism of RNase A. Comparison of the 18 O KIEs on the 2′O nucleophile, 5′O leaving group, and nonbridging phosphoryl oxygens for RNase A to values observed for hydronium- or hydroxide-catalyzed reactions indicate a late anionic transition state. Molecular dynamics simulations using an anionic phosphorane transition state mimic suggest that H-bonding by protonated His12 and Lys41 stabilizes the transition state by neutralizing the negative charge on the nonbridging phosphoryl oxygens. Quantum mechanical calculations consistent with the experimental KIEs indicate that expulsion of the 5′O remains an integral feature of the rate-limiting step both on and off the enzyme. Electrostatic interactions with positively charged amino acid site chains (His12/Lys41), together with proton transfer from His119, render departure of the 5′O less advanced compared with the solution reaction and stabilize charge buildup in the transition state. The ability to obtain a chemically detailed description of 2′- O -transphosphorylation transition states provides an opportunity to advance our understanding of biological catalysis significantly by determining how the catalytic modes and active site environments of phosphoryl transferases influence transition state structure.

To better understand the interactions between catalysts and transition states during RNA strand cleavage, primary (18)O kinetic isotope effects (KIEs) and solvent D(2)O isotope effects were measured to probe the mechanism of base-catalyzed 2'-O-transphosphorylation of the RNA dinucleotide 5'-UpG-3'. The observed (18)O KIEs for the nucleophilic 2'-O and in the 5'-O leaving group at pH 14 are both large relative to reactions of phosphodiesters with good leaving groups, indicating that the reaction catalyzed by hydroxide has a transition state (TS) with advanced phosphorus-oxygen bond fission to the leaving group ((18)k(LG) = 1.034 +/- 0.004) and phosphorus-nucleophile bond formation ((18)k(NUC) = 0.984 +/- 0.004). A breakpoint in the pH dependence of the 2'-O-transphosphorylation rate to a pH independent phase above pH 13 has been attributed to the pK(a) of the 2'-OH nucleophile. A smaller nucleophile KIE is observed at pH 12 ((18)k(NUC) = 0.995 +/- 0.004) that is interpreted as the combined effect of the equilibrium isotope effect (ca. 1.02) on deprotonation of the 2'-hydroxyl nucleophile and the intrinsic KIE on the nucleophilic addition step (ca. 0.981). An alternative mechanism in which the hydroxide ion acts as a general base is considered unlikely given the lack of a solvent deuterium isotope effect above the breakpoint in the pH versus rate profile. These results represent the first direct analysis of the transition state for RNA strand cleavage. The primary (18)O KIE results and the lack of a kinetic solvent deuterium isotope effect together provide strong evidence for a late transition state and 2'-O nucleophile activation by specific base catalysis.