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.     

In vitro selected GUAA tetraloop-binding receptors with structural plasticity and evolvability towards natural RNA structural modules

GNRA tetraloop-binding receptor interactions are key components in the macromolecular assembly of a variety of functional RNAs. In nature, there is an apparent bias for GAAA/11nt receptor and GYRA/helix interactions, with the former interaction being thermodynamically more stable than the latter. While past in vitro selections allowed isolation of novel GGAA and GUGA receptors, we report herein an in vitro selection that revealed several novel classes of specific GUAA receptors with binding affinities comparable to those from natural GAAA/11nt interactions. These GUAA receptors have structural homology with double-locked bulge RNA modules naturally occurring in ribosomal RNAs. They display mutational robustness that enables exploration of the sequence/phenotypic space associated to GNRA/receptor interactions through epistasis. Their thermodynamic self-assembly fitness landscape is characterized by a rugged neutral network with possible evolutionary trajectories toward natural GNRA/receptor interactions. High throughput sequencing analysis revealed synergetic mutations located away from the tertiary interactions that positively contribute to assembly fitness. Our study suggests that the repertoire of GNRA/receptor interactions is much larger than initially thought from the analysis of natural stable RNA molecules and also provides clues for their evolution towards natural GNRA/receptors.     

Comprehensive features of natural and in vitro selected GNRA tetraloop-binding receptors

Specific recognitions of GNRA tetraloops by small helical receptors are among the most widespread long-range packing interactions in large ribozymes. However, in contrast to GYRA and GAAA tetraloops, very few GNRA/receptor interactions have yet been identified to involve GGAA tetraloops in nature. A novel in vitro selection scheme based on a rigid self-assembling tectoRNA scaffold designed for isolation of intermolecular interactions with A-minor motifs has yielded new GGAA tetraloop-binding receptors with affinity in the nanomolar range. One of the selected receptors is a novel 12 nt RNA motif, (CCUGUG ... AUCUGG), that recognizes GGAA tetraloop hairpin with a remarkable specificity and affinity. Its physical and chemical characteristics are comparable to those of the well-studied '11nt' GAAA tetraloop receptor motif. A second less specific motif (CCCAGCCC ... GAUAGGG) binds GGRA tetraloops and appears to be related to group IC3 tetraloop receptors. Mutational, thermodynamic and comparative structural analysis suggests that natural and in vitro selected GNRA receptors can essentially be grouped in two major classes of GNRA binders. New insights about the evolution, recognition and structural modularity of GNRA and A-minor RNA-RNA interactions are proposed.     

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.     

Group I introns and gnra tetraloops: remnants of ‘The RNA world’?

GNRA tetraloops, found in high frequency in natural RNAs, make loop-receptor interactions, stabilizing the tertiary structure of Group I introns, a class of small RNAs. Analyzing 230 Group I introns, to study the distribution and sequence pattern of the GNRA tetraloops, we suggest that these features reflect the ancestral nature of these catalytic molecules, in a prebiotic RNA world. The adenosine rich GNRA tetraloops would have interacted with each other through long range RNA-RNA interactions to form higher order structures forming potential sites that render the propensity for the short RNAs to bind to metal ions from the prebiotic pool, aiding them to act as metalloenzymes.     

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.     

Intermolecular interaction between a branching ribozyme and associated homing endonuclease mRNA

RNA tertiary interactions involving docking of GNRA (N; any base; R; purine) hairpin loops into helical stem structures on other regions of the same RNA are one of the most common RNA tertiary interactions. In this study, we investigated a tertiary association between a GAAA hairpin tetraloop in a small branching ribozyme (DiGIR1) and a receptor motif (HEG P1 motif) present in a hairpin structure on a separate mRNA molecule. DiGIR1 generates a 2', 5' lariat cap at the 5' end of its downstream homing endonuclease mRNA by catalysing a self-cleavage branching reaction at an internal processing site. Upon release, the 5' end of the mRNA forms a distinct hairpin structure termed HEG P1. Our biochemical data, in concert with molecular 3D modelling, provide experimental support for an intermolecular tetraloop receptor interaction between the L9 GAAA in DiGIR1 and a GNRA tetraloop receptor-like motif (UCUAAG-CAAGA) found within the HEG P1. The biological role of this interaction appears to be linked to the homing endonuclease expression by promoting post-cleavage release of the lariat capped mRNA. These findings add to our understanding of how protein-coding genes embedded in nuclear ribosomal DNA are expressed in eukaryotes and controlled by ribozymes.     

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.     

