Hierarchy of RNA functional dynamics

NMR characterization of RNA small molecule interactions

Exciting discoveries of naturally occurring ligand-sensing and disease-linked noncoding RNAs have promoted significant interests in understanding RNA-small molecule interactions. NMR spectroscopy is a powerful tool for characterizing intermolecular interactions. In this review, we describe protocols and approaches for applying NMR spectroscopy to investigate interactions between RNA and small molecules. We review protocols for RNA sample preparation, methods for identifying RNA-binding small molecules, approaches for mapping RNA-small molecule interactions, determining complex structures, and characterizing binding kinetics. We hope this review will provide a guideline to streamline NMR applications in studying RNA-small molecule interactions, facilitating both basic mechanistic understandings of RNA functions and translational efforts in developing RNA-targeted therapeutics.

Understanding the Contributions of Conformational Changes, Thermodynamics, and Kinetics of RNA-Small Molecule Interactions

The implication of RNA in multiple cellular processes beyond protein coding has revitalized interest in the development of small molecules for therapeutically targeting RNA and for further probing its cellular biology. However, the process of rationally designing such small molecule probes is hampered by the paucity of information about fundamental molecular recognition principles of RNA. In this Review, we summarize two important and often underappreciated aspects of RNA-small molecule recognition: RNA conformational dynamics and the biophysical properties of interactions of small molecules with RNA, specifically thermodynamics and kinetics. While conformational flexibility is often said to impede RNA ligand development, the ability of small molecules to influence the RNA conformational landscape can have a significant effect on the cellular functions of RNA. An analysis of the conformational landscape of RNA and the interactions of individual conformations with ligands can thus guide the development of new small molecule probes, which needs to be investigated further. Additionally, while it is common practice to quantify the binding affinities ( K_a or K_d) of small molecules for biomacromolecules as a measure of their activity, further biophysical characterization of their interaction can provide a deeper understanding. Studies that focus on the thermodynamic and kinetic parameters for interaction between RNA and ligands are next discussed. Finally, this Review provides the reader with a perspective on how such in-depth analysis of biophysical characteristics of the interaction of RNA and small molecules can impact our understanding of these interactions and how they will benefit the future design of small molecule probes.

Coming Together: RNAs and Proteins Assemble under the Single-Molecule Fluorescence Microscope

RNAs, across their numerous classes, often work in concert with proteins in RNA-protein complexes (RNPs) to execute critical cellular functions. Ensemble-averaging methods have been instrumental in revealing many important aspects of these RNA-protein interactions, yet are insufficiently sensitive to much of the dynamics at the heart of RNP function. Single-molecule fluorescence microscopy (SMFM) offers complementary, versatile tools to probe RNP conformational and compositional changes in detail. In this review, we first outline the basic principles of SMFM as applied to RNPs, describing key considerations for labeling, imaging, and quantitative analysis. We then sample applications of in vitro and in vivo single-molecule visualization using the case studies of pre-messenger RNA (mRNA) splicing and RNA silencing, respectively. After discussing specific insights single-molecule fluorescence methods have yielded, we briefly review recent developments in the field and highlight areas of anticipated growth.

3D based on 2D: Calculating helix angles and stacking patterns using forgi 2.0, an RNA Python library centered on secondary structure elements

We present forgi, a Python library to analyze the tertiary structure of RNA secondary structure elements. Our representation of an RNA molecule is centered on secondary structure elements (stems, bulges and loops). By fitting a cylinder to the helix axis, these elements are carried over into a coarse-grained 3D structure representation. Integration with Biopython allows for handling of all-atom 3D information. forgi can deal with a variety of file formats including dotbracket strings, PDB and MMCIF files. We can handle modified residues, missing residues, cofold and multifold structures as well as nucleotide numbers starting at arbitrary positions. We apply this library to the study of stacking helices in junctions and pseudoknots and investigate how far stacking helices in solved experimental structures can divert from coaxial geometries.

3D based on 2D: Calculating helix angles and stacking patterns using forgi 2.0, an RNA Python library centered on secondary structure elements

We present forgi, a Python library to analyze the tertiary structure of RNA secondary structure elements. Our representation of an RNA molecule is centered on secondary structure elements (stems, bulges and loops). By fitting a cylinder to the helix axis, these elements are carried over into a coarse-grained 3D structure representation. Integration with Biopython allows for handling of all-atom 3D information. forgi can deal with a variety of file formats including dotbracket strings, PDB and MMCIF files. We can handle modified residues, missing residues, cofold and multifold structures as well as nucleotide numbers starting at arbitrary positions. We apply this library to the study of stacking helices in junctions and pseudo knots and investigate how far stacking helices in solved experimental structures can divert from coaxial geometries.

Predicting Site-Binding Modes of Ions and Water to Nucleic Acids Using Molecular Solvation Theory

Site binding of ions and water shapes nucleic acids folding, dynamics, and biological function, complementing the more diffuse, nonspecific "territorial" ion binding. Unlike territorial binding, prediction of site-specific binding to nucleic acids remains an unsolved challenge in computational biophysics. This work presents a new toolset based on the 3D-RISM molecular solvation theory and topological analysis that predicts cation and water site binding to nucleic acids. 3D-RISM is shown to accurately capture alkali cations and water binding to the central channel, transversal loops, and grooves of the Oxytricha nova's telomeres' G-quadruplex ( Oxy-GQ), in agreement with high-resolution crystallographic data. To improve the computed cation occupancy along the Oxy-GQ central channel, it was necessary to refine and validate new cation-oxygen parameters using structural and thermodynamic data available for crown ethers and ion channels. This single set of parameters that describes both localized and delocalized binding to various biological systems is used to gain insight into cation occupancy along the Oxy-GQ channel under various salt conditions. The paper concludes with prospects for extending the method to predict divalent cation binding to nucleic acids. This work advances the forefront of theoretical methods able to provide predictive insight into ion atmosphere effects on nucleic acids function.

