Coaxially stacked RNA helices in the catalytic center of the Tetrahymena ribozyme

Crystal structure of a cap-independent translation enhancer RNA

In eukaryotic messenger RNAs, the 5' cap structure binds to the translation initiation factor 4E to facilitate early stages of translation. Although many plant viruses lack the 5' cap structure, some contain cap-independent translation elements (CITEs) in their 3' untranslated region. The PTE (Panicum mosaic virus translation element) class of CITEs contains a G-rich asymmetric bulge and a C-rich helical junction that were proposed to interact via formation of a pseudoknot. SHAPE analysis of PTE homologs reveals a highly reactive guanosine residue within the G-rich region proposed to mediate eukaryotic initiation factor 4E (eIF4E) recognition. Here we have obtained the crystal structure of the PTE from Pea enation mosaic virus 2 (PEMV2) RNA in complex with our structural chaperone, Fab BL3-6. The structure reveals that the G-rich and C-rich regions interact through a complex network of interactions distinct from those expected for a pseudoknot. The motif, which contains a short parallel duplex, provides a structural mechanism for how the guanosine is extruded from the core stack to enable eIF4E recognition. Homologous PTE elements harbor a G-rich bulge and a three-way junction and exhibit covariation at crucial positions, suggesting that the PEMV2 tertiary architecture is conserved among these homologs.

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway

The past decade has seen a wealth of 3D structural information about complex structured RNAs and identification of functional intermediates. Nevertheless, developing a complete and predictive understanding of the folding and function of these RNAs in biology will require connection of individual rate and equilibrium constants to structural changes that occur in individual folding steps and further relating these steps to the properties and behavior of isolated, simplified systems. To accomplish these goals we used the considerable structural knowledge of the folded, unfolded, and intermediate states of P4-P6 RNA. We enumerated structural states and possible folding transitions and determined rate and equilibrium constants for the transitions between these states using single-molecule FRET with a series of mutant P4-P6 variants. Comparisons with simplified constructs containing an isolated tertiary contact suggest that a given tertiary interaction has a stereotyped rate for breaking that may help identify structural transitions within complex RNAs and simplify the prediction of folding kinetics and thermodynamics for structured RNAs from their parts. The preferred folding pathway involves initial formation of the proximal tertiary contact. However, this preference was only ∼10 fold and could be reversed by a single point mutation, indicating that a model akin to a protein-folding contact order model will not suffice to describe RNA folding. Instead, our results suggest a strong analogy with a modified RNA diffusion-collision model in which tertiary elements within preformed secondary structures collide, with the success of these collisions dependent on whether the tertiary elements are in their rare binding-competent conformations.

Visualizing the formation of an RNA folding intermediate through a fast highly modular secondary structure switch

Intermediates play important roles in RNA folding but can be difficult to characterize when short-lived or not significantly populated. By combining (15)N relaxation dispersion NMR with chemical probing, we visualized a fast (kex=k1+k-1≈423 s(-1)) secondary structural switch directed towards a low-populated (∼3%) partially folded intermediate in tertiary folding of the P5abc subdomain of the 'Tetrahymena' group I intron ribozyme. The secondary structure switch changes the base-pairing register across the P5c hairpin, creating a native-like structure, and occurs at rates of more than two orders of magnitude faster than tertiary folding. The switch occurs robustly in the absence of tertiary interactions, Mg(2+) or even when the hairpin is excised from the three-way junction. Fast, highly modular secondary structural switches may be quite common during RNA tertiary folding where they may help smoothen the folding landscape by allowing folding to proceed efficiently via additional pathways.

Promiscuous RNA binding by Polycomb repressive complex 2

Polycomb repressive complex-2 (PRC2) is a histone methyltransferase required for epigenetic silencing during development and cancer. Long non-coding RNAs (lncRNAs) recruit PRC2 to chromatin, but the general role of RNA in maintaining repressed chromatin is unknown. Here we measure the binding constant of human PRC2 to various RNAs and find comparable affinity for human lncRNAs targeted by PRC2 and irrelevant transcripts from ciliates and bacteria. PRC2 binding is size-dependent, with lower affinity for shorter RNAs. In vivo, PRC2 predominantly occupies repressed genes; PRC2 is also associated with active genes, but most of these are not regulated by PRC2. These findings support a model in which promiscuous binding of PRC2 to RNA transcripts allows it to scan for target genes that have escaped repression, leading to maintenance of the repressed state. Such RNAs may also provide a decoy for PRC2.

