Active Site Labeling of G8 in the Hairpin Ribozyme: Implications for Structure and Mechanism

Who stole the proton? Suspect general base guanine found with a smoking gun in the pistol ribozyme

The pistol ribozyme (Psr) is one among the most recently discovered classes of small nucleolytic ribozymes that catalyze site-specific RNA self-cleavage through 2'-O-transphosphorylation. The Psr contains a conserved guanine (G40) that in crystal structures is in a position suggesting it plays the role of the general base to abstract a proton from the nucleophile to activate the reaction. Although some functional data is consistent with this mechanistic role, a notable exception is 2-aminopurine (2AP) substitution which has no effect on the rate, unlike similar substitutions across other so-called "G + M" and "G + A" ribozyme classes. Herein we postulate that an alternate conserved guanine, G42, is the primary general base, and provide evidence from molecular simulations that the active site of Psr can undergo local refolding into a structure that is consistent with the common "L-platform/L-scaffold" architecture identified in G + M and G + A ribozyme classes with Psr currently the notable exception. We summarize the key currently available experimental data and present new classical and combined quantum mechanical/molecular mechanical simulation results that collectively suggest a new hypothesis. We hypothesize that there are two available catalytic pathways supported by different conformational states connected by a local refolding of the active site: (1) a primary pathway with an active site architecture aligned with the L-platform/L-scaffold framework where G42 acts as a general base, and (2) a secondary pathway with the crystallographic active site architecture where G40 acts as a general base. We go on to make several experimentally testable predictions, and suggest specific experiments that would ultimately bring closure to the mystery as to "who stole the proton in the pistol ribozyme?".

Deciphering the Self-Cleavage Reaction Mechanism of Hairpin Ribozyme

Hairpin ribozyme catalyzes the reversible self-cleavage of phosphodiester bonds which plays prominent roles in key biological processes involving RNAs. Despite impressive advances on ribozymatic self-cleavage, critical aspects of its molecular reaction mechanism remain controversially debated. Here, we generate and analyze the multidimensional free energy landscape that underlies the reaction using extensive QM/MM metadynamics simulations to investigate in detail the full self-cleavage mechanism. This allows us to answer several pertinent yet controversial questions concerning activation of the 2'-OH group, the mechanistic role of water molecules present in the active site, and the full reaction pathway including the structures of transition states and intermediates. Importantly, we find that a sufficiently unrestricted reaction subspace must be mapped using accelerated sampling methods in order to compute the underlying free energy landscape. It is shown that lower-dimensional sampling where the bond formation and cleavage steps are coupled does not allow the system to sufficiently explore the landscape. On the basis of a three-dimensional free energy surface spanned by flexible generalized coordinates, we find that 2'-OH is indirectly activated by adjacent G8 nucleobase in conjunction with stabilizing H-bonding involving water. This allows the proton of the 2'-OH group to directly migrate toward the 5'-leaving group via a nonbridging oxygen of the phosphodiester link. At variance with similar enzymatic processes where water wires connected to protonable side chains of the protein matrix act as transient proton shuttles, no such de/reprotonation events of water molecules are found to be involved in this ribozymatic transesterification. Overall, our results support an acid-catalyzed reaction mechanism where A38 nucleobase directly acts as an acid whereas G8, in stark contrast, participates only indirectly via stabilizing the nascent nucleophile for subsequent attack. Moreover, we conclude that self-cleavage of hairpin ribozyme follows an A_N + D_N two-step associative pathway where the rate-determining step is the cleavage of the phosphodiester bond. These results provide a major advancement in our understanding of the unique catalytic mechanism of hairpin ribozyme which will fruitfully impact on the design of synthetic ribozymes.

Electrophilic activity-based RNA probes reveal a self-alkylating RNA for RNA labeling

Probes that form covalent bonds with RNA molecules on the basis of their chemical reactivity would advance our ability to study the transcriptome. We developed a set of electrophilic activity-based RNA probes designed to react with unusually nucleophilic RNAs. We used these probes to identify reactive genome-encoded RNAs, resulting in the discovery of a 42-nt catalytic RNA from an archaebacterium that reacts with a 2,3-disubstituted epoxide at N7 of a specific guanosine. Detailed characterization of the catalytic RNA revealed the structural requirements for reactivity. We developed this catalytic RNA into a general tool to selectively conjugate a small molecule to an RNA of interest. This strategy enabled up to 500-fold enrichment of target RNA from total mammalian RNA or from cell lysate. We demonstrated the utility of this approach by selectively capturing proteins in yeast cell lysate that bind the ASH1 mRNA.