Deconstructing internal ribosome entry site elements: an update of structural motifs and functional divergences

Beyond the general cap-dependent translation initiation, eukaryotic organisms use alternative mechanisms to initiate protein synthesis. Internal ribosome entry site (IRES) elements are cis-acting RNA regions that promote internal initiation of translation using a cap-independent mechanism. However, their lack of primary sequence and secondary RNA structure conservation, as well as the diversity of host factor requirement to recruit the ribosomal subunits, suggest distinct types of IRES elements. In spite of this heterogeneity, conserved motifs preserve sequences impacting on RNA structure and RNA–protein interactions important for IRES-driven translation. This conservation brings the question of whether IRES elements could consist of basic building blocks, which upon evolutionary selection result in functional elements with different properties. Although RNA-binding proteins (RBPs) perform a crucial role in the assembly of ribonucleoprotein complexes, the versatility and plasticity of RNA molecules, together with their high flexibility and dynamism, determines formation of macromolecular complexes in response to different signals. These properties rely on the presence of short RNA motifs, which operate as modular entities, and suggest that decomposition of IRES elements in short modules could help to understand the different mechanisms driven by these regulatory elements. Here we will review evidence suggesting that model IRES elements consist of the combination of short modules, providing sites of interaction for ribosome subunits, eIFs and RBPs, with implications for definition of criteria to identify novel IRES-like elements genome wide.

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

Modeling Structure, Stability, and Flexibility of Double-Stranded RNAs in Salt Solutions

Double-stranded (ds) RNAs play essential roles in many processes of cell metabolism. The knowledge of three-dimensional (3D) structure, stability, and flexibility of dsRNAs in salt solutions is important for understanding their biological functions. In this work, we further developed our previously proposed coarse-grained model to predict 3D structure, stability, and flexibility for dsRNAs in monovalent and divalent ion solutions through involving an implicit structure-based electrostatic potential. The model can make reliable predictions for 3D structures of extensive dsRNAs with/without bulge/internal loops from their sequences, and the involvement of the structure-based electrostatic potential and corresponding ion condition can improve the predictions for 3D structures of dsRNAs in ion solutions. Furthermore, the model can make good predictions for thermal stability for extensive dsRNAs over the wide range of monovalent/divalent ion concentrations, and our analyses show that the thermally unfolding pathway of dsRNA is generally dependent on its length as well as its sequence. In addition, the model was employed to examine the salt-dependent flexibility of a dsRNA helix, and the calculated salt-dependent persistence lengths are in good accordance with experiments.

Structural Flexibility Enables Alternative Maturation, ARGONAUTE Sorting and Activities of miR168, a Global Gene Silencing Regulator in Plants

In eukaryotes, the RNase-III Dicer often produces length/sequence microRNA (miRNA) variants, called "isomiRs", owing to intrinsic structural/sequence determinants of the miRNA precursors (pre-miRNAs). In this study, we combined biophysics, genetics and biochemistry approaches to study Arabidopsis miR168, the key feedback regulator of central plant silencing effector protein ARGONAUTE1 (AGO1). We identified a motif conserved among plant pre-miR168 orthologs, which enables flexible internal base-pairing underlying at least three metastable structural configurations. These configurations promote alternative, accurate Dicer cleavage events generating length and structural isomiR168 variants with distinctive AGO sorting properties and modes of action. Among these isomiR168s, a duplex with a 22-nt guide strand exhibits strikingly preferential affinity for AGO10, the closest AGO1 paralog. The 22-nt miR168-AGO10 complex antagonizes AGO1 accumulation in part via "transitive RNAi", a silencing-amplification process, to maintain appropriate AGO1 cellular homeostasis. Furthermore, we found that the tombusviral P19 silencing-suppressor protein displays markedly weaker affinity for the 22-nt form among its isomiR168 cargoes, thereby promoting AGO10-directed suppression of AGO1-mediated antiviral silencing. Taken together, these findings indicate that structural flexibility, a previously overlooked property of pre-miRNAs, considerably increases the versatility and regulatory potential of individual MIRNA genes, and that some pathogens might have evolved the capacity or mechanisms to usurp this property.

The nearest neighbor and next nearest neighbor effects on the thermodynamic and kinetic properties of RNA base pair

The thermodynamic and kinetic parameters of an RNA base pair with different nearest and next nearest neighbors were obtained through long-time molecular dynamics simulation of the opening-closing switch process of the base pair near its melting temperature. The results indicate that thermodynamic parameters of GC base pair are dependent on the nearest neighbor base pair, and the next nearest neighbor base pair has little effect, which validated the nearest-neighbor model. The closing and opening rates of the GC base pair also showed nearest neighbor dependences. At certain temperature, the closing and opening rates of the GC pair with nearest neighbor AU is larger than that with the nearest neighbor GC, and the next nearest neighbor plays little role. The free energy landscape of the GC base pair with the nearest neighbor GC is rougher than that with nearest neighbor AU.

How Messenger RNA and Nascent Chain Sequences Regulate Translation Elongation

Translation elongation is a highly coordinated, multistep, multifactor process that ensures accurate and efficient addition of amino acids to a growing nascent-peptide chain encoded in the sequence of translated messenger RNA (mRNA). Although translation elongation is heavily regulated by external factors, there is clear evidence that mRNA and nascent-peptide sequences control elongation dynamics, determining both the sequence and structure of synthesized proteins. Advances in methods have driven experiments that revealed the basic mechanisms of elongation as well as the mechanisms of regulation by mRNA and nascent-peptide sequences. In this review, we highlight how mRNA and nascent-peptide elements manipulate the translation machinery to alter the dynamics and pathway of elongation.

