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

Simulation-Guided Conformational Space Exploration to Assess Reactive Conformations of a Ribozyme

Self-splicing ribozymes are small ribonucleic acid (RNA) enzymes that catalyze their own cleavage through a transphosphoesterification reaction. While this process is involved in some specific steps of viral RNA replication and splicing, it is also of importance in the context of the (putative) first autocatalytic RNA-based systems that could have preceded the emergence of modern life. The uncatalyzed phosphoester bond formation is thermodynamically very unfavorable, and many experimental studies have focused on understanding the molecular features of catalysis in these ribozymes. However, chemical reaction paths are short-lived and not easily characterized by experimental approaches, so molecular simulation approaches appear as an ideal tool to unveil the molecular details of the reaction. Here, we focus on the model hairpin ribozyme. We show that identifying a relevant initial conformation for reactivity studies, which is frequently overlooked in mixed quantum-classical studies that predominantly concentrate on the chemical reaction itself, can be highly challenging. These challenges stem from limitations in both available experimental structures (which are chemically altered to prevent self-cleavage) and the accuracy of force fields, together with the necessity for comprehensive sampling. We show that molecular dynamics simulations, combined with extensive conformational phase space exploration with Hamiltonian replica-exchange simulations, enable us to characterize the relevant conformational basins of the minimal hairpin ribozyme in the ligated state prior to self-cleavage. We find that what is usually considered a canonical reactive conformation with active site geometries and hydrogen-bond patterns that are optimal for the addition-elimination reaction with general acid/general base catalysis is metastable and only marginally populated. The thermodynamically stable conformation appears to be consistent with the expectations of a mechanism that does not require the direct participation of ribozyme residues in the reaction. While these observations may suffer from forcefield inaccuracies, all investigated forcefields lead to the same conclusions upon proper sampling, contrasting with previous investigations on shorter timescales suggesting that at least one reparametrization of the Amber99 forcefield allowed to stabilize aligned active site conformations. Our study demonstrates that identifying the most pertinent reactant state conformation holds equal importance alongside the accurate determination of the thermodynamics and kinetics of the chemical steps of the reaction.

Scaling Catalytic Contributions of Small Self-Cleaving Ribozymes

Nucleolytic ribozymes utilize general acid-base catalysis to perform phosphodiester cleavage. In most ribozyme classes, a conserved active site guanosine is positioned to act as general base, thereby activating the 2'-OH group to attack the scissile phosphate (γ-catalysis). Here, we present an atomic mutagenesis study for the pistol ribozyme class. Strikingly, "general base knockout" by replacement of the guanine N1 atom by carbon results in only 2.7-fold decreased rate. Therefore, the common view that γ-catalysis critically depends on the N1 moiety becomes challenged. For pistol ribozymes we found that γ-catalysis is subordinate in overall catalysis, made up by two other catalytic factors (α and δ). Our approach allows scaling of the different catalytic contributions (α, β, γ, δ) with unprecedented precision and paves the way for a thorough mechanistic understanding of nucleolytic ribozymes with active site guanines.

Scaling Catalytic Contributions of Small Self-Cleaving Ribozymes

Nucleolytic ribozymes utilize general acid-base catalysis to perform phosphodiester cleavage. In most ribozyme classes, a conserved active site guanosine is positioned to act as general base, thereby activating the 2'-OH group to attack the scissile phosphate (γ-catalysis). Here, we present an atomic mutagenesis study for the pistol ribozyme class. Strikingly, "general base knockout" by replacement of the guanine N1 atom by carbon results in only 2.7-fold decreased rate. Therefore, the common view that γ-catalysis critically depends on the N1 moiety becomes challenged. For pistol ribozymes we found that γ-catalysis is subordinate in overall catalysis, made up by two other catalytic factors (α and δ). Our approach allows scaling of the different catalytic contributions (α, β, γ, δ) with unprecedented precision and paves the way for a thorough mechanistic understanding of nucleolytic ribozymes with active site guanines.