2'-Fluoro substituents can mimic native 2'-hydroxyls within structured RNA

The ability of fluorine in a C-F bond to act as a hydrogen bond acceptor is controversial. To test such ability in complex RNA macromolecules, we have replaced native 2'-OH groups with 2'-F and 2'-H groups in two related systems, the Tetrahymena group I ribozyme and the ΔC209 P4-P6 RNA domain. In three cases the introduced 2'-F mimics the native 2'-OH group, suggesting that the fluorine atom can accept a hydrogen bond. In each of these cases the native hydroxyl group interacts with a purine exocyclic amine. Our results give insight about the properties of organofluorine and suggest a possible general biochemical signature for tertiary interactions between 2'-hydroxyl groups and exocyclic amino groups within RNA.

Turning limited experimental information into 3D models of RNA

Our understanding of RNA functions in the cell is evolving rapidly. As for proteins, the detailed three-dimensional (3D) structure of RNA is often key to understanding its function. Although crystallography and nuclear magnetic resonance (NMR) can determine the atomic coordinates of some RNA structures, many 3D structures present technical challenges that make these methods difficult to apply. The great flexibility of RNA, its charged backbone, dearth of specific surface features, and propensity for kinetic traps all conspire with its long folding time, to challenge in silico methods for physics-based folding. On the other hand, base-pairing interactions (either in runs to form helices or isolated tertiary contacts) and motifs are often available from relatively low-cost experiments or informatics analyses. We present RNABuilder, a novel code that uses internal coordinate mechanics to satisfy user-specified base pairing and steric forces under chemical constraints. The code recapitulates the topology and characteristic L-shape of tRNA and obtains an accurate noncrystallographic structure of the Tetrahymena ribozyme P4/P6 domain. The algorithm scales nearly linearly with molecule size, opening the door to the modeling of significantly larger structures.

Generation and development of RNA ligase ribozymes with modular architecture through "design and selection"

In vitro selection with long random RNA libraries has been used as a powerful method to generate novel functional RNAs, although it often requires laborious structural analysis of isolated RNA molecules. Rational RNA design is an attractive alternative to avoid this laborious step, but rational design of catalytic modules is still a challenging task. A hybrid strategy of in vitro selection and rational design has been proposed. With this strategy termed “design and selection,” new ribozymes can be generated through installation of catalytic modules onto RNA scaffolds with defined 3D structures. This approach, the concept of which was inspired by the modular architecture of naturally occurring ribozymes, allows prediction of the overall architectures of the resulting ribozymes, and the structural modularity of the resulting ribozymes allows modification of their structures and functions. In this review, we summarize the design, generation, properties, and engineering of four classes of ligase ribozyme generated by design and selection.

Structural inference of native and partially folded RNA by high-throughput contact mapping

The biological behaviors of ribozymes, riboswitches, and numerous other functional RNA molecules are critically dependent on their tertiary folding and their ability to sample multiple functional states. The conformational heterogeneity and partially folded nature of most of these states has rendered their characterization by high-resolution structural approaches difficult or even intractable. Here we introduce a method to rapidly infer the tertiary helical arrangements of large RNA molecules in their native and non-native solution states. Multiplexed hydroxyl radical (.OH) cleavage analysis (MOHCA) enables the high-throughput detection of numerous pairs of contacting residues via random incorporation of radical cleavage agents followed by two-dimensional gel electrophoresis. We validated this technology by recapitulating the unfolded and native states of a well studied model RNA, the P4-P6 domain of the Tetrahymena ribozyme, at subhelical resolution. We then applied MOHCA to a recently discovered third state of the P4-P6 RNA that is stabilized by high concentrations of monovalent salt and whose partial order precludes conventional techniques for structure determination. The three-dimensional portrait of a compact, non-native RNA state reveals a well ordered subset of native tertiary contacts, in contrast to the dynamic but otherwise similar molten globule states of proteins. With its applicability to nearly any solution state, we expect MOHCA to be a powerful tool for illuminating the many functional structures of large RNA molecules and RNA/protein complexes.

Quantitative analysis of protein-RNA interactions by gel mobility shift

The gel mobility shift assay is routinely used to visualize protein-RNA interactions. Its power resides in the ability to resolve free from bound RNA with high resolution in a gel matrix. We review the quantitative application of this approach to elucidate thermodynamic properties of protein-RNA complexes. Assay designs for titration, competition, and stoichiometry experiments are presented for two unrelated model complexes.