Extending the Applicability of Exact Nuclear Overhauser Enhancements to Large Proteins and RNA

Distance-dependent nuclear Overhauser enhancements (NOEs) are one of the most popular and important experimental restraints for calculating NMR structures. Despite this, they are mostly employed as semiquantitative upper distance bounds, and this discards the wealth of information that is encoded in the cross-relaxation rate constant. Information that is lost includes exact distances between protons and dynamics that occur on the sub-millisecond timescale. Our recently introduced exact measurement of the NOE (eNOE) requires little additional experimental effort relative to other NMR observables. So far, we have used eNOEs to calculate multistate ensembles of proteins up to approximately 150 residues. Here, we briefly revisit eNOE methodology and present two new directions for the use of eNOEs: applications to large proteins and RNA.

A guide to large-scale RNA sample preparation

RNA is becoming more important as an increasing number of functions, both regulatory and enzymatic, are being discovered on a daily basis. As the RNA boom has just begun, most techniques are still in development and changes occur frequently. To understand RNA functions, revealing the structure of RNA is of utmost importance, which requires sample preparation. We review the latest methods to produce and purify a variation of RNA molecules for different purposes with the main focus on structural biology and biophysics. We present a guide aimed at identifying the most suitable method for your RNA and your biological question and highlighting the advantages of different methods. Graphical abstract In this review we present different methods for large-scale production and purification of RNAs for structural and biophysical studies.

The intrinsic flexibility of the aptamer targeting the ribosomal protein S8 is a key factor for the molecular recognition

BACKGROUND

Aptamers are RNA/DNA biomolecules representing an emerging class of protein interactors and regulators. Despite the growing interest in these molecules, current understanding of chemical-physical basis of their target recognition is limited. Recently, the characterization of the aptamer targeting the protein-S8 has suggested that flexibility plays important functional roles. We investigated the structural versatility of the S8-aptamer by molecular dynamics simulations.

METHODS

Five different simulations have been conducted by varying starting structures and temperatures.

RESULTS

The simulation of S8-aptamer complex provides a dynamic view of the contacts occurring at the complex interface. The simulation of the aptamer in ligand-free state indicates that its central region is intrinsically endowed with a remarkable flexibility. Nevertheless, none of the trajectory structures adopts the structure observed in the S8-aptamer complex. The aptamer ligand-bound is very rigid in the simulation carried out at 300 K. A structural transition of this state, providing insights into the aptamer-protein recognition process, is observed in a simulation carried out at 400 K. These data indicate that a key event in the binding is linked to the widening of the central region of the aptamer. Particularly relevant is switch of the A26 base from its ligand-free state to a location that allows the G13-C28 base-pairing.

CONCLUSIONS

Intrinsic flexibility of the aptamer is essential for partner recognition. Present data indicate that S8 recognizes the aptamer through an induced-fit rather than a population-shift mechanism.

GENERAL SIGNIFICANCE

The present study provides deeper understanding of the structural basis of the structural versatility of aptamers.

New algorithms to represent complex pseudoknotted RNA structures in dot-bracket notation

Motivation

Understanding the formation, architecture and roles of pseudoknots in RNA structures are one of the most difficult challenges in RNA computational biology and structural bioinformatics. Methods predicting pseudoknots typically perform this with poor accuracy, often despite experimental data incorporation. Existing bioinformatic approaches differ in terms of pseudoknots' recognition and revealing their nature. A few ways of pseudoknot classification exist, most common ones refer to a genus or order. Following the latter one, we propose new algorithms that identify pseudoknots in RNA structure provided in BPSEQ format, determine their order and encode in dot-bracket-letter notation. The proposed encoding aims to illustrate the hierarchy of RNA folding.

Results

New algorithms are based on dynamic programming and hybrid (combining exhaustive search and random walk) approaches. They evolved from elementary algorithm implemented within the workflow of RNA FRABASE 1.0, our database of RNA structure fragments. They use different scoring functions to rank dissimilar dot-bracket representations of RNA structure. Computational experiments show an advantage of new methods over the others, especially for large RNA structures.

Availability and implementation

Presented algorithms have been implemented as new functionality of RNApdbee webserver and are ready to use at http://rnapdbee.cs.put.poznan.pl.

Contact

mszachniuk@cs.put.poznan.pl.

Supplementary information

Supplementary data are available at Bioinformatics online.

Nucleic acids: function and potential for abiogenesis

The emergence of functional cooperation between the three main classes of biomolecules - nucleic acids, peptides and lipids - defines life at the molecular level. However, how such mutually interdependent molecular systems emerged from prebiotic chemistry remains a mystery. A key hypothesis, formulated by Crick, Orgel and Woese over 40 year ago, posits that early life must have been simpler. Specifically, it proposed that an early primordial biology lacked proteins and DNA but instead relied on RNA as the key biopolymer responsible not just for genetic information storage and propagation, but also for catalysis, i.e. metabolism. Indeed, there is compelling evidence for such an 'RNA world', notably in the structure of the ribosome as a likely molecular fossil from that time. Nevertheless, one might justifiably ask whether RNA alone would be up to the task. From a purely chemical perspective, RNA is a molecule of rather uniform composition with all four bases comprising organic heterocycles of similar size and comparable polarity and pK a values. Thus, RNA molecules cover a much narrower range of steric, electronic and physicochemical properties than, e.g. the 20 amino acid side-chains of proteins. Herein we will examine the functional potential of RNA (and other nucleic acids) with respect to self-replication, catalysis and assembly into simple protocellular entities.