Evidence for the role of active site residues in the hairpin ribozyme from molecular simulations along the reaction path

The hairpin ribozyme accelerates a phosphoryl transfer reaction without catalytic participation of divalent metal ions. Residues A38 and G8 have been implicated as playing roles in general acid and base catalysis, respectively. Here we explore the structure and dynamics of key active site residues using more than 1 μs of molecular dynamics simulations of the hairpin ribozyme at different stages along the catalytic pathway. Analysis of results indicates hydrogen bond interactions between the nucleophile and proR nonbridging oxygen are correlated with active inline attack conformations. Further, the simulation results suggest a possible alternative role for G8 to promote inline fitness and facilitate activation of the nucleophile by hydrogen bonding, although this does not necessarily exclude an additional role as a general base. Finally, we suggest that substitution of G8 with N7- or N3-deazaguanosine which have elevated pK_a values, both with and without thio modifications at the 5′ leaving group position, would provide valuable insight into the specific role of G8 in catalysis.

Cobalt(III)hexaammine-dependent photocrosslinks in the hairpin ribozyme

We have utilized the hairpin ribozyme, an RNA enzyme whose structure has been solved by high-resolution methods, to develop a new tool for mapping nucleobase-stacking interactions and potential metal-binding sites in RNA molecules. This tool involves the photoactivation of a specifically bound cobalt(III)hexaammine molecule at wavelengths corresponding to excitation of the metal ion complex only; no base excitation is involved. The photoexcitation initiates a process which strongly promotes the formation of a novel covalent bond or crosslink between one base (termed the "first base"), which is close in space to the excited cobalt(III)hexaammine complex, and another base upon which the first base is closely stacked. These crosslinked species can be isolated and sequenced; their activities can be analyzed to ensure that the crosslinked structures represent an active conformation of the molecule. We have shown that, as in electron transfer in DNA, several criteria must be met to result in the successful formation of these crosslinks. These include the appropriate oxidation potential of the first donor base, the stacking and close interaction of the two donor bases involved in the crosslink, and the binding of a specific cobalt(III)hexaammine molecule to the first donor base. Additionally, we have determined that this crosslinking is pH-sensitive, although the cause of this sensitivity remains unknown. This tool has proven useful in the past for the analysis of the hairpin ribozyme folded structure, and has been applied to identify potential metal-binding sites on the hairpin and extended hammerhead ribozymes.

Do the hairpin and VS ribozymes share a common catalytic mechanism based on general acid-base catalysis? A critical assessment of available experimental data

The active centers of the hairpin and VS ribozymes are both generated by the interaction of two internal loops, and both ribozymes use guanine and adenine nucleobases to accelerate cleavage and ligation reactions. The centers are topologically equivalent and the relative positioning of key elements the same. There is good evidence that the cleavage reaction of the VS ribozyme is catalyzed by the guanine (G638) acting as general base and the adenine (A756) as general acid. We now critically evaluate the experimental mechanistic evidence for the hairpin ribozyme. We conclude that all the available data are fully consistent with a major contribution to catalysis by general acid-base catalysis involving the adenine (A38) and guanine (G8). It appears that the two ribozymes are mechanistically equivalent.

Investigation of the catalytic mechanism of a synthetic DNAzyme with protein-like functionality: an RNaseA mimic?

The protein enzyme ribonuclease A (RNaseA) cleaves RNA with catalytic perfection, although with little sequence specificity, by a divalent metal ion (M(2+))-independent mechanism in which a pair of imidazoles provides general acid and base catalysis, while a cationic amine provides electrostatic stabilization of the transition state. Synthetic imitation of this remarkable organo-catalyst ("RNaseA mimicry") has been a longstanding goal in biomimetic chemistry. The 9(25)-11 DNAzyme contains synthetically modified nucleotides presenting both imidazole and cationic amine side chains, and catalyzes RNA cleavage with turnover in the absence of M(2+) similarly to RNaseA. Nevertheless, the catalytic roles, if any, of the "protein-like" functional groups have not been defined, and hence the question remains whether 9(25)-11 engages any of these functionalities to mimic aspects of the mechanism of RNaseA. To address this question, we report a mechanistic investigation of 9(25)-11 catalysis wherein we have employed a variety of experiments, such as DNAzyme functional group deletion, mechanism-based affinity labeling, and bridging and nonbridging phosphorothioate substitution of the scissile phosphate. Several striking parallels exist between the results presented here for 9(25)-11 and the results of analogous experiments applied previously to RNaseA. Specifically, our results implicate two particular imidazoles in general acid and base catalysis and suggest that a specific cationic amine stabilizes the transition state via diastereoselective interaction with the scissile phosphate. Overall, 9(25)-11 appears to meet the minimal criteria of an RNaseA mimic; this demonstrates how added synthetic functionality can expand the mechanistic repertoire available to a synthetic DNA-based catalyst.

