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.     

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.     

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.     

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

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

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.     

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.     

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 GANC tetraloop: a novel motif in the group IIC intron structure

Tetraloops are a common building block for RNA tertiary structure, and most tetraloops fall into one of three well-characterized classes: GNRA, UNCG, and CUYG. Here, we present the sequence and structure of a fourth highly conserved class of tetraloop that occurs only within the ζ-ζ′ interaction of group IIC introns. This GANC tetraloop was identified, along with an unusual cognate receptor, in the crystal structure of the group IIC intron and through phylogenetic analysis of intron RNA sequence alignments. Unlike conventional tetraloop-receptor interactions, which are stabilized by extensive hydrogen-bonding interactions, the GANC-receptor interaction is limited to a single base stack between the conserved adenosine of the tetraloop and a single purine of the receptor, which consists of a one- to three-nucleotide bulge and does not contain an A-platform. Unlike GNRA tetraloops, the GANC tetraloop forms a sharp angle relative to the adjacent helix, bending by approximately 45° toward the major groove side of the helix. These structural attributes allow GANC tetraloops to fit precisely within the group IIC intron core, thereby demonstrating that structural motifs can adapt to function in a specific niche.     

Rules for RNA recognition of GNRA tetraloops deduced by in vitro selection: comparison with in vivo evolution.

Terminal loops with a GNRA consensus sequence are a prominent feature of large self-assembling RNA molecules. In order to investigate tertiary interactions involving GNRA loops, we have devised an in vitro selection system derived from a group I ribozyme. Two selections, destined to isolate RNA sequences that would recognize two of the most widespread loops (GUGA and GAAA), yielded variants of previously identified receptors for those loops, and also some yet unrecognized, high-affinity binders with novel specificities towards members of the GNRA family. By taking advantage of available crystal structures, we have attempted to rationalize these results in terms of RNA-RNA contacts and to expose some of the structural principles that govern GNRA loop-mediated tertiary interactions; the role of loop nucleotide 2 in ensuring specific recognition by receptors is emphasized. More generally, comparison of the products of in vitro and natural selection is shown to provide insights into the mechanisms underlying the in vivo evolution of self-assembling RNA molecules.     

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.     

Structural and evolutionary classification of G/U wobble basepairs in the ribosome

We present a comprehensive structural, evolutionary and molecular dynamics (MD) study of the G/U wobble basepairs in the ribosome based on high-resolution crystal structures, including the recent Escherichia coli structure. These basepairs are classified according to their tertiary interactions, and sequence conservation at their positions is determined. G/U basepairs participating in tertiary interactions are more conserved than those lacking any interactions. Specific interactions occurring in the G/U shallow groove pocket--like packing interactions (P-interactions) and some phosphate backbone interactions (phosphate-in-pocket interactions)--lead to higher G/U conservation than others. Two salient cases of unique phylogenetic compensation are discovered. First, a P-interaction is conserved through a series of compensatory mutations involving all four participating nucleotides to preserve or restore the G/U in the optimal orientation. Second, a G/U basepair forming a P-interaction and another one forming a phosphate-in-pocket interaction are replaced by GNRA loops that maintain similar tertiary contacts. MD simulations were carried out on eight P-interactions. The specific GU/CG signature of this interaction observed in structure and sequence analysis was rationalized, and can now be used for improving sequence alignments.     

Membrane Insertion of Marginally Hydrophobic Transmembrane Helices Depends on Sequence Context

In mammalian cells, most integral membrane proteins are initially inserted into the endoplasmic reticulum membrane by the so-called Sec61 translocon. However, recent predictions suggest that many transmembrane helices (TMHs) in multispanning membrane proteins are not sufficiently hydrophobic to be recognized as such by the translocon. In this study, we have screened 16 marginally hydrophobic TMHs from membrane proteins of known three-dimensional structure. Indeed, most of these TMHs do not insert efficiently into the endoplasmic reticulum membrane by themselves. To test if loops or TMHs immediately upstream or downstream of a marginally hydrophobic helix might influence the insertion efficiency, insertion of marginally hydrophobic helices was also studied in the presence of their neighboring loops and helices. The results show that flanking loops and nearest-neighbor TMHs are sufficient to ensure the insertion of many marginally hydrophobic helices. However, for at least two of the marginally hydrophobic helices, the local interactions are not enough, indicating that post-insertional rearrangements are involved in the folding of these proteins.     

