Mutational analysis of conserved nucleotides in a self-splicing group I intron

Optimization of RNA Aptamer SELEX Methods: Improved Aptamer Transcript 3'-End Homogeneity, PAGE Purification Yield, and Target-Bound Aptamer RNA Recovery

Since its inception in the early 1990s, SELEX remains the gold standard for discovering RNA aptamers specific for proteins and small molecules. The SELEX process has undergone countless modifications and now encompasses a breadth of innovative selection schemes to pare an aptamer library toward target-specific aptamers. Common to all these RNA aptamer SELEX processes are the steps for the preparation of DNA template and in vitro transcription of aptamer RNA. These steps have remained mostly unchanged over the past three decades and would benefit from optimization. We focused on three key areas: improving the homogeneity of in vitro transcribed aptamer RNA, increasing the efficiency of in vitro transcribed aptamer RNA purification by PAGE, and improving the quality of target-bound aptamer RNA recovered during SELEX. Together, these optimizations contribute toward a more efficient SELEX process and are applicable to both protein-based and cell-based RNA aptamer selections.

Electrochemical Biosensor Based on Well-Dispersed Boron Nitride Colloidal Nanoparticles and DNA Aptamers for Ultrasensitive Detection of Carbendazim

A selective electrochemical biosensor was developed for detecting carbendazim (CBZ) based on well-dispersed colloidal boron nitride (BN) nanocrystals and gold nanoparticles (Au NPs). BN was synthesized by "solvent cutting" to modify a glassy carbon electrode (GCE), and Au NPs were then electrodeposited. A single-stranded oligonucleotide with methylene blue (MB) was modified to the electrode surface through gold-sulfur bonds. A double-stranded DNA was formed in the presence of an aptamer. The aptamer chain can specifically bind to the target CBZ. When the aptamer binds to CBZ, the electroactive substance MB labeled at one end of the complementary chain can effectively contact the electrode surface. Detection of CBZ is realized by simultaneously monitoring the MB signal enhancement. The CBZ concentration was determined in a wide linearity range from 0.1 ng mL^-1 to 100 μg mL^-1, with a low detection limit of 0.019 ng mL^-1. This biosensor exhibited excellent selectivity and acceptable repeatability and was applied in cucumber, kiwifruit, and water samples with good recoveries, demonstrating that the strategy has remarkable potential and offers a good platform for CBZ detection.

A Phylogenetic Approach to Structural Variation in Organization of Nuclear Group I Introns and Their Ribozymes

Nuclear group I introns are restricted to the ribosomal DNA locus where they interrupt genes for small subunit and large subunit ribosomal RNAs at conserved sites in some eukaryotic microorganisms. Here, the myxomycete protists are a frequent source of nuclear group I introns due to their unique life strategy and a billion years of separate evolution. The ribosomal DNA of the myxomycete Mucilago crustacea was investigated and found to contain seven group I introns, including a direct repeat-containing intron at insertion site S1389 in the small subunit ribosomal RNA gene. We collected, analyzed, and compared 72 S1389 group IC1 introns representing diverse myxomycete taxa. The consensus secondary structure revealed a conserved ribozyme core, but with surprising sequence variations in the guanosine binding site in segment P7. Some S1389 introns harbored large extension sequences in the peripheral region of segment P9 containing direct repeat arrays. These repeats contained up to 52 copies of a putative internal guide sequence motif. Other S1389 introns harbored homing endonuclease genes in segment P1 encoding His-Cys proteins. Homing endonuclease genes were further interrupted by small spliceosomal introns that have to be removed in order to generate the open reading frames. Phylogenetic analyses of S1389 intron and host gene indicated both vertical and horizontal intron transfer during evolution, and revealed sporadic appearances of direct repeats, homing endonuclease genes, and guanosine binding site variants among the myxomycete taxa.

Heterodimerization of Group I Ribozymes Enabling Exon Recombination through Pairs of Cooperative trans-Splicing Reactions

Group I (GI) self-splicing ribozymes are attractive tools for biotechnology and synthetic biology. Several trans-splicing and related reactions based on GI ribozymes have been developed for the purpose of recombining their target mRNA sequences. By combining trans-splicing systems with rational modular engineering of GI ribozymes it was possible to achieve more complex editing of target RNA sequences. In this study we have developed a cooperative trans-splicing system through rational modular engineering with use of dimeric GI ribozymes derived from the Tetrahymena group I intron ribozyme. The resulting pairs of ribozymes exhibited catalytic activity depending on their selective dimerization. Rational modular redesign as performed in this study would facilitate the development of sophisticated regulation of double or multiple trans-splicing reactions in a cooperative manner.

Use of a Fluorescent Aptamer RNA as an Exonic Sequence to Analyze Self-Splicing Ability of aGroup I Intron from Structured RNAs

Group I self-splicing intron constitutes an important class of functional RNA molecules that can promote chemical transformation. Although the fundamental mechanism of the auto-excision from its precursor RNA has been established, convenient assay systems for its splicing activity are still useful for a further understanding of its detailed mechanism and of its application. Because some host RNA sequences, to which group I introns inserted form stable three-dimensional (3D) structures, the effects of the 3D structures of exonic elements on the splicing efficiency of group I introns are important but not a fully investigated issue. We developed an assay system for group I intron self-splicing by employing a fluorescent aptamer RNA (spinach RNA) as a model exonic sequence inserted by the Tetrahymena group I intron. We investigated self-splicing of the intron from spinach RNA, serving as a model exonic sequence with a 3D structure.

Engineering a family of synthetic splicing ribozymes

