Probing general base catalysis in the hammerhead ribozyme

Recent structural and computational studies have shed new light on the catalytic mechanism and active site structure of the RNA cleaving hammerhead ribozyme. Consequently, specific ribozyme functional groups have been hypothesized to be directly involved in general/acid base catalysis. In order to test this hypothesis, we have developed an affinity label to identify the functional general base in the S. mansoni hammerhead ribozyme. The ribozyme was reacted with a substrate analogue bearing a 2'-bromoacetamide group in place of the nucleophilic 2'-hydroxyl group which would normally be deprotonated by a general base. The electrophilic 2'-bromoacetamide group is poised to alkylate the general base, which is subsequently identified by footprinting analysis. Herein, we demonstrate alkylation of N1 of G12 in the hammerhead ribozyme in a pH and [Mg(2+)] dependent manner that is consistent with the native cleavage reaction. These results provide substantial evidence that deprotonated N1 of G12 functions directly as a general base in the hammerhead ribozyme; moreover, our experiments provide evidence that the pKa of G12 is perturbed downward in the context of the active site structure. We also observed other pH-independent alkylations, which do not appear to reflect the catalytic mechanism, but offer further insight into ribozyme conformation and structure.     

Quantum mechanical/molecular mechanical simulation study of the mechanism of hairpin ribozyme catalysis

The molecular mechanism of hairpin ribozyme catalysis is studied with molecular dynamics simulations using a combined quantum mechanical and molecular mechanical (QM/MM) potential with a recently developed semiempirical AM1/d-PhoT model for phosphoryl transfer reactions. Simulations are used to derive one- and two-dimensional potentials of mean force to examine specific reaction paths and assess the feasibility of proposed general acid and base mechanisms. Density-functional calculations of truncated active site models provide complementary insight to the simulation results. Key factors utilized by the hairpin ribozyme to enhance the rate of transphosphorylation are presented, and the roles of A38 and G8 as general acid and base catalysts are discussed. The computational results are consistent with available experimental data, provide support for a general acid/base mechanism played by functional groups on the nucleobases, and offer important insight into the ability of RNA to act as a catalyst without explicit participation by divalent metal ions.     

Solvent structure and hammerhead ribozyme catalysis

Although the hammerhead ribozyme is regarded as a prototype for understanding RNA catalysis, the mechanistic roles of associated metal ions and water molecules in the cleavage reaction remain controversial. We have investigated the catalytic potential of observed divalent metal ions and water molecules bound to a 2 A structure of the full-length hammerhead ribozyme by using X-ray crystallography in combination with molecular dynamics simulations. A single Mn(2+) is observed to bind directly to the A9 phosphate in the active site, accompanying a hydrogen-bond network involving a well-ordered water molecule spanning N1 of G12 (the general base) and 2'-O of G8 (previously implicated in general acid catalysis) that we propose, based on molecular dynamics calculations, facilitates proton transfer in the cleavage reaction. Phosphate-bridging metal interactions and other mechanistic hypotheses are also tested with this approach.     

An active-site guanine participates in glmS ribozyme catalysis in its protonated state

Active-site guanines that occupy similar positions have been proposed to serve as general base catalysts in hammerhead, hairpin, and glmS ribozymes, but no specific roles for these guanines have been demonstrated conclusively. Structural studies place G33(N1) of the glmS ribozyme of Bacillus anthracis within hydrogen-bonding distance of the 2'-OH nucleophile. Apparent pK(a) values determined from the pH dependence of cleavage kinetics for wild-type and mutant glmS ribozymes do not support a role for G33, or any other active-site guanine, in general base catalysis. Furthermore, discrepancies between apparent pK(a) values obtained from functional assays and microscopic pK(a) values obtained from pH-fluorescence profiles with ribozymes containing a fluorescent guanosine analogue, 8-azaguanosine, at position 33 suggest that the pH-dependent step in catalysis does not involve G33 deprotonation. These results point to an alternative model in which G33(N1) in its neutral, protonated form donates a hydrogen bond to stabilize the transition state.     

Atom‐Specific Mutagenesis Reveals Structural and Catalytic Roles for an Active‐Site Adenosine and Hydrated Mg2+ in Pistol Ribozymes

