Concerted kinetic folding of a multidomain ribozyme with a disrupted loop-receptor interaction

General Strategies for RNA X-ray Crystallography

An extremely small proportion of the X-ray crystal structures deposited in the Protein Data Bank are of RNA or RNA-protein complexes. This is due to three main obstacles to the successful determination of RNA structure: (1) low yields of pure, properly folded RNA; (2) difficulty creating crystal contacts due to low sequence diversity; and (3) limited methods for phasing. Various approaches have been developed to address these obstacles, such as native RNA purification, engineered crystallization modules, and incorporation of proteins to assist in phasing. In this review, we will discuss these strategies and provide examples of how they are used in practice.

DEAD-box protein CYT-19 is activated by exposed helices in a group I intron RNA

DEAD-box proteins are nonprocessive RNA helicases and can function as RNA chaperones, but the mechanisms of their chaperone activity remain incompletely understood. The Neurospora crassa DEAD-box protein CYT-19 is a mitochondrial RNA chaperone that promotes group I intron splicing and has been shown to resolve misfolded group I intron structures, allowing them to refold. Building on previous results, here we use a series of tertiary contact mutants of the Tetrahymena group I intron ribozyme to demonstrate that the efficiency of CYT-19-mediated unfolding of the ribozyme is tightly linked to global RNA tertiary stability. Efficient unfolding of destabilized ribozyme variants is accompanied by increased ATPase activity of CYT-19, suggesting that destabilized ribozymes provide more productive interaction opportunities. The strongest ATPase stimulation occurs with a ribozyme that lacks all five tertiary contacts and does not form a compact structure, and small-angle X-ray scattering indicates that ATPase activity tracks with ribozyme compactness. Further, deletion of three helices that are prominently exposed in the folded structure decreases the ATPase stimulation by the folded ribozyme. Together, these results lead to a model in which CYT-19, and likely related DEAD-box proteins, rearranges complex RNA structures by preferentially interacting with and unwinding exposed RNA secondary structure. Importantly, this mechanism could bias DEAD-box proteins to act on misfolded RNAs and ribonucleoproteins, which are likely to be less compact and more dynamic than their native counterparts.

Mapping the kinetic barriers of a Large RNA molecule's folding landscape

The folding of linear polymers into discrete three-dimensional structures is often required for biological function. The formation of long-lived intermediates is a hallmark of the folding of large RNA molecules due to the ruggedness of their energy landscapes. The precise thermodynamic nature of the barriers (whether enthalpic or entropic) that leads to intermediate formation is still poorly characterized in large structured RNA molecules. A classic approach to analyzing kinetic barriers are temperature dependent studies analyzed with Eyring's transition state theory. We applied Eyring's theory to time-resolved hydroxyl radical (•OH) footprinting kinetics progress curves collected at eight temperature from 21.5 °C to 51 °C to characterize the thermodynamic nature of folding intermediate formation for the Mg(2+)-mediated folding of the Tetrahymena thermophila group I ribozyme. A common kinetic model configuration describes this RNA folding reaction over the entire temperature range studied consisting of primary (fast) transitions to misfolded intermediates followed by much slower secondary transitions, consistent with previous studies. Eyring analysis reveals that the primary transitions are moderate in magnitude and primarily enthalpic in nature. In contrast, the secondary transitions are daunting in magnitude and entropic in nature. The entropic character of the secondary transitions is consistent with structural rearrangement of the intermediate species to the final folded form. This segregation of kinetic control reveals distinctly different molecular mechanisms during the two stages of RNA folding and documents the importance of entropic barriers to defining rugged RNA folding landscapes.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

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(2+)-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.

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.

RNA molecules with conserved catalytic cores but variable peripheries fold along unique energetically optimized pathways

