RNA tectonics: towards RNA design

TectoRNA: modular assembly units for the construction of RNA nano-objects

Structural information on complex biological RNA molecules can be exploited to design tectoRNAs or artificial modular RNA units that can self-assemble through tertiary interactions thereby forming nanoscale RNA objects. The selective interactions of hairpin tetraloops with their receptors can be used to mediate tectoRNA assembly. Here we report on the modulation of the specificity and the strength of tectoRNA assembly (in the nanomolar to micromolar range) by variation of the length of the RNA subunits, the nature of their interacting motifs and the degree of flexibility of linker regions incorporated into the molecules. The association is also dependent on the concentration of magnesium. Monitoring of tectoRNA assembly by lead(II) cleavage protection indicates that some degree of structural flexibility is required for optimal binding. With tectoRNAs one can compare the binding affinities of different tertiary motifs and quantify the strength of individual interactions. Furthermore, in analogy to the synthons used in organic chemistry to synthesize more complex organic compounds, tectoRNAs form the basic assembly units for constructing complex RNA structures on the nanometer scale. Thus, tectoRNA provides a means for constructing molecular scaffoldings that organize functional modules in three-dimensional space for a wide range of applications.

Plasticity, evolvability, and modularity in RNA

RNA folding from sequences into secondary structures is a simple yet powerful, biophysically grounded model of a genotype-phenotype map in which concepts like plasticity, evolvability, epistasis, and modularity can not only be precisely defined and statistically measured but also reveal simultaneous and profoundly non-independent effects of natural selection. Molecular plasticity is viewed here as the capacity of an RNA sequence to assume a variety of energetically favorable shapes by equilibrating among them at constant temperature. Through simulations based on experimental designs, we study the dynamics of a population of RNA molecules that evolve toward a predefined target shape in a constant environment. Each shape in the plastic repertoire of a sequence contributes to the overall fitness of the sequence in proportion to the time the sequence spends in that shape. Plasticity is costly, since the more shapes a sequence can assume, the less time it spends in any one of them. Unsurprisingly, selection leads to a reduction of plasticity (environmental canalization). The most striking observation, however, is the simultaneous slow-down and eventual halting of the evolutionary process. The reduction of plasticity entails genetic canalization, that is, a dramatic loss of variability (and hence a loss of evolvability) to the point of lock-in. The causal bridge between environmental canalization and genetic canalization is provided by a correlation between the set of shapes in the plastic repertoire of a sequence and the set of dominant (minimum free energy) shapes in its genetic neighborhood. This statistical property of the RNA genotype-phenotype map, which we call plastogenetic congruence, traps populations in regions where most genetic variation is phenotypically neutral. We call this phenomenon neutral confinement. Analytical models of neutral confinement, made tractable by the assumption of perfect plastogenetic congruence, formally connect mutation rate, the topography of phenotype space, and evolvability. These models identify three mutational regimes: that corresponding to neutral confinement, an exploration threshold corresponding to a breakdown of neutral confinement with the simultaneous persistence of the dominant phenotype, and a classic error threshold corresponding to the loss of the dominant phenotype. In a final step, we analyze the structural properties of canalized phenotypes. The reduction of plasticity leads to extreme modularity, which we analyze from several perspectives: thermophysical (melting--the RNA version of a norm of reaction), kinetic (folding pathways--the RNA version of development), and genetic (transposability--the insensitivity to genetic context). The model thereby suggests a possible evolutionary origin of modularity as a side effect of environmental canalization.

New developments in structure determination of pseudoknots

Recently, several high-resolution structures of-RNA pseudoknots have become available. Here we review the progress in this area. The majority of the structures obtained belong to the classical or H-type pseudoknot family. The most complicated pseudoknot structure elucidated so far is the Hepatitis Delta Virus ribozyme, which forms a nested double pseudoknot. In particular, the structure-function relationships of the H-type pseudoknots involved in translational frameshifting have received much attention. All molecules considered show interesting new structural motifs.

Measuring dissociation constants of RNA and aminoglycoside antibiotics by electrospray ionization mass spectrometry

