Modeling RNA tertiary structure motifs by graph-grammars

A new approach, graph-grammars, to encode RNA tertiary structure patterns is introduced and exemplified with the classical sarcin-ricin motif. The sarcin-ricin motif is found in the stem of the crucial ribosomal loop E (also referred to as the sarcin-ricin loop), which is sensitive to the alpha-sarcin and ricin toxins. Here, we generate a graph-grammar for the sarcin-ricin motif and apply it to derive putative sequences that would fold in this motif. The biological relevance of the derived sequences is confirmed by a comparison with those found in known sarcin-ricin sites in an alignment of over 800 bacterial 23S ribosomal RNAs. The comparison raised alternative alignments in few sarcin-ricin sites, which were assessed using tertiary structure predictions and 3D modeling. The sarcin-ricin motif graph-grammar was built with indivisible nucleotide interaction cycles that were recently observed in structured RNAs. A comparison of the sequences and 3D structures of each cycle that constitute the sarcin-ricin motif gave us additional insights about RNA sequence-structure relationships. In particular, this analysis revealed the sequence space of an RNA motif depends on a structural context that goes beyond the single base pairing and base-stacking interactions.     

Enhanced base-pair opening in the adenine tract of a RNA double helix

Proton exchange and nuclear magnetic resonance spectroscopy are being used to characterize the kinetics and energetics of base-pair opening in two nucleic acid double helices. One is the RNA duplex 5'-r(GCGAUAAAAAGGCC)-3'/5'-r(GGCCUUUUUAUCGC)-3', which contains a central tract of five AU base pairs. The other is the homologous DNA duplex with a central tract of five AT base pairs. The rates and the equilibrium constants of the opening reaction of each base pair are measured from the dependence of the exchange rates of imino protons on ammonia concentration, at 10 °C. The results reveal that the tract of AU base pairs in the RNA duplex differs from the homologous tract of AT base pairs in DNA in several ways. The rates of opening of AU base pairs in RNA are high and increase progressively along the tract, reaching their largest values at the 3'-end of the tract. In contrast, the opening rates of AT base pairs in DNA are much lower than those of AU base pairs. Within the tract, the largest opening rate is observed for the AT base pair at the 5'-end of the tract. These differences in opening kinetics are paralleled by differences in the stabilities of individual base pairs. All AU base pairs in the RNA are less stable than the AT base pairs in the DNA. The presence of the tract enhances these differences by increasing the stability of AT base pairs in DNA while decreasing the stability of AU base pairs in RNA. Due to these divergent trends, along the tracts, the AU base pairs become progressively less stable than AT base pairs. These findings demonstrate that tracts of AU base pairs in RNA have specific dynamic and energetic signatures that distinguish them from similar tracts of AT base pairs in DNA.     

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.     

RNAMotifScan: automatic identification of RNA structural motifs using secondary structural alignment