Replication fork regression in repetitive DNAs

Among several different types of repetitive sequences found in the human genome, this study has examined the telomeric repeat, necessary for the protection of chromosome termini, and the disease-associated triplet repeat (CTG)·(CAG)_n. Evidence suggests that replication of both types of repeats is problematic and that a contributing factor is the repetitive nature of the DNA itself. Here we have used electron microscopy to investigate DNA structures formed at replication forks on large model DNAs containing these repeat sequences, in an attempt to elucidate the contributory effect that these repetitive DNAs may have on their replication. Visualization of the DNA revealed that there is a high propensity for a paused replication fork to spontaneously regress when moving through repetitive DNAs, and that this results in a four-way chickenfoot intermediate that could present a significant block to replication in vivo, possibly leading to unwanted recombination events, amplifications or deletions.

Evaluation of the conformational switch model for alfalfa mosaic virus RNA replication

Key elements of the conformational switch model describing regulation of alfalfa mosaic virus (AMV) replication (R. C. Olsthoorn, S. Mertens, F. T. Brederode, and J. F. Bol, EMBO J. 18:4856-4864, 1999) have been tested using biochemical assays and functional studies in nontransgenic protoplasts. Although comparative sequence analysis suggests that the 3' untranslated regions of AMV and ilarvirus RNAs have the potential to fold into pseudoknots, we were unable to confirm that a proposed pseudoknot forms or has a functional role in regulating coat protein-RNA binding or viral RNA replication. Published work has suggested that the pseudoknot is part of a tRNA-like structure (TLS); however, we argue that the canonical sequence and functional features that define the TLS are absent. We suggest here that the absence of the TLS correlates directly with the distinctive requirement for coat protein to activate replication in these viruses. Experimental data are evidence that elevated magnesium concentrations proposed to stabilize the pseudoknot structure do not block coat protein binding. Additionally, covarying nucleotide changes proposed to reestablish pseudoknot pairings do not rescue replication. Furthermore, as described in the accompanying paper (L. M. Guogas, S. M. Laforest, and L. Gehrke, J. Virol. 79:5752-5761, 2005), coat protein is not, by definition, inhibitory to minus-strand RNA synthesis. Rather, the activation of viral RNA replication by coat protein is shown to be concentration dependent. We describe the 3' organization model as an alternate model of AMV replication that offers an improved fit to the available data.

A preformed compact ribosome-binding domain in the cricket paralysis-like virus IRES RNAs

The internal ribosome site RNA of the cricket paralysis-like viruses (CrPV-like) binds directly to the ribosome, assembling the translation machinery without initiation factors. This mechanism does not require initiator tRNA, and translation starts from a non-AUG codon. A wealth of biochemical data has yielded a working model for this process, but the three-dimensional structure and biophysical characteristics of the unbound CrPV-like IRES RNAs are largely unexplored. Here, we demonstrate that the CrPV-like IRESes prefold into a two-part structure in the presence of magnesium ions. The largest part is a prefolded compact RNA domain that shares folding and structural characteristics with other compactly folded RNAs such as group I intron RNAs and RNase P RNA. Chemical probing reveals that the CrPV-like IRES' compact domain contains RNA helices that are packed tightly enough to exclude solvent, and analytical ultracentrifugation indicates a large change in the shape of the IRES upon folding. Formation of this compact domain is necessary for binding of the 40S subunit, and the structural organization of the unbound IRES RNA is consistent with the hypothesis that the IRES is functionally and structurally preorganized before ribosome binding.

Crystal structure of a group I intron splicing intermediate

A recently reported crystal structure of an intact bacterial group I self-splicing intron in complex with both its exons provided the first molecular view into the mechanism of RNA splicing. This intron structure, which was trapped in the state prior to the exon ligation reaction, also reveals the architecture of a complex RNA fold. The majority of the intron is contained within three internally stacked, but sequence discontinuous, helical domains. Here the tertiary hydrogen bonding and stacking interactions between the domains, and the single-stranded joiner segments that bridge between them, are fully described. Features of the structure include: (1) A pseudoknot belt that circumscribes the molecule at its longitudinal midpoint; (2) two tetraloop-tetraloop receptor motifs at the peripheral edges of the structure; (3) an extensive minor groove triplex between the paired and joiner segments, P6-J6/6a and P3-J3/4, which provides the major interaction interface between the intron's two primary domains (P4-P6 and P3-P9.0); (4) a six-nucleotide J8/7 single stranded element that adopts a mu-shaped structure and twists through the active site, making critical contacts to all three helical domains; and (5) an extensive base stacking architecture that realizes 90% of all possible stacking interactions. The intron structure was validated by hydroxyl radical footprinting, where strong correlation was observed between experimental and predicted solvent accessibility. Models of the pre-first and pre-second steps of intron splicing are proposed with full-sized tRNA exons. They suggest that the tRNA undergoes substantial angular motion relative to the intron between the two steps of splicing.

