An obligate intermediate along the slow folding pathway of a group II intron ribozyme

Most RNA molecules collapse rapidly and reach the native state through a pathway that contains numerous traps and unproductive intermediates. The D135 group II intron ribozyme is unusual in that it can fold slowly and directly to the native state, despite its large size and structural complexity. Here we use hydroxyl radical footprinting and native gel analysis to monitor the timescale of tertiary structure collapse and to detect the presence of obligate intermediates along the folding pathway of D135. We find that structural collapse and native folding of Domain 1 precede assembly of the entire ribozyme, indicating that D1 contains an on-pathway intermediate to folding of the D135 ribozyme. Subsequent docking of Domains 3 and 5, for which D1 provides a preorganized scaffold, appears to be very fast and independent of one another. In contrast to other RNAs, the D135 ribozyme undergoes slow tertiary collapse to a compacted state, with a rate constant that is also limited by the formation D1. These findings provide a new paradigm for RNA folding and they underscore the diversity of RNA biophysical behaviors.     

Secondary structure content of the HDV ribozyme in 95% formamide.

The Hepatitis Delta Virus (HDV) ribozyme self-cleaving activity in 20 M formamide solutions is unique. Does this catalytic activity result from the conservation of its tertiary structure in 20 M formamide? We followed the ribozyme structure in formamide solutions by monitoring the amount of bound Ethidium Bromide (EB). We were able to measure the quantity of dye bound using time-resolved fluorescence spectroscopy, as an estimate of the ribozyme double helical content. This method, calibrated by using oligonucleotides with defined tertiary structure and denaturing solvents, parallels NMR and UV measurements as a function of temperature. Measurements with the HDV ribozyme lead to three conclusions: (a) both the precursor and product RNAs are structured to 24 M (95% w/w) formamide or 4 M H2O solutions which is equivalent to 4 M H2O; (b) the HDV ribozyme is the only RNA sequence investigated in this study that retains so much structure in formamide; and (c) DNA analogs of formamide resistant HDV ribozyme sequences lose their structure at less than 15 M formamide. Thus, the structural integrity of the HDV ribozyme is an intrinsic property of the RNA molecule and its sequence.     

An alternative route for the folding of large RNAs: apparent two-state folding by a group II intron ribozyme

Despite a growing literature on the folding of RNA, our understanding of tertiary folding in large RNAs derives from studies on a small set of molecular examples, with primary focus on group I introns and RNase P RNA. To broaden the scope of RNA folding models and to better understand group II intron function, we have examined the tertiary folding of a ribozyme (D135) that is derived from the self-splicing ai5gamma intron from yeast mitochondria. The D135 ribozyme folds homogeneously and cooperatively into a compact, well-defined tertiary structure that includes all regions critical for active-site organization and substrate recognition. When D135 was treated with increasing concentrations of Mg(2+) and then subjected to hydroxyl radical footprinting, similar Mg(2+) dependencies were seen for internalization of all regions of the molecule, suggesting a highly cooperative folding behavior. In this work, we show that global folding and compaction of the molecule have the same magnesium dependence as the local folding previously observed. Furthermore, urea denaturation studies indicate highly cooperative unfolding of the ribozyme that is governed by thermodynamic parameters similar to those for forward folding. In fact, D135 folds homogeneously and cooperatively from the unfolded state to its native, active structure, thereby demonstrating functional reversibility in RNA folding. Taken together, the data are consistent with two-state folding of the D135 ribozyme, which is surprising given the size and multi-domain structure of the RNA. The findings establish that the accumulation of stable intermediates prior to formation of the native state is not a universal feature of RNA folding and that there is an alternative paradigm in which the folding landscape is relatively smooth, lacking rugged features that obstruct folding to the native state.     

Disparate HDV ribozyme crystal structures represent intermediates on a rugged free-energy landscape

