Effects of osmolytes on stable UUCG tetraloops and their preference for a CG closing base pair

Osmolytes have the potential to affect the stability of secondary structure motifs and alter preferences for conserved nucleic acid sequences in the cell. To contribute to the understanding of the in vivo function of RNA we observed the effects of different classes of osmolytes on the UNCG tetraloop motif. UNCG tetraloops are the most common and stable of the RNA tetraloops and are nucleation sites for RNA folding. They also have a significant thermodynamic preference for a CG closing base pair. The thermal denaturation of model hairpins containing UUCG loops was monitored using UV-Vis spectroscopy in the presence of osmolytes with different chemical properties. Interestingly, all of the osmolytes tested destabilized the hairpins, but all had little effect on the thermodynamic preference for a CG base pair, except for polyethylene glycol (PEG) 200. PEG 200 destabilized the loop with the CG closing base pair relative to the loop with a GC closing base pair. The destabilization was linear with increasing concentrations of PEG 200, and the slope of this relationship was not perturbed by changes in the hairpin stem outside of the closing pair. This result suggests that in the presence of PEG 200, the UUCG loop with a GC closing base pair may retain some preferential interactions with the cosolute that are lost in the presence of the CG closing base pair. These results reveal that relatively small structural changes may influence how osmolytes tune the stability, and thus the function of a secondary structure motif in vivo.     

Three-dimensional motifs from the SCOR, structural classification of RNA database: extruded strands, base triples, tetraloops and U-turns

Release 2.0.1 of the Structural Classification of RNA (SCOR) database, http://scor.lbl.gov, contains a classification of the internal and hairpin loops in a comprehensive collection of 497 NMR and X-ray RNA structures. This report discusses findings of the classification that have not been reported previously. The SCOR database contains multiple examples of a newly described RNA motif, the extruded helical single strand. Internal loop base triples are classified in SCOR according to their three-dimensional context. These internal loop triples contain several examples of a frequently found motif, the minor groove AGC triple. SCOR also presents the predominant and alternate conformations of hairpin loops, as shown in the most well represented tetraloops, with consensus sequences GNRA, UNCG and ANYA. The ubiquity of the GNRA hairpin turn motif is illustrated by its presence in complex internal loops.     

The application of cluster analysis in the intercomparison of loop structures in RNA

We have developed a computational approach for the comparison and classification of RNA loop structures. Hairpin or interior loops identified in atomic resolution RNA structures were intercompared by conformational matching. The root-mean-square deviation (RMSD) values between all pairs of RNA fragments of interest, even if from different molecules, are calculated. Subsequently, cluster analysis is performed on the resulting matrix of RMSD distances using the unweighted pair group method with arithmetic mean (UPGMA). The cluster analysis objectively reveals groups of folds that resemble one another. To demonstrate the utility of the approach, a comprehensive analysis of all the terminal hairpin tetraloops that have been observed in 15 RNA structures that have been determined by X-ray crystallography was undertaken. The method found major clusters corresponding to the well-known GNRA and UNCG types. In addition, two tetraloops with the unusual primary sequence UMAC (M is A or C) were successfully assigned to the GNRA cluster. Larger loop structures were also examined and the clustering results confirmed the occurrence of variations of the GNRA and UNCG tetraloops in these loops and provided a systematic means for locating them. Nineteen examples of larger loops that closely resemble either the GNRA or UNCG tetraloop were found in the large ribosomal RNAs. When the clustering approach was extended to include all structures in the SCOR database, novel relationships were detected including one between the ANYA motif and a less common folding of the GAAA tetraloop sequence.     

Selection for thermodynamically stable DNA tetraloops using temperature gradient gel electrophoresis reveals four motifs: d(cGNNAg), d(cGNABg),d(cCNNGg), and d(gCNNGc)

