Tetraloop Motif Research Articles

Sequence-Dependent Melting and Refolding Dynamics of RNA UNCG Tetraloops Using Temperature-Jump/Drop Infrared Spectroscopy.

Time-resolved temperature-jump/drop infrared (IR) spectroscopy has been used to measure the impact of stem base sequence on the melting and refolding dynamics of ribonucleic acid (RNA) tetraloops. A series of three 12-nucleotide RNA hairpin sequences were studied, each featuring a UACG tetraloop motif and a double-stranded stem containing four base pairs. In each case, the stem comprised three GC pairs plus a single AU base pair inserted at the closing point of the loop (RNAloop), in the middle of the stem (RNAmid), or at the stem terminus (RNAend). Results from analogous DNA tetraloop (TACG) sequences were also obtained. Inclusion of AU or AT base pairs in the stem leads to faster melting of the stem-loop structure compared to a stem sequence featuring four GC base pairs while refolding times were found to be slower, consistent with a general reduction in stem-loop stability caused by the AU/AT pair. Independent measurement of the dynamic timescales for melting and refolding of ring vibrational modes of guanine (GR) and adenine (AR) provided position-specific insight into hairpin dynamics. The GR-derived data showed that DNA sequences melted more quickly (0.5 ± 0.1 to 0.7 ± 0.1 μs at 70 °C) than analogous RNA sequences (4.3 ± 0.4 to 4.4 ± 0.3 μs at 70 °C). Position-sensitive data from the AR modes suggests that DNA hairpins begin melting from the terminal end of the stem toward the loop while RNA sequences begin melting from the loop. Refolding timescales for both RNA and DNA hairpins were found to be similar (250 ± 50 μs at 70 °C) except for RNAend and DNAloop which refolded much more slowly (746 ± 36 and 430 ± 31 μs, respectively), showing that the refolding pathway is significantly impaired by the placement of AU/AT pairs at different points in the stem. We conclude that conformational changes of analogous pairs of RNA and DNA tetraloops proceed by different mechanisms.

A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States.

An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. We isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.

Incorporating a Thiophosphate Modification into a Common RNA Tetraloop Motif Causes an Unanticipated Stability Boost

GNRA (N = A, C, G, or U; R = A or G) tetraloops are common RNA secondary structural motifs and feature a phosphate stacked atop a nucleobase. The rRNA sarcin/ricin loop (SRL) is capped by GApGA, and the phosphate p stacks on G. We recently found that regiospecific incorporation of a single dithiophosphate (PS2) but not a monothiophosphate (PSO) instead of phosphate in the backbone of RNA aptamers dramatically increases the binding affinity for their targets. In the RNA:thrombin complex, the key contribution to the 1000-fold tighter binding stems from an edge-on contact between PS2 and a phenylalanine ring. Here we investigated the consequences of replacing the SRL phosphate engaged in a face-on interaction with guanine with either PS2 or PSO for stability. We found that PS2···G and Rp-PSO···G contacts stabilize modified SRLs compared to the parent loop to unexpected levels: up to 6.3 °C in melting temperature Tm and -4.7 kcal/mol in ΔΔG°. Crystal structures demonstrate that the vertical distance to guanine for the closest sulfur is just 0.05 Å longer on average compared to that of oxygen despite the larger van der Waals radius of the former (1.80 Å for S vs 1.52 Å for O). The higher stability is enthalpy-based, and the negative charge as assessed by a neutral methylphosphonate modification plays only a minor role. Quantum mechanical/molecular mechanical calculations are supportive of favorable dispersion attraction interactions by sulfur making the dominant contribution. A stacking interaction between phosphate and guanine (SRL) or uracil (U-turn) is also found in newly classified RNA tetraloop families besides GNRA.

Effects of osmolytes and macromolecular crowders on stable GAAA tetraloops and their preference for a CG closing base pair.