AMOEBA Polarizable Atomic Multipole Force Field for Nucleic Acids

The AMOEBA polarizable atomic multipole force field for nucleic acids is presented. Valence and electrostatic parameters were determined from high-level quantum mechanical data, including structures, conformational energy, and electrostatic potentials, of nucleotide model compounds. Previously derived parameters for the phosphate group and nucleobases were incorporated. A total of over 35 μs of condensed-phase molecular dynamics simulations of DNA and RNA molecules in aqueous solution and crystal lattice were performed to validate and refine the force field. The solution and/or crystal structures of DNA B-form duplexes, RNA duplexes, and hairpins were captured with an average root-mean-squared deviation from NMR structures below or around 2.0 Å. Structural details, such as base pairing and stacking, sugar puckering, backbone and χ-torsion angles, groove geometries, and crystal packing interfaces, agreed well with NMR and/or X-ray. The interconversion between A- and B-form DNAs was observed in ethanol-water mixtures at 328 K. Crystal lattices of B- and Z-form DNA and A-form RNA were examined with simulations. For the RNA tetraloop, single strand tetramers, and HIV TAR with 29 residues, the simulated conformational states, 3 J-coupling, nuclear Overhauser effect, and residual dipolar coupling data were compared with NMR results. Starting from a totally unstacked/unfolding state, the rCAAU tetranucleotide was folded into A-form-like structures during ∼1 μs molecular dynamics simulations.

Ribosome-dependent conformational flexibility changes and RNA dynamics of IRES domains revealed by differential SHAPE

Internal ribosome entry site (IRES) elements are RNA regions that recruit the translation machinery internally. Here we investigated the conformational changes and RNA dynamics of a picornavirus IRES upon incubation with distinct ribosomal fractions. Differential SHAPE analysis of the free RNA showed that nucleotides reaching the final conformation on long timescales were placed at domains 4 and 5, while candidates for long-range interactions were located in domain 3. Salt-washed ribosomes induced a fast RNA local flexibility modification of domains 2 and 3, while ribosome-associated factors changed domains 4 and 5. Consistent with this, modeling of the three-dimensional RNA structure indicated that incubation of the IRES with native ribosomes induced a local rearrangement of the apical region of domain 3, and a reorientation of domains 4 and 5. Furthermore, specific motifs within domains 2 and 3 showed a decreased flexibility upon incubation with ribosomal subunits in vitro, and presence of the IRES enhanced mRNA association to the ribosomal subunits in whole cell lysates. The finding that RNA modules can provide direct IRES-ribosome interaction suggests that linking these motifs to additional sequences able to recruit trans-acting factors could be useful to design synthetic IRESs with novel activities.

Chemical synthesis and NMR spectroscopy of long stable isotope labelled RNA

We showcase the high potential of the 2'-cyanoethoxymethyl (CEM) methodology to synthesize RNAs with naturally occurring modified residues carrying stable isotope (SI) labels for NMR spectroscopic applications. The method was applied to synthesize RNAs with sizes ranging between 60 to 80 nucleotides. The presented approach gives the possibility to selectively modify larger RNAs (>60 nucleotides) with atom-specifically ¹³C/¹⁵N-labelled building blocks. The method harbors the unique potential to address structural as well as dynamic features of these RNAs with NMR spectroscopy but also using other biophysical methods, such as mass spectrometry (MS), or small angle neutron/X-ray scattering (SANS, SAXS).

Exploring biomolecular energy landscapes

The potential energy landscape perspective provides both a conceptual and a computational framework for predicting, understanding and designing molecular properties. In this Feature Article, we highlight some recent advances that greatly facilitate structure prediction and analysis of global thermodynamics and kinetics in proteins and nucleic acids. The geometry optimisation procedures, on which these calculations are based, can be accelerated significantly using local rigidification of selected degrees of freedom, and through implementations on graphics processing units. Results of progressive local rigidification are first summarised for trpzip1, including a systematic analysis of the heat capacity and rearrangement rates. Benchmarks for all the essential optimisation procedures are then provided for a variety of proteins. Applications are then illustrated from a study of how mutation affects the energy landscape for a coiled-coil protein, and for transitions in helix morphology for a DNA duplex. Both systems exhibit an intrinsically multifunnel landscape, with the potential to act as biomolecular switches.

Probing the excited-state chemical shifts and exchange parameters by nitrogen-decoupled amide proton chemical exchange saturation transfer (HNdec-CEST)

CEST-NMR spectroscopy is a powerful tool for probing the conformational dynamics of macromolecules. We present a HN^dec-CEST experiment that simplifies the relaxation matrix, reduces fitting parameters, and enhances signal resolution. Importantly, fitting of HN^dec-CEST profiles enables robust extraction of exchange rates as well as excited-state chemical shifts and populations.

RNA force field with accuracy comparable to state-of-the-art protein force fields

The complex and often highly dynamic 3D structures of RNA molecules are central to their diverse cellular functions. Molecular dynamics (MD) simulations have played a major role in characterizing the structure and dynamics of proteins, but the physical models (“force fields”) used for simulating nucleic acids are substantially less accurate overall than those used in protein simulations, creating a major challenge for MD studies of RNA. Here, we report an RNA force field capable of describing the structural and thermodynamic properties of RNA molecules with accuracy comparable to state-of-the-art protein force fields. This force field should facilitate the use of MD simulation as a tool for the study of biologically significant RNA molecules and protein–RNA complexes.

