Crystal structure of a phage Twort group I ribozyme–product complex

Structure of a tyrosyl-tRNA synthetase splicing factor bound to a group I intron RNA

The 'RNA world' hypothesis holds that during evolution the structural and enzymatic functions initially served by RNA were assumed by proteins, leading to the latter's domination of biological catalysis. This progression can still be seen in modern biology, where ribozymes, such as the ribosome and RNase P, have evolved into protein-dependent RNA catalysts ('RNPzymes'). Similarly, group I introns use RNA-catalysed splicing reactions, but many function as RNPzymes bound to proteins that stabilize their catalytically active RNA structure. One such protein, the Neurospora crassa mitochondrial tyrosyl-tRNA synthetase (TyrRS; CYT-18), is bifunctional and both aminoacylates mitochondrial tRNA(Tyr) and promotes the splicing of mitochondrial group I introns. Here we determine a 4.5-A co-crystal structure of the Twort orf142-I2 group I intron ribozyme bound to splicing-active, carboxy-terminally truncated CYT-18. The structure shows that the group I intron binds across the two subunits of the homodimeric protein with a newly evolved RNA-binding surface distinct from that which binds tRNA(Tyr). This RNA binding surface provides an extended scaffold for the phosphodiester backbone of the conserved catalytic core of the intron RNA, allowing the protein to promote the splicing of a wide variety of group I introns. The group I intron-binding surface includes three small insertions and additional structural adaptations relative to non-splicing bacterial TyrRSs, indicating a multistep adaptation for splicing function. The co-crystal structure provides insight into how CYT-18 promotes group I intron splicing, how it evolved to have this function, and how proteins could have incrementally replaced RNA structures during the transition from an RNA world to an RNP world.

Molecular modelling of the GIR1 branching ribozyme gives new insight into evolution of structurally related ribozymes

Twin-ribozyme introns contain a branching ribozyme (GIR1) followed by a homing endonuclease (HE) encoding sequence embedded in a peripheral domain of a group I splicing ribozyme (GIR2). GIR1 catalyses the formation of a lariat with 3 nt in the loop, which caps the HE mRNA. GIR1 is structurally related to group I ribozymes raising the question about how two closely related ribozymes can carry out very different reactions. Modelling of GIR1 based on new biochemical and mutational data shows an extended substrate domain containing a GoU pair distinct from the nucleophilic residue that dock onto a catalytic core showing a different topology from that of group I ribozymes. The differences include a core J8/7 region that has been reduced and is complemented by residues from the pre-lariat fold. These findings provide the basis for an evolutionary mechanism that accounts for the change from group I splicing ribozyme to the branching GIR1 architecture. Such an evolutionary mechanism can be applied to other large RNAs such as the ribonuclease P.

Structure-function relationships in RNA and RNP enzymes: recent advances

The structural biology of ribozymes and ribonucleoprotein (RNP) enzymes is now sufficiently advanced that a true dialogue between structural and functional studies is possible. In this review, we consider three important systems in which an integration of structural and biochemical data has recently led to major advances in mechanistic understanding. In the hammerhead ribozyme, application-driven biochemical studies led to the discovery of a key structural interaction that had been omitted from previously-studied constructs. A new crystal structure of the resulting, tertiary-stabilized hammerhead has resolved a remarkable number of longstanding paradoxes in the structure-function relationship of this ribozyme. In the Group I intron ribozyme, a flurry of high-resolution structures has largely confirmed, but in some cases refined or challenged, a detailed model of a metalloenzyme active site that had previously been derived by meticulous quantitative metal ion rescue experiments. Finally, for the peptidyl transferase center of the ribosome, recent biochemical and chemical results motivated by the pioneering crystal structures have suggested a picture of a catalytic mechanism dominated by proximity and orientation effects and substrate-assisted catalysis. These results refocus attention on catalysis as a property of the integrated RNP machinery as a whole, as opposed to a narrow concern with the RNA functional groups in immediate contact with the reactive center.

RNA catalysis: ribozymes, ribosomes, and riboswitches

The catalytic mechanisms employed by RNA are chemically more diverse than initially suspected. Divalent metal ions, nucleobases, ribosyl hydroxyl groups, and even functional groups on metabolic cofactors all contribute to the various strategies employed by RNA enzymes. This catalytic breadth raises intriguing evolutionary questions about how RNA lost its biological role in some cases, but not in others, and what catalytic roles RNA might still be playing in biology.

Searching genomes for ribozymes and riboswitches

A discussion of experimental approaches and theoretical difficulties in the identification of ribozymes with novel catalytic functions.

New regulatory RNAs with complex structures have recently been discovered, among them the first catalytic riboswitch, a gene-regulatory RNA sequence with catalytic activity. Here we discuss some of the experimental approaches and theoretical difficulties attached to the identification of new ribozymes in genomes.

Modulation of individual steps in group I intron catalysis by a peripheral metal ion

