New loop-loop tertiary interactions in self-splicing introns of subgroup IC and ID: a complete 3D model of the Tetrahymena thermophila ribozyme

The long-range P3 helix of the Tetrahymena ribozyme is disrupted during folding between the native and misfolded conformations

RNAs are prone to misfolding, but how misfolded structures are formed and resolved remains incompletely understood. The Tetrahymena group I intron ribozyme folds in vitro to a long-lived misfolded conformation (M) that includes extensive native structure but is proposed to differ in topology from the native state (N). A leading model predicts that exchange of the topologies requires unwinding of the long-range, core helix P3, despite the presence of P3 in both conformations. To test this model, we constructed 16 mutations to strengthen or weaken P3. Catalytic activity and in-line probing showed that nearly all of the mutants form the M state before folding to N. The P3-weakening mutations accelerated refolding from M (3- to 30-fold) and the P3-strengthening mutations slowed refolding (6- to 1400-fold), suggesting that P3 indeed unwinds transiently. Upon depletion of Mg(2+), the mutations had analogous effects on unfolding from N to intermediates that subsequently fold to M. The magnitudes for the P3-weakening mutations were larger than in refolding from M, and small-angle X-ray scattering showed that the ribozyme expands rapidly to intermediates from which P3 is disrupted subsequently. These results are consistent with previous results indicating unfolding of native peripheral structure during refolding from M, which probably permits rearrangement of the core. Together, our results demonstrate that exchange of the native and misfolded conformations requires loss of a core helix in addition to peripheral structure. Further, the results strongly suggest that misfolding arises from a topological error within the ribozyme core, and a specific topology is proposed.

RNA 3D structure prediction by using a coarse-grained model and experimental data

RNAs form complex secondary and three-dimensional structures, and their biological functions highly rely on their structures and dynamics. Here we developed a general coarse-grained framework for RNA 3D structure prediction. A new, hybrid coarse-grained model that explicitly describes the electrostatics and hydrogen-bond interactions has been constructed based on experimental structural statistics. With the simulated annealing simulation protocol, several RNAs of less than 30-nt were folded to within 4.0 Å of the native structures. In addition, with limited restraints on Watson-Crick basepairing based on the data from NMR spectroscopy and small-angle X-ray scattering (SAXS) information, the current model was able to characterize the complex tertiary structures of large size RNAs, such as 5S ribosome and U2/U6 snRNA. We also demonstrated that the pseudoknot structure was better captured when the coordinating Mg(2+) cations and limited basepairing restraints were included. The accuracy of our model has been compared favorably with other RNA structure prediction methods presented in the previous study of RNA-Puzzles. Therefore the coarse-grained model presented here offers a unique approach for accurate prediction and modeling of RNA structures.

The GA-minor submotif as a case study of RNA modularity, prediction, and design

Complex natural RNAs such as the ribosome, group I and group II introns, and RNase P exemplify the fact that three-dimensional (3D) RNA structures are highly modular and hierarchical in nature. Tertiary RNA folding typically takes advantage of a rather limited set of recurrent structural motifs that are responsible for controlling bends or stacks between adjacent helices. Herein, the GA minor and related structural motifs are presented as a case study to highlight several structural and folding principles, to gain further insight into the structural evolution of naturally occurring RNAs, as well as to assist the rational design of artificial RNAs.

Toward a molecular understanding of RNA remodeling by DEAD-box proteins

DEAD-box proteins are superfamily 2 helicases that function in all aspects of RNA metabolism. They employ ATP binding and hydrolysis to generate tight, yet regulated RNA binding, which is used to unwind short RNA helices non-processively and promote structural transitions of RNA and RNA-protein substrates. In the last few years, substantial progress has been made toward a detailed, quantitative understanding of the structural and biochemical properties of DEAD-box proteins. Concurrently, progress has been made toward a physical understanding of the RNA rearrangements and folding steps that are accelerated by DEAD-box proteins in model systems. Here, we review the recent progress on both of these fronts, focusing on the mitochondrial DEAD-box proteins Mss116 and CYT-19 and their mechanisms in promoting the splicing of group I and group II introns.

The right angle (RA) motif: a prevalent ribosomal RNA structural pattern found in group I introns

The right angle (RA) motif, previously identified in the ribosome and used as a structural module for nano-construction, is a recurrent structural motif of 13 nucleotides that establishes a 90° bend between two adjacent helices. Comparative sequence analysis was used to explore the sequence space of the RA motif within ribosomal RNAs in order to define its canonical sequence space signature. We investigated the sequence constraints associated with the RA signature using several artificial self-assembly systems. Thermodynamic and topological investigations of sequence variants associated with the RA motif in both minimal and expanded structural contexts reveal that the presence of a helix at the 3' end of the RA motif increases the thermodynamic stability and rigidity of the resulting three-helix junction domain. A search for the RA in naturally occurring RNAs as well as its experimental characterization led to the identification of the RA in groups IC1 and ID intron ribozymes, where it is suggested to play an integral role in stabilizing peripheral structural domains. The present study exemplifies the need of empirical analysis of RNA structural motifs for facilitating the rational design and structure prediction of RNAs.

A two-piece derivative of a group I intron RNA as a platform for designing self-assembling RNA templates to promote Peptide ligation

