A new reagent for in vivo structure probing of RNA G and U residues that improves RNA structure prediction alone and combined with DMS.

A key to understanding the roles of RNA in regulating gene expression is knowing their structures in vivo. One way to obtain this information is through probing the structures of RNA with chemicals. To probe RNA structure directly in cells, membrane-permeable reagents that modify the Watson-Crick (WC) face of unpaired nucleotides can be used. Although dimethyl sulfate (DMS) has led to substantial insight into RNA structure, it has limited nucleotide specificity in vivo, with WC face reactivity only at adenine (A) and cytosine (C) at neutral pH. The reagent 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide (EDC) was recently shown to modify the WC face of guanine (G) and uracil (U). Although useful at lower concentrations in experiments that measure chemical modifications by reverse transcription (RT) stops, at higher concentrations necessary for detection by mutational profiling (MaP), EDC treatment leads to degradation of RNA. Here, we demonstrate EDC-stimulated degradation of RNA in Gram-negative and Gram-positive bacteria. In an attempt to overcome these limitations, we developed a new carbodiimide reagent, 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide methiodide (ETC), which we show specifically modifies unpaired Gs and Us in vivo without substantial degradation of RNA. We establish ETC as a probe for MaP and optimize the RT conditions and computational analysis in Escherichia coli Importantly, we demonstrate the utility of ETC as a probe for improving RNA structure prediction both alone and with DMS.

Rock, scissors, paper: How RNA structure informs function.

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

Transcriptome-wide probing reveals RNA thermometers that regulate translation of glycerol permease genes in Bacillus subtilis.

RNA structure regulates bacterial gene expression by several distinct mechanisms via environmental and cellular stimuli, one of which is temperature. While some genome-wide studies have focused on heat shock treatments and the subsequent transcriptomic changes, soil bacteria are less likely to experience such rapid and extreme temperature changes. Though RNA thermometers (RNATs) have been found in 5' untranslated leader regions (5' UTRs) of heat shock and virulence-associated genes, this RNA-controlled mechanism could regulate other genes as well. Using Structure-seq2 and the chemical probe dimethyl sulfate (DMS) at four growth temperatures ranging from 23°C to 42°C, we captured a dynamic response of the Bacillus subtilis transcriptome to temperature. Our transcriptome-wide results show RNA structural changes across all four temperatures and reveal nonmonotonic reactivity trends with increasing temperature. Then, focusing on subregions likely to contain regulatory RNAs, we examined 5' UTRs to identify large, local reactivity changes. This approach led to the discovery of RNATs that control the expression of glpF (glycerol permease) and glpT (glycerol-3-phosphate permease); expression of both genes increased with increased temperature. Results with mutant RNATs indicate that both genes are controlled at the translational level. Increased import of glycerols at high temperatures could provide thermoprotection to proteins.

In vivo structure probing of RNA in Archaea: novel insights into the ribosome structure of Methanosarcina acetivorans.

Structure probing combined with next-generation sequencing (NGS) has provided novel insights into RNA structure-function relationships. To date, such studies have focused largely on bacteria and eukaryotes, with little attention given to the third domain of life, archaea. Furthermore, functional RNAs have not been extensively studied in archaea, leaving open questions about RNA structure and function within this domain of life. With archaeal species being diverse and having many similarities to both bacteria and eukaryotes, the archaea domain has the potential to be an evolutionary bridge. In this study, we introduce a method for probing RNA structure in vivo in the archaea domain of life. We investigated the structure of ribosomal RNA (rRNA) from Methanosarcina acetivorans, a well-studied anaerobic archaeal species, grown with either methanol or acetate. After probing the RNA in vivo with dimethyl sulfate (DMS), Structure-seq2 libraries were generated, sequenced, and analyzed. We mapped the reactivity of DMS onto the secondary structure of the ribosome, which we determined independently with comparative analysis, and confirmed the accuracy of DMS probing in M. acetivorans Accessibility of the rRNA to DMS in the two carbon sources was found to be quite similar, although some differences were found. Overall, this study establishes the Structure-seq2 pipeline in the archaea domain of life and informs about ribosomal structure within M. acetivorans.