Enzymes are complex macromolecules that catalyze chemical reactions at their active sites. Important information about catalytic interactions is commonly gathered by perturbation or mutation of active site residues that directly contact substrates. However, active sites are engaged in intricate networks of interactions within the overall structure of the macromolecule, and there is a growing body of evidence about the importance of peripheral interactions in the precise structural organization of the active site. Here, we use functional studies, in conjunction with published structural information, to determine the effect of perturbation of a peripheral metal ion binding site on catalysis in a well-characterized catalytic RNA, the Tetrahymena thermophila group I ribozyme. We perturbed the metal ion binding site by site-specifically introducing a phosphorothioate substitution in the ribozyme's backbone, replacing the native ligands (the pro-R (P) oxygen atoms at positions 307 and 308) with sulfur atoms. Our data reveal that these perturbations affect several reaction steps, including the chemical step, despite the absence of direct contacts of this metal ion with the atoms involved in the chemical transformation. As structural probing with hydroxyl radicals did not reveal significant change in the three-dimensional structure upon phosphorothioate substitution, the effects are likely transmitted through local, rather subtle conformational rearrangements. Addition of Cd(2+), a thiophilic metal ion, rescues some reaction steps but has deleterious effects on other steps. These results suggest that native interactions in the active site may have been aligned by the naturally occurring peripheral residues and interactions to optimize the overall catalytic cycle.

Ribozymes, riboswitches and beyond: regulation of gene expression without proteins

Although various functions of RNA are carried out in conjunction with proteins, some catalytic RNAs, or ribozymes, which contribute to a range of cellular processes, require little or no assistance from proteins. Furthermore, the discovery of metabolite-sensing riboswitches and other types of RNA sensors has revealed RNA-based mechanisms that cells use to regulate gene expression in response to internal and external changes. Structural studies have shown how these RNAs can carry out a range of functions. In addition, the contribution of ribozymes and riboswitches to gene expression is being revealed as far more widespread than was previously appreciated. These findings have implications for understanding how cellular functions might have evolved from RNA-based origins.

Cyanobacterial ribosomal RNA genes with multiple, endonuclease-encoding group I introns

BACKGROUND

Group I introns are one of the four major classes of introns as defined by their distinct splicing mechanisms. Because they catalyze their own removal from precursor transcripts, group I introns are referred to as autocatalytic introns. Group I introns are common in fungal and protist nuclear ribosomal RNA genes and in organellar genomes. In contrast, they are rare in all other organisms and genomes, including bacteria.

RESULTS

Here we report five group I introns, each containing a LAGLIDADG homing endonuclease gene (HEG), in large subunit (LSU) rRNA genes of cyanobacteria. Three of the introns are located in the LSU gene of Synechococcus sp. C9, and the other two are in the LSU gene of Synechococcus lividus strain C1. Phylogenetic analyses show that these introns and their HEGs are closely related to introns and HEGs located at homologous insertion sites in organellar and bacterial rDNA genes. We also present a compilation of group I introns with homing endonuclease genes in bacteria.

CONCLUSION

We have discovered multiple HEG-containing group I introns in a single bacterial gene. To our knowledge, these are the first cases of multiple group I introns in the same bacterial gene (multiple group I introns have been reported in at least one phage gene and one prophage gene). The HEGs each contain one copy of the LAGLIDADG motif and presumably function as homodimers. Phylogenetic analysis, in conjunction with their patchy taxonomic distribution, suggests that these intron-HEG elements have been transferred horizontally among organelles and bacteria. However, the mode of transfer and the nature of the biological connections among the intron-containing organisms are unknown.

FR3D: finding local and composite recurrent structural motifs in RNA 3D structures

New methods are described for finding recurrent three-dimensional (3D) motifs in RNA atomic-resolution structures. Recurrent RNA 3D motifs are sets of RNA nucleotides with similar spatial arrangements. They can be local or composite. Local motifs comprise nucleotides that occur in the same hairpin or internal loop. Composite motifs comprise nucleotides belonging to three or more different RNA strand segments or molecules. We use a base-centered approach to construct efficient, yet exhaustive search procedures using geometric, symbolic, or mixed representations of RNA structure that we implement in a suite of MATLAB programs, "Find RNA 3D" (FR3D). The first modules of FR3D preprocess structure files to classify base-pair and -stacking interactions. Each base is represented geometrically by the position of its glycosidic nitrogen in 3D space and by the rotation matrix that describes its orientation with respect to a common frame. Base-pairing and base-stacking interactions are calculated from the base geometries and are represented symbolically according to the Leontis/Westhof basepairing classification, extended to include base-stacking. These data are stored and used to organize motif searches. For geometric searches, the user supplies the 3D structure of a query motif which FR3D uses to find and score geometrically similar candidate motifs, without regard to the sequential position of their nucleotides in the RNA chain or the identity of their bases. To score and rank candidate motifs, FR3D calculates a geometric discrepancy by rigidly rotating candidates to align optimally with the query motif and then comparing the relative orientations of the corresponding bases in the query and candidate motifs. Given the growing size of the RNA structure database, it is impossible to explicitly compute the discrepancy for all conceivable candidate motifs, even for motifs with less than ten nucleotides. The screening algorithm that we describe finds all candidate motifs whose geometric discrepancy with respect to the query motif falls below a user-specified cutoff discrepancy. This technique can be applied to RMSD searches. Candidate motifs identified geometrically may be further screened symbolically to identify those that contain particular basepair types or base-stacking arrangements or that conform to sequence continuity or nucleotide identity constraints. Purely symbolic searches for motifs containing user-defined sequence, continuity and interaction constraints have also been implemented. We demonstrate that FR3D finds all occurrences, both local and composite and with nucleotide substitutions, of sarcin/ricin and kink-turn motifs in the 23S and 5S ribosomal RNA 3D structures of the H. marismortui 50S ribosomal subunit and assigns the lowest discrepancy scores to bona fide examples of these motifs. The search algorithms have been optimized for speed to allow users to search the non-redundant RNA 3D structure database on a personal computer in a matter of minutes.

Structural metals in the group I intron: a ribozyme with a multiple metal ion core