Functional and kinetic constraints must be efficiently balanced during the folding process of all biopolymers. To understand how homologous RNA molecules with different global architectures fold into a common core structure we determined, under identical conditions, the folding mechanisms of three phylogenetically divergent group I intron ribozymes. These ribozymes share a conserved functional core defined by topologically equivalent tertiary motifs but differ in their primary sequence, size, and structural complexity. Time-resolved hydroxyl radical probing of the backbone solvent accessible surface and catalytic activity measurements integrated with structural-kinetic modeling reveal that each ribozyme adopts a unique strategy to attain the conserved functional fold. The folding rates are not dictated by the size or the overall structural complexity, but rather by the strength of the constituent tertiary motifs which, in turn, govern the structure, stability, and lifetime of the folding intermediates. A fundamental general principle of RNA folding emerges from this study: The dominant folding flux always proceeds through an optimally structured kinetic intermediate that has sufficient stability to act as a nucleating scaffold while retaining enough conformational freedom to avoid kinetic trapping. Our results also suggest a potential role of naturally selected peripheral A-minor interactions in balancing RNA structural stability with folding efficiency.

Enhanced specificity against misfolding in a thermostable mutant of the Tetrahymena ribozyme

Structured RNAs encode native conformations that are more stable than the vast ensembles of alternative conformations, but how this specificity is evolved is incompletely understood. Here we show that a variant of the Tetrahymena group I intron ribozyme that was generated previously by in vitro selection for enhanced thermostability also displays modestly enhanced specificity against a stable misfolded structure that is globally similar to the native state, despite the absence of selective pressure to increase the energy gap between these structures. The enhanced specificity for native folding arises from mutations in two nucleotides that are close together in space in the native structure, and additional experiments show that these two mutations do not affect the stability of the misfolded conformation relative to the largely unstructured transition state ensemble for interconversion between the native and misfolded conformers. Thus, they selectively stabilize the native state, presumably by strengthening a local tertiary contact network that cannot form in the misfolded conformation. The stabilization is larger in the presence of the peripheral element P5abc, suggesting that cooperative tertiary structure formation plays a key role in the enhanced stability. The increased specificity in the absence of explicit selection suggests that the large energy gap in the wild-type RNA may have arisen analogously, a consequence of selective pressure for stability of the functional structure. More generally, the structural rigidity and intricate networks of contacts in structured RNAs may allow them to evolve substantial structural specificity without explicit negative selection, even against closely related alternative structures.

Taming free energy landscapes with RNA chaperones

Many non-coding RNAs fold into complex three-dimensional structures, yet the self-assembly of RNA structure is hampered by mispairing, weak tertiary interactions, electrostatic barriers, and the frequent requirement that the 5' and 3' ends of the transcript interact. This rugged free energy landscape for RNA folding means that some RNA molecules in a population rapidly form their native structure, while many others become kinetically trapped in misfolded conformations. Transient binding of RNA chaperone proteins destabilize misfolded intermediates and lower the transition states between conformations, producing a smoother landscape that increases the rate of folding and the probability that a molecule will find the native structure. DEAD-box proteins couple the chemical potential of ATP hydrolysis with repetitive cycles of RNA binding and release, expanding the range of conditions under which they can refold RNA structures.

Compact intermediates in RNA folding

Large noncoding RNAs fold into their biologically functional structures via compact yet disordered intermediates, which couple the stable secondary structure of the RNA with the emerging tertiary fold. The specificity of the collapse transition, which coincides with the assembly of helical domains, depends on RNA sequence and counterions. It determines the specificity of the folding pathways and the magnitude of the free energy barriers to the ensuing search for the native conformation. By coupling helix assembly with nascent tertiary interactions, compact folding intermediates in RNA also play a crucial role in ligand binding and RNA-protein recognition.

Multiple unfolding events during native folding of the Tetrahymena group I ribozyme

Despite the ubiquitous nature of misfolded intermediates in RNA folding, little is known about their physical properties or the folding transitions that allow them to continue folding productively. Folding of the Tetrahymena group I ribozyme includes sequential accumulation of two intermediates, termed I(trap) and misfolded (M). Here, we probe the structure and folding transition of I(trap) and compare them to those of M. Hydroxyl radical and dimethyl sulfate footprinting show that both I(trap) and M are extensively structured and crudely resemble the native RNA. However, regions of the core P3-P8 domain are more exposed to solvent in I(trap) than in M. I(trap) rearranges to continue folding nearly 1000-fold faster than M, and urea accelerates folding of I(trap) much less than M. Thus, the rate-limiting transition from I(trap) requires a smaller increase in exposed surface. Mutations that disrupt peripheral tertiary contacts give large and nearly uniform increases in re-folding of M, whereas the same mutations give at most modest increases in folding from I(trap). Intriguingly, mutations within the peripheral element P5abc give 5- to 10-fold accelerations in escape from I(trap), whereas ablation of P13, which lies on the opposite surface in the native structure, near the P3-P8 domain, has no effect. Thus, the unfolding required from I(trap) appears to be local, whereas the unfolding of M appears to be global. Further, the modest effects from several mutations suggest that there are multiple pathways for escape from I(trap) and that escape is aided by loosening nearby native structural constraints, presumably to facilitate local movements of nucleotides or segments that have not formed native contacts. Overall, these and prior results suggest a model in which the global architecture and peripheral interactions of the RNA are achieved relatively early in folding. Multiple folding and re-folding events occur on the predominant pathway to the native state, with increasing native core interactions and cooperativity as folding progresses.

