RNA Electrostatics: How Ribozymes Engineer Active Sites to Enable Catalysis

Electrostatic interactions are fundamental to RNA structure and function, and intimately influenced by solvation and the ion atmosphere. RNA enzymes, or ribozymes, are catalytic RNAs that are able to enhance reaction rates over a million-fold, despite having only a limited repertoire of building blocks and available set of chemical functional groups. Ribozyme active sites usually occur at junctions where negatively charged helices come together, and in many cases leverage this strained electrostatic environment to recruit metal ions in solution that can assist in catalysis. Similar strategies have been implicated in related artificially engineered DNA enzymes. Herein, we apply Poisson-Boltzmann, 3D-RISM, and molecular simulations to study a set of metal-dependent small self-cleaving ribozymes (hammerhead, pistol, and Varkud satellite) as well as an artificially engineered DNAzyme (8-17) to examine electrostatic features and their relation to the recruitment of monovalent and divalent metal ions important for activity. We examine several fundamental roles for these ions that include: (1) structural integrity of the catalytically active state, (2) pKa tuning of residues involved in acid-base catalysis, and (3) direct electrostatic stabilization of the transition state via Lewis acid catalysis. Taken together, these examples demonstrate how RNA electrostatics orchestrates the site-specific and territorial binding of metal ions to play important roles in catalysis.

Evidence for a Catalytic Strategy to Promote Nucleophile Activation in Metal-Dependent RNA-Cleaving Ribozymes and 8-17 DNAzyme

An unique catalytic strategy was recently reported for the glmS ribozyme [Bingaman et al., Nat. Chem. Biol.2017, 13, 439-445] that involves promotion of productive hydrogen bonding of the O2' nucleophile to facilitate its activation. We provide broad evidence of this strategy in the hammerhead, pistol, and VS ribozymes and 8-17 DNAzyme, enabled by a functionally important divalent metal ion that interacts with the scissile phosphate and disrupts nonproductive competitive hydrogen bonding with the O2' nucleophile. This strategy, designated tertiary gamma (3°γ) catalysis, illustrates an additional role for divalent ions in ribozyme catalysis.

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.

Mode of action of RNase BN/RNase Z on tRNA precursors: RNase BN does not remove the CCA sequence from tRNA

RNase BN, the Escherichia coli homolog of RNase Z, was previously shown to act as both a distributive exoribonuclease and an endoribonuclease on model RNA substrates and to be inhibited by the presence of a 3'-terminal CCA sequence. Here, we examined the mode of action of RNase BN on bacteriophage and bacterial tRNA precursors, particularly in light of a recent report suggesting that RNase BN removes CCA sequences (Takaku, H., and Nashimoto, M. (2008) Genes Cells 13, 1087-1097). We show that purified RNase BN can process both CCA-less and CCA-containing tRNA precursors. On CCA-less precursors, RNase BN cleaved endonucleolytically after the discriminator nucleotide to allow subsequent CCA addition. On CCA-containing precursors, RNase BN acted as either an exoribonuclease or endoribonuclease depending on the nature of the added divalent cation. Addition of Co(2+) resulted in higher activity and predominantly exoribonucleolytic activity, whereas in the presence of Mg(2+), RNase BN was primarily an endoribonuclease. In no case was any evidence obtained for removal of the CCA sequence. Certain tRNA precursors were extremely poor substrates under any conditions tested. These findings provide important information on the ability of RNase BN to process tRNA precursors and help explain the known physiological properties of this enzyme. In addition, they call into question the removal of CCA sequences by RNase BN.

Formation of an active site in trans by interaction of two complete Varkud Satellite ribozymes

The complete VS ribozyme comprises seven helical segments, connected by three three-way RNA junctions. In the presence of Mg(2+) ions, cleavage occurs within the internal loop of helix I. This requires the participation of a guanine (G638) within the helix I loop, and a remote adenine (A756) within an internal loop of helix VI. Previous structural studies have suggested that helix I docks into the fold of the remaining part of the ribozyme, bringing A756 and G638 close to the scissile phosphate to allow the cleavage reaction to proceed. We show here that while either A756C or G638A individually exhibit very low cleavage activity, a mixture of the two variants leads to cleavage of the A756C RNA, but not the G638A RNA. The rate of cleavage depends on the concentration of the VS G638A RNA, as expected for a bimolecular interaction. This regaining of cleavage activity by complementation indicates that helix I of one VS RNA can interact with another VS RNA molecule to generate a functional active site in trans.