Metal ions play key roles in the folding and function for many structured RNAs, including group I introns. We determined the X-ray crystal structure of the Azoarcus bacterial group I intron in complex with its 5' and 3' exons. In addition to 222 nucleotides of RNA, the model includes 18 Mg(2+) and K(+) ions. Five of the metals bind within 12 A of the scissile phosphate and coordinate the majority of the oxygen atoms biochemically implicated in conserved metal-RNA interactions. The metals are buried deep within the structure and form a multiple metal ion core that is critical to group I intron structure and function. Eight metal ions bind in other conserved regions of the intron structure, and the remaining five interact with peripheral structural elements. Each of the 18 metals mediates tertiary interactions, facilitates local bends in the sugar-phosphate backbone or binds in the major groove of helices. The group I intron has a rich history of biochemical efforts aimed to identify RNA-metal ion interactions. The structural data are correlated to the biochemical results to further understand the role of metal ions in group I intron structure and function.

Ribozymes

The structural molecular biology of ribozymes took another great leap forward during the past two years. Before ribozymes were discovered in the early 1980s, all enzymes were thought to be proteins. No detailed structural information on ribozymes became available until 1994. Now, within the past two years, near atomic resolution crystal structures are available for almost all of the known ribozymes. The latest additions include ribonuclease P, group I intron structures, the ribosome (the peptidyl transferase appears to be a ribozyme) and several smaller ribozymes, including a Diels-Alderase, the glmS ribozyme and a new hammerhead ribozyme structure that reconciles 12 years of discord. Although not all ribozymes are metalloenzymes, acid-base catalysis appears to be a universal property shared by all ribozymes as well as many of their protein cousins.

A Second Divalent Metal Ion in the Group II Intron Reaction Center

Group II introns are mobile genetic elements that have been implicated as agents of genetic diversity, and serve as important model systems for investigating RNA catalysis and pre-mRNA splicing. In the absence of an atomic-resolution structure of the intron, detailed understanding of its catalytic mechanism has remained elusive. Previous identification of a divalent metal ion stabilizing the leaving group in both splicing steps suggested that the group II intron may employ a "two-metal ion" mechanism, a catalytic strategy used by a number of protein phosphoester transfer enzymes. Using metal rescue experiments, we now reveal the presence of a second metal ion required for nucleophile activation in the exon-ligation step of group II intron splicing. Coupled with biochemical and structural evidence of at least two metal ions at the group I intron reaction center, these results suggest a mechanistic paradigm for describing catalysis by large ribozymes.

Thermodynamics and kinetics for base-pair opening in the P1 duplex of the Tetrahymena group I ribozyme

The thermodynamics and kinetics for base-pair opening of the P1 duplex of the Tetrahymena group I ribozyme were studied by NMR hydrogen exchange experiments. The apparent equilibrium constants for base pair opening were measured for most of the imino protons in the P1 duplex using the base catalysts NH₃, HPO42− or TRIS. These equilibrium constants were also measured for several modified P1 duplexes, and the C-2·G23 base pair was the most stable base pair in all the duplexes. The conserved U-1·G22 base pair is required for activity of the ribozyme and the data here show that this wobble base pair destabilizes neighboring base pairs on only one side of the wobble. A 2′-OMe modification on the U-3 residue stabilized its own base pair but had little effect on the neighboring base pairs. Three base pairs, U-1·G22, C-2·G23 and A2·U21 showed unusual equilibrium constants for opening and possible implications of the opening thermodynamics of these base pairs on the undocking rates of the P1 helix with catalytic core are discussed.

Local RNA structural changes induced by crystallization are revealed by SHAPE

We present a simple approach to locate sites that undergo conformational changes upon crystallization by comparative structural mapping of the same RNA in three different environments. As a proof of principle, we probed the readily crystallized P4-P6DeltaC209 domain from the Tetrahymena thermophila group I intron in a native solution, in a solution mimicking the crystallization drop, and in crystals. We chose the selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry, which monitors the flexibility and the conformation of each nucleotide. First, SHAPE successfully revealed the structural changes that occur during the crystallization process. Specifically, 64% of the nucleotides implicated in packing contacts and present in the portion of the molecule analyzed were identified. Second, reactivity differences for some of these nucleotides were already observed in the crystallization solution, suggesting that the crystallization buffer locked down a particular structure that was favorable to crystal formation. Third, the probing of a known structure extends our understanding of the structural basis for the SHAPE reaction by suggesting that reactivity is enhanced by a C2'-endo sugar pucker. Furthermore, by identifying local conformational changes of the RNA that take place during crystallization, SHAPE could be combined with the in vitro selection of stable mutants to rationalize the design of RNA candidates for crystallization.

RNABC: forward kinematics to reduce all-atom steric clashes in RNA backbone

Although accurate details in RNA structure are of great importance for understanding RNA function, the backbone conformation is difficult to determine, and most existing RNA structures show serious steric clashes (>or= 0.4 A overlap) when hydrogen atoms are taken into account. We have developed a program called RNABC (RNA Backbone Correction) that performs local perturbations to search for alternative conformations that avoid those steric clashes or other local geometry problems. Its input is an all-atom coordinate file for an RNA crystal structure (usually from the MolProbity web service), with problem areas specified. RNABC rebuilds a suite (the unit from sugar to sugar) by anchoring the phosphorus and base positions, which are clearest in crystallographic electron density, and reconstructing the other atoms using forward kinematics. Geometric parameters are constrained within user-specified tolerance of canonical or original values, and torsion angles are constrained to ranges defined through empirical database analyses. Several optimizations reduce the time required to search the many possible conformations. The output results are clustered and presented to the user, who can choose whether to accept one of the alternative conformations. Two test evaluations show the effectiveness of RNABC, first on the S-motifs from 42 RNA structures, and second on the worst problem suites (clusters of bad clashes, or serious sugar pucker outliers) in 25 unrelated RNA structures. Among the 101 S-motifs, 88 had diagnosed problems, and RNABC produced clash-free conformations with acceptable geometry for 71 of those (about 80%). For the 154 worst problem suites, RNABC proposed alternative conformations for 72. All but 8 of those were judged acceptable after examining electron density (where available) and local conformation. Thus, even for these worst cases, nearly half the time RNABC suggested corrections suitable to initiate further crystallographic refinement. The program is available from http://kinemage.biochem.duke.edu .