A complex assembly landscape for the 30S ribosomal subunit

The ribosome is a complex macromolecular machine responsible for protein synthesis in the cell. It consists of two subunits, each of which contains both RNA and protein components. Ribosome assembly is subject to intricate regulatory control and is aided by a multitude of assembly factors in vivo, but can also be carried out in vitro. The details of the assembly process remain unknown even in the face of atomic structures of the entire ribosome and after more than three decades of research. Some of the earliest research on ribosome assembly produced the Nomura assembly map of the small subunit, revealing a hierarchy of protein binding dependencies for the 20 proteins involved and suggesting the possibility of a single intermediate. Recent work using a combination of RNA footprinting and pulse-chase quantitative mass spectrometry paints a picture of small subunit assembly as a dynamic and varied landscape, with sequential and hierarchical RNA folding and protein binding events finally converging on complete subunits. Proteins generally lock tightly into place in a 5' to 3' direction along the ribosomal RNA, stabilizing transient RNA conformations, while RNA folding and the early stages of protein binding are initiated from multiple locations along the length of the RNA.

Energy barriers, pathways, and dynamics during folding of large, multidomain RNAs

Large, multidomain RNA molecules are generally thought to fold following multiple pathways down rugged landscapes populated with intermediates and traps. A challenge to understanding RNA folding reactions is the complex relationships that exist between the structure of the RNA and its folding landscape. The identification of intermediate species that populate folding landscapes and characterization of elements of their structures are the key components to solving the RNA folding problem. This review explores recent studies that characterize the dominant pathways by which RNA folds, structural and dynamic features of intermediates that populate the folding landscape, and the energy barriers that separate the distinct steps of the folding process.

Direct measurement of tertiary contact cooperativity in RNA folding

All structured biological macromolecules must overcome the thermodynamic folding problem to populate a unique functional state among a vast ensemble of unfolded and alternate conformations. The exploration of cooperativity in protein folding has helped reveal and distinguish the underlying mechanistic solutions to this folding problem. Analogous dissections of RNA tertiary stability remain elusive, however, despite the central biological importance of folded RNA molecules and the potential to reveal fundamental properties of structured macromolecules via comparisons of protein and RNA folding. We report a direct quantitative measure of tertiary contact cooperativity in a folded RNA. We precisely measured the stability of an independently folding P4-P6 domain from the Tetrahymena thermophila group I intron by single molecule fluorescence resonance energy transfer (smFRET). Using wild-type and mutant RNAs, we found that cooperativity between the two tertiary contacts enhances P4-P6 stability by 3.2 +/- 0.2 kcal/mol.

Tertiary interactions determine the accuracy of RNA folding

RNAs must fold into unique three-dimensional structures to function in the cell, but how each polynucleotide finds its native structure is not understood. To investigate whether the stability of the tertiary structure determines the speed and accuracy of RNA folding, docking of a tetraloop with its receptor in a bacterial group I ribozyme was perturbed by site-directed mutagenesis. Disruption of the tetraloop or its receptor destabilizes tertiary interactions throughout the ribozyme by 2-3 kcal/mol, demonstrating that tertiary interactions form cooperatively in the transition from a native-like intermediate to the native state. Nondenaturing PAGE and RNase T1 digestion showed that base pairs form less homogeneously in the mutant RNAs during the transition from the unfolded state to the intermediate. Thus, tertiary interactions between helices bias the ensemble of secondary structures toward native-like conformations. Time-resolved hydroxyl radical footprinting showed that the wild-type ribozyme folds completely within 5-20 ms. By contrast, only 40-60% of a tetraloop mutant ribozyme folds in 30-40 ms, with the remainder folding in 30-200 s via nonnative intermediates. Therefore, destabilization of tetraloop-receptor docking introduces an alternate folding pathway in the otherwise smooth energy landscape of the wild-type ribozyme. Our results show that stable tertiary structure increases the flux through folding pathways that lead directly and rapidly to the native structure.