Multicomponent RNA-peptide complexes are attractive from the viewpoint of artificial design of functional biomacromolecular systems. We have developed self-folding and self-assembling RNAs that serve as templates to assist chemical ligation between two reactive peptides with RNA-binding capabilities. The design principle of previous templates, however, can be applied only to limited classes of RNA-binding peptides. In this study, we employed a two-piece derivative of a group I intron RNA from the Tetrahymena large subunit ribosomal RNA (LSU rRNA) as a platform for new template RNAs. In this group I intron-based self-assembling platform, modules for the recognition of substrate peptides can be installed independently from modules holding the platform structure. The new self-assembling platform allows us to expand the repertoire of substrate peptides in template RNA design.

Cooperative tertiary interaction network guides RNA folding

Noncoding RNAs form unique 3D structures, which perform many regulatory functions. To understand how RNAs fold uniquely despite a small number of tertiary interaction motifs, we mutated the major tertiary interactions in a group I ribozyme by single-base substitutions. The resulting perturbations to the folding energy landscape were measured using SAXS, ribozyme activity, hydroxyl radical footprinting, and native PAGE. Double- and triple-mutant cycles show that most tertiary interactions have a small effect on the stability of the native state. Instead, the formation of core and peripheral structural motifs is cooperatively linked in near-native folding intermediates, and this cooperativity depends on the native helix orientation. The emergence of a cooperative interaction network at an early stage of folding suppresses nonnative structures and guides the search for the native state. We suggest that cooperativity in noncoding RNAs arose from natural selection of architectures conducive to forming a unique, stable fold.

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.

DEAD-box proteins as RNA helicases and chaperones

DEAD-box proteins are ubiquitous in RNA-mediated processes and function by coupling cycles of ATP binding and hydrolysis to changes in affinity for single-stranded RNA. Many DEAD-box proteins use this basic mechanism as the foundation for a version of RNA helicase activity, efficiently separating the strands of short RNA duplexes in a process that involves little or no translocation. This activity, coupled with mechanisms to direct different DEAD-box proteins to their physiological substrates, allows them to promote RNA folding steps and rearrangements and to accelerate remodeling of RNA–protein complexes. This review will describe the properties of DEAD-box proteins as RNA helicases and the current understanding of how the energy from ATPase activity is used to drive the separation of RNA duplex strands. It will then describe how the basic biochemical properties allow some DEAD-box proteins to function as chaperones by promoting RNA folding reactions, with a focus on the self-splicing group I and group II intron RNAs.

Computational prediction of efficient splice sites for trans-splicing ribozymes

Group I introns have been engineered into trans-splicing ribozymes capable of replacing the 3'-terminal portion of an external mRNA with their own 3'-exon. Although this design makes trans-splicing ribozymes potentially useful for therapeutic application, their trans-splicing efficiency is usually too low for medical use. One factor that strongly influences trans-splicing efficiency is the position of the target splice site on the mRNA substrate. Viable splice sites are currently determined using a biochemical trans-tagging assay. Here, we propose a rapid and inexpensive alternative approach to identify efficient splice sites. This approach involves the computation of the binding free energies between ribozyme and mRNA substrate. We found that the computed binding free energies correlate well with the trans-splicing efficiency experimentally determined at 18 different splice sites on the mRNA of chloramphenicol acetyl transferase. In contrast, our results from the trans-tagging assay correlate less well with measured trans-splicing efficiency. The computed free energy components suggest that splice site efficiency depends on the following secondary structure rearrangements: hybridization of the ribozyme's internal guide sequence (IGS) with mRNA substrate (most important), unfolding of substrate proximal to the splice site, and release of the IGS from the 3'-exon (least important). The proposed computational approach can also be extended to fulfill additional design requirements of efficient trans-splicing ribozymes, such as the optimization of 3'-exon and extended guide sequences.

A role for a single-stranded junction in RNA binding and specificity by the Tetrahymena group I ribozyme

We have investigated the role of a single-stranded RNA junction, J1/2, that connects the substrate-containing P1 duplex to the remainder of the Tetrahymena group I ribozyme. Single-turnover kinetics, fluorescence anisotropy, and single-molecule fluorescence resonance energy transfer studies of a series of J1/2 mutants were used to probe the sequence dependence of the catalytic activity, the P1 dynamics, and the thermodynamics of docking of the P1 duplex into the ribozyme's catalytic core. We found that A29, the center A of three adenosine residues in J1/2, contributes 2 orders of magnitude to the overall ribozyme activity, and double-mutant cycles suggested that J1/2 stabilizes the docked state of P1 over the undocked state via a tertiary interaction involving A29 and the first base pair in helix P2 of the ribozyme, A31·U56. Comparative sequence analysis of this group I intron subclass suggests that the A29 interaction sets one end of a molecular ruler whose other end specifies the 5'-splice site and that this molecular ruler is conserved among a subclass of group I introns related to the Tetrahymena intron. Our results reveal substantial functional effects from a seemingly simple single-stranded RNA junction and suggest that junction sequences may evolve rapidly to provide important interactions in functional RNAs.

Attenuation of loop-receptor interactions with pseudoknot formation