Electrospray ionization mass spectrometry (ESI-MS) has been used to determine the dissociation constants (K(D)s) and binding stoichiometry for tobramycin and paromomycin with a 27-nucleotide RNA construct representing the A-site of the 16S ribosomal RNA. K(D) values determined by holding the ligand concentration fixed are compared with K(D) values derived by holding the RNA target concentration fixed. Additionally, the effect of solution conditions such as the amount of organic solvent present and the amount of salt present in the solution on the K(D) measurement is investigated. It is shown that the preferred method for determining dissociation constants using ESI-MS is holding the RNA target concentration fixed below the expected K(D) and titrating the ligand. K(D) measurements should also be carried out at as high as possible salt concentration to minimize nonspecific binding due primarily to electrostatic interactions. For tobramycin, two nonequivalent binding sites were found with K(D1) = 352 nM and K(D2) = 9 microM. For paromomycin, there is only one binding site with K(D) = 52 nM.

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.

A complex ligase ribozyme evolved in vitro from a group I ribozyme domain

Like most proteins, complex RNA molecules often are modular objects made up of distinct structural and functional domains. The component domains of a protein can associate in alternative combinations to form molecules with different functions. These observations raise the possibility that complex RNAs also can be assembled from preexisting structural and functional domains. To test this hypothesis, an in vitro evolution procedure was used to isolate a previously undescribed class of complex ligase ribozymes, starting from a pool of 10(16) different RNA molecules that contained a constant region derived from a large structural domain that occurs within self-splicing group I ribozymes. Attached to this constant region were three hypervariable regions, totaling 85 nucleotides, that gave rise to the catalytic motif within the evolved catalysts. The ligase ribozymes catalyze formation of a 3',5'-phosphodiester linkage between adjacent template-bound oligonucleotides, one bearing a 3' hydroxyl and the other a 5' triphosphate. Ligation occurs in the context of a Watson-Crick duplex, with a catalytic rate of 0.26 min(-1) under optimal conditions. The constant region is essential for catalytic activity and appears to retain the tertiary structure of the group I ribozyme. This work demonstrates that complex RNA molecules, like their protein counterparts, can share common structural domains while exhibiting distinct catalytic functions.

The effect of long-range loop-loop interactions on folding of the Tetrahymena self-splicing RNA

Bas?e-pairing between the terminal loops of helices P2.1 and P9.1a (P13) and P2 and P5c (P14) stabilize the folded structure of the Tetrahymena group I intron. Using native gel electrophoresis to analyze the folding kinetics of a natural pre-RNA containing the Tetrahymena intron, we show that P13 and P14 are the only native loop-loop interactions among six possible combinations. Other base-pairing interactions of the loop sequences stabilize misfolded and inactive pre-RNAs. Mismatches in P13 or P14 raised the midpoints and decreased the cooperativity of the Mg(2+)-dependent eqXuilibrium folding transitions. Although some mutations in P13 resulted in slightly higher folding rates, others led to slower folding compared to the wild-type, suggesting that P13 promotes formation of P3 and P7. In contrast, mismatches in P14 increased the rate of folding, suggesting that base-pairing between P5c and P2 stabilizes intermediates in which the catalytic core is misfolded. Although the peripheral helices stabilize the native structure of the catalytic core, our results show that formation of long-range interactions, and competition between correct and incorrect loop-loop base-pairs, decrease the rate at which the active pre-RNA structure is assembled.

The distribution of RNA motifs in natural sequences

Functional analysis of genome sequences has largely ignored RNA genes and their structures. We introduce here the notion of 'ribonomics' to describe the search for the distribution of and eventually the determination of the physiological roles of these RNA structures found in the sequence databases. The utility of this approach is illustrated here by the identification in the GenBank database of RNA motifs having known binding or chemical activity. The frequency of these motifs indicates that most have originated from evolutionary drift and are selectively neutral. On the other hand, their distribution among species and their location within genes suggest that the destiny of these motifs may be more elaborate. For example, the hammerhead motif has a skewed organismal presence, is phylogenetically stable and recent work on a schistosome version confirms its in vivo biological activity. The under-representation of the valine-binding motif and the Rev-binding element in GenBank hints at a detrimental effect on cell growth or viability. Data on the presence and the location of these motifs may provide critical guidance in the design of experiments directed towards the understanding and the manipulation of RNA complexes and activities in vivo.

