Increased Ribozyme Activity in Crowded Solutions

Noncoding RNAs must function in the crowded environment of the cell. Previous small-angle x-ray scattering experiments showed that molecular crowders stabilize the structure of the Azoarcus group I ribozyme, allowing the ribozyme to fold at low physiological Mg²⁺ concentrations. Here, we used an RNA cleavage assay to show that the PEG and Ficoll crowder molecules increased the biochemical activity of the ribozyme, whereas sucrose did not. Crowding lowered the Mg²⁺ threshold at which activity was detected and increased total RNA cleavage at high Mg²⁺ concentrations sufficient to fold the RNA in crowded or dilute solution. After correcting for solution viscosity, the observed reaction rate was proportional to the fraction of active ribozyme. We conclude that molecular crowders stabilize the native ribozyme and favor the active structure relative to compact inactive folding intermediates.     

Metal Ion Dependence of Cooperative Collapse Transitions in RNA

Positively charged counterions drive RNA molecules into compact configurations that lead to their biologically active structures. To understand how the valence and size of the cations influences the collapse transition in RNA, small-angle X-ray scattering was used to follow the decrease in the radius of gyration (R_g) of the Azoarcus and Tetrahymena ribozymes in different cations. Small, multivalent cations induced the collapse of both ribozymes more efficiently than did monovalent ions. Thus, the cooperativity of the collapse transition depends on the counterion charge density. Singular value decomposition of the scattering curves showed that folding of the smaller and more thermostable Azoarcus ribozyme is well described by two components, whereas collapse of the larger Tetrahymena ribozyme involves at least one intermediate. The ion-dependent persistence length, extracted from the distance distribution of the scattering vectors, shows that the Azoarcus ribozyme is less flexible at the midpoint of transition in low-charge-density ions than in high-charge-density ions. We conclude that the formation of sequence-specific tertiary interactions in the Azoarcus ribozyme overlaps with neutralization of the phosphate charge, while tertiary folding of the Tetrahymena ribozyme requires additional counterions. Thus, the stability of the RNA structure determines its sensitivity to the valence and size of the counterions.     

A kinetic and thermodynamic framework for the Azoarcus group I ribozyme reaction

Determination of quantitative thermodynamic and kinetic frameworks for ribozymes derived from the Azoarcus group I intron and comparisons to their well-studied analogs from the Tetrahymena group I intron reveal similarities and differences between these RNAs. The guanosine (G) substrate binds to the Azoarcus and Tetrahymena ribozymes with similar equilibrium binding constants and similar very slow association rate constants. These and additional literature observations support a model in which the free ribozyme is not conformationally competent to bind G and in which the probability of assuming the binding-competent state is determined by tertiary interactions of peripheral elements. As proposed previously, the slow binding of guanosine may play a role in the specificity of group I intron self-splicing, and slow binding may be used analogously in other biological processes. The internal equilibrium between ribozyme-bound substrates and products is similar for these ribozymes, but the Azoarcus ribozyme does not display the coupling in the binding of substrates that is observed with the Tetrahymena ribozyme, suggesting that local preorganization of the active site and rearrangements within the active site upon substrate binding are different for these ribozymes. Our results also confirm the much greater tertiary binding energy of the 5′-splice site analog with the Azoarcus ribozyme, binding energy that presumably compensates for the fewer base-pairing interactions to allow the 5′-exon intermediate in self splicing to remain bound subsequent to 5′-exon cleavage and prior to exon ligation. Most generally, these frameworks provide a foundation for design and interpretation of experiments investigating fundamental properties of these and other structured RNAs.     

Ion Condensation onto Ribozyme Is Site Specific and Fold Dependent

