Kinetic and secondary structure analysis of Naegleria andersoni GIR1, a group I ribozyme whose putative biological function is site-specific hydrolysis

Molecular characterization of a new member of the lariat capping twin-ribozyme introns

BACKGROUND

Twin-ribozyme introns represent a complex class of mobile group I introns that harbour a lariat capping (LC) ribozyme and a homing endonuclease gene embedded in a conventional self-splicing group I ribozyme (GIR2). Twin-ribozyme introns have so far been confined to nucleolar DNA in Naegleria amoeboflagellates and the myxomycete Didymium iridis.

RESULTS

We characterize structural organization, catalytic properties and molecular evolution of a new twin-ribozyme intron in Allovahlkampfia (Heterolobosea). The intron contains two ribozyme domains with different functions in ribosomal RNA splicing and homing endonuclease mRNA maturation. We found Allovahlkampfia GIR2 to be a typical group IC1 splicing ribozyme responsible for addition of the exogenous guanosine cofactor (exoG), exon ligation and circularization of intron RNA. The Allovahlkampfia LC ribozyme, by contrast, represents an efficient self-cleaving ribozyme that generates a small 2',5' lariat cap at the 5' end of the homing endonuclease mRNA, and thus contributes to intron mobility.

CONCLUSIONS

The discovery of a twin-ribozyme intron in a member of Heterolobosea expands the distribution pattern of LC ribozymes. We identify a putative regulatory RNA element (AP2.1) in the Allovahlkampfia LC ribozyme that involves homing endonuclease mRNA coding sequences as an important structural component.

The Naegleria genome: a free-living microbial eukaryote lends unique insights into core eukaryotic cell biology

Naegleria gruberi, a free-living protist, has long been treasured as a model for basal body and flagellar assembly due to its ability to differentiate from crawling amoebae into swimming flagellates. The full genome sequence of Naegleria gruberi has recently been used to estimate gene families ancestral to all eukaryotes and to identify novel aspects of Naegleria biology, including likely facultative anaerobic metabolism, extensive signaling cascades, and evidence for sexuality. Distinctive features of the Naegleria genome and nuclear biology provide unique perspectives for comparative cell biology, including cell division, RNA processing and nucleolar assembly. We highlight here exciting new and novel aspects of Naegleria biology identified through genomic analysis.

Experimental analyses of the chemical dynamics of ribozyme catalysis

Most ribozymes in Nature catalyze alcoholysis or hydrolysis of RNA phosphodiester bonds. Studies of the corresponding non-enzymatic reactions reveal a complex mechanistic landscape allowing for a variety of transition states and both concerted and stepwise mechanisms. High-resolution structures, incisive biochemical studies and computer simulations are providing glimpses into how ribozyme catalyzed reactions traverse this landscape. However, direct experimental tests of mechanistic detail at the chemical level are not easily achieved. Kinetic isotope effects (KIEs) probe directly the differences in the vibrational 'environment' of the atoms undergoing chemical transformation on going from the ground state to the transition state. Thus, KIEs can in principle provide direct information about transition state bonding and so may be instrumental in evaluating possible transition states for ribozyme catalyzed reactions. Understanding charge distribution in the transition state may help resolve how rate acceleration is accomplished and perhaps the similarities and differences in how RNA and protein active sites operate. Several barriers to successful application of KIE analysis to ribozymes have recently been overcome, and new chemical details are beginning to emerge.

Molecular modelling of the GIR1 branching ribozyme gives new insight into evolution of structurally related ribozymes

Twin-ribozyme introns contain a branching ribozyme (GIR1) followed by a homing endonuclease (HE) encoding sequence embedded in a peripheral domain of a group I splicing ribozyme (GIR2). GIR1 catalyses the formation of a lariat with 3 nt in the loop, which caps the HE mRNA. GIR1 is structurally related to group I ribozymes raising the question about how two closely related ribozymes can carry out very different reactions. Modelling of GIR1 based on new biochemical and mutational data shows an extended substrate domain containing a GoU pair distinct from the nucleophilic residue that dock onto a catalytic core showing a different topology from that of group I ribozymes. The differences include a core J8/7 region that has been reduced and is complemented by residues from the pre-lariat fold. These findings provide the basis for an evolutionary mechanism that accounts for the change from group I splicing ribozyme to the branching GIR1 architecture. Such an evolutionary mechanism can be applied to other large RNAs such as the ribonuclease P.

