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.     

Dual roles for the Mss116 cofactor during splicing of the ai5γ group II intron

The autocatalytic group II intron ai5γ from Saccharomyces cerevisiae self-splices under high-salt conditions in vitro, but requires the assistance of the DEAD-box protein Mss116 in vivo and under near-physiological conditions in vitro. Here, we show that Mss116 influences the folding mechanism in several ways. By comparing intron precursor RNAs with long (∼300 nt) and short (∼20 nt) exons, we observe that long exon sequences are a major obstacle for self-splicing in vitro. Kinetic analysis indicates that Mss116 not only mitigates the inhibitory effects of long exons, but also assists folding of the intron core. Moreover, a mutation in conserved Motif III that impairs unwinding activity (SAT → AAA) only affects the construct with long exons, suggesting helicase unwinding during exon unfolding, but not in intron folding. Strong parallels between Mss116 and the related protein Cyt-19 from Neurospora crassa suggest that these proteins form a subclass of DEAD-box proteins that possess a versatile repertoire of diverse activities for resolving the folding problems of large RNAs.     

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^ΔP5abc). 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^ΔP5abc 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.     

Interlocked RNA circle formation by a self-splicing yeast mitochondrial group I intron

RNA containing the aI3 group I intron of the yeast mitochondrial gene encoding cytochrome oxidase subunit I shows self-splicing in vitro. The excised intron, comprising 1514 nucleotides, is partially split into an upstream portion, containing the intronic reading frame, and a downstream portion, containing the typical group I conserved sequence elements. Full-length intron RNA and intron parts occur in linear and circular form. In the transesterification reactions leading to circle formation, only the guanosine nucleotide added during splicing is removed. Reincubation of isolated, complete circular intron RNA under self-splicing conditions leads to formation of free subintronic RNA circles. Under similar conditions, purified linear intron RNA gives rise to a number of circular and linear products, one of which consists of interlocked subintronic RNA circles. These observations suggest that the intron RNA possesses a dynamic structure in which subtle alterations in folding result in the formation of RNA products with different topology.     

Characterization of products derived from self-splicing of intron aI5 alpha which is located in the mitochondrial COX I gene of Saccharomyces cerevisiae.

We have characterized the in vitro self-splicing of intron aI5 alpha containing precursor RNA from the yeast mitochondrial gene coding for cytochrome oxidase subunit I. This intron follows the rules for group I self-splicing introns and all the characteristic products have been identified. In addition we have detected abnormal RNA products with features that indicate that the self-splicing behaviour of this intron is more complex. Two intron circles are formed by use of a major and minor intron-internal site for circle closure. A cryptic 5'-splice site located in the 3' exon results in guanosine nucleotide mediated opening at a position 30 nt downstream of the normal 3' splice site. The reactions can all be explained on the basis of the "splice guide" model proposed by Davies et al (1982 Nature 300 719-724). Although the sequence motifs at cyclization and splice sites occur more often in this intron, only some of them are allowed to interact with the internal guide sequence, suggesting that both primary structure and spatial folding of the RNA are involved in formation of productive reaction sites.     

Antibiotic-induced Oligomerisation of Group I Intron RNA

Antibiotics act as inhibitors of various biological processes. Here we demonstrate that some tuberactinomycins, hitherto known as inhibitors of prokaryotic protein synthesis and of group I intron self-splicing, have a modulatory effect on group I intron RNAs. The linear intron, which is excised during the self-splicing process, is still an active molecule capable of performing an intramolecular transesterification resulting in a circular molecule. However, in the presence of sub-inhibitory concentrations of tuberactinomycins, the intron reacts intermolecularly leading to the formation of linear head-to-tail intron-oligomers. The antibiotic stimulates the intron to reactin transinstead ofin cis. The phage T4-derivedtdintron uses the same sites for oligomerisation as for circularisation. Gel-retardation experiments demonstrate that the intron RNA forms non-covalent complexes in the presence of the antibiotic. It might be envisaged that the role of these peptide antibiotics is to bridge RNA molecules mediating RNA – RNA interactions and thus enabling their reaction. The tuberactinomycins are further able to induce the interaction of heterologous introns. The ligation of the T4 phage-derivedtdintron with theTetrahymenarRNA intron is very efficient, resulting in molecules composed of two introns derived from different species. Thetdintron attacks theTetraymenaintron at various sites, which are located within double-stranded regions. These observations suggest that small molecules like these basic peptide antibiotics could have mediated RNA–RNA interactions in a pre-protein era.     

Self-splicing of a group II intron in yeast mitochondria: dependence on 5'' exon sequences.

It has been previously suggested that self-splicing of group II introns starts with a nucleophilic attack of the 2' OH group from the branchpoint adenosine on the 5' splice junction. To investigate the sequences governing the specificity of this attack, a series of Bal31 nuclease deletion mutants was constructed in which progressively larger amounts of 5' exon have been removed starting from its 5' end. The ability of mutant RNAs to carry out self-splicing in vitro was studied. Involvement of 5' exon sequences in self-splicing activity is indicated by the fact that a mutant in which as many as 18 nucleotides of 5' exon remain is seriously disturbed in splicing, while larger deletions eliminate splicing entirely. Mutants containing a truncated 5' exon form aberrant RNAs. One of these is a 425-nucleotide RNA containing the 5' exon as well as sequences of the 5' part of the intron. Its 3' end maps at position 374 of the 887-nucleotide intron. The other is a less abundant lariat RNA probably originating from the remainder of the intron linked to the 3' exon. We interpret this large dependence of reactivity of the intron on 5' exon and adjoining intron sequences as evidence for base-pairing interactions between the exon and parts of the intron, leading to an RNA folding necessary for splicing. Possible folding models are discussed.     

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

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

