Comparison of crystal structure interactions and thermodynamics for stabilizing mutations in the Tetrahymena ribozyme

Although general mechanisms of RNA folding and catalysis have been elucidated, little is known about how ribozymes achieve structural stability at high temperature. A previous in vitro evolution experiment identified a small number of mutations that significantly increase the thermostability of the tertiary structure of the Tetrahymena ribozyme. Because we also determined the crystal structure of this thermostable ribozyme, we have for the first time the opportunity to compare the structural interactions and thermodynamic contributions of individual nucleotides in a ribozyme. We investigated the contribution of five mutations to thermostability by using temperature gradient gel electrophoresis. Unlike the case with several well-studied proteins, the effects of individual mutations on thermostability of this RNA were highly context dependent. The three most important mutations for thermostability were actually destabilizing in the wild-type background. A269G and A304G contributed to stability only when present as a pair, consistent with their proximity in the ribozyme structure. In an evolutionary context, this work supports and extends the idea that one advantage of protein enzyme systems over an RNA world is the ability of proteins to accumulate stabilizing single-site mutations, whereas RNA may often require much rarer double mutations to improve the stability of both its tertiary and secondary structures.     

Maximizing RNA folding rates: a balancing act

Large ribozymes typically require very long times to refold into their active conformation in vitro, because the RNA is easily trapped in metastable misfolded structures. Theoretical models show that the probability of misfolding is reduced when local and long-range interactions in the RNA are balanced. Using the folding kinetics of the Tetrahymena ribozyme as an example, we propose that folding rates are maximized when the free energies of forming independent domains are similar to each other. A prediction is that the folding pathway of the ribozyme can be reversed by inverting the relative stability of the tertiary domains. This result suggests strategies for optimizing ribozyme sequences for therapeutics and structural studies.     

The tolerance to exchanges of the Watson Crick base pair in the hammerhead ribozyme core is determined by surrounding elements

Tertiary interacting elements are important features of functional RNA molecules, for example, in all small nucleolytic ribozymes. The recent crystal structure of a tertiary stabilized type I hammerhead ribozyme revealed a conventional Watson-Crick base pair in the catalytic core, formed between nucleotides C3 and G8. We show that any Watson-Crick base pair between these positions retains cleavage competence in two type III ribozymes. In the Arabidopsis thaliana sequence, only moderate differences in cleavage rates are observed for the different base pairs, while the peach latent mosaic viroid (PLMVd) ribozyme exhibits a preference for a pyrimidine at position 3 and a purine at position 8. To understand these differences, we created a series of chimeric ribozymes in which we swapped sequence elements that surround the catalytic core. The kinetic characterization of the resulting ribozymes revealed that the tertiary interacting loop sequences of the PLMVd ribozyme are sufficient to induce the preference for Y3-R8 base pairs in the A. thaliana hammerhead ribozyme. In contrast to this, only when the entire stem-loops I and II of the A. thaliana sequences are grafted on the PLMVd ribozyme is any Watson-Crick base pair similarly tolerated. The data provide evidence for a complex interplay of secondary and tertiary structure elements that lead, mediated by long-range effects, to an individual modulation of the local structure in the catalytic core of different hammerhead ribozymes.     

Structure and function of regulatory RNA elements: Ribozymes that regulate gene expression

Since their discovery in the 1980s, it has gradually become apparent that there are several functional classes of naturally occurring ribozymes. These include ribozymes that mediate RNA splicing (the Group I and Group II introns, and possibly the RNA components of the spliceosome), RNA processing ribozymes (RNase P, which cleaves precursor tRNAs and other structural RNA precursors), the peptidyl transferase center of the ribosome, and small, self-cleaving genomic ribozymes (including the hammerhead, hairpin, HDV and VS ribozymes). The most recently discovered functional class of ribozymes include those that are embedded in the untranslated regions of mature mRNAs that regulate the gene's translational expression. These include the prokaryotic glmS ribozyme, a bacterial riboswitch, and a variant of the hammerhead ribozyme, which has been found embedded in mammalian mRNAs. With the discovery of a mammalian riboswitch ribozyme, the question of how an embedded hammerhead ribozyme's switching mechanism works becomes a compelling question. Recent structural results suggest several possibilities.     

The solution structure of the VS ribozyme active site loop reveals a dynamic "hot-spot"