Characterization of the mechanical unfolding of RNA pseudoknots

The pseudoknot is an important RNA structural element that provides an excellent model system for studying the contributions of tertiary interactions to RNA stability and to folding kinetics. RNA pseudoknots are also of interest because of their key role in the control of ribosomal frameshifting by viral RNAs. Their mechanical properties are directly relevant to their unfolding by ribosomes during translation. We have used optical tweezers to study the kinetics and thermodynamics of mechanical unfolding and refolding of single RNA molecules. Here we describe the unfolding of the frameshifting pseudoknot from infectious bronchitis virus (IBV), three constituent hairpins, and three mutants of the IBV pseudoknot. All four pseudoknots cause −1 programmed ribosomal frameshifting. We have measured the free energies and rates of mechanical unfolding and refolding of the four frameshifting pseudoknots. Our results show that the IBV pseudoknot requires a higher force than its corresponding hairpins to unfold. Furthermore, its rate of unfolding changes little with increasing force, in contrast with the rate of hairpin unfolding. The presence of Mg²⁺ significantly increases the kinetic barriers to unfolding the IBV pseudoknot, but has only a minor effect on the hairpin unfolding. The greater mechanical stability of pseudoknots compared to hairpins, and their kinetic insensitivity to force supports the hypothesis that −1 frameshifting depends on the difficulty of unfolding the mRNA.

A new RNA branching activity: the GIR1 ribozyme

The formation of lariat intermediates during the first step of splicing of group II introns and spliceosomal introns is a well-studied fundamental reaction in molecular biology. Apart from this prominent example, there are surprisingly few occurrences of branched nucleotides or even 2',5'-phosphodiester bonds in biology. We recently described a new ribozyme, the GIR1 branching ribozyme, which catalyzes the formation of a tiny lariat that caps an mRNA. This new example together with work on artificial branching ribozymes and deoxyribozymes shows that branching is facile and points to the possibility that branching reactions could be more prevalent than previously recognized.

An mRNA is capped by a 2', 5' lariat catalyzed by a group I-like ribozyme

Twin-ribozyme introns are formed by two ribozymes belonging to the group I family and occur in some ribosomal RNA transcripts. The group I-like ribozyme, GIR1, liberates the 5' end of a homing endonuclease messenger RNA in the slime mold Didymium iridis. We demonstrate that this cleavage occurs by a transesterification reaction with the joining of the first and the third nucleotide of the messenger by a 2',5'-phosphodiester linkage. Thus, a group I-like ribozyme catalyzes an RNA branching reaction similar to the first step of splicing in group II introns and spliceosomal introns. The resulting short lariat, by forming a protective 5' cap, might have been useful in a primitive RNA world.

Hydrolytic cleavage by a group I intron ribozyme is dependent on RNA structures not important for splicing

DiGIR2 is the group I splicing-ribozyme of the mobile twin-ribozyme intron Dir.S956-1, present in Didymium nuclear ribosomal DNA. DiGIR2 is responsible for intron excision, exon ligation, 3'-splice site hydrolysis, and full-length intron RNA circle formation. We recently reported that DiGIR2 splicing (intron excision and exon ligation) competes with hydrolysis and subsequent full-length intron circularization. Here we present experimental evidence that hydrolysis at the 3'-splice site in DiGIR2 is dependent on structural elements within the P9 subdomain not involved in splicing. Whereas the GCGA tetra-loop in P9b was found to be important in hydrolytic cleavage, probably due to tertiary RNA-RNA interactions, the P9.2 hairpin structure was found to be essential for hydrolysis. The most important positions in P9.2 include three adenosines in the terminal loop (L9.2) and a consensus kink-turn motif in the proximal stem. We suggest that the L9.2 adenosines and the kink-motif represent key regulatory elements in the splicing and hydrolytic reaction pathways.

DiGIR1 and NaGIR1: naturally occurring group I-like ribozymes with unique core organization and evolved biological role

