Dynamic changes in mRNA nucleocytoplasmic localization in the nitrate response of Arabidopsis roots.

Nitrate is a nutrient and signal that regulates gene expression. The nitrate response has been extensively characterized at the organism, organ, and cell-type-specific levels, but intracellular mRNA dynamics remain unexplored. To characterize nuclear and cytoplasmic transcriptome dynamics in response to nitrate, we performed a time-course expression analysis after nitrate treatment in isolated nuclei, cytoplasm, and whole roots. We identified 402 differentially localized transcripts (DLTs) in response to nitrate treatment. Induced DLT genes showed rapid and transient recruitment of the RNA polymerase II, together with an increase in the mRNA turnover rates. DLTs code for genes involved in metabolic processes, localization, and response to stimulus indicating DLTs include genes with relevant functions for the nitrate response that have not been previously identified. Using single-molecule RNA FISH, we observed early nuclear accumulation of the NITRATE REDUCTASE 1 (NIA1) transcripts in their transcription sites. We found that transcription of NIA1, a gene showing delayed cytoplasmic accumulation, is rapidly and transiently activated; however, its transcripts become unstable when they reach the cytoplasm. Our study reveals the dynamic localization of mRNAs between the nucleus and cytoplasm as an emerging feature in the temporal control of gene expression in response to nitrate treatment in Arabidopsis roots.

Genome-wide quantification of RNA flow across subcellular compartments reveals determinants of the mammalian transcript life cycle.

Dissecting the regulatory mechanisms controlling mammalian transcripts from production to degradation requires quantitative measurements of mRNA flow across the cell. We developed subcellular TimeLapse-seq to measure the rates at which RNAs are released from chromatin, exported from the nucleus, loaded onto polysomes, and degraded within the nucleus and cytoplasm in human and mouse cells. These rates varied substantially, yet transcripts from genes with related functions or targeted by the same transcription factors and RNA-binding proteins flowed across subcellular compartments with similar kinetics. Verifying these associations uncovered a link between DDX3X and nuclear export. For hundreds of RNA metabolism genes, most transcripts with retained introns were degraded by the nuclear exosome, while the remaining molecules were exported with stable cytoplasmic lifespans. Transcripts residing on chromatin for longer had extended poly(A) tails, whereas the reverse was observed for cytoplasmic mRNAs. Finally, machine learning identified molecular features that predicted the diverse life cycles of mRNAs.

Probing RNA structure and dynamics using nanopore and next generation sequencing.

It has become increasingly evident that the structures RNAs adopt are conformationally dynamic; the various structured states that RNAs sample govern their interactions with other nucleic acids, proteins, and ligands to regulate a myriad of biological processes. Although several biophysical approaches have been developed and used to study the dynamic landscape of structured RNAs, technical limitations have limited their application to all classes of RNA due to variable size and flexibility. Recent advances combining chemical probing experiments with next-generation- and direct sequencing have emerged as an alternative approach to exploring the conformational dynamics of RNA. In this review, we provide a methodological overview of the sequencing-based techniques used to study RNA conformational dynamics. We discuss how different techniques have enabled us to better understand the propensity of RNAs from a variety of different classes to sample multiple conformational states. Finally, we present examples of the ways these techniques have reshaped how we think about RNA structure.

Multiple transcriptome analyses reveal mouse testis developmental dynamics.

