Functional involvement of G8 in the hairpin ribozyme cleavage 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?".

The role of Na+ in catalysis by the 8-17 DNAzyme

The 8-17 DNAzyme is the most studied deoxyribozyme in terms of its molecular mechanism; hence it has become a model system to understand the basis behind DNA catalysis. New functional studies and the recent attainment of high-resolution X-ray structures, in addition to theoretical calculations have offered a great opportunity to gain a broader comprehension of its mechanism; however many aspects are unclear yet, especially regarding the precise role of metal ions in catalysis. Recently, molecular dynamics simulations have suggested for the first time a specific and dynamical participation of Na⁺ in the mechanism through the reaction pathway, besides the roles proposed for divalent metal cofactors. Herein, we present experimental evidence of a cooperative role of the monovalent cation Na⁺ in catalysis that is in line with these theoretical suggestions. Our findings show a clear influence of the concentration of Na⁺ on the activity of the 8-17 DNAzyme when Pb²⁺ is used as the cofactor. Interestingly, this effect is not noticed with Mg²⁺, indicating a particular contribution of the monovalent ion to catalysis that would operate preferentially with Pb²⁺. We have also found that Na⁺ affects the pK_a of the general base and the general acid, indicating its influence on general acid-base catalysis, already identified as part of the mechanism of the 8-17 DNAzyme. Finally, our results emphasize the need to consider Na⁺ carefully in the design and analysis of functional studies of catalytic DNAs and its possible specific role in their mechanisms.

Structural Insights Into Tautomeric Dynamics in Nucleic Acids and in Antiviral Nucleoside Analogs