Protein-Facilitated Folding of Group II Intron Ribozymes

Multiple studies hypothesize that DEAD-box proteins facilitate folding of the ai5γ group II intron. However, these conclusions are generally inferred from splicing kinetics, and not from direct monitoring of DEAD-box protein-facilitated folding of the intron. Using native gel electrophoresis and dimethyl sulfate structural probing, we monitored Mss-116-facilitated folding of ai5γ intron ribozymes and a catalytically active self-splicing RNA containing full-length intron and short exons. We found that the protein directly stimulates folding of these RNAs by accelerating formation of the compact near-native state. This process occurs in an ATP-independent manner, although ATP is required for the protein turnover. As Mss 116 binds RNA nonspecifically, most binding events do not result in the formation of the compact state, and ATP is required for the protein to dissociate from such nonproductive complexes and rebind the unfolded RNA. Results obtained from experiments at different concentrations of magnesium ions suggest that Mss 116 stimulates folding of ai5γ ribozymes by promoting the formation of unstable folding intermediates, which is then followed by a cascade of folding events resulting in the formation of the compact near-native state. Dimethyl sulfate probing results suggest that the compact state formed in the presence of the protein is identical to the near-native state formed more slowly in its absence. Our results also indicate that Mss 116 does not stabilize the native state of the ribozyme, but that such stabilization results from binding of attached exons.     

Fluorescence of covalently attached pyrene as a general RNA folding probe

Fluorescence techniques are commonly and powerfully applied to monitor biomolecular folding. In a limited fashion, the fluorescence emission intensity of covalently attached pyrene has been used as a reporter of RNA conformational changes. Here, we pursue two goals: we examine the relationship between tether identity and fluorescence response, and we determine the general utility of pyrene fluorescence to monitor RNA folding. The P4–P6 domain of the Tetrahymena group I intron RNA was systematically modified at multiple nucleotide positions with pyrene derivatives that provide a range of tether lengths and compositions between the RNA and chromophore. Certain tethers typically lead to a superior fluorescence signal upon RNA folding, as demonstrated by equilibrium titrations with Mg²⁺. In addition, useful fluorescence responses were obtained with pyrene placed at several nucleotide positions dispersed throughout P4–P6. This suggests that monitoring of tertiary folding by fluorescence of covalently attached pyrene will be generally applicable to structured RNA molecules.     

Influence of RNA structural stability on the RNA chaperone activity of the Escherichia coli protein StpA

Proteins with RNA chaperone activity are able to promote folding of RNA molecules by loosening their structure. This RNA unfolding activity is beneficial when resolving misfolded RNA conformations, but could be detrimental to RNAs with low thermodynamic stability. In order to test this idea, we constructed various RNAs with different structural stabilities derived from the thymidylate synthase (td) group I intron and measured the effect of StpA, an Escherichia coli protein with RNA chaperone activity, on their splicing activity in vivo and in vitro. While StpA promotes splicing of the wild-type td intron and of mutants with wild-type-like stability, splicing of mutants with a lower structural stability is reduced in the presence of StpA. In contrast, splicing of an intron mutant, which is not destabilized but which displays a reduced population of correctly folded RNAs, is promoted by StpA. The sensitivity of an RNA towards StpA correlates with its structural stability. By lowering the temperature to 25°C, a temperature at which the structure of these mutants becomes more stable, StpA is again able to stimulate splicing. These observations clearly suggest that the structural stability of an RNA determines whether the RNA chaperone activity of StpA is beneficial to folding.     

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

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

A molecular beacon assay for monitoring RNA splicing

Small molecule targeting of self-splicing RNAs like group I and II introns has been limited in part by the lack of a universal high-throughput screening platform for studies of splicing inhibition and kinetics. Here, we present the development of a molecular beacon assay for monitoring the accumulation of spliced exons during RNA splicing reactions. In this case, we applied it to the autocatalyzed reaction of the H.c.LSU group II intron found in the mitochondria of the pathogenic dimorphic fungus Histoplasma capsulatum. We find that a molecular beacon with the loop length of 18 nucleotides selectively recognizes ligated exons formed during self-splicing and exhibits high fluorescent signal upon binding of its target. We demonstrate that the fluorescent assay using molecular beacons can be successfully applied to kinetic characterization of the splicing reaction and determination of inhibition constants for small molecules. The results presented herein offer support for a molecular beacon approach to identifying small molecule inhibitors of intron splicing.     

Two major tertiary folding transitions of the Tetrahymena catalytic RNA.

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

The brace for a growing scaffold: mss116 protein promotes RNA folding by stabilizing an early assembly intermediate

The ai5γ group II intron requires a protein cofactor to facilitate native folding in the cell. Yeast protein Mss116 greatly accelerates intron folding under near-physiological conditions both in vivo and in vitro. Although the effect of Mss116 on the kinetics of ai5γ ribozyme folding and catalysis has been extensively studied, the precise structural role and interaction sites of Mss116 have been elusive. Using Nucleotide Analog Interference Mapping to study the folding of splicing precursor constructs, we have identified specific intron functional groups that participate in Mss116-facilitated folding and we have determined their role in the folding mechanism. The data indicate that Mss116 stabilizes an early, obligate folding intermediate within intron domain 1, thereby laying the foundation for productive folding to the native state. In addition, the data reveal an important role for the IBS2 exon sequence and for the terminus of domain 6, during the folding of self-splicing group IIB intron constructs.     

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

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