Cellular Small Molecules Contribute to Twister Ribozyme Catalysis

The number of self-cleaving small ribozymes has increased sharply in recent years. Advances have been made in describing these ribozymes in terms of four catalytic strategies: α describes in-line attack, β describes neutralization of the nonbridging oxygens, γ describes activation of the nucleophile, and δ describes stabilization of the leaving group. Current literature presents the rapid self-cleavage of the twister ribozyme in terms of all four of these classic catalytic strategies. Herein, we describe the nonspecific contribution of small molecules to ribozyme catalysis. At biological pH, the rate of the wild-type twister ribozyme is enhanced up to 5-fold in the presence of moderate buffer concentrations, similar to the 3-5-fold effects reported previously for buffer catalysis for protein enzymes. We observe this catalytic enhancement not only with standard laboratory buffers, but also with diverse biological small molecules, including imidazole, amino acids, and amino sugars. Brønsted plots suggest that small molecules assist in proton transfer, most likely with δ catalysis. Cellular small molecules provide a simple way to overcome the limited functional diversity of RNA and have the potential to participate in the catalytic mechanisms of many ribozymes in vivo.

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.

Structural basis for the fast self-cleavage reaction catalyzed by the twister ribozyme

Twister is a recently discovered RNA motif that is estimated to have one of the fastest known catalytic rates of any naturally occurring small self-cleaving ribozyme. We determined the 4.1-Å resolution crystal structure of a twister sequence from an organism that has not been cultured in isolation, and it shows an ordered scissile phosphate and nucleotide 5' to the cleavage site. A second crystal structure of twister from Orzyza sativa determined at 3.1-Å resolution exhibits a disordered scissile phosphate and nucleotide 5' to the cleavage site. The core of twister is stabilized by base pairing, a large network of stacking interactions, and two pseudoknots. We observe three nucleotides that appear to mediate catalysis: a guanosine that we propose deprotonates the 2'-hydroxyl of the nucleotide 5' to the cleavage site and a conserved adenosine. We suggest the adenosine neutralizes the negative charge on a nonbridging phosphate oxygen atom at the cleavage site. The active site also positions the labile linkage for in-line nucleophilic attack, and thus twister appears to simultaneously use three strategies proposed for small self-cleaving ribozymes. The twister crystal structures (i) show its global structure, (ii) demonstrate the significance of the double pseudoknot fold, (iii) provide a possible hypothesis for enhanced catalysis, and (iv) illuminate the roles of all 10 highly conserved nucleotides of twister that participate in the formation of its small and stable catalytic pocket.

Comparison of ab Initio, DFT, and Semiempirical QM/MM Approaches for Description of Catalytic Mechanism of Hairpin Ribozyme

We have analyzed the capability of state-of-the-art multiscale computational approaches to provide atomic-resolution electronic structure insights into possible catalytic scenarios of the hairpin ribozyme by evaluating potential and free energy surfaces of the reactions by various hybrid QM/MM methods. The hairpin ribozyme is a unique catalytic RNA that achieves rate acceleration similar to other small self-cleaving ribozymes but without direct metal ion participation. Guanine 8 (G8) and adenine 38 (A38) have been identified as the catalytically essential nucleobases. However, their exact catalytic roles are still being investigated. In line with the available experimental data, we considered two reaction scenarios involving protonated A38H(+) as a general acid which is further assisted by either canonical G8 or deprotonated G8(-) forms. We used the spin-component scaled Møller-Plesset (SCS-MP2) method at the complete basis set limit as the reference method. The semiempirical AM1/d-PhoT and SCC-DFTBPR methods provided acceptable activation barriers with respect to the SCS-MP2 data but predicted significantly different reaction pathways. DFT functionals (BLYP and MPW1K) yielded the same reaction pathway as the SCS-MP2 method. The activation barriers were slightly underestimated by the GGA BLYP functional, although with accuracy comparable to the semiempirical methods. The SCS-MP2 method and hybrid MPW1K functional gave activation barriers that were closest to those derived from experimentally measured rate constants.

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.

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.

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.