Controlling RNA splicing opens up possibilities for the synthetic biologist. The Tetrahymena ribozyme is a model group I self-splicing ribozyme that has been shown to be useful in synthetic circuits. To create additional splicing ribozymes that can function in synthetic circuits, we generated synthetic ribozyme variants by rationally mutating the Tetrahymena ribozyme. We present an alignment visualization for the ribozyme termed as structure information diagram that is similar to a sequence logo but with alignment data mapped on to secondary structure information. Using the alignment data and known biochemical information about the Tetrahymena ribozyme, we designed synthetic ribozymes with different primary sequences without altering the secondary structure. One synthetic ribozyme with 110 nt mutated retained 12% splicing efficiency in vivo. The results indicate that our biochemical understanding of the ribozyme is accurate enough to engineer a family of active splicing ribozymes with similar secondary structure but different primary sequences.

The active site of the ribosome is composed of two layers of conserved nucleotides with distinct roles in peptide bond formation and peptide release

Peptide bond formation and peptide release are catalyzed in the active site of the large subunit of the ribosome where universally conserved nucleotides surround the CCA ends of the peptidyl- and aminoacyl-tRNA substrates. Here, we describe the use of an affinity-tagging system for the purification of mutant ribosomes and analysis of four universally conserved nucleotides in the innermost layer of the active site: A2451, U2506, U2585, and A2602. While pre-steady-state kinetic analysis of the peptidyl transferase activity of the mutant ribosomes reveals substantially reduced rates of peptide bond formation using the minimal substrate puromycin, their rates of peptide bond formation are unaffected when the substrates are intact aminoacyl-tRNAs. These mutant ribosomes do, however, display substantial defects in peptide release. These results reveal a view of the catalytic center in which an inner shell of conserved nucleotides is pivotal for peptide release, while an outer shell is responsible for promoting peptide bond formation.

Large phenotype jumps in biomolecular evolution

By defining the phenotype of a biopolymer by its active three-dimensional shape, and its genotype by its primary sequence, we propose a model that predicts and characterizes the statistical distribution of a population of biopolymers with a specific phenotype that originated from a given genotypic sequence by a single mutational event. Depending on the ratio g(0) that characterizes the spread of potential energies of the mutated population with respect to temperature, three different statistical regimes have been identified. We suggest that biopolymers found in nature are in a critical regime with g(0) approximately 1-6, corresponding to a broad, but not too broad, phenotypic distribution resembling a truncated Lévy flight. Thus the biopolymer phenotype can be considerably modified in just a few mutations. The proposed model is in good agreement with the experimental distribution of activities determined for a population of single mutants of a group-I ribozyme.

Evolution in vitro: analysis of a lineage of ribozymes

BACKGROUND

Catalytic RNAs, or ribozymes, possessing both a genotype and a phenotype, are ideal molecules for evolution experiments in vitro. A large, heterogeneous pool of RNAs can be subjected to multiple rounds of selection, amplification and mutation, leading to the development of variants that have some desired phenotype. Such experiments allow the investigator to correlate specific genetic changes with quantifiable alterations of the catalytic properties of the RNA. In addition, patterns of evolutionary change can be discerned through a detailed examination of the genotypic composition of the evolving RNA population.

RESULTS

Beginning with a pool of 10(13) variants of the Tetrahymena ribozyme, we carried out in vitro evolution experiments that led to the generation of ribozymes with the ability to cleave an RNA substrate in the presence of Ca2+ ions, an activity that does not exist for the wild-type molecule. Over the course of 12 generations, a seven-error variant emerged that has substantial Ca(2+)-dependent RNA-cleavage activity. Advantageous mutations increased in frequency in the population according to three distinct dynamics--logarithmic, linear and transient. Through a comparative analysis of 31 individual variants, we infer how certain mutations influence the catalytic properties of the ribozyme.

CONCLUSIONS

In vitro evolution experiments make it possible to elucidate important aspects of both evolutionary biology and structural biochemistry on a reasonable short time scale.

Rearrangement of a stable RNA secondary structure during VS ribozyme catalysis

The Neurospora VS ribozyme recognizes and cleaves a substrate RNA that contains a GC-rich stem loop. In contrast to most RNA secondary structures that are stable during tertiary or quaternary folding, this substrate undergoes extensive ribozyme-induced rearrangement in the presence of magnesium in which the base pairings of at least seven of the ten nucleotides in the stem are changed. This conformational switch is essential for catalytic activity with the wild-type substrate and creates a metal-binding secondary structure motif near the cleavage site. Base pair rearrangement is accompanied by bulging a cytosine from the middle of the stem, indicating that ribozymes may perform base flipping, an activity previously observed only with protein enzymes that modify DNA.

Characterization of P8 and J8/7 elements in the conserved core of the tetrahymena group I intron ribozyme

The universally conserved core region in the group I intron ribozymes is responsible for its catalytic activity. The structural elements in this region have been known to organize the active site of this class of ribozymes. However, it has been unclear whether all elements are requisite or some elements are dispensable for conducting the catalysis. To investigate the necessity of these elements in the catalysis, we prepared and examined a series of mutants having a nick or deletion in these elements. In this report, we show that two elements, P8 and 5' portion of J8/7, are nonessential for activity.

A core folding model for catalysis by the hammerhead ribozyme accounts for its extraordinary sensitivity to abasic mutations