RNA tetraloops can recognize receptors to mediate long-range interactions in stable natural RNAs. In vitro selected GNRA tetraloop/receptor interactions are usually more ‘G/C-rich’ than their ‘A/U-rich’ natural counterparts. They are not as widespread in nature despite comparable biophysical and chemical properties. Moreover, while AA, AC and GU dinucleotide platforms occur in natural GAAA/11 nt receptors, the AA platform is somewhat preferred to the others. The apparent preference for ‘A/U-rich’ GNRA/receptor interactions in nature might stem from an evolutionary adaptation to avoid folding traps at the level of the larger molecular context. To provide evidences in favor of this hypothesis, several riboswitches based on natural and artificial GNRA receptors were investigated in vitro for their ability to prevent inter-molecular GNRA/receptor interactions by trapping the receptor sequence into an alternative intra-molecular pseudoknot. Extent of attenuation determined by native gel-shift assays and co-transcriptional assembly is correlated to the G/C content of the GNRA receptor. Our results shed light on the structural evolution of natural long-range interactions and provide design principles for RNA-based attenuator devices to be used in synthetic biology and RNA nanobiotechnology.

A two-dimensional mutate-and-map strategy for non-coding RNA structure

Non-coding RNAs fold into precise base-pairing patterns to carry out critical roles in genetic regulation and protein synthesis, but determining RNA structure remains difficult. Here, we show that coupling systematic mutagenesis with high-throughput chemical mapping enables accurate base-pair inference of domains from ribosomal RNA, ribozymes and riboswitches. For a six-RNA benchmark that has challenged previous chemical/computational methods, this 'mutate-and-map' strategy gives secondary structures that are in agreement with crystallography (helix error rates, 2%), including a blind test on a double-glycine riboswitch. Through modelling of partially ordered states, the method enables the first test of an interdomain helix-swap hypothesis for ligand-binding cooperativity in a glycine riboswitch. Finally, the data report on tertiary contacts within non-coding RNAs, and coupling to the Rosetta/FARFAR algorithm gives nucleotide-resolution three-dimensional models (helix root-mean-squared deviation, 5.7 Å) of an adenine riboswitch. These results establish a promising two-dimensional chemical strategy for inferring the secondary and tertiary structures that underlie non-coding RNA behaviour.

Structure-function analysis from the outside in: long-range tertiary contacts in RNA exhibit distinct catalytic roles

The conserved catalytic core of the Tetrahymena group I ribozyme is encircled by peripheral elements. We have conducted a detailed structure-function study of the five long-range tertiary contacts that fasten these distal elements together. Mutational ablation of each of the tertiary contacts destabilizes the folded ribozyme, indicating a role of the peripheral elements in overall stability. Once folded, three of the five tertiary contact mutants exhibit defects in overall catalysis that range from 20- to 100-fold. These and the subsequent results indicate that the structural ring of peripheral elements does not act as a unitary element; rather, individual connections have distinct roles as further revealed by kinetic and thermodynamic dissection of the individual reaction steps. Ablation of P14 or the metal ion core/metal ion core receptor (MC/MCR) destabilizes docking of the substrate-containing P1 helix into tertiary interactions with the ribozyme's conserved core. In contrast, ablation of the L9/P5 contact weakens binding of the guanosine nucleophile by slowing its association, without affecting P1 docking. The P13 and tetraloop/tetraloop receptor (TL/TLR) mutations had little functional effect and small, local structural changes, as revealed by hydroxyl radical footprinting, whereas the P14, MC/MCR, and L9/P5 mutants show structural changes distal from the mutation site. These changes extended into regions of the catalytic core involved in docking or guanosine binding. Thus, distinct allosteric pathways couple the long-range tertiary contacts to functional sites within the conserved core. This modular functional specialization may represent a fundamental strategy in RNA structure-function interrelationships.

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.

Occurrence and characteristics of group 1 introns found at three different positions within the 28S ribosomal RNA gene of the dematiaceous Phialophora verrucosa: phylogenetic and secondary structural implications

BACKGROUND

Group 1 introns (ribozymes) are among the most ancient and have the broadest phylogenetic distribution among the known self-splicing ribozymes. Fungi are known to be rich in rDNA group 1 introns. In the present study, five sequences of the 28S ribosomal RNA gene (rDNA) regions of pathogenic dematiaceous Phialophora verrucosa were analyzed using PCR by site-specific primers and were found to have three insertions, termed intron-F, G and H, at three positions of the gene. We investigated the distribution of group 1 introns in this fungus by surveying 34 strains of P. verrucosa and seven strains of Phialophora americana as the allied species.

RESULTS

Intron-F's (inserted at L798 position) were found in 88% of P. verrucosa strains, while intron-G's (inserted at L1921) at 12% and intron-H's (inserted at L2563) at 18%. There was some correlation between intron distribution and geographic location. In addition, we confirmed that the three kinds of introns are group 1 introns from results of BLAST search, alignment analysis and Reverse Transcriptase-Polymerase Chain Reaction (RT-PCR). Prediction of secondary structures and phylogenetic analysis of intron sequences identified introns-F and G as belonging to subgroup IC1. In addition, intron-H was identified as IE.

CONCLUSION

The three intron insertions and their insertion position in the 28S rDNA allowed the characterization of the clinical and environmental isolates of P. verrucosa and P. americana into five genotypes. All subgroups of introns-F and G and intron-H were characterized and observed for the first time in both species.

Intermolecular interaction between a branching ribozyme and associated homing endonuclease mRNA