The pH dependence of hairpin ribozyme catalysis reflects ionization of an active site adenine

Understanding how self-cleaving ribozymes mediate catalysis is crucial in light of compelling evidence that human and bacterial gene expression can be regulated through RNA self-cleavage. The hairpin ribozyme catalyzes reversible phosphodiester bond cleavage through a mechanism that does not require divalent metal cations. Previous structural and biochemical evidence implicated the amidine group of an active site adenosine, A38, in a pH-dependent step in catalysis. We developed a way to determine microscopic pK(a) values in active ribozymes based on the pH-dependent fluorescence of 8-azaadenosine (8azaA). We compared the microscopic pK(a) for ionization of 8azaA at position 38 with the apparent pK(a) for the self-cleavage reaction in a fully functional hairpin ribozyme with a unique 8azaA at position 38. Microscopic and apparent pK(a) values were virtually the same, evidence that A38 protonation accounts for the decrease in catalytic activity with decreasing pH. These results implicate the neutral unprotonated form of A38 in a transition state that involves formation of the 5'-oxygen-phosphorus bond.

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.

The glmS ribozyme: use of a small molecule coenzyme by a gene-regulatory RNA

The glmS ribozyme is the first known example of a natural ribozyme that has evolved to require binding of an exogenous small molecule for activity. In Gram-positive bacteria, this RNA domain is part of the messenger RNA (mRNA) encoding the essential enzyme that synthesizes glucosamine-6-phosphate (GlcN6P). When present at physiologic concentration, this small molecule binds to the glmS ribozyme and uncovers a latent self-cleavage activity that ultimately leads to degradation of the mRNA. Biochemical and structural studies reveal that the RNA adopts a rigid fold stabilized by three pseudoknots and the packing of a peripheral domain against the ribozyme core. GlcN6P binding to this pre-organized RNA does not induce conformational changes; rather, the small molecule functions as a coenzyme, providing a catalytically essential amine group to the active site. The ribozyme is not a passive player, however. Active site functional groups are essential for catalysis, even in the presence of GlcN6P. In addition to being a superb experimental system with which to analyze how RNA catalysts can exploit small molecule coenzymes to broaden their chemical versatility, the presence of the glmS ribozyme in numerous pathogenic bacteria make this RNA an attractive target for the development of new antibiotics and antibacterial strategies.

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.

Small self-cleaving ribozymes

The hammerhead, hairpin, hepatitis delta virus (HDV), Varkud Satellite (VS), and glmS ribozymes catalyze sequence-specific intramolecular cleavage of RNA. They range between 50 and 150 nucleotides in length, and are known as the "small self-cleaving ribozymes." Except for the glmS ribozyme that functions as a riboswitch in Gram-positive bacteria, they were originally discovered as domains of satellite RNAs. However, recent studies show that several of them are broadly distributed in genomes of organisms from many phyla. Each of these ribozymes has a unique overall architecture and active site organization. Crystal structures have revealed how RNA active sites can bind preferentially to the transition state of a reaction, whereas mechanistic studies have shown that nucleobases can efficiently perform general acid-base and electrostatic catalysis. This versatility explains the abundance of ribozymes in contemporary organisms and also supports a role for catalytic RNAs early in evolution.

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.

Direct Raman measurement of an elevated base pKa in the active site of a small ribozyme in a precatalytic conformation