The testes are the organs of gamete production and testosterone synthesis. Up to date, no model system is available for mammalian testicular development, and only few studies have characterized the mouse testis transcriptome from no more than three postnatal ages. To describe the transcriptome landscape of the developing mouse testis and identify the potential molecular mechanisms underlying testis maturation, we examined multiple RNA-seq data of mouse testes from 3-week-old (puberty) to 11-week-old (adult). Sperm cells appeared as expected in 5-week-old mouse testis, suggesting the proper sample collection. The principal components analysis revealed the genes from 3w to 4w clustered away from other timepoints, indicating they may be the important nodes for testicular development. The pairwise comparisons at two adjacent timepoints identified 7,612 differentially expressed genes (DEGs), resulting in 58 unique mRNA expression patterns. Enrichment analysis identified functions in tissue morphogenesis (3-4w), regulation of peptidase activity (4-5w), spermatogenesis (7-8w), and antigen processing (10-11w), suggesting distinct functions in different developmental periods. 50 hub genes and 10 gene cluster modules were identified in the testis maturation process by protein-protein interaction (PPI) network analysis, and the miRNA-lncRNA-mRNA, miRNA-circRNA-mRNA and miRNA-circRNA-lncRNA-mRNA competing endogenous RNA (ceRNA) networks were constructed. The results suggest that testis maturation is a complex developmental process modulated by various molecules, and that some potential RNA-RNA interactions may be involved in specific developmental stages. In summary, this study provides an update on the molecular basis of testis development, which may help to understand the molecular mechanisms of mouse testis development and provide guidance for mouse reproduction.

RNA folding kinetics control riboswitch sensitivity in vivo.

Riboswitches are ligand-responsive gene-regulatory RNA elements that perform key roles in maintaining cellular homeostasis. Understanding how riboswitch sensitivity is controlled is critical to understanding how highly conserved aptamer domains are deployed in a variety of contexts with different sensitivity demands. Here we uncover new roles by which RNA folding dynamics control riboswitch sensitivity in cells. By investigating the Clostridium beijerinckii pfl ZTP riboswitch, we identify multiple mechanistic routes of altering expression platform sequence and structure to slow RNA folding, all of which enhance riboswitch sensitivity. Applying these methods to riboswitches with diverse aptamer architectures that regulate transcription and translation with ON and OFF logic demonstrates the generality of our findings, indicating that any riboswitch that operates in a kinetic regime can be sensitized by slowing expression platform folding. Comparison of the most sensitized versions of these switches to equilibrium aptamer:ligand dissociation constants suggests a limit to the sensitivities achievable by kinetic RNA switches. Our results add to the growing suite of knowledge and approaches that can be used to rationally program cotranscriptional RNA folding for biotechnology applications, and suggest general RNA folding principles for understanding dynamic RNA systems in other areas of biology.

Coordination of rhythmic RNA synthesis and degradation orchestrates 24- and 12-h RNA expression patterns in mouse fibroblasts.

Circadian RNA expression is essential to ultimately regulate a plethora of downstream rhythmic biochemical, physiological, and behavioral processes. Both transcriptional and posttranscriptional mechanisms are considered important to drive rhythmic RNA expression; however, the extent to which each regulatory process contributes to the rhythmic RNA expression remains controversial. To systematically address this, we monitored RNA dynamics using metabolic RNA labeling technology during a circadian cycle in mouse fibroblasts. We find that rhythmic RNA synthesis is the primary contributor of 24-h RNA rhythms, while rhythmic degradation is more important for 12-h RNA rhythms. These rhythms were predominantly regulated by Bmal1 and/or the core clock mechanism, and the interplay between rhythmic synthesis and degradation has a significant impact in shaping rhythmic RNA expression patterns. Interestingly, core clock RNAs are regulated by multiple rhythmic processes and have the highest amplitude of synthesis and degradation, presumably critical to sustain robust rhythmicity of cell-autonomous circadian rhythms. Our study yields invaluable insights into the temporal dynamics of both 24- and 12-h RNA rhythms in mouse fibroblasts.

Enhanced TROSY Effect in [2-19 F, 2-13 C] Adenosine and ATP Analogs Facilitates NMR Spectroscopy of Very Large Biological RNAs in Solution.

