UUCG Tetraloop Research Articles

NMR analysis of a loop-bulge structure of UUCGA pentaloop

Although structures of many RNA loops, such as GNRA and UNCG tetraloops, were well known, it is still possible to find more RNA structures. In the present study, solution structure of an RNA fragment having UUCGA pentaloop was analyzed by NMR spectroscopy. It was found that the UUCG tetraloop is formed and the adenosine residue at the 3ʹ side of the tetraloop is bulged out. The characteristic motif of the loop-bulge structure has also been found in other RNAs including CUUGU and CUGGC pentaloops. Along with the recently found T-hairpin structure with a UUUGAUU loop, in which UUUGA pentaloop and UU bulge are formed, the loop-bulge structures can be categorized as an RNA motif and it may be called as the integrated structure loop, I-loop.

Deep Boosted Molecular Dynamics: Accelerating Molecular Simulations with Gaussian Boost Potentials Generated Using Probabilistic Bayesian Deep Neural Network.

We have developed a new deep boosted molecular dynamics (DBMD) method. Probabilistic Bayesian neural network models were implemented to construct boost potentials that exhibit Gaussian distribution with minimized anharmonicity, thereby allowing for accurate energetic reweighting and enhanced sampling of molecular simulations. DBMD was demonstrated on model systems of alanine dipeptide and the fast-folding protein and RNA structures. For alanine dipeptide, 30 ns DBMD simulations captured up to 83-125 times more backbone dihedral transitions than 1 μs conventional molecular dynamics (cMD) simulations and were able to accurately reproduce the original free energy profiles. Moreover, DBMD sampled multiple folding and unfolding events within 300 ns simulations of the chignolin model protein and identified low-energy conformational states comparable to previous simulation findings. Finally, DBMD captured a general folding pathway of three hairpin RNAs with the GCAA, GAAA, and UUCG tetraloops. Based on a deep learning neural network, DBMD provides a powerful and generally applicable approach to boosting biomolecular simulations. DBMD is available with open source in OpenMM at https://github.com/MiaoLab20/DBMD/.

Structural insights into recognition of SL4, the UUCG stem-loop, of human U1 snRNA by the ubiquitin-like domain, including the C-terminal tail in the SF3A1 subunit of U2 snRNP.

The pre-spliceosomal complex involves interactions between U1 and U2 snRNPs, where a ubiquitin-like domain (ULD) of SF3A1, a component of U2 snRNP, binds to the stem-loop 4 (SL4; the UUCG tetraloop) of U1 snRNA in U1 snRNP. Here, we reported the 1.80 Å crystal structure of human SF3A1 ULD (ULDSF3A1) complexed with SL4. The structural part of ULDSF3A1 (res. 704-785) adopts a typical β-grasp fold with a topology of β1-β2-α1-310a-β3-β4-310b-β5, closely resembling that of ubiquitin, except for the length and structure of the β1/β2 loop. A patch on the surface formed by three ULDSF3A1-specific residues, Lys756 (β3), Phe763 (β4), and Lys765 (following β4), contacts the canonical UUCG tetraloop structure. In contrast, the directly following C-terminal tail composed of 786KERGGRKK793 was essentially stretched. The main or side chains of all the residues interacted with the major groove of the stem helix; the RGG residues adopted a peculiar conformation for RNA recognition. These findings were confirmed by mutational studies using bio-layer interferometry. Collectively, a unique combination of the β-grasp fold and the C-terminal tail constituting ULDSF3A1 is required for the SL4-specific binding. This interaction mode also suggests that putative post-translational modifications, including ubiquitination in ULDSF3A1, directly inhibit SL4 binding.

Structures and Dynamics of DNA Mini-Dumbbells Are Force Field Dependent.