Catalytic RNA molecules can achieve rate acceleration by shifting base pK(a) values toward neutrality. Prior evidence has suggested that base A38 of the hairpin ribozyme plays an important role in phosphoryl transfer, possibly functioning as a general acid, or by orienting a specific water molecule for proton transfer. To address the role of A38, we used Raman spectroscopy to measure directly the pK(a) of the N1-imino moiety in the context of hairpin ribozyme crystals representative of a "precatalytic" conformation. The results revealed that the pK(a) of A38 is shifted to 5.46 +/- 0.05 relative to 3.68 +/- 0.06 derived from a reference solution of the nucleotide AMP. The elevated pK(a) correlates well with the first titration point of the macroscopic pH-rate profile of the hairpin ribozyme in solution and strongly supports A38 as a general acid catalyst in bond scission. The results confirm that A38 is protonated before the transition state, which would promote phosphorane development. Overall, the results establish a cogent structure-function paradigm that expands our understanding of how RNA structure can enhance nucleobase reactivity to catalyze biological reactions.

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.

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.

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.

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.

Focus on function: single molecule RNA enzymology

The ability of RNA to catalyze chemical reactions was first demonstrated 25 years ago with the discovery that group I introns and RNase P function as RNA enzymes (ribozymes). Several additional ribozymes were subsequently identified, most notably the ribosome, followed by intense mechanistic studies. More recently, the introduction of single molecule tools has dissected the kinetic steps of several ribozymes in unprecedented detail and has revealed surprising heterogeneity not evident from ensemble approaches. Still, many fundamental questions of how RNA enzymes work at the molecular level remain unanswered. This review surveys the current status of our understanding of RNA catalysis at the single molecule level and discusses the existing challenges and opportunities in developing suitable assays.

Self-association of adenine-dependent hairpin ribozymes

Hairpin ribozymes are flexible molecules that catalyse reversible self-cleavage after the docking of two independently folded internal loops, A and B. The activities, self-association and structures in solution of two 85 base adenine-dependent hairpin ribozymes (ADHR1 and ADHR2) were studied by native gel electrophoresis, analytical centrifugation, and small angle neutron scattering. Bi-molecular RNA interactions such as linear-linear, loop-loop, loop-linear or kissing interactions have been found to be important in the control of various biological functions, and hairpin loops present rich potential for establishing both intra- and intermolecular interactions through standard Watson-Crick base pairing or non-canonical interactions. Similar results were obtained for ADHR1 and ADHR2. At room temperature, they indicated end-to-end self-association of the ribozymes in rod-like structures with a cross-section corresponding to two double strands side-by-side. Dimers, which predominate at low concentration ( approximately 0.1 mg/ml), associate into longer rods, with increasing concentration ( approximately 1 mg/ml). Above 65 degrees C, the dimers and rods dissociated into compact monomers, with a radius of gyration similar to that of tRNA (about 70 bases). The dimers were non-active for catalysis, which suggests that dimer formation, probably by preventing the correct docking of loops A and B, could act as an inhibition mechanism for the regulation of hairpin ribozyme catalysis.

A chemo-genetic approach for the study of nucleobase participation in nucleolytic ribozymes

A novel chemo-genetic approach for the analysis of general acid-base catalysis by nucleobases in ribozymes is reviewed. This involves substitution of a C-nucleoside with imidazole in place of a natural nucleobase. The Varkud satellite ribozyme in which the nucleobase at the critical 756 position has been replaced by imidazole is active in both cleavage and ligation reactions. Similarly, a modified hairpin ribozyme with the nucleobase at position 8 substituted by imidazole is active in cleavage and ligation reactions. Although the rates are lower than those of the natural ribozymes, they are significantly greater than other variants at these positions. The dependence of the hairpin ribozyme reaction rates on pH has been studied. Both cleavage and ligation reactions display a bell-shaped pH dependence, consistent with general acid-base catalysis involving the nucleotide at position 8.

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.

Functional Analysis of Hairpin Ribozyme Active Site Architecture