Osmolytes and macromolecular crowders have the potential to influence the stability of secondary structure motifs and alter preferences for conserved nucleic acid sequences in vivo. To further understand the cellular function of RNA we observed the effects of a model osmolyte, polyethylene glycol (PEG) 200, and a model macromolecular crowding agent, PEG 8000, on the GAAA tetraloop motif. GAAA tetraloops are conserved, stable tetraloops, and are critical participants in RNA tertiary structure. They also have a thermodynamic preference for a CG closing base pair. The thermal denaturation of model hairpins containing GAAA loops was monitored using UV-Vis spectroscopy in the presence and absence of PEG 200 or PEG 8000. Both of the cosolutes tested influenced the thermodynamic preference for a CG base pair by destabilizing the loop with a CG closing base pair relative to the loop with a GC closing base pair. This result also extended to a related DNA triloop, which provides further evidence that the interactions between the loop and closing base pair are identical for the d(GCA) triloop and the GAAA tetraloop. Our results suggest that in the presence of model PEG molecules, loops with a GC closing base pair may retain some preferential interactions with the cosolutes that are lost in the presence of the CG closing base pair. These results reveal that relatively small structural changes could influence how neutral cosolutes tune the stability and function of secondary structure motifs in vivo.

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

ABSTRACTOsmolytes 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.

Synthesis, Thermodynamic Properties, and Crystal Structure of RNA Oligonucleotides Containing 5-Hydroxymethylcytosine.

5-Hydroxymethylcytosine (hm5C) is an RNA modification that has attracted significant interest because of the finding that RNA hydroxymethylation can favor mRNA translation. For insight into the mechanistic details of hm5C function to be obtained, the availability of RNAs containing this modification at defined positions that can be used for in vitro studies is highly desirable. In this work, we present an eight-step route to 5-hydroxymethylcytidine (hm5rC) phosphoramidite for solid-phase synthesis of modified RNA oligonucleotides. Furthermore, we examined the effects of hm5rC on RNA duplex stability and its impact on structure formation using the sarcin-ricin loop (SRL) motif. Thermal denaturation experiments revealed that hm5rC increases RNA duplex stability. By contrast, when cytosine within an UNCG tetraloop motif was replaced by hm5rC, the thermodynamic stability of the corresponding hairpin fold was attenuated. Importantly, incorporation of hm5rC into the SRL motif resulted in an RNA crystal structure at 0.85 Å resolution. Besides changes in the hydration pattern at the site of modification, a slight opening of the hm5rC-G pair compared to the unmodified C-G in the native structure was revealed.

On the conformational stability of the smallest RNA kissing complexes maintained through two G·C base pairs

Two identical 5′GACG3′ tetra-loop motifs with different stem sequences (called H2 and H3) are found in the 5′ end region of Moloney Murine Leukemia Virus (MMLV) genomic RNA. They play important roles in RNA dimerization and encapsidation through two identical tetra-loops (5′GACG3′) forming a loop-to-loop kissing complex, the smallest RNA kissing complex ever found in nature. We examined the effects of a loop-closing base pair as well as a stem sequence on the conformational stability of the kissing complex. UV melting analysis and gel electrophoresis were performed on eight RNA sequences mimicking the H2 and H3 hairpin tetra-loops with variation in loop-closing base pairs. Our results show that changing the loop-closing base pair from the wildtype (5′A·U3′ for H3, 5′U·A3′ for H2) to 5′G·C3’/5′C·G3′ has significant effect on the stability of the kissing complexes: the substitution to 5′C·G3′ significantly decreases both thermal and mechanical stability, while switching to the 5′G·C3′ significantly increases the mechanical stability only. The kissing complexes with the wildtype loop-closing base pairs (5′A·U3′ for H3 and 5′U·A3′ for H2) show different stability when attached to a different stem sequence (H2 stem vs. H3 stem). This suggests that not only the loop-closing base pair itself, but also the stem sequence, affects the conformational stability of the RNA kissing complex.

Simultaneous NMR characterisation of multiple minima in the free energy landscape of an RNA UUCG tetraloop.