Recent studies have shown that RNA structural motifs play essential roles in RNA folding and interaction with other molecules. Computational identification and analysis of RNA structural motifs remains a challenging task. Existing motif identification methods based on 3D structure may not properly compare motifs with high structural variations. Other structural motif identification methods consider only nested canonical base-pairing structures and cannot be used to identify complex RNA structural motifs that often consist of various non-canonical base pairs due to uncommon hydrogen bond interactions. In this article, we present a novel RNA structural alignment method for RNA structural motif identification, RNAMotifScan, which takes into consideration the isosteric (both canonical and non-canonical) base pairs and multi-pairings in RNA structural motifs. The utility and accuracy of RNAMotifScan is demonstrated by searching for kink-turn, C-loop, sarcin-ricin, reverse kink-turn and E-loop motifs against a 23S rRNA (PDBid: 1S72), which is well characterized for the occurrences of these motifs. Finally, we search these motifs against the RNA structures in the entire Protein Data Bank and the abundances of them are estimated. RNAMotifScan is freely available at our supplementary website (http://genome.ucf.edu/RNAMotifScan).     

High base pair opening rates in tracts of GC base pairs

Sequence-dependent structural features of the DNA double helix have a strong influence on the base pair opening dynamics. Here we report a detailed study of the kinetics of base pair breathing in tracts of GC base pairs in DNA duplexes derived from 1H NMR measurements of the imino proton exchange rates upon titration with the exchange catalyst ammonia. In the limit of infinite exchange catalyst concentration, the exchange times of the guanine imino protons of the GC tracts extrapolate to much shorter base pair lifetimes than commonly observed for isolated GC base pairs. The base pair lifetimes in the GC tracts are below 5 ms for almost all of the base pairs. The unusually rapid base pair opening dynamics of GC tracts are in striking contrast to the behavior of AT tracts, where very long base pair lifetimes are observed. The implication of these findings for the structural principles governing spontaneous helix opening as well as the DNA-binding specificity of the cytosine-5-methyltransferases, where flipping of the cytosine base has been observed, are discussed.     

Proton exchange and base pair opening in a DNA triple helix

Nuclear magnetic resonance spectroscopy has been used to characterize opening reactions and stabilities of individual base pairs in two related DNA structures. The first is the triplex structure formed by the DNA 31-mer 5'-AGAGAGAACCCCTTCTCTCTTTTTCTCTCTT-3'. The structure belongs to the YRY (or parallel) family of triple helices. The second structure is the hairpin double helix formed by the DNA 20-mer 5'-AGAGAGAACCCCTTCTCTCT-3' and corresponds to the duplex part of the YRY triplex. The rates of exchange of imino protons with solvent in the two structures have been measured by magnetization transfer from water and by real-time exchange at 10 degrees C in 100 mM NaCl and 5 mM MgCl2 at pH 5.5 and in the presence of two exchange catalysts. The results indicate that the exchange of imino protons in protonated cytosines is most likely limited by the opening of Hoogsteen C+G base pairs. The base pair opening parameters estimated from imino proton exchange rates suggest that the stability of individual Hoogsteen base pairs in the DNA triplex is comparable to that of Watson-Crick base pairs in double-helical DNA. In the triplex structure, the exchange rates of imino protons in Watson-Crick base pairs are up to 5000-fold lower than those in double-helical DNA. This result suggests that formation of the triplex structure enhances the stability of Watson-Crick base pairs by up to 5 kcal/mol. This stabilization depends on the specific location of each triad in the triplex structure.     

Assignment and modeling of the Rev Response Element RNA bound to a Rev peptide using 13C-heteronuclear NMR

Summary The Rev Response Element (RRE) RNA-Rev protein interaction is important for regulation of gene expression in the human immunodeficiency virus. A model system for this interaction, which includes stem IIB of the RRE RNA and an arginine-rich peptide from the RNA-binding domain of Rev, was studied using multidimensional heteronuclear NMR. Assignment of the RNA when bound to the peptide was obtained from NMR experiments utilizing uniformly and specifically ¹³C-labeled RNA. Isotopic filtering experiments on the specifically labeled RNA enabled unambiguous assignment of unusual nonsequential NOE patterns present in the internal loop of the RRE. A three-dimensional model of the RNA in the complex was obtained using restrained molecular dynamics calculations. The internal loop contains two purine-purine base pairs, which are stacked to form one continuous helix flanked by two A-form regions. The formation of a G-G base pair in the internal loop requires an unusual structure of the phosphate backbone. This structural feature is consistent with mutational data as being important for the binding of Rev to the RRE. The G-G base pair may play an important role in opening the normally narrow major groove of A-form RNA to permit binding of the Rev basic domain.     

Isoalloxazine derivatives promote photocleavage of natural RNAs at G.U base pairs embedded within helices.

We have recently shown that isoalloxazine derivatives are able to photocleave RNA specifically at G.U base pairs embedded within a helical stack. The reaction involves the selective molecular recognition of G.U base pairs by the isoalloxazine ring and the removal of one nucleoside downstream of the uracil residue. Divalent metal ions are absolutely required for cleavage. Here we extend our studies to complex natural RNA molecules with known secondary and tertiary structures, such as tRNAs and a group I intron (td). G.U pairs were cleaved in accordance with the phylogenetically and experimentally derived secondary and tertiary structures. Tandem G.U pairs or certain G.U pairs located at a helix extremity were not affected. These new cleavage data, together with the RNA crystal structure, allowed us to perform molecular dynamics simulations to provide a structural basis for the observed specificity. We present a stable structural model for the ternary complex of the G. U-containing helical stack, the isoalloxazine molecule and a metal ion. This model provides significant new insight into several aspects of the cleavage phenomenon, mechanism and specificity for G. U pairs. Our study shows that in large natural RNAs a secondary structure motif made of an unusual base pair can be recognized and cleaved with high specificity by a low molecular weight molecule. This photocleavage reaction thus opens up the possibility of probing the accessibility of G.U base pairs, which are endowed with specific structural and functional roles in numerous structured and catalytic RNAs and interactions of RNA with proteins, in folded RNAs.     

The snRNP 15.5K protein folds its cognate K-turn RNA: a combined theoretical and biochemical study

The human 15.5K protein binds to the 5' stem-loop of U4 snRNA, promotes the assembly of the spliceosomal U4/U6 snRNP, and is required for the recruitment of the 61K protein and the 20/60/90K protein complex to the U4 snRNA. In the crystallographic structure of the 15.5K-U4 snRNA complex, the conformation of the RNA corresponds to the family of kink-turn (K-turn) structural motifs. We simulated the complex and the free RNA, showing how the protein binding and the intrinsic flexibility contribute to the RNA folding process. We found that the RNA is significantly more flexible in the absence of the 15.5K protein. Conformational transitions such as the interconversion between alternative purine stacking schemes, the loss of G-A base pairs, and the opening of the K-turn occur only in the free RNA. Furthermore, the stability of one canonical G-C base pair is influenced both by the binding of the 15.5K protein and the nature of the adjacent structural element in the RNA. We performed chemical RNA modification experiments and observed that the free RNA lacks secondary structure elements, a result in excellent agreement with the simulations. Based on these observations, we propose a protein-assisted RNA folding mechanism in which the RNA intrinsic flexibility functions as a catalyst.     

RNA-protein interaction. An analysis with RNA oligonucleotides of the recognition by alpha-sarcin of a ribosomal domain critical for function

alpha-Sarcin is a cytotoxic protein that inactivates ribosomes by hydrolyzing a single phosphodiester bond on the 3' side of G-4325 in eukaryotic 28 S rRNA. We have examined the requirements for the recognition by alpha-sarcin of this domain using a synthetic oligoribonucleotide (35-mer) that reproduces the sequence and, we presume, the secondary structure (a stem, a bulged nucleotide, and a loop) at the site of modification. The wild type structure and a large number of variants were transcribed in vitro from synthetic DNA templates with phage T7 RNA polymerase. Recognition of the substrate is strongly favored by a G at the position that corresponds to 4325. There is an absolute requirement for a helical stem; however, it can be reduced from the 7 base pairs in the natural structure to 3 without loss of specificity. The nature of the base pairs in the stem modifies but does not abolish recognition; whereas, the bulged nucleotide does not contribute to identification. Cleavage is materially affected by altering the nucleotides in the universal sequence surrounding G-4325 and changing the position in the loop of the tetranucleotide GAG(sarcin)A leads to loss of recognition by the toxin. We propose that the alpha-sarcin domain RNA participates in elongation factor catalyzed binding of aminoacyl-tRNA and of translocation; that translocation is driven by transitions in the structure of the alpha-sarcin domain RNA initiated by the binding of the factors, or the hydrolysis of GTP, or both; and that to toxin inactivates the ribosomes by preventing this transition.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

The G x U wobble base pair. A fundamental building block of RNA structure crucial to RNA function in diverse biological systems

The G x U wobble base pair is a fundamental unit of RNA secondary structure that is present in nearly every class of RNA from organisms of all three phylogenetic domains. It has comparable thermodynamic stability to Watson-Crick base pairs and is nearly isomorphic to them. Therefore, it often substitutes for G x C or A x U base pairs. The G x U wobble base pair also has unique chemical, structural, dynamic and ligand-binding properties, which can only be partially mimicked by Watson-Crick base pairs or other mispairs. These features mark sites containing G x U pairs for recognition by proteins and other RNAs and allow the wobble pair to play essential functional roles in a remarkably wide range of biological processes.     

Ribozymes correctly cleave a model substrate and endogenous RNA in vivo

The alpha-sarcin domain of 28 S RNA in Xenopus oocytes is attacked by several catalytic toxins (e.g. alpha-sarcin and ricin) that abolish protein synthesis. We synthesized 6 ribozymes targeted to the alpha-sarcin domain and to an oligoribonucleotide (34-mer) that mimics this domain. Sarcin ribozyme 5 (SR5) efficiently cleaved after the CUC site in the synthetic 34-mer in vitro at 50 degrees C. SR5 also cut the same site when both substrate and ribozyme were coinjected or injected separately into oocytes at 18 degrees C. Correct cleavage in vivo was shown by isolating and sequencing the large cleavage fragment. The cleavage reaction appeared to function equally well in the oocyte nucleus and cytoplasm. SR5 also correctly cleaved endogenous 28 S RNA in oocytes, although cutting was much less efficient than with alpha-sarcin. We therefore demonstrated that a ribozyme specifically cuts both a model substrate and a cellular RNA in vivo. Earlier work showed that certain injected deoxyoligonucleotides complementary to the alpha-sarcin region abolish protein synthesis. Oocyte protein synthesis was also abolished by an SR5 containing a single G----U substitution that inactivates RNA catalysis, indicating that SR5's translational suppression was perhaps due to antisense function rather than ribozyme cleavage.     

Motif prediction in ribosomal RNAs Lessons and prospects for automated motif prediction in homologous RNA molecules

The traditional way to infer RNA secondary structure involves an iterative process of alignment and evaluation of covariation statistics between all positions possibly involved in basepairing. Watson-Crick basepairs typically show covariations that score well when examples of two or more possible basepairs occur. This is not necessarily the case for non-Watson-Crick basepairing geometries. For example, for sheared (trans Hoogsteen/Sugar edge) pairs, one base is highly conserved (always A or mostly A with some C or U), while the other can vary (G or A and sometimes C and U as well). RNA motifs consist of ordered, stacked arrays of non-Watson-Crick basepairs that in the secondary structure representation form hairpin or internal loops, multi-stem junctions, and even pseudoknots. Although RNA motifs occur recurrently and contribute in a modular fashion to RNA architecture, it is usually not apparent which bases interact and whether it is by edge-to-edge H-bonding or solely by stacking interactions. Using a modular sequence-analysis approach, recurrent motifs related to the sarcin-ricin loop of 23S RNA and to loop E from 5S RNA were predicted in universally conserved regions of the large ribosomal RNAs (16S- and 23S-like) before the publication of high-resolution, atomic-level structures of representative examples of 16S and 23S rRNA molecules in their native contexts. This provides the opportunity to evaluate the predictive power of motif-level sequence analysis, with the goal of automating the process for predicting RNA motifs in genomic sequences. The process of inferring structure from sequence by constructing accurate alignments is a circular one. The crucial link that allows a productive iteration of motif modeling and realignment is the comparison of the sequence variations for each putative pair with the corresponding isostericity matrix to determine which basepairs are consistent both with the sequence and the geometrical data.     

Effect of single-residue bulges on RNA double-helical structures: crystallographic database analysis and molecular dynamics simulation studies

Asymmetric bulge loop motifs are widely dispersed in all types of functional RNAs. They are frequently occurring structural motifs in folded RNA structures and appear commonly in pre-microRNA and ribosomes, where they are involved in specific RNA–RNA and RNA–protein interactions. It is therefore necessary to understand such motifs from a structural point of view. We analyzed all available RNA structures and identified quite a few fragments of double helices that contain bulges. We found that these discontinuities often introduce kinks into the double helices, which also affects the stacking overlap between the base pairs across the irregularity. In order to understand the influence of these bulges on stability and flexibility, we carried out molecular dynamics simulations of three different single-residue bulge-containing RNA helices using the CHARMM36 force field. The structural variability at the junctions of RNA bulges is expected to differ from that in continuous double-helical stretches. The structural features of the junction region were observed to vary noticeably depending on the orientation of the bulge residue. When the base of the bulge residue is looped out, the RNA stretch behaves like a standard long A-form RNA double helix, whereas the entire RNA behaves differently when the base of the bulge residue is intercalated between base pairs inside the RNA stem. Such single-base intercalation was found to introduce a permanent kink into the composite double helix, which could be a recognition element for Dicer during the maturation of miRNA.     

Modular RNA architecture revealed by computational analysis of existing pseudoknots and ribosomal RNAs

Modular architecture is a hallmark of RNA structures, implying structural, and possibly functional, similarity among existing RNAs. To systematically delineate the existence of smaller topologies within larger structures, we develop and apply an efficient RNA secondary structure comparison algorithm using a newly developed two-dimensional RNA graphical representation. Our survey of similarity among 14 pseudoknots and subtopologies within ribosomal RNAs (rRNAs) uncovers eight pairs of structurally related pseudoknots with non-random sequence matches and reveals modular units in rRNAs. Significantly, three structurally related pseudoknot pairs have functional similarities not previously known: one pair involves the 3′ end of brome mosaic virus genomic RNA (PKB134) and the alternative hammerhead ribozyme pseudoknot (PKB173), both of which are replicase templates for viral RNA replication; the second pair involves structural elements for translation initiation and ribosome recruitment found in the viral internal ribosome entry site (PKB223) and the V4 domain of 18S rRNA (PKB205); the third pair involves 18S rRNA (PKB205) and viral tRNA-like pseudoknot (PKB134), which probably recruits ribosomes via structural mimicry and base complementarity. Additionally, we quantify the modularity of 16S and 23S rRNAs by showing that RNA motifs can be constructed from at least 210 building blocks. Interestingly, we find that the 5S rRNA and two tree modules within 16S and 23S rRNAs have similar topologies and tertiary shapes. These modules can be applied to design novel RNA motifs via build-up-like procedures for constructing sequences and folds.     

A three-dimensional RNA motif in Potato spindle tuber viroid mediates trafficking from palisade mesophyll to spongy mesophyll in Nicotiana benthamiana

Cell-to-cell trafficking of RNA is an emerging biological principle that integrates systemic gene regulation, viral infection, antiviral response, and cell-to-cell communication. A key mechanistic question is how an RNA is specifically selected for trafficking from one type of cell into another type. Here, we report the identification of an RNA motif in Potato spindle tuber viroid (PSTVd) required for trafficking from palisade mesophyll to spongy mesophyll in Nicotiana benthamiana leaves. This motif, called loop 6, has the sequence 5'-CGA-3'...5'-GAC-3' flanked on both sides by cis Watson-Crick G/C and G/U wobble base pairs. We present a three-dimensional (3D) structural model of loop 6 that specifies all non-Watson-Crick base pair interactions, derived by isostericity-based sequence comparisons with 3D RNA motifs from the RNA x-ray crystal structure database. The model is supported by available chemical modification patterns, natural sequence conservation/variations in PSTVd isolates and related species, and functional characterization of all possible mutants for each of the loop 6 base pairs. Our findings and approaches have broad implications for studying the 3D RNA structural motifs mediating trafficking of diverse RNA species across specific cellular boundaries and for studying the structure-function relationships of RNA motifs in other biological processes.     

Frequency and isostericity of RNA base pairs

Most of the hairpin, internal and junction loops that appear single-stranded in standard RNA secondary structures form recurrent 3D motifs, where non-Watson–Crick base pairs play a central role. Non-Watson–Crick base pairs also play crucial roles in tertiary contacts in structured RNA molecules. We previously classified RNA base pairs geometrically so as to group together those base pairs that are structurally similar (isosteric) and therefore able to substitute for each other by mutation without disrupting the 3D structure. Here, we introduce a quantitative measure of base pair isostericity, the IsoDiscrepancy Index (IDI), to more accurately determine which base pair substitutions can potentially occur in conserved motifs. We extract and classify base pairs from a reduced-redundancy set of RNA 3D structures from the Protein Data Bank (PDB) and calculate centroids (exemplars) for each base combination and geometric base pair type (family). We use the exemplars and IDI values to update our online Basepair Catalog and the Isostericity Matrices (IM) for each base pair family. From the database of base pairs observed in 3D structures we derive base pair occurrence frequencies for each of the 12 geometric base pair families. In order to improve the statistics from the 3D structures, we also derive base pair occurrence frequencies from rRNA sequence alignments.