Molecular dynamics suggest multifunctionality of an adenine imino group in acid-base catalysis of the hairpin ribozyme

Despite numerous structural and biochemical investigations, the catalytic mechanism of hairpin ribozyme self-cleavage remains elusive. To gain insight into the coupling of active site dynamics with activity of this small catalytic RNA, we analyzed a total of approximately 300 ns of molecular dynamics (MD) simulations. Our simulations predict improved global stability for an in vitro selected "gain of function" mutation, which is validated by native gel electrophoretic mobility shift assay. We observe that active site nucleobases and water molecules stabilize a geometry favorable to catalysis through a dynamic hydrogen bonding network. Simulations in which A38 is unprotonated show its N1 move into close proximity of the active site 2'-OH, indicating that A38 may act as a general base during cleavage, a role that has generally been discounted due to the longer distances observed in crystal structures involving inactivating substrate analogs. By contrast, simulations in which N1 of A38 is protonated place N1 in close proximity to the 5'-oxygen leaving group, which supports the proposal that A38 serves as a general acid. In analogy to protein enzymes, we discuss a plausible mechanism in which A38 acts bifunctionally and shuttles a proton directly from the 2'-OH to the 5'-oxygen. Furthermore, our simulations suggest an important role for protonation of N1 of A38 in promoting a favorable geometry similar to that observed in transition-state analog crystal structures, and support previously proposed roles of A38, G8, and long residency water molecules in transition-state stabilization.

A self-cleaving DNA enzyme modified with amines, guanidines and imidazoles operates independently of divalent metal cations (M2+)

The selection of modified DNAzymes represents an important endeavor in expanding the chemical and catalytic properties of catalytic nucleic acids. Few examples of such exist and to date, there is no example where three different modified bases have been simultaneously incorporated for catalytic activity. Herein, dCTP, dATP and dUTP bearing, respectively, a cationic amine, an imidazole and a cationic guanidine, were enzymatically polymerized on a DNA template for the selection of a highly functionalized DNAzyme, called DNAzyme 9-86, that catalyzed (M(2+))-independent self-cleavage under physiological conditions at a single ribo(cytosine)phosphodiester linkage with a rate constant of (0.134 +/- 0.026) min(-1). A pH rate profile analysis revealed pK(a)'s of 7.4 and 8.1, consistent with both general acid and base catalysis. The presence of guanidinium cations permits cleavage at significantly higher temperatures than previously observed for DNAzymes with only amines and imidazoles. Qualitatively, DNAzyme 9-86 presents an unprecedented ensemble of synthetic functionalities while quantitatively it expresses one of the highest reported values for any self-cleaving nucleic acid when investigated under M(2+)-free conditions at 37 degrees C.

Synthesis of a cross-linked branched polymer network in the interior of a protein cage

A goal of biomimetic chemistry is to use the hierarchical architecture inherent in biological systems to guide the synthesis of functional three-dimensional structures. Viruses and other highly symmetrical protein cage architectures provide defined scaffolds to initiate hierarchical structure assembly. Here we demonstrate that a cross-linked branched polymer can be initiated and synthesized within the interior cavity of a protein cage architecture. Creating this polymer network allows for the spatial control of pendant reactive sites and dramatically increases the stability of the cage architecture. This material was generated by the sequential coupling of multifunctional monomers using click chemistry to create a branched cross-linked polymer network. Analysis of polymer growth by mass spectrometry demonstrated that the polymer was initiated at the interior surface of the cage at genetically introduced cysteine reactive sites. The polymer grew as expected to generation 2.5 where it was limited by the size constraints of the cavity. The polymer network was fully cross-linked across protein subunits that make up the cage and extended the thermal stability for the cage to at least 120 degrees C. The introduced reactive centers were shown to be active and their number density increased with increasing generation. This synthetic approach provides a new avenue for creating defined polymer networks, spatially constrained by a biological template.

A rugged free energy landscape separates multiple functional RNA folds throughout denaturation