The highly charged RNA molecules, with each phosphate carrying a single negative charge, cannot fold into well-defined architectures with tertiary interactions in the absence of ions. For ribozymes, divalent cations are known to be more efficient than monovalent ions in driving them to a compact state, although Mg²⁺ ions are needed for catalytic activities. Therefore, how ions interact with RNA is relevant in understanding RNA folding. It is often thought that most of the ions are territorially and nonspecifically bound to the RNA, as predicted by the counterion condensation theory. Here, we show using simulations of Azoarcus ribozyme, based on an accurate coarse-grained three-site interaction model with explicit divalent and monovalent cations, that ion condensation is highly specific and depends on the nucleotide position. The regions with high coordination between the phosphate groups and the divalent cations are discernible even at very low Mg²⁺ concentrations when the ribozyme does not form tertiary interactions. Surprisingly, these regions also contain the secondary structural elements that nucleate subsequently in the self-assembly of RNA, implying that ion condensation is determined by the architecture of the folded state. These results are in sharp contrast to interactions of ions (monovalent and divalent) with rigid charged rods, in which ion condensation is uniform and position independent. The differences are explained in terms of the dramatic nonmonotonic shape fluctuations in the ribozyme as it folds with increasing Mg²⁺ or Ca²⁺ concentration.     

Structural Rearrangements Linked to Global Folding Pathways of the Azoarcus Group I Ribozyme

Stable RNAs must fold into specific three-dimensional structures to be biologically active, yet many RNAs form metastable structures that compete with the native state. Our previous time-resolved footprinting experiments showed that Azoarcus group I ribozyme forms its tertiary structure rapidly (τ < 30 ms) without becoming significantly trapped in kinetic intermediates. Here, we use stopped-flow fluorescence spectroscopy to probe the global folding kinetics of a ribozyme containing 2-aminopurine in the loop of P9. The modified ribozyme was catalytically active and exhibited two equilibrium folding transitions centered at 0.3 and 1.6 mM Mg²⁺, consistent with previous results. Stopped-flow fluorescence revealed four kinetic folding transitions with observed rate constants of 100, 34, 1, and 0.1 s^− 1 at 37 °C. From comparison with time-resolved Fe(II)-ethylenediaminetetraacetic acid footprinting of the modified ribozyme under the same conditions, these folding transitions were assigned to formation of the I_C intermediate, tertiary folding and docking of the nicked P9 tetraloop, reorganization of the P3 pseudoknot, and refolding of nonnative conformers, respectively. The footprinting results show that 50-60% of the modified ribozyme folds in less than 30 ms, while the rest of the RNA population undergoes slow structural rearrangements that control the global folding rate. The results show how small perturbations to the structure of the RNA, such as a nick in P9, populate kinetic folding intermediates that are not observed in the natural ribozyme.     

Mg2+ Impacts the Twister Ribozyme through Push-Pull Stabilization of Nonsequential Phosphate Pairs

RNA molecules perform a variety of biological functions for which the correct three-dimensional structure is essential, including as ribozymes where they catalyze chemical reactions. Metal ions, especially Mg²⁺, neutralize these negatively charged nucleic acids and specifically stabilize RNA tertiary structures as well as impact the folding landscape of RNAs as they assume their tertiary structures. Specific binding sites of Mg²⁺ in folded conformations of RNA have been studied extensively; however, the full range of interactions of the ion with compact intermediates and unfolded states of RNA is challenging to investigate, and the atomic details of the mechanism by which the ion facilitates tertiary structure formation is not fully known. Here, umbrella sampling combined with oscillating chemical potential Grand Canonical Monte Carlo/molecular dynamics simulations are used to capture the energetics and atomic-level details of Mg²⁺-RNA interactions that occur along an unfolding pathway of the Twister ribozyme. The free energy profiles reveal stabilization of partially unfolded states by Mg²⁺, as observed in unfolding experiments, with this stabilization being due to increased sampling of simultaneous interactions of Mg²⁺ with two or more nonsequential phosphate groups. Notably, these results indicate a push-pull mechanism in which the Mg²⁺-RNA interactions actually lead to destabilization of specific nonsequential phosphate-phosphate interactions (i.e., pushed apart), whereas other interactions are stabilized (i.e., pulled together), a balance that stabilizes unfolded states and facilitates the folding of Twister, including the formation of hydrogen bonds associated with the tertiary structure. This study establishes a better understanding of how Mg²⁺-ion interactions contribute to RNA structural properties and stability.     