The VS ribozyme is the largest ribozyme in its class and is also the least structurally characterized thus far. The current working model of the VS ribozyme locates the active site in stem-loop VI. The solution structure of this active site loop was determined using high resolution NMR spectroscopy. The structure reveals that the ground-state conformation of the active site differs significantly from that determined previously from chemical structure probing and mutational analysis of the ribozyme in its active conformation, which contains several looped out bases. In contrast, the base-pairing scheme found for the isolated loop contains three mismatched base-pairs: an A+-C, a G-U wobble, and a sheared G-A base-pair and no looped out bases. Dynamics observed within the active site loop provide insight into the mechanism by which the RNA can rearrange its secondary structure into an "activated" conformation prior to cleavage. These findings lend support to the idea that RNA secondary structure is more fluid than once believed and that a better understanding of structure and dynamic features of ribozymes is required to unravel the intricacies of their catalytic abilities.     

Identification of genes that function in the TNF-alpha-mediated apoptotic pathway using randomized hybrid ribozyme libraries

Now that the sequences of many genomes are available, methods are required for the rapid identification of functional genes. We describe here a simple system for the isolation of genes that function in the tumor necrosis factor-alpha (TNF-alpha)-mediated pathway of apoptosis, using RNA helicase-associated ribozyme libraries with randomized substrate-binding arms. Because target-site accessibility considerably limits the effective use of intracellular ribozymes, the effectiveness of a conventional ribozyme library has been low. To overcome this obstacle, we attached to ribozymes an RNA motif (poly(A)-tail) able to interact with endogenous RNA helicase(s) so that the resulting helicase-attached, hybrid ribozymes can more easily attack target sites regardless of their secondary or tertiary structures. When the phenotype of cells changes upon introduction of a ribozyme library, genes responsible for these changes may be identified by sequencing the active ribozyme clones. In the case of TNF-alpha-mediated apoptosis, when a ribozyme library was introduced into MCF-7 cells, surviving clones were completely or partially resistant to TNF-alpha-induced apoptosis. We identified many pro-apoptotic genes and partial sequences of previously uncharacterized genes using this method. Our gene discovery system should be generally applicable to the identification of functional genes in various systems.     

Rearrangement of a stable RNA secondary structure during VS ribozyme catalysis

The Neurospora VS ribozyme recognizes and cleaves a substrate RNA that contains a GC-rich stem loop. In contrast to most RNA secondary structures that are stable during tertiary or quaternary folding, this substrate undergoes extensive ribozyme-induced rearrangement in the presence of magnesium in which the base pairings of at least seven of the ten nucleotides in the stem are changed. This conformational switch is essential for catalytic activity with the wild-type substrate and creates a metal-binding secondary structure motif near the cleavage site. Base pair rearrangement is accompanied by bulging a cytosine from the middle of the stem, indicating that ribozymes may perform base flipping, an activity previously observed only with protein enzymes that modify DNA.     

Tertiary structure stabilization promotes hairpin ribozyme ligation.

The hairpin ribozyme catalyzes a reversible RNA cleavage reaction that participates in processing intermediates of viral satellite RNA replication in plants. A minimal hairpin ribozyme consists of two helix-loop-helix segments. These segments associate noncoaxially in the active folded structure in a way that brings catalytically important loop nucleotides into close proximity. The hairpin ribozyme in the satellite RNA of Tobacco Ringspot Virus assembles in the context of a four-way helical junction. Recent physical characterization of hairpin ribozyme structures using fluorescence resonance energy transfer demonstrated enhanced stability of the folded structure in the context of a four-way helical junction compared to minimal hairpin ribozyme variants. Analysis of the functional consequences of this modification of the helical junction has revealed two changes in the hairpin ribozyme kinetic mechanism. First, ribozymes with a four-way helical junction bind 3' cleavage products with much higher affinity than minimal hairpin ribozymes, evidence that tertiary interactions within the folded structure contribute to product binding energy. Second, the balance between ligation and cleavage shifts in favor of ligation. The enhanced ligation activity of hairpin ribozymes that contain a four-way helical junction supports the notion that tertiary structure stability is a major determinant of the hairpin ribozyme proficiency as a ligase and illustrates the link between RNA structure and biological function.     

Fast formation of the P3-P7 pseudoknot: a strategy for efficient folding of the catalytically active ribozyme

Formation of the P3-P7 pseudoknot structure, the core of group I ribozymes, requires long-range base pairing. Study of the Tetrahymena ribozyme appreciates the hierarchical folding of the large, multidomain RNA, in which the P3-P7 core folds significantly slower than do the other domains. Here we explored the formation of the P3-P7 pseudoknot of the Candida ribozyme that has been reported to concertedly fold to the catalytically active structure with a rate constant of 2 min(-1). We demonstrate that pseudoknot formation occurs during the rapid ribozyme compaction, coincident with formation of many tertiary interactions of the ribozyme. A low physiological concentration of magnesium (1.5 mM) is sufficient to fully support the pseudoknot formation. The presence of nonnative intermediates containing an unfolded P3-P7 region is evident. However, catalysis-based analysis shows these nonnative intermediates are stable and fail to convert to the catalytically active structure, suggesting that rapid pseudoknot formation is essential for folding of the active ribozyme. Interestingly, RNAstructure predicts no stable Alt P3 structure for the Candida ribozyme, but two stable Alt P3s for the Tetrahymena ribozyme, explaining the dramatic difference in folding of the P3-P7 core of these two ribozymes. We propose that rapid formation of the P3-P7 pseudoknot represents a folding strategy ensuring efficient production of the catalytically active structure of group I ribozymes, which sheds new light on the mechanism of effective ribozyme folding in vivo.     

