RNA Structural Dynamics Research Articles

Monitoring Real-Time Temperature Dynamics of a Short RNA Hairpin Using Förster Resonance Energy Transfer and Circular Dichroism.

RNA molecules play crucial roles in gene expression regulation and cellular signaling, and these functions are governed by the formation of RNA secondary and tertiary structures. These structures are highly dynamic and subject to rapid changes in response to environmental cues, temperature in particular. Thermosensitive RNA secondary structures have been harnessed by multiple organisms to survey their temperature environment and to adjust gene expression accordingly. It is thus highly desirable to observe RNA structural changes in real time over a range of temperatures. Multiple approaches have been developed to study structural dynamics, but many of these require extensive processing of the RNA, large amounts of RNA input, and/or cannot be applied under physiological conditions. Here, we describe the use of a dually fluorescently labeled RNA oligonucleotide (containing a predicted hairpin structure) to monitor subtle RNA structural dynamics in vitro by Förster resonance energy transfer (FRET) and circular dichroism (CD) spectroscopy. These approaches can be employed under physiologically relevant conditions over a range of temperatures and with RNA concentrations as low as 200nM; they enable us to observe RNA structural dynamics in real time and to correlate these dynamics with changes in biological processes such as translation.

Structural Analysis of the Hepatitis B Virus RNA Encapsidation Signal ε by Selective 2'-Hydroxyl Acylation Analyzed by Primer Extension (SHAPE).

RNA structure is crucial for RNA function, including in viral cis-elements such as the hepatitis B virus (HBV) RNA encapsidation signal ε. Interacting with the viral polymerase ε mediates packaging of the pregenomic (pg) RNA into capsids, initiation of reverse transcription, and it affects the mRNA functions of pgRNA. As free RNA, the 61-nucleotide (nt) ε sequence adopts a bipartite stem-loop structure with a central bulge and an apical loop. Due to stable Watson-Crick base pairing, this was already predicted by early RNA folding programs and confirmed by classical enzymatic and chemical structure probing. A newer, high-resolution probing technique exploits the selective acylation of solvent-accessible 2'-hydroxyls in the RNA backbone by electrophilic compounds such as 2-methylnicotinic acid imidazolide (NAI), followed by mapping of the modified sites by primer extension. This SHAPE principle has meanwhile been extended to numerous applications. Here we provide a basic protocol for NAI-based SHAPE of isolated HBV ε RNA which already provided insights into the impact of mutations, and preliminarily, of polymerase binding on the RNA structural dynamics. While the focus is on NAI modification, we also briefly cover target RNA preparation by in vitro transcription, primer extension using a radiolabeled primer, and analysis of the resulting cDNAs by denaturing polyacrylamide gelelectrophoresis (PAGE). Given the high tolerance of SHAPE chemistry to different conditions, including applicability in live cells, we expect this technique to greatly facilitate deciphering the conformational dynamics underlying the various functions of the ε element, especially in concert with the recently solved three-dimensional structure of the free RNA.

Visualizing a two-state conformational ensemble in stem-loop 3 of the transcriptional regulator 7SK RNA.

Structural plasticity is integral to RNA function; however, there are currently few methods to quantitatively resolve RNAs that have multiple structural states. NMR spectroscopy is a powerful approach for resolving conformational ensembles but is size-limited. Chemical probing is well-suited for large RNAs but provides limited structural and kinetics information. Here, we integrate the two approaches to visualize a two-state conformational ensemble for the central stem-loop 3 (SL3) of 7SK RNA, a critical element for 7SK RNA function in transcription regulation. We find that the SL3 distal end exchanges between two equally populated yet structurally distinct states in both isolated SL3 constructs and full-length 7SK RNA. We rationally designed constructs that lock SL3 into a single state and demonstrate that both chemical probing and NMR data fit to a linear combination of the two states. Comparison of vertebrate 7SK RNA sequences shows either or both states are highly conserved. These results provide new insights into 7SK RNA structural dynamics and demonstrate the utility of integrating chemical probing with NMR spectroscopy to gain quantitative insights into RNA conformational ensembles.

Cellular roadmaps of viroid infection.

