Visualizing the solvent-inaccessible core of a group II intron ribozyme

Group II introns are well recognized for their remarkable catalytic capabilities, but little is known about their three-dimensional structures. In order to obtain a global view of an active enzyme, hydroxyl radical cleavage was used to define the solvent accessibility along the backbone of a ribozyme derived from group II intron ai5gamma. These studies show that a highly homogeneous ribozyme population folds into a catalytically compact structure with an extensively internalized catalytic core. In parallel, a model of the intron core was built based on known tertiary contacts. Although constructed independently of the footprinting data, the model implicates the same elements for involvement in the catalytic core of the intron.     

Elimination of a group II intron from a plastid gene causes a mutant phenotype

Group II introns are found in bacteria and cell organelles (plastids, mitochondria) and are thought to represent the evolutionary ancestors of spliceosomal introns. It is generally believed that group II introns are selfish genetic elements that do not have any function. Here, we have scrutinized this assumption by analyzing two group II introns that interrupt a plastid gene (ycf3) involved in photosystem assembly. Using stable transformation of the plastid genome, we have generated mutant plants that lack either intron 1 or intron 2 or both. Interestingly, the deletion of intron 1 caused a strong mutant phenotype. We show that the mutants are deficient in photosystem I and that this deficiency is directly related to impaired ycf3 function. We further show that, upon deletion of intron 1, the splicing of intron 2 is strongly inhibited. Our data demonstrate that (i) the loss of a group II intron is not necessarily phenotypically neutral and (ii) the splicing of one intron can depend on the presence of another.     

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

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

Occurrence of matK in a trnK group II intron in charophyte green algae and phylogeny of the Characeae

A group II intron containing the matK gene, which encodes a splicing-associated maturase, was found in the trnK (lysine tRNA) exon in the chloroplast genome of the six extant genera of green algae in the family Characeae, which among green algae are the sister group to embryophytes (land plants). The characean trnK intron (～2.5 kilobases [kb]) and matK ORF (～1.5 kb) are comparable in size to the intron and ORF of land plants, in which they are similarly found inserted in the trnK exon. Domain X, a sequence of conserved amino acid residues within matK, occurs in the Characeae. Phylogenetic analysis using maximum likelihood (GTR + I + gamma likelihood model) and parsimony (branch and bound search) yielded one tree with high bootstrap support for all branches. The matK tree was congruent with the rbcL tree for the same taxa. The number and proportion of informative sites was higher in matK (501, 31% of matK sequence) compared to rbcL (122, 10%). Characeae branch lengths were on average more than five times longer for matK compared to rbcL and provided better resolution within the Characeae. These findings along with recent genomic analyses demonstrate that the intron and matK invaded the chloroplast genome of green algae prior to the evolution of land plants.     

Guiding ribozyme cleavage through motif recognition: the mechanism of cleavage site selection by a group ii intron ribozyme