RNA molecules in solution tend to undergo structural fluctuations of relatively large amplitude and to populate a range of different conformations some of which with low populations. It is still very challenging, however, to characterise the structures of these low populated states and to understand their functional roles. In the present study, we address this problem by using NMR residual dipolar couplings (RDCs) as structural restraints in replica-averaged metadynamics (RAM) simulations. By applying this approach to a 14-mer RNA hairpin containing the prototypical UUCG tetraloop motif, we show that it is possible to construct the free energy landscape of this RNA molecule. This free energy landscapes reveals the surprisingly rich dynamics of the UUCG tetraloop and identifies the multiple substates that exist in equilibrium owing to thermal fluctuations. The approach that we present is general and can be applied to the study of the free energy landscapes of other RNA or RNA-protein systems.

Highly sampled tetranucleotide and tetraloop motifs enable evaluation of common RNA force fields.

Recent modifications and improvements to standard nucleic acid force fields have attempted to fix problems and issues that have been observed as longer timescale simulations have become routine. Although previous work has shown the ability to fold the UUCG stem–loop structure, until now no group has attempted to quantify the performance of current force fields using highly converged structural populations of the tetraloop conformational ensemble. In this study, we report the use of multiple independent sets of multidimensional replica exchange molecular dynamics (M-REMD) simulations with different initial conditions to generate well-converged conformational ensembles for the tetranucleotides r(GACC) and r(CCCC), as well as the larger UUCG tetraloop motif. By generating what is to our knowledge the most complete RNA structure ensembles reported to date for these systems, we remove the coupling between force field errors and errors due to incomplete sampling, providing a comprehensive comparison between current top-performing MD force fields for RNA. Of the RNA force fields tested in this study, none demonstrate the ability to correctly identify the most thermodynamically stable structure for all three systems. We discuss the deficiencies present in each potential function and suggest areas where improvements can be made. The results imply that although “short” (nsec-μsec timescale) simulations may stay close to their respective experimental structures and may well reproduce experimental observables, inevitably the current force fields will populate alternative incorrect structures that are more stable than those observed via experiment.

Sense overlapping transcripts in IS1341-type transposase genes are functional non-coding RNAs in archaea

The existence of sense overlapping transcripts that share regulatory and coding information in the same genomic sequence shows an additional level of prokaryotic gene expression complexity. Here we report the discovery of ncRNAs associated with IS1341-type transposase (tnpB) genes, at the 3'-end of such elements, with examples in archaea and bacteria. Focusing on the model haloarchaeon Halobacterium salinarum NRC-1, we show the existence of sense overlapping transcripts (sotRNAs) for all its IS1341-type transposases. Publicly available transcriptome compendium show condition-dependent differential regulation between sotRNAs and their cognate genes. These sotRNAs allowed us to find a UUCA tetraloop motif that is present in other archaea (ncRNA family HgcC) and in a H. salinarum intergenic ncRNA derived from a palindrome associated transposable elements (PATE). Overexpression of one sotRNA and the PATE-derived RNA harboring the tetraloop motif improved H. salinarum growth, indicating that these ncRNAs are functional.

Tospovirus ambisense genomic RNA segments use almost complete repertoire of stable tetraloops in the intergenic region

The intergenic regions of the ambisense RNA segments of viruses from the Tospovirus genus form large extended RNA structures that regulate virus replication. Using comparative structure analysis, we show the presence of conserved alternative conformations at the apical parts of these structures. In one conformation, a branched Y-shape, the 5'-proximal hairpin arms are mostly capped by exceptionally stable tetraloop motifs. The tetraloop hairpins are folded in both virus and virus-complementary sense RNAs, and different tetraloops can functionally replace each other. Folding simulations show that the branched Y-shape structures can undergo a conformational transition to alternative extended rod-like conformations. Functional importance of both alternatives is supported by nucleotide covariations. The balanced equilibrium between alternative structures is evidenced by native gel electrophoresis of mutant RNA transcripts with shifted equilibria. The tetraloops play a role in the stability and dynamics of structures but may also be recognized by proteins involved in translation and/or replication.