Large RNAs are central to cellular functions, but characterizing such RNAs remains challenging by solution NMR. We present two labeling technologies based on [2-19 F, 2-13 C]-adenosine, which allow the incorporation of aromatic 19 F-13 C spin pairs. The labels when coupled with the transverse relaxation optimized spectroscopy (TROSY) enable us to probe RNAs comprising up to 124 nucleotides. With our new [2-19 F, 2-13 C]-adenosine-phosphoramidite, all resonances of the human hepatitis B virus epsilon RNA could be readily assigned. With [2-19 F, 2-13 C]-adenosine triphosphate, the 124 nt pre-miR-17-NPSL1-RNA was produced via in vitro transcription and the TROSY spectrum of this 40 kDa [2-19 F, 2-13 C]-A-labeled RNA featured sharper resonances than the [2-1 H, 2-13 C]-A sample. The mutual cancelation of the chemical-shift-anisotropy and the dipole-dipole-components of TROSY-resonances leads to narrow linewidths over a wide range of molecular weights. With the synthesis of a non-hydrolysable [2-19 F, 2-13 C]-adenosine-triphosphate, we facilitate the probing of co-factor binding in kinase complexes and NMR-based inhibitor binding studies in such systems. Our labels allow a straightforward assignment for larger RNAs via a divide-and-conquer/mutational approach. The new [2-19 F, 2-13 C]-adenosine precursors are a valuable addition to the RNA NMR toolbox and will allow the study of large RNAs/RNA protein complexes in vitro and in cells.

Dynamic RNA synthetic biology: new principles, practices and potential.

An increased appreciation of the role of RNA dynamics in governing RNA function is ushering in a new wave of dynamic RNA synthetic biology. Here, we review recent advances in engineering dynamic RNA systems across the molecular, circuit and cellular scales for important societal-scale applications in environmental and human health, and bioproduction. For each scale, we introduce the core concepts of dynamic RNA folding and function at that scale, and then discuss technologies incorporating these concepts, covering new approaches to engineering riboswitches, ribozymes, RNA origami, RNA strand displacement circuits, biomaterials, biomolecular condensates, extracellular vesicles and synthetic cells. Considering the dynamic nature of RNA within the engineering design process promises to spark the next wave of innovation that will expand the scope and impact of RNA biotechnologies.

Allosteric activation of RhlB by RNase E induces partial duplex opening in substrate RNA.

The E. coli DEAD-Box helicase RhlB is responsible for ATP-dependent unwinding of structured mRNA to facilitate RNA degradation by the protein complex degradosome. The allosteric interaction with complex partner RNase E is necessary to stimulate both, RhlB's ATPase and RNA unwinding activity to levels comparable with other DEAD-Box helicases. However, the structural changes of the helicase RhlB induced by binding of RNase E have not been characterized and how those lead to increased reaction rates has remained unclear. We investigated the origin of this activation for RNA substrates with different topologies. Using NMR spectroscopy and an RNA centered approach, we could show that RNase E binding increases the affinity of RhlB towards a subset of RNA substrates, which leads to increased ATP turnover rates. Most strikingly, our studies revealed that in presence of RNase E (694-790) RhlB induces a conformational change in an RNA duplex with 5'- overhang even in absence of ATP, leading to partial duplex opening. Those results indicate a unique and novel activation mode of RhlB among DEAD-Box helicases, as ATP binding is thought to be an essential prerequisite for RNA unwinding.

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 no 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 conservation of both states, suggesting functional importance. 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.

Coordination of rhythmic RNA synthesis and degradation orchestrates 24-hour and 12-hour RNA expression patterns in mouse fibroblasts.