Introducing abasic nucleotides at each of 13 positions in the conserved core of the hammerhead ribozyme causes a large decrease in the extent of catalysis [Peracchi, A., et al. (1996) Proc. Natl. Acad. Sci. U.S.A. 93, 11522]. This extreme sensitivity to structural defects is in contrast to the behavior of protein enzymes and larger ribozymes. Several additional differences in the behavior of the hammerhead relative to that of protein enzymes and larger ribozymes are described herein. The deleterious effects of the abasic mutations are not relieved by lowering the temperature, by increasing the concentration of monovalent or divalent metal ions, or by adding polyamines, in contrast to effects observed with protein enzymes and large RNA enzymes. In addition, the abasic mutations do not significantly weaken substrate binding. These results and previous observations are all accounted for by a "core folding" model in which the stable ground state structure of the hammerhead ribozyme complexed with the substrate is a partially folded state that must undergo an additional folding event to achieve its catalytic conformation. We propose that the peculiar behavior of the hammerhead arises because the limited structural interconnections in a small RNA enzyme do not allow the ground state to stably adopt the catalytic conformation; within the globally folded catalytic conformation, limited structural interconnections may further impair catalysis by hampering the precise alignment of active site functional groups. This behavior represents a basic manifestation of the well-recognized interconnection between folding and catalysis.

A preorganized active site in the crystal structure of the Tetrahymena ribozyme

Group I introns possess a single active site that catalyzes the two sequential reactions of self-splicing. An RNA comprising the two domains of the Tetrahymena thermophila group I intron catalytic core retains activity, and the 5.0 angstrom crystal structure of this 247-nucleotide ribozyme is now described. Close packing of the two domains forms a shallow cleft capable of binding the short helix that contains the 5' splice site. The helix that provides the binding site for the guanosine substrate deviates significantly from A-form geometry, providing a tight binding pocket. The binding pockets for both the 5' splice site helix and guanosine are formed and oriented in the absence of these substrates. Thus, this large ribozyme is largely preorganized for catalysis, much like a globular protein enzyme.

Smaller, faster ribozymes reveal the catalytic core of Neurospora VS RNA

We have investigated the structural requirements for cis-cleavage of the VS ribozyme by designing deletions, substitutions, and circular permutations based on the secondary structure model. Four of the six helices predicted in the model have been shortened, resulting in self-cleaving RNAs of only 121 to 126 nucleotides. Remarkably, the shorter ribozymes exhibit a 30 to 40-fold faster cis-cleavage rate. The increase in activity results from disrupting an inhibitory helix whose 5' side contains bases upstream of the cleavage site, and from constructing a circular permutation that tethers the helix containing the cleavage site to a shortened version of the rest of the ribozyme. The non-essential regions identified by the deletions map to the periphery of a recently proposed structure model, revealing a central ribozyme core that contains the essential structural elements required for activity of the VS ribozyme.

Crystals by design: a strategy for crystallization of a ribozyme derived from the Tetrahymena group I intron

Recently, the 2.8 A crystal structure of one domain of the self-splicing Tetrahymena group I intron was reported. Although it revealed much about RNA tertiary interactions, it contained only half of the active site. We have now designed a series of larger molecules that contain about 70% of the intron and all of the catalytic core. These RNAs were efficient in cleavage of a substrate RNA, consisting of the approximately 100 nucleotides from the 5' end of the intron, at a site corresponding to the 5' splice site. A sparse matrix was designed specifically for large RNAs and used to screen for preliminary crystallization conditions. Of the six RNAs initially tested, five were crystallized in this initial trial. Two of these crystals were further examined. The first diffracted X-rays to only approximately 16 A resolution, even when the crystal were very large. The second diffracted as high as 3.5 A, but the crystals were twinned and therefore unusable for structural studies. Site-specific mutagenesis was performed on the latter RNA to disrupt interactions that might have been responsible for the twinning. One of these mutant RNAs produced large, single, diffraction-quality crystals. The crystals belong to the tetragonal space group P42212 and have large unit cell dimensions, a=b=178 A and c=199 A. Thus, by variation of both sequence elements and crystallization conditions, crystals of a 247 nucleotide catalytic RNA were obtained.

RNA tertiary structure mediation by adenosine platforms

The crystal structure of a group I intron domain reveals an unexpected motif that mediates both intra- and intermolecular interactions. At three separate locations in the 160-nucleotide domain, adjacent adenosines in the sequence lie side-by-side and form a pseudo-base pair within a helix. This adenosine platform opens the minor groove for base stacking or base pairing with nucleotides from a noncontiguous RNA strand. The platform motif has a distinctive chemical modification signature that may enable its detection in other structured RNAs. The ability of this motif to facilitate higher order folding provides one explanation for the abundance of adenosine residues in internal loops of many RNAs.

The Cbp2 protein suppresses splice site mutations in a group I intron

The Cbp2 protein facilitates the folding of a group I intron in the COB pre-mRNA of yeast mitochondria. Based on its ability to suppress mutations affecting the auto-catalytic reaction, the protein appears to play a role in the selection of splice sites. Adding Cbp2 did not overcome the effects of mutations in P1 whose primary effect was on the first step of splicing. In contrast, most mutations affecting the ligation of exons were suppressed in vitro by Cbp2. These included mutations in P1, P9.0 and P10. In fact, a mutant transcript lacking both P9.0 and P10 ligated efficiently in the presence of Cbp2. P9.0 and P10 mutations also reduced the rate of cleavage at the 5' splice junction, and this effect was only partially mitigated by adding Cbp2. A competitive secondary structure near the 3' splice junction blocked Cbp2-stimulated splicing, but this mutation could be suppressed by co-transcriptional splicing in the presence of Cbp2. Our data underscore the importance of the interaction between the 5' and 3' splice junctions in group I introns and suggest that nucleotide-nucleotide interactions that stabilize the structure of group I introns can be superceded by protein-RNA interactions.

A synthetic model for triple-helical domains in self-splicing group I introns studied by ultraviolet and circular dichroism spectroscopy