The hepatitis delta virus (HDV) ribozyme is a member of the class of small, self-cleaving catalytic RNAs found in a wide range of genomes from HDV to human. Both pre- and post-catalysis (precursor and product) crystal structures of the cis-acting genomic HDV ribozyme have been determined. These structures, together with extensive solution probing, have suggested that a significant conformational change accompanies catalysis. A recent crystal structure of a trans-acting precursor, obtained at low pH and by molecular replacement from the previous product conformation, conforms to the product, raising the possibility that it represents an activated conformer past the conformational change. Here, using fluorescence resonance energy transfer (FRET), we discovered that cleavage of this ribozyme at physiological pH is accompanied by a structural lengthening in magnitude comparable to previous trans-acting HDV ribozymes. Conformational heterogeneity observed by FRET in solution appears to have been removed upon crystallization. Analysis of a total of 1.8 µsec of molecular dynamics (MD) simulations showed that the crystallographically unresolved cleavage site conformation is likely correctly modeled after the hammerhead ribozyme, but that crystal contacts and the removal of several 2′-oxygens near the scissile phosphate compromise catalytic in-line fitness. A cis-acting version of the ribozyme exhibits a more dynamic active site, while a G-1 residue upstream of the scissile phosphate favors poor fitness, allowing us to rationalize corresponding changes in catalytic activity. Based on these data, we propose that the available crystal structures of the HDV ribozyme represent intermediates on an overall rugged RNA folding free-energy landscape.     

Thermodynamics and kinetics of the hairpin ribozyme from atomistic folding/unfolding simulations

We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active “docked” conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.     

Terbium-mediated footprinting probes a catalytic conformational switch in the antigenomic hepatitis delta virus ribozyme

The two forms of the hepatitis delta virus ribozyme are derived from the genomic and antigenomic RNA strands of the human hepatitis delta virus (HDV), where they serve a crucial role in pathogen replication by catalyzing site-specific self-cleavage reactions. The HDV ribozyme requires divalent metal ions for formation of its tertiary structure, consisting of a tight double-nested pseudoknot, and for efficient self- (or cis-) cleavage. Comparison of recently solved crystal structures of the cleavage precursor and 3' product indicates that a significant conformational switch is required for catalysis by the genomic HDV ribozyme. Here, we have used the lanthanide metal ion terbium(III) to footprint the precursor and product solution structures of the cis-acting antigenomic HDV ribozyme. Inhibitory Tb(3+) binds with high affinity to similar sites on RNA as Mg(2+) and subsequently promotes slow backbone scission. We find subtle, yet significant differences in the terbium(III) footprinting pattern between the precursor and product forms of the antigenomic HDV ribozyme, consistent with differences in conformation as observed in the crystal structures of the genomic ribozyme. In addition, UV melting profiles provide evidence for a less tight tertiary structure in the precursor. In both the precursor and product we observe high-affinity terbium(III) binding sites in joining sequence J4/2 (Tb(1/2) approximately 4 microM) and loop L3, which are key structural components forming the catalytic core of the HDV ribozyme, as well as in several single-stranded regions such as J1/2 and the L4 tetraloop (Tb(1/2) approximately 50 microM). Sensitized luminescence spectroscopy confirms that there are at least two affinity classes of Tb(3+) binding sites. Our results thus demonstrate that a significant conformational change accompanies catalysis in the antigenomic HDV ribozyme in solution, similar to the catalytic conformational switch observed in crystals of the genomic form, and that structural and perhaps catalytic metal ions bind close to the catalytic core.     

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.     

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.     

Two major tertiary folding transitions of the Tetrahymena catalytic RNA.

The L-21 Tetrahymena ribozyme, an RNA molecule with sequence-specific endoribonuclease activity derived from a self-splicing group I intron, provides a model system for studying the RNA folding problem. A 160 nucleotide, independently folding domain of tertiary structure (the P4-P6 domain) comprises about half of the ribozyme. We now apply Fe(II)-EDTA cleavage to mutants of the ribozyme to explore the role of individual structural elements in tertiary folding of the RNA at equilibrium. Deletion of peripheral elements near the 3' end of the ribozyme destabilizes a region of the catalytic core (P3-P7) without altering the folding of the P4-P6 domain. Three different mutations within the P4-P6 domain that destabilize its folding also shift the folding of the P3-P7 region of the catalytic core to higher MgCl2 concentrations. We conclude that the role of the extended P4-P6 domain and of the 3'-terminal peripheral elements is at least in part to stabilize the catalytic core. The organization of RNA into independently folding domains of tertiary structure may be common in large RNAs, including ribosomal RNAs. Furthermore, the observation of domain-domain interactions in a catalytic RNA supports the feasibility of a primitive spliceosome without any proteins.     

The poly(A) site sequence in HDV RNA alters both extent and rate of self-cleavage of the antigenomic ribozyme