The dynamic mechanisms by which RNAs acquire biologically functional structures are of increasing importance to the rapidly expanding fields of RNA therapeutics and biotechnology. Large energy barriers separating misfolded and functional states arising from alternate base pairing are a well-appreciated characteristic of RNA. In contrast, it is typically assumed that functionally folded RNA occupies a single native basin of attraction that is free of deeply dividing energy barriers (ergodic hypothesis). This assumption is widely used as an implicit basis to interpret experimental ensemble-averaged data. Here, we develop an experimental approach to isolate persistent sub-populations of a small RNA enzyme and show by single molecule fluorescence resonance energy transfer (smFRET), biochemical probing and high-resolution mass spectrometry that commitment to one of several catalytically active folds occurs unexpectedly high on the RNA folding energy landscape, resulting in partially irreversible folding. Our experiments reveal the retention of molecular heterogeneity following the complete loss of all native secondary and tertiary structure. Our results demonstrate a surprising longevity of molecular heterogeneity and advance our current understanding beyond that of non-functional misfolds of RNA kinetically trapped on a rugged folding-free energy landscape.

Structural effects of nucleobase variations at key active site residue Ade38 in the hairpin ribozyme

The hairpin ribozyme requires functional groups from Ade38 to achieve efficient bond cleavage or ligation. To identify molecular features that contribute to catalysis, structures of position 38 base variants 2,6-diaminopurine (DAP), 2-aminopurine (AP), cytosine (Cyt), and guanine (Gua) were determined between 2.2 and 2.8 A resolution. For each variant, two substrate modifications were compared: (1) a 2'-O-methyl-substituent at Ade-1 was used in lieu of the nucleophile to mimic the precatalytic state, and (2) a 3'-deoxy-2',5'-phosphodiester linkage between Ade-1 and Gua+1 was used to mimic a reaction-intermediate conformation. While the global fold of each variant remained intact, the results revealed the importance of Ade38 N1 and N6 groups. Absence of N6 resulting from AP38 coincided with failure to localize the precatalytic scissile phosphate. Cyt38 severely impaired catalysis in a prior study, and its structures here indicated an anti base conformation that sequesters the imino moiety from the scissile bond. Gua38 was shown to be even more deleterious to activity. Although the precatalytic structure was nominally affected, the reaction-intermediate conformation indicated a severe electrostatic clash between the Gua38 keto oxygen and the pro-Rp oxygen of the scissile bond. Overall, position 38 modifications solved in the presence of 2'-OMe Ade-1 deviated from in-line geometry, whereas variants with a 2',5' linkage exhibited S-turn destabilization, as well as base conformational changes from syn to anti. These findings demonstrate the importance of the Ade38 Watson-Crick face in attaining a reaction-intermediate state and the sensitivity of the RNA fold to restructuring when electrostatic and shape features fail to complement.

An important role of G638 in the cis-cleavage reaction of the Neurospora VS ribozyme revealed by a novel nucleotide analog incorporation method

We describe a chemical coupling procedure that allows joining of two RNAs, one of which contains a site-specific base analog substitution, in the absence of divalent ions. This method allows incorporation of nucleotide analogs at specific positions even into large, cis-cleaving ribozymes. Using this method we have studied the effects of substitution of G638 in the cleavage site loop of the VS ribozyme with a variety of purine analogs having different functional groups and pK(a) values. Cleavage rate versus pH profiles combined with kinetic solvent isotope experiments indicate an important role for G638 in proton transfer during the rate-limiting step of the cis-cleavage reaction.

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.

A comparison of vanadate to a 2'-5' linkage at the active site of a small ribozyme suggests a role for water in transition-state stabilization

The potential for water to participate in RNA catalyzed reactions has been the topic of several recent studies. Here, we report crystals of a minimal, hinged hairpin ribozyme in complex with the transition-state analog vanadate at 2.05 A resolution. Waters are present in the active site and are discussed in light of existing views of catalytic strategies employed by the hairpin ribozyme. A second structure harboring a 2',5'-phosphodiester linkage at the site of cleavage was also solved at 2.35 A resolution and corroborates the assignment of active site waters in the structure containing vanadate. A comparison of the two structures reveals that the 2',5' structure adopts a conformation that resembles the reaction intermediate in terms of (1) the positioning of its nonbridging oxygens and (2) the covalent attachment of the 2'-O nucleophile with the scissile G+1 phosphorus. The 2',5'-linked structure was then overlaid with scissile bonds of other small ribozymes including the glmS metabolite-sensing riboswitch and the hammerhead ribozyme, and suggests the potential of the 2',5' linkage to elicit a reaction-intermediate conformation without the need to form metalloenzyme complexes. The hairpin ribozyme structures presented here also suggest how water molecules bound at each of the nonbridging oxygens of G+1 may electrostatically stabilize the transition state in a manner that supplements nucleobase functional groups. Such coordination has not been reported for small ribozymes, but is consistent with the structures of protein enzymes. Overall, this work establishes significant parallels between the RNA and protein enzyme worlds.