Structural studies were performed on synthetic oligonucleotides with sequences corresponding to the P4/P6 and J3/4, J6/7 regions of the self-splicing group I intron of the bacteriophage T4 nrdB pre-mRNA, which correspond to the proposed triple-helical domain in the Tetrahymena thermophila intron. A 23-mer RNA was synthesized as a mixed ribo-deoxyribo oligonucleotide, modeling an expected base-paired region of P4 along with the J3/4 and P6 (5'-end bases of P6) regions. strand modeling the 3'-end bases of P6 and J6/7 regions, with which a triple helix may form, was synthesized as a pure oligoribonucleotide (7-mer RNA). The interactions of these oligonucleotides have been characterized by UV and circular dichroism (CD) spectroscopy. The results show that the 23-mer RNA forms a stable hairpin modeling the P4 base-paired region. Triple helix association between the 23-mer RNA hairpin and the 7-mer RNA single strand was detected by CD in the presence of Mg2+ (>5mM) but not in presence of a monovalent cation like Na+ (up to 500 mM). Studies on selected variants of both 7-mer and 23-mer RNAs were carried out. The results show that for the association of the two partner strands not only the formation of P6 helix but also triplet interactions between two strands are required. The association of the two strands in general follow a pattern predicted by comparative sequence analysis. Parallel studies on pure oligoribonucleotides having base sequence corresponding to those of the oligoribonucleotides showed no evidence of association under similar conditions, which could indicate that the 2'-hydroxyl groups of the riboses might play an important role in hydrogen bonding to form the required nucleoside triples. Molecular modeling studies on the proposed "plaited triple helix" formed by the association of the 23-mer RNA hairpin and 7-mer RNA single strand showed that the structure is sterically and energetically feasible.

Conformation of an RNA molecule that models the P4/P6 junction for group I introns

We present a three-dimensional structure of a 34-nucleotide RNA molecule determined by NMR spectroscopy. The molecule was designed to form a junction between two double-helical stems whose sequence was based on the P4/P6 domain from group I introns. There are 5' and 3' single-strand overhangs at the junctions of the stems. Contrary to our expectations, we found that the 3' end of the molecule is placed in the minor and not the major groove of the P4 helix. As a result of tertiary contacts and stacking interactions from nucleotides in the 3' end, the junction helices are rotated in a left-handed fashion and do not stack coaxially. This conformation is highly dependent on the presence of single-stranded nucleotides at the 3' overhang. When the 3' end is removed, the molecule assumes a radically different structure with 5' end in the minor groove of the P6 helix and overall right-handed rotation between the stems. Only one nucleotide at the 3' end is sufficient to change the geometry of the junction.

Cotranscriptional splicing of a group I intron is facilitated by the Cbp2 protein

The nuclear CBP2 gene encodes a protein essential for the splicing of a mitochondrial group I intron in Saccharomyces cerevisiae. This intron (bI5) is spliced autocatalytically in the presence of high concentrations of magnesium and monovalent salt but requires the Cbp2 protein for splicing under physiological conditions. Addition of Cbp2 during RNA synthesis permitted cotranscriptional splicing. Splicing did not occur in the transcription buffer in the absence of synthesis. The Cbp2 protein appeared to modify the folding of the intron during RNA synthesis: pause sites for RNA polymerase were altered in the presence of the protein, and some mutant transcripts that did not splice after transcription did so during transcription in the presence of Cbp2. Cotranscriptional splicing also reduced hydrolysis at the 3' splice junction. These results suggest that Cbp2 modulates the sequential folding of the ribozyme during its synthesis. In addition, splicing during transcription led to an increase in RNA synthesis with both T7 RNA polymerase and mitochondrial RNA polymerase, implying a functional coupling between transcription and splicing.

Mutations in the mitochondrial split gene COXI are preferentially located in exons: a mapping study of 170 mutants

We have analysed the precise location of a large number (170) of mutations affecting the structural gene for subunit I of the cytochrome c oxidase complex. This gene, COXI, is 12.9 kb long and the major part of the sequence (i.e. 11.3 kb) is composed of introns. Several conclusions can be drawn from this study: (1) A significant proportion (84/170) of the mutations cannot be assigned to a single position within the gene by deletion mapping, in spite of clearly being located in it. These mutations are probably large deletions or multiple mutations. (2) Four mutants carry distant double mutations, which have been individually localized. (3) Eighty-two mutants have lesions that are restricted to very short regions of the gene and we therefore conclude that they are most probably due to single hits; amongst these single mutations, 41 are unambiguously located in exons and 28 in introns. This result implies that, at least in this particular split gene, the probability of selection of a mutant phenotype in an exon is, on the average, 13.3 times greater than in an intron, in spite of the existence, within most of these introns, of open reading frames specifying intronic proteins. The evolutionary significance and biological implications of these results are discussed.

Suppression of mutations in the core of the Tetrahymena ribozyme by spermidine, ethanol and by substrate stabilization

We have previously described a collection of mutations in conserved residues of the core of the Tetrahymena self-splicing intron. Most of these single base substitutions have less than 10% of the activity of their parental intron derivative [Couture, S., et al., (1990) J. Mol. Biol., 215, 345-358]. We examined the effect of two agents known to stabilize RNA structure, spermidine and ethanol, on the activity of many of these mutant RNAs. In the presence of either 5 mM spermidine or 20% ethanol most substitution mutations were partially or completely suppressed. These conditions also increased the temperature optima of both wild-type and mutant ribozymes. In addition, we find that mutations are also suppressed by a high concentration of GTP, a substrate in the reaction which is bound specifically by the intron. Thus we observe a general suppression of mutations in an RNA enzyme (ribozyme) by spermidine, ethanol and by substrate stabilization. These results are consistent with the idea that most mutations destabilize the folded structure of the ribozyme and can be suppressed by any of a variety of stabilizing influences.

Kinetic intermediates in RNA folding

The folding pathways of large, highly structured RNA molecules are largely unexplored. Insight into both the kinetics of folding and the presence of intermediates was provided in a study of the Mg(2+)-induced folding of the Tetrahymena ribozyme by hybridization of complementary oligodeoxynucleotide probes. This RNA folds via a complex mechanism involving both Mg(2+)-dependent and Mg(2+)-independent steps. A hierarchical model for the folding pathway is proposed in which formation of one helical domain (P4-P6) precedes that of a second helical domain (P3-P7). The overall rate-limiting step is formation of P3-P7, and takes place with an observed rate constant of 0.72 +/- 0.14 minute-1. The folding mechanism of large RNAs appears similar to that of many multidomain proteins in that formation of independently stable substructures precedes their association into the final conformation.

A non-directed, hydroxylamine-generated suppressor mutation in the P3 pairing region of the bacteriophage T4 td intron partially restores self-splicing capability

Hydroxylamine (HA) mutagenesis of an HA-induced splicing-defective bacteriophage T4 td intron mutant with a mutation in the intron P3 RNA pairing region was used to generate pseudorevertants. Because HA can only cause GC to AT transitions, the original mutant (H104A) could not undergo true reversion, yet the compensatory mutation on the opposite side of the P3 helix, which was complementary to the original H104A mutation, could occur. A pseudorevertant was isolated that contained both the original H104A mutation and the compensatory mutation HS9. By phenotypic and molecular genetic criteria, this double mutant (H104A-HS9) was shown to be able to undergo significant RNA splicing, thus confirming the existence and functional importance of the long-range P3 pairing region in this phage intron. The second-site suppressor mutation (HS9) was isolated by phage cross and also exhibited some self-splicing ability. A correlation exists between the strength of P3 helix Watson-Crick base pairing and the apparent level of splicing when wild-type, H104A, HS9, and H104A-HS9 are compared. This suggests that the primary role of the P3 RNA pairing region in the T4 td intron is structural in contributing to the critical RNA secondary structure.

Deletion of P9 and stem-loop structures downstream from the catalytic core affects both 5' and 3' splicing activities in a group-I intron

The P9 stem-loop is one of the conserved structural elements found in all group-I introns. Using two deletion mutants in this region of the Tetrahymena thermophilia large ribosomal subunit intron, we show that removal of the P9 element, either alone, or together with the non-conserved downstream P9.1 and P9.2 elements, results in an intron incapable of the first step of the splicing reaction at a low concentration of Mg2+. The mutant introns also require high concentrations of Mg2+ for the second step in splicing, as well as hydrolysis reactions, suggesting that P9, as well as P9.1 and P9.2, are important structural elements in the final folded form of the intron. In addition, RNase-T1-mediated-structure-probing experiments demonstrated that the loss of P9, P9.1 and P9.2 changes the structural context of the region binding the 5' splice site. The deletions lead to less efficient recognition of the 3' splice site and an accumulation of unligated exons. These observations support the view that the P9, P9.1 and P9.2 stem-loops play an important role in the binding of the 3' splice site.

Representation of the secondary and tertiary structure of group I introns

Group I introns, which are widespread in nature, carry out RNA self-splicing. The secondary structure common to these introns was for the most part established a decade ago. Information about their higher order structure has been derived from a range of experimental approaches, comparative sequence analysis, and molecular modelling. This information now provides the basis for a new two-dimensional structural diagram that more accurately represents the domain organization and orientation of helices within the intron, the coaxial stacking of certain helices, and the proximity of key nucleotides in three-dimensional space. It is hoped that this format will facilitate the detailed comparison of group I intron structures.

A three-dimensional model of the Rev-binding element of HIV-1 derived from analyses of aptamers

Coordinated variations in the sequence of the Rev-binding element of HIV-1, identified by in vitro genetic selections, have been used as distance and conformational constraints for molecular modelling. Three-dimensional models of the wild-type Rev-binding element and several, evolved RNA ligands (aptamers) have been constructed. These models demonstrate that non-Watson-Crick pairings open the major groove allowing access of an alpha-helical peptide from Rev, and explain why some selected RNA sequences can bind Rev more tightly than others.

A group-I self-splicing intron in the nuclear small subunit rRNA-encoding gene of the green alga, Chlorella ellipsoidea C-87

We report the presence of a 442-bp group-I self-splicing intron in the nuclear small subunit (SSU) rRNA-encoding gene (rDNA) of the unicellular green alga, Chlorella ellipsoidea C-87 (C. saccharophila 211-1a). The intron was found to be inserted at a position within the highly conserved helix 48 that was close to the 3' terminus of the SSU rRNA. The position was exactly the same as previously identified for the Pneumocystis carinii intron. A secondary structure model for the C. ellipsoidea intron contained all P1-P10 motifs of the group-I introns. Although the overall secondary structure of the C. ellipsoidea intron was substantially different from that of the intron in the nuclear large subunit rDNA of Tetrahymena thermophila, the nucleotide (nt) sequences constituting the catalytic core were strikingly conserved between the two; only three of 48 nt were different. The C. ellipsoidea intron was autocatalytically excised from the transcript in vitro via the group-I mechanism under somewhat unique conditions. No SSU rDNA intron was found in six other Chlorella species, including C. fusca var. vacuolata, C. kessleri, C. minutissima, C. protothecoides, C. sorokiniana and C. vulgaris.

In vitro genetic analysis of the hinge region between helical elements P5-P4-P6 and P7-P3-P8 in the sunY group I self-splicing intron

Modeling of the group I intron RNA suggests that its catalytic core is primarily composed of two extended structural elements (stacked helices P5-P4-P6 and P7-P3-P8) whose relative orientation is partially determined by base-triple interactions between paired regions P4 and P6, and single-stranded joining regions J6/7 and J3/4, respectively. In vitro genetic selection was used to isolate functional sequence variants of the proposed triple helical domain of the sunY intron. Comparative sequence analysis of the selected variants provided supporting evidence for the two previously established base-triples between P4 and J6/7 and provided the first experimental evidence for an interaction between P6(1) and J3/4(3). Sequence covariations also indicated that a simple relationship exists between the length of a single-stranded joining region, J3/4, and the identity of a particular base-pair, P4(1). Selected variants based on a core structure with an extra nucleotide inserted in J3/4 revealed two different responses to this structural perturbation: a base-triple interaction and an intrahelical bulged pyrimidine. Chemical modification analysis supported the existence of these alternative structures. The function of this region of the ribozyme can therefore be fulfilled by at least three different structures.

Selection and design of high-affinity RNA ligands for HIV-1 Rev

We have used in vitro selection to isolate minimal, high-affinity RNA ligands for the Rev protein of HIV-1. Sequence analysis reveals that the tightest binding aptamers exhibit some similarity to a Rev-binding element (RBE) localized within the Rev-responsive element (RRE), but also contain novel sequence and structural motifs. A short helical stem and bulged nucleotides (nt) CUC ... UYGAG that have no counterpart in the wild-type (wt) element contribute to high-affinity binding. We have designed and synthesized a short (37 nt) RNA molecule that incorporates this motif; this RNA ligand has from three- to fivefold tighter binding than the full-length wt element, and up to 16-fold tighter than minimal wt RBEs. A guanosine:guanosine pairing that is postulated to occur in the wt element has been altered to other base pairings in the context of our optimized minimal element. RNAs that contain non-Watson-Crick base pairings, that can be modeled as isosteric to the wt G:G pair, bind Rev up to 160-fold tighter than elements that contain canonical Watson-Crick pairings or non-isosteric mismatches. These results support the hypothesis that Rev recognizes structural features associated with a non-Watson-Crick base pair.

Selective optimization of the Rev-binding element of HIV-1

RNA molecules that can bind to the Rev protein of HIV-1 have been isolated from random sequence nucleic acid pools based on a minimal Rev-binding element (RBE) found within the Rev Responsive Element (RRE). While the selected sequences are related to the wild-type element, they also contain substitutions that allow them to bind Rev up to 10-fold better in vitro. A hypothesized homopurine pairing at G48:G71 is generally replaced by A48:A71; the occasional selection of C48:A71 suggests that R71 may be in a syn conformation. These data support the structural model for the RBE originally proposed by Bartel et al. (1). Additional interactions with the Rev protein are promoted by the sequence CUC ... UYGAG, found in one class of high-affinity aptamers, but absent from the wild-type element. Within each class of aptamers different residues and substructures covary with one another to generate optimal Rev-binding surfaces. The interdependencies of different nucleotide substitutions suggest structural models for both the wild-type RBE and the selected high-affinity aptamers.

Movement of the guide sequence during RNA catalysis by a group I ribozyme

Ribozymes derived from the self-splicing pre-ribosomal RNA of Tetrahymena act as sequence-specific endonucleases. The reaction involves binding an RNA or DNA substrate by base pairing to the internal guide sequence (IGS) to form helix P1. Site-specific photo-crosslinking localized the 5' end of the IGS in helix P1 to the vicinity of conserved bases between helices P4 and P5, supporting a major feature of the Michel-Westhof three-dimensional structure model. The crosslinked ribozyme retained catalytic activity. When not base-paired, the IGS was still specifically crosslinked, but the major site was 37 A distant from the reactive site in the experimentally supported three-dimensional model. The data indicate that a substantial induced-fit conformational change accompanies P1 formation, and they provide a physical basis for understanding the transport of oligonucleotides to the catalytic core of the ribozyme. The ability of RNA to orchestrate large-scale conformational changes may help explain why the ribosome and the spliceosome are RNA-based machines.

Evolution in vitro of an RNA enzyme with altered metal dependence

The Tetrahymena group I ribozyme catalyses a sequence-specific phosphodiester cleavage reaction on an external RNA oligonucleotide substrate in the presence of a divalent metal cation cofactor. This reaction proceeds readily with either Mg2+ or Mn2+, but no detectable reaction has been reported when other divalent cations are used as the sole cofactor. Cations such as Ca2+, Sr2+ and Ba2+ can stabilize the correct folded conformation of the ribozyme, thereby partially alleviating the Mg2+ or Mn2+ requirement. But catalysis by the ribozyme involves coordination of either Mg2+ or Mn2+ at the active site, resulting in an overall requirement for one of these two cations. Here we use an in vitro evolution process to obtain variants of the Tetrahymena ribozyme that are capable of cleaving an RNA substrate in reaction mixtures containing Ca2+ as the divalent cation. These findings extend the range of different chemical environments available to RNA enzymes and illustrate the power of in vitro evolution in generating macromolecular catalysts with desired properties.

A base-triple structural domain in RNA

An oligonucleotide modeled on a proposed base-triple domain of the Tetrahymena group I intron has been characterized by NMR. The oligonucleotide contains two double-helix regions with adjacent single-stranded nucleotides. The NMR data show that the two helices stack coaxially, although the rotation between the two helices is approximately twice as large as the rotation between normal base pairs. The rotation between the two helices allows the single-stranded nucleotides to form U.U.G and A.G.C base triples in the minor groove. The A.G.C base triple contains a hydrogen bond between the adenine N1 and a 2'-hydroxyl in the minor groove of the G.C pair. A similar hydrogen bond between an adenine and a 2'-hydroxyl in transfer RNA suggests that this could be a recurring tertiary interaction in RNA.

Mutational analysis of the yeast U2 snRNA suggests a structural similarity to the catalytic core of group I introns

We have used an in vitro reconstitution system to determine the effects of a large number of mutations in the highly conserved 5' terminal domain of the yeast U2 snRNA on pre-mRNA splicing. Whereas many mutations have little or no functional consequence, base substitutions in two regions were found to have drastic effects on pre-mRNA splicing. A previously unrecognized function for the U2 snRNA in the second step of splicing was found by alteration of the absolutely conserved sequence AGA upstream of the branch point recognition sequence. The effects of these mutations suggest the formation of a structure involving the U2 snRNA similar to the guanosine-binding site found in the catalytic core of group I introns.

A region of group I introns that contains universally conserved residues but is not essential for self-splicing

The catalytic core of the self-splicing group I intron RNAs is composed of six paired regions together with their connecting sequences; these are thought to form two elongated domains, with paired regions P5, P4, and P6 aligned along one axis and P8, P3, and P7 along the other. Most of the very highly conserved residues of the group I introns lie in or near P7, but two occur in L4, the internal loop connecting P4 and P5. It is generally believed that such bases are conserved because they are essential for splicing. Mutants were created in a member of each of the two major subclasses of group I introns, in which P5, L4, and the distal portion of P4 were deleted. Splicing activity was still detected in these mutants, albeit substantially weakened; splicing was accurate and occurred by the normal group I mechanism, with addition of a guanosine molecule to the intron. Thus the deleted region, containing two universally conserved bases, is not essential but facilitates splicing. Another reaction characteristic of group I introns, hydrolysis of the 3' splice site, was less severely affected by the deletions. The results are discussed in terms of the prevailing three-dimensional model for the core structure of the group I introns.

Structure and evolution of myxomycete nuclear group I introns: a model for horizontal transfer by intron homing

We have examined five nuclear group I introns, located at three different positions in the large subunit ribosomal RNA (LSU rRNA) gene of the two myxomycete species, Didymium iridis and Physarum polycephalum. Structural models of intron RNAs, including secondary and tertiary interactions, are proposed. This analysis revealed that the Physarum intron 2 contains an unusual core region that lacks the P8 segment, as well as several of the base-triples known to be conserved among group I introns. Structural and evolutionary comparisons suggest that the corresponding introns 1 and 2 were present in a common ancestor of Didymium and Physarum, and that the five introns in LSU rRNA genes of these myxomycetes were acquired in three different events. Evolutionary relationships, inferred from the sequence analysis of several different nuclear group I introns and the ribosomal RNA genes of the intron-harbouring organisms, strongly support horizontal transfer of introns in the course of evolution. We propose a model that may explain how myxomycetes in natural environments obtained their nuclear group I introns.

Directed evolution of an RNA enzyme

An in vitro evolution procedure was used to obtain RNA enzymes with a particular catalytic function. A population of 10(13) variants of the Tetrahymena ribozyme, a group I ribozyme that catalyzes sequence-specific cleavage of RNA via a phosphoester transfer mechanism, was generated. This enzyme has a limited ability to cleave DNA under conditions of high temperature or high MgCl2 concentration, or both. A selection constraint was imposed on the population of ribozyme variants such that only those individuals that carried out DNA cleavage under physiologic conditions were amplified to produce "progeny" ribozymes. Mutations were introduced during amplification to maintain heterogeneity in the population. This process was repeated for ten successive generations, resulting in enhanced (100 times) DNA cleavage activity.

Mutational evidence for competition between the P1 and the P10 helices of a mitochondrial group I intron

A guanosine to cytosine transversion at position 2 of the fifth intron of the mitochondrial gene COB blocks the ligation step of splicing. This mutation prevents the formation of a base pair within the P1 helix of this group I intron--the RNA duplex formed between the 3' end of the upstream exon and the internal guide sequence. The mutation also reduces the rate of the first step of splicing (guanosine addition at the 5' splice junction) while stimulating hydrolysis at the 3' intron-exon boundary. Consequently, the ligation of exons is blocked because the 3' exon is removed prior to cleavage at the 5' splice junction. The lesion can be suppressed by second-site mutations that preserve the potential for base-pairing at this position. Because the P1 duplex and the P10 duplex (between the guide sequence and the 3' exon) overlap at the affected pairings represent alternative structures that do not, indeed cannot, form simultaneously.