The ribozyme self-cleavage site in the antigenomic sequence of hepatitis delta virus (HDV) RNA is 33-nt downstream of the poly(A) site for the delta antigen mRNA. An HDV antigenomic ribozyme precursor RNA that included the upstream poly(A) processing site was used to test the hypothesis that nonribozyme sequence near the poly(A) site could affect ribozyme activity. Relative to ribozyme precursor without the extra upstream sequences, the kinetic profile for self-cleavage of the longer precursor was altered in two ways. First, only half of the precursor RNA self-cleaved. The cleaved fraction could be increased or decreased with mutations in the upstream sequence. These mutations, which were predicted to alter the relative stability of competing secondary structures within the precursor, changed the distribution of alternative RNA structures that are resolved in native-gel electrophoresis. Second, the active fraction cleaved with an observed rate constant that was higher than that of the ribozyme without the upstream sequences. Moreover, the higher rate constants occurred at lower, near-physiological, divalent metal ion concentrations (1–2 mM). Modulation of ribozyme activity, through competing alternative structures, could be part of a mechanism that allows a biologically significant choice between maturation of the mRNA and processing of replication intermediates.     

Design of a highly reactive HDV ribozyme sequence uncovers facilitation of RNA folding by alternative pairings and physiological ionic strength

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA that resides in the HDV genome and regulates its replication. The native fold of the ribozyme is complex, having two pseudoknots. Earlier work implicated four non-native pairings in slowing pseudoknot formation: Alt 1, Alt 2, Alt 3, and Alt P1. The goal of the present work was design of a kinetically simplified and maximally reactive construct for in vitro mechanistic and structural studies. The initial approach chosen was site-directed mutagenesis in which known alternative pairings were destabilized while leaving the catalytic core intact. Based on prior studies, the G11C/U27Delta double mutant was prepared. However, biphasic kinetics and antisense oligonucleotide response trends opposite those of the well-studied G11C mutant were observed suggesting that new alternative pairings with multiple registers, termed Alt X and Alt Y, had been created. Enzymatic structure mapping of oligonucleotide models supported this notion. This led to a model wherein Alt 2 and the phylogenetically conserved Alt 3 act as "folding guides", facilitating folding of the major population of the RNA molecules by hindering formation of the Alt X and Alt Y registers. Attempts to eliminate the strongest of the Alt X pairings by rational design of a quadruple mutant only resulted in more complex kinetic behavior. In an effort to simultaneously destabilize multiple alternative pairings, studies were carried out on G11C/U27Delta in the presence of urea or increased monovalent ion concentration. Inclusion of physiological ionic strength allowed the goal of monophasic, fast-folding (kobs approximately 60 min(-1)) kinetics to be realized. To account for this, a model is developed wherein Na+, which destabilizes secondary and tertiary structures in the presence of Mg2+, facilitates native folding by destabilizing the multiple alternative secondary structures with a higher-order dependence.     

Folding path of P5abc RNA involves direct coupling of secondary and tertiary structures

Folding mechanisms in which secondary structures are stabilized through the formation of tertiary interactions are well documented in protein folding but challenge the folding hierarchy normally assumed for RNA. However, it is increasingly clear that RNA could fold by a similar mechanism. P5abc, a small independently folding tertiary domain of the Tetrahymena thermophila group I ribozyme, is known to fold by a secondary structure rearrangement involving helix P5c. However, the extent of this rearrangement and the precise stage of folding that triggers it are unknown. We use experiments and simulations to show that the P5c helix switches to the native secondary structure late in the folding pathway and is directly coupled to the formation of tertiary interactions in the A-rich bulge. P5c mutations show that the switch in P5c is not rate-determining and suggest that non-native interactions in P5c aid folding rather than impede it. Our study illustrates that despite significant differences in the building blocks of proteins and RNA, there may be common ways in which they self-assemble.     

Monitoring of an RNA Multistep Folding Pathway by Isothermal Titration Calorimetry

Isothermal titration calorimetry was used to monitor the energetic landscape of a catalytic RNA, specifically that of the hepatitis delta virus ribozyme. Using mutants that isolated various tertiary interactions, the thermodynamic parameters of several ribozyme-substrate intermediates were determined. The results shed light on the impact of several tertiary interactions on the global structure of the ribozyme. In addition, the data indicate that the formation of the P1.1 pseudoknot is the limiting step of the molecular mechanism. Last, as illustrated here, isothermal titration calorimetry appears to be a method of choice for the elucidation of an RNA's folding pathway.     