Insight into the functional versatility of RNA through model-making with applications to data fitting

The RNA World hypothesis posits that life emerged from self-replicating RNA molecules. For any single biopolymer to be the basis for life, it must both store information and perform diverse functions. It is well known that RNA can store information. Advances in recent years have revealed that RNA can exhibit remarkable functional versatility as well. In an effort to judge the functional versatility of RNA and thereby the plausibility that RNA was at one point the basis for life, a statistical chemical approach is adopted. Essential biological functions are reduced to simple molecular models in a minimalist, biopolymer-independent fashion. The models dictate requisite states, populations of states, and physical and chemical changes occurring between the states. Next, equations are derived from the models, which lead to complex phenomenological constants such as observed and functional constants that are defined in terms of familiar elementary chemical descriptors: intrinsic rate constants, microscopic ligand equilibrium constants, secondary structure stability, and ligand concentration. Using these equations, simulations of functional behavior are performed. These functional models provide practical frameworks for fitting and organizing real data on functional RNAs such as ribozymes and riboswitches. At the same time, the models allow the suitability of RNA as a basis for life to be judged. We conclude that RNA, while inferior to extant proteins in most, but not all, functional respects, may be more versatile than proteins, performing a wider range of elementary biological functions at a tolerable level. Inspection of the functional models and various RNA structures uncovers several surprising ways in which the nucleobases can conspire to afford chemical catalysis and evolvability. These models support the plausibility that RNA, or a closely related informational biopolymer, could serve as the basis for a fairly simple form of life.

An allosteric-feedback mechanism for protein-assisted group I intron splicing

The I-AniI maturase facilitates self-splicing of a mitochondrial group I intron in Aspergillus nidulans. Binding occurs in at least two steps: first, a specific but labile encounter complex rapidly forms and then this intermediate is slowly resolved into a native, catalytically active RNA/protein complex. Here we probe the structure of the RNA throughout the assembly pathway. Although inherently unstable, the intron core, when bound by I-AniI, undergoes rapid folding to a near-native state in the encounter complex. The next transition includes the slow destabilization and docking into the core of the peripheral stacked helix that contains the 5' splice site. Mutational analyses confirm that both transitions are important for native complex formation. We propose that protein-driven destabilization and docking of the peripheral stacked helix lead to subtle changes in the I-AniI binding site that facilitate native complex formation. These results support an allosteric-feedback mechanism of RNA-protein recognition in which proteins engaged in an intermediate complex can influence RNA structure far from their binding sites. The linkage of these changes to stable binding ensures that the protein and RNA do not get sequestered in nonfunctional complexes.

Structural protein analysis of the polyvalent staphylococcal bacteriophage 812

Phage 812 is a polyvalent phage with a very broad host range in the genus Staphylococcus, which makes it a suitable candidate for phage therapy of staphylococcal infections. This proteomic study, combining the results of both 1-DE and 2-DE followed by PMF, led to the identification of 24 virion proteins. Twenty new proteins, not yet identified by proteome analysis of closely related staphylococcal phages K and G1 were identified using this approach. Fifteen proteins were assigned unambiguously to the head-tail genome module; the remaining nine proteins are encoded by genes of the left or right arms of the phage genome. As expected, the most abundant proteins in the electrophoretic patterns are the major capsid protein, the major tail sheath protein and proteins identical to ORF 50 and ORF 95 of phage K, although their function is only putative. Identification of these 20 new proteins contributes substantially to a detailed characterization of phage virions, knowledge of which is necessary for rational phage therapy.

The interaction networks of structured RNAs

All pairwise interactions occurring between bases which could be detected in three-dimensional structures of crystallized RNA molecules are annotated on new planar diagrams. The diagrams attempt to map the underlying complex networks of base–base interactions and, especially, they aim at conveying key relationships between helical domains: co-axial stacking, bending and all Watson–Crick as well as non-Watson–Crick base pairs. Although such wiring diagrams cannot replace full stereographic images for correct spatial understanding and representation, they reveal structural similarities as well as the conserved patterns and distances between motifs which are present within the interaction networks of folded RNAs of similar or unrelated functions. Finally, the diagrams could help devising methods for meaningfully transforming RNA structures into graphs amenable to network analysis.

Unusual metal specificity and structure of the group I ribozyme from Chlamydomonas reinhardtii 23S rRNA