Circadian RNA expression is essential to ultimately regulate a plethora of downstream rhythmic biochemical, physiological, and behavioral processes. Both transcriptional and post-transcriptional mechanisms are considered important to drive rhythmic RNA expression, however, the extent to which each regulatory process contributes to the rhythmic RNA expression remains controversial. To systematically address this, we monitored RNA dynamics using metabolic RNA labeling technology during a circadian cycle in mouse fibroblasts. We find that rhythmic RNA synthesis is the primary contributor of 24 hr RNA rhythms, while rhythmic degradation is more important for 12 hr RNA rhythms. These rhythms were predominantly regulated by Bmal1 and/or the core clock mechanism, and interplay between rhythmic synthesis and degradation has a significant impact in shaping rhythmic RNA expression patterns. Interestingly, core clock RNAs are regulated by multiple rhythmic processes and have the highest amplitude of synthesis and degradation, presumably critical to sustain robust rhythmicity of cell-autonomous circadian rhythms. Our study yields invaluable insights into the temporal dynamics of both 24 hr and 12 hr RNA rhythms in mouse fibroblasts.

Live cell imaging of DNA and RNA with fluorescent signal amplification and background reduction techniques.

Illuminating DNA and RNA dynamics in live cell can elucidate their life cycle and related biochemical activities. Various protocols have been developed for labeling the regions of interest in DNA and RNA molecules with different types of fluorescent probes. For example, CRISPR-based techniques have been extensively used for imaging genomic loci. However, some DNA and RNA molecules can still be difficult to tag and observe dynamically, such as genomic loci in non-repetitive regions. In this review, we will discuss the toolbox of techniques and methodologies that have been developed for imaging DNA and RNA. We will also introduce optimized systems that provide enhanced signal intensity or low background fluorescence for those difficult-to-tag molecules. These strategies can provide new insights for researchers when designing and using techniques to visualize DNA or RNA molecules.

Multi-Site Conformational Exchange in the Synthetic Neomycin-Sensing Riboswitch Studied by 19 F NMR.

The synthetic neomycin-sensing riboswitch interacts with its cognate ligand neomycin as well as with the related antibiotics ribostamycin and paromomycin. Binding of these aminoglycosides induces a very similar ground state structure in the RNA, however, only neomycin can efficiently repress translation initiation. The molecular origin of these differences has been traced back to differences in the dynamics of the ligand:riboswitch complexes. Here, we combine five complementary fluorine based NMR methods to accurately quantify seconds to microseconds dynamics in the three riboswitch complexes. Our data reveal complex exchange processes with up to four structurally different states. We interpret our findings in a model that shows an interplay between different chemical groups in the antibiotics and specific bases in the riboswitch. More generally, our data underscore the potential of 19 F NMR methods to characterize complex exchange processes with multiple excited states.

a6A-seq: N 6-allyladenosine-based cellular messenger RNA metabolic labelling and sequencing.

The integration of RNA metabolic labelling by nucleoside analogues with high-throughput RNA sequencing has been harnessed to study RNA dynamics. The immunoprecipitation purification or chemical pulldown technique is generally required to enrich the analogue-labelled RNAs. Here we developed an a6A-seq method, which takes advantage of N6-allyladenosine (a6A) metabolic labelling on cellular mRNAs and profiles them in an immunoprecipitation-free and mutation-based manner. a6A plays a role as a chemical sequencing tag in that the iodination of a6A in mRNAs results in 1,N 6-cyclized adenosine (cyc-A), which induces base misincorporation during RNA reverse transcription, thus making a6A-labelled mRNAs detectable by sequencing. A nucleic acid melting assay was utilized to investigate why cyc-A prefers to be paired with guanine. a6A-seq was utilized to study cellular gene expression changes under a methionine-free stress condition. Compared with regular RNA-seq, a6A-seq could more sensitively detect the change of mRNA production over a time scale. The experiment of a6A-containing mRNA immunoprecipitation followed by qPCR successfully validated the high-throughput a6A-seq data. Together, our results show a6A-seq is an effective tool to study RNA dynamics.

An overview of structural approaches to study therapeutic RNAs.