Helix-length compensation studies reveal the adaptability of the VS ribozyme architecture

Compensatory mutations in RNA are generally regarded as those that maintain base pairing, and their identification forms the basis of phylogenetic predictions of RNA secondary structure. However, other types of compensatory mutations can provide higher-order structural and evolutionary information. Here, we present a helix-length compensation study for investigating structure-function relationships in RNA. The approach is demonstrated for stem-loop I and stem-loop V of the Neurospora VS ribozyme, which form a kissing-loop interaction important for substrate recognition. To rapidly characterize the substrate specificity (k(cat)/K(M)) of several substrate/ribozyme pairs, a procedure was established for simultaneous kinetic characterization of multiple substrates. Several active substrate/ribozyme pairs were identified, indicating the presence of limited substrate promiscuity for stem Ib variants and helix-length compensation between stems Ib and V. 3D models of the I/V interaction were generated that are compatible with the kinetic data. These models further illustrate the adaptability of the VS ribozyme architecture for substrate cleavage and provide global structural information on the I/V kissing-loop interaction. By exploring higher-order compensatory mutations in RNA our approach brings a deeper understanding of the adaptability of RNA structure, while opening new avenues for RNA research.     

Ribozymes which cleave arenavirus RNAs: identification of susceptible target sites and inhibition by target site secondary structure.

The development of safe and effective antiviral agents has been a slow process, largely because of the difficulty in distinguishing between virus and host functions; materials toxic to the virus are frequently harmful also to the host in which the agent resides. Recently, techniques which target nucleic acid sequences as a means of reducing gene expression have emerged. This antisense armamentarium includes ribozymes, RNA enzymes which cleave other RNA molecules in a sequence-specific manner. We wish to assess the ability of ribozymes to control animal virus infection. Reasoning that the viruses most vulnerable to ribozyme intervention will be those whose complete life cycle is based on RNA (with no DNA stage), we have begun to develop ribozymes directed toward lymphocytic choriomeningitis virus (LCMV), the prototype of the arenavirus family. Using ribozymes of the hammerhead variety, we have identified several sites on the LCMV genome which can be efficiently cleaved in trans. The efficiency of cleavage is site dependent, and we demonstrate that secondary structure at the target site can abolish ribozyme cleavage. Computer-assisted analysis indicates that much of the LCMV genome may be involved in base pairing, which may render it similarly resistant to ribozyme attack. The few remaining open regions of LCMV lack a GUC target site, on which most studies to date have relied. Here we show that AUC, CUC, and AUU are alternative sites which can be cleaved by trans-acting ribozymes. This finding is important given the aforementioned restriction of available sites, imposed by secondary structure.     

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.     

Additional roles of a peripheral loop-loop interaction in the Neurospora VS ribozyme

Many RNAs contain tertiary interactions that contribute to folding the RNA into its functional 3D structure. In the VS ribozyme, a tertiary loop-loop kissing interaction involving stem-loops I and V is also required to rearrange the secondary structure of stem-loop I such that nucleotides at the base of stem I, which contains the cleavage-ligation site, can adopt the conformation required for activity. In the current work, we have used mutants that constitutively adopt the catalytically permissive conformation to search for additional roles of the kissing interaction in vitro. Using mutations that disrupt or restore the kissing interaction, we find that the kissing interaction contributes ~1000-fold enhancement to the rates of cleavage and ligation. Large Mg(2+)-dependent effects on equilibrium were also observed: in the presence of the kissing interaction cleavage is favored >10-fold at micromolar concentrations of Mg(2+); whereas ligation is favored >10-fold at millimolar concentrations of Mg(2+). In the absence of the kissing interaction cleavage exceeds ligation at all concentrations of Mg(2+). These data provide evidence that the kissing interaction strongly affects the observed cleavage and ligation rate constants and the cleavage-ligation equilibrium of the ribozyme.     

Concordant exploration of the kinetics of RNA folding from global and local perspectives