Group I intron ribozymes require cations for folding and catalysis, and the current literature indicates that a number of cations can promote folding, but only Mg2+ and Mn2+ support both processes. However, some group I introns are active only with Mg2+, e.g. three of the five group I introns in Chlamydomonas reinhardtii. We have investigated one of these ribozymes, an intron from the 23S LSU rRNA gene of Chlamydomonas reinhardtii (Cr.LSU), by determining if the inhibition by Mn2+ involves catalysis, folding, or both. Kinetic analysis of guanosine-dependent cleavage by a Cr.LSU ribozyme, 23S.5 Delta Gb, that lacks the 3' exon and intron-terminal G shows that Mn2+ does not affect guanosine binding or catalysis, but instead promotes misfolding of the ribozyme. Surprisingly, ribozyme misfolding induced by Mn2+ is highly cooperative, with a Hill coefficient larger than that of native folding induced by Mg2+. At lower Mn2+ concentrations, metal inhibition is largely alleviated by the guanosine cosubstrate (GMP). The concentration dependence of guanosine cosubstrate-induced folding suggests that it functions by interacting with the G binding site, perhaps by displacing an inhibitory Mn2+. Because of these and other properties of Cr.LSU, the tertiary structure of the intron from 23S.5 Delta Gb was examined using Fe2+-EDTA cleavage. The ground-state structure shows evidence of an unusually open ribozyme core: the catalytic P3-P7 domain and the nucleotides that connect it to the P4-P5-P6 domain are exposed to solvent. The implications of this structure for the in vitro and in vivo properties of this intron ribozyme are discussed.

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

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

Structure of the S-adenosylmethionine riboswitch regulatory mRNA element

Riboswitches are cis-acting genetic regulatory elements found in the 5'-untranslated regions of messenger RNAs that control gene expression through their ability to bind small molecule metabolites directly. Regulation occurs through the interplay of two domains of the RNA: an aptamer domain that responds to intracellular metabolite concentrations and an expression platform that uses two mutually exclusive secondary structures to direct a decision-making process. In Gram-positive bacteria such as Bacillus species, riboswitches control the expression of more than 2% of all genes through their ability to respond to a diverse set of metabolites including amino acids, nucleobases and protein cofactors. Here we report the 2.9-angstroms resolution crystal structure of an S-adenosylmethionine (SAM)-responsive riboswitch from Thermoanaerobacter tengcongensis complexed with S-adenosylmethionine, an RNA element that controls the expression of several genes involved in sulphur and methionine metabolism. This RNA folds into a complex three-dimensional architecture that recognizes almost every functional group of the ligand through a combination of direct and indirect readout mechanisms. Ligand binding induces the formation of a series of tertiary interactions with one of the helices, serving as a communication link between the aptamer and expression platform domains.

The ARTS web server for aligning RNA tertiary structures

RNA molecules with common structural features may share similar functional properties. Structural comparison of RNAs and detection of common substructures is, thus, a highly important task. Nevertheless, the current available tools in the RNA community provide only a partial solution, since they either work at the 2D level or are suitable for detecting predefined or local contiguous tertiary motifs only. Here, we describe a web server built around ARTS, a method for aligning tertiary structures of nucleic acids (both RNA and DNA). ARTS receives a pair of 3D nucleic acid structures and searches for a priori unknown common substructures. The search is truly 3D and irrespective of the order of the nucleotides on the chain. The identified common substructures can be large global folds with hundreds and even thousands of nucleotides as well as small local motifs with at least two successive base pairs. The method is highly efficient and has been used to conduct an all-against-all comparison of all the RNA structures in the Protein Data Bank. The web server together with a software package for download are freely accessible at http://bioinfo3d.cs.tau.ac.il/ARTS.

Mg2+ mimicry in the promotion of group I ribozyme activities by aminoglycoside antibiotics

Aminoglycoside antibiotics inhibit several types of ribozymes, including group I introns, by displacing critical Mg2+ ions. However, they stimulate activity of the small hairpin ribozyme. We show here that aminoglycosides promote self-splicing of the Cr.psbA2 group I intron at subthreshold Mg2+ concentrations. Neomycin is the most effective of the aminoglycosides tested; it stimulates splicing of Cr.psbA2 at micromolar concentrations, and, in this respect, is >100-fold more effective than spermidine. At optimal Mg2+ for Cr.psbA2 splicing, these drugs, especially kanamycin B and tobramycin, promote GTP attack at the 3' splice-site. Kinetic analysis suggests that this is due to an alternatively folded state of the ribozyme that is induced, or stabilized, by aminoglycosides. A similar effect is observed at high Mg2+ concentrations. Comparing the effects of structurally related aminoglycosides indicates that splicing promotion is more sensitive to drug structure than misfolding and occurs at lower drug concentrations. These data show that aminoglycosides can promote biochemical activities of a large ribozyme by acting as a Mg2+ mimic. The results also underscore the functional diversity of group I introns in nature.

RNA splicing: group I intron crystal structures reveal the basis of splice site selection and metal ion catalysis

The group I intron has served as a model for RNA catalysis since its discovery 25 years ago. Four recently determined high-resolution crystal structures complement extensive biochemical studies on this system. Structures of the Azoarcus, Tetrahymena and bacteriophage Twort group I introns mimic different states of the splicing or ribozyme reaction pathway and provide information on splice site selection and metal ion catalysis. The 5'-splice site is selected by formation of a conserved G.U wobble pair between the 5'-exon terminus and the intron. The 3'-splice site is identified through stacking of three base triples, in which the middle triple contains the conserved terminal nucleotide of the intron, OmegaG. The structures support a two-metal-ion mechanism for group I intron splicing that might have corollaries to group II intron and pre-mRNA splicing by the spliceosome.

Structure of ribonuclease P--a universal ribozyme