AA.AG@helix.ends: A:A and A:G base-pairs at the ends of 16 S and 23 S rRNA helices.

This study reveals that AA and AG oppositions occur frequently at the ends of helices in RNA crystal and NMR structures in the PDB database and in the 16 S and 23 S rRNA comparative structure models, with the G usually 3' to the helix for the AG oppositions. In addition, these oppositions are frequently base-paired and usually in the sheared conformation, although other conformations are present in NMR and crystal structures. These A:A and A:G base-pairs are present in a variety of structural environments, including GNRA tetraloops, E and E-like loops, interfaced between two helices that are coaxially stacked, tandem G:A base-pairs, U-turns, and adenosine platforms. Finally, given structural studies that reveal conformational rearrangements occurring in regions of the RNA with AA and AG oppositions at the ends of helices, we suggest that these conformationally unique helix extensions might be associated with functionally important structural rearrangements.     

Isolation and characterization of a family of stable RNA tetraloops with the motif YNMG that participate in tertiary interactions

RNA is known to fold into a variety of structural elements, many of which have sufficient sequence complexity to make the thermodynamic study of each possible variant impractical. We previously reported a method for isolating stable and unstable RNA sequences from combinatorial libraries using temperature gradient gel electrophoresis (TGGE). This method was used herein to analyze a six-nucleotide RNA hairpin loop library. Three rounds of in vitro selection were performed using TGGE, and unusually stable RNAs were identified by cloning and sequencing. Known stable tetraloops were found, including sequences belonging to the UNCG motif closed by a CG base pair, and the CUUG motif closed by a GC base pair. In addition, unknown tetraloops were found that were nearly as stable as cUNCGg, including sequences related through substitution of the U with a C (Y), the C with an A (M), or both. These substitutions allow hydrogen bonding and stacking interactions in the UNCG loop to be maintained. Thermodynamic analysis of YNMG and variant loops confirmed optimal stability with Y at position 1 and M at position 3. Similarity in structure and stability among YNMG loops was further supported by deoxyribose substitution, CD, and NMR experiments. A conserved tertiary interaction in 16S rRNA exists between a YAMG loop at position 343 and two adenines in the loop at position 159 (Escherichia coli numbering). NMR and functional group substitution experiments suggest that YNAG loops in particular have enhanced flexibility, which allows the tertiary interaction to be maintained with diverse loop sequences at position 159. Taken together, these results support the existence of an extended family of UNCG-like tetraloops with the motif cYNMGg that are thermodynamically stable and structurally similar and can engage in tertiary interactions in large RNA molecules.     

Promoting RNA helical stacking via A-minor junctions

RNA molecules take advantage of prevalent structural motifs to fold and assemble into well-defined 3D architectures. The A-minor junction is a class of RNA motifs that specifically controls coaxial stacking of helices in natural RNAs. A sensitive self-assembling supra-molecular system was used as an assay to compare several natural and previously unidentified A-minor junctions by native polyacrylamide gel electrophoresis and atomic force microscopy. This class of modular motifs follows a topological rule that can accommodate a variety of interchangeable A-minor interactions with distinct local structural motifs. Overall, two different types of A-minor junctions can be distinguished based on their functional self-assembling behavior: one group makes use of triloops or GNRA and GNRA-like loops assembling with helices, while the other takes advantage of more complex tertiary receptors specific for the loop to gain higher stability. This study demonstrates how different structural motifs of RNA can contribute to the formation of topologically equivalent helical stacks. It also exemplifies the need of classifying RNA motifs based on their tertiary structural features rather than secondary structural features. The A-minor junction rule can be used to facilitate tertiary structure prediction of RNAs and rational design of RNA parts for nanobiotechnology and synthetic biology.     

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.     