Sequential folding of the genomic ribozyme of the hepatitis delta virus: structural analysis of RNA transcription intermediates

The structures of the model oligoribonucleotides that mimic the consecutive stages in the transcription of genomic HDV ribozyme have been analyzed by the Pb(2+)-induced cleavage method, partial digestion with specific nucleases and chemical probing. In the transcription intermediates, the P1 and P4 helical segments are found to be present in the final folded forms in which they exist in the full-length transcript. However, the region corresponding to the central hairpin forms another thermodynamically stable hairpin structure. Its correct folding requires the presence of a ribozyme 3'-terminal sequence and the formation of helix P2. This confirms the ribozyme structure of the pseudoknot type and points to the crucial role of helix P2 in its overall folding. Moreover, we show that the J4/2 region can be specifically cleaved in the presence of selected divalent metal ions in the full-length transcript, but not in a shorter one lacking six 3'-terminal nucleotides, which cannot form the pseudoknotted structure. Thus, a particular RNA conformation around that cleavage site is required for specific hydrolysis, and the J4/2 region seems to be involved in the formation of a general metal ion binding site. Recently, it has been proposed that, in the antigenomic ribozyme, a four nucleotide sequence within the J1/2 region may contribute to the folding pathway, being part of a mechanism responsible for controlling ribozyme cleavage activity. Our study shows that in the genomic ribozyme the central hairpin region may contribute to a similar mechanism, providing a barrier to the formation of an active structure in the ribozyme folding pathway.

Design of allosteric hammerhead ribozymes activated by ligand-induced structure stabilization

BACKGROUND

Ribozymes can function as allosteric enzymes that undergo a conformational change upon ligand binding to a site other than the active site. Although allosteric ribozymes are not known to exist in nature, nucleic acids appear to be well suited to display such advanced forms of kinetic control. Current research explores the mechanisms of allosteric ribozymes as well as the strategies and methods that can be used to create new controllable enzymes.

RESULTS

In this study, we exploit the modular nature of certain functional RNAs to engineer allosteric ribozymes that are activated by flavin mononucleotide (FMN) or theophylline. By joining an FMN- or theophylline-binding domain to a hammerhead ribozyme by different stem II elements, we have identified a minimal connective bridge comprised of a G.U wobble pair that is responsive to ligand binding. Binding of FMN or theophylline to its allosteric site induces a conformational change in the RNA that stabilizes the wobble pair and ultimately favors the active form of the catalytic core. These ligand-sensitive ribozymes exhibit rate enhancements of more than 100-fold in the presence of FMN and of approximately 40-fold in the presence of theophylline.

CONCLUSIONS

An adaptive strategy for modular rational design has proven to be an effective approach to the engineering of novel allosteric ribozymes. This strategy was used to create allosteric ribozymes that function by a mechanism involving ligand-induced structure stabilization. Conceivably, similar engineering strategies and allosteric mechanisms could be used to create a variety of novel allosteric ribozymes that function with other effector molecules.

Oligonucleotide aptamers that recognize small molecules

Nucleic acid receptors ('aptamers'), which recognize a large variety of organic molecules of low molecular weight, have been isolated from combinatorial nucleic acid libraries by in vitro selection methods. Structural studies of nucleic acid-small molecule complexes provide insight into both the principles of molecular recognition by this class of biopolymers and the architecture of tertiary motifs in nucleic acid folding. Aptamers that recognize small molecules are increasingly applied as tools in molecular biology, from the detection of oxidative damage in DNA to conditional gene expression and from their use as modules for the engineering of allosteric ribozymes to biosensors.

Engineering precision RNA molecular switches

Ligand-specific molecular switches composed of RNA were created by coupling preexisting catalytic and receptor domains via structural bridges. Binding of ligand to the receptor triggers a conformational change within the bridge, and this structural reorganization dictates the activity of the adjoining ribozyme. The modular nature of these tripartite constructs makes possible the rapid construction of precision RNA molecular switches that trigger only in the presence of their corresponding ligand. By using similar enzyme engineering strategies, new RNA switches can be made to operate as designer molecular sensors or as a new class of genetic control elements.

Stepping through an RNA structure: A novel approach to conformational analysis