Ribonuclease P (RNase P) is one of only two known universal ribozymes and was one of the first ribozymes to be discovered. It is involved in RNA processing, in particular the 5' maturation of tRNA. Unlike most other natural ribozymes, it recognizes and cleaves its substrate in trans. RNase P is a ribonucleoprotein complex containing one RNA subunit and as few as one protein subunit. It has been shown that, in bacteria and in some archaea, the RNA subunit alone can support catalysis. The structure and function of bacterial RNase P RNA have been studied extensively, but the detailed catalytic mechanism is not yet fully understood. Recently, structures of one of the structural domains and of the entire RNA component of RNase P from two different bacteria have been described. These structures provide the first atomic-level information on the structural assembly of the RNA component, and the regions involved in substrate recognition and catalysis. Comparison of these structures reveals a highly conserved core that comprises two universally conserved structural modules. Interestingly, the same structural core can be found in the context of different scaffolds.

Monovalent cations use multiple mechanisms to resolve ribozyme misfolding

Recent efforts have been made to unravel the independent roles of monovalent cations in RNA folding, primarily using the Tetrahymena ribozyme as a model. Here we report how monovalent cations impact the folding of the Candida ribozyme. Interestingly, this ribozyme requires an order of magnitude less monovalent cations (Na+ and Tris+) to commit to a new folding starting state in which the J3/4:P6 base triple is partially formed and mispairing in the L2.1 and L6 terminal loops is resolved. When Mg2+-induced ribozyme folding proceeded on the same energy landscape, the altered starting state led to a rapid assembly of the correct ribozyme core and a fivefold to 10-fold increase in the ribozyme activity. Moreover, when the ribozyme folding was started from a misfolding-prone state, high millimolar concentrations of monovalent cations moderately elevated the ribozyme activity by efficiently resolving the misfolding of a peripheral element, P5abc.

A model system for investigating lineshape/structure correlations in RNA site-directed spin labeling

In RNA site-directed spin labeling (SDSL) studies, structural and dynamic information at the individual RNA nucleotide level is derived from the observed electron paramagnetic resonance spectrum of a covalently attached nitroxide. A systematic approach for RNA SDSL is to establish a library that categorizes observed spectral lineshapes based on known RNA structures, thus enabling lineshape-based structure identification at any RNA site. To establish the first RNA SDSL library, selective secondary structure elements have been systematically engineered into a model RNA. Nitroxide lineshapes reporting features specific to each element were obtained utilizing a new avidin-tethering scheme for suppressing spectral effects due to uniform RNA tumbling. The data demonstrated two key features required for a SDSL library with a predicting power: (i) spectral divergence--distinctive lineshape for different elements; and (ii) spectral convergence--similar lineshape for the same element in different contexts. This sets the foundation for further RNA SDSL library development.

Comparison of crystal structure interactions and thermodynamics for stabilizing mutations in the Tetrahymena ribozyme

Although general mechanisms of RNA folding and catalysis have been elucidated, little is known about how ribozymes achieve structural stability at high temperature. A previous in vitro evolution experiment identified a small number of mutations that significantly increase the thermostability of the tertiary structure of the Tetrahymena ribozyme. Because we also determined the crystal structure of this thermostable ribozyme, we have for the first time the opportunity to compare the structural interactions and thermodynamic contributions of individual nucleotides in a ribozyme. We investigated the contribution of five mutations to thermostability by using temperature gradient gel electrophoresis. Unlike the case with several well-studied proteins, the effects of individual mutations on thermostability of this RNA were highly context dependent. The three most important mutations for thermostability were actually destabilizing in the wild-type background. A269G and A304G contributed to stability only when present as a pair, consistent with their proximity in the ribozyme structure. In an evolutionary context, this work supports and extends the idea that one advantage of protein enzyme systems over an RNA world is the ability of proteins to accumulate stabilizing single-site mutations, whereas RNA may often require much rarer double mutations to improve the stability of both its tertiary and secondary structures.

Nucleic acid visualization with UCSF Chimera

With the increase in the number of large, 3D, high-resolution nucleic acid structures, particularly of the 30S and 50S ribosomal subunits and the intact bacterial ribosome, advancements in the visualization of nucleic acid structural features are essential. Large molecular structures are complicated and detailed, and one goal of visualization software is to allow the user to simplify the display of some features and accent others. We describe an extension to the UCSF Chimera molecular visualization system for the purpose of displaying and highlighting nucleic acid characteristics, including a new representation of sugar pucker, several options for abstraction of base geometries that emphasize stacking and base pairing, and an adaptation of the ribbon backbone to accommodate the nucleic acid backbone. Molecules are displayed and manipulated interactively, allowing the user to change the representations as desired for small molecules, proteins and nucleic acids. This software is available as part of the UCSF Chimera molecular visualization system and thus is integrated with a suite of existing tools for molecular graphics.

Topology of three-way junctions in folded RNAs

The three-way junctions contained in X-ray structures of folded RNAs have been compiled and analyzed. Three-way junctions with two helices approximately coaxially stacked can be divided into three main families depending on the relative lengths of the segments linking the three Watson-Crick helices. Each family has topological characteristics with some conservation in the non-Watson-Crick pairs within the linking segments as well as in the types of contacts between the segments and the helices. The most populated family presents tertiary interactions between two helices as well as extensive shallow/minor groove contacts between a linking segment and the third helix. On the basis of the lengths of the linking segments, some guidelines could be deduced for choosing a topology for a three-way junction on the basis of a secondary structure. Examples and prediction bas'ed on those rules are discussed.

Atomic level architecture of group I introns revealed