Flexible nucleic acid structures can be challenging to accurately resolve with currently available experimental structural determination techniques. As an alternative, molecular dynamics (MD) simulations can provide a window into understanding the unique dynamics and population distributions of these biomolecules. Previously, molecular dynamics simulations of noncanonical (non-duplex) nucleic acids have proven difficult to accurately model. With a new influx of improved nucleic acid force fields, achieving an in-depth understanding of the dynamics of flexible nucleic acid structures may be achievable. In this project, currently available nucleic acid force fields are evaluated using a flexible yet stable model system: the DNA mini-dumbbell. Prior to MD simulations, nuclear magnetic resonance (NMR) re-refinement was accomplished using improved refinement techniques in explicit solvent to yield DNA mini-dumbbell structures with better agreement between the newly determined PDB snapshots, with the NMR data itself, as well as the unrestrained simulation data. Starting from newly determined structures, a total aggregate of over 800 μs of production data between 2 DNA mini-dumbbell sequences and 8 force fields was collected to compare to these newly refined structures. The force fields tested spanned from traditional Amber force fields: bsc0, bsc1, OL15, and OL21 to Charmm force fields: Charmm36 and the Drude polarizable force field, as well as force fields from independent developers: Tumuc1 and CuFix/NBFix. The results indicated slight variations not only between the different force fields but also between the sequences as well. Given our previous experiences with high populations of potentially anomalous structures in RNA UUCG tetraloops and in various tetranucleotides, we expected the mini-dumbbell system to be challenging to accurately model. Surprisingly, many of the recently developed force fields generated structures in good agreement with experiments. Yet, each of the force fields provided a different distribution of potentially anomalous structures.

Toward Convergence in Folding Simulations of RNA Tetraloops: Comparison of Enhanced Sampling Techniques and Effects of Force Field Modifications.

Atomistic molecular dynamics simulations represent an established technique for investigation of RNA structural dynamics. Despite continuous development, contemporary RNA simulations still suffer from suboptimal accuracy of empirical potentials (force fields, ffs) and sampling limitations. Development of efficient enhanced sampling techniques is important for two reasons. First, they allow us to overcome the sampling limitations, and second, they can be used to quantify ff imbalances provided they reach a sufficient convergence. Here, we study two RNA tetraloops (TLs), namely the GAGA and UUCG motifs. We perform extensive folding simulations and calculate folding free energies (ΔGfold°) with the aim to compare different enhanced sampling techniques and to test several modifications of the nonbonded terms extending the AMBER OL3 RNA ff. We demonstrate that replica-exchange solute tempering (REST2) simulations with 12-16 replicas do not show any sign of convergence even when extended to a timescale of 120 μs per replica. However, the combination of REST2 with well-tempered metadynamics (ST-MetaD) achieves good convergence on a timescale of 5-10 μs per replica, improving the sampling efficiency by at least 2 orders of magnitude. Effects of ff modifications on ΔGfold° energies were initially explored by the reweighting approach and then validated by new simulations. We tested several manually prepared variants of the gHBfix potential which improve stability of the native state of both TLs by ∼2 kcal/mol. This is sufficient to conveniently stabilize the folded GAGA TL while the UUCG TL still remains under-stabilized. Appropriate adjustment of van der Waals parameters for C-H···O5' base-phosphate interaction may further stabilize the native states of both TLs by ∼0.6 kcal/mol.

Sequence-specific RNA recognition by an RGG motif connects U1 and U2 snRNP for spliceosome assembly

In mammals, the structural basis for the interaction between U1 and U2 small nuclear ribonucleoproteins (snRNPs) during the early steps of splicing is still elusive. The binding of the ubiquitin-like (UBL) domain of SF3A1 to the stem-loop 4 of U1 snRNP (U1-SL4) contributes to this interaction. Here, we determined the 3D structure of the complex between the UBL of SF3A1 and U1-SL4 RNA. Our crystallography, NMR spectroscopy, and cross-linking mass spectrometry data show that SF3A1-UBL recognizes, sequence specifically, the GCG/CGC RNA stem and the apical UUCG tetraloop of U1-SL4. In vitro and invivo mutational analyses support the observed intermolecular contacts and demonstrate that the carboxyl-terminal arginine-glycine-glycine-arginine (RGGR) motif of SF3A1-UBL binds sequence specifically by inserting into the RNA major groove. Thus, the characterization of the SF3A1-UBL/U1-SL4 complex expands the repertoire of RNA binding domains and reveals the capacity of RGG/RG motifs to bind RNA in a sequence-specific manner.

UUCG RNA Tetraloop as a Formidable Force-Field Challenge for MD Simulations

