Structure and dynamics of an RNA tetraloop: a joint molecular dynamics and NMR study

Biophysical characterization and design of a minimal version of the Hoechst RNA aptamer

RNA aptamers are oligonucleotides, selected through Systematic Evolution of Ligands by EXponential Enrichment (SELEX), that can bind to specific target molecules with high affinity. One such molecule is the RNA aptamer that binds to a blue-fluorescent Hoechst dye that was modified with bulky t-Bu groups to prevent non-specific binding to DNA. This aptamer has potential for biosensor applications; however, limited information is available regarding its conformation, molecular interactions with the ligand, and binding mechanism. The study presented here aims to biophysically characterize the Hoechst RNA aptamer when complexed with the t-Bu Hoechst dye and to further optimize the RNA sequence by designing and synthesizing new sequence variants. Each variant aptamer-t-Bu Hoechst complex was evaluated through a combination of fluorescence emission, native polyacrylamide gel electrophoresis, fluorescence titration, and isothermal titration calorimetry experiments. The results were used to design a minimal version of the aptamer consisting of only 21 nucleotides. The performed study also describes a more efficient method for synthesizing the t-Bu Hoechst dye derivative. Understanding the biophysical properties of the t-Bu Hoechst dye-RNA complex lays the foundation for nuclear magnetic resonance spectroscopy studies and its potential development as a building block for an aptamer-based biosensor that can be used in medical, environmental or laboratory settings.

Investigating the Intramolecular Competition of Different RNA Binding Motifs for Neomycin B by Native Top-Down Mass Spectrometry

The ongoing search for small molecule drugs that target ribonucleic acids (RNA) is complicated by a limited understanding of the principles that govern RNA-small molecule interactions. Here we have used stoichiometry-resolved native top-down mass spectrometry (MS) to study the binding of neomycin B to small model hairpin RNAs, an unstructured RNA, and a viral RNA construct. For 15-22 nt model RNAs with hairpin structure, we found that neomycin B binding to hairpin loops relies on interactions with both the nucleobases and the 2'-OH groups, and that a simple 5' or 3' overhang can introduce an additional binding motif. For a 47 nt RNA construct derived from stem IA of the human immunodeficiency virus 1 (HIV-1) rev response element (RRE) RNA, native top-down MS identified four different binding motifs, of which the purine-rich internal loop showed the highest affinity for neomycin B. Stoichiometry-resolved binding site mapping by native top-down MS allows for a new perspective on binding specificity, and has the potential to reveal unexpected principles of small molecule binding to RNA.

Integrated NMR/Molecular Dynamics Determination of the Ensemble Conformation of a Thermodynamically Stable CUUG RNA Tetraloop

Both experimental and theoretical structure determinations of RNAs have remained challenging due to the intrinsic dynamics of RNAs. We report here an integrated nuclear magnetic resonance/molecular dynamics (NMR/MD) structure determination approach to describe the dynamic structure of the CUUG tetraloop. We show that the tetraloop undergoes substantial dynamics, leading to averaging of the experimental data. These dynamics are particularly linked to the temperature-dependent presence of a hydrogen bond within the tetraloop. Interpreting the NMR data by a single structure represents the low-temperature structure well but fails to capture all conformational states occurring at a higher temperature. We integrate MD simulations, starting from structures of CUUG tetraloops within the Protein Data Bank, with an extensive set of NMR data, and provide a structural ensemble that describes the dynamic nature of the tetraloop and the experimental NMR data well. We thus show that one of the most stable and frequently found RNA tetraloops displays substantial dynamics, warranting such an integrated structural approach.

Molecular mechanism of ATP and RNA binding to Zika virus NS3 helicase and identification of repurposed drugs using molecular dynamics simulations