Twenty-two years after their discovery as ribozymes, the self-splicing group I introns are finally disclosing their architecture at the atomic level. The crystal structures of three group I introns solved at moderately high resolution (3.1-3.8A) reveal a remarkably conserved catalytic core bound to the metal ions required for activity. The structure of the core is stabilized by an intron-specific set of long-range interactions that involves peripheral elements. Group I intron structures thus provide much awaited and extremely valuable snapshots of how these ribozymes coordinate substrate binding and catalysis.

Structure, folding and mechanisms of ribozymes

The past two years have seen exciting developments in RNA catalysis. A completely new ribozyme (possibly two) has come along and several new structures have been determined, including three different group I intron species. Although the origins of catalysis remain incompletely understood, a significant convergence of views has happened in the past year, together with the discovery of new super-fast ribozymes. There is persuasive evidence of general acid-base chemistry in nucleolytic ribozymes, whereas catalysis of peptidyl transfer in the ribosome seems to result largely from orientation and proximity effects. Lastly, important new folding-enhancing elements have been discovered.

RNA structure: the long and the short of it

The database of RNA structure has grown tremendously since the crystal structure analyses of ribosomal subunits in 2000–2001. During the past year, the trend toward determining the structure of large, complex biological RNAs has accelerated, with the analysis of three intact group I introns, A- and B-type ribonuclease P RNAs, a riboswitch–substrate complex and other structures. The growing database of RNA structures, coupled with efforts directed at the standardization of nomenclature and classification of motifs, has resulted in the identification and characterization of numerous RNA secondary and tertiary structure motifs. Because a large proportion of RNA structure can now be shown to be composed of these recurring structural motifs, a view of RNA as a modular structure built from a combination of these building blocks and tertiary linkers is beginning to emerge. At the same time, however, more detailed analysis of water, metal, ligand and protein binding to RNA is revealing the effect of these moieties on folding and structure formation. The balance between the views of RNA structure either as strictly a construct of preformed building blocks linked in a limited number of ways or as a flexible polymer assuming a global fold influenced by its environment will be the focus of current and future RNA structural biology.

Structural Evidence for a Two-Metal-Ion Mechanism of Group I Intron Splicing

We report the 3.4 angstrom crystal structure of a catalytically active group I intron splicing intermediate containing the complete intron, both exons, the scissile phosphate, and all of the functional groups implicated in catalytic metal ion coordination, including the 2'-OH of the terminal guanosine. This structure suggests that, like protein phosphoryltransferases, an RNA phosphoryltransferase can use a two-metal-ion mechanism. Two Mg2+ ions are positioned 3.9 angstroms apart and are directly coordinated by all six of the biochemically predicted ligands. The evolutionary convergence of RNA and protein active sites on the same inorganic architecture highlights the intrinsic chemical capacity of the two-metal-ion catalytic mechanism for phosphoryl transfer.

RNA Structure: Reading the Ribosome

The crystal structures of the ribosome and its subunits have increased the amount of information about RNA structure by about two orders of magnitude. This is leading to an understanding of the principles of RNA folding and of the molecular interactions that underlie the functional capabilities of the ribosome and other RNA systems. Nearly all of the possible types of RNA tertiary interactions have been found in ribosomal RNA. One of these, an abundant tertiary structural motif called the A-minor interaction, has been shown to participate in both aminoacyl-transfer RNA selection and in peptidyl transferase; it may also play an important role in the structural dynamics of the ribosome.

Site-directed spin labeling studies reveal solution conformational changes in a GAAA tetraloop receptor upon Mg(2+)-dependent docking of a GAAA tetraloop

The Mg(2+)-dependent GAAA tetraloop interaction with its 11 nucleotide receptor is one of the most frequently occurring long-range tertiary interactions in RNAs. To explore conformational changes in the receptor during tetraloop docking, nitroxide spin labels were attached at each of four uridine bases, one at a time, within an RNA molecule containing the receptor sequence. In the presence of Mg2+ and the tetraloop, the electron paramagnetic resonance (EPR) spectrum of one of the labeled bases reflected a large increase in mobility, indicating unstacking of the base upon tetraloop docking. This provides direct evidence that base unstacking is an intrinsic feature of the solution tetraloop-receptor complex formed in the presence of Mg2+. Additional evidence suggests that in solution the bound receptor conformation is similar to that observed in the crystal structure of a group I intron ribozyme domain. In Mg2+ alone, a receptor conformation with an unstacked base was not detectable, suggesting that this conformation is of higher standard state free energy than that of the free receptor. This leads to the conclusion that the extensive RNA-RNA interactions observed in the crystal structure of the tetraloop-receptor complex provide larger interaction energy than the measured apparent affinity between the tetraloop and the free receptor. This is compatible with a high specificity of the tetraloop-receptor interaction.

Evolution from DNA to RNA recognition by the bI3 LAGLIDADG maturase

LAGLIDADG endonucleases bind across adjacent major grooves via a saddle-shaped surface and catalyze DNA cleavage. Some LAGLIDADG proteins, called maturases, facilitate splicing by group I introns, raising the issue of how a DNA-binding protein and an RNA have evolved to function together. In this report, crystallographic analysis shows that the global architecture of the bI3 maturase is unchanged from its DNA-binding homologs; in contrast, the endonuclease active site, dispensable for splicing facilitation, is efficiently compromised by a lysine residue replacing essential catalytic groups. Biochemical experiments show that the maturase binds a peripheral RNA domain 50 A from the splicing active site, exemplifying long-distance structural communication in a ribonucleoprotein complex. The bI3 maturase nucleic acid recognition saddle interacts at the RNA minor groove; thus, evolution from DNA to RNA function has been mediated by a switch from major to minor groove interaction.

