Molecular Dynamics Study of Intrinsic Stability in Six RNA Terminal Loop Motifs

Abstract

Single stranded RNA molecules can assume a wide range of tertiary structures beyond the canonical A-form double helix. Certain sequences, termed motifs, are more common than a random distribution would suggest. The existence of such motifs can be rationalized in structural terms. In this study, we have investigated the intrinsic structural stability of RNA terminal loop motifs using multiple MD simulations in explicit water. Representative loops were chosen from the major tetraloop motifs, including also the U-turn motif. Not all loops retain their folded starting structure, but lowering the temperature to 277 K, or adding adjacent base pairs from the stem to which the motif is attached, helps stabilizing the folded loop structure.     

Automated classification of RNA 3D motifs and the RNA 3D Motif Atlas

The analysis of atomic-resolution RNA three-dimensional (3D) structures reveals that many internal and hairpin loops are modular, recurrent, and structured by conserved non-Watson–Crick base pairs. Structurally similar loops define RNA 3D motifs that are conserved in homologous RNA molecules, but can also occur at nonhomologous sites in diverse RNAs, and which often vary in sequence. To further our understanding of RNA motif structure and sequence variability and to provide a useful resource for structure modeling and prediction, we present a new method for automated classification of internal and hairpin loop RNA 3D motifs and a new online database called the RNA 3D Motif Atlas. To classify the motif instances, a representative set of internal and hairpin loops is automatically extracted from a nonredundant list of RNA-containing PDB files. Their structures are compared geometrically, all-against-all, using the FR3D program suite. The loops are clustered into motif groups, taking into account geometric similarity and structural annotations and making allowance for a variable number of bulged bases. The automated procedure that we have implemented identifies all hairpin and internal loop motifs previously described in the literature. All motif instances and motif groups are assigned unique and stable identifiers and are made available in the RNA 3D Motif Atlas (http://rna.bgsu.edu/motifs), which is automatically updated every four weeks. The RNA 3D Motif Atlas provides an interactive user interface for exploring motif diversity and tools for programmatic data access.     

Role of the closing base pair for d(GCA) hairpin stability: free energy analysis and folding simulations

Hairpin loops belong to the most important structural motifs in folded nucleic acids. The d(GNA) sequence in DNA can form very stable trinucleotide hairpin loops depending, however, strongly on the closing base pair. Replica-exchange molecular dynamics (REMD) were employed to study hairpin folding of two DNA sequences, d(gcGCAgc) and d(cgGCAcg), with the same central loop motif but different closing base pairs starting from single-stranded structures. In both cases, conformations of the most populated conformational cluster at the lowest temperature showed close agreement with available experimental structures. For the loop sequence with the less stable G:C closing base pair, an alternative loop topology accumulated as second most populated conformational state indicating a possible loop structural heterogeneity. Comparative-free energy simulations on induced loop unfolding indicated higher stability of the loop with a C:G closing base pair by ~3 kcal mol(-1) (compared to a G:C closing base pair) in very good agreement with experiment. The comparative energetic analysis of sampled unfolded, intermediate and folded conformational states identified electrostatic and packing interactions as the main contributions to the closing base pair dependence of the d(GCA) loop stability.     

The human T-cell lymphotropic virus type-I dimerization initiation site forms a hairpin loop, unlike previously characterized retroviral dimerization motifs

The formation of genomic RNA dimers during the retroviral life cycle is essential for optimal viral replication and infectivity. The sequences and RNA structures responsible for this interaction are located in the untranslated 5' leader RNA, along with other cis-acting signals. Dimer formation occurs by specific interaction between identical structural motifs. It is believed that an initial kissing hairpin forms following self-recognition by autocomplementary RNA loops, leading to formation of an extended stable duplex. The dimerization initiation site (DIS) of the deltaretrovirus human T-cell lymphotropic virus type-I (HTLV-I) has been previously localized to a 14-nucleotide sequence predicted to contain an RNA stem loop. Biochemical probing of the monomeric RNA structure using RNAse T1, RNAse V1, RNAse U2, lead acetate, and dimethyl sulfate has led to the generation of the first structural map of the HTLV-I DIS. A comprehensive data set of individual nucleotide modifications reveals that the structural motif responsible for HTLV-I RNA dimerization forms a trinucleotide RNA loop, unlike any previously characterized retroviral dimerization motif. Molecular modeling demonstrates that this can be formed by an unusual C:synG base pair closing the loop. Comparative phylogeny indicates that such a motif may also exist in other deltaretroviruses.     

Clustering RNA structural motifs in ribosomal RNAs using secondary structural alignment

RNA structural motifs are the building blocks of the complex RNA architecture. Identification of non-coding RNA structural motifs is a critical step towards understanding of their structures and functionalities. In this article, we present a clustering approach for de novo RNA structural motif identification. We applied our approach on a data set containing 5S, 16S and 23S rRNAs and rediscovered many known motifs including GNRA tetraloop, kink-turn, C-loop, sarcin-ricin, reverse kink-turn, hook-turn, E-loop and tandem-sheared motifs, with higher accuracy than the state-of-the-art clustering method. We also identified a number of potential novel instances of GNRA tetraloop, kink-turn, sarcin-ricin and tandem-sheared motifs. More importantly, several novel structural motif families have been revealed by our clustering analysis. We identified a highly asymmetric bulge loop motif that resembles the rope sling. We also found an internal loop motif that can significantly increase the twist of the helix. Finally, we discovered a subfamily of hexaloop motif, which has significantly different geometry comparing to the currently known hexaloop motif. Our discoveries presented in this article have largely increased current knowledge of RNA structural motifs.     

Secondary structure model of the Mason-Pfizer monkey virus 5' leader sequence: identification of a structural motif common to a variety of retroviruses.

A stable secondary structure model is presented for the region 3' of the primer-binding site to 130 bases into the gag sequence of the prototype type D retrovirus Mason-Pfizer monkey virus. Using biochemical probing of RNA from this region in association with free energy minimization, we have identified a stem-loop structure in the region, which from other studies has been shown to be important for genomic RNA encapsidation. The structure involves a highly stable stem of five G-C pairs terminating in a heptaloop. Comparison of the Mason-Pfizer monkey virus structure with one predicted for squirrel monkey retrovirus demonstrates an identical stem and a common ACC motif in the loop. Free energy studies of the secondary structure of the 5' regions of eight other retroviruses predict stem loops which have similar GAYC motifs. We believe this may represent a common structural and sequence motif which among other functions may be involved in genomic RNA packaging in these viruses.     

Ni2+-binding RNA motifs with an asymmetric purine-rich internal loop and a G-A base pair.

RNA molecules with high affinity for immobilized Ni2+ were isolated from an RNA pool with 50 randomized positions by in vitro selection-amplification. The selected RNAs preferentially bind Ni2+ and Co2+ over other cations from first series transition metals. Conserved structure motifs, comprising about 15 nt, were identified that are likely to represent the Ni2+ binding sites. Two conserved motifs contain an asymmetric purine-rich internal loop and probably a mismatch G-A base pair. The structure of one of these motifs was studied with proton NMR spectroscopy and formation of the G-A pair at the junction of helix and internal loop was demonstrated. Using Ni2+ as a paramagnetic probe, a divalent metal ion binding site near this G-A base pair was identified. Ni2+ ions bound to this motif exert a specific stabilization effect. We propose that small asymmetric purine-rich loops that contain a G-A interaction may represent a divalent metal ion binding site in RNA.     

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.     

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.     

A novel loop-loop recognition motif in the yeast ribosomal protein L30 autoregulatory RNA complex

The yeast Saccharomyces cerevisiae ribosomal protein L30 negatively autoregulates its production by binding to a helix-loop-helix structure formed in its pre-mRNA and its mRNA. A three-dimensional solution structure of the L30 protein in complex with its regulatory RNA has been solved using NMR spectroscopy. In the complex, the helix-loop-helix RNA adopts a sharply bent conformation at the internal loop region. Unusual RNA features include a purine stack, a reverse Hoogsteen base pair (G11anti-G56syn) and highly distorted backbones. The L30 protein is folded in a three-layer alpha/beta/alpha sandwich topology, and three loops at one end of the sandwich make base-specific contacts with the RNA internal loop. The protein-RNA binding interface is divided into two clusters, including hydrophobic and aromatic stacking interactions centering around G56, and base-specific hydrogen-bonding contacts to A57, G58 and G10-U60 wobble base pair. Both the protein and the RNA exhibit a partially induced fit for binding, where loops in the protein and the internal loop in the RNA become more ordered upon complex formation. The specific interactions formed between loops on L30 and the internal loop on the mRNA constitute a novel loop-loop recognition motif where an intimate RNA-protein interface is formed between regions on both molecules that lack regular secondary structure.     

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.     

NMR structure of a 4 x 4 nucleotide RNA internal loop from an R2 retrotransposon: identification of a three purine-purine sheared pair motif and comparison to MC-SYM predictions

The NMR solution structure is reported of a duplex, 5'GUGAAGCCCGU/3'UCACAGGAGGC, containing a 4 × 4 nucleotide internal loop from an R2 retrotransposon RNA. The loop contains three sheared purine-purine pairs and reveals a structural element found in other RNAs, which we refer to as the 3RRs motif. Optical melting measurements of the thermodynamics of the duplex indicate that the internal loop is 1.6 kcal/mol more stable at 37°C than predicted. The results identify the 3RRs motif as a common structural element that can facilitate prediction of 3D structure. Known examples include internal loops having the pairings: 5'GAA/3'AGG, 5'GAG/3'AGG, 5'GAA/3'AAG, and 5'AAG/3'AGG. The structural information is compared with predictions made with the MC-Sym program.     

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.     

Non-Watson Crick base pairs might stabilize RNA structural motifs in ribozymes — a comparative study of group-I intron structures

In recent decades studies on RNA structure and function have gained significance due to discoveries on diversified functions of RNA. A common element for RNA secondary structure formed by series of non-Watson/Watson Crick base pairs, internal loops and pseudoknots have been the highlighting feature of recent structural determination of RNAs. The recent crystal structure of group-I introns has demonstrated that these might constitute RNA structural motifs in ribozymes, playing a crucial role in their enzymatic activity. To understand the functional significance of these non-canonical base pairs in catalytic RNA, we analysed the sequences of group-I introns from nuclear genes. The results suggest that they might form the building blocks of folded RNA motifs which are crucial to the catalytic activity of the ribozyme. The conservation of these, as observed from divergent organisms, argues for the presence of non-canonical base pairs as an important requisite for the structure and enzymatic property of ribozymes by enabling them to carry out functions such as replication, polymerase activity etc. in primordial conditions in the absence of proteins.     

RNA aptamers that specifically bind to a 16S ribosomal RNA decoding region construct

RNA–RNA recognition is a critical process in controlling many key biological events, such as translation and ribozyme functions. The recognition process governing RNA–RNA interactions can involve complementary Watson–Crick (WC) base pair binding, or can involve binding through tertiary structural interaction. Hence, it is of interest to determine which of the RNA–RNA binding events might emerge through an in vitro selection process. The A-site of the 16S rRNA decoding region was chosen as the target, both because it possesses several different RNA structural motifs, and because it is the rRNA site where codon/anticodon recognition occurs requiring recognition of both mRNA and tRNA. It is shown here that a single family of RNA molecules can be readily selected from two different sizes of RNA library. The tightest binding aptamer to the A-site 16S rRNA construct, 109.2-3, has its consensus sequences confined to a stem–loop region, which contains three nucleotides complementary to three of the four nucleotides in the stem–loop region of the A-site 16S rRNA. Point mutations on each of the three nucleotides on the stem–loop of the aptamer abolish its binding capacity. These studies suggest that the RNA aptamer 109.2-3 interacts with the simple 27 nt A-site decoding region of 16S rRNA through their respective stem–loops. The most probable mode of interaction is through complementary WC base pairing, commonly referred to as a loop–loop ‘kissing’ motif. High affinity binding to the other structural motifs in the decoding region were not observed.     

Nuclear magnetic resonance structure of the Varkud satellite ribozyme stem-loop V RNA and magnesium-ion binding from chemical-shift mapping

An important step in the substrate recognition of the Neurospora Varkud Satellite (VS) ribozyme is the formation of a magnesium-dependent loop/loop interaction between the terminal loops of stem-loops I and V. We have studied the structure of stem-loop V by nuclear magnetic resonance spectroscopy and shown that it adopts a U-turn conformation, a common motif found in RNA. Structural comparisons indicate that the U-turn of stem-loop V fulfills some but not all of the structural characteristics found in canonical U-turn structures. This U-turn conformation exposes the Watson-Crick faces of the bases within stem-loop V (G697, A698, and C699) and makes them accessible for interaction with stem-loop I. Using chemical-shift mapping, we show that magnesium ions interact with the loop of the isolated stem-loop V and induce a conformational change that may be important for interaction with stem-loop I. This study expands our understanding of the role of U-turn motifs in RNA structure and function and provides insights into the mechanism of substrate recognition in the VS ribozyme.     

RNAMotif, an RNA secondary structure definition and search algorithm

RNA molecules fold into characteristic secondary and tertiary structures that account for their diverse functional activities. Many of these RNA structures are assembled from a collection of RNA structural motifs. These basic building blocks are used repeatedly, and in various combinations, to form different RNA types and define their unique structural and functional properties. Identification of recurring RNA structural motifs will therefore enhance our understanding of RNA structure and help associate elements of RNA structure with functional and regulatory elements. Our goal was to develop a computer program that can describe an RNA structural element of any complexity and then search any nucleotide sequence database, including the complete prokaryotic and eukaryotic genomes, for these structural elements. Here we describe in detail a new computational motif search algorithm, RNAMotif, and demonstrate its utility with some motif search examples. RNAMotif differs from other motif search tools in two important aspects: first, the structure definition language is more flexible and can specify any type of base–base interaction; second, RNAMotif provides a user controlled scoring section that can be used to add capabilities that patterns alone cannot provide.     

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.