Congenital Zika virus syndrome has caused a public health emergency of international concern. So far, there are no drugs available to prevent or treat the infection caused by Zika virus. The Zika virus NS3 helicase is a potential protein target for drug discovery due to its vital role in viral genome replication. NS3 helicase unwinds the viral RNA to enable the reproduction of the viral genome by the NS5 protein. NS3 helicase has two crucial binding sites; the ATP binding site and the RNA binding site. We used molecular docking and molecular dynamics (MD) simulations to study the structural behavior of Zika virus NS3 helicase in its apo form and in the presence of ATP, single-stranded RNA, and both ATP-RNA to understand their potential implications in NS3 helicase activity. Further, we have carried out virtual screening of FDA approved drugs, followed by molecular docking to identify the ATP-competitive hit molecules as probable Zika virus NS3 helicase inhibitors. The MD simulations trajectories were analyzed using normal mode analysis and principal component analysis that reveals fluctuations in the R-loop. These findings aid in understanding the molecular mechanisms of the simultaneous binding of ATP and RNA, and guide the design and discovery of new inhibitors of the Zika virus NS3 helicase as a promising drug target to treat the Zika virus infection. Communicated by Ramaswamy H. Sarma.

Hofmeister Series for Metal-Cation-RNA Interactions: The Interplay of Binding Affinity and Exchange Kinetics

A large variety of physicochemical properties involving RNA depends on the type of metal cation present in solution. In order to gain microscopic insight into the origin of these ion specific effects, we apply molecular dynamics simulations to describe the interactions of metal cations and RNA. For the three most common ion binding sites on RNA, we calculate the binding affinities and exchange rates of eight different mono- and divalent metal cations. Our results reveal that binding sites involving phosphate groups preferentially bind metal cations with high charge density (such as Mg²⁺) in inner-sphere conformations while binding sites involving N7 or O6 atoms preferentially bind cations with low charge density (such as K⁺). The binding affinity therefore follows a direct Hofmeister series at the backbone but is reversed at the nucleobases leading to a high selectivity of ion binding sites on RNA. In addition, the exchange rates for cation binding cover almost 5 orders of magnitude, leading to a vastly different time scale for the lifetimes of contact pairs. Taken together, the site-specific binding affinities and the specific lifetime of contact pairs provide the microscopic explanation of ion specific effects observed in a wide variety of macroscopic RNA properties. Finally, combining the results from atomistic simulations with extended Poisson-Boltzmann theory allows us to predict the distribution of metal cations around double-stranded RNA at finite concentrations and to reproduce the results of ion counting experiments with good accuracy.

FSATOOL: A useful tool to do the conformational sampling and trajectory analysis work for biomolecules

Reliable conformational sampling and trajectory analysis are always important to the study of the folding or binding mechanisms of biomolecules. Generally, one has to prepare many complicated parameters and follow a lot of steps to obtain the final data. The whole process is too complicated to new users. In this article, we provide a convenient and user-friendly tool that is compatible to AMBER, called fast sampling and analysis tool (FSATOOL). FSATOOL has some useful features. First and the most important, the whole work is extremely simplified into two steps, one is the fast sampling procedure and the other is the trajectory analysis procedure. Second, it contains several powerful sampling methods for the simulation on graphics process unit, including our previous mixing replica exchange molecular dynamics method. The method combines the advantages of the biased and unbiased simulations. Finally, it extracts the dominant transition pathways automatically from the folding network by Markov state model. Users do not need to do the tedious intermediate steps by hand. To illustrate the usage of FSATOOL in practice, we perform one simulation for a RNA hairpin in explicit solvent. All the results are presented. © 2019 Wiley Periodicals, Inc.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Resolving sugar puckers in RNA excited states exposes slow modes of repuckering dynamics

Recent studies have shown that RNAs exist in dynamic equilibrium with short-lived low-abundance 'excited states' that form by reshuffling base pairs in and around non-canonical motifs. These conformational states are proposed to be rich in non-canonical motifs and to play roles in the folding and regulatory functions of non-coding RNAs but their structure proves difficult to characterize given their transient nature. Here, we describe an approach for determining sugar pucker conformation in RNA excited states through nuclear magnetic resonance measurements of C1΄ and C4΄ rotating frame spin relaxation (R1ρ) in uniformly 13C/15N labeled RNA samples. Application to HIV-1 TAR exposed slow modes of sugar repuckering dynamics at the μs and ms timescale accompanying transitions between non-helical (C2΄-endo) to helical (C3΄-endo) conformations during formation of two distinct excited states. In contrast, we did not obtain any evidence for slow sugar repuckering dynamics for nucleotides in a variety of structural contexts that do not undergo non-helical to helical transitions. Our results outline a route for significantly improving the conformational characterization of RNA excited states and suggest that slow modes of repuckering dynamics gated by transient changes in secondary structure are quite common in RNA.