The mechanism by which group II introns cleave the correct phosphodiester linkage was investigated by studying the reaction of mutant substrates with a ribozyme derived from intron ai5gamma. While fidelity was found to be quite high in most cases, a single mutation on the substrate (+1C) resulted in a dramatic loss of fidelity. When this mutation was combined with a second mutation that induces a bulge in the exon binding site 1/intron binding site 1 (EBS1/IBS1) duplex, the base-pairing register of the EBS1/IBS1 duplex was shifted and the cleavage site moved to a downstream position on the substrate. Conversely, when mismatches were incorporated at the EBS1/IBS1 terminus, the duplex was effectively truncated and cleavage occurred at an upstream site. Taken together, these data demonstrate that the cleavage site of a group II intron ribozyme can be tuned at will by manipulating the thermodynamic stability and structure of the EBS1/IBS1 pairing. The results are consistent with a model in which the cleavage site is not designated through recognition of specific nucleotides (such as the 5'-terminal residue of EBS1). Instead, the ribozyme detects a structure at the junction between single and double-stranded residues on the bound substrate. This finding explains the puzzling lack of phylogenetic conservation in ribozyme and substrate sequences near group II intron target sites.     

An obligate intermediate along the slow folding pathway of a group II intron ribozyme

Most RNA molecules collapse rapidly and reach the native state through a pathway that contains numerous traps and unproductive intermediates. The D135 group II intron ribozyme is unusual in that it can fold slowly and directly to the native state, despite its large size and structural complexity. Here we use hydroxyl radical footprinting and native gel analysis to monitor the timescale of tertiary structure collapse and to detect the presence of obligate intermediates along the folding pathway of D135. We find that structural collapse and native folding of Domain 1 precede assembly of the entire ribozyme, indicating that D1 contains an on-pathway intermediate to folding of the D135 ribozyme. Subsequent docking of Domains 3 and 5, for which D1 provides a preorganized scaffold, appears to be very fast and independent of one another. In contrast to other RNAs, the D135 ribozyme undergoes slow tertiary collapse to a compacted state, with a rate constant that is also limited by the formation D1. These findings provide a new paradigm for RNA folding and they underscore the diversity of RNA biophysical behaviors.     

A folding control element for tertiary collapse of a group II intron ribozyme

Ribozymes derived from the group II intron ai5gamma collapse to a compact intermediate, folding to the native state through a slow, direct pathway that is unperturbed by kinetic traps. Molecular collapse of ribozyme D135 requires high magnesium concentrations and is thought to involve a structural element in domain 1 (D1). We used nucleotide analog interference mapping, in combination with nondenaturing gel electrophoresis, to identify RNA substructures and functional groups that are essential for D135 tertiary collapse. This revealed that the most crucial atoms for compaction are located within a small section of D1 that includes the kappa and zeta elements. This small substructure controls specific collapse of the molecule and, in later steps of the folding pathway, it forms the docking site for catalytic D5. In this way, the stage is set for proper active site formation during the earliest steps of ribozyme folding.     

Origin and evolution of the chloroplast trnK (matK) intron: a model for evolution of group II intron RNA structures

The trnK intron of plants encodes the matK open reading frame (ORF), which has been used extensively as a phylogenetic marker for classification of plants. Here we examined the evolution of the trnK intron itself as a model for group II intron evolution in plants. Representative trnK intron sequences were compiled from species spanning algae to angiosperms, and four introns were newly sequenced. Phylogenetic analyses showed that the matK ORFs belong to the ML (mitochondrial-like) subclass of group II intron ORFs, indicating that they were derived from a mobile group II intron of the class. RNA structures of the introns were folded and analyzed, which revealed progressive RNA structural deviations and degenerations throughout plant evolution. The data support a model in which plant organellar group II introns were derived from bacterial-like introns that had "standard" RNA structures and were competent for self-splicing and mobility and that subsequently the ribozyme structures degenerated to ultimately become dependent upon host-splicing factors. We propose that the patterns of RNA structure evolution seen for the trnK intron will apply to the other group II introns in plants.     

Optimizing the substrate specificity of a group I intron ribozyme

Group I ribozymes can repair mutant RNAs via trans-splicing. Unfortunately, substrate specificity is quite low for the trans-splicing reaction catalyzed by the group I ribozyme from Tetrahymenathermophila. We have used a systematic approach based on biochemical knowledge of the function of the Tetrahymena ribozyme to optimize its ability to discriminate against nonspecific substrates in vitro. Ribozyme derivatives that combine a mutation which indirectly slows down the rate of the chemical cleavage step by weakening guanosine binding with additional mutations that weaken substrate binding have greatly enhanced specificity with short oligonucleotide substrates and an mRNA fragment derived from the p53 gene. Moreover, compared to the wild-type ribozyme, reaction of a more specific ribozyme with targeted substrates is much less sensitive to the presence of nonspecific RNA competitors. These results demonstrate how a detailed understanding of the biochemistry of a catalytic RNA can facilitate the design of customized ribozymes with improved properties for therapeutic applications.     

Splice site selection and role of the lariat in a group II intron.

The structural elements involved in 5' and 3' splice site (SS) selection in a group II intron were analyzed. While 5' SS selection appears to be defined by only one element, the EBS1-IBS1 pairing, four distinct structural components contribute to 3' SS selection, one of which being analogous to the "internal guide sequence" described for group I introns. Moreover, some of the mutants analyzed during this study induce efficient 5' SS hydrolysis and suggest how 5' SS transesterification is selected against hydrolysis. Finally, the lariat structure was found to accelerate both steps of splicing, suggesting that it "locks" the ribozyme in an active configuration.     

Restriction for gene insertion within the Lactococcus lactis Ll.LtrB group II intron

The Ll.LtrB intron, from the low G+C gram-positive bacterium Lactococcus lactis, was the first bacterial group II intron shown to splice and mobilize in vivo. The detailed retrohoming and retrotransposition pathways of Ll.LtrB were studied in both L. lactis and Escherichia coli. This bacterial retroelement has many features that would make it a good gene delivery vector. Here we report that the mobility efficiency of Ll.LtrB expressing LtrA in trans is only slightly affected by the insertion of fragments <100 nucleotides within the loop region of domain IV. In contrast, Ll.LtrB mobility efficiency is drastically decreased by the insertion of foreign sequences >1 kb. We demonstrate that the inhibitory effect caused by the addition of expression cassettes on Ll.LtrB mobility efficiency is not sequence specific, and not due to the expression, or the toxicity, of the cargo genes. Using genetic screens, we demonstrate that in order to maintain intron mobility, the loop region of domain IV, more specifically domain IVb, is by far the best region to insert foreign sequences within Ll.LtrB. Poisoned primer extension and Northern blot analyses reveal that Ll.LtrB constructs harboring cargo sequences splice less efficiently, and show a significant reduction in lariat accumulation in L. lactis. This suggests that cargo-containing Ll.LtrB variants are less stable. These results reveal the potential, yet limitations, of the Ll.LtrB group II intron to be used as a gene delivery vector, and validate the random insertion approach described in this study to create cargo-containing Ll.LtrB variants that are mobile.     

Constraints on the evolution of plastid introns: the group II intron in the gene encoding tRNA-Val(UAC).

The evolution of the group II intron in the plastid gene encoding tRNA(Val)UAC (trnV) from seven plant taxa was studied by aligning secondary and other structural features. Levels of evolutionary divergence between six angiosperms and a liverwort, Marchantia polymorpha, were compared for the six domains commonly demonstrated for group II introns and were shown to be statistically heterogeneous. Evolutionary rates varied substantially among various domains and other features. Domain II showed the highest evolutionary rate, approaching the synonymous substitution rate reported for cpDNA-encoded genes, while domain VI and the helix and loop region bearing EBS1 evolved at rates similar to those for nonsynonymous substitutions of a number of cpDNA-encoded genes. The minimum free-energy structure of domain I varied among the seven taxa, suggesting that possible protein-RNA or tertiary interactions are important for intron processing.     

Fluorescence and solution NMR study of the active site of a 160-kDa group II intron ribozyme

We have reconstructed the group II intron from Pylaiella littoralis (PL) into a hydrolytic ribozyme, comprising domains 1-3 (D123) connected in cis plus domain 5 (D5) supplied in trans that efficiently cleaves spliced exon substrates. Using a novel gel-based fluorescence assay and nuclear magnetic resonance (NMR) spectroscopy, we monitored the direct binding of D5 to D123, characterized the kinetics of the spliced exon hydrolysis reaction (which is mechanistically analogous to the reverse of the second catalytic step of splicing), and identified the binding surface of D123 on D5. This PL ribozyme acts as an RNA endonuclease even at low monovalent (100 mM KCl) and divalent ion concentrations (1-10 mM MgCl(2)). This is in contrast to other group II intron ribozyme systems that require high levels of salt, making NMR analysis problematic. D5 binds tightly to D123 with a K(d) of 650 +/- 250 nM, a K(m) of approximately 300 nM, and a K(cat) of 0.02 min(-1) under single turnover conditions. Within the approximately 160-kDa D123-D5 binary complex, site-specific binding to D123 leads to dramatic chemical shift perturbation of residues localized to the tetraloop and internal bulge within D5, suggesting a structural switch model for D5-assisted splicing. This minimal ribozyme thus recapitulates the essential features of the reverse of the second catalytic step and represents a well-behaved system for ongoing high-resolution structural work to complement folding and catalytic functional studies.     

Folding of the group I intron ribozyme from the 26S rRNA gene of Candida albicans

Preincubation of the group I intron Ca.LSU from Candida albicans at 37°C in the absence of divalent cations results in partial folding of this intron. This is indicated by increased resistance to T1 ribonuclease cleavage of many G residues in most local helices, including P4-P6, as well as the non-local helix P7, where the G binding site is located. These changes correlate with increased gel mobility and activation of catalysis by precursor RNA containing this intron after preincubation. The presence of divalent cations or spermidine during preincubation results in formation of the predicted helices, as indicated by protection of additional G residues. However, addition of these cations during preincubation of the precursor RNA alters its gel mobility and eliminates the preincubation activation of precursor RNA seen in the absence of cations. These results suggest that, in the presence of divalent cations or spermidine, Ca.LSU folds into a more ordered, stable but misfolded conformation that is less able to convert into the catalytically active form than the ribozyme preincubated without cations. These results indicate that, like the group I intron of Tetrahymena, multiple folding pathways exist for Ca.LSU. However, it appears that the role cations play in the multiple folding pathways leading to the catalytically active form may differ between folding of these two group I introns.     

A map of the binding site for catalytic domain 5 in the core of a group II intron ribozyme.

Group II introns are ribozymes with a complex tertiary architecture that is of great interest as a model for RNA folding. Domain 5 (D5) is a highly conserved region of the intron that is considered one of the most critical structures in the catalytic core. Despite its central importance, the means by which D5 interacts with other core elements is unclear. To obtain a map of potential interaction sites, dimethyl sulfate was used to footprint regions of the intron that are involved in D5 binding. These studies were complemented by measurements of D5 binding to a series of truncated intron derivatives. In this way, the minimal region of the intron required for strong D5 association was defined and the sites most likely to represent thermodynamically significant positions of tertiary contact were identified. These studies show that ground-state D5 binding is mediated by tertiary contacts to specific regions of D1, including a tetraloop receptor and an adjacent three-way junction. In contrast, D2 and D3 are not found to stabilize D5 association. These data highlight the significance of D1-D5 interactions and will facilitate the identification of specific tertiary contacts between them.