RNA tertiary interactions involving docking of GNRA (N; any base; R; purine) hairpin loops into helical stem structures on other regions of the same RNA are one of the most common RNA tertiary interactions. In this study, we investigated a tertiary association between a GAAA hairpin tetraloop in a small branching ribozyme (DiGIR1) and a receptor motif (HEG P1 motif) present in a hairpin structure on a separate mRNA molecule. DiGIR1 generates a 2', 5' lariat cap at the 5' end of its downstream homing endonuclease mRNA by catalysing a self-cleavage branching reaction at an internal processing site. Upon release, the 5' end of the mRNA forms a distinct hairpin structure termed HEG P1. Our biochemical data, in concert with molecular 3D modelling, provide experimental support for an intermolecular tetraloop receptor interaction between the L9 GAAA in DiGIR1 and a GNRA tetraloop receptor-like motif (UCUAAG-CAAGA) found within the HEG P1. The biological role of this interaction appears to be linked to the homing endonuclease expression by promoting post-cleavage release of the lariat capped mRNA. These findings add to our understanding of how protein-coding genes embedded in nuclear ribosomal DNA are expressed in eukaryotes and controlled by ribozymes.

A mutate-and-map strategy accurately infers the base pairs of a 35-nucleotide model RNA

We present a rapid experimental strategy for inferring base pairs in structured RNAs via an information-rich extension of classic chemical mapping approaches. The mutate-and-map method, previously applied to a DNA/RNA helix, systematically searches for single mutations that enhance the chemical accessibility of base-pairing partners distant in sequence. To test this strategy for structured RNAs, we have carried out mutate-and-map measurements for a 35-nt hairpin, called the MedLoop RNA, embedded within an 80-nt sequence. We demonstrate the synthesis of all 105 single mutants of the MedLoop RNA sequence and present high-throughput DMS, CMCT, and SHAPE modification measurements for this library at single-nucleotide resolution. The resulting two-dimensional data reveal visually clear, punctate features corresponding to RNA base pair interactions as well as more complex features; these signals can be qualitatively rationalized by comparison to secondary structure predictions. Finally, we present an automated, sequence-blind analysis that permits the confident identification of nine of the 10 MedLoop RNA base pairs at single-nucleotide resolution, while discriminating against all 1460 false-positive base pairs. These results establish the accuracy and information content of the mutate-and-map strategy and support its feasibility for rapidly characterizing the base-pairing patterns of larger and more complex RNA systems.

Predicting and modeling RNA architecture

A general approach for modeling the architecture of large and structured RNA molecules is described. The method exploits the modularity and the hierarchical folding of RNA architecture that is viewed as the assembly of preformed double-stranded helices defined by Watson-Crick base pairs and RNA modules maintained by non-Watson-Crick base pairs. Despite the extensive molecular neutrality observed in RNA structures, specificity in RNA folding is achieved through global constraints like lengths of helices, coaxiality of helical stacks, and structures adopted at the junctions of helices. The Assemble integrated suite of computer tools allows for sequence and structure analysis as well as interactive modeling by homology or ab initio assembly with possibilities for fitting within electronic density maps. The local key role of non-Watson-Crick pairs guides RNA architecture formation and offers metrics for assessing the accuracy of three-dimensional models in a more useful way than usual root mean square deviation (RMSD) values.

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.

A folding "framework structure" of Tetrahymena group I intron

We have published the dynamic extended folding (DEF) method, which is a RNA secondary structure prediction approach-to simulate the in vivo RNA co-transcriptional folding process. In order to verify the reliability of the method, we selected the X-ray-determined Tetrahymena group I intron as a sample to construct the framework of its folding secondary structure. Our prediction coincides well with the secondary structure predicted by T.R. Cech and the X-ray diffraction crystal structure determined by Lehnert V. Our results show that the DEF framework structure of Tetrahymena group I intron reflects its function sites in a concise and straightforward manner, and the scope of the simulation was expanded.

Roles of DEAD-box proteins in RNA and RNP Folding

RNAs and RNA-protein complexes (RNPs) traverse rugged energy landscapes as they fold to their native structures, and many continue to undergo conformational rearrangements as they function. Due to the inherent stability of local RNA structure, proteins are required to assist with RNA conformational transitions during initial folding and in exchange between functional structures. DEAD-box proteins are superfamily 2 RNA helicases that are ubiquitously involved in RNA-mediated processes. Some of these proteins use an ATP-dependent cycle of conformational changes to disrupt RNA structure nonprocessively, accelerating structural transitions of RNAs and RNPs in a manner that bears a strong resemblance to the activities of certain groups of protein chaperones. This review summarizes recent work using model substrates and tractable self-splicing intron RNAs, which has given new insights into how DEAD-box proteins promote RNA folding steps and conformational transitions, and it summarizes recent progress in identifying sites and mechanisms of DEAD-box protein activity within more complex cellular targets.

Coarse-grained model for simulation of RNA three-dimensional structures