RNAs provide considerable opportunities as therapeutic agent to expand the plethora of classical therapeutic targets, from extracellular and surface proteins to intracellular nucleic acids and its regulators, in a wide range of diseases. RNA versatility can be exploited to recognize cell types, perform cell therapy, and develop new vaccine classes. Therapeutic RNAs (aptamers, antisense nucleotides, siRNA, miRNA, mRNA and CRISPR-Cas9) can modulate or induce protein expression, inhibit molecular interactions, achieve genome editing as well as exon-skipping. A common RNA thread, which makes it very promising for therapeutic applications, is its structure, flexibility, and binding specificity. Moreover, RNA displays peculiar structural plasticity compared to proteins as well as to DNA. Here we summarize the recent advances and applications of therapeutic RNAs, and the experimental and computational methods to analyze their structure, by biophysical techniques (liquid-state NMR, scattering, reactivity, and computational simulations), with a focus on dynamic and flexibility aspects and to binding analysis. This will provide insights on the currently available RNA therapeutic applications and on the best techniques to evaluate its dynamics and reactivity.

Phosphomimetic S207D Lysyl-tRNA Synthetase Binds HIV-1 5'UTR in an Open Conformation and Increases RNA Dynamics.

Interactions between lysyl-tRNA synthetase (LysRS) and HIV-1 Gag facilitate selective packaging of the HIV-1 reverse transcription primer, tRNALys3. During HIV-1 infection, LysRS is phosphorylated at S207, released from a multi-aminoacyl-tRNA synthetase complex and packaged into progeny virions. LysRS is critical for proper targeting of tRNALys3 to the primer-binding site (PBS) by specifically binding a PBS-adjacent tRNA-like element (TLE), which promotes release of the tRNA proximal to the PBS. However, whether LysRS phosphorylation plays a role in this process remains unknown. Here, we used a combination of binding assays, RNA chemical probing, and small-angle X-ray scattering to show that both wild-type (WT) and a phosphomimetic S207D LysRS mutant bind similarly to the HIV-1 genomic RNA (gRNA) 5'UTR via direct interactions with the TLE and stem loop 1 (SL1) and have a modest preference for binding dimeric gRNA. Unlike WT, S207D LysRS bound in an open conformation and increased the dynamics of both the PBS region and SL1. A new working model is proposed wherein a dimeric phosphorylated LysRS/tRNA complex binds to a gRNA dimer to facilitate tRNA primer release and placement onto the PBS. Future anti-viral strategies that prevent this host factor-gRNA interaction are envisioned.

The dynamic, motile and deformative properties of RNA nanoparticles facilitate the third milestone of drug development.

Besides mRNA, rRNA, and tRNA, cells contain many other noncoding RNA that display critical roles in the regulation of cellular functions. Human genome sequencing revealed that the majority of non-protein-coding DNA actually codes for non-coding RNAs. The dynamic nature of RNA results in its motile and deformative behavior. These conformational transitions such as the change of base-pairing, breathing within complemented strands, and pseudoknot formation at the 2D level as well as the induced-fit and conformational capture at the 3D level are important for their biological functions including regulation, translation, and catalysis. The dynamic, motile and catalytic activity has led to a belief that RNA is the origin of life. We have recently reported that the deformative property of RNA nanoparticles enhances their penetration through the leaky blood vessel of cancers which leads to highly efficient tumor accumulation. This special deformative property also enables RNA nanoparticles to pass the glomerulus, overcoming the filtration size limit, resulting in fast renal excretion and rapid body clearance, thus low or no toxicity. The biodistribution of RNA nanoparticles can be further improved by the incorporation of ligands for cancer targeting. In addition to the favorable biodistribution profiles, RNA nanoparticles possess other properties including self-assembly, negative charge, programmability, and multivalency; making it a great material for pharmaceutical applications. The intrinsic negative charge of RNA nanoparticles decreases the toxicity of drugs by preventing nonspecific binding to the negative charged cell membrane and enhancing the solubility of hydrophobic drugs. The polyvalent property of RNA nanoparticles allows the multi-functionalization which can apply to overcome drug resistance. This review focuses on the summary of these unique properties of RNA nanoparticles, which describes the mechanism of RNA dynamic, motile and deformative properties, and elucidates and prepares to welcome the RNA therapeutics as the third milestone in pharmaceutical drug development.