Viroids are single-stranded circular noncoding RNAs that infect plants. According to the International Committee on Taxonomy of Viruses, there are 44 viroids known to date. Notably, more than 20 000 distinct viroid-like RNA sequences have recently been identified in existing sequencing datasets, suggesting an unprecedented complexity in biological roles of viroids and viroid-like RNAs. Interestingly, a human pathogen, hepatitis delta virus (HDV), also replicates via a rolling circle mechanism like viroids. Therefore, knowledge of viroid infection is informative for research on HDV and other viroid-like RNAs reported from various organisms. Here, we summarize recent advancements in understanding viroid shuttling among subcellular compartments for completing replication cycles, emphasizing regulatory roles of RNA motifs and structural dynamics in diverse biological processes. We also compare the knowledge of viroid intracellular trafficking with known pathways governing cellular RNA movement in cells. Future investigations on regulatory RNA structures and cognate factors in regulating viroid subcellular trafficking and replication will likely provide new insights into RNA structure-function relationships and facilitate the development of strategies controlling RNA localization and function in cells.

Biocontrol potential of lipopeptides produced by Paenibacillus polymyxa AF01 against Neoscytalidium dimidiatum in pitaya.

Pitaya canker, caused by Neoscytalidium dimidiatum, is one of the most important fungal diseases that cause significant losses in production. To replace chemical pesticides, the use of biocontrol strains to manage plant diseases has been the focus of research. In this study, the bacterial strain AF01, identified as Paenibacillus polymyxa, exhibited significant antifungal effects against N. dimidiatum and four other pitaya fungal pathogens. The strain P. polymyxa AF01 produces 13 fusaricidins, which directly inhibit mycelial growth, spore germination and germ tube elongation by causing the membrane integrity and cell ultrastructure to incur irreversible damage. Pot experiment and yield test confirmed that AF01 provided preservative effects by reducing the disease index. In comparison to the untreated control groups, RNA-seq data showed that P. polymyxa AF01 selectively blocked some transcription and translation processes and inhibited RNA and DNA structural dynamics, energy production and conversion, and signal transduction, particularly cell wall biosynthesis, changes in membrane permeability, and impairment of protein biosynthesis. Thus, P. polymyxa AF01 could be potentially useful as a suitable biocontrol agent for pitaya canker.

The Evaluation of SHAPE-MaP RNA Structure Probing Protocols Reveals a Novel Role of Mn2+ in the Detection of 2'-OH Adducts.

Chemical probing, for decades, has been one of the most popular tools for studying the secondary structure of RNA molecules. Recently, protocols for simultaneous analysis of multiple RNAs have been developed, enabling in vivo transcriptome-wide interrogation of the RNA structure dynamics. One of the most popular methods is the selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP). In this study, we describe the evaluation of this protocol by addressing the influence of the reverse transcription enzymes, buffer conditions, and chemical probes on the properties of the cDNA library and the quality of mutational profiling-derived structural signals. Our results reveal a SuperScript IV (SSIV) reverse transcriptase as a more efficient enzyme for mutational profiling of SHAPE adducts and shed new light on the role of Mn2+ cations in the modulation of SSIV readthrough efficiency.

Identification of metal ion-binding sites in RNA structures using deep learning method.

Metal ion is an indispensable factor for the proper folding, structural stability and functioning of RNA molecules. However, it is very difficult for experimental methods to detect them in RNAs. With the increase of experimentally resolved RNA structures, it becomes possible to identify the metal ion-binding sites in RNA structures through in-silico methods. Here, we propose an approach called Metal3DRNA to identify the binding sites of the most common metal ions (Mg2+, Na+ and K+) in RNA structures by using a three-dimensional convolutional neural network model. The negative samples, screened out based on the analysis for binding surroundings of metal ions, are more like positive ones than the randomly selected ones, which are beneficial to a powerful predictor construction. The microenvironments of the spatial distributions of C, O, N and P atoms around a sample are extracted as features. Metal3DRNA shows a promising prediction power, generally surpassing the state-of-the-art methods FEATURE and MetalionRNA. Finally, utilizing the visualization method, we inspect the contributions of nucleotide atoms to the classification in several cases, which provides a visualization that helps to comprehend the model. The method will be helpful for RNA structure prediction and dynamics simulation study. Availability and implementation: The source code is available at https://github.com/ChunhuaLiLab/Metal3DRNA.

RNA structure mediated thermoregulation: What can we learn from plants?

RNA molecules have the capacity to form a multitude of distinct secondary and tertiary structures, but only the most energetically favorable conformations are adopted at any given time. Formation of such structures strongly depends on the environment and consequently, these structures are highly dynamic and may refold as their surroundings change. Temperature is one of the most direct physical parameters that influence RNA structure dynamics, and in turn, thermosensitive RNA structures can be harnessed by a cell to perceive and respond to its temperature environment. Indeed, many thermosensitive RNA structures with biological function have been identified in prokaryotic organisms, but for a long time such structures remained elusive in eukaryotes. Recent discoveries, however, reveal that thermosensitive RNA structures are also found in plants, where they affect RNA stability, pre-mRNA splicing and translation efficiency in a temperature-dependent manner. In this minireview, we provide a short overview of thermosensitive RNA structures in prokaryotes and eukaryotes, highlight recent advances made in identifying such structures in plants and discuss their similarities and differences to established prokaryotic RNA thermosensors.