In vivo-like nearest neighbor parameters improve prediction of fractional RNA base-pairing in cells.

We conducted a thermodynamic analysis of RNA stability in Eco80 artificial cytoplasm, which mimics in vivo conditions, and compared it to transcriptome-wide probing of mRNA. Eco80 contains 80% of Escherichia coli metabolites, with biological concentrations of metal ions, including 2 mM free Mg2+ and 29 mM metabolite-chelated Mg2+. Fluorescence-detected binding isotherms (FDBI) were used to conduct a thermodynamic analysis of 24 RNA helices and found that these helices, which have an average stability of -12.3 kcal/mol, are less stable by ΔΔGo37 ∼1 kcal/mol. The FDBI data was used to determine a set of Watson-Crick free energy nearest neighbor parameters (NNPs), which revealed that Eco80 reduces the stability of three NNPs. This information was used to adjust the NN model using the RNAstructure package. The in vivo-like adjustments have minimal effects on the prediction of RNA secondary structures determined in vitro and in silico, but markedly improve prediction of fractional RNA base pairing in E. coli, as benchmarked with our in vivo DMS and EDC RNA chemical probing data. In summary, our thermodynamic and chemical probing analyses of RNA helices indicate that RNA secondary structures are less stable in cells than in artificially stable in vitro buffer conditions.

Genome-wide analysis of the in vivo tRNA structurome reveals RNA structural and modification dynamics under heat stress.

RNA structure plays roles in myriad cellular events including transcription, translation, and RNA processing. Genome-wide analyses of RNA secondary structure in vivo by chemical probing have revealed critical structural features of mRNAs and long ncRNAs. Here, we examine the in vivo secondary structure of a small RNA class, tRNAs. Study of tRNA structure is challenging because tRNAs are heavily modified and strongly structured. We introduce "tRNA structure-seq," a new workflow that accurately determines in vivo secondary structures of tRNA. The workflow combines dimethyl sulfate (DMS) probing, ultra-processive RT, and mutational profiling (MaP), which provides mutations opposite DMS and natural modifications thereby allowing multiple modifications to be identified in a single read. We applied tRNA structure-seq to E. coli under control and stress conditions. A leading folding algorithm predicts E. coli tRNA structures with only ∼80% average accuracy from sequence alone. Strikingly, tRNA structure-seq, by providing experimental restraints, improves structure prediction under in vivo conditions to ∼95% accuracy, with more than 14 tRNAs predicted completely correctly. tRNA structure-seq also quantifies the relative levels of tRNAs and their natural modifications at single nucleotide resolution, as validated by LC-MS/MS. Our application of tRNA structure-seq yields insights into tRNA structure in living cells, revealing that it is not immutable but has dynamics, with partial unfolding of secondary and tertiary tRNA structure under heat stress that is correlated with a loss of tRNA abundance. This method is applicable to other small RNAs, including those with natural modifications and highly structured regions.

Flanking Sequence Cotranscriptionally Regulates Twister Ribozyme Activity.

Small nucleolytic ribozymes are RNAs that cleave their own phosphodiester backbone. While proteinaceous enzymes are regulated by a variety of known mechanisms, methods of regulation for ribozymes remain unclear. Twister is one ribozyme class for which many structural and catalytic properties have been elucidated. However, few studies have analyzed the activity of twister ribozymes in the context of a native flanking sequence, even though ribozymes as transcribed in nature do not exist in isolation. Interactions between the ribozyme and its neighboring sequences can induce conformational changes that inhibit self-cleavage, providing a regulatory mechanism that could naturally determine ribozyme activity in vivo and in synthetic applications. To date, eight twister ribozymes have been identified within the staple crop rice (Oryza sativa). Herein, we select several twister ribozymes from rice and show that they are differentially regulated by their flanking sequence using published RNA-seq data sets, structure probing, and cotranscriptional cleavage assays. We found that the Osa 1-2 ribozyme does not interact with its flanking sequences. However, sequences flanking the Osa 1-3 and Osa 1-8 ribozymes form inactive conformations, referred to here as "ribozymogens", that attenuate ribozyme self-cleavage activity. For the Osa 1-3 ribozyme, we show that activity can be rescued upon addition of a complementary antisense oligonucleotide, suggesting ribozymogens can be controlled via external signals. In all, our data provide a plausible mechanism wherein flanking sequence differentially regulates ribozyme activity in vivo. More broadly, the ability to regulate ribozyme behavior locally has potential applications in control of gene expression and synthetic biology.