Drawing from the growing database of complex three-dimensional RNA structures, a systematic method has been developed for classifying and analyzing the variety of conformations adopted by nucleic acids. This method is based on the development of a reduced representation for nucleic acid backbone conformation, simplifying the formidable eight-dimensional problem that has long complicated nucleic acid conformational analysis. Two pseudotorsion angles (eta and theta) have been defined, based on the selection of two appropriate pivot points along the RNA backbone, P and C4'. These pseudotorsions, together with a complete library of conventional torsion angles, can be calculated for any RNA structure or all-atom model using a new program called AMIGOS. Having computed eta and theta pseudotorsions for each position on an RNA molecule, they can be represented on a two-dimensional plot similar to the phi-phi plots that have traditionally been used for protein conformational analysis. Like a Ramachandran plot, clusters of residues appear at discrete regions on an eta-theta plot. Nucleotides within these clusters share conformational properties, often belonging to the same type of structural motif such as A-platforms, sheared tandem purine-purine pairs and GNRA tetraloops. An eta-theta plot provides a two-dimensional representation of the conformational properties of an entire RNA molecule, facilitating rapid analysis of structural features. In addition to the utility of eta-theta plots for intuitive visualization of conformational space, the pseudotorsional convention described here should significantly simplify approaches to macromolecular modeling of RNA structure.

Defining functional groups, core structural features and inter-domain tertiary contacts essential for group II intron self-splicing: a NAIM analysis