The hairpin ribozyme is a small catalytic motif found in plant satellite RNAs where it catalyzes a reversible self-cleavage reaction during processing of replication intermediates. Crystallographic studies of hairpin ribozymes have provided high resolution views of the RNA functional groups that comprise the active site and stimulated biochemical studies that probed the contributions of nucleobase functional groups to catalytic chemistry. The dramatic loss of activity that results from perturbation of active site architecture points to the importance of positioning and orientation in catalytic rate acceleration. The current study focuses on the network of noncovalent interactions that align nucleophilic and leaving group oxygens in the orientation required for the S(N)2-type reaction mechanism and orient the active site nucleobases near the reactive phosphate to facilitate catalytic chemistry. Nucleotide modifications that alter or eliminate individual hydrogen bonding partners had different effects on the activation barrier to catalysis, the stability of ribozyme complexes in the ground state, and the internal equilibrium between cleavage and ligation of bound products. Furthermore, substitution of hydrogen bond donors and acceptors with seemingly equivalent pairs sometimes had very different functional consequences. These biochemical analyses augment high resolution structural information to provide insights into the functional significance of active site architecture.

Calculation of pKas in RNA: on the structural origins and functional roles of protonated nucleotides

pK(a) calculations based on the Poisson-Boltzmann equation have been widely used to study proteins and, more recently, DNA. However, much less attention has been paid to the calculation of pK(a) shifts in RNA. There is accumulating evidence that protonated nucleotides can stabilize RNA structure and participate in enzyme catalysis within ribozymes. Here, we calculate the pK(a) shifts of nucleotides in RNA structures using numerical solutions to the Poisson-Boltzmann equation. We find that significant shifts are predicted for several nucleotides in two catalytic RNAs, the hairpin ribozyme and the hepatitis delta virus ribozyme, and that the shifts are likely to be related to their functions. We explore how different structural environments shift the pK(a)s of nucleotides from their solution values. RNA structures appear to use two basic strategies to shift pK(a)s: (a) the formation of compact structural motifs with structurally-conserved, electrostatic interactions; and (b) the arrangement of the phosphodiester backbone to focus negative electrostatic potential in specific regions.

Structural basis of glmS ribozyme activation by glucosamine-6-phosphate

The glmS ribozyme is the only natural catalytic RNA known to require a small-molecule activator for catalysis. This catalytic RNA functions as a riboswitch, with activator-dependent RNA cleavage regulating glmS messenger RNA expression. We report crystal structures of the glmS ribozyme in precleavage states that are unliganded or bound to the competitive inhibitor glucose-6-phosphate and in the postcleavage state. All structures superimpose closely, revealing a remarkably rigid RNA that contains a preformed active and coenzyme-binding site. Unlike other riboswitches, the glmS ribozyme binds its activator in an open, solvent-accessible pocket. Our structures suggest that the amine group of the glmS ribozyme-bound coenzyme performs general acid-base and electrostatic catalysis.

Fluorescent N2,N3-epsilon-adenine nucleoside and nucleotide probes: synthesis, spectroscopic properties, and biochemical evaluation

N1,N(6)-ethenoadenine, epsilon-A, nucleos(t)ides have been previously applied as fluorescent probes in numerous biochemical systems. However, these epsilon-A analogues lack the H-bonding capability of adenine. To improve the fluorescence characteristics while preserving the H-bonding pattern required for molecular recognition, we designed a novel probe: N(2),N3-etheno-adenosine, (N(2),N3-epsilon-A). Here, we describe four novel syntheses of the target epsilon-nucleoside and related analogues. These methods are short, facile, and provide the product regiospecifically. In addition, we report the absorption and emission spectra of N(2),N3-epsilon-A and the dependence of the spectral features on the pH and polarity of the medium. Specifically, maximum emission of N(2),N3-epsilon-A in water is observed at 420 nm (phi=0.03, excitation at 290 nm). The biochemical relevance of the new probe was evaluated with respect to the P2Y(1) receptor and NTPDases 1 and 2. N(2),N3-epsilon-ATP was found to be almost equipotent with ATP at the P2Y(1) receptor and was hydrolyzed by NTPDases 1 and 2 at about 80 % of the rate of ATP. Furthermore, protein binding does not seem to shift the fluorescence of N(2),N3-epsilon-ATP. Based on the fluorescence and full recognition by ATP-binding proteins, we propose N(2),N3-epsilon-ATP and related nucleo(s)tides as unique probes for the investigation of adenine nucleo(s)tide-binding proteins as well as for other biochemical applications.