Generation of a catalytic module on a self-folding RNA

It is theoretically possible to obtain a catalytic site of an artificial ribozyme from a random sequence consisting of a limited numbers of nucleotides. However, this strategy has been inadequately explored. Here, we report an in vitro selection technique that exploits modular construction of a structurally constrained RNA to acquire a catalytic site for RNA ligation from a short random sequence. To practice the selection, a sequence of 30 nucleotides was located close to the putative reaction site in a derivative of a naturally occurring self-folding RNA whose crystal structure is known. RNAs whose activity depended on the starting three-dimensional structure were selected with 3'-5' ligation specificity, indicating that the strategy can be used to acquire a variety of catalytic sites and other functional RNA modules.

RNA crystallization

RNA molecules may be crystallized using variations of the methods developed for protein crystallography. As the technology has become available to synthesize and purify RNA molecules in the quantities and with the quality that is required for crystallography, the field of RNA structure has exploded. The first consideration when crystallizing an RNA is the sequence, which may be varied in a rational way to enhance crystallizability or prevent formation of alternate structures. Once a sequence has been designed, the RNA may be synthesized chemically by solid-state synthesis or it may be produced enzymatically using RNA polymerase and an appropriate DNA template. Purification of milligram quantities of RNA can be accomplished by HPLC or gel electrophoresis. As with proteins, crystallization of RNA is usually accomplished by vapor diffusion techniques. There are several considerations that are either unique to RNA crystallization or more important for RNA crystallization. Techniques for design, synthesis, purification, and crystallization of RNAs will be reviewed here.

Design, construction, and analysis of a novel class of self-folding RNA

RNA can play multiple biological roles through use of its three-dimensional (3-D) structures. Recent advances in RNA structural biology have revealed that complex RNA 3D structures are assemblages of double-stranded helices with a variety of tertiary structural motifs. By employing RNA tertiary structural motifs together with the helices, we designed a novel class of self-folding RNA. In RNA composed of three helices (P1, P2, and P3), P1 interacts with P3 via a tetraloop-receptor interaction and P2 forms consecutive base-triples. Two designed RNAs of this class were prepared and their folding properties indicate that they form defined tertiary structures as designed. These RNAs may be used as modular units for constructing artificial ribozymes or nanometer-scale materials.

A small structural element, Pc-J5/5a, plays dual roles in a group IC1 intron RNA

The P4-P6 domain of group IC1 intron ribozymes such as that of the Tetrahymena autonomously folds into a hairpin-shaped structure in which the J5/5a region serves as a hinge. Phylogenetic comparisons of these IC1 introns suggested that the J5/5a region (termed Pc-J5/5a motif) in a subclass of IC1 introns such as the one from Pneumocystis carinii functions not only as a hinge but also as a receptor for a GAAA-tetraloop. We investigated the role of this motif by transplanting the structural unit, Pc-J5/5a motif, of Pneumocystis carinii into the P4-P6 domain of the Tetrahymena intron. The results showed that the Pc-J5/5a motif binds to a GAAA loop with high affinity and also facilitates the bending of the Tetrahymena P4-P6 domain more positively than the original J5/5a region.

Determining the influence of structure on hybridization using oligonucleotide arrays

We have studied the effects of structure on nucleic acid heteroduplex formation by analyzing hybridization of tRNAphe to a complete set of complementary oligonucleotides, ranging from single nucleotides to dodecanucleotides. The analysis points to features in tRNA that determine heteroduplex yield. All heteroduplexes that give high yield include both double-stranded stems as well as single-stranded regions. Bases in the single-stranded regions are stacked onto the stems, and heteroduplexes terminate at potential interfaces for coaxial stacking. Heteroduplex formation is disfavored by sharp turns or a lack of helical order in single-stranded regions, competition from bases displaced from a stem, and stable tertiary interactions. The study is relevant to duplex formation on oligonucleotide microarrays and to antisense technologies.