RNA misfolding and the action of chaperones

RNA folds to a myriad of three-dimensional structures and performs an equally diverse set of functions. The ability of RNA to fold and function in vivo is all the more remarkable because, in vitro, RNA has been shown to have a strong propensity to adopt misfolded, non-functional conformations. A principal factor underlying the dominance of RNA misfolding is that local RNA structure can be quite stable even in the absence of enforcing global tertiary structure. This property allows non-native structure to persist, and it also allows native structure to form and stabilize non-native contacts or non-native topology. In recent years it has become clear that one of the central reasons for the apparent disconnect between the capabilities of RNA in vivo and its in vitro folding properties is the presence of RNA chaperones, which facilitate conformational transitions of RNA and therefore mitigate the deleterious effects of RNA misfolding. Over the past two decades, it has been demonstrated that several classes of non-specific RNA binding proteins possess profound RNA chaperone activity in vitro and when overexpressed in vivo, and at least some of these proteins appear to function as chaperones in vivo. More recently, it has been shown that certain DExD/H-box proteins function as general chaperones to facilitate folding of group I and group II introns. These proteins are RNA-dependent ATPases and have RNA helicase activity, and are proposed to function by using energy from ATP binding and hydrolysis to disrupt RNA structure and/or to displace proteins from RNA-protein complexes. This review outlines experimental studies that have led to our current understanding of the range of misfolded RNA structures, the physical origins of RNA misfolding, and the functions and mechanisms of putative RNA chaperone proteins.

Kinetic redistribution of native and misfolded RNAs by a DEAD-box chaperone

DExD/H-box proteins are ubiquitously involved in RNA-mediated processes and use ATP to accelerate conformational changes in RNA. However, their mechanisms of action, and what determines which RNA species are targeted, are not well understood. Here we show that the DExD/H-box protein CYT-19, a general RNA chaperone, mediates ATP-dependent unfolding of both the native conformation and a long-lived misfolded conformation of a group I catalytic RNA with efficiencies that depend on the stabilities of the RNA species but not on specific structural features. CYT-19 then allows the RNA to refold, changing the distribution from equilibrium to kinetic control. Because misfolding is favoured kinetically, conditions that allow unfolding of the native RNA yield large increases in the population of misfolded species. Our results suggest that DExD/H-box proteins act with sufficient breadth and efficiency to allow structured RNAs to populate a wider range of conformations than would be present at equilibrium. Thus, RNAs may face selective pressure to stabilize their active conformations relative to inactive ones to avoid significant redistribution by DExD/H-box proteins. Conversely, RNAs whose functions depend on forming multiple conformations may rely on DExD/H-box proteins to increase the populations of less stable conformations, thereby increasing their overall efficiencies.

Distinct contribution of electrostatics, initial conformational ensemble, and macromolecular stability in RNA folding

We distinguish the contribution of the electrostatic environment, initial conformational ensemble, and macromolecular stability on the folding mechanism of a large RNA using a combination of time-resolved "Fast Fenton" hydroxyl radical footprinting and exhaustive kinetic modeling. This integrated approach allows us to define the folding landscape of the L-21 Tetrahymena thermophila group I intron structurally and kinetically from its earliest steps with unprecedented accuracy. Distinct parallel pathways leading the RNA to its native form upon its Mg(2+)-induced folding are observed. The structures of the intermediates populating the pathways are not affected by variation of the concentration and type of background monovalent ions (electrostatic environment) but are altered by a mutation that destabilizes one domain of the ribozyme. Experiments starting from different conformational ensembles but folding under identical conditions show that whereas the electrostatic environment modulates molecular flux through different pathways, the initial conformational ensemble determines the partitioning of the flux. This study showcases a robust approach for the development of kinetic models from collections of local structural probes.

Deletion of the P5abc peripheral element accelerates early and late folding steps of the Tetrahymena group I ribozyme