Telomerase RNA is more than a DNA template

The addition of telomeric DNA to chromosome ends is an essential cellular activity that compensates for the loss of genomic DNA that is due to the inability of the conventional DNA replication apparatus to duplicate the entire chromosome. The telomerase reverse transcriptase and its associated RNA bind to the very end of the telomere via a sequence in the RNA and specific protein-protein interactions. Telomerase RNA also provides the template for addition of new telomeric repeats by the reverse-transcriptase protein subunit. In addition to the template, there are 3 other conserved regions in telomerase RNA that are essential for normal telomerase activity. Here we briefly review the conserved core regions of telomerase RNA and then focus on a recent study in fission yeast that determined the function of another conserved region in telomerase RNA called the Stem Terminus Element (STE). (1) The STE is distant from the templating core of telomerase in both the linear and RNA secondary structure, but, nonetheless, affects the fidelity of telomere sequence addition and, in turn, the ability of telomere binding proteins to bind and protect chromosome ends. We will discuss possible mechanisms of STE action and the suitability of the STE as an anti-cancer target.

Elucidating the Relation between Internal Motions and Dihedral Angles in an RNA Hairpin Using Molecular Dynamics

Molecular dynamics simulations were performed to characterize the internal motions of the ribonucleic acid apical stem loop of human hepatitis B virus. The NMR relaxation rates calculated directly from the trajectory are in good agreement with the experiment. Calculated order parameters follow the experimental pattern. Order parameters lower than 0.8 are observed for nucleotides that are weakly hydrogen bonded to their base pair partner, unpaired, or part of the loop. These residues show slow decay of the internal correlation functions of their base and sugar C-H vectors. Concerted motions around backbone dihedral angles influence the amplitude of motion of the sugar and base C-H vectors. The order parameters for base C-H vectors are also affected by the fluctuation of the glycosidic dihedral angle.

What NMR Relaxation Can Tell Us about the Internal Motion of an RNA Hairpin: A Molecular Dynamics Simulation Study

Classical molecular dynamics simulations of a 14-mer UUCG RNA hairpin are performed to study its conformational dynamics and corresponding NMR relaxation parameters. The direct calculation of the relaxation rates from the trajectory yields good agreement with experiment, indicating the validity of the theoretical model. Various ways to provide a link between theory and experiment are considered, including the "model-free approach" of Lipari and Szabo and Gaussian axial fluctuation model of Brüschweilwer. It is studied if the underlying assumptions of these approaches are satisfied in the case of a flexible RNA hairpin. Being consistent with the analysis of the NMR experiments, Lipari-Szabo fits of the first 100 or 1000 ps of the internal correlation functions lead to a nice agreement between calculated and experimental order parameters and internal correlation times. Finally, the relation between NMR order parameters and the underlying internal motion of the RNA hairpin is discussed in detail. A principal component analysis reveals that the principal motions of the molecule account only partially for the measured NMR order parameters, because the latter are insensitive to internal dynamics occurring on a nanosecond time scale due to molecular tumbling.

Influence of Nitroxide Spin Labels on RNA Structure: A Molecular Dynamics Simulation Study

Pulsed electron double resonance (PELDOR) experiments on oligonucleotides provide a distance ruler that allows the measurement of nanometer distances accurately. The technique requires attachment of nitroxide spin labels to the nucleotides, which may possibly perturb its conformation. To study to what extent nitroxide spin labels may affect RNA structure, all-atom molecular dynamics simulations in explicit solvent are performed for six double-labeled RNA duplexes. A new parametrization of the force field for the nitroxide spin label is developed, which leads to intramolecular distances that are in good agreement with experimental results. Comparison of the results for spin-labeled and unlabeled RNA reveals that the conformational effect of the spin label depends significantly on whether the spin label is attached to the major or the minor groove of RNA. While major-groove spin labeling may to some extent affect the conformation of nearby base pairs, minor-groove spin labeling has the advantage of mostly preserving the RNA conformation.