NMR elucidation of the role of Mg2+ in the structure and stability of the conserved RNA motifs of the EMCV IRES element

The three dimensional solution structures of a highly conserved 16mer RNA, endowed with a classic 'GNRA' tetraloop motif, and a 17mer RNA containing a cytosine-rich heptaloop which is predicted to be a potential receptor for the former RNA, of the I-domain of Encephalomyocarditis virus IRES have been determined by NMR spectroscopy. As Mg(2+) plays an important role in the activity of the IRES, the corresponding NMR structures of the Mg(2+) bound RNA complexes have also been determined. These RNA NMR structures, 16mer (21 constraints per residue), 16mer RNA/Mg(2+) (21 constraints per residue), 17mer (17 constraints per residue) and 17mer RNA/Mg(2+) (16 constraints per residue), were calculated to a high degree of precision with low RMSDs and low clash scores represent, to the best our knowledge, the first structures of a type II picornavirus IRES. Conformational analysis of the average structure showed that the RNAs and their Mg(2+) complexes adopt characteristic A-form helical structures, stabilised by canonical and non-canonical interactions in both the stem and loop regions. The GCGA tetraloop of the 16mer folds into a standard GNRA conformation, with the structural role of A550 being in the form of a G547.A550 sheared base-pair made up of two hydrogen bonds. Further, the previously uncharacterised AACCCCA heptaloop present in the 17mer forms a compact tertiary loop motif, held together by strong π-π interactions. Analysis of the NMR structures demonstrates that the role of Mg(2+) is principally to confer enhanced stability to the RNAs whereby the tetra- and heptaloops can achieve optimum conformation for any RNA-RNA interactions which are crucial for understanding the structure-function relationship of the IRES.

High resolution de‐novo folding of hyperstable RNA tetraloops using molecular dynamics simulations

RNA's biochemical versatility arises from its ability to form tertiary structures, which are governed by the intrinsic flexibility of single‐stranded loops, weak stacking interactions, and non‐canonical hydrogen bonding patterns. However, in contrast to the numerous documented successes in protein folding, there has been no successful reports to‐date of Angstrom‐resolution folding of RNA from the unfolded state. This is caused by intrinsic errors in the force‐field parameters themselves, which have not been significantly improved in 20 years.We present preliminary results of an improved RNA force‐field which has been extensively calibrated against vapor‐pressure osmometry, ultrasonic absorption, NMR, and circular dichroism experiments on short oligonucleotides. Using replica‐exchange molecular dynamics, we are able to reversibly fold all three hyperstable tetraloop motifs (UUCG, GCAA, and CUUG) from the unfolded state. All of the unique non‐canonical interactions observed in high‐resolution structures are recapitulated by the folding simulations with an all‐atom RMSD of 1 to 3 Å. By accurately capturing the balance between the key driving forces for RNA tertiary folding, we anticipate that large‐scale RNA tertiary folding via molecular dynamics simulations will become feasible in the near future.

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.

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.

Molecular analysis of the role of IRES stem-loop V in replicative capacities and translation efficiencies of Coxsackievirus B3 mutants

Coxsackievirus B3 (CVB3) is a principal viral cause of acute myocarditis in humans and has been implicated in the pathogenesis of dilated cardiomyopathy. The natural genetic determinants of cardiovirulence for CVB3 have not been identified, although using strains engineered in the laboratory, it has been demonstrated elsewhere that, for several wild-type CB3 strains, the primary molecular determinant of cardiovirulence phenotype localizes to the 5' nontranslated region (5'NTR) and capsid. Stable RNA tetraloop motifs are found frequently in biologically active RNAs. These motifs carry out a wide variety of functions in RNA folding, in RNA-RNA and RNA-protein interactions. A great deal of knowledge about the structures and functions of tetraloop motifs has accumulated largely due to intensive theoretical, biochemical, and biophysical studies on one most frequently occurring family of tetraloop sequences, namely, the GNRA sequence, especially the GNAA sequence conserved in all enteroviruses. Here in this study, through construction of CVB3 chimeric mutants, the predicted stem loop (SL) V within the 5'NTR has been identified as important in determining viral cardiovirulence. Replication assays in HeLa cell monolayers revealed that wild-type CVB3 virus and two of the six mutants constructed here grow efficiently, whereas other mutant viruses replicate poorly. Furthermore, the in vitro translation products from these mutants and wild-type CVB3, demonstrated that the two mutants who replicate efficiently, translated at relatively equivalent amount than the wild-type. However, other mutants demonstrated a low efficiency in their production of protein when translated in a Rabbit Reticulocytes Lysats.

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

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.