The P5abc peripheral element stabilizes the Tetrahymena group I ribozyme and enhances its catalytic activity. Despite its beneficial effects on the native structure, prior studies have shown that early formation of P5abc structure during folding can slow later folding steps. Here we use a P5abc deletion variant E(deltaP5abc) to systematically probe the role of P5abc throughout tertiary folding. Time-resolved hydroxyl radical footprinting shows that E(deltaP5abc) forms its earliest stable tertiary structure on the millisecond time scale, approximately 5-fold faster than the wild-type ribozyme, and stable structure spreads throughout E(deltaP5abc) in seconds. Nevertheless, activity measurements show that the earliest detectable formation of native E(deltaP5abc) ribozyme is much slower (approximately 0.6 min(-1)), in a manner similar to that of the wild type. Also similar, only a small fraction of E(deltaP5abc) attains the native state on this time scale under standard conditions at 25 degrees C, whereas the remainder misfolds; footprinting experiments show that the misfolded conformer shares structural features with the long-lived misfolded conformer of the wild-type ribozyme. Thus, P5abc does not have a large overall effect on the rate-limiting step(s) along this pathway. However, once misfolded, E(deltaP5abc) refolds to the native state 80-fold faster than the wild-type ribozyme and is less accelerated by urea, indicating that P5abc stabilizes the misfolded structure relative to the less-ordered transition state for refolding. Together, the results suggest that, under these conditions, even the earliest tertiary folding intermediates of the wild-type ribozyme represent misfolded species and that P5abc is principally a liability during the tertiary folding process.

Nonspecific binding to structured RNA and preferential unwinding of an exposed helix by the CYT-19 protein, a DEAD-box RNA chaperone

We explore the interactions of CYT-19, a DExD/H-box protein that functions in folding of group I RNAs, with a well characterized misfolded species of the Tetrahymena ribozyme. Consistent with its function, CYT-19 accelerates refolding of the misfolded RNA to its native state. Unexpectedly, CYT-19 performs another reaction much more efficiently; it unwinds the 6-bp P1 duplex formed between the ribozyme and its oligonucleotide substrate. Furthermore, CYT-19 performs this reaction 50-fold more efficiently than it unwinds the same duplex free in solution, suggesting that it forms additional interactions with the ribozyme, most likely using a distinct RNA binding site from the one responsible for unwinding. This site can apparently bind double-stranded RNA, as attachment of a simple duplex adjacent to P1 recapitulates much of the activation provided by the ribozyme. Unwinding the native P1 duplex does not accelerate refolding of the misfolded ribozyme, implying that CYT-19 can disrupt multiple contacts on the RNA, consistent with its function in folding of multiple RNAs. Further experiments showed that the P1 duplex unwinding activity is virtually the same whether the ribozyme is misfolded or native but is abrogated by formation of tertiary contacts between the P1 duplex and the body of the ribozyme. Together these results suggest a mechanism for CYT-19 and other general DExD/H-box RNA chaperones in which the proteins bind to structured RNAs and efficiently unwind loosely associated duplexes, which biases the proteins to disrupt nonnative base pairs and gives the liberated strands an opportunity to refold.

The paradoxical behavior of a highly structured misfolded intermediate in RNA folding

Like many structured RNAs, the Tetrahymena group I ribozyme is prone to misfolding. Here we probe a long-lived misfolded species, referred to as M, and uncover paradoxical aspects of its structure and folding. Previous work indicated that a non-native local secondary structure, termed alt P3, led to formation of M during folding in vitro. Surprisingly, hydroxyl radical footprinting, fluorescence measurements with site-specifically incorporated 2-aminopurine, and functional assays indicate that the native P3, not alt P3, is present in the M state. The paradoxical behavior of alt P3 presumably arises because alt P3 biases folding toward M, but, after commitment to this folding pathway and before formation of M, alt P3 is replaced by P3. Further, structural and functional probes demonstrate that the misfolded ribozyme contains extensive native structure, with only local differences between the two states, and the misfolded structure even possesses partial catalytic activity. Despite the similarity of these structures, re-folding of M to the native state is very slow and is strongly accelerated by urea, Na+, and increased temperature and strongly impeded by Mg2+ and the presence of native peripheral contacts. The paradoxical observations of extensive native structure within the misfolded species but slow conversion of this species to the native state are readily reconciled by a model in which the misfolded state is a topological isomer of the native state, and computational results support the feasibility of this model. We speculate that the complex topology of RNA secondary structures and the inherent rigidity of RNA helices render kinetic traps due to topological isomers considerably more common for RNA than for proteins.