Tertiary structure around the guanosine-binding site of the Tetrahymena ribozyme

A cleavage reagent directed to the active site of the Tetrahymena catalytic RNA was synthesized by derivatization of the guanosine substrate with a metal chelator. When complexed with iron(II), this reagent cleaved the RNA in five regions. Cleavage at adenosine 207, which is far from the guanosine-binding site in the primary and secondary structure, provides a constraint for the higher order folding of the RNA. This cleavage site constitutes physical evidence for a key feature of the Michel-Westhof model. Targeting a reactive entity to a specific site should be generally useful for determining proximity within folded RNA molecules or ribonucleoprotein complexes.

Circle reopening in the Tetrahymena ribozyme resembles site-specific hydrolysis at the 3' splice site

The Tetrahymena intron, after splicing from its flanking exons, can mediate its own circularization. This is followed by site-specific hydrolysis of the phosphodiester bond formed during the circularization reaction. The structural components involved in recognition of this bond for hydrolysis have not been established. We have made base substitutions to the P9.0 pairing and at the 3'-terminal guanosine residue (G414) of the intron to investigate their effects on circle formation and reopening. We have found that disruption of either P9.0 pairing or binding of the terminal nucleotide result in the formation of a large circle, C-413:5E23 from precursor RNA molecules that have undergone hydrolysis at the 3' splice site. This circle is formed at the phosphodiester bond of the 5'-terminal guanosine residue of the upstream exon via nucleophilic attack by the 3'-terminal nucleotide of the intron. The large circle is novel since it can reopen eight bases downstream from the original circularization junction at a site resembling the normal 3' splice site, restoring a guanosine to the 3' terminus and re-establishing P9.0 pairing. The new 3' terminus of the intron is capable of recircularization at any of the three normal wild-type sites. We conclude that both P9.0 and the 3'-terminal guanosine residue are required for the selection of the phosphodiester bond hydrolysed during circle reopening.