Single-Molecule Fluorescence Studies of RNA: A Decade's Progress

Over the past decade, single-molecule fluorescence studies have elucidated the structure-function relationship of RNA molecules. The real-time observation of individual RNAs by single-molecule fluorescence has unveiled the dynamic behavior of complex RNA systems in unprecedented detail, revealing the presence of transient intermediate states and their kinetic pathways. This review provides an overview of how single-molecule fluorescence has been used to explore the dynamics of RNA folding and catalysis.

Identification and characterization of a short 2'-3' bond-forming ribozyme

A large number of natural and artificial ribozymes have been isolated since the demonstration of the catalytic potential of RNA, with the majority of these catalyzing phosphate hydrolysis or transesterification reactions. Here, we describe and characterize an extremely short ribozyme that catalyzes the positionally specific transesterification that produces a 2'-3' phosphodiester bond between itself and a branch substrate provided in trans, cleaving itself internally in the process. Although this ribozyme was originally derived from constructs based on snRNAs, its minimal catalytic motif contains essentially no snRNA sequence and the reaction it catalyzes is not directly related to either step of pre-mRNA splicing. Our data have implications for the intrinsic reactivity of the large amount of RNA sequence space known to be transcribed in nature and for the validity and utility of the use of protein-free systems to study pre-mRNA splicing.

A critical assessment of the utility of protein-free splicing systems

U2 and U6 snRNAs form part of the catalytic spliceosome and represent strong candidates for components of its active site. Over the past decade it has become clear that these snRNAs are capable of catalyzing several different chemical reactions, leading to the widespread conclusion that the spliceosome is a ribozyme. Here, we discuss the advances in both protein-free and fully spliceosomal systems that would be required to conclude that the reactions observed to be catalyzed by protein-free snRNAs are related to splicing and question the reliability of snRNA-only systems as tools for mechanistic splicing research.

Metal ion specificities for folding and cleavage activity in the Schistosoma hammerhead ribozyme

The effects of various metal ions on cleavage activity and global folding have been studied in the extended Schistosoma hammerhead ribozyme. Fluorescence resonance energy transfer was used to probe global folding as a function of various monovalent and divalent metal ions in this ribozyme. The divalent metals ions Ca(2+), Mg(2+), Mn(2+), and Sr(2+) have a relatively small variation (less than sixfold) in their ability to globally fold the hammerhead ribozyme, which contrasts with the very large difference (>10,000-fold) in apparent rate constants for cleavage for these divalent metal ions in single-turnover kinetic experiments. There is still a very large range (>4600-fold) in the apparent rate constants for cleavage for these divalent metal ions measured in high salt (2 M NaCl) conditions where the ribozyme is globally folded. These results demonstrate that the identity of the divalent metal ion has little effect on global folding of the Schistosoma hammerhead ribozyme, whereas it has a very large effect on the cleavage kinetics. Mechanisms by which the identity of the divalent metal ion can have such a large effect on cleavage activity in the Schistosoma hammerhead ribozyme are discussed.

The complete VS ribozyme in solution studied by small-angle X-ray scattering

We have used small-angle X-ray solution scattering to obtain ab initio shape reconstructions of the complete VS ribozyme. The ribozyme occupies an electron density envelope with an irregular shape, into which helical sections have been fitted. The ribozyme is built around a core comprising a near-coaxial stack of three helices, organized by two three-way helical junctions. An additional three-way junction formed by an auxiliary helix directs the substrate stem-loop, juxtaposing the cleavage site with an internal loop to create the active complex. This is consistent with the current view of the probable mechanism of trans-esterification in which adenine and guanine nucleobases contributed by the interacting loops combine in general acid-base catalysis.

A glimpse into the active site of a group II intron and maybe the spliceosome, too