Perturbation of the hierarchical folding of a large RNA by the destabilization of its Scaffold's tertiary structure

The P4-P6 domain serves as a scaffold against which the periphery and catalytic core organize and fold during Mg2+-mediated folding of the Tetrahymena thermophila ribozyme. The most prominent structural motif of the P4-P6 domain is the tetraloop-tetraloop receptor interaction which "clamps" the distal parts of its hairpin-like structure. Destabilization of the tertiary structure of the P4-P6 domain by perturbation of the tetraloop-tetraloop receptor interaction alters the Mg2+-mediated folding pathway. The folding hierarchy of P5c approximately P4-P6 > periphery > catalytic core that is a striking attribute of the folding of the wild-type RNA is abolished. The initial steps in folding of the mutant RNA are > or =50-fold faster than those of the wild-type ribozyme with the earliest observed tertiary contacts forming around regions known to specifically bind Mg2+. The interaction between the mutant tetraloop and the tetraloop receptor appears coincidently with slowly forming catalytic core tertiary contacts. Thus, the stability conferred upon the P4-P6 domain by the tetraloop-tetraloop receptor interaction dictates the preferred folding pathway by stabilizing an early intermediate. A sub-denaturing concentration of urea diminishes the early barrier to folding the wild-type ribozyme along with complex effects on the subsequent steps of folding the wild-type and mutant RNA.

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.

Structural specificity conferred by a group I RNA peripheral element

Like proteins, structured RNAs must specify a native conformation that is more stable than all other possible conformations. Local structure is much more stable for RNA than for protein, so it is likely that the principal challenge for RNA is to stabilize the native structure relative to misfolded and partially folded intermediates rather than unfolded structures. Many structured RNAs contain peripheral structural elements, which surround the core elements. Although it is clear that peripheral elements stabilize structure within RNAs that contain them, it has not yet been explored whether they specifically stabilize the native states relative to alternative folds. A two-piece version of the group I intron RNA from Tetrahymena is used here to show that the peripheral element P5abc binds to the native conformation of the rest of the RNA 50,000 times more tightly than it binds to a long-lived misfolded conformation. Thus, P5abc stabilizes the native conformation by approximately 6 kcal/mol relative to this misfolded conformation. Further, activity measurements show that for the RNA lacking P5abc, the native conformation is only marginally preferred over the misfolded conformation (<0.5 kcal/mol), indicating that the peripheral structure of this RNA is required to achieve a significant thermodynamic preference for the native state. Such "structural specificity" may be a general function of RNA peripheral domains.

Predicting the secondary structures and tertiary interactions of 211 group I introns in IE subgroup

The large number of currently available group I intron sequences in the public databases provides opportunity for studying this large family of structurally complex catalytic RNA by large-scale comparative sequence analysis. In this study, the detailed secondary structures of 211 group I introns in the IE subgroup were manually predicted. The secondary structure-favored alignments showed that IE introns contain 14 conserved stems. The P13 stem formed by long-range base-pairing between P2.1 and P9.1 is conserved among IE introns. Sequence variations in the conserved core divide IE introns into three distinct minor subgroups, namely IE1, IE2 and IE3. Co-variation of the peripheral structural motifs with core sequences supports that the peripheral elements function in assisting the core structure folding. Interestingly, host-specific structural motifs were found in IE2 introns inserted at S516 position. Competitive base-pairing is found to be conserved at the junctions of all long-range paired regions, suggesting a possible mechanism of establishing long-range base-pairing during large RNA folding. These findings extend our knowledge of IE introns, indicating that comparative analysis can be a very good complement for deepening our understanding of RNA structure and function in the genomic era.

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.

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(2+) at physiological concentrations, with both showing an Mg(2+) 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(2+); neither slower kinetics nor obvious Mg(2+) 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.

Multiple monovalent ion-dependent pathways for the folding of the L-21 Tetrahymena thermophila ribozyme