An axial binding site in the Tetrahymena precursor RNA

Previous studies allow the construction of three distinct models of the binding of G and arginine within the active site of the Tetrahymena self-splicing preribosomal precursor RNA. These models (base triple, axial I and axial II) are now distinguished by measurements on the specificity of RNAs with nucleotide substitutions at positions spanning the site. Because the semi-conserved unpaired nucleotide 263 has no effect on substrate or inhibitor selection by the Tetrahymena RNA we conclude that the axial I model is improbable. In contrast, data with substituted RNAs and nucleoside analogs suggest that nucleotide 265 makes a hydrogen bond with the substrate. Accordingly the active site appears axial because substrate contacts exist at more than one nucleotide on the 5' side of the P7 helix. The effects of this hydrogen bond are observable in cases where the donor or acceptor is on the RNA, whether nucleotide 265 is a purine or pyrimidine, or whether nucleotide 265 is mispaired, wobble paired or normally paired. This pattern is consistent with the axial II model. Molecular dynamics and energy minimization calculations lead to the same conclusions as these site-directed substitutions; the base triple and axial I models are unstable dynamically. Under thermal agitation, the third model site (axial II) is transformed to a related, but more stable structure, axial III. The axial III active site is characterized by the extrusion of the conserved bulged base 263 from the P7 helix, a semi-pocket for G base formed by stacking of nucleotide 262, and formation of all bonds to the G base originally proposed for both the base triple and axial II sites. Because of these hydrogen bonds the axial III site is also consistent with data on enzymatic specificity. The axial III model indicates an unforeseen capacity for pocket formation within the groove of an RNA helix, suggests that the site may be unusually flexible, and bears on a hypothesis concerning the origin of the genetic code.