Hairpins play important roles in the function of DNA, forming cruciforms and affecting processes such as replication and recombination. Temperature gradient gel electrophoresis (TGGE) and in vitro selection have been used to isolate thermodynamically stable DNA hairpins from a six-nucleotide random library. The TGGE-selection process was optimized such that known stable DNA tetraloops were recovered, and the selection appears to be exhaustive. In the selection, four families of exceptionally stable DNA loops were identified: d(cGNNAg), d(cGNABg), d(cCNNGg), and d(gCNNGc). (Lowercase denotes the closing base pair; N = A, C, G, or T; and B = C, G, or T.) It appears that the known stable d(cGNAg) triloop motif can be embedded into a tetraloop, with the extra nucleotide inserted into either the middle of the loop, d(cGNNAg), or at the 3'-end of the loop, d(cGNABg). For d(cGNNAg) and d(cGNABg), a CG closing base pair was strongly preferred over a GC, with DeltaDeltaG degrees (37) approximately 2 kcal/mol. Members of the two families, d(cCNNGg) and d(gCNNGc), are similar in stability. The loop sequences and closing base pairs identified for exceptionally stable DNA tetraloops show many similarities to those known for exceptionally stable RNA tetraloops. These data provide an expanded set of thermodynamic rules for the formation of tetraloops in DNA.     

Structural organization of a viral IRES depends on the integrity of the GNRA motif

Little is known about the tertiary structure of internal ribosome entry site (IRES) elements. The central domain of foot-and-mouth disease (FMDV) IRES, named 3 or I, contains a conserved GNRA motif, essential for IRES activity. We have combined functional analysis with RNA probing to define its structural organization. We have found that a UNCG motif does not functionally substitute the GNRA motif; moreover, binding of synthetic GNRA stem-loops to domain 3 was significantly reduced in RNAs bearing UCCG or GUAG substitutions. The apical region of domain 3 consists of a four-way junction where residues of the GNRA tetraloop are responsible for the organization of the adjacent stem-loops, as deduced from ribonucleases and dimethyl sulfate accessibility. A single A-to-G substitution in the fourth position of this motif led to a strong RNA reorganization, affecting several nucleotides away in the secondary structure of domain 3. The study of mutants bearing UNCG or GUAG tetraloops revealed lack of protection to chemical attack in native RNA at specific nucleotides relative to the parental GUAA, suggesting that the GNRA motif dictates the organization and stability of domain 3. This effect is likely mediated by the interaction with distant residues. Therefore, the GNRA motif plays a crucial role in the organization of IRES structure with important consequences on activity.     

A thermodynamic study of unusually stable RNA and DNA hairpins.

About 70% of the RNA tetra-loop sequences identified in ribosomal RNAs from different organisms fall into either (UNCG) or (GNRA) families (where N = A, C, G, or U; and R = A or G). RNA hairpins with these loop sequences form unusually stable tetra-loop structures. We have studied the RNA hairpin GGAC(UUCG)GUCC and several sequence variants to determine the effect of changing the loop sequence and the loop-closing base pair on the thermodynamic stability of (UNCG) tetra-loops. The hairpin GGAG(CUUG)CUCC with the conserved loop G(CUUG)C was also unusually stable. We have determined melting temperatures (Tm), and obtained thermodynamic parameters for DNA hairpins with sequences analogous to stable RNA hairpins with (UNCG), C(GNRA)G, C(GAUA)G, and G(CUUG)C loops. DNA hairpins with (TTCG), (dUdUCG), and related sequences in the loop, unlike their RNA counterparts, did not form unusually stable hairpins. However, DNA hairpins with the consensus loop sequence C(GNRA)G were very stable compared to hairpins with C(TTTT)G or C(AAAA)G loops. The C(GATA)G and G(CTTG)C loops were also extra stable. The relative stabilities of the unusually stable DNA hairpins are similar to those observed for their RNA analogs.     

A conserved motif in group IC3 introns is a new class of GNRA receptor

Terminal tetraloops consisting of GNRA sequences are often found in biologically active large RNAs. The loops appear to contribute towards the organization of higher order RNA structures by forming specific tertiary interactions with their receptors. Group IC3 introns which possess a GAAA loop in the L2 region often have a phylogenetically conserved motif in their P8 domains. In this report, we show that this conserved motif stands as a new class of receptor that distinguishes the sequences of GNRA loops less stringently than previously known receptors. The motif can functionally substitute an 11 nt motif receptor in the Tetrahymena ribozyme. Its structural and functional similarity to one class of synthetic receptors obtained from in vitro selection is observed.     

A GCUA tetranucleotide loop found in the poliovirus oriL by in vivo SELEX (un)expectedly forms a YNMG-like structure: Extending the YNMG family with GYYA