The X-ray crystal structure of an excised group II self-splicing intron was recently solved by the Pyle group. Here we review some of the notable features of this structure and what they may tell us about the catalytic active site of the group II ribozyme and potentially the spliceosome. The new structure validates the central role of domain V in both the structure and catalytic function of the ribozyme and resolves several outstanding puzzles raised by previous biochemical, genetic and structural studies. While lacking both exons as well as the cleavage sites and nucleophiles, the structure reveals how a network of tertiary interactions can position two divalent metal ions in a configuration that is ideal for catalysis.

Nucleotide synthetase ribozymes may have emerged first in the RNA world

Though the "RNA world" hypothesis has gained a central role in ideas concerning the origin of life, the scenario concerning its emergence remains uncertain. It has been speculated that the first scene may have been the emergence of a template-dependent RNA synthetase ribozyme, which catalyzed its own replication: thus, "RNA replicase." However, the speculation remains uncertain, primarily because of the large sequence length requirement of such a replicase and the lack of a convincing mechanism to ensure its self-favoring features. Instead, we propose a nucleotide synthetase ribozyme as an alternative candidate, especially considering recent experimental evidence suggesting the possibility of effective nonenzymatic template-directed synthesis of RNA. A computer simulation was conducted to support our proposal. The conditions for the emergence of the nucleotide synthetase ribozyme are discussed, based on dynamic analysis on a computer. We suggest the template-dependent RNA synthetase ribozyme emerged later, perhaps after the emergence of protocells.

Leadzyme formed in vivo interferes with tobacco mosaic virus infection in Nicotiana tabacum

We developed a new method for inhibiting tobacco mosaic virus infection in tobacco plants based on specific RNA hydrolysis induced by a leadzyme. We identified a leadzyme substrate target sequence in genomic tobacco mosaic virus RNA and designed a 16-mer oligoribonucleotide capable of forming a specific leadzyme motif with a five-nucleotide catalytic loop. The synthetic 16-mer RNA was applied with nontoxic, catalytic amount of lead to infected tobacco leaves. We observed inhibition of tobacco mosaic virus infection in tobacco leaves in vivo due to specific tobacco mosaic virus RNA cleavage effected by leadzyme. A significant reduction in tobacco mosaic virus accumulation was observed even when the leadzyme was applied up to 2 h after inoculation of leaves with tobacco mosaic virus. This process, called leadzyme interference, is determined by specific recognition and cleavage of the target site by the RNA catalytic strand in the presence of Pb(2+).

Nucleobase catalysis in the hairpin ribozyme

RNA catalysis is important in the processing and translation of RNA molecules, yet the mechanisms of catalysis are still unclear in most cases. We have studied the role of nucleobase catalysis in the hairpin ribozyme, where the scissile phosphate is juxtaposed between guanine and adenine bases. We show that a modified ribozyme in which guanine 8 has been substituted by an imidazole base is active in both cleavage and ligation, with ligation rates 10-fold faster than cleavage. The rates of both reactions exhibit bell-shaped dependence on pH, with pK(a) values of 5.7 +/- 0.1 and 7.7 +/- 0.1 for cleavage and 6.1 +/- 0.3 and 6.9 +/- 0.3 for ligation. The data provide good evidence for general acid-base catalysis by the nucleobases.

Structure of a self-splicing group II intron catalytic effector domain 5: parallels with spliceosomal U6 RNA

Domain 5 (D5) is absolutely required for all catalytic functions of group II introns. Here we describe the solution NMR structure, electrostatic calculations, and detailed magnesium ion-binding surface of D5 RNA from the Pylaiella littoralis large ribosomal RNA intron (D5-PL). The overall structure consists of a hairpin capped by a GNRA tetraloop. The stem is divided into lower and upper helices of 8 and 5 bp, respectively, separated by an internal bulge. The D5-PL internal bulge nucleotides stack into the helical junction, resulting in a coupling between the bulge A25 and the closing base pair (G8-C27) of the lower helix. Comparison of the D5-PL structure to previously reported related structures indicates that our structure is most similar, in the helical regions, to the crystal structure of D5 from yeast Ai5gamma (D5-Ai5gamma) and the NMR structure of the U6 snRNA stem-loop region. Our structure differs in many respects from both the NMR and X-ray structures of D5-Ai5gamma in the bulge region. Electrostatic calculations and NMR chemical shift perturbation analyses reveal magnesium ion-binding sites in the tetraloop, internal bulge, and the AGC triad in the lower stem. Our results suggest that the structure, electrostatic environment, and the magnesium ion-binding sites within the tetraloop, bulge, and triad regions are conserved features of the splicing machinery of both the group II introns and the spliceosome that are likely key for catalytic function.