Mutations in a nonconserved sequence of the Tetrahymena ribozyme increase activity and specificity

The RNA substrate-binding site of the Tetrahymena ribozyme is connected to the catalytic core by the joining region J1/2. Although J1/2 is not conserved among group I introns, small insertions or deletions in this sequence have dramatic effects, enhancing the turnover number and sequence specificity of ribozyme-catalyzed RNA cleavage. Measurements of rate constants for individual steps in the reaction have revealed the basis of these improvements. Ironically, the higher turnover and specificity both result from decreased affinity for RNA, rather than better cleavage. These results provide evidence that the nonconserved J1/2 sequence positions the RNA substrate to optimize tertiary interactions and ensure cleavage at the position corresponding to the 5' splice site. The wild-type RNA is well adapted to its biological function, and its limitations in multiple turnover can be corrected by mutation.

A self-splicing group I intron in the nuclear pre-rRNA of the green alga, Ankistrodesmus stipitatus

The nuclear small subunit ribosomal RNA gene of the unicellular green alga Ankistrodesmus stipitatus contains a group I intron, the first of its kind to be found in the nucleus of a member of the plant kingdom. The intron RNA closely resembles the group I intron found in the large subunit rRNA precursor of Tetrahymena thermophila, differing by only eight nucleotides of 48 in the catalytic core and having the same peripheral secondary structure elements. The Ankistrodesmus RNA self-splices in vitro, yielding the typical group I intron splicing intermediates and products. Unlike the Tetrahymena intron, however, splicing is accelerated by high concentrations of monovalent cations and is rate-limited by the exon ligation step. This system provides an opportunity to understand how limited changes in intron sequence and structure alter the properties of an RNA catalytic center.