The cloverleaf structure in the 5'-untranslated region of enterovirus RNA that regulates viral RNA replication contains an evolutionarily conserved YNMG tetraloop closed by a Y-G base pair. This loop is believed to interact specifically with the viral protease 3C. To further characterize the specificity of this interaction, the tetraloop and two flanking base pairs of the poliovirus RNA were randomized, and viable viral clones were obtained using in vivo SELEX. Among many different mutants with the canonical YNMG sequences to be described elsewhere, a large-plaque-forming clone contained a deviating uGCUAg sequence. The NMR structure of a small hairpin capped with uGCUAg that we present here shows that the GCUA tetraloop adopts a novel fold, which is highly similar to that of the YNMG tetraloop with common stacking properties and hydrogen-bond interactions including an unusual syn conformation of the adenosine. Thermodynamic studies show moderate stabilities of hairpins with canonical YNMG and the novel GCUA loops, which, together with the similarity of spatial structures, illustrates that the tetraloop structure itself is crucial for the RNA-protein interaction required for the viral replication. A re-evaluation of the ribosomal secondary structure database reveals a hairpin containing a GCUA loop, which covaries with YNMG and is involved in a tertiary interaction, and in the 50S ribosomal subunit from Haloarcula marismortui the structurally comparable apex of stem-loop 35a is a recognition site for protein L2. These observations show a more general occurrence and importance of the so-far unrecognized GYYA hairpin loops.     

Thermodynamic parameters for loop formation in RNA and DNA hairpin tetraloops.

We determined the melting temperatures (Tm) and thermodynamic parameters of 15 RNA and 19 DNA hairpins at 1 M NaCl, 0.01 M sodium phosphate, 0.1 mM EDTA, at pH 7. All these hairpins have loops of four bases, the most common loop size in 16S and 23S ribosomal RNAs. The RNA hairpins varied in loop sequence, loop-closing base pair (A.U, C.G, or G.C), base sequence of the stem, and stem size (four or five base pairs). The DNA hairpins varied in loop sequence, loop-closing base pair (C.G, or G.C), and base sequence of the four base-pair stem. Thermodynamic properties of a hairpin may be represented by nearest-neighbor interactions of the stem plus contributions from the loop. Thus, we obtained thermodynamic parameters for the formation of RNA and DNA tetraloops. For the tetraloops we studied, a free energy of loop formation (at 37 degrees C) of about +3 kcal/mol is most common for either RNA or DNA. There are extra stable loops with delta G degrees 37 near +1 kcal/mol, but the sequences are not necessarily the same for RNA and DNA. The closing base pair is also important; changing from C.G to G.C lowered the stability of several tetraloops in both RNA and DNA. These values will be useful in predicting RNA and DNA secondary structures.     

Isolation and characterization of thermodynamically stable and unstable RNA hairpins from a triloop combinatorial library

Hairpins are the most common elements of RNA secondary structure, playing important roles in RNA tertiary architecture and forming protein binding sites.Triloops are common in a variety of naturally occurring RNA hairpins, but little is known about their thermodynamic stability. Reported here are the sequences and thermodynamic parameters for a variety of stable and unstable triloop hairpins. Temperature gradient gel electrophoresis (TGGE) can be used to separate a simple RNA combinatorial library based on thermal stability [Bevilacqua, J. M., and Bevilacqua, P. C. (1998) Biochemistry 45, 15877-15884]. Here we introduce the application of TGGE to separating and analyzing a complex RNA combinatorial library based on thermal stability, using an RNA triloop library. Several rounds of in vitro selection of an RNA triloop library were carried out using TGGE, and preferences for exceptionally stable and unstable closing base pairs and loop sequences were identified. For stable hairpins, the most common closing base pair is CG, and U-rich loop sequences are preferred. Closing base pairs of GC and UA result in moderately stable hairpins when combined with a stable loop sequence. For unstable hairpins, the most common closing base pairs are AU and UG, and U-rich loop sequences are no longer preferred. In general, the contributions of the closing base pair and loop sequence to overall hairpin stability appear to be additive. Thermodynamic parameters for individual hairpins determined by UV melting are generally consistent with outcomes from selection experiments, with hairpins containing a CG closing base pair having a DeltaDeltaG degrees (37) 2.1-2.5 kcal/mol more favorable than hairpins with other closing base pairs. Sequences and thermodynamic rules for triloop hairpins should aid in RNA structure prediction and determination of whether naturally occurring triloop hairpins are thermodynamically stable.     