In vitro selection of RNAs that undergo autolytic cleavage with Pb2+.

An in vitro selection method has been developed to obtain RNA molecules that specifically undergo autolytic cleavage reactions by Pb2+ ion. The method utilizes a circular RNA intermediate which is regenerated following the cleavage reaction to allow amplification and multiple cycles of selection. Pb2+ is known to catalyze a specific cleavage reaction between U17 and G18 of yeast tRNA(Phe). Starting from pools of RNA molecules which have a random distribution of sequences at nine or ten selected positions in the sequence of yeast tRNA(Phe), we have isolated many RNA molecules that undergo rapid and specific self-cleavage with Pb2+ at a variety of different sites. Terminal truncation experiments suggest that most of these self-cleaving RNA molecules do not fold like tRNA. However, two of the variants are cleaved rapidly with Pb2+ at U17 even though they lack the highly conserved nucleotides G18 and G19. Both specific mutations and terminal truncation experiments suggest that the D and T loops of these two variants interact in a manner similar to that of tRNA(Phe) despite the absence of the G18U55 and G19C56 tertiary interactions. A model for an alternate tertiary interaction involving a U17U55 pair is presented. This model may be relevant to the structure of about 100 mitochondrial tRNAs that also lack G18 and G19. The selection method presented here can be directly applied to isolate catalytic RNAs that undergo cleavage in the presence of other metal ions, modified nucleotides, or sequence-specific nucleases.     

In vitro selection and analysis of RNA aptamer recognize arginine-rich motif (ARM) model peptide on a QCM.

To study RNA-peptide interactions, we performed an in vitro selection of RNA on a 27 MHz quartz-crystal microbalance (QCM) on which a simple R5 helix peptide was immobilized as a model of N peptide from bacteriophade lambda. The consensus sequences including a GNRA tetraloop were obtained from a random RNA pool after the 7th cycle selection.     

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.     

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.     

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.     

The GAAA tetraloop-receptor interaction contributes differentially to folding thermodynamics and kinetics for the P4-P6 RNA domain