Genome-Wide Analysis of RNA Secondary Structure.

Single-stranded RNA molecules fold into extraordinarily complicated secondary and tertiary structures as a result of intramolecular base pairing. In vivo, these RNA structures are not static. Instead, they are remodeled in response to changes in the prevailing physicochemical environment of the cell and as a result of intermolecular base pairing and interactions with RNA-binding proteins. Remarkable technical advances now allow us to probe RNA secondary structure at single-nucleotide resolution and genome-wide, both in vitro and in vivo. These data sets provide new glimpses into the RNA universe. Analyses of RNA structuromes in HIV, yeast, Arabidopsis, and mammalian cells and tissues have revealed regulatory effects of RNA structure on messenger RNA (mRNA) polyadenylation, splicing, translation, and turnover. Application of new methods for genome-wide identification of mRNA modifications, particularly methylation and pseudouridylation, has shown that the RNA "epitranscriptome" both influences and is influenced by RNA structure. In this review, we describe newly developed genome-wide RNA structure-probing methods and synthesize the information emerging from their application.

Biological solution conditions and flanking sequence modulate LLPS of RNA G-quadruplex structures.

Guanine-rich regions of DNA or RNA can form structures with two or more consecutive G-quartets called G-quadruplexes (GQ). Recent studies reveal the potential for these structures to aggregate in vitro. Here, we report effects of in vivo concentrations of additives-amino acids, nucleotides, and crowding agents-on the structure and solution behavior of RNAs containing GQ-forming sequences. We found that cytosine nucleotides destabilize a model GQ structure at biological salt concentrations, while free amino acids and other nucleotides do not do so to a substantial degree. We also report that the tendency of folded GQs to form droplets or to aggregate depends on the nature of flanking sequence and the presence of additives. Notably, in the presence of biological amounts of polyamines, flanking regions on the 5'-end of the RNA drive more droplet-like phase separation, while flanking regions on the 3'-end, as well as both the 5'- and 3'-ends, induce more condensed, granular structures. Finally, we provide an example of a biological sequence in the presence of polyamines and show that crowders such as PEG and dextran can selectively cause its phase separation. These findings have implications for the participation of GQS in LLPS in vivo.

RNA multimerization as an organizing force for liquid-liquid phase separation.

RNA interactions are exceptionally strong and highly redundant. As such, nearly any two RNAs have the potential to interact with one another over relatively short stretches, especially at high RNA concentrations. This is especially true for pairs of RNAs that do not form strong self-structure. Such phenomena can drive liquid-liquid phase separation, either solely from RNA-RNA interactions in the presence of divalent or organic cations, or in concert with proteins. RNA interactions can drive multimerization of RNA strands via both base-pairing and tertiary interactions. In this article, we explore the tendency of RNA to form stable monomers, dimers, and higher order structures as a function of RNA length and sequence through a focus on the intrinsic thermodynamic, kinetic, and structural properties of RNA. The principles we discuss are independent of any specific type of biomolecular condensate, and thus widely applicable. We also speculate how external conditions experienced by living organisms can influence the formation of nonmembranous compartments, again focusing on the physical and structural properties of RNA. Plants, in particular, are subject to diverse abiotic stresses including extreme temperatures, drought, and salinity. These stresses and the cellular responses to them, including changes in the concentrations of small molecules such as polyamines, salts, and compatible solutes, have the potential to regulate condensate formation by melting or strengthening base-pairing. Reversible condensate formation, perhaps including regulation by circadian rhythms, could impact biological processes in plants, and other organisms.

RNA sequence and structure control assembly and function of RNA condensates.