Structural Insights Reveal the Dynamics of the Repeating r(CAG) Transcript Found in Huntington's Disease (HD) and Spinocerebellar Ataxias (SCAs)

In humans, neurodegenerative disorders such as Huntington's disease (HD) and many spinocerebellar ataxias (SCAs) have been found to be associated with CAG trinucleotide repeat expansion. An important RNA-mediated mechanism that causes these diseases involves the binding of the splicing regulator protein MBNL1 (Muscleblind-like 1 protein) to expanded r(CAG) repeats. Moreover, mutant huntingtin protein translated from expanded r(CAG) also yields toxic effects. To discern the role of mutant RNA in these diseases, it is essential to gather information about its structure. Detailed insight into the different structures and conformations adopted by these mutant transcripts is vital for developing therapeutics targeting them. Here, we report the crystal structure of an RNA model with a r(CAG) motif, which is complemented by an NMR-based solution structure obtained from restrained Molecular Dynamics (rMD) simulation studies. Crystal structure data of the RNA model resolved at 2.3 Å reveals non-canonical pairing of adenine in 5´-CAG/3´-GAC motif samples in different syn and anti conformations. The overall RNA structure has helical parameters intermediate to the A- and B-forms of nucleic acids due to the global widening of major grooves and base-pair preferences near internal AA loops. The comprehension of structural behaviour by studying the spectral features and the dynamics also supports the flexible nature of the r(CAG) motif.

Theory of nonrigid rotational motion applied to NMR relaxation in RNA

Solution NMR spectroscopy can elucidate many features of the structure and dynamics of macromolecules, yet relaxation measurements, the most common source of experimental information on dynamics, can sample only certain ranges of dynamic rates. A complete characterization of motion of a macromolecule thus requires the introduction of complementary experimental approaches. Solid-state NMR spectroscopy successfully probes the time scale of nanoseconds to microseconds, a dynamic window where solution NMR results have been deficient, and probes conditions where the averaging effects of rotational diffusion of the molecule are absent. Combining the results of the two distinct techniques within a single framework provides greater insight into dynamics, but this task requires the common interpretation of results recorded under very different experimental conditions. Herein, we provide a unified description of dynamics that is robust to the presence of large-scale conformational exchange, where the diffusion tensor of the molecule varies on a time scale comparable to rotational diffusion in solution. We apply this methodology to the HIV-1 TAR RNA molecule, where conformational rearrangements are both substantial and functionally important. The formalism described herein is of greater generality than earlier combined solid-state/solution NMR interpretations, if detailed molecular structures are available, and can offer a more complete description of RNA dynamics than either solution or solid-state NMR spectroscopy alone.

Nonequilibrium NMR Methods for Monitoring Protein and RNA Folding

The review introduces to time-resolved NMR spectroscopic investigations of the kinetics of protein and RNA folding. The description of the experimental investigations is discussed in the context of possible kinetic folding pathways showing the extent of information that can be gained from the various kinetic experiments. The review introduces to four different methods to initiate folding reactions in connection with time-resolved NMR experiments and discusses examples of refolding of the model proteinα-lactalbumin and of bistable RNAs.

Identifying and searching for conserved RNA localisation signals

RNA localisation is an important mode of delivering proteins to their site of function. Cis-acting signals within the RNAs, which can be thought of as zip-codes, determine the site of localisation. There are few examples of fully characterised RNA signals, but the signals are thought to be defined through a combination of primary, secondary, and tertiary structures. In this chapter, we describe a selection of computational methods for predicting RNA secondary structure, identifying localisation signals, and searching for similar localisation signals on a genome-wide scale. The chapter is aimed at the biologist rather than presenting the details of each of the individual methods.