Molecular Simulations Matching Denaturation Experiments for N6-Methyladenosine.

Post-transcriptional modifications are crucial for RNA function and can affect its structure and dynamics. Force-field-based classical molecular dynamics simulations are a fundamental tool to characterize biomolecular dynamics, and their application to RNA is flourishing. Here, we show that the set of force-field parameters for N6-methyladenosine (m6A) developed for the commonly used AMBER force field does not reproduce duplex denaturation experiments and, specifically, cannot be used to describe both paired and unpaired states. Then, we use reweighting techniques to derive new parameters matching available experimental data. The resulting force field can be used to properly describe paired and unpaired m6A in both syn and anti conformation, which thus opens the way to the use of molecular simulations to investigate the effects of N6 methylations on RNA structural dynamics.

Automatic Learning of Hydrogen-Bond Fixes in the AMBER RNA Force Field.

The capability of current force fields to reproduce RNA structural dynamics is limited. Several methods have been developed to take advantage of experimental data in order to enforce agreement with experiments. Here, we extend an existing framework which allows arbitrarily chosen force-field correction terms to be fitted by quantification of the discrepancy between observables back-calculated from simulation and corresponding experiments. We apply a robust regularization protocol to avoid overfitting and additionally introduce and compare a number of different regularization strategies, namely, L1, L2, Kish size, relative Kish size, and relative entropy penalties. The training set includes a GACC tetramer as well as more challenging systems, namely, gcGAGAgc and gcUUCGgc RNA tetraloops. Specific intramolecular hydrogen bonds in the AMBER RNA force field are corrected with automatically determined parameters that we call gHBfixopt. A validation involving a separate simulation of a system present in the training set (gcUUCGgc) and new systems not seen during training (CAAU and UUUU tetramers) displays improvements regarding the native population of the tetraloop as well as good agreement with NMR experiments for tetramers when using the new parameters. Then, we simulate folded RNAs (a kink–turn and L1 stalk rRNA) including hydrogen bond types not sufficiently present in the training set. This allows a final modification of the parameter set which is named gHBfix21 and is suggested to be applicable to a wider range of RNA systems.

Genome-wide analysis of the in vivo tRNA structurome reveals RNA structural and modification dynamics under heat stress

RNA structure plays roles in myriad cellular events including transcription, translation, and RNA processing. Genome-wide analyses of RNA secondary structure invivo by chemical probing have revealed critical structural features of mRNAs and long ncRNAs. Here, we examine the invivo secondary structure of a small RNA class, tRNAs. Study of tRNA structure is challenging because tRNAs are heavily modified and strongly structured. We introduce "tRNA structure-seq," a new workflow that accurately determines invivo secondary structures of tRNA. The workflow combines dimethyl sulfate (DMS) probing, ultra-processive RT, and mutational profiling (MaP), which provides mutations opposite DMS and natural modifications thereby allowing multiple modifications to be identified in a single read. We applied tRNA structure-seq to E. coli under control and stress conditions. A leading folding algorithm predicts E. coli tRNA structures with only ∼80% average accuracy from sequence alone. Strikingly, tRNA structure-seq, by providing experimental restraints, improves structure prediction under invivo conditions to ∼95% accuracy, with more than 14 tRNAs predicted completely correctly. tRNA structure-seq also quantifies the relative levels of tRNAs and their natural modifications at single nucleotide resolution, as validated by LC-MS/MS. Our application of tRNA structure-seq yields insights into tRNA structure in living cells, revealing that it is not immutable but has dynamics, with partial unfolding of secondary and tertiary tRNA structure under heat stress that is correlated with a loss of tRNA abundance. This method is applicable to other small RNAs, including those with natural modifications and highly structured regions.

3’ Terminal end‐labeling of structured RNA by Klenow DNA polymerase

