Visualization of a tertiary structural domain of the Tetrahymena group I intron by electron microscopy

High-throughput determination of RNA structure by proximity ligation

We present an unbiased method to globally resolve RNA structures through pairwise contact measurements between interacting regions. RNA Proximity Ligation (RPL) uses proximity ligation of native RNA followed by deep sequencing to yield chimeric reads with ligation junctions in the vicinity of structurally proximate bases. We apply RPL in both baker's yeast (Saccharomyces cerevisiae) and human cells and generate contact probability maps for ribosomal and other abundant RNAs, including yeast snoRNAs, the RNA subunit of the signal recognition particle, and the yeast U2 spliceosomal RNA homolog. RPL measurements correlate with established secondary structures for these RNA molecules, including stem-loop structures and long-range pseudoknots. We anticipate that RPL will complement the current repertoire of computational and experimental approaches in enabling the high-throughput determination of secondary and tertiary RNA structures.

Many ways to loop DNA

In the 1960s, I developed methods for directly visualizing DNA and DNA-protein complexes using an electron microscope. This made it possible to examine the shape of DNA and to visualize proteins as they fold and loop DNA. Early applications included the first visualization of true nucleosomes and linkers and the demonstration that repeating tracts of adenines can cause a curvature in DNA. The binding of DNA repair proteins, including p53 and BRCA2, has been visualized at three- and four-way junctions in DNA. The trombone model of DNA replication was directly verified, and the looping of DNA at telomeres was discovered.

Visualizing large RNA molecules in solution

Single-stranded RNAs (ssRNAs) longer than a few hundred nucleotides do not have a unique structure in solution. Their equilibrium properties therefore reflect the average of an ensemble of structures. We use cryo-electron microscopy to image projections of individual long ssRNA molecules and characterize the anisotropy of their ensembles in solution. A flattened prolate volume is found to best represent the shapes of these ensembles. The measured sizes and anisotropies are in good agreement with complementary determinations using small-angle X-ray scattering and coarse-grained molecular dynamics simulations. A long viral ssRNA is compared with shorter noncoding transcripts to demonstrate that prolate geometry and flatness are generic properties independent of sequence length and origin. The anisotropy persists under physiological as well as low-ionic-strength conditions, revealing a direct correlation between secondary structure asymmetry and 3D shape and size. We discuss the physical origin of the generic anisotropy and its biological implications.

Structure of long telomeric RNA transcripts: the G-rich RNA forms a compact repeating structure containing G-quartets

Recent studies have identified RNA transcripts arising from mammalian telomeres with the transcript of the C-rich strand (r(5'-UUAGGG-3')(n)) being much more abundant than transcripts of the G-rich strand. Here we used transmission electron microscopy, CD, and nuclease digestion to investigate the structure of approximately 640-nucleotide (nt) RNA transcripts of the C-rich and G-rich strands of mammalian telomeric DNA. The CD spectrum of the C-rich RNA in low salt (10 mm KCl) or high salt (100 mm KCl) was typical of mixed sequence RNA, whereas the CD spectrum for the G-rich RNA differed with changes characteristic of parallel G-quadruplexes at the higher salt concentration. Electron microscopy visualization of the C-rich RNA revealed relatively extended unstructured molecules 59.7 +/- 17.8 nm in length and with a width consistent with single-stranded RNA following metal coating. In contrast, the G-rich RNA was observed as round particles and short, thick rods with the rods being most prevalent in high salt conditions and absent in low salt. The rods were 22.7 +/- 4.8 nm in length and 7.6 nm in width. Digestion of the G-rich RNA with T1 RNA nuclease revealed a ladder of bands whose sizes were integral multiples of 24 nt plus a 4-nt overhang. These observations suggest a model in which G-rich telomeric RNA folds into chains of particles each consisting of four (UUAGGG) repeats stabilized by parallel G-quartets and joined by UUA linkers. These chains further condense to form short rods and round particles.

Perturbed folding kinetics of circularly permuted RNAs with altered topology

The folding pathway of the Tetrahymena ribozyme correlates inversely with the sequence distance between native interactions, or contact order. The rapidly folding P4-P6 domain has a low contact order, while the slowly folding P3-P7 region has a high contact order. To examine the role of topology and contact order in RNA folding, we screened for circular permutants of the ribozyme that retain catalytic activity. Permutants beginning in the P4-P6 domain fold 5 to 20 times more slowly than the wild-type ribozyme. By contrast, 50% of a permuted RNA that disjoins a non-native interaction in P3 folds tenfold faster than the wild-type ribozyme. Hence, the probability of rapidly folding to the native state depends on the topology of tertiary domains.