Chemical rescue, multiple ionizable groups, and general acid-base catalysis in the HDV genomic ribozyme

In the ribozyme from the hepatitis delta virus (HDV) genomic strand RNA, a cytosine side chain is proposed to facilitate proton transfer in the transition state of the reaction and, thus, act as a general acid-base catalyst. Mutation of this active-site cytosine (C75) reduced RNA cleavage rates by as much as one million-fold, but addition of exogenous cytosine and certain nucleobase or imidazole analogs can partially rescue activity in these mutants. However, pH-rate profiles for the rescued reactions were bell shaped, and only one leg of the pH-rate curve could be attributed to ionization of the exogenous nucleobase or buffer. When a second potential ionizable nucleobase (C41) was removed, one leg of the bell-shaped curve was eliminated in the chemical-rescue reaction. With this construct, the apparent pK(a) determined from the pH-rate profile correlated with the solution pK(a) of the buffer, and the contribution of the buffer to the rate enhancement could be directly evaluated in a free-energy or Brønsted plot. The free-energy relationship between the acid dissociation constant of the buffer and the rate constant for cleavage (Brønsted value, beta, = approximately 0.5) was consistent with a mechanism in which the buffer acted as a general acid-base catalyst. These data support the hypothesis that cytosine 75, in the intact ribozyme, acts as a general acid-base catalyst.

Water in the active site of an all-RNA hairpin ribozyme and effects of Gua8 base variants on the geometry of phosphoryl transfer

The hairpin ribozyme requires functional group contributions from G8 to assist in phosphodiester bond cleavage. Previously, replacement of G8 by a series of nucleobase variants showed little effect on interdomain docking, but a 3-250-fold effect on catalysis. To identify G8 features that contribute to catalysis within the hairpin ribozyme active site, structures for five base variants were determined by X-ray crystallography in a resolution range between 2.3 and 2.7 A. For comparison, a native all-RNA "G8" hairpin ribozyme structure was refined to 2.05 A resolution. The native structure revealed a scissile bond angle (tau) of 158 degrees, which is close to the requisite 180 degrees "in-line" geometry. Mutations G8(inosine), G8(diaminopurine), G8(aminopurine), G8(adenosine), and G8(uridine) folded properly, but exhibited nonideal scissile bond geometries (tau ranging from 118 degrees to 93 degrees) that paralleled their diminished solution activities. A superposition ensemble of all structures, including a previously described hairpin ribozyme-vanadate complex, indicated the scissile bond can adopt a variety of conformations resulting from perturbation of the chemical environment and provided a rationale for how the exocyclic amine of nucleobase 8 promotes productive, in-line geometry. Changes at position 8 also caused variations in the A-1 sugar pucker. In this regard, variants A8 and U8 appeared to represent nonproductive ground states in which their 2'-OH groups mimicked the pro-R, nonbridging oxygen of the vanadate transition-state complex. Finally, the results indicated that ordered water molecules bind near the 2'-hydroxyl of A-1, lending support to the hypothesis that solvent may play an important role in the reaction.

Nucleotide analogues to investigate RNA structure and function

RNA plays an essential cellular role in nearly every aspect of the transmission and expression of genetic information, including regulatory roles that have significance for cellular development. Access to RNA bearing synthetic modifications has allowed biological chemists to probe deep into the inner workings of cellular processes. Here, we describe recent advances in harnessing the power of nucleotide analogues to obtain mechanistic and biological insights into RNA structure, function and dynamics.

Structure, folding and mechanisms of ribozymes