Mechanisms of covalent self-assembly of the Azoarcus ribozyme from four fragment oligonucleotides

RNA oligomers of length 40–60 nt can self-assemble into covalent versions of the Azoarcus group I intron ribozyme. This process requires a series of recombination reactions in which the internal guide sequence of a nascent catalytic complex makes specific interactions with a complement triplet, CAU, in the oligomers. However, if the CAU were mutated, promiscuous self-assembly may be possible, lessening the dependence on a particular set of oligomer sequences. Here, we assayed whether oligomers containing mutations in the CAU triplet could still self-construct Azoarcus ribozymes. The mutations CAC, CAG, CUU and GAU all inhibited self-assembly to some degree, but did not block it completely in 100 mM MgCl₂. Oligomers containing the CAC mutation retained the most self-assembly activity, while those containing GAU retained the least, indicating that mutations more 5′ in this triplet are the most deleterious. Self-assembly systems containing additional mutant locations were progressively less functional. Analyses of properly self-assembled ribozymes revealed that, of two recombination mechanisms possible for self-assembly, termed ‘tF2’ and ‘R2F2’, the simpler one-step ‘tF2’ mechanism is utilized when mutations exist. These data suggest that self-assembling systems are more facile than previously believed, and have relevance to the origin of complex ribozymes during the RNA World.     

Molecular beacons as probes of RNA unfolding under native conditions

Hybridization of fluorescent molecular beacons provides real-time detection of RNA secondary structure with high specificity. We used molecular beacons to measure folding and unfolding rates of the Tetrahymena group I ribozyme under native conditions. A molecular beacon targeted against 15 nt in the 5′ strand of the P3 helix specifically hybridized with misfolded forms of the ribozyme, without invading the native tertiary structure. The beacon associated with the misfolded ribozyme 300 times more slowly than with an unstructured oligonucleotide containing the same target sequence, suggesting that the misfolded ribozyme core remains structured in the absence of Mg²⁺. The rate of beacon hybridization under native conditions revealed a linear relationship between the free energy of unfolding and Mg²⁺ concentration. A small fraction of the RNA population unfolded very rapidly, suggesting parallel unfolding in one step or through misfolded intermediates.     

A class I ligase ribozyme with reduced Mg2+ dependence: Selection,sequence analysis,and identification of functional tertiary interactions

The class I ligase was among the first ribozymes to have been isolated from random sequences and represents the catalytic core of several RNA-directed RNA polymerase ribozymes. The ligase is also notable for its catalytic efficiency and structural complexity. Here, we report an improved version of this ribozyme, arising from selection that targeted the kinetics of the chemical step. Compared with the parent ribozyme, the improved ligase achieves a modest increase in rate enhancement under the selective conditions and shows a sharp reduction in [Mg²⁺] dependence. Analysis of the sequences and kinetics of successful clones suggests which mutations play the greatest part in these improvements. Moreover, backbone and nucleobase interference maps of the parent and improved ligase ribozymes complement the newly solved crystal structure of the improved ligase to identify the functionally significant interactions underlying the catalytic ability and structural complexity of the ligase ribozyme.     

Conformational heterogeneity at position U37 of an all-RNA hairpin ribozyme with implications for metal binding and the catalytic structure of the S-turn

The hairpin ribozyme is an RNA enzyme that performs site-specific phosphodiester bond cleavage between nucleotides A-1 and G+1 within its cognate substrate. Previous functional studies revealed that the minimal hairpin ribozyme exhibited "gain-of-function" cleavage properties resulting from U39C or U39 to propyl linker (C3) modifications. Furthermore, each "mutant" displayed different magnesium-dependence in its activity. To investigate the molecular basis for these gain-of-function variants, crystal structures of minimal, junctionless hairpin ribozymes were solved in native (U39), and mutant U39C and U39(C3) forms. The results revealed an overall molecular architecture comprising two docked internal loop domains folded into a wishbone shape, whose tertiary interface forms a sequestered active site. All three minimal hairpin ribozymes bound Co(NH(3))(6)(3+) at G21/A40, the E-loop/S-turn boundary. The native structure also showed that U37 of the S-turn adopts both sequestered and exposed conformations that differ by a maximum displacement of 13 A. In the sequestered form, the U37 base packs against G36, and its 2'-hydroxyl group forms a water mediated hydrogen bond to O4' of G+1. These interactions were not observed in previous four-way-junction hairpin ribozyme structures due to crystal contacts with the U1A splicing protein. Interestingly, the U39C and U39(C3) mutations shifted the equilibrium conformation of U37 into the sequestered form through formation of new hydrogen bonds in the S-turn, proximal to the essential nucleotide A38. A comparison of all three new structures has implications for the catalytically relevant conformation of the S-turn and suggests a rationale for the distinctive metal dependence of each mutant.     