A test of the model to predict unusually stable RNA hairpin loop stability

To investigate the accuracy of a model [Giese et al., 1998, Biochemistry37:1094-1100 and Mathews et al., 1999, JMol Biol 288:911-940] that predicts the stability of RNA hairpin loops, optical melting studies were conducted on sets of hairpins previously determined to have unusually stable thermodynamic parameters. Included were the tetraloops GNRA and UNCG (where N is any nucleotide and R is a purine), hexaloops with UU first mismatches, and the hairpin loop of the iron responsive element, CAGUGC. The experimental values for the GNRA loops are in excellent agreement (deltaG degrees 37 within 0.2 kcal/mol and melting temperature (TM) within 4 degrees C) with the values predicted by the model. When the UNCG hairpin loops are treated as tetraloops, and a bonus of 0.8 kcal/mol included in the prediction to account for the extra stable first mismatch (UG), the measured and predicted values are also in good agreement (deltaG degrees 37 within 0.7 kcal/mol and TM within 3 degrees C). Six hairpins with unusually stable UU first mismatches also gave good agreement with the predictions (deltaG degrees 37 within 0.5 kcal/mol and TM within 8 degrees C), except for hairpins closed by wobble base pairs. For these hairpins, exclusion of the additional stabilization term for UU first mismatches improved the prediction (AG degrees 37 within 0.1 kcal/mol and TM within 3 degrees C). Hairpins with the iron-responsive element loop were not predicted well by the model, as measured deltaG degrees 37 values were at least 1 kcal/mol greater than predicted.     

Common and distinctive features of GNRA tetraloops based on a GUAA tetraloop structure at 1.4 A resolution

GNRA tetraloops (N is A, C, G, or U; R is A or G) are basic building blocks of RNA structure that often interact with proteins or other RNA structural elements. Understanding sequence-dependent structural variation among different GNRA tetraloops is an important step toward elucidating the molecular basis of specific GNRA tetraloop recognition by proteins and RNAs. Details of the geometry and hydration of this motif have been based on high-resolution crystallographic structures of the GRRA subset of tetraloops; less is known about the GYRA subset (Y is C or U). We report here the structure of a GUAA tetraloop determined to 1.4 A resolution to better define these details and any distinctive features of GYRA tetraloops. The tetraloop is part of a 27-nt structure that mimics the universal sarcin/ricin loop from Escherichia coli 23S ribosomal RNA in which a GUAA tetraloop replaces the conserved GAGA tetraloop. The adenosines of the GUAA tetraloop form an intermolecular contact that is a commonplace RNA tertiary interaction called an A-minor motif. This is the first structure to reveal in great detail the geometry and hydration of a GUAA tetraloop and an A-minor motif. Comparison of tetraloop structures shows a common backbone geometry for each of the eight possible tetraloop sequences and suggests a common hydration. After backbone atom superposition, equivalent bases from different tetraloops unexpectedly depart from coplanarity by as much as 48 degrees. This variation displaces the functional groups of tetraloops implicated in protein and RNA binding, providing a recognition feature.     

稳定的RNA发夹结构及其生物学功能

以UNCG、GNRA、CUUG(N=A、U、C或G,R=G或A)为端环能够形成稳定的、保守的发夹结构。高分辨率的溶液结构、晶体结构和计算机模拟等方法从原子水平上解析了这些发夹特殊的结构特征。在体内,它们发挥着重要的生物学功能:在折叠过程中作为折叠的起始位置帮助组织RNA分子正确折叠;与核酸受体结合参与三级相互作用;与蛋白质发生相互作用;阻止逆转录酶的延伸等等。另外,由于C(UUCG)G发夹极其稳定的特征,在体外RNA分子的实验测定中它还是稳定核酸结构的理想工具。这些稳定的发夹广泛分布于体内rRNA、催化RNA和非编码mRNA中。但在对人类编码区mRNA结构特征的研究当中,却未发现C(UUCG)G发夹。     

NMR structure of the 3' stem-loop from human U4 snRNA