Combining Coarse-Grained Simulations and Single Molecule Analysis Reveals a Three-State Folding Model of the Guanidine-II Riboswitch.

Riboswitch RNAs regulate gene expression by conformational changes induced by environmental conditions and specific ligand binding. The guanidine-II riboswitch is proposed to bind the small molecule guanidinium and to subsequently form a kissing loop interaction between the P1 and P2 hairpins. While an interaction was shown for isolated hairpins in crystallization and electron paramagnetic resonance experiments, an intrastrand kissing loop formation has not been demonstrated. Here, we report the first evidence of this interaction in cis in a ligand and Mg2+ dependent manner. Using single-molecule FRET spectroscopy and detailed structural information from coarse-grained simulations, we observe and characterize three interconvertible states representing an open and kissing loop conformation as well as a novel Mg2+ dependent state for the guanidine-II riboswitch from E. coli. The results further substantiate the proposed switching mechanism and provide detailed insight into the regulation mechanism for the guanidine-II riboswitch class. Combining single molecule experiments and coarse-grained simulations therefore provides a promising perspective in resolving the conformational changes induced by environmental conditions and to yield molecular insights into RNA regulation.

Estimating RNA dynamics using one time point for one sample in a single-pulse metabolic labeling experiment.

BACKGROUND: Over the past decade, experimental procedures such as metabolic labeling for determining RNA turnover rates at the transcriptome-wide scale have been widely adopted and are now turning to single cell measurements. Several computational methods to estimate RNA synthesis, processing and degradation rates from such experiments have been suggested, but they all require several RNA sequencing samples. Here we present a method that can estimate those three rates from a single sample. METHODS: Our method relies on the analytical solution to the Zeisel model of RNA dynamics. It was validated on metabolic labeling experiments performed on mouse embryonic stem cells. Resulting degradation rates were compared both to previously published rates on the same system and to a state-of-the-art method applied to the same data. RESULTS: Our method is computationally efficient and outputs rates that correlate well with previously published data sets. Using it on a single sample, we were able to reproduce the observation that dynamic biological processes tend to involve genes with higher metabolic rates, while stable processes involve genes with lower rates. This supports the hypothesis that cells control not only the mRNA steady-state abundance, but also its responsiveness, i.e., how fast steady state is reached. Moreover, degradation rates obtained with our method compare favourably with the other tested method. CONCLUSIONS: In addition to saving experimental work and computational time, estimating rates for a single sample has several advantages. It does not require an error-prone normalization across samples and enables the use of replicates to estimate uncertainty and assess sample quality. Finally the method and theoretical results described here are general enough to be useful in other contexts such as nucleotide conversion methods and single cell metabolic labeling experiments.

Thoughts on how to think (and talk) about RNA structure.

Recent events have pushed RNA research into the spotlight. Continued discoveries of RNA with unexpected diverse functions in healthy and diseased cells, such as the role of RNA as both the source and countermeasure to a severe acute respiratory syndrome coronavirus 2 infection, are igniting a new passion for understanding this functionally and structurally versatile molecule. Although RNA structure is key to function, many foundational characteristics of RNA structure are misunderstood, and the default state of RNA is often thought of and depicted as a single floppy strand. The purpose of this perspective is to help adjust mental models, equipping the community to better use the fundamental aspects of RNA structural information in new mechanistic models, enhance experimental design to test these models, and refine data interpretation. We discuss six core observations focused on the inherent nature of RNA structure and how to incorporate these characteristics to better understand RNA structure. We also offer some ideas for future efforts to make validated RNA structural information available and readily used by all researchers.