Requirement for canonical base pairing in the short pseudoknot structure of genomic hepatitis delta virus ribozyme

The tertiary structure of the 3′-cleaved product of the genomic hepatitis delta virus (HDV) ribozyme was solved by X-ray crystallographic analysis. In this structure, three single-stranded regions (SSrA, -B and -C) interact intricately with one another via hydrogen bonds between nucleotide bases, phosphate oxygens and 2′-OHs to form a nested double pseudoknot structure. Among these interactions, two Watson–Crick (W–C) base pairs, 726G–710C and 727G–709C, that form between SSrA and SSrC (P1.1) seem to be especially important for compact folding. To characterize the importance of these base pairs, ribozymes were subjected to in vitro selection from a pool of RNA molecules randomly substituted at positions 709, 710, 726 and 727. The results establish the importance of the two W–C base pairs for activity, although some mutants are active with one G–C base pair. In addition, the kinetic parameters were analyzed in all 16 combinations with two canonical base pairs. Comparison of variant ribozymes with the wild-type ribozyme reveals that the difference in reaction rates for these variants (ΔΔG^‡) is not simply accounted for by the differences in the stability of P1.1 (ΔΔG⁰₃₇). The role played by Mg²⁺ ions in formation of the P1.1 structure is also discussed.     

Concerted folding of a Candida ribozyme into the catalytically active structure posterior to a rapid RNA compaction

Folding of the major population of Tetrahymena intron RNA into the catalytically active structure is trapped in a slow pathway. In this report, folding of Candida albicans intron was investigated using the trans-acting Ca.L-11 ribozyme as a model. We demonstrated that both the catalytic activity (k_obs) and compact folding equilibrium of Ca.L-11 are strongly dependent on Mg²⁺ at physiological concentrations, with both showing an Mg²⁺ Hill coefficient of 3. Formation of the compact structure of Ca.L-11 is shown to occur very rapidly, on a subsecond time scale similar to that of RNase T1 cleavage. Most of the ribozyme RNA population folds into the catalytically active structure with a rate constant of 2 min^–1 at 10 mM Mg²⁺; neither slower kinetics nor obvious Mg²⁺ inhibition is observed. These results suggest that folding of the Ca.L-11 ribozyme is initiated by a rapid magnesium-dependent RNA compaction, which is followed by a slower searching for the native contacts to form the catalytically active structure without interference from the long-lived trapped states. This model thus provides an ideal system to address a range of interesting aspects of RNA folding, such as conformational searching, ion binding and the role of productive intermediates.     

Catalytic cleavage of cis- and trans-acting antigenomic delta ribozymes in the presence of various divalent metal ions