Molecular dynamics (MD) simulation has become a powerful tool for characterizing at an atomic level of detail the conformational changes undergone by proteins. The application of such simulations to RNA structures, however, has proven more challenging, due in large part to the fact that the physical models (“force fields”) available for MD simulations of RNA molecules are substantially less accurate in many respects than those currently available for proteins. Here, we introduce an extensive revision of a widely used RNA force field in which the parameters have been modified, based on quantum mechanical calculations and existing experimental information, to more accurately reflect the fundamental forces that stabilize RNA structures. We evaluate these revised parameters through long-timescale MD simulations of a set of RNA molecules that covers a wide range of structural complexity, including single-stranded RNAs, RNA duplexes, RNA hairpins, and riboswitches. The structural and thermodynamic properties measured in these simulations exhibited dramatically improved agreement with experimentally determined values. Based on the comparisons we performed, this RNA force field appears to achieve a level of accuracy comparable to that of state-of-the-art protein force fields, thus significantly advancing the utility of MD simulation as a tool for elucidating the structural dynamics and function of RNA molecules and RNA-containing biological assemblies.

Structural basis for the dimerization of Nab2 generated by RNA binding provides insight into its contribution to both poly(A) tail length determination and transcript compaction in Saccharomyces cerevisiae

In Saccharomyces cerevisiae generation of export-competent mRNPs terminates the nuclear phase of the gene expression pathway and facilitates transport to the cytoplasm for translation. Nab2 functions in this process to control both mRNP compaction that facilitates movement through nuclear pore complexes and the length of transcript poly(A) tails. Nab2 has a modular structure that includes seven CCCH Zn fingers that bind to A-rich RNAs and fingers 5–7 are critical for these functions. Here, we demonstrate, using both biophysical and structural methods, that binding A₁₁G RNA induces dimerization of Zn fingers 5–7 mediated by the novel spatial arrangement of the fingers promoting each RNA chain binding two protein chains. The dimerization of Nab2 induced by RNA binding provides a basis for understanding its function in both poly(A) tail length regulation and in the compaction of mature transcripts to facilitate nuclear export.

Distributed biotin-streptavidin transcription roadblocks for mapping cotranscriptional RNA folding

RNA folding during transcription directs an order of folding that can determine RNA structure and function. However, the experimental study of cotranscriptional RNA folding has been limited by the lack of easily approachable methods that can interrogate nascent RNA structure at nucleotide resolution. To address this, we previously developed cotranscriptional selective 2΄-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq) to simultaneously probe all intermediate RNA transcripts during transcription by stalling elongation complexes at catalytically dead EcoRIE111Q roadblocks. While effective, the distribution of elongation complexes using EcoRIE111Q requires laborious PCR using many different oligonucleotides for each sequence analyzed. Here, we improve the broad applicability of cotranscriptional SHAPE-Seq by developing a sequence-independent biotin-streptavidin (SAv) roadblocking strategy that simplifies the preparation of roadblocking DNA templates. We first determine the properties of biotin-SAv roadblocks. We then show that randomly distributed biotin-SAv roadblocks can be used in cotranscriptional SHAPE-Seq experiments to identify the same RNA structural transitions related to a riboswitch decision-making process that we previously identified using EcoRIE111Q. Lastly, we find that EcoRIE111Q maps nascent RNA structure to specific transcript lengths more precisely than biotin-SAv and propose guidelines to leverage the complementary strengths of each transcription roadblock in cotranscriptional SHAPE-Seq.

RNA structure prediction: from 2D to 3D

We summarize different levels of RNA structure prediction, from classical 2D structure to extended secondary structure and motif-based research toward 3D structure prediction of RNA. We outline the importance of classical secondary structure during all those levels of structure prediction.

Synthesis and incorporation of 13C-labeled DNA building blocks to probe structural dynamics of DNA by NMR

We report the synthesis of atom-specifically 13C-modified building blocks that can be incorporated into DNA via solid phase synthesis to facilitate investigations on structural and dynamic features via NMR spectroscopy. In detail, 6-13C-modified pyrimidine and 8-13C purine DNA phosphoramidites were synthesized and incorporated into a polypurine tract DNA/RNA hybrid duplex to showcase the facile resonance assignment using site-specific labeling. We also addressed micro- to millisecond dynamics in the mini-cTAR DNA. This DNA is involved in the HIV replication cycle and our data points toward an exchange process in the lower stem of the hairpin that is up-regulated in the presence of the HIV-1 nucleocapsid protein 7. As another example, we picked a G-quadruplex that was earlier shown to exist in two folds. Using site-specific 8-13C-2'deoxyguanosine labeling we were able to verify the slow exchange between the two forms on the chemical shift time scale. In a real-time NMR experiment the re-equilibration of the fold distribution after a T-jump could be monitored yielding a rate of 0.012 min-1. Finally, we used 13C-ZZ-exchange spectroscopy to characterize the kinetics between two stacked X-conformers of a Holliday junction mimic. At 25°C, the refolding process was found to occur at a forward rate constant of 3.1 s-1 and with a backward rate constant of 10.6 s-1.

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.

Inter-helical conformational preferences of HIV-1 TAR-RNA from maximum occurrence analysis of NMR data and molecular dynamics simulations