A comprehensive characterization of a group IB intron and its encoded maturase reveals that protein-assisted splicing requires an almost intact intron RNA

The group I intron (AnCOB) of the mitochondrial apocytochrome b gene from Aspergillus nidulans encodes a bi-functional maturase protein that is also a DNA endonuclease. Although the AnCOB intron self-splices, the encoded maturase protein greatly facilitates splicing, in part, by stabilizing RNA tertiary structure. To determine their role in self-splicing and in protein-assisted splicing, several peripheral RNA sub-domains in the 313 nucleotide intron were deleted (P2, P9, P9.1) or truncated (P5ab, P6a). The sequence in two helices (P2 and P9) was also inverted. Except for P9, the deleted regions are not highly conserved among group I introns and are often dispensable for catalytic activity. Nevertheless, despite the very tight binding of AnCOB RNA to the maturase and the high activity of the bimolecular complex (the rate of 5' splice-site cleavage was >20 min(-1) with guanosine as the cofactor), the intron was surprisingly sensitive to these modifications. Several mutations inactivated splicing completely and virtually all impaired splicing to varying degrees. Mutants containing comparatively small deletions in various regions of the intron significantly decreased binding affinity (generally >10(4)-fold), indicating that none of the domains that remained constitutes the primary recognition site of the maturase. The data argue that tight binding requires tertiary interactions that can be maintained by only a relatively intact intron RNA, and that the binding mechanism of the maturase differs from those of two other well-characterized group I intron splicing factors, CYT-18 and Cpb2. A model is proposed in which the protein promotes widespread cooperative folding of an RNA lacking extensive initial tertiary structure.

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.

Gold-tagged RNA-A probe for macromolecular assemblies

Ribonucleic acids (RNAs) play a key role in many fundamental life processes. These polymers are often found complexed with proteins in extremely large particles whose molecular mass may reach several millions of daltons (e.g., ribosomes, spliceosomes, and viruses). Structural studies of such RNA-protein complexes should help elucidate their mode of action. For the structural analyses of many macromolecular assemblies, electron microscopy (EM) has served an instrumental role. However, localization by EM of RNA within biological complexes is not yet a straightforward undertaking. Here we describe a methodology for the covalent tagging of RNA molecules with gold clusters, thereby enabling their direct visualization by microscopical methods. Our strategy involves transcription in vitro of RNAs that carry free thiol groups, using ribonucleoside triphosphate analogs containing a substituent with a terminal thiol group on their heterocyclic ring. This synthesis is followed by coupling of gold clusters to the thiolated transcript through a maleimido group. Visualization of such gold-tagged RNAs by transmission electron microscopy showed spots of gold clusters, with a diameter of 1-2 nm, arranged at nearly regular distances on an imaginary curve that presumably corresponds to the RNA chain. This assignment was corroborated by atomic force microscopy that exhibited images of RNA chains in which knob-like structures, whose height corresponds to the diameter of the gold clusters, were clearly seen. This study demonstrates the potential use of nucleic acids that are covalently labeled with gold clusters for the structural characterization of protein-RNA complexes.

Visualization of double-stranded RNAs from the myotonic dystrophy protein kinase gene and interactions with CUG-binding protein

Myotonic dystrophy (DM) is associated with a (CTG) (n) triplet repeat expansion in the 3'-untranslated region of the myotonic dystrophy protein kinase (DMPK) gene. Using electron microscopy, we visualized large RNAs containing up to 130 CUG repeats and studied the binding of purified CUG-binding protein (CUG-BP) to these RNAs. Electron microscopic examination revealed perfect double-stranded (ds)RNA segments whose lengths were that expected for duplex RNA. The RNA dominant mutation model for DM pathogenesis predicts that the expansion mutation acts at the RNA level by forming long dsRNAs that sequester certain RNA-binding proteins. To test this model, we examined the subcellular distribution and RNA-binding properties of CUG-BP. While previous studies have demonstrated that mutant DMPK transcripts accumu-late in nuclear foci, the localization pattern of CUG-BP in both normal and DM cells was similar. Although CUG-BP in nuclear extracts preferentially photocrosslinked to DMPK transcripts, this binding was not proportional to (CUG) (n) repeat size. Moreover, CUG-BP localized to the base of the RNA hairpin and not along the stem, as visualized by electron micro-scopy. These results provide the first visual evidence that the DM expansion forms an RNA hairpin structure and suggest that CUG-BP is unlikely to be a sequestered factor.

Magnesium-dependent folding of self-splicing RNA: exploring the link between cooperativity, thermodynamics, and kinetics

Folding of the Tetrahymena self-splicing RNA into its active conformation involves a set of discrete intermediate states. The Mg2+-dependent equilibrium transition from the intermediates to the native structure is more cooperative than the formation of the intermediates from the unfolded states. We show that the degree of cooperativity is linked to the free energy of each transition and that the rate of the slow transition from the intermediates to the native state decreases exponentially with increasing Mg2+ concentration. Monovalent salts, which stabilize the folded RNA nonspecifically, induce states that fold in less than 30 s after Mg2+ is added to the RNA. A simple model is proposed that predicts the folding kinetics from the Mg2+-dependent change in the relative stabilities of the intermediate and native states.

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.

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.

Flexibility of RNA

One of the fundamental properties of the RNA helix is its intrinsic resistance to bend- or twist-deformations. Results of a variety of physical measurements point to a persistence length of 700-800 A for double-stranded RNA in the presence of magnesium cations, approximately 1.5-2.0-fold larger than the corresponding value for DNA. Although helix flexibility represents an important, quantifiable measure of the forces of interaction within the helix, it must also be considered in describing conformational variation of nonhelix elements (e.g. internal loops, branches), since the latter always reflect the properties of the flanking helices; that is, such elements are never completely rigid. For one important element of tertiary structure, namely, the core of yeast tRNAPhe, the above consideration has led to the conclusion that the core is not substantially more flexible than an equivalent length of pure helix.

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.

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.

The kinetic folding pathway of the Tetrahymena ribozyme reveals possible similarities between RNA and protein folding

We have probed the nature of the individual kinetic steps in the folding of the Tetrahymena ribozyme by studying the folding kinetics of mutant ribozymes. After rapid formation of the first structural subdomain, a slow step precedes stable formation of the second subdomain. The two central helices of the second subdomain form in an interdependent manner, and this structural subunit therefore also constitutes a kinetic folding unit. The slow folding step includes formation of tertiary interactions in a triple-helical scaffold that orients the two subdomains of the RNA. The rapid and early formation of short range secondary structure, the hierarchical formation of kinetic folding units corresponding to structural subdomains, and the formation of tertiary interactions between subdomains late during the folding process appear to be common features of the folding mechanism for both RNA and proteins.

The P9.1-P9.2 peripheral extension helps guide folding of the Tetrahymena ribozyme

We have previously proposed a hierarchical model for the folding mechanism of the Tetrahymena ribozyme that may illustrate general features of the folding pathways of large RNAs. While the role of elements in the conserved catalytic core of this ribozyme during the folding process is beginning to emerge, the participation of non-conserved peripheral extensions in the kinetic folding mechanism has not yet been addressed. We now show that the 3'-terminal P9.1-P9.2 extension of the Tetrahymena ribozyme plays an important role during the folding process and appears to guide formation of the catalytic core.

Primary and secondary structure analyses of the rDNA group-I introns of the Zygnematales (Charophyta)

The Zygnematales (Charophyta) contain a group-I intron (subgroupIC1) within their nuclear-encoded small subunit ribosomal DNA (SSU rDNA) coding region. This intron, which is inserted after position 1506 (relative to the SSU rDNA of Escherichia coli), is proposed to have been vertically inherited since the origin of the Zygnematales approximately 350-400 million years ago. Primary and secondary structure analyses were carried out to model group-I intron evolution in the Zygnematales. Secondary structure analyses support genetic data regarding sequence conservation within regions known to be functionally important for in vitro self-splicing of group-I introns. Comparisons of zygnematalean group-I intron secondary structures also provided some new insights into sequences that may have important roles in in vivo RNA splicing. Sequence analyses showed that sequence divergence rates and the nucleotide compositions of introns and coding regions within any one taxon varied widely, suggesting that the "1506" group-I introns and rDNA coding regions in the Zygnematales evolve independently.

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.

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.

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

Coaxial stacking of helical elements is a determinant of three-dimensional structure in RNA. In the catalytic center of the Tetrahymena group I intron, helices P4 and P6 are part of a tertiary structural domain that folds independently of the remainder of the intron. When P4 and P6 were fused with a phosphodiester linkage, the resulting RNA retained the detailed tertiary interactions characteristic of the native P4-P6 domain and even required lower magnesium ion concentrations for folding. These results indicate that P4 and P6 are coaxial in the P4-P6 domain and, therefore, in the native ribozyme. Helix fusion could provide a general method for identifying pairs of coaxially stacked helices in biological RNA molecules.

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.