Myotonic dystrophy: molecular windows on a complex etiology

Myotonic dystrophy (DM) is the most common form of adult onset muscular dystrophy, with an incidence of approximately 1 in 8500 adults. DM is caused by an expanded number of trinucleotide repeats in the 3'-untranslated region (UTR) of a cAMP-dependent protein kinase (DM protein kinase, DMPK). Although a large number of transgenic animals have been generated with different gene constructions and knock-outs, none of them faithfully recapitulates the multisystemic and often severe phenotype seen in human patients. The transgenic data suggest that myotonic dystrophy is not caused simply by a biochemical deficiency or abnormality in the DM kinase gene product. Emerging studies suggest that two novel pathogenetic mechanisms may play a role in the disease: the expanded repeats appear to cause haploinsufficiency of a neighboring homeobox gene and also abnormal DMPK RNA appears to have a detrimental effect on RNA homeostasis. The complex, multisystemic phenotype may reflect an underlying multifaceted molecular pathophysiology: the facial dysmorphology may be due to pattern defects caused by haploinsufficiency of the homeobox gene, while the muscle disease and endocrine abnormalities may be due to both altered RNA metabolism and deficiency of the cAMP DMPK protein.

RNA seeing double: close-packing of helices in RNA tertiary structure

Structured RNA molecules play essential roles in RNA processing, chromosome maintenance and protein biosynthesis. RNA necessarily uses different strategies than proteins for folding and assembly of complex architectures. The RNA-folding problem is largely an issue of helical packing: how does RNA organize and pack short, double-helical segments to produce active sites and recognition motifs for proteins? Noncanonical base pairs, metal ions and 2'-hydroxyl groups are key elements in RNA higher-order structure formation.

Time-resolved synchrotron X-ray "footprinting", a new approach to the study of nucleic acid structure and function: application to protein-DNA interactions and RNA folding

Hydroxyl radicals (.OH) can cleave the phosphodiester backbone of nucleic acids and are valuable reagents in the study of nucleic acid structure and protein-nucleic acid interactions. Irradiation of solutions by high flux "white light" X-ray beams based on bending magnet beamlines at the National Synchrotron Light Source (NSLS) yields sufficient concentrations of .OH so that quantitative nuclease protection ("footprinting") studies of DNA and RNA can be conducted with a duration of exposure in the range of 50 to 100 ms. The sensitivity of DNA and RNA to X-ray mediated .OH cleavage is equivalent. Both nucleic acids are completely protected from synchrotron X-ray induced cleavage by the presence of thiourea in the sample solution, demonstrating that cleavage is suppressed by a free radical scavenger. The utility of this time-dependent approach to footprinting is demonstrated with a synchrotron X-ray footprint of a protein-DNA complex and by a time-resolved footprinting analysis of the Mg(2+)-dependent folding of the Tetrahymena thermophilia L-21 ScaI ribozyme RNA. Equilibrium titrations reveal differences among the ribozyme domains in the cooperativity of Mg(2+)-dependent .OH protection. RNA .OH protection progress curves were obtained for several regions of the ribozyme over timescales of 30 seconds to several minutes. Progress curves ranging from > or = 3.5 to 0.4 min-1 were obtained for the P4-P6 and P5 sub-domains and the P3-P7 domain, respectively. The .OH protection progress curves have been correlated with the available biochemical, structural and modeling data to generate a model of the ribozyme folding pathway. Rate differences observed for specific regions within domains provide evidence for steps in the folding pathway not previously observed. Synchrotron X-ray footprinting is a new approach of general applicability for the study of time-resolved structural changes of nucleic acid conformation and protein-nucleic acid complexes.

A common RNA structural motif involved in the internal initiation of translation of cellular mRNAs