The NMR structure of the 3' stem-loop (3'SL) from human U4 snRNA was determined to gain insight into the structural basis for conservation of this stem-loop sequence from vertebrates. 3'SL sequences from human, rat, mouse and chicken U4 snRNA each consist of a 7 bp stem capped by a UACG tetraloop. No high resolution structure has previously been reported for a UACG tetraloop. The UACG tetraloop portion of the 3'SL was especially well defined by the NMR data, with a total of 92 NOE-derived restraints (about 15 per residue), including 48 inter-residue restraints (about 8 per residue) for the tetraloop and closing C-G base pair. Distance restraints were derived from NOESY spectra using MARDIGRAS with random error analysis. Refinement of the 20mer RNA hairpin structure was carried out using the programs DYANA and miniCarlo. In the UACG tetraloop, U and G formed a base pair stabilized by two hydrogen bonds, one between the 2'-hydroxyl proton of U and carbonyl oxygen of G, another between the imino proton of G and carbonyl oxygen O2 of U. In addition, the amino group of C formed a hydrogen bond with the phosphate oxygen of A. G adopted a syn orientation about the glycosidic bond, while the sugar puckers of A and C were either C2'-endo or flexible. The conformation of the UACG tetraloop was, overall, similar to that previously reported for UUCG tetraloops, another member of the UNCG class of tetraloops. The presence of an A, rather than a U, at the variable position, however, presents a distinct surface for interaction of the 3'SL tetraloop with either RNA or protein residues that may stabilize interactions important for active spliceosome formation. Such tertiary interactions may explain the conservation of the UACG tetraloop motif in 3'SL sequences from U4 snRNA in vertebrates.     

Predicting U-turns in ribosomal RNA with comparative sequence analysis

The U-turn is a well-known RNA motif characterized by a sharp reversal of the RNA backbone following a single-stranded uridine base. In experimentally determined U-turn motifs, the nucleotides 3' to the turn are frequently involved in tertiary interactions, rendering this motif particularly attractive in RNA modeling and functional studies. The U-turn signature is composed of an UNR sequence pattern flanked by a Y:Y, Y:A (Y=pyrimidine) or G:A base juxtaposition. We have identified 33 potential UNR-type U-turns and 25 related GNRA-type U-turns in a large set of aligned 16 S and 23 S rRNA sequences. U-turn candidates occur in hairpin loops (34 times) as well as in internal and multi-stem loops (24 times). These are classified into ten families based on loop type, sequence pattern (UNR or GNRA) and the nature of the closing base juxtaposition. In 13 cases, the bases on the 3' side of the turn, or on the immediate 5' side, are involved in tertiary covariations, making these sites strong candidates for tertiary interactions.     

UNAC tetraloops: to what extent do they mimic GNRA tetraloops?

The structures of four small RNAs each containing a different version of the UNAC loop were determined in solution using NMR spectroscopy and restrained molecular dynamics. The UMAC tetraloops (where M is A or C) exhibited a typical GNRA fold including at least one hydrogen bond between the first U and fourth C. In contrast, UGAC and UUAC tetraloops have a different orientation of the first and fourth residues, such that they do not closely mimic the GNRA fold. Although the UMAC tetraloops are excellent structural mimics of the GNRA tetraloop backbone, sequence comparisons typically do not reveal co‐variation between the two loop types. The limited covariation is attributed to differences in the location of potential hydrogen bond donors and acceptors as a result of the replacement of the terminal A of GNRA with C in the UMAC version. Thus, UMAC loops do not readily form the common GNRA tetraloop‐receptor interaction. The loop at positions 863‐866 in E. coli 16S ribosomal RNA appears to be a major exception. However, in this case the GNRA loop does not in fact engage in the usual base to backbone tertiary interactions. In summary, UMAC loops are not just an alternative sequence version of the GNRA loop family, but instead they expand the types of interactions, or lack thereof, that are possible. From a synthetic biology perspective their inclusion in an artificial RNA may allow the establishment of a stable loop structure while minimizing unwanted long range interactions or permitting alternative long‐range interactions. © 2012 Wiley Periodicals, Inc. Biopolymers 97: 617–628, 2012.     

Co-conservation of rRNA tetraloop sequences and helix length suggests involvement of the tetraloops in higher-order interactions