Synchrotron hydroxyl radical (*OH) footprinting is a technique that monitors the local changes in solvent accessibility of the RNA backbone on milliseconds to minutes time-scales. The Mg(2+)-dependent folding of the L-21 Sca 1 Tetrahymena thermophila ribozyme has been followed using this technique at an elevated concentration of monovalent ion (200 mM NaCl) and as a function of the initial annealing conditions and substrate. Previous studies conducted at low concentrations of monovalent ion displayed sequential folding of the P4-P6 domain, the peripheral helices and the catalytic core, with each protection displaying monophasic kinetics. For ribozyme annealed in buffer containing 200 mM NaCl and folded by the addition of 10 mM MgCl(2), multiple kinetic phases are observed for *OH protections throughout the ribozyme. The independently folding P4-P6 domain is the first to fold with its protections displaying 50-90% burst phase amplitudes. That the folding of P4-P6 within the ribozyme does not display the 100% burst phase of isolated P4-P6 at 200 mM NaCl shows that interactions with the remainder of the ribozyme impede this domain's folding. In addition, *OH protections constituting each side of a tertiary contact are not coincident in some cases, consistent with the formation of transient non-native interactions. While the peripheral contacts and triple helical scaffold exhibit substantial burst phases, the slowest protection to appear is J8/7 in the catalytic core, which displays a minimal burst amplitude and whose formation is coincident with the recovery of catalytic activity. The number of kinetic phases as well as their amplitudes and rates are different when the ribozyme is annealed in low-salt buffer and folded by the concomitant addition of monovalent and divalent cations. Annealed substrate changes the partitioning of the ribozyme among the multiple folding populations. These results provide a map of the early steps in the ribozyme's folding landscape and the degree to which the preferred pathways are dependent upon the initial reaction conditions.

Mispaired P3 region in the hierarchical folding pathway of the Tetrahymena ribozyme

BACKGROUND

The Tetrahymena group I ribozyme folds into a complex three-dimensional structure for performing catalytic reactions. The catalysis depends on its catalytic core consisting of two helical domains, P4-P6 and P3-P7, connected by single stranded regions. In the folding process, most of this ribozyme folds in a hierarchical manner in which a kinetically stable intermediate determines the overall folding rate.

RESULTS

Although the nature of this intermediate has not yet been elucidated, a mispaired P3 stem (alt-P3) appears a likely candidate. To examine the effects of the alt-P3 structure on the kinetic and thermodynamic properties of the active structure of the ribozyme or its P3-P7 domain formation, we prepared and analysed variant ribozymes in which relative stabilities of the original P3 and alt-P3 structure were altered systematically.

CONCLUSION

The results indicate that the alt-P3 structure is not the major rate-limiting factor in the folding process.

Relationship between the self-splicing activity and the solidity of the master domain of the Tetrahymena group I ribozyme

The highly conserved P3-P7 domain of the Group I intron ribozymes is known to contain essential elements, such as the binding site for the cofactor guanosine, required for conducting the splicing reaction. We investigated the domain of the Tetrahymena intron ribozyme and its variants in order to clarify the relationship between its stability and function. We found that the destabilization of the P3-P7 domain facilitates the active structure formation at high magnesium ion concentrations where the formation is retarded for the wild type. The destabilized domain also increases K(GTP)(m) although this can be compensated by increasing the concentration of Mg(2+), indicating that the stable domain is required for establishing a tight guanosine binding site. The results suggest that the stability of the domain affects the rate-limiting step in the RNA folding pathway and also regulates the efficiency of the splicing reaction.

Productive folding to the native state by a group II intron ribozyme

Group II introns are large catalytic RNA molecules that fold into compact structures essential for the catalysis of splicing and intron mobility reactions. Despite a growing body of information on the folded state of group II introns at equilibrium, there is currently no information on the folding pathway and little information on the ionic requirements for folding. Folding isotherms were determined by hydroxyl radical footprinting for the 32 individual protections that are distributed throughout a group II intron ribozyme derived from intron ai5gamma. The isotherms span a similar range of Mg(2+) concentrations and share a similar index of cooperativity. Time-resolved hydroxyl radical footprinting studies show that all regions of the ribozyme fold slowly and with remarkable synchrony into a single catalytically active structure at a rate comparable to those of other ribozymes studied thus far. The rate constants for the formation of tertiary contacts and recovery of catalytic activity are identical within experimental error. Catalytic activity analyses in the presence of urea provide no evidence that the slow folding of the ai5gamma intron is attributable to the presence of unproductive kinetic traps along the folding pathway. Taken together, the data suggest that the rate-limiting step for folding of group II intron ai5gamma occurs early along the reaction pathway. We propose that this behavior resembles protein folding that is limited in rate by high contact order, or the need to form key tertiary interactions from partners that are located far apart in the primary or secondary structure.