Slow exchange model of nonrigid rotational motion in RNA for combined solid-state and solution NMR studies

Functional RNA molecules are conformationally dynamic and sample a multitude of dynamic modes over a wide range of frequencies. Thus, a comprehensive description of RNA dynamics requires the inclusion of a broad range of motions across multiple dynamic rates which must be derived from multiple spectroscopies. Here we describe a slow conformational exchange theoretical approach to combining the description of local motions in RNA that occur in the nanosecond to microsecond window and are detected by solid-state NMR with nonrigid rotational motion of the HIV-1 transactivation response element (TAR) RNA in solution as observed by solution NMR. This theoretical model unifies the experimental results generated by solution and solid-state NMR and provides a comprehensive view of the dynamics of HIV-1 TAR RNA, a well-known paradigm of an RNA where function requires extensive conformational rearrangements. This methodology provides a quantitative atomic level view of the amplitudes and rates of the local and collective displacements of the TAR RNA molecule and provides directly motional parameters for the conformational capture hypothesis of this classical RNA-ligand interaction.

Loop and stem dynamics during RNA hairpin folding and unfolding

2-Aminopurine (2AP) is a fluorescent adenine analog that probes mainly base stacking in nucleic acids. We labeled the loop or the stem of the RNA hairpin gacUACGguc with 2AP to study folding thermodynamics and kinetics at both loci. Thermal melts and fast laser temperature jumps detected by 2AP fluorescence monitored the stability and folding/unfolding kinetics. The observed thermodynamic and kinetic traces of the stem and loop mutants, though strikingly different at a first glance, can be fitted to the same free-energy landscape. The differences between the two probe locations arise because base stacking decreases upon unfolding in the stem, whereas it increases in the loop. We conclude that 2AP is a conservative adenine substitution for mapping out the contributions of different RNA structural elements to the overall folding process. Molecular dynamics (MD) totaling 0.6 μsec were performed to look at the conformations populated by the RNA at different temperatures. The combined experimental data, and MD simulations lead us to propose a minimal four-state free-energy landscape for the RNA hairpin. Analysis of this landscape shows that a sequential folding model is a good approximation for the full folding dynamics. The frayed state formed initially from the native state is a heterogeneous ensemble of structures whose stem is frayed either from the end or from the loop.

Network models for molecular kinetics and their initial applications to human health

Molecular kinetics underlies all biological phenomena and, like many other biological processes, may best be understood in terms of networks. These networks, called Markov state models (MSMs), are typically built from physical simulations. Thus, they are capable of quantitative prediction of experiments and can also provide an intuition for complex conformational changes. Their primary application has been to protein folding; however, these technologies and the insights they yield are transferable. For example, MSMs have already proved useful in understanding human diseases, such as protein misfolding and aggregation in Alzheimer's disease.

Referencing strategy for the direct comparison of nuclear magnetic resonance and molecular dynamics motional parameters in RNA

Nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations are both techniques that can be used to characterize the structural dynamics of biomolecules and their underlying time scales. Comparison of relaxation parameters obtained through each methodology allows for cross validation of techniques and for complementarity in the analysis of dynamics. Here we present a combined NMR/MD study of the dynamics of HIV-1 transactivation response (TAR) RNA. We compute relaxation constants (R(1), R(2), and NOE) and model-free parameters (S(2) and tau) from a 65 ns molecular dynamics (MD) trajectory and compare them with the respective parameters measured in a domain-elongation NMR experiment. Using the elongated domain as the frame of reference for all computed parameters allows for a direct comparison between experiment and simulation. We see good agreement for many parameters and gain further insight into the nature of the local and global dynamics of TAR, which are found to be quite complex, spanning multiple time scales. For the few cases where agreement is poor, comparison of the dynamical parameters provides insight into the limits of each technique. We suggest a frequency-matching procedure that yields an upper bound for the time scale of dynamics to which the NMR relaxation experiment is sensitive.

Fast folding of an RNA tetraloop on a rugged energy landscape detected by a stacking-sensitive probe