The group I-like ribozyme GIR1 is a unique example of a naturally occurring ribozyme with an evolved biological function. GIR1 generates the 5'-end of a nucleolar encoded messenger RNA involved in intron mobility. GIR1 is found as a cis-cleaving ribozyme within two very different rDNA group I introns (twin-ribozyme introns) in distantly related organisms. The Didymium GIR1 (DiGIR1) and Naegleria GIR1 (NaGIR1) share fundamental features in structural organization and reactivity, and display significant differences when compared to the related group I splicing ribozymes. GIR1 lacks the characteristic P1 segment present in all group I splicing ribozymes, it has a novel core organization, and it catalyses two site-specific hydrolytic cleavages rather than splicing. DiGIR1 and NaGIR1 appear to have originated from eubacterial group I introns in order to fulfil a common biological challenge: the expression of a protein encoding gene in a nucleolar context.

A rapid and simple nonradioactive method for in vitro testing of ribozyme activity

Ribozymes that target specific messenger RNA transcripts are powerful tools in the emerging fields of functional genomics, proteomics, and metabolomics. We have found that successful in vitro testing greatly increases the likelihood of producing ribozymes with good efficacy in living cells. A rapid and simple nonradioactive method for systematic in vitro testing of ribozyme-cleaving activity is reported. Ribozymes are synthesized enzymatically from double-stranded DNA (dsDNA) oligonucleotides without vector cloning. Substrate target DNA template is cloned into a vector flanked with SP6 and T7 promoters at multiple cloning sites that permit colorimetric screening and ampicillin selection, enhancing the efficiency of the cloning procedure. Ribozyme cleavage products are satisfactorily resolved on 2.0% NuSieve 3:1 agarose (FMC Products, Rockland, ME)/formaldehyde gels by electrophoresis. This method avoids the preparation of polyacrylamide gels. Using this procedure, the ribozyme, target substrate RNA, and ribozyme cleavage products are all easily detected by ethidium bromide staining. Resolution and detection are fast and simple, eliminating the need for either polyacrylamide gel analysis or radiolabeling. The use of RNase inhibitors in the assays is also assessed and discussed.

Flanking sequences with an essential role in hydrolysis of a self-cleaving group I-like ribozyme

DiGIR1 is a group I-like ribozyme derived from the mobile twin ribozyme group I intron DiSSU1 in the nuclear ribosomal DNA of the myxomycete Didymium iridis. This ribozyme is responsible for intron RNA processing in vitro and in vivo at two internal sites close to the 5'-end of the intron endo-nuclease open reading frame and is a unique example of a group I ribozyme with an evolved biological function. DiGIR1 is the smallest functional group I ribozyme known from nature and has an unusual core organization including the 6 bp P15 pseudoknot. Here we report results of functional and structural analyses that identify RNA elements critical for hydrolysis outside the DiGIR1 ribozyme core moiety. Results from deletion analysis, disruption/compensation mutagenesis and RNA structure probing analysis all support the existence of two new segments, named P2 and P2.1, involved in the hydrolysis of DiGIR1. Significant decreases in the hydrolysis rate, k (obs), were observed in disruption mutants involving both segments. These effects were restored by compensatory base pairing mutants. The possible role of P2 is to tether the ribozyme core, whereas P2.1 appears to be more directly involved in catalysis.

Structure, stability and function of RNA pseudoknots involved in stimulating ribosomal frameshifting

Programmed -1 ribosomal frameshifting has become the subject of increasing interest over the last several years, due in part to the ubiquitous nature of this translational recoding mechanism in pathogenic animal and plant viruses. All cis-acting frameshift signals encoded in mRNAs are minimally composed of two functional elements: a heptanucleotide "slippery sequence" conforming to the general form X XXY YYZ, followed by an RNA structural element, usually an H-type RNA pseudoknot, positioned an optimal number of nucleotides (5 to 9) downstream. The slippery sequence itself promotes a low level ( approximately 1 %) of frameshifting; however, downstream pseudoknots stimulate this process significantly, in some cases up to 30 to 50 %. Although the precise molecular mechanism of stimulation of frameshifting remains poorly understood, significant advances have been made in our knowledge of the three-dimensional structures, thermodynamics of folding, and functional determinants of stimulatory RNA pseudoknots derived from the study of several well-characterized frameshift signals. These studies are summarized here and provide new insights into the structural requirements and mechanism of programmed -1 ribosomal frameshifting.

Identification of the minimal essential RNA sequences responsible for site-specific targeting of the Leishmania RNA virus 1-4 capsid endoribonuclease