Tetraloops with the generic sequence GNRA are commonly found in RNA secondary structure, and interactions of such tetraloops with "receptors" elsewhere in RNA play important roles in RNA structure and folding. However, the contributions of tetraloop-receptor interactions specifically to the kinetics of RNA tertiary folding, rather than the thermodynamics of maintaining tertiary structure once folded, have not been reported. Here we investigate the role of the key GAAA tetraloop-receptor motif in folding of the P4-P6 domain of the Tetrahymena group I intron RNA. Insertions of one or more nucleotides into the tetraloop significantly disrupt the thermodynamics of tertiary folding; single-nucleotide insertions shift the folding free energy by 2-4 kcal/mol (DeltaDeltaG(o)'). The folding kinetics of several modified P4-P6 domains were determined by stopped-flow fluorescence spectroscopy, using an internally incorporated pyrene residue as the chromophore. In contrast to the thermodynamic results, the kinetics of Mg(2+)-induced folding were barely affected by the tetraloop modifications, with a DeltaDeltaG(++) of 0.2-0.4 kcal/mol and a Phi value (ratio of the kinetic and thermodynamic contributions) of <0.1. These data indicate an early transition state for folding of P4-P6 with respect to forming the tetraloop-receptor contact, consistent with previous results for modifications elsewhere in P4-P6. We conclude that the GAAA tetraloop-receptor motif contributes little to the stabilization of the transition state for Mg(2+)-induced P4-P6 folding. Rather, the tetraloop-receptor motif acts to clamp the RNA once folding has occurred. This is the first report to correlate the kinetic and thermodynamic contributions of an important RNA tertiary motif, the GNRA tetraloop-receptor. The results are related to possible models for the Mg(2+)-induced folding of the P4-P6 RNA, including a model invoking rapid nonspecific electrostatic collapse.     

GNRA/receptor interacting modules: versatile modular units for natural and artificial RNA architectures

Interactions between GNRA tetraloops and their receptors are found frequently as modular units in various types of naturally occurring structured RNAs. Due to their functional importance, GNRA/receptor interactions have been studied extensively with regard to their 3D structures and biochemical and biophysical properties. Artificial non-natural GNRA/receptor modules have also been generated not only to obtain a better understanding of this class of motifs in natural RNA structures but also for application of these modular units to the design and construction of artificial RNA structures that can be used as platforms to generate functional RNAs applicable for nanobiotechnology. In this review, we present a survey of structures, functions, and analyses as well as artificial generation and application of GNRA/receptor interacting modules.     

A comparative study on two GNRA-tetraloop receptors: 11-nt and IC3 motifs

Natural RNAs often contain terminal loops consisting of GNRA (N=A, G, C, U; R=A, G) and their receptors, which bind to the loops via long-range RNA-RNA interactions. Among several known receptors, two characteristic structural elements have been identified that are termed the 11-nt motif (CCUAAG-UAUGG) and IC3 motif (CCCUAAC-GAGGG). These two motifs that share a similar secondary structure have been shown to exhibit distinctively different binding specificities. The 11-nt motif recognizes a GAAA loop with highest specificity among the known receptors, whereas the IC3 motif distinguishes GAAA from other GNRA loops less stringently than any other receptors. To identify the elements in the receptors that determine the binding specificity, a series of chimeric receptors derived from the two motifs were prepared and their properties were examined. We identified characteristic base-pairs and a particular U residue in the receptors as such elements by means of a gel mobility shift assay that evaluates the degree of the tetraloop-receptor interaction. The relationship between the elements and the specificity is discussed together with a model that describes a possible evolutional linkage between the two receptors.     

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

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

Solution structure of an ATP-binding RNA aptamer reveals a novel fold.

In vitro selection has been used to isolate several RNA aptamers that bind specifically to biological cofactors. A well-characterized example in the ATP-binding RNA aptamer family, which contains a conserved 11-base loop opposite a bulged G and flanked by regions of double-stranded RNA. The nucleotides in the consensus sequence provide a binding pocket for ATP (or AMP), which binds with a Kd in the micromolar range. Here we present the three-dimensional solution structure of a 36-nucleotide ATP-binding RNA aptamer complexed with AMP, determined from NMR-derived distance and dihedral angle restraints. The conserved loop and bulged G form a novel compact, folded structure around the AMP. The backbone tracing of the loop nucleotides can be described by a Greek zeta (zeta). Consecutive loop nucleotides G, A, A form a U-turn at the bottom of the zeta, and interact with the AMP to form a structure similar to a GNRA tetraloop, with AMP standing in for the final A. Two asymmetric G. G base pairs close the stems flanking the internal loop. Mutated aptamers support the existence of the tertiary interactions within the consensus nucleotides and with the AMP found in the calculated structures.     

Equilibrium conformational dynamics in an RNA tetraloop from massively parallel molecular dynamics

Conformational equilibrium within the ubiquitous GNRA tetraloop motif was simulated at the ensemble level, including 10 000 independent all-atom molecular dynamics trajectories totaling over 110 µs of simulation time. This robust sampling reveals a highly dynamic structure comprised of 15 conformational microstates. We assemble a Markov model that includes transitions ranging from the nanosecond to microsecond timescales and is dominated by six key loop conformations that contribute to fluctuations around the native state. Mining of the Protein Data Bank provides an abundance of structures in which GNRA tetraloops participate in tertiary contact formation. Most predominantly observed in the experimental data are interactions of the native loop structure within the minor groove of adjacent helical regions. Additionally, a second trend is observed in which the tetraloop assumes non-native conformations while participating in multiple tertiary contacts, in some cases involving multiple possible loop conformations. This tetraloop flexibility can act to counterbalance the energetic penalty associated with assuming non-native loop structures in forming tertiary contacts. The GNRA motif has thus evolved not only to readily participate in simple tertiary interactions involving native loop structure, but also to easily adapt tetraloop secondary conformation in order to participate in larger, more complex tertiary interactions.