We investigate the microsecond-timescale kinetics of the RNA hairpin ga*cUUCGguc. The fluorescent nucleotide 2-aminopurine (a*) reports mainly on base stacking. Ten kinetic traces and the temperature denaturation curve are globally fitted to four-state models of the free-energy surface. In the best-fitting sequential model, the hairpin unfolds over successively larger barriers in at least three stages: stem fraying and increased base-stacking fluctuations; concerted loss of hydrogen bonding and partial unstacking; and additional unstacking of single strands at the highest temperatures. Parallel and trap models also provide adequate fits: such pathways probably also play a role in the complete free-energy surface of the hairpin. To interpret the model states structurally, 200 ns of molecular dynamics, including six temperature-jump simulations, were run. Although the sampling is by no means comprehensive, five different states were identified using hydrogen bonding and base stacking as reaction coordinates. The four to five states required to explain the experiments or simulations set a lower limit on the complexity of this small RNA hairpin's energy landscape.

RNA phosphodiester backbone dynamics of a perdeuterated cUUCGg tetraloop RNA from phosphorus-31 NMR relaxation analysis

We have analyzed the relaxation properties of all (31)P nuclei in an RNA cUUCGg tetraloop model hairpin at proton magnetic field strengths of 300, 600 and 900 MHz in solution. Significant H, P dipolar contributions to R (1) and R (2) relaxation are observed in a protonated RNA sample at 600 MHz. These contributions can be suppressed using a perdeuterated RNA sample. In order to interpret the (31)P relaxation data (R (1), R (2)), we measured the (31)P chemical shift anisotropy (CSA) by solid-state NMR spectroscopy under various salt and hydration conditions. A value of 178.5 ppm for the (31)P CSA in the static state (S (2) = 1) could be determined. In order to obtain information about fast time scale dynamics we performed a modelfree analysis on the basis of our relaxation data. The results show that subnanosecond dynamics detected around the phosphodiester backbone are more pronounced than the dynamics detected for the ribofuranosyl and nucleobase moieties of the individual nucleotides (Duchardt and Schwalbe, J Biomol NMR 32:295-308, 2005; Ferner et al., Nucleic Acids Res 36:1928-1940, 2008). Furthermore, the dynamics of the individual phosphate groups seem to be correlated to the 5' neighbouring nucleobases.

A new computational approach for mechanical folding kinetics of RNA hairpins

Based on an ensemble of kinetically accessible conformations, we propose a new analytical model for RNA folding kinetics. The model gives populational kinetics, kinetic rates, transition states, and pathways from the rate matrix. Applications of the new kinetic model to mechanical folding of RNA hairpins such as trans-activation-responsive RNA reveal distinct kinetic behaviors in different force regimes, from zero force to forces much stronger than the critical force for the folding-unfolding transition. In the absence of force or a low force, folding can be initiated (nucleated) at any position by forming the first base stack and there exist many pathways for the folding process. In contrast, for a higher force, the folding/unfolding would predominantly proceed along a single zipping/unzipping pathway. Studies for different hairpin-forming sequences indicate that depending on the nucleotide sequence, a kinetic intermediate can emerge in the low force regime but disappear in high force regime, and a new kinetic intermediate, which is absent in the low and high force regimes, can emerge in the medium force range. Variations of the force lead to changes in folding cooperativity and rate-limiting steps. The predicted network of pathways for trans-activation-responsive RNA suggests two parallel dominant pathways. The rate-limiting folding steps (at f = 8 pN) are the formation of specific basepairs that are 2-4 basepairs away from the loop. At a higher force (f = 11 pN), the folding rate is controlled by the formation of the bulge loop. The predicted rates and transition states are in good agreement with the experimental data for a broad force regime.

NMR and MD studies of the temperature-dependent dynamics of RNA YNMG-tetraloops