The 5'-non-translated regions (5'NTR) of human immunoglobulin heavy chain binding protein (BiP), Antennapedia (Antp) ofDrosophilaand human fibroblast growth factor 2 (FGF-2) mRNAs are reported to mediate translation initiation by an internal ribosome binding mechanism. In this study, we investigate predicted features of the higher order structures folded in these 5'NTR sequences. Statistical analyses of RNA folding detected a 92 nt unusual folding region (UFR) from 129 to 220, close to the initiator AUG in the BiP mRNA. Details of the structural analyses show that the UFR forms a Y-type stem-loop structure with an additional stem-loop in the 3'-end resembling the common structure core found in the internal ribosome entry site (IRES) elements of picornavirus. The Y-type structural motif is also conserved among a number of divergent BiP mRNAs. We also find two RNA elements in the 5'-leader sequence of human FGF-2. The first RNA element (96 nt) is 2 nt upstream of the first CUG start codon located in the reported IRES element of human FGF-2. The second (107 nt) is immediately upstream of the authentic initiator AUG of the main open reading frame. Intriguingly, the folded RNA structural motif in the two RNA elements is conserved in other members of FGF family and shares the same structural features as that found in the 5'NTR of divergent BiP mRNAs. We suggest that the common RNA structural motif conserved in the diverse BiP and FGF-2 mRNAs has a general function in the internal ribosome binding mechanism of cellular mRNAs.

Hierarchy and dynamics of RNA folding

The evidence showing that the self-assembly of complex RNAs occurs in discrete transitions, each relating to the folding of sub-systems of increasing size and complexity starting from a state with most of the secondary structure, is reviewed. The reciprocal influence of the concentration of magnesium ions and nucleotide mutations on tertiary structure is analyzed. Several observations demonstrate that detrimental mutations can be rescued by high magnesium concentrations, while stabilizing mutations lead to a lesser dependence on magnesium ion concentration. Recent data point to the central controlling and monitoring roles of RNA-binding proteins that can bind to the different folding stages, either before full establishment of the secondary structure or at the molten globule state before the cooperative transition to the final three-dimensional structure.

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.

A curved RNA helix incorporating an internal loop with G.A and A.A non-Watson-Crick base pairing

The crystal structure of the RNA dodecamer 5'-GGCC(GAAA)GGCC-3' has been determined from x-ray diffraction data to 2.3-A resolution. In the crystal, these oligomers form double helices around twofold symmetry axes. Four consecutive non-Watson-Crick base pairs make up an internal loop in the middle of the duplex, including sheared G.A pairs and novel asymmetric A.A pairs. This internal loop sequence produces a significant curvature and narrowing of the double helix. The helix is curved by 34 degrees from end to end and the diameter is narrowed by 24% in the internal loop. A Mn2+ ion is bound directly to the N7 of the first guanine in the Watson-Crick region following the internal loop and the phosphate of the preceding residue. This Mn2+ location corresponds to a metal binding site observed in the hammerhead catalytic RNA.

Length changes in the joining segment between domains 5 and 6 of a group II intron inhibit self-splicing and alter 3' splice site selection

Domain 5 (D5) and domain 6 (D6) are adjacent folded hairpin substructures of self-splicing group II introns that appear to interact within the active ribozyme. Here we describe the effects of changing the length of the 3-nucleotide segment joining D5 to D6 [called J(56)3] on the splicing reactions of intron 5 gamma of the COXI gene of yeast mitochondrial DNA. Shortened variants J(56)0 and J(56)1 were defective in vitro for branching, and the second splicing step was performed inefficiently and inaccurately. The lengthened variant J(56)5 had a milder defect-splicing occurred at a reduced rate but with correct branching and a mostly accurate 3' splice junction choice. Yeast mitochondria were transformed with the J(56)5 allele, and the resulting yeast strain was respiration deficient because of ineffective aI5 gamma splicing. Respiration-competent revertants were recovered, and in one type a single joiner nucleotide was deleted while in the other type a nucleotide of D6 was deleted. Although these revertants still showed partial splicing blocks in vivo and in vitro, including a substantial defect in the second step of splicing, both spliced accurately in vivo. These results establish that a 3-nucleotide J(56) is optimal for this intron, especially for the accuracy of 3' splice junction selection, and indicate that D5 and D6 are probably not coaxially stacked.

Crystal structure of a group I ribozyme domain: principles of RNA packing

Group I self-splicing introns catalyze their own excision from precursor RNAs by way of a two-step transesterification reaction. The catalytic core of these ribozymes is formed by two structural domains. The 2.8-angstrom crystal structure of one of these, the P4-P6 domain of the Tetrahymena thermophila intron, is described. In the 160-nucleotide domain, a sharp bend allows stacked helices of the conserved core to pack alongside helices of an adjacent region. Two specific long-range interactions clamp the two halves of the domain together: a two-Mg2+-coordinated adenosine-rich corkscrew plugs into the minor groove of a helix, and a GAAA hairpin loop binds to a conserved 11-nucleotide internal loop. Metal- and ribose-mediated backbone contacts further stabilize the close side-by-side helical packing. The structure indicates the extent of RNA packing required for the function of large ribozymes, the spliceosome, and the ribosome.