Intracellular condensates formed through liquid-liquid phase separation (LLPS) primarily contain proteins and RNA. Recent evidence points to major contributions of RNA self-assembly in the formation of intracellular condensates. As the majority of previous studies on LLPS have focused on protein biochemistry, effects of biological RNAs on LLPS remain largely unexplored. In this study, we investigate the effects of crowding, metal ions, and RNA structure on formation of RNA condensates lacking proteins. Using bacterial riboswitches as a model system, we first demonstrate that LLPS of RNA is promoted by molecular crowding, as evidenced by formation of RNA droplets in the presence of polyethylene glycol (PEG 8K). Crowders are not essential for LLPS, however. Elevated Mg2+ concentrations promote LLPS of specific riboswitches without PEG. Calculations identify key RNA structural and sequence elements that potentiate the formation of PEG-free condensates; these calculations are corroborated by key wet-bench experiments. Based on this, we implement structure-guided design to generate condensates with novel functions including ligand binding. Finally, we show that RNA condensates help protect their RNA components from degradation by nucleases, suggesting potential biological roles for such higher-order RNA assemblies in controlling gene expression through RNA stability. By utilizing both natural and artificial RNAs, our study provides mechanistic insight into the contributions of intrinsic RNA properties and extrinsic environmental conditions to the formation and regulation of condensates comprised of RNAs.

The Metabolome Weakens RNA Thermodynamic Stability and Strengthens RNA Chemical Stability.

We examined the complex network of interactions among RNA, the metabolome, and divalent Mg2+ under conditions that mimic the Escherichia coli cytoplasm. We determined Mg2+ binding constants for the top 15 E. coli metabolites, comprising 80% of the metabolome by concentration at physiological pH and monovalent ion concentrations. These data were used to inform the development of an artificial cytoplasm that mimics in vivo E. coli conditions, which we term "Eco80". We empirically determined that the mixture of E. coli metabolites in Eco80 approximated single-site binding behavior toward Mg2+ in the biologically relevant free Mg2+ range of ∼0.5 to 3 mM Mg2+, using a Mg2+-sensitive fluorescent dye. Effects of Eco80 conditions on the thermodynamic stability, chemical stability, structure, and catalysis of RNA were examined. We found that Eco80 conditions lead to opposing effects on the thermodynamic and chemical stabilities of RNA. In particular, the thermodynamic stability of RNA helices was weakened by 0.69 ± 0.12 kcal/mol, while the chemical stability was enhanced ∼2-fold, which can be understood using the speciation of Mg2+ between weak and strong Mg2+-metabolite complexes in Eco80. Overall, the use of Eco80 reflects RNA function in vivo and enhances the biological relevance of mechanistic studies of RNA.

Tissue-specific changes in the RNA structurome mediate salinity response in Arabidopsis.

Little is known concerning the effects of abiotic factors on in vivo RNA structures. We applied Structure-seq to assess the in vivo mRNA structuromes of Arabidopsis thaliana under salinity stress, which negatively impacts agriculture. Structure-seq utilizes dimethyl sulfate reactivity to identify As and Cs that lack base-pairing or protection. Salt stress refolded transcripts differentially in root versus shoot, evincing tissue specificity of the structurome. Both tissues exhibited an inverse correlation between salt stress-induced changes in transcript reactivity and changes in abundance, with stress-related mRNAs showing particular structural dynamism. This inverse correlation is more pronounced in mRNAs wherein the mean reactivity of the 5'UTR, CDS, and 3'UTR concertedly change under salinity stress, suggesting increased susceptibility to abundance control mechanisms in transcripts exhibiting this phenomenon, which we name "concordancy." Concordant salinity-induced increases in reactivity were notably observed in photosynthesis genes, thereby implicating mRNA structural loss in the well-known depression of photosynthesis by salt stress. Overall, changes in secondary structure appear to impact mRNA abundance, molding the functional specificity of the transcriptome under stress.

Structure-seq2 probing of RNA structure upon amino acid starvation reveals both known and novel RNA switches in Bacillus subtilis.

