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.

An in vivo selection method to optimize trans-splicing ribozymes

Group I intron ribozymes can repair mutated mRNAs by replacing the 3'-terminal portion of the mRNA with their own 3'-exon. This trans-splicing reaction has the potential to treat genetic disorders and to selectively kill cancer cells or virus-infected cells. However, these ribozymes have not yet been used in therapy, partially due to a low in vivo trans-splicing efficiency. Previous strategies to improve the trans-splicing efficiencies focused on designing and testing individual ribozyme constructs. Here we describe a method that selects the most efficient ribozymes from millions of ribozyme variants. This method uses an in vivo rescue assay where the mRNA of an inactivated antibiotic resistance gene is repaired by trans-splicing group I intron ribozymes. Bacterial cells that express efficient trans-splicing ribozymes are able to grow on medium containing the antibiotic chloramphenicol. We randomized a 5'-terminal sequence of the Tetrahymena thermophila group I intron and screened a library with 9 × 10⁶ ribozyme variants for the best trans-splicing activity. The resulting ribozymes showed increased trans-splicing efficiency and help the design of efficient trans-splicing ribozymes for different sequence contexts. This in vivo selection method can now be used to optimize any sequence in trans-splicing ribozymes.

Plasticity of the RNA kink turn structural motif

The kink turn (K-turn) is an RNA structural motif found in many biologically significant RNAs. While most examples of the K-turn have a similar fold, the crystal structure of the Azoarcus group I intron revealed a novel RNA conformation, a reverse kink turn bent in the direction opposite that of a consensus K-turn. The reverse K-turn is bent toward the major grooves rather than the minor grooves of the flanking helices, yet the sequence differs from the K-turn consensus by only a single nucleotide. Here we demonstrate that the reverse bend direction is not solely defined by internal sequence elements, but is instead affected by structural elements external to the K-turn. It bends toward the major groove under the direction of a tetraloop-tetraloop receptor. The ability of one sequence to form two distinct structures demonstrates the inherent plasticity of the K-turn sequence. Such plasticity suggests that the K-turn is not a primary element in RNA folding, but instead is shaped by other structural elements within the RNA or ribonucleoprotein assembly.

A conformational switch in the DiGIR1 ribozyme involved in release and folding of the downstream I-DirI mRNA

DiGIR1 is a group I-like cleavage ribozyme found as a structural domain within a nuclear twin-ribozyme group I intron. DiGIR1 catalyzes cleavage by branching at an Internal Processing Site (IPS) leading to formation of a lariat cap at the 5'-end of the 3'-cleavage product. The 3'-cleavage product is subsequently processed into an mRNA encoding a homing endonuclease. By analysis of combinations of 5'- and 3'-deletions, we identify a hairpin in the 5'-UTR of the mRNA (HEG P1) that is formed by conformational switching following cleavage. The formation of HEG P1 inhibits the reversal of the branching reaction, thus giving it directionality. Furthermore, the release of the mRNA is a consequence of branching rather than hydrolytic cleavage. A model is put forward that explains the release of the I-DirI mRNA with a lariat cap and a structured 5'-UTR as a direct consequence of the DiGIR1 branching reaction. The role of HEG P1 in GIR1 branching is reminiscent of that of hairpin P-1 in splicing of the Tetrahymena rRNA group I intron and illustrates a general principle in RNA-directed RNA processing.

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.

RNA kink turns to the left and to the right

A helix-loop-helix within the group I intron has most of the canonical sequence elements of a kink turn (K-turn), yet it bends in the opposite direction. The reverse K-turn kinks toward the major rather than the minor grooves of the flanking helices. This suggests that there are two distinct subclasses of tertiary structures that a K-turn secondary structure can adopt. The final structure may be specified by external factors, such as protein binding or the tertiary structural context, rather than the intrinsic conformation of the RNA.

A base triple in the Tetrahymena group I core affects the reaction equilibrium via a threshold effect