Explicit solvent atomistic molecular dynamics (MD) simulations represent an established technique to study structural dynamics of RNA molecules and an important complement for diverse experimental methods. However, performance of molecular mechanical (MM) force fields (ff's) remains far from satisfactory even after decades of development, as apparent from a problematic structural description of some important RNA motifs. Actually, some of the smallest RNA molecules belong to the most challenging systems for MD simulations and, among them, the UUCG tetraloop is saliently difficult. We report a detailed analysis of UUCG MD simulations, depicting the sequence of events leading to the loss of the UUCG native state during MD simulations. The total amount of MD simulation data analyzed in this work is close to 1.3 ms. We identify molecular interactions, backbone conformations, and substates that are involved in the process. Then, we unravel specific ff deficiencies using diverse quantum mechanical/molecular mechanical (QM/MM) and QM calculations. Comparison between the MM and QM methods shows discrepancies in the description of the 5'-flanking phosphate moiety and both signature sugar-base interactions. Our work indicates that poor behavior of the UUCG tetraloop in simulations is a complex issue that cannot be attributed to one dominant and straightforwardly correctable factor. Instead, there is a concerted effect of multiple ff inaccuracies that are coupled and amplifying each other. We attempted to improve the simulation behavior by some carefully tailored interventions, but the results were still far from satisfactory, underlying the difficulties in development of accurate nucleic acid ff's.

Reducing the measurement time of exact NOEs by non-uniform sampling.

We have previously reported on the measurement of exact NOEs (eNOEs), which yield a wealth of additional information in comparison to conventional NOEs. We have used these eNOEs in a variety of applications, including calculating high-resolution structures of proteins and RNA molecules. The collection of eNOEs is challenging, however, due to the need to measure a NOESY buildup series consisting of typically four NOESY spectra with varying mixing times in a single measurement session. While the 2D version can be completed in a few days, a fully sampled 3D-NOESY buildup series can take 10days or more to acquire. This can be both expensive as well as problematic in the case of samples that are not stable over such a long period of time. One potential method to significantly decrease the required measurement time of eNOEs is to use non-uniform sampling (NUS) to decrease the number of points measured in the indirect dimensions. The effect of NUS on the extremely tight distance restraints extracted from eNOEs may be very pronounced. Therefore, we investigated the fidelity of eNOEs measured from three test cases at decreasing NUS densities: the 18.4kDa protein human Pin1, the 4.1kDa WW domain of Pin1 (both in 3D), and a 4.6kDa 14mer RNA UUCG tetraloop (2D). Our results show that NUS imparted negligible error on the eNOE distances derived from good quality data down to 10% sampling for all three cases, but there is a noticeable decrease in the eNOE yield that is dependent upon the underlying sparsity, and thus complexity, of the sample. For Pin1, this transition occurred at roughly 40% while for the WW domain and the UUCG tetraloop it occurred at lower NUS densities of 20% and 10%, respectively. We rationalized these numbers through reconstruction simulations under various conditions. The extent of this loss depends upon the number of scans taken as well as the number of peaks to be reconstructed. Based on these findings, we have created guidelines for choosing an optimal NUS density depending on the number of peaks needed to be reconstructed in the densest region of a 2D or 3D NOESY spectrum.

Integrating NMR and simulations reveals motions in the UUCG tetraloop.

We provide an atomic-level description of the structure and dynamics of the UUCG RNA stem–loop by combining molecular dynamics simulations with experimental data. The integration of simulations with exact nuclear Overhauser enhancements data allowed us to characterize two distinct states of this molecule. The most stable conformation corresponds to the consensus three-dimensional structure. The second state is characterized by the absence of the peculiar non-Watson–Crick interactions in the loop region. By using machine learning techniques we identify a set of experimental measurements that are most sensitive to the presence of non-native states. We find that although our MD ensemble, as well as the consensus UUCG tetraloop structures, are in good agreement with experiments, there are remaining discrepancies. Together, our results show that (i) the MD simulation overstabilize a non-native loop conformation, (ii) eNOE data support its presence with a population of ≈10% and (iii) the structural interpretation of experimental data for dynamic RNAs is highly complex, even for a simple model system such as the UUCG tetraloop.

Integrating NMR and Simulations Reveals Motions in the UUCG Tetraloop

Integrating NMR and Simulations Reveals Motions in the UUCG Tetraloop

Conformational adaptation of UNCG loops upon crowding.

If the A-form helix is the major structural motif found in RNA, the loops that cap them constitute the second most important family of motifs. Among those, two are overrepresented, GNRA and UNCG tetraloops. Recent surveys of RNA structures deposited in the PDB show that GNRA and UNCG tetraloops can adopt tertiary folds that are very different from their canonical conformations, characterized by the presence of a U-turn of a Z-turn, respectively. Crystallographic data from both a lariat-capping (LC) ribozyme and a group II intron ribozyme reveal that a given UUCG tetraloop can adopt a distinct fold depending on its structural environment. Specifically, when the crystal packing applies relaxed constraints on the loop, the canonical Z-turn conformation is observed. In contrast, a highly packed environment induces “squashing” of the tetraloop by distorting its sugar-phosphate backbone in a specific way that expels the first and fourth nucleobases out of the loop, and falls in van der Waals distance of the last base pair of the helix, taking the place of the pair formed between the first and fourth residues in Z-turn loops. The biological relevance of our observations is supported by the presence of similarly deformed loops in the highly packed environment of the ribosome and in a complex between a dsRNA and a RNase III. The finding that Z-turn loops change conformation under higher molecular packing suggests that, in addition to their demonstrated role in stabilizing RNA folding, they may contribute to the three-dimensional structure of RNA by mediating tertiary interactions with distal residues.

High-resolution small RNA structures from exact nuclear Overhauser enhancement measurements without additional restraints

RNA not only translates the genetic code into proteins, but also carries out important cellular functions. Understanding such functions requires knowledge of the structure and dynamics at atomic resolution. Almost half of the published RNA structures have been solved by nuclear magnetic resonance (NMR). However, as a result of severe resonance overlap and low proton density, high-resolution RNA structures are rarely obtained from nuclear Overhauser enhancement (NOE) data alone. Instead, additional semi-empirical restraints and labor-intensive techniques are required for structural averages, while there are only a few experimentally derived ensembles representing dynamics. Here we show that our exact NOE (eNOE) based structure determination protocol is able to define a 14-mer UUCG tetraloop structure at high resolution without other restraints. Additionally, we use eNOEs to calculate a two-state structure, which samples its conformational space. The protocol may open an avenue to obtain high-resolution structures of small RNA of unprecedented accuracy with moderate experimental efforts.

Automated Force-Field Parametrization Guided by Multisystem Ensemble Averages

RNA structure and dynamics play a fundamental role in many cellular processes such as gene expression inhibition, splicing and catalysis. Molecular dynamics is a computational tool that can be used to investigate RNA structure and dynamics at atomistic resolution. However, its capability to predict and explain experimental data is limited by the accuracy of the employed potential energy functions, also known as force fields. Recent works have shown that state-of-the-art force fields could predict unphysical RNA conformations that are not in agreement with experiments. The emerging strategy to overcome these limitations is to complement molecular dynamics with experimental data included as restraints. In a recent work, we suggested a maximum-entropy based method to enforce solution experiments in molecular dynamics simulations. Importantly, this approach reduces the risk of overfitting by simultaneously adapting the force-field corrections to multiple systems. We here push this idea further and develop a general scheme to fit arbitrary force-field parameters given a set of ensemble averages. Such quantities can range from NMR data, such as 3J couplings or NOE, to native state populations. The key feature of our method is the possibility to concurrently combine ensemble averages from multiple systems into a unique error function to be minimized, drastically enhancing corrections’ transferability. The method is applied to the difficult case of GAGA and UUCG tetraloops for which we are able to maximize their native state population by refining torsional potentials alone.

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.

Noncanonical α/γ Backbone Conformations in RNA and the Accuracy of Their Description by the AMBER Force Field.

The sugar-phosphate backbone of RNA can exist in diverse rotameric substates, giving RNA molecules enormous conformational variability. The most frequent noncanonical backbone conformation in RNA is α/γ = t/t, which is derived from the canonical backbone by a crankshaft motion and largely preserves the standard geometry of the RNA duplex. A similar conformation also exists in DNA, where it has been extensively studied and shown to be involved in DNA-protein interactions. However, the function of the α/γ = t/t conformation in RNA is poorly understood. Here, we present molecular dynamics simulations of several prototypical RNA structures obtained from X-ray and NMR experiments, including canonical and mismatched RNA duplexes, UUCG and GAGA tetraloops, Loop E, the sarcin-ricin loop, a parallel guanine quadruplex, and a viral pseudoknot. The stability of various noncanonical α/γ backbone conformations was analyzed with two AMBER force fields, ff99bsc0χOL3 and ff99bsc0χOL3 with the recent εζOL1 and βOL1 corrections for DNA. Although some α/γ substates were stable with seemingly well-described equilibria, many were unstable in our simulations. Notably, the most frequent noncanonical conformer α/γ = t/t was unstable in both tested force fields. Possible reasons for this instability are discussed. Our work reveals a potentially important artifact in RNA force fields and highlights a need for further force field refinement.

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.

Structural fidelity and NMR relaxation analysis in a prototype RNA hairpin.

RNA hairpins are widespread and very stable motifs that contribute decisively to RNA folding and biological function. The GTP1G2C3A4C5U6U7C8G9G10U11G12C13C14 construct (with a central UUCG tetraloop) has been extensively studied by solution NMR, and offers and excellent opportunity to evaluate the structure and dynamical description afforded by molecular dynamics (MD) simulations. Here, we compare average structural parameters and NMR relaxation rates estimated from a series of multiple independent explicit solvent MD simulations using the two most recent RNA AMBER force fields (ff99 and ff10). Predicted overall tumbling times are ∼20% faster than those inferred from analysis of NMR data and follow the same trend when temperature and ionic strength is varied. The Watson–Crick stem and the “canonical” UUCG loop structure are maintained in most simulations including the characteristic syn conformation along the glycosidic bond of G9, although some key hydrogen bonds in the loop are partially disrupted. Our analysis pinpoints G9–G10 backbone conformations as a locus of discrepancies between experiment and simulation. In general the results for the more recent force-field parameters (ff10) are closer to experiment than those for the older ones (ff99). This work provides a comprehensive and detailed comparison of state of the art MD simulations against a wide variety of solution NMR measurements.

Molecular dynamics study of the controlled destabilization of an RNA hairpin structure by a covalently attached azobenzene switch.

As shown in recent experimental studies, photoswitches like azobenzene can act as efficient regulators of the folding and unfolding of DNA and RNA duplexes. Here we explore the details of the conformational changes induced by azobenzene attachment, focusing upon a small 14-mer RNA hairpin structure. The azobenzene chromophore is covalently bound to the stem region adjacent to a UUCG tetraloop which is known to represent a particularly stable structure. Since the characteristic time scale of conformational changes exceeds the nanosecond scale (and by far exceeds the ultrafast time scale of trans-to-cis photoswitching), equilibrium simulations using enhanced sampling by replica exchange molecular dynamics (REMD) are employed to investigate the influence of trans versus cis azobenzene attachment on the stability of the hairpin. We report on the analysis of fluctuations and conformational landscapes, along with calculations of relative melting temperatures. The simulations are found to reproduce certain experimentally predicted trends for azobenzene-modified RNA; in particular, both trans and cis conformers have a destabilizing effect. This effect is significantly enhanced for the cis conformer, even though the latter tends to flip out of the double-stranded stem region.

RNA does the folding dance of twist, turn, stack

RNA molecules have long been known to adopt unusual conformations. The first crystal structure of an RNA was that of yeast tRNAphe (1, 2), which provided a portent of what was to come. With more structures available now from X-ray crystallography and NMR, that promise has been met beyond expectations. The RNAs that are folded in molecular dynamics simulations by Chen and Garcia (3) are tetraloops (RNA hairpin loops of only four nucleotides) that are among the most unusual small structures known today. Although or because they are unusual with atypical thermal stability, these loops are also ubiquitous in longer RNA strands [rUUCG (4); rGNRA and rCUUG (5)]; they were first identified as the most abundant hairpins in prokaryotic rRNA. The function of these loops varies: the GNRA tetraloops (where N is any base and R is a purine) is often found in RNA:RNA tertiary interactions to anchor a global fold. The UUCG tetraloop is solitary, with no known interactions with RNA or proteins: it is thought to nucleate folding and so direct in vivo RNA folding patterns. These little hairpins have been the objects of considerable in vitro experimental structural characterization by NMR and X-ray crystallography, and are found to consistently form a single unique structure. The rUUCG tetraloop in particular has been subjected to many thermodynamic investigations to measure the sources of its unusual stability (6). The rCUUG hairpin is structurally the most flexible, because it tucks its first U into the minor groove but its second U is relatively unconstrained. Among the three tetraloops, CUUG has a propensity for significant conformational exchange (7). Although the three tetraloops used in the Chen and Garcia (3) computational folding experiments are all unusually stable, they achieve their stability by very different means. Their structural diversity makes them excellent candidates for folding studies, because no one mechanism dominates their ability to adopt their unique structures.