Detecting conformational heterogeneity in biological macromolecules is a key for the understanding of their biological function. We here provide a comparison between two independent approaches to assess conformational heterogeneity: molecular dynamics simulations, performed without inclusion of any experimental data, and maximum occurrence (MaxOcc) distribution over the topologically available conformational space. The latter only reflects the extent of the averaging and identifies regions which are most compliant with the experimentally measured NMR Residual Dipolar Couplings (RDCs). The analysis was performed for the HIV-1 TAR RNA, consisting of two helical domains connected by a flexible bulge junction, for which four sets of RDCs were available as well as an 8.2 μs all-atom molecular dynamics simulation. A sample and select approach was previously applied to extract from the molecular dynamics trajectory conformational ensembles in agreement with the four sets of RDCs. The MaxOcc analysis performed here identifies the most likely sampled region in the conformational space of the system which, strikingly, overlaps well with the structures independently sampled in the molecular dynamics calculations and even better with the RDC selected ensemble.

Cotranscriptional folding of a riboswitch at nucleotide resolution

RNAs can begin to fold immediately as they emerge from RNA polymerase. During cotranscriptional folding, interactions between nascent RNAs and ligands are able to direct the formation of alternative RNA structures, a feature exploited by noncoding RNAs called riboswitches to make gene-regulatory decisions. Despite their importance, cotranscriptional folding pathways have yet to be uncovered with sufficient resolution to reveal how cotranscriptional folding governs RNA structure and function. To access cotranscriptional folding at nucleotide resolution, we extended selective 2'-hydroxyl acylation analyzed by primer-extension sequencing (SHAPE-seq) to measure structural information of nascent RNAs during transcription. Using cotranscriptional SHAPE-seq, we determined how the cotranscriptional folding pathway of the Bacillus cereus crcB fluoride riboswitch undergoes a ligand-dependent bifurcation that delays or promotes terminator formation via a series of coordinated structural transitions. Our results directly link cotranscriptional RNA folding to a genetic decision and establish a framework for cotranscriptional analysis of RNA structure at nucleotide resolution.

Excited States of Nucleic Acids Probed by Proton Relaxation Dispersion NMR Spectroscopy

In this work an improved stable isotope labeling protocol for nucleic acids is introduced. The novel building blocks eliminate/minimize homonuclear (13) C and (1) H scalar couplings thus allowing proton relaxation dispersion (RD) experiments to report accurately on the chemical exchange of nucleic acids. Using site-specific (2) H and (13) C labeling, spin topologies are introduced into DNA and RNA that make (1) H relaxation dispersion experiments applicable in a straightforward manner. The novel RNA/DNA building blocks were successfully incorporated into two nucleic acids. The A-site RNA was previously shown to undergo a two site exchange process in the micro- to millisecond time regime. Using proton relaxation dispersion experiments the exchange parameters determined earlier could be recapitulated, thus validating the proposed approach. We further investigated the dynamics of the cTAR DNA, a DNA transcript that is involved in the viral replication cycle of HIV-1. Again, an exchange process could be characterized and quantified. This shows the general applicablility of the novel labeling scheme for (1) H RD experiments of nucleic acids.

Shortening the HIV-1 TAR RNA Bulge by a Single Nucleotide Preserves Motional Modes over a Broad Range of Time Scales

Helix-junction-helix (HJH) motifs are flexible building blocks of RNA architecture that help define the orientation and dynamics of helical domains. They are also frequently involved in adaptive recognition of proteins and small molecules and in the formation of tertiary contacts. Here, we use a battery of nuclear magnetic resonance techniques to examine how deleting a single bulge residue (C24) from the human immunodeficiency virus type 1 (HIV-1) transactivation response element (TAR) trinucleotide bulge (U23-C24-U25) affects dynamics over a broad range of time scales. Shortening the bulge has an effect on picosecond-to-nanosecond interhelical and local bulge dynamics similar to that casued by increasing the Mg(2+) and Na(+) concentration, whereby a preexisting two-state equilibrium in TAR is shifted away from a bent flexible conformation toward a coaxial conformation, in which all three bulge residues are flipped out and flexible. Surprisingly, the point deletion minimally affects microsecond-to-millisecond conformational exchange directed toward two low-populated and short-lived excited conformational states that form through reshuffling of bases pairs throughout TAR. The mutant does, however, adopt a slightly different excited conformational state on the millisecond time scale, in which U23 is intrahelical, mimicking the expected conformation of residue C24 in the excited conformational state of wild-type TAR. Thus, minor changes in HJH topology preserve motional modes in RNA occurring over the picosecond-to-millisecond time scales but alter the relative populations of the sampled states or cause subtle changes in their conformational features.

Thirty-five years of research into ribozymes and nucleic acid catalysis: where do we stand today?

Since the discovery of the first catalytic RNA in 1981, the field of ribozyme research has developed from the discovery of catalytic RNA motifs in nature and the elucidation of their structures and catalytic mechanisms, into a field of engineering and design towards application in diagnostics, molecular biology and medicine. Owing to the development of powerful protocols for selection of nucleic acid catalysts with a desired functionality from random libraries, the spectrum of nucleic acid supported reactions has greatly enlarged, and importantly, ribozymes have been accompanied by DNAzymes. Current areas of research are the engineering of allosteric ribozymes for artificial regulation of gene expression, the design of ribozymes and DNAzymes for medicinal and environmental diagnostics, and the demonstration of RNA world relevant ribozyme activities. In addition, new catalytic motifs or novel genomic locations of known motifs continue to be discovered in all branches of life by the help of high-throughput bioinformatic approaches. Understanding the biological role of the catalytic RNA motifs widely distributed in diverse genetic contexts belongs to the big challenges of future RNA research.

Nucleic acid polymeric properties and electrostatics: Directly comparing theory and simulation with experiment