Group II introns are self-splicing RNA molecules that are of considerable interest as ribozymes, mobile genetic elements and examples of folded RNA. Although these introns are among the most common ribozymes, little is known about the chemical and structural determinants for their reactivity. By using nucleotide analog interference mapping (NAIM), it has been possible to identify the nucleotide functional groups (Rp phosphoryls, 2'-hydroxyls, guanosine exocyclic amines, adenosine N7 and N6) that are most important for composing the catalytic core of the intron. The majority of interference effects occur in clusters located within the two catalytically essential Domains 1 and 5 (D1 and D5). Collectively, the NAIM results indicate that key tetraloop-receptor interactions display a specific chemical signature, that the epsilon-epsilon' interaction includes an elaborate array of additional features and that one of the most important core structures is an uncharacterized three-way junction in D1. By combining NAIM with site-directed mutagenesis, a new tertiary interaction, kappa-kappa', was identified between this region and the most catalytically important section of D5, adjacent to the AGC triad in stem 1. Together with the known zeta-zeta' interaction, kappa-kappa' anchors D5 firmly into the D1 scaffold, thereby presenting chemically essential D5 functionalities for participation in catalysis.

MANIP: an interactive tool for modelling RNA

Large RNA structures can be viewed as assemblies of smaller units or modules that are usually clearly identified (helices, hairpin loops, other recurrent motifs, etc.). We have developed a program, MANIP, which allows the rapid assembly of separate motifs (each with a specified sequence) into a complex three-dimensional architecture. The already determined modules are present in a database from which they can be extracted with the appropriate sequence. Their assembly is performed in real time on the computer screen with buttons and dials that command rotation and translation of any chosen fragment with respect to the chosen pivot, or that generate all possible variations of any torsion angle within a specified segment either in the 5' or in the 3' direction. The possible in-built manipulations follow the general stereochemical rules of RNA structure. MANIP automatically recognizes and displays the allowed and nonallowed hydrogen bonds between the residues. The program is interfaced with a rapid and automatic online refinement tool of partial or full assemblies, NUCLIN-NUCLSQ. The refinement protocol incorporates canonical as well as noncanonical base pairing constraints together with restraints imposed by covalent geometry, stereochemistry, and van der Waals contacts. The computer package runs on UNIX Silicon Graphics workstations and is written in C with OpenGL and X11/Motif libraries.

Derivation of the three-dimensional architecture of bacterial ribonuclease P RNAs from comparative sequence analysis

The secondary structure of bacterial RNase P RNA, a ribozyme responsible for the maturation of the 5' end of tRNAs, is well established on the basis of sequence comparison analysis. RNase P RNA secondary structures fall into two types, A and B, which share a common core formed by the assembly of two main folding domains, but differ in their peripheral elements.A revised alignment of 137 available sequences reveals new covariations allowing for the refinement of both types of secondary structures. Phylogenetic evidence is thus provided for the extension of stems P11, P14, P19, P10.1 and P15.1 through further canonical base-pairs or GAellipsisGA mismatches. These refinements led in turn to a new organization of the catalytic core, with coaxial stackings of helices P2 and P19 as well as P1 and P4. New inter-domain tertiary interactions involve loop L9 and helix P1 and loop L8 with helix P4. These features were incorporated into atomic-scale 3D models of RNase P RNA for representatives of each structural type, namely Escherichia coli and Bacillus subtilis. In each model, the juxtaposition of the core helices creates a cradle onto which the pre-tRNA substrate binds with most evolutionarily conserved residues converging towards the cleavage site. The inner cores of both types are stabilized similarly, albeit by different peripheral elements, emphasizing the modular and hierarchical organisation of the architecture of RNase P RNAs. Similarities are thus apparent between the type A modules, P16/P17/P6 and P13/P14, and their type B analogs, P5.1/P15.1 and P10. 1/P10.1a, respectively. Other noteworthy features of these models include compactness and good agreement with published crosslinking data.

Emerging themes in RNA folding

RNAs, like proteins, readily form specific structures adapted for ligand binding and catalysis. Since they are composed of completely different chemical building blocks, however, RNAs and proteins necessarily use distinct strategies to assemble complex architectures. While burial of hydrophobic residues drives protein folding, the hydrophobic effect in RNA contributes primarily to the formation of secondary structure. To form tertiary structure, RNA must overcome electrostatic repulsions from the phosphate backbone. How do negatively charged double helices pack together to produce catalytic centers and ligand binding surfaces? Here, we review our understanding of the principles that underlie RNA folding based on the structural information currently available.

ESSA: an integrated and interactive computer tool for analysing RNA secondary structure

With ESSA, we propose an approach of RNA secondary structure analysis based on extensive viewing within a friendly graphical interface. This computer program is organized around the display of folding models produced by two complementary methods suitable to draw long RNA molecules. Any feature of interest can be managed directly on the display and highlighted by a rich combination of colours and symbols with emphasis given to structural probe accessibilities. ESSA also includes a word searching procedure allowing easy visual identification of structural features even complex and degenerated. Analysis functions make it possible to calculate the thermodynamic stability of any part of a folding using several models and compare homologous aligned RNA both in primary and secondary structure. The predictive capacities of ESSA which brings together the experimental, thermodynamic and comparative methods, are increased by coupling it with a program dedicated to RNA folding prediction based on constraints management and propagation. The potentialities of ESSA are illustrated by the identification of a possible tertiary motif in the LSU rRNA and the visualization of a pseudoknot in S15 mRNA.

RNA as a drug target

Historically, the pharmaceutical industry has focused on proteins, rather than nucleic acids, as drug targets. But recent advances in the fields of RNA synthesis, structure determination and therapeutic target identification make the systematic exploitation of RNA as a drug target a realistic goal.

Protein-facilitated RNA folding

In the absence of protein collaborators, both simple and complex RNAs often misfold or are unfolded. Biologically important RNAs solve their folding problem, in part, using the assistance of chaperone and cofactor proteins. Recent work emphasizes several rules for RNA-protein complexes: formation involves induced fit; many large RNAs fold slowly; and ribonucleoprotein assembly requires multiple steps. Finally, protein binding can introduce thermodynamic side effects.

Rational design of allosteric ribozymes

BACKGROUND

Efficient operation of cellular processes relies on the strict control that each cell exerts over its metabolic pathways. Some protein enzymes are subject to allosteric regulation, in which binding sites located apart from the enzyme's active site can specifically recognize effector molecules and alter the catalytic rate of the enzyme via conformational changes. Although RNA also performs chemical reactions, no ribozymes are known to operate as true allosteric enzymes in biological systems. It has recently been established that small-molecule receptors can readily be made of RNA, as demonstrated by the in vitro selection of various RNA aptamers that can specifically bind corresponding ligand molecules. We set out to examine whether the catalytic activity of an existing ribozyme could be brought under the control of an effector molecule by designing conjoined aptamer-ribozyme complexes.

RESULTS

By joining an ATP-binding RNA to a self-cleaving ribozyme, we have created the first example of an allosteric ribozyme that has a catalytic rate that can be controlled by ATP. A 180-fold reduction in rate is observed upon addition of either adenosine or ATP, but no inhibition is detected in the presence of dATP or other nucleoside triphosphates. Mutations in the aptamer domain that are expected to eliminate ATP binding or that increase the distance between aptamer and ribozyme domains result in a loss of ATP-specific allosteric control. Using a similar design approach, allosteric hammerhead ribozymes that are activated in the presence of ATP were created and another ribozyme that can be controlled by theophylline was created.

CONCLUSIONS

The catalytic features of these conjoined aptamer-ribozyme constructs demonstrate that catalytic RNAs can also be subject to allosteric regulation-a key feature of certain protein enzymes. Moreover, by using simple rational design strategies, it is now possible to engineer new catalytic polynucleotides which have rates that can be tightly and specifically controlled by small effector molecules.

Context dependent RNA-RNA recognition in a three-dimensional model of the 16S rRNA core

A 3-D model of the core of the 16S rRNA of Escherichia coli containing 328 residues has been built in the protein map derived from neutron scattering data with the help of all the available phylogenetic, biochemical, and cross-linking data. The three pseudoknots of the 16S-core cluster, through the arrangement of complex three-, four- and five-way junctions, around the neck and at the subunit interface. The roles in assembly, initiation or elongation of the three pseudoknots in ribosomal dynamics are emphasized. The 530-loop, localized on the periphery of the 30S particle, could be built with and without a pseudoknot independently of the state of the particle. The pseudoknot of the central domain controls the dynamics of an helix connected to the subunit interface which could trigger some mechanism during translation. The process of the model construction is compatible with a folding scenario in which the 5'-terminal pseudoknot controls the assembly of the central junction and the subsequent folding of the 3'-major domain. The modelling, together with the phylogenetic analysis and the experimental data, point to several potential RNA-RNA contacts which depend on the structural and sequence context in which they occur.

The New World of ribozymes

The number of RNA molecules that have novel catalytic activities has dramatically increased during the past two years. This ribozymic boom is not due to the discovery of additional examples of natural ribozymes but rather to the development of artificial ribozymes isolated by in vitro selection and evolution techniques. The structural and functional complexities of these artificial ribozymes, however, do not match those of the larger natural ribozymes. The understanding of both RNA structure and catalysis performed by natural and artificial ribozymes paves the way for the creation of RNA molecules that are able to efficiently catalyze more complex reactions.

New loop-loop tertiary interactions in self-splicing introns of subgroup IC and ID: a complete 3D model of the Tetrahymena thermophila ribozyme

BACKGROUND

Group I introns self-splice via two consecutive trans-esterification reactions in the presence of guanosine cofactor and magnesium ions. Comparative sequence analysis has established that a catalytic core of about 120 nucleotides is conserved in all known group I introns. This core is generally not sufficient for activity, however, and most self-splicing group I introns require non-conserved peripheral elements to stabilize the complete three-dimensional (3D) structure. The physico-chemical properties of group I introns make them excellent systems for unraveling the structural basis of the RNA-RNA interactions responsible for promoting the self-assembly of complex RNAs.

RESULTS

We present phylogenetic and experimental evidence for the existence of three additional tertiary base pairings between hairpin loops within peripheral components of subgroup IC1 and ID introns. Each of these new long range interactions, called P13, P14 and P16, involves a terminal loop located in domain 2. Although domains 2 of IC and ID introns share very strong sequence similarity, their terminal loops interact with domains 5 and 9 (subgroup IC1) and domain 6 (subgroup ID). Based on these tertiary contacts, comparative sequence analysis, and published experimental results such as Fe(II)-EDTA protection patterns, we propose 3D models for two entire group I introns, the subgroup IC1 intron in the large ribosomal precursor RNA of Tetrahymena thermophila and the SdCob.1 subgroup ID intron found in the cytochrome b gene of Saccharomyces douglasii.

CONCLUSIONS

Three-dimensional models of group I introns belonging to four different subgroups are now available. They all emphasize the modular and hierarchical organization of the architecture of group I introns and the widespread use of base-pairings between terminal hairpin loops for stabilizing the folded and active structures of large and complex RNA molecules.