RNA labeling is a key step in studying RNA localization, folding, and structural dynamics using fluorescence based assays. Current methods of fluorophore attachment to long (>100 nucleotides) and structured RNAs are hampered by low labeling yields and limitations in attachment position. These challenges have led to an interest in developing alternative methods that would enable the site‐specific labeling of long and structured RNAs. To this end, we have developed a site‐specific, chemo‐enzymatic method of fluorophore incorporation to RNAs. This novel method uses Klenow enzyme, a DNA polymerase, to add modified nucleotides to RNAs 3’ end in a templated manner. To facilitate this modified nucleotide addition in structured RNAs, we designed helper DNA oligonucleotides to disrupt any pre‐existing RNA intramolecular structure. Using this method, we have achieved efficient 3’ terminal end‐labeling of chemically synthesized and in‐vitro transcribed RNAs of different lengths with cyanine dyes.This method will be routinely used in fluorescence‐based assays that investigate the structure‐function relationships of RNAs involved in cancer and HIV progression.

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

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

Discovery of a large-scale, cell-state-responsive allosteric switch in the 7SK RNA using DANCE-MaP

7SK is a conserved noncoding RNA that regulates transcription by sequestering the transcription factor P-TEFb. 7SK function entails complex changes in RNA structure, but characterizing RNA dynamics in cells remains an unsolved challenge. We developed a single-molecule chemical probing strategy, DANCE-MaP (deconvolution and annotation of ribonucleic conformational ensembles), that defines per-nucleotide reactivity, direct base pairing interactions, tertiary interactions, and thermodynamic populations for each state in RNA structural ensembles from a single experiment. DANCE-MaP reveals that 7SK RNA encodes a large-scale structural switch that couples dissolution of the P-TEFb binding site to structural remodeling at distal release factor binding sites. The 7SK structural equilibrium shifts in response to cell growth and stress and can be targeted to modulate expression of P-TEFbresponsive genes. Our study reveals that RNA structural dynamics underlie 7SK function as an integrator of diverse cellular signals to control transcription and establishes the power of DANCE-MaP to define RNA dynamics in cells.

Visualizing RNA Structures by SAXS-Driven MD Simulations.

The biological role of biomolecules is intimately linked to their structural dynamics. Experimental or computational techniques alone are often insufficient to determine accurate structural ensembles in atomic detail. We use all-atom molecular dynamics (MD) simulations and couple it to small-angle X-ray scattering (SAXS) experiments to resolve the structural dynamics of RNA molecules. To accomplish this task, we utilize a set of re-weighting and biasing techniques tailored for RNA molecules. To showcase our approach, we study two RNA molecules: a riboswitch that shows structural variations upon ligand binding, and a two-way junction RNA that displays structural heterogeneity and sensitivity to salt conditions. Integration of MD simulations and experiments allows the accurate construction of conformational ensembles of RNA molecules. We observe a dynamic change of the SAM-I riboswitch conformations depending on its binding partners. The binding of SAM and Mg2+ cations stabilizes the compact state. The absence of Mg2+ or SAM leads to the loss of tertiary contacts, resulting in a dramatic expansion of the riboswitch conformations. The sensitivity of RNA structures to the ionic strength demonstrates itself in the helix junction helix (HJH). The HJH shows non-monotonic compaction as the ionic strength increases. The physics-based picture derived from the experimentally guided MD simulations allows biophysical characterization of RNA molecules. All in all, SAXS-guided MD simulations offer great prospects for studying RNA structural dynamics.

Genome-wide RNA structure changes during human neurogenesis modulate gene regulatory networks

The distribution, dynamics, and function of RNA structures in human development are under-explored. Here, we systematically assayed RNA structural dynamics and their relationship with gene expression, translation, and decay during human neurogenesis. We observed that the human ESC transcriptome is globally more structurally accessible than differentiated cells and undergoes extensive RNA structure changes, particularly in the 3' UTR. Additionally, RNA structure changes during differentiation are associated with translation and decay. We observed that RBP and miRNA binding is associated with RNA structural changes during early neuronal differentiation, and splicing is associated during later neuronal differentiation. Furthermore, our analysis suggests that RBPs are major factors in structure remodeling and co-regulate additional RBPs and miRNAs through structure. We demonstrated an example of this by showing that PUM2-induced structure changes on LIN28A enable miR-30 binding. This study deepens our understanding of the widespread and complex role of RNA-based gene regulation during human development.

The Role of RNA Secondary Structure in Regulation of Gene Expression in Bacteria.