The pistol RNA motif represents a new class of self‐cleaving ribozymes of yet unknown biological function. Our recent crystal structure of a pre‐catalytic state of this RNA shows guanosine G40 and adenosine A32 close to the G53–U54 cleavage site. While the N1 of G40 is within 3.4 Å of the modeled G53 2′‐OH group that attacks the scissile phosphate, thus suggesting a direct role in general acid–base catalysis, the function of A32 is less clear. We present evidence from atom‐specific mutagenesis that neither the N1 nor N3 base positions of A32 are involved in catalysis. By contrast, the ribose 2′‐OH of A32 seems crucial for the proper positioning of G40 through a H‐bond network that involves G42 as a bridging unit between A32 and G40. We also found that disruption of the inner‐sphere coordination of the active‐site Mg²⁺ cation to N7 of G33 makes the ribozyme drastically slower. A mechanistic proposal is suggested, with A32 playing a structural role and hydrated Mg²⁺ playing a catalytic role in cleavage.     

Rescue of an abasic hairpin ribozyme by cationic nucleobases: evidence for a novel mechanism of RNA catalysis

The hairpin ribozyme catalyzes a reversible phosphodiester cleavage reaction. We examined the roles of conserved nucleobases in catalysis using an abasic ribozyme rescue strategy. Loss of the active site G8 nucleobase reduced the cleavage rate constant by 350-fold while loss of A9 and A10 nucleobases reduced activity less than 10-fold. Certain heterocyclic amines restored partial activity when provided in solution to the variant lacking G8. Heterocyclic amines that were capable of rescue shared the exocyclic amine and cyclic amide in common with the Watson-Crick hydrogen bonding face of guanine. In contrast to the shallow pH dependence of unmodified ribozyme activity, rescue activity increased sharply with decreasing pH. These results support a novel model for RNA catalysis in which a cationic nucleobase contributes electrostatic stabilization to negative charge developing in the transition state.     

An unusual pH-independent and metal-ion-independent mechanism for hairpin ribozyme catalysis

Background: Hairpin ribozymes (RNA enzymes) catalyze the same chemical reaction as ribonuclease A and yet RNAs do not usually have functional groups analogous to the catalytically essential histidine and lysine sidechains of protein ribonucleases. Some RNA enzymes appear to recruit metal ions to act as Lewis acids in charge stabilization and metal-bound hydroxide for general base catalysis, but it has been reported that the hairpin ribozyme functions in the presence of metal ion chelators. This led us to investigate whether the hairpin ribozyme exploits a metal-ion-independent catalytic strategy.Results: Substitution of sulfur for nonbridging oxygens of the reactive phosphate of the hairpin ribozyme has small, stereospecific and metal-ionindependent effects on cleavage and ligation mediated by this ribozyme. Cobalt hexammine, an exchange-inert metal complex, supports full hairpin ribozyme activity, and the ribozyme's catalytic rate constants display only a shallow dependence on pH.Conclusions: Direct metal ion coordination to phosphate oxygens is not essential for hairpin ribozyme catalysis and metal-bound hydroxide does not serve as the general base in this catalysis. Several models might account for the unusual pH and metal ion independence: hairpin cleavage and ligation might be limited by a slow conformational change; a pH-independent or metalcation-independent chemical step, such as breaking the 5′ oxygen-phosphorus bond, might be rate determining; or finally, functional groups within the ribozyme might participate directly in catalytic chemistry. Whichever the case, the hairpin ribozyme appears to employ a unique strategy for RNA catalysis.     

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.     

Crucial Roles of Two Hydrated Mg2+ Ions in Reaction Catalysis of the Pistol Ribozyme

Pistol ribozymes constitute a new class of small self‐cleaving RNAs. Crystal structures have been solved, providing three‐dimensional snapshots along the reaction coordinate of pistol phosphodiester cleavage, corresponding to the pre‐catalytic state, a vanadate mimic of the transition state, and the product. The results led to the proposed underlying chemical mechanism. Importantly, a hydrated Mg²⁺ ion remains innersphere‐coordinated to N7 of G33 in all three states, and is consistent with its likely role as acid in general acid base catalysis (δ and β catalysis). Strikingly, the new structures shed light on a second hydrated Mg²⁺ ion that approaches the scissile phosphate from its binding site in the pre‐cleavage state to reach out for water‐mediated hydrogen bonding in the cyclophosphate product. The major role of the second Mg²⁺ ion appears to be the stabilization of product conformation. This study delivers a mechanistic understanding of ribozyme‐catalyzed backbone cleavage.     

Density functional theory investigation on the mechanism of the hepatitis delta virus ribozyme