DNA (2'-deoxyribonucleic acid) and RNA (ribonucleic acid) play diverse functional roles in biology and disease. Despite being comprised primarily of only four cognate nucleobases, nucleic acids can adopt complex three-dimensional structures, and RNA in particular, can catalyze biochemical reactions to regulate a wide variety of biological processes. Such chemical versatility is due in part to the phenomenon of nucleobase tautomerism, whereby the bases can adopt multiple, yet distinct isomeric forms, known as tautomers. For nucleobases, tautomers refer to structural isomers that differ from one another by the position of protons. By altering the position of protons on nucleobases, many of which play critical roles for hydrogen bonding and base pairing interactions, tautomerism has profound effects on the biochemical processes involving nucleic acids. For example, the transient formation of minor tautomers during replication could generate spontaneous mutations. These mutations could arise from the stabilization of mismatches, in the active site of polymerases, in conformations involving minor tautomers that are indistinguishable from canonical base pairs. In this review, we discuss the evidence for tautomerism in DNA, and its consequences to the fidelity of DNA replication. Also reviewed are RNA systems, such as the riboswitches and self-cleaving ribozymes, in which tautomerism plays a functional role in ligand recognition and catalysis, respectively. We also discuss tautomeric nucleoside analogs that are efficacious as antiviral drug candidates such as molnupiravir for coronaviruses and KP1212 for HIV. The antiviral efficacy of these analogs is due, in part, to their ability to exist in multiple tautomeric forms and induce mutations in the replicating viral genomes. From a technical standpoint, minor tautomers of nucleobases are challenging to identify directly because they are rare and interconvert on a fast, millisecond to nanosecond, time scale. Nevertheless, many approaches including biochemical, structural, computational and spectroscopic methods have been developed to study tautomeric dynamics in RNA and DNA systems, and in antiviral nucleoside analogs. An overview of these methods and their applications is included here.

Investigation of the pKa of the Nucleophilic O2' of the Hairpin Ribozyme

Small ribozymes cleave their RNA phosphodiester backbone by catalyzing a transphosphorylation reaction wherein a specific O2' functions as the nucleophile. While deprotonation of this alcohol through its acidification would increase its nucleophilicity, little is known about the pK_a of this O2' in small ribozymes, in part because high pK_a's are not readily accessible experimentally. Herein, we turn to molecular dynamics to calculate the pK_a of the nucleophilic O2' in the hairpin ribozyme and to study interactions within the active site that may impact its value. We estimate the pK_a of the nucleophilic O2' in the wild-type hairpin ribozyme to be 18.5 ± 0.8, which is higher than the reference compound, and identify a correlation between proper positioning of the O2' for nucleophilic attack and elevation of its pK_a. We find that monovalent ions may play a role in depression of the O2' pK_a, while the exocyclic amine appears to be important for organizing the ribozyme active site. Overall, this study suggests that the pK_a of the O2' is raised in the ground state and lowers during the course of the reaction owing to positioning and metal ion interactions.

Identification of over 200-fold more hairpin ribozymes than previously known in diverse circular RNAs

Self-cleaving ribozymes are catalytic RNAs that cut themselves at a specific inter-nucleotide linkage. They serve as a model of RNA catalysis, and as an important tool in biotechnology. For most of the nine known structural classes of self-cleaving ribozymes, at least hundreds of examples are known, and some are present in multiple domains of life. By contrast, only four unique examples of the hairpin ribozyme class are known, despite its discovery in 1986. We bioinformatically predicted 941 unique hairpin ribozymes of a different permuted form from the four previously known hairpin ribozymes, and experimentally confirmed several diverse predictions. These results profoundly expand the number of natural hairpin ribozymes, enabling biochemical analysis based on natural sequences, and suggest that a distinct permuted form is more biologically relevant. Moreover, all novel hairpins were discovered in metatranscriptomes. They apparently reside in RNA molecules that vary both in size—from 381 to 5170 nucleotides—and in protein content. The RNA molecules likely replicate as circular single-stranded RNAs, and potentially provide a dramatic increase in diversity of such RNAs. Moreover, these organisms have eluded previous attempts to isolate RNA viruses from metatranscriptomes—suggesting a significant untapped universe of viruses or other organisms hidden within metatranscriptome sequences.

Hydrated metal ion as a general acid in the catalytic mechanism of the 8-17 DNAzyme

The RNA-cleaving 8-17 DNAzyme, which is a metalloenzyme that depends on divalent metal ions for its function, is the most studied catalytic DNA in terms of its mechanism. By the end of 2017, a report of the crystal structure of the enzyme-substrate complex in the presence of Pb²⁺ probed some of the previous findings and opened new questions, especially around the participation of the metal ion in the catalytic mechanism and the promiscuity exhibited by the enzyme in terms of the metal cofactor required for catalysis. In this article we explore the role of the divalent metal ion in the mechanism of the 8-17 DNAzyme as a general acid, by measuring the influence of pH over the activity of a slower variant of the enzyme in the presence of Pb²⁺. We replaced G14, which has been identified as a general base in the mechanism of the enzyme, by the unnatural analog 2-aminopurine, with a lower pK_a value of the N1 group. With this approach, we obtained a bell-shaped pH-rate profile with experimental pK_a values of 5.4 and 7.0. Comparing these results with previous pH-rate profiles in the presence of Mg²⁺, our findings suggest the stabilization of the 5'-O leaving group by the hydrated metal ion acting as a general acid, in addition to the activation of the 2'-OH nucleophile by the general base G14.

Insights into DNA catalysis from structural and functional studies of the 8-17 DNAzyme

The review examines functional knowledge gathered over two decades of research on the 8-17 DNAzyme, focusing on three aspects: the structural requirements for catalysis, the role of metal ions and the participation of general acid-base catalysis.

Novel ribozymes: discovery, catalytic mechanisms, and the quest to understand biological function

Small endonucleolytic ribozymes promote the self-cleavage of their own phosphodiester backbone at a specific linkage. The structures of and the reactions catalysed by members of individual families have been studied in great detail in the past decades. In recent years, bioinformatics studies have uncovered a considerable number of new examples of known catalytic RNA motifs. Importantly, entirely novel ribozyme classes were also discovered, for most of which both structural and biochemical information became rapidly available. However, for the majority of the new ribozymes, which are found in the genomes of a variety of species, a biological function remains elusive. Here, we concentrate on the different approaches to find catalytic RNA motifs in sequence databases. We summarize the emerging principles of RNA catalysis as observed for small endonucleolytic ribozymes. Finally, we address the biological functions of those ribozymes, where relevant information is available and common themes on their cellular activities are emerging. We conclude by speculating on the possibility that the identification and characterization of proteins that we hypothesize to be endogenously associated with catalytic RNA might help in answering the ever-present question of the biological function of the growing number of genomically encoded, small endonucleolytic ribozymes.

Comparison of the Structures and Mechanisms of the Pistol and Hammerhead Ribozymes

Comparison of the secondary and three-dimensional structures of the hammerhead and pistol ribozymes reveals many close similarities, so in this work we have asked if they are mechanistically identical. We have determined a new crystal structure of the pistol ribozyme and have shown that G40 acts as general base in the cleavage reaction. The conformation in the active site ensures an in-line attack of the O2' nucleophile, and the conformation at the scissile phosphate and the position of the general base are closely similar to those in the hammerhead ribozyme. However, the two ribozymes differ in the nature of the general acid. 2'-Amino substitution experiments indicate that the general acid of the hammerhead ribozyme is the O2' of G8, while that of the pistol ribozyme is a hydrated metal ion. The two ribozymes are related but mechanistically distinct.

Gibbs/MCMC Sampling for Multiple RNA Interaction with Sub-optimal Solutions

Multiple RNA interaction can be modeled as a problem in combinatorial optimization, where the "optimal" structure is driven by an energy-minimization-like algorithm. However, the actual structure may not be optimal in this computational sense. Moreover, it is not necessarily unique. Therefore, alternative sub-optimal solutions are needed to cover the biological ground. We present a combinatorial formulation for the Multiple RNA Interaction problem with approximation algorithms to handle various interaction patterns, which when combined with Gibbs sampling and MCMC (Markov Chain Monte Carlo), can efficiently generate a reasonable number of optimal and sub-optimal solutions. When viable structures are far from an optimal solution, exploring dependence among different parts of the interaction can increase their score and boost their candidacy for the sampling algorithm. By clustering the solutions, we identify few representatives that are distinct enough to suggest possible alternative structures.

Evidence of a General Acid-Base Catalysis Mechanism in the 8-17 DNAzyme

DNAzymes are catalytic DNA molecules that can perform a variety of reactions. Although advances have been made in obtaining DNAzymes via in vitro selection and many of them have been developed into sensors and imaging agents for metal ions, bacteria, and other molecules, the structural features responsible for these enzymatic reactions are still not well understood. Previous studies of the 8-17 DNAzyme have suggested conserved guanines close to the phosphodiester transfer site may play a role in the catalytic reaction. To identify the specific guanine and functional group of the guanine responsible for the reaction, we herein report the effects of replacing G1.1 and G14 (G; p Ka,N1 = 9.4) with analogues with a different p Ka at the N1 position, such as inosine (G14I; p Ka,N1 = 8.7), 2,6-diaminopurine (G14diAP; p Ka,N1 = 5.6), and 2-aminopurine (G14AP; p Ka,N1 = 3.8) on pH-dependent reaction rates. A comparison of the pH dependence of the reaction rates of these DNAzymes demonstrated that G14 in the bulge loop next to the cleavage site, is involved in proton transfer at the catalytic site. In contrast, we did not find any evidence of G1.1 being involved in acid-base catalysis. These results support general acid-base catalysis as a feasible strategy used in DNA catalysis, as in RNA and protein enzymes.

Elucidation of Catalytic Strategies of Small Nucleolytic Ribozymes From Comparative Analysis of Active Sites

A number of small, self-cleaving ribozyme classes have been identified including the hammerhead, hairpin, hepatitis delta virus (HDV), Varkud satellite (VS), glmS, twister, hatchet, pistol, and twister sister ribozymes. Within the active sites of these ribozymes, myriad functional groups contribute to catalysis. There has been extensive structure-function analysis of individual ribozymes, but the extent to which catalytic devices are shared across different ribozyme classes is unclear. As such, emergent catalytic principles for ribozymes may await discovery. Identification of conserved catalytic devices can deepen our understanding of RNA catalysis specifically and of enzymic catalysis generally. To probe similarities and differences amongst ribozyme classes, active sites from more than 80 high-resolution crystal structures of self-cleaving ribozymes were compared computationally. We identify commonalities amongst ribozyme classes pertaining to four classic catalytic devices: deprotonation of the 2'OH nucleophile (γ), neutralization of the non-bridging oxygens of the scissile phosphate (β), neutralization of the O5' leaving group (δ), and in-line nucleophilic attack (α). In addition, we uncover conservation of two catalytic devices, each of which centers on the activation of the 2'OH nucleophile by a guanine: one to acidify the 2'OH by hydrogen bond donation to it (γ') and one to acidify the 2'OH by releasing it from non-productive interactions by competitive hydrogen bonding (γ''). Our findings reveal that the amidine functionalities of G, A, and C are especially important for these strategies, and help explain absence of U at ribozyme active sites. The identified γ' and γ'' catalytic strategies help unify the catalytic strategies shared amongst catalytic RNAs and may be important for large ribozymes, as well as protein enzymes that act on nucleic acids.

How RNA acts as a nuclease: some mechanistic comparisons in the nucleolytic ribozymes

Recent structural and mechanistic studies have shed considerable light on the catalytic mechanisms of nucleolytic ribozymes. The discovery of several new ribozymes in this class has now allowed comparisons to be made, and the beginnings of mechanistic groupings to emerge.

The structure of a nucleolytic ribozyme that employs a catalytic metal ion

The TS ribozyme (originally called "twister sister") is a catalytic RNA. We present a crystal structure of the ribozyme in a pre-reactive conformation. Two co-axial helical stacks are organized by a three-way junction and two tertiary contacts. Five divalent metal ions are directly coordinated to RNA ligands, making important contributions to the RNA architecture. The scissile phosphate lies in a quasihelical loop region that is organized by a network of hydrogen bonding. A divalent metal ion is directly bound to the nucleobase 5' to the scissile phosphate, with an inner-sphere water molecule positioned to interact with the O2' nucleophile. The rate of ribozyme cleavage correlated in a log-linear manner with divalent metal ion pKa, consistent with proton transfer in the transition state, and we propose that the bound metal ion is a likely general base for the cleavage reaction. Our data indicate that the TS ribozyme functions predominantly as a metalloenzyme.

Structural and Biochemical Properties of Novel Self-Cleaving Ribozymes

Fourteen well-defined ribozyme classes have been identified to date, among which nine are site-specific self-cleaving ribozymes. Very recently, small self-cleaving ribozymes have attracted renewed interest in their structure, biochemistry, and biological function since the discovery, during the last three years, of four novel ribozymes, termed twister, twister sister, pistol, and hatchet. In this review, we mainly address the structure, biochemistry, and catalytic mechanism of the novel ribozymes. They are characterized by distinct active site architectures and divergent, but similar, biochemical properties. The cleavage activities of the ribozymes are highly dependent upon divalent cations, pH, and base-specific mutations, which can cause changes in the nucleotide arrangement and/or electrostatic potential around the cleavage site. It is most likely that a guanine and adenine in close proximity of the cleavage site are involved in general acid-base catalysis. In addition, metal ions appear to play a structural rather than catalytic role although some of their crystal structures have shown a direct metal ion coordination to a non-bridging phosphate oxygen at the cleavage site. Collectively, the structural and biochemical data of the four newest ribozymes could contribute to advance our mechanistic understanding of how self-cleaving ribozymes accomplish their efficient site-specific RNA cleavages.

The Novel Chemical Mechanism of the Twister Ribozyme

We describe the multifactorial origins of catalysis by the twister ribozyme. We provide evidence that the adenine immediately 3' to the scissile phosphate (A1) acts as a general acid. Substitution of ring nitrogen atoms indicates that very unusually the N3 of A1 is the proton donor to the oxyanion leaving group. A1 is accommodated in a specific binding pocket that raises its pKa toward neutrality, juxtaposes its N3 with the O5' to be protonated, and helps create the in-line trajectory required for nucleophilic attack. A1 performs general acid catalysis while G33 acts as a general base. A 100-fold stereospecific phosphorothioate effect at the scissile phosphate is consistent with a significant stabilization of the transition state by the ribozyme, and functional group substitution at G33 indicates that its exocyclic N2 interacts directly with the scissile phosphate. A model of the ribozyme active site is proposed that accommodates these catalytic strategies.

Chemistry and Biology of Self-Cleaving Ribozymes

Self-cleaving ribozymes were discovered 30 years ago, but their biological distribution and catalytic mechanisms are only beginning to be defined. Each ribozyme family is defined by a distinct structure, with unique active sites accelerating the same transesterification reaction across the families. Biochemical studies show that general acid-base catalysis is the most common mechanism of self-cleavage, but metal ions and metabolites can be used as cofactors. Ribozymes have been discovered in highly diverse genomic contexts throughout nature, from viroids to vertebrates. Their biological roles include self-scission during rolling-circle replication of RNA genomes, co-transcriptional processing of retrotransposons, and metabolite-dependent gene expression regulation in bacteria. Other examples, including highly conserved mammalian ribozymes, suggest that many new biological roles are yet to be discovered.

Crystal structure of the Varkud satellite ribozyme

Varkud Satellite (VS) ribozyme mediates rolling circle replication of a plasmid found in the Neurospora mitochondria. We report crystal structures of this ribozyme at 3.1Å resolution, revealing an intertwined dimer formed by an exchange of substrate helices. Within each protomer, an arrangement of three-way helical junctions organizes seven helices into a global fold that creates a docking site for the substrate helix of the other protomer, resulting in the formation of two active sites in trans. This mode of RNA-RNA association resembles the process of domain swapping in proteins and has implications for RNA regulation and evolution. Within each active site, adenine and guanine nucleobases abut the scissile phosphate, poised to serve direct roles in catalysis. Similarities to the active sites of the hairpin and hammerhead ribozymes highlight the functional significance of active site features, underscore the ability of RNA to access functional architectures from distant regions of sequence space, and suggest convergent evolution.

Role of tautomerism in RNA biochemistry

Heterocyclic nucleic acid bases and their analogs can adopt multiple tautomeric forms due to the presence of multiple solvent-exchangeable protons. In DNA, spontaneous formation of minor tautomers has been speculated to contribute to mutagenic mispairings during DNA replication, whereas in RNA, minor tautomeric forms have been proposed to enhance the structural and functional diversity of RNA enzymes and aptamers. This review summarizes the role of tautomerism in RNA biochemistry, specifically focusing on the role of tautomerism in catalysis of small self-cleaving ribozymes and recognition of ligand analogs by riboswitches. Considering that the presence of multiple tautomers of nucleic acid bases is a rare occurrence, and that tautomers typically interconvert on a fast time scale, methods for studying rapid tautomerism in the context of nucleic acids under biologically relevant aqueous conditions are also discussed.

Arginine as a general acid catalyst in serine recombinase-mediated DNA cleavage

Members of the serine family of site-specific DNA recombinases use an unusual constellation of amino acids to catalyze the formation and resolution of a covalent protein-DNA intermediate. A recent high resolution structure of the catalytic domain of Sin, a particularly well characterized family member, provided a detailed view of the catalytic site. To determine how the enzyme might protonate and stabilize the 3'O leaving group in the strand cleavage reaction, we examined how replacing this oxygen with a sulfur affected the cleavage rate by WT and mutant enzymes. To facilitate direct comparison of the cleavage rates, key experiments used suicide substrates that prevented religation after cleavage. The catalytic defect associated with mutation of one of six highly conserved arginine residues, Arg-69 in Sin, was partially rescued by a 3' phosphorothiolate substrate. We conclude that Arg-69 has an important role in stabilizing the 3'O leaving group and is the prime candidate for the general acid that protonates the 3'O, in good agreement with the position it occupies in the high resolution structure of the active site of Sin.

A transition-state interaction shifts nucleobase ionization toward neutrality to facilitate small ribozyme catalysis

One mechanism by which ribozymes can accelerate biological reactions is by adopting folds that favorably perturb nucleobase ionization. Herein we used Raman crystallography to directly measure pK(a) values for the Ade38 N1 imino group of a hairpin ribozyme in distinct conformational states. A transition-state analogue gave a pK(a) value of 6.27 ± 0.05, which agrees strikingly well with values measured by pH-rate analyses. To identify the chemical attributes that contribute to the shifted pK(a), we determined crystal structures of hairpin ribozyme variants containing single-atom substitutions at the active site and measured their respective Ade38 N1 pK(a) values. This approach led to the identification of a single interaction in the transition-state conformation that elevates the base pK(a) > 0.8 log unit relative to the precatalytic state. The agreement of the microscopic and macroscopic pK(a) values and the accompanying structural analysis supports a mechanism in which Ade38 N1(H)+ functions as a general acid in phosphodiester bond cleavage. Overall the results quantify the contribution of a single electrostatic interaction to base ionization, which has broad relevance for understanding how RNA structure can control chemical reactivity.

General acid-base catalysis mediated by nucleobases in the hairpin ribozyme

The catalytic mechanism by which the hairpin ribozyme accelerates cleavage or ligation of the phosphodiester backbone of RNA has been incompletely understood. There is experimental evidence for an important role for an adenine (A38) and a guanine (G8), and it has been proposed that these act in general acid-base catalysis. In this work we show that a large reduction in cleavage rate on substitution of A38 by purine (A38P) can be reversed by replacement of the 5'-oxygen atom at the scissile phosphate by sulfur (5'-PS), which is a much better leaving group. This is consistent with A38 acting as the general acid in the unmodified ribozyme. The rate of cleavage of the 5'-PS substrate by the A38P ribozyme increases with pH log-linearly, indicative of a requirement for a deprotonated base with a relatively high pK(a). On substitution of G8 by diaminopurine, the 5'-PS substrate cleavage rate at first increases with pH and then remains at a plateau, exhibiting an apparent pK(a) consistent with this nucleotide acting in general base catalysis. Alternative explanations for the pH dependence of hairpin ribozyme reactivity are discussed, from which we conclude that general acid-base catalysis by A38 and G8 is the simplest and most probable explanation consistent with all the experimental data.

QM/MM studies of hairpin ribozyme self-cleavage suggest the feasibility of multiple competing reaction mechanisms

The hairpin ribozyme is a prominent member of small ribozymes since it does not require metal ions to achieve catalysis. Guanine 8 (G8) and adenine 38 (A38) have been identified as key participants in self-cleavage and -ligation. We have carried out hybrid quantum-mechanical/molecular mechanical (QM/MM) calculations to evaluate the energy along several putative reaction pathways. The error of our DFT description of the QM region was tested and shown to be ~1 kcal/mol. We find that self-cleavage of the hairpin ribozyme may follow several competing microscopic reaction mechanisms, all with calculated activation barriers in good agreement with those from experiment (20-21 kcal/mol). The initial nucleophilic attack of the A-1(2'-OH) group on the scissile phosphate is predicted to be rate-limiting in all these mechanisms. An unprotonated G8(-) (together with A38H(+)) yields a feasible activation barrier (20.4 kcal/mol). Proton transfer to a nonbridging phosphate oxygen also leads to feasible reaction pathways. Finally, our calculations consider thio-substitutions of one or both nonbridging oxygens of the scissile phosphate and predict that they have only a negligible effect on the reaction barrier, as observed experimentally.

Modeling the RNA 2'OH activation: possible roles of metal ion and nucleobase as catalysts in self-cleaving ribozymes

The RNA 2'OH activation as taking place in the first chemical step of self-cleaving ribozymes is studied theoretically by DFT and MP2 methods using a continuum solvation model (CPCM). The reaction of proton transfer is studied in the presence of two kinds of catalysts: a fully hydrated metal ion (Mg(2+)) or partially hydrated nucleobase (guanine), taken separately or together leading to three different modes of activation. The metal ion is either directly bound (inner-sphere) or indirectly bound (outer-sphere) to the 2'OH group and a hydroxide ion acts as a general or specific base; the nucleobase is taken in anionic or in neutral enol-tautomeric forms playing itself the role of general base. The presence of a close metal ion (outer-sphere) lowers the pK(a) value of the 2'OH group by several log units in both metal-ion and nuleobase catalysis. The direct metal coordination to the 2'OH group (inner-sphere) further stabilizes the developing negative charge on the nucleophile. The switching from the inner-sphere to the outer-sphere coordination appears to be driven by the energy cost for reorganizing the first coordination shell rather than by the electrostatic repulsion between the ligands. The metal-ion catalysis is more effective with a specific base in the dianionic mechanism. On the other hand, the nucleobase catalysis is more effective in the monoanionic mechanism and in the presence of a metal ion acting as a cofactor through nonspecific electrostatic interactions. The results establish a baseline to study the possible roles of metal and nucleobase catalysts and their environment in more realistic models for self-cleaving ribozymes.

Catalysis by the nucleolytic ribozymes

The nucleolytic ribozymes use general acid-base catalysis to contribute significantly to their rate enhancement. The VS (Varkud satellite) ribozyme uses a guanine and an adenine nucleobase as general base and acid respectively in the cleavage reaction. The hairpin ribozyme is probably closely similar, while the remaining nucleolytic ribozymes provide some interesting contrasts.

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.

Nucleobase-mediated general acid-base catalysis in the Varkud satellite ribozyme

Existing evidence suggests that the Varkud satellite (VS) ribozyme accelerates the cleavage of a specific phosphodiester bond using general acid-base catalysis. The key functionalities are the nucleobases of adenine 756 in helix VI of the ribozyme, and guanine 638 in the substrate stem loop. This results in a bell-shaped dependence of reaction rate on pH, corresponding to groups with pK(a) = 5.2 and 8.4. However, it is not possible from those data to determine which nucleobase is the acid, and which the base. We have therefore made substrates in which the 5' oxygen of the scissile phosphate is replaced by sulfur. This labilizes the leaving group, removing the requirement for general acid catalysis. This substitution restores full activity to the highly impaired A756G ribozyme, consistent with general acid catalysis by A756 in the unmodified ribozyme. The pH dependence of the cleavage of the phosphorothiolate-modified substrates is consistent with general base catalysis by nucleobase at position 638. We conclude that cleavage of the substrate by the VS ribozyme is catalyzed by deprotonation of the 2'-O nucleophile by G638 and protonation of the 5'-O leaving group by A756.

Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme

The hairpin ribozyme is a prominent member of the group of small catalytic RNAs (RNA enzymes or ribozymes) because it does not require metal ions to achieve catalysis. Biochemical and structural data have implicated guanine 8 (G8) and adenine 38 (A38) as catalytic participants in cleavage and ligation catalyzed by the hairpin ribozyme, yet their exact role in catalysis remains disputed. To gain insight into dynamics in the active site of a minimal self-cleaving hairpin ribozyme, we have performed extensive classical, explicit-solvent molecular dynamics (MD) simulations on time scales of 50-150 ns. Starting from the available X-ray crystal structures, we investigated the structural impact of the protonation states of G8 and A38, and the inactivating A-1(2'-methoxy) substitution employed in crystallography. Our simulations reveal that a canonical G8 agrees well with the crystal structures while a deprotonated G8 profoundly distorts the active site. Thus MD simulations do not support a straightforward participation of the deprotonated G8 in catalysis. By comparison, the G8 enol tautomer is structurally well tolerated, causing only local rearrangements in the active site. Furthermore, a protonated A38H(+) is more consistent with the crystallography data than a canonical A38. The simulations thus support the notion that A38H(+) is the dominant form in the crystals, grown at pH 6. In most simulations, the canonical A38 departs from the scissile phosphate and substantially perturbs the structures of the active site and S-turn. Yet, we occasionally also observe formation of a stable A-1(2'-OH)...A38(N1) hydrogen bond, which documents the ability of the ribozyme to form this hydrogen bond, consistent with a potential role of A38 as general base catalyst. The presence of this hydrogen bond is, however, incompatible with the expected in-line attack angle necessary for self-cleavage, requiring a rapid transition of the deprotonated 2'-oxyanion to a position more favorable for in-line attack after proton transfer from A-1(2'-OH) to A38(N1). The simulations revealed a potential force field artifact, occasional but irreversible formation of "ladder-like", underwound A-RNA structure in one of the external helices. Although it does not affect the catalytic center of the hairpin ribozyme, further studies are under way to better assess possible influence of such force field behavior on long RNA simulations.

Analysis of RNA folding by native polyacrylamide gel electrophoresis

Polyacrylamide gel electrophoresis under native conditions (native PAGE) is a well-established and versatile method for probing nucleic acid conformation and nucleic acid-protein interactions. Native PAGE has been used to measure RNA folding equilibria and kinetics under a wide variety of conditions. Advantages of this method are its adaptability, absolute determination of reaction endpoints, and direct analysis of conformational hetereogeneity within a sample. Native PAGE is also useful for resolving ligand-induced structural changes.

Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM

Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.

Exploring ribozyme conformational changes with X-ray crystallography

Relating three-dimensional fold to function is a central challenge in RNA structural biology. Toward this goal, X-ray crystallography has long been considered the "gold standard" for structure determinations at atomic resolution, although NMR spectroscopy has become a powerhouse in this arena as well. In the area of dynamics, NMR remains the dominant technique to probe the magnitude and timescales of molecular motion. Although the latter area remains largely unassailable by conventional crystallographic methods, inroads have been made on proteins using Laue radiation on timescales of ms to ns. Proposed 'fourth generation' radiation sources, such as free-electron X-ray lasers, promise ps- to fs-timescale resolution, and credible evidence is emerging that supports the feasibility of single molecule imaging. At present however, the preponderance of RNA structural information has been derived from timescale and motion insensitive crystallographic techniques. Importantly, developments in computing, automation and high-flux synchrotron sources have propelled the rapidity of 'conventional' RNA crystal structure determinations to timeframes of hours once a suitable set of phases is obtained. With a sufficient number of crystal structures, it is possible to create a structural ensemble that can provide insight into global and local molecular motion characteristics that are relevant to biological function. Here we describe techniques to explore conformational changes in the hairpin ribozyme, a representative non-protein-coding RNA catalyst. The approaches discussed include: (i) construct choice and design using prior knowledge to improve X-ray diffraction; (ii) recognition of long-range conformational changes and (iii) use of single-base or single-atom changes to create ensembles. The methods are broadly applicable to other RNA systems.

Comparative enzymology and structural biology of RNA self-cleavage

Self-cleaving hammerhead, hairpin, hepatitis delta virus, and glmS ribozymes comprise a family of small catalytic RNA motifs that catalyze the same reversible phosphodiester cleavage reaction, but each motif adopts a unique structure and displays a unique array of biochemical properties. Recent structural, biochemical, and biophysical studies of these self-cleaving RNAs have begun to reveal how active site nucleotides exploit general acid-base catalysis, electrostatic stabilization, substrate destabilization, and positioning and orientation to reduce the free energy barrier to catalysis. Insights into the variety of catalytic strategies available to these model RNA enzymes are likely to have important implications for understanding more complex RNA-catalyzed reactions fundamental to RNA processing and protein synthesis.

Identification of an imino group indispensable for cleavage by a small ribozyme

The hairpin ribozyme is a small, noncoding RNA (ncRNA) that catalyzes a site-specific phosphodiester bond cleavage reaction. Prior biochemical and structural analyses pinpointed the amidine moiety of base Ade38 as a key functional group in catalysis, but base changes designed to probe function resulted in localized misfolding of the active site. To define the requirements for chemical activity using a conservative modification, we synthesized and incorporated N1-deazaadenosine into the full-length ribozyme construct. This single-atom variant severely impairs activity, although the active-site fold remains intact in the accompanying crystal structures. The results demonstrate the essentiality of the imino moiety as well as the importance of its interaction with the substrate in the precatalytic and transition-state conformations. This work demonstrates the efficacy of single-atom approaches in the analysis of ncRNA structure-function relationships.

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.

Direct measurement of the ionization state of an essential guanine in the hairpin ribozyme

Active site guanines are critical for self-cleavage reactions of several ribozymes, but their precise functions in catalysis are unclear. To learn whether protonated or deprotonated forms of guanine predominate in the active site, microscopic pKa values were determined for ionization of 8-azaguanosine substituted for G8 in the active site of a fully functional hairpin ribozyme in order to determine microscopic pKa values for 8-azaguanine deprotonation from the pH dependence of fluorescence. Microscopic pKa values above 9 for deprotonation of 8-azaguanine in the active site were about 3 units higher than apparent pKa values determined from the pH dependence of self-cleavage kinetics. Thus, the increase in activity with increasing pH does not correlate with deprotonation of G8, and most of G8 is protonated at neutral pH. These results do not exclude a role in proton transfer, but a simple interpretation is that G8 functions in the protonated form, perhaps by donating hydrogen bonds.

The Small Ribozymes: Common and Diverse Features Observed through the FRET Lens

The hammerhead, hairpin, HDV, VS and glmS ribozymes are the five known, naturally occurring catalytic RNAs classified as the "small ribozymes". They share common reaction chemistry in cleaving their own backbone by phosphodiester transfer, but are diverse in their secondary and tertiary structures, indicating that Nature has found at least five independent solutions to a common chemical task. Fluorescence resonance energy transfer (FRET) has been extensively used to detect conformational changes in these ribozymes and dissect their reaction pathways. Common and diverse features are beginning to emerge that, by extension, highlight general biophysical properties of non-protein coding RNAs.

RNA intramolecular dynamics by single-molecule FRET

Investigation of single RNA molecules using fluorescence resonance energy transfer (FRET) is a powerful approach to investigate dynamic and thermodynamic aspects of the folding process of a given RNA. Its application requires interdisciplinary work from the fields of chemistry, biochemistry, and physics. The present work gives detailed instructions on the synthesis of RNA molecules labeled with two fluorescent dyes interacting by FRET, as well as on their investigation by single-molecule fluorescence spectroscopy.

Mechanistic characterization of the HDV genomic ribozyme: solvent isotope effects and proton inventories in the absence of divalent metal ions support C75 as the general acid

The hepatitis delta virus (HDV) ribozyme uses the nucleobase C75 and a hydrated Mg(2+) ion as the general acid-base catalysts in phosphodiester bond cleavage at physiological salt. A mechanistic framework has been advanced that involves one Mg(2+)-independent and two Mg(2+)-dependent channels. The rate-pH profile for wild-type (WT) ribozyme in the Mg(2+)-free channel is inverted relative to the fully Mg(2+)-dependent channel, with each having a near-neutral pKa. Inversion of the rate-pH profile was used as the crux of a mechanistic argument that C75 serves as general acid both in the presence and absence of Mg(2+). However, subsequent studies on a double mutant (DM) ribozyme suggested that the pKa observed for WT in the absence of Mg(2+) arises from ionization of C41, a structural nucleobase. To investigate this further, we acquired rate-pH/pD profiles and proton inventories for WT and DM in the absence of Mg(2+). Corrections were made for effects of ionic strength on hydrogen ion activity and pH meter readings. Results are accommodated by a model wherein the Mg(2+)-free pKa observed for WT arises from ionization of C75, and DM reactivity is compromised by protonation of C41. The Brønsted base appears to be water or hydroxide ion depending on pH. The observed pKa's are related to salt-dependent pH titrations of a model oligonucleotide, as well as electrostatic calculations, which support the local environment for C75 in the absence of Mg(2+) being similar to that in the presence of Mg(2+) and impervious to bulk ions. Accordingly, the catalytic role of C75 as the general acid does not appear to depend on divalent ions or the identity of the Brønsted base.

The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid-base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK (a) values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK (a) values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid-base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid-base catalysis indicate that k (1), the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK (a). The inference that k (1) is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK (a) is that of the general acid and the higher pK (a) is that of the general base.

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.

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.

Electrostatic interactions in the hairpin ribozyme account for the majority of the rate acceleration without chemical participation by nucleobases

Molecular dynamics simulations using a combined quantum mechanical/molecular mechanical potential are used to determine the two-dimensional free energy profiles for the mechanism of RNA transphosphorylation in solution and catalyzed by the hairpin ribozyme. A mechanism is explored whereby the reaction proceeds without explicit chemical participation by conserved nucleobases in the active site. The ribozyme lowers the overall free energy barrier by up to 16 kcal/mol, accounting for the majority of the observed rate enhancement. The barrier reduction in this mechanism is achieved mainly by the electrostatic environment provided by the ribozyme without recruitment of active site nucleobases as acid or base catalysts. The results establish a baseline mechanism that invokes only the solvation and specific hydrogen-bonding interactions present in the ribozyme active site and provide a departure point for the exploration of alternate mechanisms where nucleobases play an active chemical role.

Catalytic strategies of self-cleaving ribozymes

[Structure: see text]. Five naturally occurring nucleolytic ribozymes have been identified: the hammerhead, hairpin, glmS, hepatitis delta virus (HDV), and Varkud satellite (VS) ribozymes. All of these RNA enzymes catalyze self-scission of the RNA backbone using a chemical mechanism equivalent to that of RNase A. RNase A uses four basic strategies to promote this reaction: geometric constraints, activation of the nucleophile, transition-state stabilization, and leaving group protonation. In this Account, we discuss the current thinking on how nucleolytic ribozymes harness RNase A's four sources of catalytic power. The geometry of the phosphodiester cleavage reaction constrains the nucleotides flanking the scissile phosphate so that they are unstacked from a canonical A-form helix and thus require alternative stabilization. Crystal structures and mutational analysis reveal that cross-strand base pairing, along with unconventional stacking and tertiary hydrogen-bonding interactions, work to stabilize the splayed conformation in nucleolytic ribozymes. Deprotonation of the 2'-OH nucleophile greatly increases its nucleophilicity in the strand scission reaction. Crystal structures of the hammerhead, hairpin, and glmS ribozymes reveal the N1 of a G residue within hydrogen-bonding distance of the 2'-OH. In each case, this residue has also been shown to be important for catalysis. In the HDV ribozyme, a hydrated magnesium has been implicated as the general base. Catalysis by the VS ribozyme requires both an A and a G, but the precise role of either has not been elucidated. Enzymes can lower the energy of a chemical reaction by binding more tightly to the transition state than to the ground states. Comparison of the hairpin ground- and transition-state mimic structures reveal greater hydrogen bonding to the transition-state mimic structure, suggesting transition-state stabilization as a possible catalytic strategy. However, the hydrogen-bonding pattern in the glmS ribozyme transition-state mimic structure and the ground-state structures are equivalent. Protonation of the 5'-O leaving group by a variety of functional groups can promote the cleavage reaction. In the HDV ribozyme, the general acid is a conserved C residue. In the hairpin ribozyme, a G residue has been implicated in protonation of the leaving group. An A in the hammerhead ribozyme probably plays a similar role. In the glmS ribozyme, an exogenous cofactor may provide the general acid. This diversity is in contrast to the relatively small number of functional groups that serve as a general base, where at least three of the nucleolytic ribozymes may use the N1 of a G.

A contact photo-cross-linking investigation of the active site of the 8-17 deoxyribozyme

The small RNA-cleaving 8-17 deoxyribozyme (DNAzyme) has been the subject of extensive mechanistic and structural investigation, including a number of recent single-molecule studies of its global folding. Little detailed insight exists, however, into this DNAzyme's active site; for instance, the identity of specific nucleotides that are proximal to or in contact with the scissile site in the substrate. Here, we report a systematic replacement of a number of bases within the magnesium-folded DNAzyme-substrate complex with thio- and halogen-substituted base analogues, which were then photochemically activated to generate contact cross-links within the complex. Mapping of the cross-links revealed a striking pattern of DNAzyme-substrate cross-links but an absence of significant intra-DNAzyme cross-links. Notably, the two nucleotides directly flanking the scissile phosphodiester cross-linked strongly with functionally important elements within the DNAzyme, the thymine of a G.T wobble base pair, a WCGR bulge loop, and a terminal AGC loop. Mutation of the wobble base pair to a G-C pair led to a significant folding instability of the DNAzyme-substrate complex. The cross-linking patterns obtained were used to generate a model for the DNAzyme's active site that had the substrate's scissile phosphodiester sandwiched between the DNAzyme's wobble thymine and its AGC and WCGR loops.

Probing RNA structural dynamics and function by fluorescence resonance energy transfer (FRET)

Biological function of RNA is often mediated by cyclic switching between several (meta-)stable arrangements of tertiary structure. Fluorophore labeling of RNA offers a unique view into these folding and conformational switching events, since a fluorescence signal is sensitive to its molecular environment and can be continuously monitored in real time to produce kinetic rate information. This unit focuses on the practical implications of using fluorescence resonance energy transfer (FRET) to probe RNA structural dynamics and function. FRET is a particularly powerful fluorescence technique since, in addition to kinetic data, it provides insights into the structural basis of a conformational rearrangement. Protocols describe how to postsynthetically label RNA for FRET and how to acquire and analyze FRET data. Support protocols describe methods for deprotecting synthetic RNA and for purifying RNA by gel electrophoresis and HPLC. Considerations for selecting appropriate RNA, fluorophores, and labeling strategies are discussed in detail in the commentary.

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.

Shared traits on the reaction coordinates of ribonuclease and an RNA enzyme

Reaction-intermediate analogs have been used to understand how phosphoryl transfer enzymes promote catalysis. Herein we report the first structure of a small ribozyme crystallized with a 3'-OH, 2',5'-linkage in lieu of the normal phosphodiester substrate. The new structure shares features of the reaction coordinate exhibited in prior ribozyme structures including a vanadate complex that mimicked the oxyphosphorane transition state. As such, the structure exhibits reaction-intermediate traits that allow direct comparison of stabilizing interactions to the 3'-OH, 2',5'-linkage contributed by the RNA enzyme and its protein counterpart, ribonuclease. Clear similarities are observed between the respective structures including hydrogen bonds to the non-bridging oxygens of the scissile phosphate. Other commonalities include carefully poised water molecules that may alleviate charge build-up in the transition state and placement of a positive charge near the leaving group. The advantages of 2',5'-linkages to investigate phosphoryl-transfer reactions are discussed, and argue for their expanded use in structural studies.