Magnesium dependence of the amplified conformational switch in the trans-acting hepatitis delta virus ribozyme

The ability of divalent metal ions to participate in both structure formation and catalytic chemistry of RNA enzymes (ribozymes) has made it difficult to separate their cause and effect in ribozyme function. For example, the recently solved crystal structures of precursor and product forms of the cis-cleaving genomic hepatitis delta virus (HDV) ribozyme show a divalent metal ion bound in the active site that is released upon catalysis due to an RNA conformational change. This conformational switch is associated with a repositioning of the catalytically involved base C75 in the active-site cleft, thus controlling catalysis. These findings confirm previous data from fluorescence resonance energy transfer (FRET) on a trans-acting form of the HDV ribozyme that found a global conformational change to accompany catalysis. Here, we further test the conformational switch model by measuring the Mg(2+) dependence of the global conformational change of the trans-acting HDV ribozyme, using circular dichroism and time-resolved FRET as complementary probes of secondary and tertiary structure formation, respectively. We observe significant differences in both structure and Mg(2+) affinity of the precursor and product forms, in the presence and absence of 300 mM Na(+) background. The precursor shortens while the product extends with increasing Mg(2+) concentration, essentially amplifying the structural differences observed in the crystal structures. In addition, the precursor has an approximately 2-fold and approximately 13-fold lower Mg(2+) affinity than the product in secondary and tertiary structure formation, respectively. We also have compared the C75 wild-type with the catalytically inactive C75U mutant and find significant differences in global structure and Mg(2+) affinity for both their precursor and product forms. Significantly, the Mg(2+) affinity of the C75 wild-type is 1.7-2.1-fold lower than that of the C75U mutant, in accord with the notion that C75 is essential for a catalytic conformational change that leads to a decrease in the local divalent metal ion affinity and release of a catalytic metal. Thus, a consistent picture emerges in which divalent metal ions and RNA functional groups are intimately intertwined in affecting structural dynamics and catalysis in the HDV ribozyme.     

RNA tertiary interactions mediate native collapse of a bacterial group I ribozyme

Large RNAs collapse into compact intermediates in the presence of counterions before folding to the native state. We previously found that collapse of a bacterial group I ribozyme correlates with the formation of helices within the ribozyme core, but occurs at Mg2+ concentrations too low to support stable tertiary structure and catalytic activity. Here, using small-angle X-ray scattering, we show that Mg2+-induced collapse is a cooperative folding transition that can be fit by a two-state model. The Mg2+ dependence of collapse is similar to the Mg2+ dependence of helix assembly measured by partial ribonuclease T1 digestion and of an unfolding transition measured by UV hypochromicity. The correspondence between multiple probes of RNA structure further supports a two-state model. A mutation that disrupts tertiary contacts between the L9 tetraloop and its helical receptor destabilized the compact state by 0.8 kcal/mol, while mutations in the central triplex were less destabilizing. These results show that native tertiary interactions stabilize the compact folding intermediates under conditions in which the RNA backbone remains accessible to solvent.     

Activation of PKR by RNA misfolding: HDV ribozyme dimers activate PKR

Protein Kinase R (PKR), the double-stranded RNA (dsRNA)-activated protein kinase, plays important roles in innate immunity. Previous studies have shown that PKR is activated by long stretches of dsRNA, RNA pseudoknots, and certain single-stranded RNAs; however, regulation of PKR by RNAs with globular tertiary structure has not been reported. In this study, the HDV ribozyme is used as a model of a mostly globular RNA. In addition to a catalytic core, the ribozyme contains a peripheral 13-bp pairing region (P4), which, upon shortening, affects neither the catalytic activity of the ribozyme nor its ability to crystallize. We report that the HDV ribozyme sequence alone can activate PKR. To elucidate the RNA structural basis for this, we prepared a number of HDV variants, including those with shortened or lengthened P4 pairing regions, with the anticipation that lengthening the P4 extension would yield a more potent activator since it would offer more base pairs of dsRNA. Surprisingly, the variant with a shortened P4 was the most potent activator. Through native gel mobility and enzymatic structure mapping experiments we implicate misfolded HDV ribozyme dimers as the PKR-activating species, and show that the shortened P4 leads to enhanced occupancy of the RNA dimer. These observations have implications for how RNA misfolding relates to innate immune response and human disease.