Characterization of a native hammerhead ribozyme derived from schistosomes

A recent re-examination of the role of the helices surrounding the conserved core of the hammerhead ribozyme has identified putative loop-loop interactions between stems I and II in native hammerhead sequences. These extended hammerhead sequences are more active at low concentrations of divalent cations than are minimal hammerheads. The loop-loop interactions are proposed to stabilize a more active conformation of the conserved core. Here, a kinetic and thermodynamic characterization of an extended hammerhead sequence derived from Schistosoma mansoni is performed. Biphasic kinetics are observed, suggesting the presence of at least two conformers, one cleaving with a fast rate and the other with a slow rate. Replacing loop II with a poly(U) sequence designed to eliminate the interaction between the two loops results in greatly diminished activity, suggesting that the loop-loop interactions do aid in forming a more active conformation. Previous studies with minimal hammerheads have shown deleterious effects of Rp-phosphorothioate substitutions at the cleavage site and 5' to A9, both of which could be rescued with Cd2+. Here, phosphorothioate modifications at the cleavage site and 5' to A9 were made in the schistosome-derived sequence. In Mg2+, both phosphorothioate substitutions decreased the overall fraction cleaved without significantly affecting the observed rate of cleavage. The addition of Cd2+ rescued cleavage in both cases, suggesting that these are still putative metal binding sites in this native sequence.

Linkage between proton binding and folding in RNA: a thermodynamic framework and its experimental application for investigating pKa shifting

Perturbation of pKa values can change the favored protonation states of the nucleobases at biological pH and thereby modulate the function of RNA and DNA molecules. In an effort to understand the driving forces for pKa shifting specific to nucleic acids, we developed a thermodynamic framework that relates proton binding to the nucleobases and the helix-coil transition. Key features that emerge from the treatment are a comprehensive description of all the actions of proton binding on RNA folding: acid and alkaline denaturation of the helix and pKa shifting in the folded state. Practical experimental approaches for measuring pKas from thermal denaturation experiments are developed. Microscopic pka values (where ka is the acid dissociation constant) for the unfolded state were determined directly by experiments on unstructured oligonucleotides, which led to a macroscopic pKa for the ensemble of unfolded states shifted toward neutrality. The formalism was then applied to pH-dependent UV melting data for model DNA oligonucleotides. Folded-state pka) values were in good agreement with the outcome of pH titrations, and the acid and alkaline denaturation regions were well described. The formalism developed here is similar to that of Draper and coworkers for Mg2+ binding to RNA, except that the unfolded state is described explicitly owing to the presence of specific proton-binding sites on the bases. A principal conclusion is that it should be possible to attain large pKa shifts by designing RNA molecules that fold cooperatively.

Kinetic analysis of ribozyme-substrate complex formation in yeast

Many RNA-mediated reactions in transcription, translation, RNA processing, and transport require assembly of RNA complexes, yet assembly pathways remain poorly understood. Assembly mechanisms can be difficult to assess in a biological context because many components interact in complex pathways and individual steps are difficult to isolate experimentally. Our previous studies of self-cleaving hairpin ribozymes showed that kinetic and equilibrium parameters measured in yeast agree well with parameters measured in vitro under ionic conditions that mimic the intracellular environment. We now report studies of intermolecular reactions with ribozyme and target sequences expressed in yeast as separate chimeric U3 snoRNAs. In this system, intracellular cleavage rates reflect the kinetics of ribozyme-substrate complex formation through annealing of base-paired helices. Second-order rate constants increased with increasing helix length for in vitro reactions with 2 mM MgCl(2) and 150 mM NaCl and in vivo but not in reactions with 10 mM MgCl(2). Thus, efficient RNA complex formation required a larger extent of complementarity in vivo than in vitro under conditions with high concentrations of divalent cations. The most efficient intracellular cleavage reactions exhibited second-order rate constants that were 15- to 30-fold below rate constants for cleavage of oligonucleotides in vitro. Careful analysis of structural features that influence cleavage efficiency points to substrate binding as the rate-determining step in the intracellular cleavage pathway. Second-order rate constants for intermolecular cleavage agree well with diffusion coefficients reported for U3 snoRNPs in vivo suggesting that complex formation between chimeric ribozyme and substrate snoRNPs in yeast nuclei is diffusion limited.