Density functional theory methods have been used to investigate the hepatitis delta virus (HDV) ribozyme and its catalyzed phosphodiester cleavage. In particular, the effects of the environment's polarity and/or specific hydrogen-bond interactions on the proton affinity of the active site cytosine's N3 ring center have been considered. In addition, the basicities of possible hydrated Mg2+ ion species were also examined. The mechanism previously proposed for the HDV ribozyme in which the active site cytosine (C75) is protonated and thus acts as an acid while the Mg2+ species acts as the complementary base was then investigated. The possible role of tautomerization of C75 is also discussed.     

A unique mechanism for RNA catalysis: the role of metal cofactors in hairpin ribozyme cleavage

Background: Ribozymes are biological catalysts that promote the hydrolysis and transesterification of phosphate diesters of RNA. They typically require divalent magnesium ions for activation, although it has proven difficult to differentiate structural from catalytic roles for the magnesium ions and to identify the molecular mechanism of catalysis. Direct inner-sphere coordination is usually invoked in the catalytic step, although there is no evidence to support the generality of such a pathway for all ribozymes.Results: We studied the catalytic pathway for the hairpin class of ribozyme. The substitutionally inert transition metal complex cobalt hexaammine [Co(NH₃)₆³⁺) was shown to be as active as Mg²⁺(aq) in promoting hairpin ribozyme activity, demonstrating that inner-sphere pathways are not used by this class of ribozyme. These results were confirmed by studies with R_P- and S_P-phosphorothioate substrate analogs which show a similar reactivity to that of the native substrate towards the magnesium-activated ribozyme. Monovalent cations enhance the activity of Co(NH₃)₆³⁺-promoted reactions, but inhibit Mg²⁺-activated catalysis, demonstrating a requirement for hydrated cations at several key sites in the ribozyme.Conclusions: These results provide clear support for a model of RNA catalysis that does not involve direct coordination of magnesium to the phosphate ester, nor activation of a bound water molecule. A mechanism in which catalysis is carried out by functional groups on the RNA ribozyme itself is possible; such functional groups are likely to have pK_a values that are appropriate for carrying out this catalysis. The metal cofactor would then serve to define the architecture of the catalytic pocket and contribute to the stabilization of transient species, as has been described earlier. Hydrolytic pathways in nucleic acid reactions are apparently more diverse than was previously thought, and the hairpin ribozyme falls into a mechanistically distinct class from the Tetrahymena and the hammerhead ribozymes.     

Elucidating the mechanism for the reduction of nitrite by copper nitrite reductase--a contribution from quantum chemical studies

Density functional methods have been applied to investigate the properties of the active site of copper-containing nitrite reductases and possible reaction mechanisms for the enzyme catalysis. The results for a model of the active site indicate that a hydroxyl intermediate is not formed during the catalytic cycle, but rather a state with a protonated nitrite bound to the reduced copper. Electron affinity calculations indicate that reduction of the T2 copper site does not occur immediately after nitrite binding. Proton affinity calculations are indicative of substantial pK(a) differences between different states of the T2 site. The calculations further suggest that the reaction does not proceed until uptake of a second proton from the bulk solution. They also indicate that Asp-92 may play both a key role as a proton donor to the substrate, and a structural role in promoting catalysis. In the D92N mutant another base, presumably a nearby histidine (His-249) may take the role as the proton donor. On the basis of these model calculations and available experimental evidence, an ordered reaction mechanism for the reduction of nitrite is suggested. An investigation of the binding modes of the nitric oxide product and the nitrite substrate to the model site has also been made, indicating that nitric oxide prefers to bind in an end-on fashion to the reduced T2 site.     

A conformationally restricted guanosine analog reveals the catalytic relevance of three structures of an RNA enzyme

Recent studies indicate that RNA function can be enhanced by the incorporation of conformationally restricted nucleotides. Herein, we use 8-bromoguanosine, a nucleotide analog with an enforced syn conformation, to elucidate the catalytic relevance of ribozyme structures. We chose to study the lead-dependent ribozyme (leadzyme) because structural models derived from NMR, crystal, and computational (MC-Sym) studies differ in which of the three active site guanosines (G7, G9, or G24) have a syn glycosidic torsion angle. Kinetic assays were carried out on 8BrG variants at these three guanosine positions. These data indicate that an 8BrG24 leadzyme is hyperactive, while 8BrG7 and 8BrG9 leadzymes have reduced activity. These findings support the computational model of the leadzyme, rather than the NMR and crystal structures, as being the most relevant to phosphodiester bond cleavage.     

Diels-Alder ribozyme catalysis: a computational approach