The accurate prediction of an RNA's three-dimensional structure from its "primary structure" will have a tremendous influence on the experimental design and its interpretation and ultimately our understanding of the many functions of RNA. This paper presents a general coarse-grained (CG) potential for modeling RNA 3-D structures. Each nucleotide is represented by five pseudo atoms, two for the backbone (one for the phosphate and another for the sugar) and three for the base to represent base-stacking interactions. The CG potential has been parametrized from statistical analysis of 688 RNA experimental structures. Molecular dynamic simulations of 15 RNA molecules with the length of 12-27 nucleotides have been performed using the CG potential, with performance comparable to that from all-atom simulations. For ~75% of systems tested, simulated annealing led to native-like structures at least once out of multiple repeated runs. Furthermore, with weak distance restraints based on the knowledge of three to five canonical Watson-Crick pairs, all 15 RNAs tested are successfully folded to within 6.5 Å of native structures using the CG potential and simulated annealing. The results reveal that with a limited secondary structure model the current CG potential can reliably predict the 3-D structures for small RNA molecules. We also explored an all-atom force field to construct atomic structures from the CG simulations.

Generation and development of RNA ligase ribozymes with modular architecture through "design and selection"

In vitro selection with long random RNA libraries has been used as a powerful method to generate novel functional RNAs, although it often requires laborious structural analysis of isolated RNA molecules. Rational RNA design is an attractive alternative to avoid this laborious step, but rational design of catalytic modules is still a challenging task. A hybrid strategy of in vitro selection and rational design has been proposed. With this strategy termed “design and selection,” new ribozymes can be generated through installation of catalytic modules onto RNA scaffolds with defined 3D structures. This approach, the concept of which was inspired by the modular architecture of naturally occurring ribozymes, allows prediction of the overall architectures of the resulting ribozymes, and the structural modularity of the resulting ribozymes allows modification of their structures and functions. In this review, we summarize the design, generation, properties, and engineering of four classes of ligase ribozyme generated by design and selection.

Computational approaches to 3D modeling of RNA

Many exciting discoveries have recently revealed the versatility of RNA and its importance in a variety of functions within the cell. Since the structural features of RNA are of major importance to their biological function, there is much interest in predicting RNA structure, either in free form or in interaction with various ligands, including proteins, metabolites and other molecules. In recent years, an increasing number of researchers have developed novel RNA algorithms for predicting RNA secondary and tertiary structures. In this review, we describe current experimental and computational advances and discuss recent ideas that are transforming the traditional view of RNA folding. To evaluate the performance of the most recent RNA 3D folding algorithms, we provide a comparative study in order to test the performance of available 3D structure prediction algorithms for an RNA data set of 43 structures of various lengths and motifs. We find that the algorithms vary widely in terms of prediction quality across different RNA lengths and topologies; most predictions have very large root mean square deviations from the experimental structure. We conclude by outlining some suggestions for future RNA folding research.

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.

Exploring RNA structure by integrative molecular modelling

RNA molecular modelling is adequate to rapidly tackle the structure of RNA molecules. With new structured RNAs constituting a central class of cellular regulators discovered every year, the need for swift and reliable modelling methods is more crucial than ever. The pragmatic method based on interactive all-atom molecular modelling relies on the observation that specific structural motifs are recurrently found in RNA sequences. Once identified by a combination of comparative sequence analysis and biochemical data, the motifs composing the secondary structure of a given RNA can be extruded in three dimensions (3D) and used as building blocks assembled manually during a bioinformatic interactive process. Comparing the models to the corresponding crystal structures has validated the method as being powerful to predict the RNA topology and architecture while being less accurate regarding the prediction of base-base interactions. These aspects as well as the necessary steps towards automation will be discussed.

RNA reprogramming and repair based on trans-splicing group I ribozymes

While many traditional gene therapy strategies attempt to deliver new copies of wild-type genes back to cells harboring the defective genes, RNA-directed strategies offer a range of novel therapeutic applications. Revision or reprogramming of mRNA is a form of gene therapy that modifies mRNA without directly changing the transcriptional regulation or the genomic gene sequence. Group I ribozymes can be engineered to act in trans by recognizing a separate RNA molecule in a sequence-specific manner, and to covalently link a new RNA sequence to this separate RNA molecule. Group I ribozymes have been shown to repair defective transcripts that cause human genetic or malignant diseases, as well as to replace transcript sequences by foreign RNA resulting in new cellular functions. This review provides an overview of current strategies using trans-splicing group I ribozymes in RNA repair and reprogramming.

Identification of antisense RNA stem-loops that inhibit RNA-protein interactions using a bacterial reporter system

Many well-characterized examples of antisense RNAs from prokaryotic systems involve hybridization of the looped regions of stem–loop RNAs, presumably due to the high thermodynamic stability of the resulting loop–loop and loop–linear interactions. In this study, the identification of RNA stem–loops that inhibit U1A protein binding to the hpII RNA through RNA–RNA interactions was attempted using a bacterial reporter system based on phage λ N-mediated antitermination. As a result, loop sequences possessing 7–8 base complementarity to the 5′ region of the boxA element important for functional antitermination complex formation, but not the U1 hpII loop, were identified. In vitro and in vivo mutational analysis strongly suggested that the selected loop sequences were binding to the boxA region, and that the structure of the antisense stem–loop was important for optimal inhibitory activity. Next, in an attempt to demonstrate the ability to inhibit the interaction between the U1A protein and the hpII RNA, the rational design of an RNA stem–loop that inhibits U1A-binding to a modified hpII was carried out. Moderate inhibitory activity was observed, showing that it is possible to design and select antisense RNA stem–loops that disrupt various types of RNA–protein interactions.

The Mrs1 splicing factor binds the bI3 group I intron at each of two tetraloop-receptor motifs