Functional identification of catalytic metal ion binding sites within RNA

The viability of living systems depends inextricably on enzymes that catalyze phosphoryl transfer reactions. For many enzymes in this class, including several ribozymes, divalent metal ions serve as obligate cofactors. Understanding how metal ions mediate catalysis requires elucidation of metal ion interactions with both the enzyme and the substrate(s). In the Tetrahymena group I intron, previous work using atomic mutagenesis and quantitative analysis of metal ion rescue behavior identified three metal ions (M_A, M_B, and M_C) that make five interactions with the ribozyme substrates in the reaction's transition state. Here, we combine substrate atomic mutagenesis with site-specific phosphorothioate substitutions in the ribozyme backbone to develop a powerful, general strategy for defining the ligands of catalytic metal ions within RNA. In applying this strategy to the Tetrahymena group I intron, we have identified the pro-S_P phosphoryl oxygen at nucleotide C262 as a ribozyme ligand for M_C. Our findings establish a direct connection between the ribozyme core and the functionally defined model of the chemical transition state, thereby extending the known set of transition-state interactions and providing information critical for the application of the recent group I intron crystallographic structures to the understanding of catalysis.

A combination of substrate atomic mutagenesis with site-specific substitutions in the ribozyme backbone allow the ligands of catalytic metal ions to be identified.

A peripheral element assembles the compact core structure essential for group I intron self-splicing

The presence of non-conserved peripheral elements in all naturally occurring group I introns underline their importance in ensuring the natural intron function. Recently, we reported that some peripheral elements are conserved in group I introns of IE subgroup. Using self-splicing activity as a readout, our initial screening revealed that one such conserved peripheral elements, P2.1, is mainly required to fold the catalytically active structure of the Candida ribozyme, an IE intron. Unexpectedly, the essential function of P2.1 resides in a sequence-conserved short stem of P2.1 but not in a long-range interaction associated with the loop of P2.1 that stabilizes the ribozyme structure. The P2.1 stem is indispensable in folding the compact ribozyme core, most probably by forming a triple helical interaction with two core helices, P3 and P6. Surprisingly, although the ribozyme lacking the P2.1 stem renders a loosely folded core and the loss of self-splicing activity requires two consecutive transesterifications, the mutant ribozyme efficiently catalyzes the first transesterification reaction. These results suggest that the intron self-splicing demands much more ordered structure than does one independent transesterification, highlighting that the universally present peripheral elements achieve their functional importance by enabling the highly ordered structure through diverse tertiary interactions.

The catalytic diversity of RNAs

The natural RNA enzymes catalyse phosphate-group transfer and peptide-bond formation. Initially, metal ions were proposed to supply the chemical versatility that nucleotides lack. In the ensuing decades, structural and mechanistic studies have substantially altered this initial viewpoint. Whereas self-splicing ribozymes clearly rely on essential metal-ion cofactors, self-cleaving ribozymes seem to use nucleotide bases for their catalytic chemistry. Despite the overall differences in chemical features, both RNA and protein enzymes use similar catalytic strategies.

Ribozyme catalysis: not different, just worse

Evolution has resoundingly favored protein enzymes over RNA-based catalysts, yet ribozymes occupy important niches in modern cell biology that include the starring role in catalysis of protein synthesis on the ribosome. Recent results from structural and biochemical studies show that natural ribozymes use an impressive range of catalytic mechanisms, beyond metalloenzyme chemistry and analogous to more chemically diverse protein enzymes. These findings make it increasingly possible to compare details of RNA- and protein-based catalysis.

Structure and assembly of group I introns

Self-splicing group I introns have served as a model for RNA catalysis and folding for over two decades. New three-dimensional structures now bring the details into view. Revelations include an unanticipated turn in the RNA backbone around the guanosine-binding pocket. Two metal ions in the active site coordinate the substrate and phosphates from all three helical domains.

Linkage between proton binding and folding in RNA: implications for RNA catalysis

Small ribozymes use their nucleobases to catalyse phosphodiester bond cleavage. The hepatitis delta virus ribozyme employs C75 as a general acid to protonate the 5′-bridging oxygen leaving group, and to accomplish this task efficiently, it shifts its pKa towards neutrality. Simulations and thermodynamic experiments implicate linkage between folding and protonation in nucleobase pKa shifting. Even small oligonucleotides are shown to fold in a highly co-operative manner, although they do so in a context-specific fashion. Linkage between protonation and co-operativity of folding may drive pKa shifting and provide for enhanced function in RNA.

Predicting the secondary structures and tertiary interactions of 211 group I introns in IE subgroup

The large number of currently available group I intron sequences in the public databases provides opportunity for studying this large family of structurally complex catalytic RNA by large-scale comparative sequence analysis. In this study, the detailed secondary structures of 211 group I introns in the IE subgroup were manually predicted. The secondary structure-favored alignments showed that IE introns contain 14 conserved stems. The P13 stem formed by long-range base-pairing between P2.1 and P9.1 is conserved among IE introns. Sequence variations in the conserved core divide IE introns into three distinct minor subgroups, namely IE1, IE2 and IE3. Co-variation of the peripheral structural motifs with core sequences supports that the peripheral elements function in assisting the core structure folding. Interestingly, host-specific structural motifs were found in IE2 introns inserted at S516 position. Competitive base-pairing is found to be conserved at the junctions of all long-range paired regions, suggesting a possible mechanism of establishing long-range base-pairing during large RNA folding. These findings extend our knowledge of IE introns, indicating that comparative analysis can be a very good complement for deepening our understanding of RNA structure and function in the genomic era.