Catalytic activity of four structural variants of the antigenomic delta ribozyme, two cis- and two trans-acting, has been compared in the presence of selected divalent metal ions that effectively support catalysis. The ribozymes differ in regions that are not directly involved in formation of the ribozyme active site: the region immediately preceding the catalytic cleavage site, the P4 stem and a stretch of the viral RNA sequence extending the minimal ribozyme sequence at its 3′-terminus. The variants show high cleavage activity in the presence of Mg²⁺, Ca²⁺ and Mn²⁺, lower with Co²⁺ and Sr²⁺ and some variants are also active with Cd²⁺ and Zn²⁺ ions. In the presence of a particular metal ion the ribozymes cleave, however with different initial rates, according to pseudo-first or higher order kinetics and to different final cleavage extents. On the other hand, relatively small differences are observed in the reactions induced by various metal ions. The cleavage of trans-acting ribozymes induced by Mg²⁺ is partially inhibited in the presence of Na⁺, spermidine and some other divalent metal ions. The inert Co(NH₃)₆³⁺ complex is unable to support catalysis, as reported earlier for the genomic ribozyme. The results are discussed in terms of the influence of structural elements peripheral to the ribozyme active site on its cleavage rate and efficiency as well as the role of metal ions in the cleavage mechanism. Some implications concerning further studies and possible applications of delta ribozymes are also considered.     

Structure-function relationships of two closely related group IC3 intron ribozymes from Azoarcus and Synechococcus pre-tRNA

The two group IC3 pre-tRNA introns from Azoarcus and Synechococcus share very analogous secondary structures. They are small group I ribozymes that possess only two peripheral domains, P2 and P9. However, the 3′-splice site hydrolysis activity of the Synechococcus ribozyme critically depends on P2 whereas that of Azoarcus does not, indicating that the structure–function relationships of the two ribozymes are strikingly different despite their structural resemblance. To identify the element(s) that determines the catalytic properties of these ribozymes, we undertook analyses of chimeric ribozymes prepared by swapping their structural elements. We found that the difference can be attributed to a small number of nucleotides within the conserved core region. Further analysis by employing in vitro selection revealed that a base triple interaction (P4bp3 × J6/7-2) is a critical element for determining activity and suggests the existence of a novel base quintuple involving the base triple P4bp5 × J8/7-5.     

Fast folding of a ribozyme by stabilizing core interactions: evidence for multiple folding pathways in RNA

Folding of the Tetrahymena ribozyme under physiological conditions in vitro is limited by slow conversion of long-lived intermediates to the active structure. These intermediates arise because the most stable domain of the ribozyme folds 10-50 times more rapidly than the core region containing helix P3. Native gel electrophoresis and time-resolved X-ray-dependent hydroxyl radical cleavage revealed that mutations that weaken peripheral interactions between domains accelerated folding fivefold, while a point mutation that stabilizes P3 enabled 80 % of the mutant RNA to reach the native conformation within 30 seconds at 22 degrees C. The P3 mutation increased the folding rate of the catalytic core as much as 50-fold, so that both domains of the ribozyme were formed at approximately the same rate. The results show that the ribozyme folds rapidly without significantly populating metastable intermediates when native interactions in the ribozyme core are stabilized relative to peripheral structural elements.     

Understanding the Role of Three-Dimensional Topology in Determining the?Folding Intermediates of Group I Introns

Many RNA molecules exert their biological function only after folding to unique three-dimensional structures. For long, noncoding RNA molecules, the complexity of finding the native topology can be a major impediment to correct folding to the biologically active structure. An RNA molecule may fold to a near-native structure but not be able to continue to the correct structure due to a topological barrier such as crossed strands or incorrectly stacked helices. Achieving the native conformation thus requires unfolding and refolding, resulting in a long-lived intermediate. We investigate the role of topology in the folding of two phylogenetically related catalytic group I introns, the Twort and Azoarcus group I ribozymes. The kinetic models describing the Mg²⁺-mediated folding of these ribozymes were previously determined by time-resolved hydroxyl (⋅OH) radical footprinting. Two intermediates formed by parallel intermediates were resolved for each RNA. These data and analytical ultracentrifugation compaction analyses are used herein to constrain coarse-grained models of these folding intermediates as we investigate the role of nonnative topology in dictating the lifetime of the intermediates. Starting from an ensemble of unfolded conformations, we folded the RNA molecules by progressively adding native constraints to subdomains of the RNA defined by the ⋅OH time-progress curves to simulate folding through the different kinetic pathways. We find that nonnative topologies (arrangement of helices) occur frequently in the folding simulations despite using only native constraints to drive the reaction, and that the initial conformation, rather than the folding pathway, is the major determinant of whether the RNA adopts nonnative topology during folding. From these analyses we conclude that biases in the initial conformation likely determine the relative flux through parallel RNA folding pathways.     