The past two years have seen exciting developments in RNA catalysis. A completely new ribozyme (possibly two) has come along and several new structures have been determined, including three different group I intron species. Although the origins of catalysis remain incompletely understood, a significant convergence of views has happened in the past year, together with the discovery of new super-fast ribozymes. There is persuasive evidence of general acid-base chemistry in nucleolytic ribozymes, whereas catalysis of peptidyl transfer in the ribosome seems to result largely from orientation and proximity effects. Lastly, important new folding-enhancing elements have been discovered.

The catalytic diversity of RNAs

The natural RNA enzymes catalyse phosphate-group transfer and peptide-bond formation. Initially, metal ions were proposed to supply the chemical versatility that nucleotides lack. In the ensuing decades, structural and mechanistic studies have substantially altered this initial viewpoint. Whereas self-splicing ribozymes clearly rely on essential metal-ion cofactors, self-cleaving ribozymes seem to use nucleotide bases for their catalytic chemistry. Despite the overall differences in chemical features, both RNA and protein enzymes use similar catalytic strategies.

Ribozyme catalysis: not different, just worse

Evolution has resoundingly favored protein enzymes over RNA-based catalysts, yet ribozymes occupy important niches in modern cell biology that include the starring role in catalysis of protein synthesis on the ribosome. Recent results from structural and biochemical studies show that natural ribozymes use an impressive range of catalytic mechanisms, beyond metalloenzyme chemistry and analogous to more chemically diverse protein enzymes. These findings make it increasingly possible to compare details of RNA- and protein-based catalysis.

General acid catalysis by the hepatitis delta virus ribozyme

Recent crystallographic and functional analyses of RNA enzymes have raised the possibility that the purine and pyrimidine nucleobases may function as general acid-base catalysts. However, this mode of nucleobase-mediated catalysis has been difficult to establish unambiguously. Here, we used a hyperactivated RNA substrate bearing a 5'-phosphorothiolate to investigate the role of a critical cytosine residue in the hepatitis delta virus ribozyme. The hyperactivated substrate specifically suppressed the deleterious effects of cytosine mutations and pH changes, thereby linking the protonation of the nucleobase to leaving-group stabilization. We conclude that the active-site cytosine provides general acid catalysis, mediating proton transfer to the leaving group through a protonated N3-imino nitrogen. These results establish a specific role for a nucleobase in a ribozyme reaction and support the proposal that RNA nucleobases may function in a manner analogous to that of catalytic histidine residues in protein enzymes.

Role of an active site adenine in hairpin ribozyme catalysis

The hairpin ribozyme is a small catalytic RNA that accelerates reversible cleavage of a phosphodiester bond. Structural and mechanistic studies suggest that divalent metals stabilize the functional structure but do not participate directly in catalysis. Instead, two active site nucleobases, G8 and A38, appear to participate in catalytic chemistry. The features of A38 that are important for active site structure and chemistry were investigated by comparing cleavage and ligation reactions of ribozyme variants with A38 modifications. An abasic substitution of A38 reduced cleavage and ligation activity by 14,000-fold and 370,000-fold, respectively, highlighting the critical role of this nucleobase in ribozyme function. Cleavage and ligation activity of unmodified ribozymes increased with increasing pH, evidence that deprotonation of some functional group with an apparent pK(a) value near 6 is important for activity. The pH-dependent transition in activity shifted by several pH units in the basic direction when A38 was substituted with an abasic residue, or with nucleobase analogs with very high or low pK(a) values that are expected to retain the same protonation state throughout the experimental pH range. Certain exogenous nucleobases that share the amidine group of adenine restored activity to abasic ribozyme variants that lack A38. The pH dependence of chemical rescue reactions also changed according to the intrinsic basicity of the rescuing nucleobase, providing further evidence that the protonation state of the N1 position of purine analogs is important for rescue activity. These results are consistent with models of the hairpin ribozyme catalytic mechanism in which interactions with A38 provide electrostatic stabilization to the transition state.