Folding of the natural hammerhead ribozyme is enhanced by interaction of auxiliary elements

It has been shown that the activity of the hammerhead ribozyme at microM magnesium ion concentrations is markedly increased by the inclusion of loops in helices I and II. We have studied the effect of such loops on the magnesium ion-induced folding of the ribozyme, using fluorescence resonance energy transfer. We find that with the loops in place, folding into the active conformation occurs in a single step, in the microM range of magnesium ion concentration. Disruption of the loop-loop interaction leads to a reversion to two-step folding, with the second stage requiring mM concentrations of magnesium ion. Sodium ions also promote the folding of the natural form of the ribozyme at high concentrations, but the folding occurs as a two-stage process. The loops clearly act as important auxiliary elements in the function of the ribozyme, permitting folding to occur efficiently under physiological conditions.

The three-dimensional architecture of the class I ligase ribozyme

The class I ligase ribozyme catalyzes a Mg(++)-dependent RNA-ligation reaction that is chemically analogous to a single step of RNA polymerization. Indeed, this ribozyme constitutes the catalytic domain of an accurate and general RNA polymerase ribozyme. The ligation reaction is also very rapid in both single- and multiple-turnover contexts and thus is informative for the study of RNA catalysis as well as RNA self-replication. Here we report the initial characterization of the three-dimensional architecture of the ligase. When the ligase folds, several segments become protected from hydroxyl-radical cleavage, indicating that the RNA adopts a compact tertiary structure. Ribozyme folding was largely, though not completely, Mg(++) dependent, with a K(1/2[Mg]) < 1 mM, and was observed over a broad temperature range (20 degrees C -50 degrees C). The hydroxyl-radical mapping, together with comparative sequence analyses and analogy to a region within 23S ribosomal RNA, were used to generate a three-dimensional model of the ribozyme. The predictive value of the model was tested and supported by a photo-cross-linking experiment.

The Varkud satellite ribozyme

The VS ribozyme is the largest nucleolytic ribozyme, for which there is no crystal structure to date. The ribozyme consists of five helical sections, organized by two three-way junctions. The global structure has been determined by solution methods, particularly FRET. The substrate stem-loop binds into a cleft formed between two helices, while making a loop-loop contact with another section of the ribozyme. The scissile phosphate makes a close contact with an internal loop (the A730 loop), the probable active site of the ribozyme. This loop contains a particularly critical nucleotide A756. Most changes to this nucleotide lead to three-orders of magnitude slower cleavage, and the Watson-Crick edge is especially important. NAIM experiments indicate that a protonated base is required at this position for the ligation reaction. A756 is thus a strong candidate for nucleobase participation in the catalytic chemistry.

Dissection of the ion-induced folding of the hammerhead ribozyme using 19F NMR

We have used (19)F NMR to analyze the metal ion-induced folding of the hammerhead ribozyme by selective incorporation of 5fluorouridine. We have studied the chemical shift and linewidths of (19)F resonances of 5-fluorouridine at the 4 and 7 positions in the ribozyme core as a function of added Mg(2+). The data fit well to a simple two-state model whereby the formation of domain 1 is induced by the noncooperative binding of Mg(2+) with an association constant in the range of 100 to 500 M(-1), depending on the concentration of monovalent ions present. The results are in excellent agreement with data reporting on changes in the global shape of the ribozyme. However, the NMR experiments exploit reporters located in the center of the RNA sections undergoing the folding transitions, thereby allowing the assignment of specific nucleotides to the separate stages. The results define the folding pathway at high resolution and provide a time scale for the first transition in the millisecond range.