Nucleic acids are biopolymers that carry genetic information and are also involved in various gene regulation functions such as gene silencing and protein translation. Because of their negatively charged backbones, nucleic acids are polyelectrolytes. To adequately understand nucleic acid folding and function, we need to properly describe its i) polymer/polyelectrolyte properties and ii) associating ion atmosphere. While various theories and simulation models have been developed to describe nucleic acids and the ions around them, many of these theories/simulations have not been well evaluated due to complexities in comparison with experiment. In this review, I discuss some recent experiments that have been strategically designed for straightforward comparison with theories and simulation models. Such data serve as excellent benchmarks to identify limitations in prevailing theories and simulation parameters.

Saturation of recognition elements blocks evolution of new tRNA identities

Understanding the principles that led to the current complexity of the genetic code is a central question in evolution. Expansion of the genetic code required the selection of new transfer RNAs (tRNAs) with specific recognition signals that allowed them to be matured, modified, aminoacylated, and processed by the ribosome without compromising the fidelity or efficiency of protein synthesis. We show that saturation of recognition signals blocks the emergence of new tRNA identities and that the rate of nucleotide substitutions in tRNAs is higher in species with fewer tRNA genes. We propose that the growth of the genetic code stalled because a limit was reached in the number of identity elements that can be effectively used in the tRNA structure.

Insights into the mechanisms of RNA secondary structure destabilization by the HIV-1 nucleocapsid protein

The mature HIV-1 nucleocapsid protein NCp7 (NC) plays a key role in reverse transcription facilitating the two obligatory strand transfers. Several properties contribute to its efficient chaperon activity: preferential binding to single-stranded regions, nucleic acid aggregation, helix destabilization, and rapid dissociation from nucleic acids. However, little is known about the relationships between these different properties, which are complicated by the ability of the protein to recognize particular HIV-1 stem-loops, such as SL1, SL2, and SL3, with high affinity and without destabilizing them. These latter properties are important in the context of genome packaging, during which NC is part of the Gag precursor. We used NMR to investigate destabilization of the full-length TAR (trans activating response element) RNA by NC, which is involved in the first strand transfer step of reverse transcription. NC was used at a low protein:nucleotide (nt) ratio of 1:59 in these experiments. NMR data for the imino protons of TAR identified most of the base pairs destabilized by NC. These base pairs were adjacent to the loops in the upper part of the TAR hairpin rather than randomly distributed. Gel retardation assays showed that conversion from the initial TAR-cTAR complex to the fully annealed form occurred much more slowly at the 1:59 ratio than at the higher ratios classically used. Nevertheless, NC significantly accelerated the formation of the initial complex at a ratio of 1:59.

Regulatory effects of cotranscriptional RNA structure formation and transitions

RNAs, which play significant roles in many fundamental biological processes of life, fold into sophisticated and precise structures. RNA folding is a dynamic and intricate process, which conformation transition of coding and noncoding RNAs form the primary elements of genetic regulation. The cellular environment contains various intrinsic and extrinsic factors that potentially affect RNA folding in vivo, and experimental and theoretical evidence increasingly indicates that the highly flexible features of the RNA structure are affected by these factors, which include the flanking sequence context, physiochemical conditions, cis RNA-RNA interactions, and RNA interactions with other molecules. Furthermore, distinct RNA structures have been identified that govern almost all steps of biological processes in cells, including transcriptional activation and termination, transcriptional mutagenesis, 5'-capping, splicing, 3'-polyadenylation, mRNA export and localization, and translation. Here, we briefly summarize the dynamic and complex features of RNA folding along with a wide variety of intrinsic and extrinsic factors that affect RNA folding. We then provide several examples to elaborate RNA structure-mediated regulation at the transcriptional and posttranscriptional levels. Finally, we illustrate the regulatory roles of RNA structure and discuss advances pertaining to RNA structure in plants. WIREs RNA 2016, 7:562-574. doi: 10.1002/wrna.1350 For further resources related to this article, please visit the WIREs website.

Conformational Rearrangements of Individual Nucleotides during RNA-Ligand Binding Are Rate-Differentiated

A pronounced rate differentiation has been found for conformational rearrangements of individual nucleobases that occur during ligand recognition of the preQ1 class-I riboswitch aptamer from Thermoanaerobacter tengcongensis. Rate measurements rely on the 2ApFold approach by analyzing the fluorescence response of riboswitch variants, each with a single, strategically positioned 2-aminopurine nucleobase substitution. Observed rate discrimination between the fastest and the slowest conformational adaption is 22-fold, with the largest rate observed for the rearrangement of a nucleoside directly at the binding site and the smallest rate observed for the 3'-unpaired nucleoside that stacks onto the pseudo-knot-closing Watson-Crick base pair. Our findings provide novel insights into how compact, prefolded RNAs that follow the induced-fit recognition mechanism adapt local structural elements in response to ligand binding on a rather broad time scale and how this process culminates in a structural signal that is responsible for efficient downregulation of ribosomal translation.

Rapid Exchange Between Free and Bound States in RNA-Dendrimer Polyplexes: Implications on the Mechanism of Delivery and Release