In a combined NMR/MD study, the temperature-dependent changes in the conformation of two members of the RNA YNMG-tetraloop motif (cUUCGg and uCACGg) have been investigated at temperatures of 298, 317 and 325 K. The two members have considerable different thermal stability and biological functions. In order to address these differences, the combined NMR/MD study was performed. The large temperature range represents a challenge for both, NMR relaxation analysis (consistent choice of effective bond length and CSA parameter) and all-atom MD simulation with explicit solvent (necessity to rescale the temperature). A convincing agreement of experiment and theory is found. Employing a principle component analysis of the MD trajectories, the conformational distribution of both hairpins at various temperatures is investigated. The ground state conformation and dynamics of the two tetraloops are indeed found to be very similar. Furthermore, both systems are initially destabilized by a loss of the stacking interactions between the first and the third nucleobase in the loop region. While the global fold is still preserved, this initiation of unfolding is already observed at 317 K for the uCACGg hairpin but at a significantly higher temperature for the cUUCGg hairpin.

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.

iRED analysis of TAR RNA reveals motional coupling, long-range correlations, and a dynamical hinge

The HIV-1 transactivation response RNA element (TAR), which is essential to the lifecycle of the virus, has been suggested, based on NMR and hydrodynamic measurements, to undergo substantial, collective, structural dynamics that are important for its function. To deal with the significant coupling between overall diffusional rotation and internal motion expected to exist in TAR, here we utilize an isotropic reorientational eigenmode dynamics analysis of simulated molecular trajectories to obtain a detailed description of TAR dynamics and an accurately quantified pattern of dynamical correlations. The analysis demonstrates the inseparability of internal and overall motional modes, confirms the existence and reveals the nature of collective domain dynamics, and additionally reveals that the hinge for these motions is centered on residues U23, C24, and C41. Results also indicate the existence of long-range communication between the loop and the core of the RNA, and between the loop and the bulge. Additionally, the isotropic reorientational eigenmode dynamics analysis explains, from a dynamical perspective, several existing biochemical mutational studies and suggests new mutations for future structural dynamics studies.

Structure and thermodynamics of a conserved U2 snRNA domain from yeast and human

The spliceosome is a dynamic ribonucleoprotein complex responsible for the removal of intron sequences from pre-messenger RNA. The highly conserved 5' end of the U2 small nuclear RNA (snRNA) makes key base-pairing interactions with the intron branch point sequence and U6 snRNA. U2 stem I, a stem-loop located in the 5' region of U2, has been implicated in spliceosome assembly and may modulate the folding of the U2 and U6 snRNAs in the spliceosome active site. Here we present the NMR structures of U2 stem I from human and Saccharomyces cerevisiae. These sequences represent the two major classes of U2 stem I, distinguished by the identity of tandem wobble pairs (UU/UU in yeast and CA/GU in human) and the presence of post-transcriptional modifications (four 2'-O-methyl groups and two pseudouracils in human). The structures reveal that the UU/UU and CA/GU tandem wobble pairs are nearly isosteric. The tandem wobble pairs separate two thermodynamically distinct regions of Watson-Crick base pairs, with the modified nucleotides in human stem I conferring a significant increase in stability. We hypothesize that the separate thermodynamic stabilities of U2 stem I exist to allow the structure to transition through different folded conformations during spliceosome assembly and catalysis.

RNA localization signals: deciphering the message with bioinformatics

mRNA localization is an important posttranscriptional method of targeting proteins to their site of function. The sorting of transcripts to their correct intracellular destination is achieved by a number of mechanisms, including selective degradation or transport by molecular motors along the cytoskeleton. In all cases, this involves mRNA localization signals, or so called zip codes, being recognized by trans-acting cellular factors. In a few cases, primary sequence motifs for RNA localization can be identified, but in general, localization signals operate at the level of secondary (2D) and tertiary (3D) structure. This inevitably means that searching for localization signal motifs is a complex task requiring specialist knowledge of bioinformatics. Furthermore, the publications describing these searching methods tend to be aimed at the bioinformatics community. In this review, we have surveyed the major tools for folding, comparing, and searching potential mRNA localization signals in transcripts or across genomes. Our aim is to provide an overview for biologists, who lack specialist computer and bioinformatics training, of the major RNA bioinformatics tools that are available. The examples provided are focused on mRNA localization signals and RNA stem-loop structures, however, these tools can be applied to the study of any RNA signals.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.