Time-resolved small-angle X-ray scattering (SAXS) with millisecond time-resolution reveals two discrete phases of global compaction upon Mg2+-mediated folding of the Tetrahymena thermophila ribozyme. Electrostatic relaxation of the RNA occurs rapidly and dominates the first phase of compaction during which the observed radius of gyration (R(g)) decreases from 75 angstroms to 55 angstroms. A further decrease in R(g) to 45 angstroms occurs in a well-defined second phase. An analysis of mutant ribozymes shows that the latter phase depends upon the formation of long-range tertiary contacts within the P4-P6 domain of the ribozyme; disruption of the three remaining long-range contacts linking the peripheral helices has no effect on the 55-45 angstroms compaction transition. A better understanding of the role of specific tertiary contacts in compaction was obtained by concordant time-resolved hydroxyl radical (OH) analyses that report local changes in the solvent accessibility of the RNA backbone. Comparison of the global and local measures of folding shows that formation of a subset of native tertiary contacts (i.e. those defining the ribozyme core) can occur within a highly compact ensemble whose R(g) is close to that of the fully folded ribozyme. Analyses of additional ribozyme mutants and reaction conditions establish the generality of the rapid formation of a partially collapsed state with little to no detectable tertiary structure. These studies directly link global RNA compaction with formation of tertiary structure as the molecule acquires its biologically active structure, and underscore the strong dependence on salt of both local and global measures of folding kinetics.     

Ionization of a critical adenosine residue in the neurospora Varkud Satellite ribozyme active site

The Varkud Satellite (VS) ribozyme catalyzes a site-specific self-cleavage reaction that generates 5'-OH and 2',3'-cyclic phosphate products. Other ribozymes that perform an equivalent reaction appear to employ ionization of an active site residue, either to neutralize the negatively charged transition state or to act as a general acid-base catalyst. To test for important base ionization events in the VS ribozyme ligation reaction, we performed nucleotide analogue interference mapping (NAIM) with a series of ionization-sensitive adenosine and cytidine analogues. A756, a catalytically critical residue located within the VS active site, was the only nucleotide throughout the VS ribozyme that displayed the pH-dependent interference pattern characteristic of functional base ionization. We observed unique rescue of 8-azaadenosine (pK(a) 2.2) and purine riboside (pK(a) 2.1) interference at A756 at reduced reaction pH, suggestive of an ionization-specific effect. These results are consistent with protonation and/or deprotonation of A756 playing a direct role in the VS ribozyme reaction mechanism. In addition, NAIM experiments identified several functional groups within the RNA that play important roles in ribozyme folding and/or catalysis. These include residues in helix II, helix VI (730 loop), the II-III-VI and III-IV-V helix junctions, and loop V.     

Trans-activation of the Tetrahymena group I intron ribozyme via a non-native RNA-RNA interaction

The peripheral P2.1 domain of the Tetrahymena group I intron ribozyme has been shown to be non-essential for splicing. We found, however, that separately prepared P2.1 RNA efficiently accelerates the 3' splice-site-specific hydrolysis reaction of a mutant ribozyme lacking both P2.1 and its upstream region in trans. We report here the unusual properties of this trans-activation. Compensatory mutational analysis revealed that non-native long-range base-pairings between the loop region of P2.1 RNA and L5c region of the mutant ribozyme are needed for the activation in spite of the fact that P2.1 forms base-pairings with P9.1 in the Tetrahymena ribozyme. The trans -activation depends on the non-native RNA-RNA interaction together with the higher order structure of P2.1 RNA. This activation is unique among the known trans-activations that utilize native tertiary interactions or RNA chaperons.     

RNA double cleavage by a hairpin-derived twin ribozyme

The hairpin ribozyme is a small catalytic RNA that catalyses reversible sequence-specific RNA hydrolysis in trans. It consists of two domains, which interact with each other by docking in an antiparallel fashion. There is a region between the two domains acting as a flexible hinge for interdomain interactions to occur. Hairpin ribozymes with reverse-joined domains have been constructed by dissecting the domains at the hinge and rejoining them in reverse order. We have used both the conventional and reverse-joined hairpin ribozymes for the design of a hairpin-derived twin ribozyme. We show that this twin ribozyme cleaves a suitable RNA substrate at two specific sites while maintaining the target specificity of the individual monoribozymes. For characterisation of the studied ribozymes we have evaluated a quantitative assay of sequence-specific ribozyme activity using fluorescently labelled RNA substrates in conjunction with an automated DNA sequencer. This assay was found to be applicable with hairpin and hairpin-derived ribozymes. The results demonstrate the potential of hairpin ribozymes for multi-target strategies of RNA cleavage and suggest the possibility for employing hairpin-derived twin ribozymes as powerful tools for RNA manipulation in vitro and in vivo.