A computational comparison of the Diels-Alder reaction of a maleimide and an anthracene in water and the active site of the ribozyme Diels-Alderase is reported. During the course of the catalyzed reaction, the maleimide is held in the hydrophobic pocket while the anthracene approaches to the maleimide through the back passage of the active site. The active site is so narrow that the anthracene has to adopt a tilted approach angle toward maleimide. The conformation of the active site changes marginally at different states of the reaction. Active site dynamics contribution to catalysis has been ruled out. The active site stabilizes the product more than the transition state (TS). The reaction coordinates of the ribozyme reaction in TS, RC1-CD1 and RC4-CD2, are 2.35 and 2.33 A, respectively, compared to 2.37 and 2.36 A in water. The approach angle of anthracene toward maleimide is twisted by 18 degrees in the TS structure of ribozyme reaction while no twisted angle is found in TS of the reaction in water. The free energy barriers for reactions in both ribozyme and water were obtained by umbrella sampling combined with SCCDFTB/MM. The calculated free energy barriers for the ribozyme and water reactions are in good agreement with the experimental values. As expected, Mulliken charges of the atoms involved in the ribozyme reaction change in a similar manner as that of the reaction in water. The proficiency of the Diels-Alder ribozyme reaction originates from the active site holding the two reactants in reactive conformations, in which the reacting atoms are brought together in van der Waals distances and reactants approach to each other at an appropriate angle.     

Coordination environment of a site-bound metal ion in the hammerhead ribozyme determined by 15N and 2H ESEEM spectroscopy

Although site-bound Mg2+ ions have been proposed to influence RNA structure and function, establishing the molecular properties of such sites has been challenging due largely to the unique electrostatic properties of the RNA biopolymer. We have previously determined that, in solution, the hammerhead ribozyme (a self-cleaving RNA) has a high-affinity metal ion binding site characterized by a K(d,app) < 10 microM for Mn2+ in 1 M NaCl and speculated that this site has functional importance in the ribozyme cleavage reaction. Here we determine both the precise location and the hydration level of Mn2+ in this site using ESEEM (electron spin-echo envelope modulation) spectroscopy. Definitive assignment of the high-affinity site to the activity-sensitive A9/G10.1 region is achieved by site-specific labeling of G10.1 with 15N guanine. The coordinated metal ion retains four water ligands as measured by 2H ESEEM spectroscopy. The results presented here show that a functionally important, specific metal binding site is uniquely populated in the hammerhead ribozyme even in a background of high ionic strength. Although it has a relatively high thermodynamic affinity, this ion remains partially hydrated and is chelated to the RNA by just two ligands.     

RNA enzymes with two small-molecule substrates

Background: The ‘RNA world’ hypothesis posits ancient organisms employing versatile catalysis by RNAs. In particular, such a metabolism would have required RNA catalysts that join small molecules. Such anabolic reactions now occur very widely, for example in phospholipid, terpene, amino acid and nucleotide synthetic pathways in modern organisms. Present RNA systems, however, do not perform such reactions using substrates that do not base pair. Here we ask whether this lack is a methodological artifact due to the practice of selection-amplification, or a fundamental property of active sites reconstructed within RNA structures.Results: Three rationally modified RNA enzymes, Iso6-G, Iso6-2G and Iso6-3G, catalyze the formation of (5′→5′) polyphosphate-linked oligonucleotides in trans. One of these, Iso6-G RNA, has a specific substrate site for a guanosine triphosphate, GTP, dGTP or ddGTP, and one nonspecific substrate site for a terminal-phosphate-containing small molecule. This ribozyme catalyzes multiple turnovers, proceeding at a constant rate. Guanosine specificity is probably not attributable to Watson-Crick base pairing.Conclusions: Ribozymes can readily bind multiple small-molecule substrates simultaneously and catalyze reactions that build up larger products, apparently independent of substrate-RNA Watson-Crick base pairing. RNA enzymes therefore parallel proteins, which often overcome the entropic difficulties of positioning multiple small substrates for catalysis of anabolic reactions. These results support the idea of a complex ancestral metabolism based on RNA catalysis.     

General base catalysis for cleavage by the active-site cytosine of the hepatitis delta virus ribozyme: QM/MM calculations establish chemical feasibility

