The hepatitis delta virus (HDV) ribozyme catalyzes a self-cleavage reaction using a combination of nucleobase and metal ion catalysis. Both divalent and monovalent ions can catalyze this reaction, although the rate is slower with monovalent ions alone. Herein, we use quantum mechanical/molecular mechanical (QM/MM) free energy simulations to investigate the mechanism of this ribozyme and to elucidate the roles of the catalytic metal ion. With Mg2+ at the catalytic site, the self-cleavage mechanism is observed to be concerted with a phosphorane-like transition state and a free energy barrier of ∼13 kcal/mol, consistent with free energy barrier values extrapolated from experimental studies. With Na+ at the catalytic site, the mechanism is observed to be sequential, passing through a phosphorane intermediate, with free energy barriers of 2–4 kcal/mol for both steps; moreover, proton transfer from the exocyclic amine of protonated C75 to the nonbridging oxygen of the scissile phosphate occurs to stabilize the phosphorane intermediate in the sequential mechanism. To explain the slower rate observed experimentally with monovalent ions, we hypothesize that the activation of the O2′ nucleophile by deprotonation and orientation is less favorable with Na+ ions than with Mg2+ ions. To explore this hypothesis, we experimentally measure the pKa of O2′ by kinetic and NMR methods and find it to be lower in the presence of divalent ions rather than only monovalent ions. The combined theoretical and experimental results indicate that the catalytic Mg2+ ion may play three key roles: assisting in the activation of the O2′ nucleophile, acidifying the general acid C75, and stabilizing the nonbridging oxygen to prevent proton transfer to it.

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G•U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg 2+ ion to participate in catalysis. Herein, we used molecular dynamics (MD) and X-ray crystallography to characterize the conformation and metal binding characteristics of this base pair in product and precleaved forms. Beginning with a crystal structure of the product form, we observed formation of the reverse G•U wobble during MD trajectories. We also demonstrated that this base pair is compatible with the diffraction data for the product-bound state. During MD trajectories of the product form, Na + ions interacted with the reverse G•U wobble in the RNA active site, and a Mg 2+ ion, introduced in certain trajectories, remained bound at this site. Beginning with a crystal structure of the precleaved form, the reverse G•U wobble with bound Mg 2+ remained intact during MD simulations. When we removed Mg 2+ from the starting precleaved structure, Na + ions interacted with the reverse G•U wobble. In support of the computational results, we observed competition between Na + and Mg 2+ in the precleaved ribozyme crystallographically. Non-linear Poisson-Boltzmann calculations revealed a negatively charged patch near the reverse G•U wobble. This anionic pocket likely serves to bind metal ions and to help shift the pK a of the catalytic nucleobase, C75. Thus, the reverse G•U wobble motif serves to organize two catalytic elements, a metal ion and catalytic nucleobase, within the active site of the HDV ribozyme.RNA is involved in many aspects of biology, where it serves both informational and functional roles (1-3). Indeed, RNA can act as a riboswitch, binding small molecules and † This project was supported by NIH grant R01GM095923 (B.L.G and P.C.B), NIH grant GM56207 (S.H.S.), instrumentation funded by the NSF through grant OCI-0821527, the Purdue University Department of Biochemistry, the Markey Center for Structural Biology and the Purdue University Center for Cancer Research (B.L.G.). * To whom correspondence should be addressed. B.L.G.: telephone (765) 496-6165; fax (765) 494-7897; barbgolden@purdue.edu. S.H.-S. telephone (814) 865-6442; fax (814) 865-2927; shs@chem.psu.edu. P.C.B. telephone (814) 863-3812; fax (814) 865-2927. pcb@chem.psu.edu. ⊥ Present Address: Department of Biochemistry, University of Illinois, 600 S. Mathews Ave., Urbana IL 61801.Supporting Information Available Procedures for MD; protonated cytosine partial charges; additional MD trajectories; metal ion movement during MD; heavy-atom RMSD plots; additional NLPB results. Also provided are tables of crystallographic data collection statistics; crystallographic Na + coordinates for the pre-cleaved ribozyme; and average distances between the reverse G•U wobble and metal ions from MD. This material is available free of charge via the Internet at http://pubs.acs.org. NIH Public Access Author ManuscriptBioche...

Metal ion and nucleobase catalysis are important for ribozyme mechanism, but the extent to which they cooperate is unclear. A crystal structure of the hepatitis delta virus (HDV) ribozyme suggested that the pro-RP oxygen at the scissile phosphate directly coordinates a catalytic Mg2+ ion and is within hydrogen bonding distance of the amine of the general acid C75. Prior studies on the genomic HDV ribozyme, however, showed neither a thio effect nor metal ion rescue using Mn2+. Here, we combine experiment and theory to explore phosphorothioate substitutions at the scissile phosphate. We report significant thio effects at the scissile phosphate and metal ion rescue with Cd2+. Reaction profiles with an SP-phosphorothioate substitution are indistinguishable from those of the unmodified substrate in the presence of Mg2+ or Cd2+, supporting that the pro-SP oxygen does not coordinate metal ions. The RP-phosphorothioate substitution, however, exhibits biphasic kinetics, with the fast-reacting phase displaying a thio effect of up to 5-fold effect and the slow-reacting phase displaying a thio effect of ~1,000-fold. Moreover, the fast- and slow-reacting phases give metal ion rescues in Cd2+ of up to 10- and 330-fold, respectively. The metal ion rescues are unconventional in that they arise from Cd2+ inhibiting the oxo substrate but not the RP substrate. This metal ion rescue suggests a direct interaction of the catalytic metal ion with the pro-RP oxygen, in line with experiments on the antigenomic HDV ribozyme. Experiments without divalent ions, with mutants that interfere with Mg2+ binding, or with C75 deleted suggest that the pro-RP oxygen plays at most a redundant role in positioning C75. Quantum mechanical/molecular mechanical (QM/MM) studies indicate that the metal ion contributes to catalysis by interacting with both the pro-RP oxygen and the nucleophilic 2’- hydroxyl, supporting the experimental findings.

The hepatitis delta virus ribozyme catalyzes an RNA cleavage reaction using a catalytic nucleobase and a divalent metal ion. The catalytic base, C75, serves as a general acid and has a pKa shifted towards neutrality. Less is known about the role of metal ions in the mechanism. A recent crystal structure of the pre-cleavage ribozyme identified a Mg2+ ion that interacts through its partial hydration sphere with the G25•U20 reverse wobble. In addition, this Mg2+ ion is in position to directly coordinate the nucleophile, the 2’-hydroxyl of U(-1), suggesting it can serve as a Lewis acid to facilitate deprotonation of the 2’-hydroxyl. To test the role of the active site Mg2+ ion, we replaced the G25•U20 reverse wobble with an isosteric A25•C20 reverse wobble. This change was found to significantly reduce the negative potential at the active site, as supported by electrostatics calculations, suggesting that active site Mg2+ binding could be adversely affected by the mutation. Kinetic analysis and molecular dynamics of the A25•C20 double mutant suggest that this variant stably folds into an active structure. However, pH-rate profiles of the double mutant are inverted relative to the profiles for wild-type ribozyme, suggesting that the A25•C20 double mutant has lost the active site metal ion. Overall, these studies support a model wherein the partially hydrated Mg2+ positioned at the G25•U20 reverse wobble is catalytic and could serve as a Lewis acid, a Brønsted base, or both to facilitate deprotonation of the nucleophile.