A combination of solution NMR, dynamic light scattering (DLS), and fluorescence quenching assays were employed to obtain insights into the dynamics and structural features of a polyplex system consisting of HIV-1 transactivation response element (TAR) and PEGylated generation 5 poly(amidoamine) dendrimer (G5-PEG). NMR chemical shift mapping and (13)C spin relaxation based dynamics measurements depict the polyplex system as a highly dynamic assembly where the RNA, with its local structure and dynamics preserved, rapidly exchanges ( Collapse

Key Words

• therapeutic rna delivery
• pamam
• pegylation
• tar
• dls
• nmr spin relaxation

Collapse

MESH Headings

• Dendrimers/chemistry
• HIV-1/chemistry
• Magnetic Resonance Spectroscopy/methods
• Polyamines/chemistry
• Polyethylene Glycols/chemistry
• RNA/chemistry
• Response Elements/genetics
• Transfection/methods

Collapse

Grants

• R01 AI066975 NIAID NIH HHS
• P50 GM103297 NIGMS NIH HHS
• R01 AI087463-01 NIAID NIH HHS
• R01 EB005028 NIBIB NIH HHS
• R01EB005028 NIBIB NIH HHS
• R01 GM089846 NIGMS NIH HHS

Collapse

Affiliation(s)

Anisha Shakya Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Casey A. Dougherty Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Yi Xue Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Hashim M. Al-Hashimi Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Mark M. Banaszak Holl Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA

Collapse

Understanding the kinetic mechanism of RNA single base pair formation

RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.

RNA as a stable polymer to build controllable and defined nanostructures for material and biomedical applications

The value of polymers is manifested in their vital use as building blocks in material and life sciences. Ribonucleic acid (RNA) is a polynucleic acid, but its polymeric nature in materials and technological applications is often overlooked due to an impression that RNA is seemingly unstable. Recent findings that certain modifications can make RNA resistant to RNase degradation while retaining its authentic folding property and biological function, and the discovery of ultra-thermostable RNA motifs have adequately addressed the concerns of RNA unstability. RNA can serve as a unique polymeric material to build varieties of nanostructures including nanoparticles, polygons, arrays, bundles, membrane, and microsponges that have potential applications in biomedical and material sciences. Since 2005, more than a thousand publications on RNA nanostructures have been published in diverse fields, indicating a remarkable increase of interest in the emerging field of RNA nanotechnology. In this review, we aim to: delineate the physical and chemical properties of polymers that can be applied to RNA; introduce the unique properties of RNA as a polymer; review the current methods for the construction of RNA nanostructures; describe its applications in material, biomedical and computer sciences; and, discuss the challenges and future prospects in this field.

Single Molecule Cluster Analysis dissects splicing pathway conformational dynamics

We report Single Molecule Cluster Analysis (SiMCAn), which utilizes hierarchical clustering of hidden Markov modeling-fitted single-molecule fluorescence resonance energy transfer (smFRET) trajectories to dissect the complex conformational dynamics of biomolecular machines. We used this method to study the conformational dynamics of a precursor mRNA during the splicing cycle as carried out by the spliceosome. By clustering common dynamic behaviors derived from selectively blocked splicing reactions, SiMCAn was able to identify the signature conformations and dynamic behaviors of multiple ATP-dependent intermediates. In addition, it identified an open conformation adopted late in splicing by a 3' splice-site mutant, invoking a mechanism for substrate proofreading. SiMCAn enables rapid interpretation of complex single-molecule behaviors and should prove useful for the comprehensive analysis of a plethora of dynamic cellular machines.

Elastic network models for RNA: a comparative assessment with molecular dynamics and SHAPE experiments

Elastic network models (ENMs) are valuable and efficient tools for characterizing the collective internal dynamics of proteins based on the knowledge of their native structures. The increasing evidence that the biological functionality of RNAs is often linked to their innate internal motions poses the question of whether ENM approaches can be successfully extended to this class of biomolecules. This issue is tackled here by considering various families of elastic networks of increasing complexity applied to a representative set of RNAs. The fluctuations predicted by the alternative ENMs are stringently validated by comparison against extensive molecular dynamics simulations and SHAPE experiments. We find that simulations and experimental data are systematically best reproduced by either an all-atom or a three-beads-per-nucleotide representation (sugar-base-phosphate), with the latter arguably providing the best balance of accuracy and computational complexity.

High-throughput cellular RNA device engineering

Methods for rapidly assessing sequence-structure-function landscapes and developing conditional gene-regulatory devices are critical to our ability to manipulate and interface with biology. We describe a framework for engineering RNA devices from preexisting aptamers that exhibit ligand-responsive ribozyme tertiary interactions. Our methodology utilizes cell sorting, high-throughput sequencing, and statistical data analyses to enable parallel measurements of the activities of hundreds of thousands of sequences from RNA device libraries in the absence and presence of ligands. Our tertiary interaction RNA devices exhibit improved performance in terms of gene silencing, activation ratio, and ligand sensitivity as compared to optimized RNA devices that rely on secondary structure changes. We apply our method to building biosensors for diverse ligands and determine consensus sequences that enable ligand-responsive tertiary interactions. These methods advance our ability to develop broadly applicable genetic tools and to elucidate understanding of the underlying sequence-structure-function relationships that empower rational design of complex biomolecules.

Conjoined use of EM and NMR in RNA structure refinement

More than 40% of the RNA structures have been determined using nuclear magnetic resonance (NMR) technique. NMR mainly provides local structural information of protons and works most effectively on relatively small biomacromolecules. Hence structural characterization of large RNAs can be difficult for NMR alone. Electron microscopy (EM) provides global shape information of macromolecules at nanometer resolution, which should be complementary to NMR for RNA structure determination. Here we developed a new energy term in Xplor-NIH against the density map obtained by EM. We conjointly used NMR and map restraints for the structure refinement of three RNA systems — U2/U6 small-nuclear RNA, genome-packing motif (Ψ^CD)₂ from Moloney murine leukemia virus, and ribosome-binding element from turnip crinkle virus. In all three systems, we showed that the incorporation of a map restraint, either experimental or generated from known PDB structure, greatly improves structural precision and accuracy. Importantly, our method does not rely on an initial model assembled from RNA duplexes, and allows full torsional freedom for each nucleotide in the torsion angle simulated annealing refinement. As increasing number of macromolecules can be characterized by both NMR and EM, the marriage between the two techniques would enable better characterization of RNA three-dimensional structures.