The hepatitis delta virus (HDV) ribozyme is an RNA motif embedded in human pathogenic HDV RNA. Previous experimental studies have established that the active-site nucleotide C75 is essential for self-cleavage of the ribozyme, although its exact catalytic role in the process remains debated. Structural data from X-ray crystallography generally indicate that C75 acts as the general base that initiates catalysis by deprotonating the 2'-OH nucleophile at the cleavage site, while a hydrated magnesium ion likely protonates the 5'-oxygen leaving group. In contrast, some mechanistic studies support the role of C75 acting as general acid and thus being protonated before the reaction. We report combined quantum chemical/molecular mechanical calculations for the C75 general base pathway, utilizing the available structural data for the wild type HDV genomic ribozyme as a starting point. Several starting configurations differing in magnesium ion placement were considered and both one-dimensional and two-dimensional potential energy surface scans were used to explore plausible reaction paths. Our calculations show that C75 is readily capable of acting as the general base, in concert with the hydrated magnesium ion as the general acid. We identify a most likely position for the magnesium ion, which also suggests it acts as a Lewis acid. The calculated energy barrier of the proposed mechanism, approximately 20 kcal/mol, would lower the reaction barrier by approximately 15 kcal/mol compared with the uncatalyzed reaction and is in good agreement with experimental data.     

The glmS ribozyme tunes the catalytically critical pK(a) of its coenzyme glucosamine-6-phosphate

The glmS ribozyme riboswitch is the first known natural catalytic RNA that employs a small-molecule cofactor. Binding of glucosamine-6-phosphate (GlcN6P) uncovers the latent self-cleavage activity of the RNA, which adopts a catalytically competent conformation that is nonetheless inactive in the absence of GlcN6P. Structural and analogue studies suggest that the amine of GlcN6P functions as a general acid-base catalyst, while its phosphate is important for binding affinity. However, the solution pK(a) of the amine is 8.06 ± 0.05, which is not optimal for proton transfer. Here we used Raman crystallography directly to determine the pK(a)'s of GlcN6P bound to the glmS ribozyme. Binding to the RNA lowers the pK(a) of the amine of GlcN6P to 7.26 ± 0.09 and raises the pK(a) of its phosphate to 6.35 ± 0.09. Remarkably, the pK(a)'s of these two functional groups are unchanged from their values for free GlcN6P (8.06 ± 0.05 and 5.98 ± 0.05, respectively) when GlcN6P binds to the catalytically inactive but structurally unperturbed G40A mutant of the ribozyme, thus implicating the ribozyme active site guanine in pK(a) tuning. This is the first demonstration that a ribozyme can tune the pK(a) of a small-molecule ligand. Moreover, the anionic glmS ribozyme in effect stabilizes the neutral amine of GlcN6P by lowering its pK(a). This is unprecedented and illustrates the chemical sophistication of ribozyme active sites.     

A hydrogen-bonding triad stabilizes the chemical transition state of a group I ribozyme

BACKGROUND: The group I intron is an RNA enzyme capable of efficiently catalyzing phosphoryl-transfer reactions. Functional groups that stabilize the chemical transition state of the cleavage reaction have been identified, but they are all located within either the 5'-exon (P1) helix or the guanosine cofactor, which are the substrates of the reaction. Functional groups within the ribozyme active site are also expected to assist in transition-state stabilization, and their role must be explored to understand the chemical basis of group I intron catalysis. RESULTS: Using nucleotide analog interference mapping and site-specific functional group substitution experiments, we demonstrate that the 2'-OH at A207, a highly conserved nucleotide in the ribozyme active site, specifically stabilizes the chemical transition state by approximately 2 kcal mol-1. The A207 2'-OH only makes its contribution when the U(-1) 2'-OH immediately adjacent to the scissile phosphate is present, suggesting that the 2'-OHs of A207 and U(-1) interact during the chemical step. CONCLUSIONS: These data support a model in which the 3'-oxyanion leaving group of the transesterification reaction is stabilized by a hydrogen-bonding triad consisting of the 2'-OH groups of U(-1) and A207 and the exocyclic amine of G22. Because all three nucleotides occur within highly conserved non-canonical base pairings, this stabilization mechanism is likely to occur throughout group I introns. Although this mechanism utilizes functional groups distinctive of RNA enzymes, it is analogous to the transition states of some protein enzymes that perform similar phosphoryl-transfer reactions.     

Nucleobase participation in ribozyme catalysis

We constructed a modified form of the VS ribozyme containing an imidazole ring in place of adenine at position 756. The novel ribozyme is active in both cleavage and ligation reactions. The reaction is efficient, although relatively slow. The results are consistent with a role for nucleobase catalysis in the catalytic mechanism of this ribozyme.