The Leishmania RNA virus 1-4 capsid protein possesses an endoribonuclease activity responsible for single-site-specific cleavage within the 450-nucleotide 5' untranslated region of its own viral RNA transcript. To characterize the minimal essential RNA determinants required for site-specific cleavage, mutated RNA transcripts were examined for susceptibility to cleavage by the virus capsid protein in an in vitro assay. Deletion analyses revealed that all determinants necessary for accurate cleavage are encoded in viral nucleotides 249 to 342. Nuclease mapping and site-specific mutagenesis of the minimal RNA sequence defined a stem-loop structure that is located 40 nucleotides upstream from the cleavage site (nucleotide 320) and that is essential for accurate RNA cleavage. Abrogation of cleavage by disruption of base pairing within the stem-loop was reversed through the introduction of complementary nucleotide substitutions that reestablished the structure. We also provide evidence that divalent cations, essential components of the cleavage reaction, stabilized the stem-loop structure in solution. That capsid-specific antiserum eliminated specific RNA cleavage provides further evidence that the virus capsid gene encodes the essential endoribonuclease activity.

Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure

An improved dynamic programming algorithm is reported for RNA secondary structure prediction by free energy minimization. Thermodynamic parameters for the stabilities of secondary structure motifs are revised to include expanded sequence dependence as revealed by recent experiments. Additional algorithmic improvements include reduced search time and storage for multibranch loop free energies and improved imposition of folding constraints. An extended database of 151,503 nt in 955 structures? determined by comparative sequence analysis was assembled to allow optimization of parameters not based on experiments and to test the accuracy of the algorithm. On average, the predicted lowest free energy structure contains 73 % of known base-pairs when domains of fewer than 700 nt are folded; this compares with 64 % accuracy for previous versions of the algorithm and parameters. For a given sequence, a set of 750 generated structures contains one structure that, on average, has 86 % of known base-pairs. Experimental constraints, derived from enzymatic and flavin mononucleotide cleavage, improve the accuracy of structure predictions.

I-NjaI, a nuclear intron-encoded homing endonuclease from Naegleria, generates a pentanucleotide 3' cleavage-overhang within a 19 base-pair partially symmetric DNA recognition site

Different species of the amoebo-flagellate Naegleria harbor optional group I introns in the nuclear ribosomal DNA that contain open reading frames. Intron proteins from Naegleria jamiesoni, Naegleria andersoni, and Naegleria italica (named I-NjaI, I-NanI and I-NitI, respectively) were expressed in Escherichia coli and found to be isoschizomeric homing endonucleases that specifically recognize and cleave intron-lacking homologous alleles of ribosomal DNA. The I-NjaI endonuclease was affinity purified, characterized in more detail, and found to generate five-nucleotide 3' staggered ends at the intron insertion site which differs from the ends generated by all other known homing endonucleases. The recognition site was delimited and found to cover an approximately 19 base-pair partially symmetric sequence spanning both the cleavage site and the intron insertion site. The palindromic feature was supported by mutational analysis of the target DNA. All single-site substitutions within the recognition sequence were cleaved by the purified I-NjaI endonuclease, but at different efficiencies. The center of symmetry and cleavage was found to be completely degenerate in specificity, which resembles that of the subclass IIW bacterial restriction enzymes.

Group I twintrons: genetic elements in myxomycete and schizopyrenid amoeboflagellate ribosomal DNAs

Protists are unicellular eukaryotes which represent a significant fraction of the global biodiversity. The myxomycete Didymium and the schizopyrenid amoeboflagellate Naegleria are distantly related protists. However, we have noted several striking similarities in life cycle, cell morphology, and ribosomal DNA organization between these organisms. Both have multicopy nuclear extrachromosomal ribosomal DNAs. Here the small subunit ribosomal RNA genes are interrupted by an optional group I twintron, a novel category among the group I introns. Group I twintrons are mobile self-splicing introns of 1.3-1.4 kb in size, with a complex organization at the RNA level. A group I twintron consists of two distinct ribozymes (catalytic RNAs) with different functions in RNA processing, and an open reading frame encoding a functional homing endonuclease--all with prospects of application as molecular tools in biotechnology. Updated RNA secondary structure models of group I twintrons, as well as an example of in vitro ribozyme activity, are presented. We suggest that the group I twintrons have been independently established in myxomycetes and schizopyrenid amoeboflagellates by horizontal gene transfer due to a combination of the phagocytotic behavior in natural environments and the extrachromosomal multicopy nature of ribosomal DNA.