Direct Revelation of Multiple Conformations in RNA by Femtosecond Dynamics

RNA structures are intrinsically dynamic, and the dynamics are important for interactions and recognition. We developed a strategy using femtosecond dynamics to quantitatively investigate the heterogeneous nature of RNA structures. The different conformations that exist in an RNA molecule are resolved based on their distinct temporal behaviors of the fluorescent decays for the incorporated fluorophore on ultrafast time scales. Using this strategy, we have probed the GNRA tetraloop motif and revealed alternative structures that have not been directly observed by other techniques. A model is proposed for the coexisting multiple conformations and the conformational dynamics for this motif.

Automated extraction and classification of RNA tertiary structure cyclic motifs

A minimum cycle basis of the tertiary structure of a large ribosomal subunit (LSU) X-ray crystal structure was analyzed. Most cycles are small, as they are composed of 3- to 5 nt, and repeated across the LSU tertiary structure. We used hierarchical clustering to quantify and classify the 4 nt cycles. One class is defined by the GNRA tetraloop motif. The inspection of the GNRA class revealed peculiar instances in sequence. First is the presence of UA, CA, UC and CC base pairs that substitute the usual sheared GA base pair. Second is the revelation of GNR(Xn)A tetraloops, where Xn is bulged out of the classical GNRA structure, and of GN/RA formed by the two strands of interior-loops. We were able to unambiguously characterize the cycle classes using base stacking and base pairing annotations. The cycles identified correspond to small and cyclic motifs that compose most of the LSU RNA tertiary structure and contribute to its thermodynamic stability. Consequently, the RNA minimum cycles could well be used as the basic elements of RNA tertiary structure prediction methods.

Residue Specific Ribose and Nucleobase Dynamics of the cUUCGg RNA Tetraloop Motif by MNMR 13C Relaxation

The dynamics of the nucleobase and the ribose moieties in a 14-nt RNA cUUCGg hairpin-loop uniformly labeled with 13C and 15N were studied by 13C spin relaxation experiments. R1, R1rho and the 13C-[1H] steady-state NOE of C6 and C1' in pyrimidine and C8 and C1' in purine residues were obtained at 298 K. The relaxation data were analyzed by the model-free formalism to yield dynamic information on timescales of pico-, nano- and milli-seconds. An axially symmetric diffusion tensor with an overall rotational correlation time tau(c) of 2.31 +/- 0.13 ns and an axial ratio of 1.35 +/- 0.02 were determined. Both findings are in agreement with hydrodynamic calculations. For the nucleobase carbons, the validity of different reported 13C chemical shift anisotropy values (Stueber, D. and Grant, D. M., 2002 J. Am. Chem. Soc. 124, 10539-10551; Fiala et al., 2000 J. Biomol. NMR 16, 291-302; Sitkoff, D. and Case, D. A., 1998 Prog. NMR Spectroscopy 32, 165-190) is discussed. The resulting dynamics are in agreement with the structural features of the cUUCGg motif in that all residues are mostly rigid (0.82 < S2 < 0.96) in both the nucleobase and the ribose moiety except for the nucleobase of U7, which is protruding into solution (S2 = 0.76). In general, ribose mobility follows nucleobase dynamics, but is less pronounced. Nucleobase dynamics resulting from the analysis of 13C relaxation rates were found to be in agreement with 15N relaxation data derived dynamic information (Akke et al., 1997 RNA 3, 702-709).