Synthesis of RNA containing inosine: analysis of the sequence requirements for the 5' splice site of the Tetrahymena group I intron

Two protected derivatives of the ribonucleoside inosine have been prepared to serve as building blocks for phosphoramidite-based synthesis of RNA. Two different synthetic routes address the unusual solubility characteristics of inosine and its derivatives. The final products of the different synthetic pathways, 5'-O-(dimethoxytrityl)-2'-O-(t-butyldimethylsiyl) inosine 3'-O-(beta-cyanoethyldiisopropylamino) phosphoramidite 5a, and O6-p-nitrophenylethyl-5'-O-(dimethoxytrityl)-2'-O-(t-butyldimethylsilyl) inosine 3'-O-(methyldiisopropylamino) phosphoramidite 5b, were chemically incorporated into short oligoribonucleotides which also contained the four standard ribonucleoside bases. The oligomers were chosen to study base-specific interactions between an RNA substrate and an RNA enzyme derived from the Group I Tetrahymena self-splicing intron. The oligomers were shown to be biochemically competent using a trans cleavage assay with the modified Tetrahymena intron. The results confirm the dependence of the catalytic activity on a wobble base pair, rather than a Watson-Crick base pair, in the helix at the 5'-splice site. Furthermore, comparison of guanosine and inosine in a wobble base pair allows one to assess the importance of the guanine 2-amino group for biological activity. The preparation of the inosine phosphoramidites adds to the repertoire of base analogues available for the study of RNA catalysis and RNA-protein interactions.