Terminal loops containing four nucleotides (tetraloops) are common in structural RNAs, and they frequently conform to one of three sequence motifs, GNRA, UNCG, or CUUG. Here we compare available sequences and secondary structures for rRNAs from bacteria, and we show that helices capped by phylogenetically conserved GNRA loops display a strong tendency to be of conserved length. The simplest interpretation of this correlation is that the conserved GNRA loops are involved in higher-order interactions, intramolecular or intermolecular, resulting in a selective pressure for maintaining the lengths of these helices. A small number of conserved UNCG loops were also found to be associated with conserved length helices, consistent with the possibility that this type of tetraloop also takes part in higher-order interactions.     

Characterization of an RNA receptor motif that recognizes a GCGA tetraloop

Tertiary interactions between a new RNA motif and RNA tetraloops were analyzed to determine whether this new motif shows preference for a GCGA tetraloop. In the structural context of a ligase ribozyme, this motif discriminated GCGA loop from 3 other tetraloops. The affinity between the GCGA loop and its receptor is strong enough to carry out the ribozyme activity.     

Chrysanthemum chlorotic mottle viroid RNA: dissection of the pathogenicity determinant and comparative fitness of symptomatic and non-symptomatic variants

Chrysanthemum chlorotic mottle viroid (CChMVd) is a small RNA (398-401nt) with hammerhead ribozymes in both polarity strands that mediate self-cleavage of the oligomeric RNA intermediates generated in a rolling-circle mechanism of replication. Within the in vivo branched RNA conformation of CChMVd, a tetraloop has been identified as a major determinant of pathogenicity. Here we present a detailed study of this tetraloop by site-directed mutagenesis, bioassay of the CChMV-cDNA clones and analysis of the resulting progenies. None of the changes introduced in the tetraloop, including its substitution by a triloop or a pentaloop, abolished infectivity. In contrast to observations for other RNAs, the thermodynamically stable GAAA tetraloop characteristic of non-symptomatic CChMVd-NS strains was not functionally interchangeable for other stable tetraloops of the UNCG family, suggesting that the sequence, rather than the structure, is the major factor governing conservation of this motif. In most cases, the changes introduced initially led to symptomless infections, which eventually evolved to be symptomatic concurrently with the prevalence in the progeny of the UUUC tetraloop characteristic of symptomatic CChMVd-S strains. Only in one case did the GAAA tetraloop emerge and eventually dominate the progeny in infected plants that were non-symptomatic. These results revealed two major fitness peaks in the tetraloop (UUUC and GAAA), whose adjacent stem was also under strong selection pressure. Co-inoculations with CChMVd-S and -NS variants showed that only when the latter was in a 100- or 1000-fold excess did the infected plants remain symptomless, confirming the higher biological fitness of the S variant and explaining the lack of symptom expression previously observed in cross-protection experiments.     

Evidence that folding of an RNA tetraloop hairpin is less cooperative than its DNA counterpart

Hairpin secondary structural elements play important roles in the folding and function of RNA and DNA molecules. Previous work from our lab on small DNA hairpin loop motifs, d(cGNAg) and d(cGNABg) (where B is C, G, or T), showed that folding is highly cooperative and obeys indirect coupling, consistent with a concerted transition. Herein, we investigate folding of the related, exceptionally stable RNA hairpin motif, r(cGNRAg) (where R is A or G). Previous NMR characterization identified a complex network of seven hydrogen bonds in this loop. We inserted three carbon (C3) spacers throughout the loop and found coupling between G1 of the loop and the CG closing base pair, similar to that found in DNA. These data support a GNRA motif being expandable at any position but before the G. Thermodynamic measurements of nucleotide-analogue-substituted oligonucleotides revealed pairwise-coupling free energies ranging from weak to strong. When coupling free energies were remeasured in the background of changes at a third site, they remained essentially unchanged even though all of the sites were coupled to each other. This type of coupling, referred to as "direct", is peculiar to the RNA loop. The data suggest that, for small stable loops, folding of RNA obeys a model with nearest-neighbor interactions, while folding of DNA follows a more concerted process in which the stabilizing interactions are linked through a conformational change. The lesser cooperativity in RNA loops may provide a more robust loop that can withstand mutations without a severe loss in stability. These differences may enhance the ability of RNA to evolve.