Mice lacking the myotonic dystrophy protein kinase develop a late onset progressive myopathy

Myotonic dystrophy (DM) is an autosomal dominant disorder resulting from the expansion of a CTG repeat in the 3' untranslated region of a putative protein kinase (DMPK). To elucidate the role of DMPK in DM pathogenesis we have developed Dmpk deficient (Dmpk-/-) mice. Dmpk-/-mice develop a late-onset, progressive skeletal myopathy that shares some pathological features with DM. Muscles from mature mice show variation in fibre size, increased fibre degeneration and fibrosis. Adult Dmpk-/-mice show ultrastructural changes in muscle and a 50% decrease in force generation compared to young mice. Our results indicate that DMPK may be necessary for the maintenance of skeletal muscle structure and function and suggest that a decrease in DMPK levels may contribute to DM pathology.

RNA tectonics: towards RNA design

Our understanding of the structural, folding and catalytic properties of RNA molecules has increased enormously in recent years. The discovery of catalytic RNA molecules by Sidney Altman and Tom Cech, the development of in vitro selection procedures, and the recent crystallizations of hammerhead ribozymes and of a large domain of an autocatalytic group 1 intron are some of the milestones that have contributed to the explosion of the RNA field. The availability of a three-dimensional model for the catalytic core of group 1 introns contributed also a heuristic drive toward the development of new techniques and approaches for unravelling RNA architecture, folding and stability. Here, we emphasize the mosaic structure of RNA and review some of the recent literature pertinent to this working framework. In the long run, RNA tectonics aims at constructing combinatorial libraries, using RNA mosaic units for creating molecules with dedicated shapes and properties.

A common structural core in the internal ribosome entry sites of picornavirus, hepatitis C virus, and pestivirus

Cap-independent translations of viral RNAs of enteroviruses and rhinoviruses, cardioviruses and aphthoviruses, hepatitis A and C viruses (HAV and HCV), and pestivirus are initiated by the direct binding of 40S ribosomal subunits to a cis-acting genetic element termed the internal ribosome entry site (IRES) or ribosome landing pad (RLP) in the 5' noncoding region (5'NCR). RNA higher ordered structure models for these IRES elements were derived by a combined approach using thermodynamic RNA folding, Monte Carlo simulation, and phylogenetic comparative analysis. The structural differences among the three groups of picornaviruses arise not only from point mutations, but also from the addition or deletion of structural domains. However, a common core can be identified in the proposed structural models of these IRES elements from enteroviruses and rhinoviruses, cardioviruses and aphthoviruses, and HAV. The common structural core identified within the picornavirus IRES is also conserved in the 5'NCR of the divergent viruses, HCV, and pestiviruses. Furthermore, the proposed structural motif shares a structural feature similar to that observed in the catalytic core of the group 1 intron. The conserved structural motif from these divergent sequences that looks like the common core region of group 1 introns is probably a crucial element involved in the IRES-dependent translation.

Protein facilitation of group I intron splicing by assembly of the catalytic core and the 5' splice site domain

The yeast mitochondrial group I intron b15 undergoes self-splicing at high Mg2+ concentrations, but requires the splicing factor CBP2 for reaction under physiological conditions. Chemical accessibility and UV cross-linking experiments now reveal that self-processing is slow because functional elements are not properly positioned in an active tertiary structure. Folding energy provided by CBP2 drives assembly of two RNA domains that comprise the catalytic core and meditates association of an approximately 100 nt 5' domain that contains the 5' splice site. Thus, the protein assembles RNA secondary structure elements into a specific three-dimensional array while the RNA provides the catalytic center. The division of labor between RNA and protein illustrated by this simple system reveals principles applicable to complex ribonucleoprotein assemblies such as the spliceosome and ribosome.

RNA folding

The number of known motifs for RNA folding and RNA tertiary organization is expanding rapidly as we learn more about the diverse biological functions of RNA. Problems in protein and RNA folding have melded in recent investigations of ribonucleoprotein folding. Theoretical and experimental models are rapidly being developed for the pathways and stabilizing forces involved in RNA folding.