Tertiary Motifs Revealed in Analyses of Higher-Order RNA Junctions

RNA junctions are secondary-structure elements formed when three or more helices come together. They are present in diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze currently solved 3D RNA junctions in terms of base-pair interactions and 3D configurations. First, we study base-pair interaction diagrams for solved RNA junctions with 5 to 10 helices and discuss common features. Second, we compare these higher-order junctions to those containing 3 or 4 helices and identify global motif patterns such as coaxial stacking and parallel and perpendicular helical configurations. These analyses show that higher-order junctions organize their helical components in parallel and helical configurations similar to lower-order junctions. Their sub-junctions also resemble local helical configurations found in three- and four-way junctions and are stabilized by similar long-range interaction preferences such as A-minor interactions. Furthermore, loop regions within junctions are high in adenine but low in cytosine, and in agreement with previous studies, we suggest that coaxial stacking between helices likely forms when the common single-stranded loop is small in size; however, other factors such as stacking interactions involving noncanonical base pairs and proteins can greatly determine or disrupt coaxial stacking. Finally, we introduce the ribo-base interactions: when combined with the along-groove packing motif, these ribo-base interactions form novel motifs involved in perpendicular helix-helix interactions. Overall, these analyses suggest recurrent tertiary motifs that stabilize junction architecture, pack helices, and help form helical configurations that occur as sub-elements of larger junction networks. The frequent occurrence of similar helical motifs suggest nature's finite and perhaps limited repertoire of RNA helical conformation preferences. More generally, studies of RNA junctions and tertiary building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

Long membrane helices and short loops predicted less accurately

Low-resolution experiments suggest that most membrane helices span over 17-25 residues and that most loops between two helices are longer than 15 residues. Both constraints have been used explicitly in the development of prediction methods. Here, we compared the largest possible sequence-unique data sets from high- and low-resolution experiments. For the high-resolution data, we found that only half of the helices fall into the expected length interval and that half of the loops were shorter than 10 residues. We compared the accuracy of detecting short loops and long helices for 28 advanced and simple prediction methods: All methods predicted short loops less accurately than longer ones. In particular, loops shorter than 7 residues appeared to be very difficult to detect by current methods. Similarly, all methods tended to be more accurate for longer than for shorter helices. However, helices with more than 32 residues were predicted less accurately than all other helices. Our findings may suggest particular strategies for improving predictions of membrane helices.     

Achieving specificity in selected and wild-type N peptide-RNA complexes: the importance of discrimination against noncognate RNA targets

The boxB RNA pentaloops from the P22 and lambda phages each adopt a GNRA tetraloop fold upon binding their cognate arginine-rich N peptides. The third loop base in P22 boxB (3-out) and the fourth in lambda boxB (4-out) are excluded to accommodate this structure. Previously, we selected a pool of lambda N sequences with random amino acids at loop contacting positions 13-22 for binding to either of these two GNRA-folded pentaloops or a canonical GNRA tetraloop and isolated a class of peptides with a new conserved arginine (R15). Here, we characterize the binding of lambda N and these R15 peptides using fluorescent titrations with 2-aminopurine labeled versions of the three GNRA-folded loops and circular dichroism spectrometry. All peptides preferentially bind the lambda boxB RNA loop. lambda N and R15 peptide specificity against the P22 loop arises from the cost of rearranging its loop into the 4-out GNRA structure. Modeling indicates that the interaction of R8 with an additional loop phosphate in the 4-out GNRA pentaloop selectively stabilizes this complex relative to the tetraloop. R15 peptides gain additional discrimination against the tetraloop because their arginine also preferentially interacts with the 4-out GNRA pentaloop phosphate backbone, whereas K14 and W18 of lambda N contribute equal affinity when binding the tetraloop. Nonspecific electrostatic interactions by basic residues near the C-termini of these peptides create significantly steeper salt dependencies in association constants for noncognate loops, aiding discrimination at high salt concentrations. Our results emphasize the importance of considering specificity against noncognate as well as nonspecific targets in the combinatorial and rational design of biopolymers capable of macromolecular recognition.