Most large ribozymes require protein cofactors in order to function efficiently. The yeast mitochondrial bI3 group I intron requires two proteins for efficient splicing, Mrs1 and the bI3 maturase. Mrs1 has evolved from DNA junction resolvases to function as an RNA cofactor for at least two group I introns; however, the RNA binding site and the mechanism by which Mrs1 facilitates splicing were unknown. Here we use high-throughput RNA structure analysis to show that Mrs1 binds a ubiquitous RNA tertiary structure motif, the GNRA tetraloop-receptor interaction, at two sites in the bI3 RNA. Mrs1 also interacts at similar tetraloop-receptor elements, as well as other structures, in the self-folding Azoarcus group I intron and in the RNase P enzyme. Thus, Mrs1 recognizes general features found in the tetraloop-receptor motif. Identification of the two Mrs1 binding sites now makes it possible to create a model of the complete six-component bI3 ribonucleoprotein. All protein cofactors bind at the periphery of the RNA such that every long-range RNA tertiary interaction is stabilized by protein binding, involving either Mrs1 or the bI3 maturase. This work emphasizes the strong evolutionary pressure to bolster RNA tertiary structure with RNA-binding interactions as seen in the ribosome, spliceosome, and other large RNA machines.

Methods of site-specific labeling of RNA with fluorescent dyes

Single molecule fluorescence techniques offer unique insights into mechanisms of conformational changes of RNA. Knowing how to make fluorescently labeled RNA molecules and understanding potential limitations of different labeling strategies is essential for successful implementation of single molecule fluorescence techniques. This chapter offers a step by step overview of the process of obtaining RNA constructs ready for single molecule measurements. Several alternative methods are described for each step, and ways of troubleshooting the most common problems, in particular, splinted RNA ligation, are suggested.

Protein roles in group I intron RNA folding: the tyrosyl-tRNA synthetase CYT-18 stabilizes the native state relative to a long-lived misfolded structure without compromising folding kinetics

The Neurospora crassa CYT-18 protein is a mitochondrial tyrosyl-tRNA synthetase that also promotes self-splicing of group I intron RNAs by stabilizing the functional structure in the conserved core. CYT-18 binds the core along the same surface as a common peripheral element, P5abc, suggesting that CYT-18 can replace P5abc functionally. In addition to stabilizing structure generally, P5abc stabilizes the native conformation of the Tetrahymena group I intron relative to a globally similar misfolded conformation that has only local differences within the core and is populated significantly at equilibrium by a ribozyme variant lacking P5abc (E(DeltaP5abc)). Here, we show that CYT-18 specifically promotes formation of the native group I intron core from this misfolded conformation. Catalytic activity assays demonstrate that CYT-18 shifts the equilibrium of E(DeltaP5abc) toward the native state by at least 35-fold, and binding assays suggest an even larger effect. Thus, similar to P5abc, CYT-18 preferentially recognizes the native core, despite the global similarity of the misfolded core and despite forming crudely similar complexes, as revealed by dimethyl sulfate footprinting. Interestingly, the effects of CYT-18 and P5abc on folding kinetics differ. Whereas P5abc inhibits refolding of the misfolded conformation by forming peripheral contacts that must break during refolding, CYT-18 does not display analogous inhibition, most likely because it relies to a greater extent on direct interactions with the core. Although CYT-18 does not encounter this RNA in vivo, our results suggest that it stabilizes its cognate group I introns relative to analogous misfolded intermediates. By specifically recognizing native structural features, CYT-18 may also interact with earlier folding intermediates to avoid RNA misfolding or to trap native contacts as they form. More generally, our results highlight the ability of a protein cofactor to stabilize a functional RNA structure specifically without incurring associated costs in RNA folding kinetics.

Metal ion dependence of cooperative collapse transitions in RNA

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

Knowledge-based instantiation of full atomic detail into coarse-grain RNA 3D structural models

Motivation: The recent development of methods for modeling RNA 3D structures using coarse-grain approaches creates a need to bridge low- and high-resolution modeling methods. Although they contain topological information, coarse-grain models lack atomic detail, which limits their utility for some applications.

Results: We have developed a method for adding full atomic detail to coarse-grain models of RNA 3D structures. Our method [Coarse to Atomic (C2A)] uses geometries observed in known RNA crystal structures. Our method rebuilds full atomic detail from ideal coarse-grain backbones taken from crystal structures to within 1.87–3.31 Å RMSD of the full atomic crystal structure. When starting from coarse-grain models generated by the modeling tool NAST, our method builds full atomic structures that are within 1.00 Å RMSD of the starting structure. The resulting full atomic structures can be used as starting points for higher resolution modeling, thus bridging high- and low-resolution approaches to modeling RNA 3D structure.

Availability: Code for the C2A method, as well as the examples discussed in this article, are freely available at www.simtk.org/home/c2a.

Contact:russ.altman@stanford.edu

Probing the dynamics of the P1 helix within the Tetrahymena group I intron