Due to the high exposition to changing environmental conditions, bacteria have developed many mechanisms enabling immediate adjustments of gene expression. In many cases, the required speed and plasticity of the response are provided by RNA-dependent regulatory mechanisms. This is possible due to the very high dynamics and flexibility of an RNA structure, which provide the necessary sensitivity and specificity for efficient sensing and transduction of environmental signals. In this review, we will discuss the current knowledge about known bacterial regulatory mechanisms which rely on RNA structure. To better understand the structure-driven modulation of gene expression, we describe the basic theory on RNA structure folding and dynamics. Next, we present examples of multiple mechanisms employed by RNA regulators in the control of bacterial transcription and translation.

Monitoring Real-time Temperature Dynamics of a Short RNA Hairpin Using Förster Resonance Energy Transfer and Circular Dichroism

RNA secondary structures are highly dynamic and subject to prompt changes in response to the environment. Temperature in particular has a strong impact on RNA structural conformation, and temperature-sensitive RNA hairpin structures have been exploited by multiple organisms to modify the rate of translation in response to temperature changes. Observing RNA structural changes in real-time over a range of temperatures is therefore highly desirable. A variety of approaches exists that probe RNA secondary structures, but many of these either require large amount and/or extensive processing of the RNA or cannot be applied under physiological conditions, rendering the observation of structural dynamics over a range of temperatures difficult. Here, we describe the use of a dually fluorescently labelled RNA oligonucleotide (containing the predicted hairpin structure) that can be used to monitor subtle RNA-structural dynamics by Förster Resonance Energy Transfer (FRET) at different temperatures with RNA concentration as low as 200 nM. FRET efficiency varies as a function of the fluorophores' distance; high efficiency can thus be correlated to a stable hairpin structure, whilst a reduction in FRET efficiency reflects a partial opening of the hairpin or a destabilisation of this structure. The same RNA sequence can also be used for Circular Dichroism spectroscopy to observe global changes of RNA secondary structure at a given temperature. The combination of these approaches allowed us to monitor RNA structural dynamics over a range of temperatures in real-time and correlate structural changes to plant biology phenotypes. Graphic abstract:Monitoring temperature-dependent RNA structural dynamics using Förster Resonance Energy Transfer (FRET).

Unraveling the structure and biological functions of RNA triple helices.

It has been nearly 63 years since the first characterization of an RNA triple helix in vitro by Gary Felsenfeld, David Davies, and Alexander Rich. An RNA triple helix consists of three strands: A Watson–Crick RNA double helix whose major‐groove establishes hydrogen bonds with the so‐called “third strand”. In the past 15 years, it has been recognized that these major‐groove RNA triple helices, like single‐stranded and double‐stranded RNA, also mediate prominent biological roles inside cells. Thus far, these triple helices are known to mediate catalysis during telomere synthesis and RNA splicing, bind to ligands and ions so that metabolite‐sensing riboswitches can regulate gene expression, and provide a clever strategy to protect the 3′ end of RNA from degradation. Because RNA triple helices play important roles in biology, there is a renewed interest in better understanding the fundamental properties of RNA triple helices and developing methods for their high‐throughput discovery. This review provides an overview of the fundamental biochemical and structural properties of major‐groove RNA triple helices, summarizes the structure and function of naturally occurring RNA triple helices, and describes prospective strategies to isolate RNA triple helices as a means to establish the “triplexome”.This article is categorized under:RNA Structure and Dynamics > RNA Structure and DynamicsRNA Structure and Dynamics > RNA Structure, Dynamics and ChemistryRNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

RNA structural dynamics regulate early embryogenesis through controlling transcriptome fate and function

BackgroundVertebrate early embryogenesis is initially directed by a set of maternal RNAs and proteins, yet the mechanisms controlling this program remain largely unknown. Recent transcriptome-wide studies on RNA structure have revealed its pervasive and crucial roles in RNA processing and functions, but whether and how RNA structure regulates the fate of the maternal transcriptome have yet to be determined.ResultsHere we establish the global map of four nucleotide-based mRNA structures by icSHAPE during zebrafish early embryogenesis. Strikingly, we observe that RNA structurally variable regions are enriched in the 3′ UTR and contain cis-regulatory elements important for maternal-to-zygotic transition (MZT). We find that the RNA-binding protein Elavl1a stabilizes maternal mRNAs by binding to the cis-elements. Conversely, RNA structure formation suppresses Elavl1a’s binding leading to the decay of its maternal targets.ConclusionsOur study finds that RNA structurally variable regions are enriched in mRNA 3′ UTRs and contain cis-regulatory elements during zebrafish early embryogenesis. We reveal that Elavl1a regulates maternal RNA stability in an RNA structure-dependent fashion. Overall, our findings reveal a broad and fundamental role of RNA structure-based regulation in vertebrate early embryogenesis.