RNA structure influences numerous processes in all organisms. In bacteria, these processes include transcription termination and attenuation, small RNA and protein binding, translation initiation, and mRNA stability, and can be regulated via metabolite availability and other stresses. Here we use Structure-seq2 to probe the in vivo RNA structurome of Bacillus subtilis grown in the presence and absence of amino acids. Our results reveal that amino acid starvation results in lower overall dimethyl sulfate (DMS) reactivity of the transcriptome, indicating enhanced protection owing to protein binding or RNA structure. Starvation-induced changes in DMS reactivity correlated inversely with transcript abundance changes. This correlation was particularly pronounced in genes associated with the stringent response and CodY regulons, which are involved in adaptation to nutritional stress, suggesting that RNA structure contributes to transcript abundance change in regulons involved in amino acid metabolism. Structure-seq2 accurately reported on four known amino acid-responsive riboswitches: T-box, SAM, glycine, and lysine riboswitches. Additionally, we discovered a transcription attenuation mechanism that reduces yfmG expression when amino acids are added to the growth medium. We also found that translation of a leader peptide (YfmH) encoded just upstream of yfmG regulates yfmG expression. Our results are consistent with a model in which a slow rate of yfmH translation caused by limitation of the amino acids encoded in YfmH prevents transcription termination in the yfmG leader region by favoring formation of an overlapping antiterminator structure. This novel RNA switch offers a way to simultaneously monitor the levels of multiple amino acids.

Single-nucleotide control of tRNA folding cooperativity under near-cellular conditions.

RNA folding is often studied by renaturing full-length RNA in vitro and tracking folding transitions. However, the intracellular transcript folds as it emerges from the RNA polymerase. Here, we investigate the folding pathways and stability of numerous late-transcriptional intermediates of yeast and Escherichia coli transfer RNAs (tRNAs). Transfer RNA is a highly regulated functional RNA that undergoes multiple steps of posttranscriptional processing and is found in very different lengths during its lifetime in the cell. The precursor transcript is extended on both the 5' and 3' ends of the cloverleaf core, and these extensions get trimmed before addition of the 3'-CCA and aminoacylation. We studied the thermodynamics and structures of the precursor tRNA and of late-transcriptional intermediates of the cloverleaf structure. We examined RNA folding at both the secondary and tertiary structural levels using multiple biochemical and biophysical approaches. Our findings suggest that perhaps nature has selected for a single-base addition to control folding to the functional 3D structure. In near-cellular conditions, yeast tRNAPhe and E. coli tRNAAla transcripts fold in a single, cooperative transition only when nearly all of the nucleotides in the cloverleaf are transcribed by indirectly enhancing folding cooperativity. Furthermore, native extensions on the 5' and 3' ends do not interfere with cooperative core folding. This highly controlled cooperative folding has implications for recognition of tRNA by processing and modification enzymes and quality control of tRNA in cells.

Functional Roles of Chelated Magnesium Ions in RNA Folding and Function.

RNA regulates myriad cellular events such as transcription, translation, and splicing. To perform these essential functions, RNA often folds into complex tertiary structures in which its negatively charged ribose-phosphate backbone interacts with metal ions. Magnesium, the most abundant divalent metal ion in cells, neutralizes the backbone, thereby playing essential roles in RNA folding and function. This has been known for more than 50 years, and there are now thousands of in vitro studies, most of which have used ≥10 mM free Mg2+ ions to achieve optimal RNA folding and function. In the cell, however, concentrations of free Mg2+ ions are much lower, with most Mg2+ ions chelated by metabolites. In this Perspective, we curate data from a number of sources to provide extensive summaries of cellular concentrations of metabolites that bind Mg2+ and to estimate cellular concentrations of metabolite-chelated Mg2+ species, in the representative prokaryotic and eukaryotic systems Escherichia coli, Saccharomyces cerevisiae, and iBMK cells. Recent research from our lab and others has uncovered the fact that such weakly chelated Mg2+ ions can enhance RNA function, including its thermodynamic stability, chemical stability, and catalysis. We also discuss how metabolite-chelated Mg2+ complexes may have played roles in the origins of life. It is clear from this analysis that bound Mg2+ should not be simply considered non-RNA-interacting and that future RNA research, as well as protein research, could benefit from considering chelated magnesium.

Long Tracts of Guanines Drive Aggregation of RNA G-Quadruplexes in the Presence of Spermine.