Previous work on group I introns has suggested that a central base triple might be more important for the first rather than the second step of self-splicing, leading to a model in which the base triple undergoes a conformational change during self-splicing. Here, we use the well-characterized L-21 ScaI ribozyme derived from the Tetrahymena group I intron to probe the effects of base-triple disruption on individual reaction steps. Consistent with previous results, reaction of a ternary complex mimicking the first chemical step in self-splicing is slowed by mutations in this base triple, whereas reaction of a ternary complex mimicking the second step of self-splicing is not. Paradoxically, mechanistic dissection of the base-triple disruption mutants indicates that active site binding is weakened uniformly for the 5'-splice site and the 5'-exon analog, mimics for the species bound in the first and second step of self-splicing. Nevertheless, the 5'-exon analog remains bound at the active site, whereas the 5'-splice site analog does not. This differential effect arises despite the uniform destabilization, because the wild-type ribozyme binds the 5'-exon analog more strongly in the active site than in the 5'-splice site analog. Thus, binding into the active site constitutes an additional barrier to reaction of the 5'-splice site analog, but not the 5'-exon analog, resulting in a reduced reaction rate constant for the first step analog, but not the second step analog. This threshold model explains the self-splicing observations without the need to invoke a conformational change involving the base triple, and underscores the importance of quantitative dissection for the interpretation of effects from mutations.

The ability to form full-length intron RNA circles is a general property of nuclear group I introns

In addition to splicing, group I intron RNA is capable of an alternative two-step processing pathway that results in the formation of full-length intron circular RNA. The circularization pathway is initiated by hydrolytic cleavage at the 3' splice site and followed by a transesterification reaction in which the intron terminal guanosine attacks the 5' splice site presented in a structure analogous to that of the first step of splicing. The products of the reactions are full-length circular intron and unligated exons. For this reason, the circularization reaction is to the benefit of the intron at the expense of the host. The circularization pathway has distinct structural requirements that differ from those of splicing and appears to be specifically suppressed in vivo. The ability to form full-length circles is found in all types of nuclear group I introns, including those from the Tetrahymena ribosomal DNA. The biological function of the full-length circles is not known, but the fact that the circles contain the entire genetic information of the intron suggests a role in intron mobility.

Effect of transcription on folding of the Tetrahymena ribozyme

Sequential formation of RNA interactions during transcription can bias the folding pathway and ultimately determine the functional state of a transcript. The kinetics of cotranscriptional folding of the Tetrahymena L-21 ribozyme was compared with refolding of full-length transcripts under the same conditions. Sequential folding after transcription by phage T7 or Escherichia coli polymerase is only twice as fast as refolding, and the yield of native RNA is the same. By contrast, a greater fraction of circularly permuted variants folded correctly at early times during transcription than during refolding. Hybridization of complementary oligonucleotides suggests that cotranscriptional folding enables a permuted RNA beginning at G303 to escape non-native interactions in P3 and P9. We propose that base pairing of upstream sequences during transcription elongation favors branched secondary structures that increase the probability of forming the native ribozyme structure.

RNA chaperone StpA loosens interactions of the tertiary structure in the td group I intron in vivo

Efficient splicing of the td group I intron in vivo is dependent on the ribosome. In the absence of translation, the pre-mRNA is trapped in nonnative-splicing-incompetent conformations. Alternatively, folding of the pre-mRNA can be promoted by the RNA chaperone StpA or by the group I intron-specific splicing factor Cyt-18. To understand the mechanism of action of RNA chaperones, we probed the impact of StpA on the structure of the td intron in vivo. Our data suggest that StpA loosens tertiary interactions. The most prominent structural change was the opening of the base triples, which are involved in the correct orientation of the two major intron core domains. In line with the destabilizing activity of StpA, splicing of mutant introns with a reduced structural stability is sensitive to StpA. In contrast, Cyt-18 strengthens tertiary contacts, thereby rescuing splicing of structurally compromised td mutants in vivo. Our data provide direct evidence for protein-induced conformational changes within catalytic RNA in vivo. Whereas StpA resolves tertiary contacts enabling the RNA to refold, Cyt-18 contributes to the overall compactness of the td intron in vivo.