RNA conformational transformations are integral to RNA's biological functions. Further, structured RNA molecules exist as a series of dynamic intermediates in the course of folding or complexation with proteins. Thus, an understanding of RNA folding and function will require deep and incisive understanding of its dynamic behavior. However, existing tools to investigate RNA dynamics are limited. Here, we introduce a powerful fluorescence polarization anisotropy approach that utilizes a rare base analogue that retains substantial fluorescence when incorporated into helices. We show that 6-methylisoxanthopterin (6-MI) can be used to follow the nanosecond dynamics of individual helices. We then use 6-MI to probe the dynamics of an individual helix, referred to as P1, within the 400nt Tetrahymena group I ribozyme. Comparisons of the dynamics of the P1 helix in wild type and mutant ribozymes and in model constructs reveal a highly immobilized docked state of the P1 helix, as expected, and a relatively mobile "open complex" or undocked state. This latter result rules out a model in which slow docking of the P1 helix into its cognate tertiary interactions arises from a stable alternatively docked conformer. The results are consistent with a model in which stacking and tertiary interactions of the A(3) tether connecting the P1 helix to the body of the ribozyme limit P1 mobility and slow its docking, and this model is supported by cross-linking results. The ability to isolate the nanosecond motions of individual helices within complex RNAs and RNA/protein complexes will be valuable in distinguishing between functional models and in discerning the fundamental behavior of important biological species.

Tertiary motifs revealed in analyses of higher-order RNA junctions

RNA junctions are secondary-structure elements formed when three or more helices come together. They are present in diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze currently solved 3D RNA junctions in terms of base-pair interactions and 3D configurations. First, we study base-pair interaction diagrams for solved RNA junctions with 5 to 10 helices and discuss common features. Second, we compare these higher-order junctions to those containing 3 or 4 helices and identify global motif patterns such as coaxial stacking and parallel and perpendicular helical configurations. These analyses show that higher-order junctions organize their helical components in parallel and helical configurations similar to lower-order junctions. Their sub-junctions also resemble local helical configurations found in three- and four-way junctions and are stabilized by similar long-range interaction preferences such as A-minor interactions. Furthermore, loop regions within junctions are high in adenine but low in cytosine, and in agreement with previous studies, we suggest that coaxial stacking between helices likely forms when the common single-stranded loop is small in size; however, other factors such as stacking interactions involving noncanonical base pairs and proteins can greatly determine or disrupt coaxial stacking. Finally, we introduce the ribo-base interactions: when combined with the along-groove packing motif, these ribo-base interactions form novel motifs involved in perpendicular helix-helix interactions. Overall, these analyses suggest recurrent tertiary motifs that stabilize junction architecture, pack helices, and help form helical configurations that occur as sub-elements of larger junction networks. The frequent occurrence of similar helical motifs suggest nature's finite and perhaps limited repertoire of RNA helical conformation preferences. More generally, studies of RNA junctions and tertiary building blocks can ultimately help in the difficult task of RNA 3D structure prediction.

Dynamic extended folding: modeling the RNA secondary structures during co-transcriptional folding

For RNA secondary structure prediction, it is an important issue that how to deal with co-transcriptional folding during the RNA synthesis in the cell. On one hand, co-transcriptional folding, leads to the correct final structure of the whole RNA molecule. On the other hand, it may form the recognition sites for the progress of the transcription. Considering the hurdles in the experimental determination of RNA folding structures, we proposed a so-called "dynamic extended folding simulation" approach. We used two human pre-mRNA samples, the first functional alpha-gene HBZ and the fifth beta-gene HBB, to "display" the co-transcriptional folding images in detail. The modeling process starts from the prediction of a 30-nucleotide (nt) sequence, then in each update 30 nts was extended, say, 1-30, 1-60, 1-90, 1-120,..., 1-1651 nts (for HBB, 1-1606 nts). We selected the RNAstructure program to predict the folding secondary structures of all the segments. We defined "hairpin" as the unit of the secondary structure and analyzed the states of such unit during the sequential dynamic extended folding processes. We found that some hairpins are "conserved", i.e., after its appearance, it always is there in the followed foldings. Some hairpins present partially in the folding segments, and some hairpins appear for only once or twice. This phenomenon vividly depicts the generation and adjusting of the temporal structural units during the co-transcriptional folding process. It is these "hairpins" that support the thermodynamically stable structure at the end of the RNA synthesis. They may also play a role in RNA splicing process and even in the folding structure of the synthesized protein.

Coarse-grained modeling of large RNA molecules with knowledge-based potentials and structural filters

Understanding the function of complex RNA molecules depends critically on understanding their structure. However, creating three-dimensional (3D) structural models of RNA remains a significant challenge. We present a protocol (the nucleic acid simulation tool [NAST]) for RNA modeling that uses an RNA-specific knowledge-based potential in a coarse-grained molecular dynamics engine to generate plausible 3D structures. We demonstrate NAST's capabilities by using only secondary structure and tertiary contact predictions to generate, cluster, and rank structures. Representative structures in the best ranking clusters averaged 8.0 +/- 0.3 A and 16.3 +/- 1.0 A RMSD for the yeast phenylalanine tRNA and the P4-P6 domain of the Tetrahymena thermophila group I intron, respectively. The coarse-grained resolution allows us to model large molecules such as the 158-residue P4-P6 or the 388-residue T. thermophila group I intron. One advantage of NAST is the ability to rank clusters of structurally similar decoys based on their compatibility with experimental data. We successfully used ideal small-angle X-ray scattering data and both ideal and experimental solvent accessibility data to select the best cluster of structures for both tRNA and P4-P6. Finally, we used NAST to build in missing loops in the crystal structures of the Azoarcus and Twort ribozymes, and to incorporate crystallographic data into the Michel-Westhof model of the T. thermophila group I intron, creating an integrated model of the entire molecule. Our software package is freely available at https://simtk.org/home/nast.