Modular construction of a tertiary RNA structure: the specificity domain of the Bacillus subtilis RNase P RNA

The structure of the specificity domain (S-domain) of the Bacillus subtilis RNase P RNA has been proposed to be composed of a core and a buttress module, analogous to the bipartite structure of the P4-P6 domain of the Tetrahymena group I ribozyme. The core module is the functional unit of the S-domain and contains the binding site for the T stem-loop of a tRNA. The buttress module provides structural stability to the core module and consists of a GA3 tetraloop and its receptor. To explicitly test the hypothesis that modular construction can describe the structure of the S-domain and is a useful RNA design strategy, we analyzed the equilibrium folding and substrate binding of three classes of S-domain mutants. Addition or deletion of a base pair in the helical linker region between the modules only modestly destabilizes the tertiary structure. tRNA binding selectivity is affected in one but not in two other mutants of this class. Elimination of the GA3 tetraloop-receptor interactions significantly destabilizes the core module and results in the loss of tRNA binding selectivity. Replacing the buttress module with that of a homologous RNase P RNA maintains the tRNA binding selectivity. Overall, we have observed that the linker regions between the two modules can tolerate moderate structural changes and that the buttress modules can be shuffled between homologous S-domains. These results suggest that it is feasible to design an RNA using a buttress module to stabilize a functional module.

Beyond kinetic traps in RNA folding

Large RNAs often have rugged folding energy landscapes that result in severe misfolding and slow folding kinetics. Several interdependent parameters that contribute to misfolding are now well understood and examples of large RNAs and ribonucleoproteins that avoid kinetic traps have been reported. These advances have facilitated the exploration of fundamental RNA folding processes that were previously inaccessible.

A deteriorated triple-helical scaffold accelerates formation of the Tetrahymena ribozyme active structure

The Tetrahymena group I ribozyme requires a hierarchical folding process to form its correct three-dimensional structure. Ribozyme activity depends on the catalytic core consisting of two domains, P4-P6 and P3-P7, connected by a triple-helical scaffold. The folding proceeds in the following order: (i) fast folding of the P4-P6 domain, (ii) slow folding of the P3-P7 domain, and (iii) structure rearrangement to form the active ribozyme structure. The third step is believed to directly determine the conformation of the active catalytic domain, but as yet the precise mechanisms remain to be elucidated. To investigate the folding kinetics of this step, we analyzed mutant ribozymes having base substitution(s) in the triple-helical scaffold and found that disruption of the scaffold at A105G results in modest slowing of the P3-P7 folding (1.9-fold) and acceleration of step (iii) by 5.9-fold. These results suggest that disruption or destabilization of the scaffold is a normal component in the formation process of the active structure of the wild type ribozyme.

Formation of a GNRA tetraloop in P5abc can disrupt an interdomain interaction in the Tetrahymena group I ribozyme

The secondary structure of a truncated P5abc subdomain (tP5abc, a 56-nucleotide RNA) of the Tetrahymena thermophila group I intron ribozyme changes when its tertiary structure forms. We have now used heteronuclear NMR spectroscopy to determine its conformation in solution. The tP5abc RNA that contains only secondary structure is extended compared with the tertiary folded form; both forms coexist in slow chemical exchange (the interconversion rate constant is slower than 1 s(-1)) in the presence of magnesium. Kinetic experiments have shown that tertiary folding of the P5abc subdomain is one of the earliest folding transitions in the group I intron ribozyme, and that it leads to a metastable misfolded intermediate. Previous mutagenesis studies suggest that formation of the extended P5abc structure described here destabilize a misfolded intermediate. This study shows that the P5abc RNA subdomain containing a GNRA tetraloop in P5c (in contrast to the five-nucleotide loop P5c in the tertiary folded ribozyme) can disrupt the base-paired interdomain (P14) interaction between P5c and P2.