Biochemical analysis of pistol self-cleaving ribozymes

Pistol RNAs are members of a distinct class of self-cleaving ribozymes that was recently discovered by using a bioinformatics search strategy. Several hundred pistol ribozymes share a consensus sequence including 10 highly conserved nucleotides and many other modestly conserved nucleotides associated with specific secondary structure features, including three base-paired stems and a pseudoknot. A representative pistol ribozyme from the bacterium Lysinibacillus sphaericus was found to promote RNA strand scission with a rate constant of ∼10 min⁻¹ under physiological Mg²⁺ and pH conditions. The reaction proceeds via the nucleophilic attack of a 2′-oxygen atom on the adjacent phosphorus center, and thus adheres to the same general catalytic mechanism of internal phosphoester transfer as found with all other classes of natural self-cleaving ribozymes discovered to date. Analyses of the kinetic characteristics and the metal ion requirements of the cleavage reaction reveal that members of this ribozyme class likely use several catalytic strategies to promote the rapid cleavage of RNA.     

Dissecting the influence of Mg2+ on 3D architecture and ligand-binding of the guanine-sensing riboswitch aptamer domain

Long-range tertiary interactions determine the three-dimensional structure of a number of metabolite-binding riboswitch RNA elements and were found to be important for their regulatory function. For the guanine-sensing riboswitch of the Bacillus subtilis xpt-pbuX operon, our previous NMR-spectroscopic studies indicated pre-formation of long-range tertiary contacts in the ligand-free state of its aptamer domain. Loss of the structural pre-organization in a mutant of this RNA (G37A/C61U) resulted in the requirement of Mg²⁺ for ligand binding. Here, we investigate structural and stability aspects of the wild-type aptamer domain (Gsw) and the G37A/C61U-mutant (Gsw^loop) of the guanine-sensing riboswitch and their Mg²⁺-induced folding characteristics to dissect the role of long-range tertiary interactions, the link between pre-formation of structural elements and ligand-binding properties and the functional stability. Destabilization of the long-range interactions as a result of the introduced mutations for Gsw^loop or the increase in temperature for both Gsw and Gsw^loop involves pronounced alterations of the conformational ensemble characteristics of the ligand-free state of the riboswitch. The increased flexibility of the conformational ensemble can, however, be compensated by Mg²⁺. We propose that reduction of conformational dynamics in remote regions of the riboswitch aptamer domain is the minimal pre-requisite to pre-organize the core region for specific ligand binding.     

Interaction between RNA molecules of a two-component trans analog of antigenomic HDV ribozyme

A study was made of the association of the RNA components forming a B:LS two-component rans analog of the antigenomic HDV ribozyme. The B:LS ribozyme differed from known trans ribozymes in the sizes and nucleotide sequences of its components (33 and 34 nt, respectively), the topology of its functional parts, and the lack of a very short cleavage product. Compared to the cis ribozyme, B:LS showed similar dependences on the reaction conditions (Mg²⁺ concentration, pH, temperature) and a similar biphasic kinetic curve of self-cleavage. The kinetic model of B:LS self-cleavage (available at www.cardio.ru/labgen/RZ_e.html) describes a possible cause of the biphasic kinetic curve as a change in the rate-limiting step of consecutive conformational transitions accompanying self-cleavage. Another possible cause is an interaction between the molecules involved in cleavage, i.e., multimerization of whole ribozyme molecules with their components or the reaction products. B:LS provides a convenient model for studying such interactions, since the mode of component binding allows generation of 1B:2LS and 2B:1LS complexes as well as complexes with the cleavage products. Nondenaturing PAGE was used to study the factors affecting association and dissociation of the ribozyme components. The possibility of interactions between the RNA components of the cis and trans ribozymes was demonstrated experimentally. It was shown that the ribozyme is capable of multimerization when LS is in excess over B and that the cleavage products are not significantly involved in this process. The results suggest intermolecular interactions for the cleavage of the natural cis ribozyme.