Architecture of a Diels-Alderase ribozyme with a preformed catalytic pocket

Artificial ribozymes catalyze a variety of chemical reactions. Their structures and reaction mechanisms are largely unknown. We have analyzed a ribozyme catalyzing Diels-Alder cycloaddition reactions by comprehensive mutation analysis and a variety of probing techniques. New tertiary interactions involving base pairs between nucleotides of the 5' terminus and a large internal loop forming a pseudoknot fold were identified. The probing data indicate a preformed tertiary structure that shows no major changes on substrate or product binding. Based on these observations, a molecular architecture featuring a Y-shaped arrangement is proposed. The tertiary structure is formed in a rather unusual way; that is, the opposite sides of the asymmetric internal loop are clamped by the four 5'-terminal nucleotides, forming two adjacent two base-pair helices. It is proposed that the catalytic pocket is formed by a wedge within one of these helices.

Role of an active site guanine in hairpin ribozyme catalysis probed by exogenous nucleobase rescue

The hairpin ribozyme is a small catalytic RNA with reversible phosphodiester cleavage activity. Biochemical and structural studies exclude a requirement for divalent metal cation cofactors and implicate one active site nucleobase in particular, G8, in the catalytic mechanism. Our previous work demonstrated that the cleavage activity that is lost when G8 is replaced by an abasic residue is restored when certain nucleobases are provided in solution. The specificity and pH dependence of exogenous nucleobase rescue were consistent with several models of the rescue mechanism, including general acid base catalysis, electrostatic stabilization of negative charge in the transition state or a requirement for protonation to facilitate exogenous nucleobase binding. Detailed analyses of exogenous nucleobase rescue for both cleavage and ligation reactions now allow us to refine models of the rescue mechanism. Activity increased with increasing pH for both unmodified ribozyme reactions and unrescued reactions of abasic variants lacking G8. This similarity in pH dependence argues against a role for G8 as a general base catalyst, because G8 deprotonation could not be responsible for the pH-dependent transition in the abasic variant. Exogenous nucleobase rescue of both cleavage and ligation activity increased with decreasing pH, arguing against a role for rescuing nucleobases in general acid catalysis, because a nucleobase that contributes general acid catalysis in the cleavage pathway should provide general base catalysis in ligation. Analysis of the concentration dependence of cytosine rescue at high and low pH demonstrated that protonation promotes catalysis within the nucleobase-bound ribozyme complex but does not stabilize nucleobase binding in the ground state. These results support an electrostatic stabilization mechanism in which exogenous nucleobase binding counters negative charge that develops in the transition state.

Single-molecule enzymology of RNA: essential functional groups impact catalysis from a distance

The hairpin ribozyme is a minimalist paradigm for studying RNA folding and function. In this enzyme, two domains dock by induced fit to form a catalytic core that mediates a specific backbone cleavage reaction. Here, we have fully dissected its reversible reaction pathway, which comprises two structural transitions (docking/undocking) and a chemistry step (cleavage/ligation), by applying a combination of single-molecule fluorescence resonance energy transfer (FRET) assays, ensemble cleavage assays, and kinetic simulations. This has allowed us to quantify the effects that modifications of essential functional groups remote from the site of catalysis have on the individual rate constants. We find that all ribozyme variants show similar fractionations into effectively noninterchanging molecule subpopulations of distinct undocking rate constants. This leads to heterogeneous cleavage activity as commonly observed for RNA enzymes. A modification at the domain junction additionally leads to heterogeneous docking. Surprisingly, most modifications not only affect docking/undocking but also significantly impact the internal chemistry rate constants over a substantial distance from the site of catalysis. We propose that a network of coupled molecular motions connects distant parts of the RNA with its reaction site, which suggests a previously undescribed analogy between RNA and protein enzymes. Our findings also have broad implications for applications such as the action of drugs and ligands distal to the active site or the engineering of allostery into RNA.