Toward predicting self-splicing and protein-facilitated splicing of group I introns

In the current era of massive discoveries of noncoding RNAs within genomes, being able to infer a function from a nucleotide sequence is of paramount interest. Although studies of individual group I introns have identified self-splicing and nonself-splicing examples, there is no overall understanding of the prevalence of self-splicing or the factors that determine it among the >2300 group I introns sequenced to date. Here, the self-splicing activities of 12 group I introns from various organisms were assayed under six reaction conditions that had been shown previously to promote RNA catalysis for different RNAs. Besides revealing that assessing self-splicing under only one condition can be misleading, this survey emphasizes that in vitro self-splicing efficiency is correlated with the GC content of the intron (>35% GC was generally conductive to self-splicing), and with the ability of the introns to form particular tertiary interactions. Addition of the Neurospora crassa CYT-18 protein activated splicing of two nonself-splicing introns, but inhibited the second step of self-splicing for two others. Together, correlations between sequence, predicted structure and splicing begin to establish rules that should facilitate our ability to predict the self-splicing activity of any group I intron from its sequence.

Splicing and evolution of an unusually small group I intron

Introns are common in the rRNA gene loci of fungal genomes, but biochemical studies to investigate splicing are rare. Here, self-splicing of a very small (67 nucleotide) group I intron is demonstrated. The PaSSU intron (located within the rRNA small subunit gene of Phialophora americana) splices in vitro under group I intron conditions. Most group I ribozymes contain pairing regions P1-P10, with a conserved G.U pair at the 5' splice site, and a G at the 3' intron border. The PaSSU intron contains only P1, P7, and P10. While it contains the G.U pair at the 5' splice, a U is found at the 3' end of the intron instead of a G. Phylogenetic analysis places it within subgroup IC1, whose members are found in the nuclear rRNA genes of fungi. The structural elements are similar to those in the centermost regions of other group I introns. Its size can be explained by a single large deletion that removed P2 through much of P9. Part of the original P9 region has assumed the function of P7. Its small size and genealogy makes it an excellent model to study RNA catalysis and evolution.

Time-resolved footprinting for the study of the structural dynamics of DNA–protein interactions

Transcription is often regulated at the level of initiation by the presence of transcription factors or nucleoid proteins or by changing concentrations of metabolites. These can influence the kinetic properties and/or structures of the intermediate RNA polymerase–DNA complexes in the pathway. Time-resolved footprinting techniques combine the high temporal resolution of a stopped-flow apparatus with the specific structural information obtained by the probing agent. Combined with a careful quantitative analysis of the evolution of the signals, this approach allows for the identification and kinetic and structural characterization of the intermediates in the pathway of DNA sequence recognition by a protein, such as a transcription factor or RNA polymerase. The combination of different probing agents is especially powerful in revealing different aspects of the conformational changes taking place at the protein–DNA interface. For example, hydroxyl radical footprinting, owing to their small size, provides a map of the solvent-accessible surface of the DNA backbone at a single nucleotide resolution; modification of the bases using potassium permanganate can reveal the accessibility of the bases when the double helix is distorted or melted; cross-linking experiments report on the formation of specific amino acid–DNA contacts, and DNase I footprinting results in a strong signal-to-noise ratio from DNA protection at the binding site and hypersensitivity at curved or kinked DNA sites. Recent developments in protein footprinting allow for the direct characterization of conformational changes of the proteins in the complex.

A three-dimensional model of a group II intron RNA and its interaction with the intron-encoded reverse transcriptase

Group II introns are self-splicing ribozymes believed to be the ancestors of spliceosomal introns. Many group II introns encode reverse transcriptases that promote both RNA splicing and intron mobility to new genomic sites. Here we used a circular permutation and crosslinking method to establish 16 intramolecular distance relationships within the mobile Lactococcus lactis Ll.LtrB-DeltaORF intron. Using these new constraints together with 13 established tertiary interactions and eight published crosslinks, we modeled a complete three-dimensional structure of the intron. We also used the circular permutation strategy to map RNA-protein interaction sites through fluorescence quenching and crosslinking assays. Our model provides a comprehensive structural framework for understanding the function of group II ribozymes, their natural structural variations, and the mechanisms by which the intron-encoded protein promotes RNA splicing and intron mobility. The model also suggests an arrangement of active site elements that may be conserved in the spliceosome.

Unwinding RNA's secrets: advances in the biology, physics, and modeling of complex RNAs

The rapid development of our understanding of the diverse biological roles fulfilled by non-coding RNA has motivated interest in the basic macromolecular behavior, structure, and function of RNA. We focus on two areas in the behavior of complex RNAs. First, we present advances in the understanding of how RNA folding is accomplished in vivo by presenting a mechanism for the action of DEAD-box proteins. Members of this family are intimately associated with almost all cellular processes involving RNA, mediating RNA structural rearrangements and chaperoning their folding. Next, we focus on advances in understanding, and characterizing the basic biophysical forces that govern the folding of complex RNAs. Ultimately we expect that a confluence and synergy between these approaches will lead to profound understanding of RNA and its biology.