G-Quadruplexes (GQs) are compact, stable structures in DNA and RNA comprised of two or more tiers of quartets whose G-rich motif of tracts of two or more G's occurs commonly within genomes and transcriptomes. While thermodynamically stable in vitro, these structures remain difficult to study in vivo. One approach to understanding GQ in vivo behavior is to test whether conditions and molecules found in cells facilitate their folding. Polyamines are biogenic polycations that interact with RNA. Among common polyamines, spermine contains the highest charge and is found in eukaryotes, making it a good candidate for association with high-charge density nucleic acid structures like GQs. Using a variety of techniques, including ultraviolet-detected thermal denaturation, circular dichroism, size exclusion chromatography, and confocal microscopy, on an array of quadruplex sequence variants, we find that eukaryotic biological concentrations of spermine induce microaggregation of three-tiered G-rich sequences, but not of purely two-tiered structures, although higher spermine concentrations induce aggregation of even these. The formation of microaggregates can also be induced by addition of as little as a single G to a two-tiered structure; moreover, they form at biological temperatures, are sensitive to salt, and can form in the presence of at least some flanking sequence. Notably, GQ aggregation is not observed under prokaryotic-like conditions of no spermine and higher NaCl concentrations. The sequence, polyamine, and salt specificity of microaggregation reported herein have implications for the formation and stability of G-rich nucleic acid aggregates in vivo and for functional roles for understudied GQ sequences with only two quadruplex tiers.

In vivo RNA structural probing of uracil and guanine base-pairing by 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC).

Many biological functions performed by RNAs arise from their in vivo structures. The structure of the same RNA can differ in vitro and in vivo owing in part to the influence of molecules ranging from protons to secondary metabolites to proteins. Chemical reagents that modify the Watson-Crick (WC) face of unprotected RNA bases report on the absence of base-pairing and so are of value to determining structures adopted by RNAs. Reagents have thus been sought that can report on the native RNA structures that prevail in living cells. Dimethyl sulfate (DMS) and glyoxal penetrate cell membranes and inform on RNA secondary structure in vivo through modification of adenine (A), cytosine (C), and guanine (G) bases. Uracil (U) bases, however, have thus far eluded characterization in vivo. Herein, we show that the water-soluble carbodiimide 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC) is capable of modifying the WC face of U and G in vivo, favoring the former nucleobase by a factor of ∼1.5, and doing so in the eukaryote rice, as well as in the Gram-negative bacterium Escherichia coli While both EDC and glyoxal target Gs, EDC reacts with Gs in their typical neutral state, while glyoxal requires Gs to populate the rare anionic state. EDC may thus be more generally useful; however, comparison of the reactivity of EDC and glyoxal may allow the identification of Gs with perturbed pKas in vivo and genome-wide. Overall, use of EDC with DMS allows in vivo probing of the base-pairing status of all four RNA bases.

Design of highly active double-pseudoknotted ribozymes: a combined computational and experimental study.

Design of RNA sequences that adopt functional folds establishes principles of RNA folding and applications in biotechnology. Inverse folding for RNAs, which allows computational design of sequences that adopt specific structures, can be utilized for unveiling RNA functions and developing genetic tools in synthetic biology. Although many algorithms for inverse RNA folding have been developed, the pseudoknot, which plays a key role in folding of ribozymes and riboswitches, is not addressed in most algorithms. For the few algorithms that attempt to predict pseudoknot-containing ribozymes, self-cleavage activity has not been tested. Herein, we design double-pseudoknot HDV ribozymes using an inverse RNA folding algorithm and test their kinetic mechanisms experimentally. More than 90% of the positively designed ribozymes possess self-cleaving activity, whereas more than 70% of negative control ribozymes, which are predicted to fold to the necessary structure but with low fidelity, do not possess it. Kinetic and mutation analyses reveal that these RNAs cleave site-specifically and with the same mechanism as the WT ribozyme. Most ribozymes react just 50- to 80-fold slower than the WT ribozyme, and this rate can be improved to near WT by modification of a junction. Thus, fast-cleaving functional ribozymes with multiple pseudoknots can be designed computationally.