High-throughput determination of RNA structures

19F NMR in RNA structural biology: exploring structures, dynamics, and small molecule interactions

RNA molecules play essential roles in numerous biological pathways, making them attractive targets for drug discovery. Despite the challenges in developing small molecules targeting RNA, the success in developing compounds that modulate RNA function underscores its therapeutic potential. 19F NMR spectroscopy has emerged as a powerful tool in structural biology and drug discovery, particularly for studying macromolecular structures and ligand interactions. As RNA continues to gain prominence as a drug target, 19F NMR is expected to play a pivotal role in advancing RNA-focused drug discovery. This review describes the diverse applications of 19F NMR in RNA biology, including its use in characterizing RNA structures, probing molecular dynamics, identifying small-molecule binders, and investigating interaction mechanisms of small-molecule ligands. By providing detailed structural and ligand binding insights, 19F NMR will facilitate the discovery of RNA-targeting therapeutics and deepen our understanding of RNA modulatory mechanisms.

A comparative analysis of mRNA enrichment strategies and guidance for improving their efficiency

The pervasive presence of ribosomal RNA (rRNA) in total RNA poses a considerable challenge to the accurate analysis of cellular transcriptomes. In this study, we comprehensively analyze strategies for Saccharomyces cerevisiae mRNA enrichment, aiming to either separate polyadenylated RNAs or selectively remove rRNAs. Our findings reveal that a single round of mRNA enrichment in recommended conditions proves insufficient for both methods, prompting the exploration of strategies to enhance their efficiency. We show that adjusting the oligo (dT) magnetic beads-to-RNA ratio leads to significant improvement, but even better results can be achieved with two rounds of mRNA enrichment. We propose experimental conditions that reduce rRNA content in total yeast RNA to less than 10%, as confirmed by capillary electrophoresis and NG sequencing. Based on the obtained data, we selected the most time- and cost-effective option of polyadenylated RNA selection in yeast total RNA. Furthermore, we demonstrate that RNA modification with SHAPE reagent (NAI) does not interfere with the optimized mRNA enrichment protocol. These insights contribute to mRNA enrichment strategies and underscore the importance of optimizing mRNA isolation methodologies for downstream analyses.

Functions and therapeutic applications of pseudouridylation

The success of using pseudouridine (Ψ) and its methylation derivative in mRNA vaccines against SARS-CoV-2 has sparked a renewed interest in this RNA modification, known as the 'fifth nucleotide' of RNA. In this Review, we discuss the emerging functions of pseudouridylation in gene regulation, focusing on how pseudouridine in mRNA, tRNA and ribosomal RNA (rRNA) regulates translation. We also discuss the effects of pseudouridylation on RNA secondary structure, pre-mRNA splicing, and in vitro mRNA stability. In addition to nuclear-genome-encoded RNAs, pseudouridine is also present in mitochondria-encoded rRNA, mRNA and tRNA, where it has different distributions and functions compared with their nuclear counterparts. We then discuss the therapeutic potential of programmable pseudouridylation and mRNA vaccine optimization through pseudouridylation. Lastly, we briefly describe the latest quantitative pseudouridine detection methods. We posit that pseudouridine is a highly promising modification that merits further epitranscriptomics investigation and therapeutic application.

Structural determinants of inverted Alu-mediated backsplicing revealed by -MaP and -JuMP

Biogenesis of circular RNA usually involves a backsplicing reaction where the downstream donor site is ligated to the upstream acceptor site by the spliceosome. For this reaction to occur, these sites must be in proximity. Inverted repeat sequences, such as Alu elements, if positioned in the upstream and downstream introns, can base pair and represent one mechanism for inducing proximity. Here, we investigate the pre-mRNA structure of the human HIPK3 gene at exon 2, which forms a circular RNA via backsplicing. We leverage multiple chemical probing approaches, including the recently developed SHAPE-JuMP (selective 2'-hydroxyl acylation analyzed by primer extension and juxtaposed merged pairs) strategy, to characterize secondary and tertiary interactions in the pre-mRNA that govern backsplicing. Our data confirm that the antisense Alu elements AluSz(-) and AluSq2(+), in the upstream and downstream introns, form a highly paired interaction. Circularization requires formation of long-range Alu-mediated base pairs but does not require the full-length AluSq2(+). In addition to confirming long-range base pairs, our SHAPE-JuMP data identified multiple long-range interactions between non-pairing nucleotides. Genome-wide analysis of inverted repeats flanking circular RNAs confirms that the presence of these elements favors circularization, but with modest predictive power. Together, our study suggests that secondary structure considerations alone do not fully explain backsplicing and that additional interactions are involved.

Stress-induced changes in endogenous TP53 mRNA 5' regulatory region

Tumor suppressor protein p53 is regulated in a number of ways, including during initiation of TP53 mRNA translation. The 5' end of TP53 mRNA contains regulatory structures that enable noncanonical initiation using mechanisms that remain poorly described. Here we analyze per-nucleotide reactivity changes in the 5' end secondary structure of TP53 mRNA under in-cell conditions using A549 human lung carcinoma cells. We first construct a cell-free secondary structure model using SHAPE reagent 5-nitroisatoic anhydride on gently extracted and deproteinated RNA. We observe previously described regulatory features of the TP53 mRNA 5' end including two motifs which we refer to as long stem-loop (LSL) and short stem-loop (SSL), respectively. We observe a domain-forming helix that groups LSL and SSL, forming a three-helix junction. Applying in-cell selective 2' hydroxyl acylation analyzed by primer extension and mutational profiling, we assess reactivity profiles with unstressed cells and with chemically induced stress conditions expected to stimulate TP53 cap-independent translation. We analyze the effects of etoposide-induced DNA damage, CoCl2-induced hypoxia, and 5' cap inhibition with 4EGI-1 treatment. Identifying stress-associated changes in the TP53 5' end may help elucidate the role of regulatory RNA structure in cap-independent translation. Using ΔSHAPE, we identify in-cell protection sites that correspond with previously described RNA-protein binding sites on the apical loops of LSL and SSL. Furthermore, we identify several other potential interaction sites, some associated with specific types of stress. Some noteworthy changes include ΔSHAPE sites proximal to the start codons, at the three-helix junction and on the domain-forming helix. We summarize potential interactions on the cell-free secondary structure model.

Long-read RNA sequencing: A transformative technology for exploring transcriptome complexity in human diseases

Long-read RNA sequencing (RNA-seq) is emerging as a powerful and versatile technology for studying human transcriptomes. By enabling the end-to-end sequencing of full-length transcripts, long-read RNA-seq opens up avenues for investigating various RNA species and features that cannot be reliably interrogated by standard short-read RNA-seq methods. In this review, we present an overview of long-read RNA-seq, delineating its strengths over short-read RNA-seq, as well as summarizing recent advances in experimental and computational approaches to boost the power of long-read-based transcriptomics. We describe a wide range of applications of long-read RNA-seq, and highlight its expanding role as a foundational technology for exploring transcriptome variations in human diseases.

Comprehensive benchmarking of large language models for RNA secondary structure prediction

In recent years, inspired by the success of large language models (LLMs) for DNA and proteins, several LLMs for RNA have also been developed. These models take massive RNA datasets as inputs and learn, in a self-supervised way, how to represent each RNA base with a semantically rich numerical vector. This is done under the hypothesis that obtaining high-quality RNA representations can enhance data-costly downstream tasks, such as the fundamental RNA secondary structure prediction problem. However, existing RNA-LLM have not been evaluated for this task in a unified experimental setup. Since they are pretrained models, assessment of their generalization capabilities on new structures is a crucial aspect. Nonetheless, this has been just partially addressed in literature. In this work we present a comprehensive experimental and comparative analysis of pretrained RNA-LLM that have been recently proposed. We evaluate the use of these representations for the secondary structure prediction task with a common deep learning architecture. The RNA-LLM were assessed with increasing generalization difficulty on benchmark datasets. Results showed that two LLMs clearly outperform the other models, and revealed significant challenges for generalization in low-homology scenarios. Moreover, in this study we provide curated benchmark datasets of increasing complexity and a unified experimental setup for this scientific endeavor. Source code and curated benchmark datasets with increasing complexity are available in the repository: https://github.com/sinc-lab/rna-llm-folding/.

Consistent features observed in structural probing data of eukaryotic RNAs

Understanding RNA structure is crucial for elucidating its regulatory mechanisms. With the recent commercialization of messenger RNA vaccines, the profound impact of RNA structure on stability and translation efficiency has become increasingly evident, underscoring the importance of understanding RNA structure. Chemical probing of RNA has emerged as a powerful technique for investigating RNA structure in living cells. This approach utilizes chemical probes that selectively react with accessible regions of RNA, and by measuring reactivity, the openness and potential of RNA for protein binding or base pairing can be inferred. Extensive experimental data generated using RNA chemical probing have significantly contributed to our understanding of RNA structure in cells. However, it is crucial to acknowledge potential biases in chemical probing data to ensure an accurate interpretation. In this study, we comprehensively analyzed transcriptome-scale RNA chemical probing data in eukaryotes and report common features. Notably, in all experiments, the number of bases modified in probing was small, the bases showing the top 10% reactivity well reflected the known secondary structure, bases with high reactivity were more likely to be exposed to solvent and low reactivity did not reflect solvent exposure, which is important information for the analysis of RNA chemical probing data.

Characterizing 3D RNA structural features from DMS reactivity

Dimethyl sulfate (DMS) chemical mapping probes RNA structure, where low reactivity is generally interpreted as Watson-Crick (WC) base pairs and high reactivity as unpaired nucleotides. Studies examining DMS reactivity of RNAs with known 3D structures have identified nucleotides that deviate from this interpretation with distinct solvent accessibility and hydrogen bonding patterns. Understanding the frequency of these outliers and their recurring structural 3D features remains incomplete. To address this knowledge gap, we systematically analyzed DMS reactivity patterns across a library of 7,500 RNA constructs containing two-way junctions with known 3D structures. We observe DMS reactivity exists on a continuum over four orders of magnitude with approximately 10% overlap in reactivity between WC and non-WC nucleotides. We find that non-WC bases with WC-like DMS protection exhibit increased hydrogen bonding and decreased solvent accessibility, whereas WC pairs exhibiting greater DMS reactivity tend to flank junctions, correlating with weaker base stacking and greater junction dynamics. Furthermore, we discover that DMS reactivity values in non-canonical pairs correlate with atomic distances and base pair geometry, enabling discrimination between different 3D conformations. These DMS reactivity patterns indicate that DMS reactivity provides atomic-scale information about RNA 3D conformations, which can be used to model RNA structures and dynamics.

RASP v2.0: an updated atlas for RNA structure probing data

RNA molecules function in numerous biological processes by folding into intricate structures. Here we present RASP v2.0, an updated database for RNA structure probing data featuring a substantially expanded collection of datasets along with enhanced online structural analysis functionalities. Compared to the previous version, RASP v2.0 includes the following improvements: (i) the number of RNA structure datasets has increased from 156 to 438, comprising 216 transcriptome-wide RNA structure datasets, 141 target-specific RNA structure datasets, and 81 RNA-RNA interaction datasets, thereby broadening species coverage from 18 to 24, (ii) a deep learning-based model has been implemented to impute missing structural signals for 59 transcriptome-wide RNA structure datasets with low structure score coverage, significantly enhancing data quality, particularly for low-abundance RNAs, (iii) three new online analysis modules have been deployed to assist RNA structure studies, including missing structure score imputation, RNA secondary and tertiary structure prediction, and RNA binding protein (RBP) binding prediction. By providing a resource of much more comprehensive RNA structure data, RASP v2.0 is poised to facilitate the exploration of RNA structure-function relationships across diverse biological processes. RASP v2.0 is freely accessible at http://rasp2.zhanglab.net/.

The conformational space of RNase P RNA in solution

RNA conformational diversity has fundamental biological roles1-5, but direct visualization of its full conformational space in solution has not been possible using traditional biophysical techniques. Using solution atomic force microscopy, a deep neural network and statistical analyses, we show that the ribonuclease P (RNase P) RNA adopts heterogeneous conformations consisting of a conformationally invariant core and highly flexible peripheral structural elements that sample a broad conformational space, with amplitudes as large as 20-60 Å in a multitude of directions, with very low net energy cost. Increasing Mg2+ drives compaction and enhances enzymatic activity, probably by narrowing the conformational space. Moreover, analyses of the correlations and anticorrelations between spatial flexibility and sequence conservation suggest that the functional roles of both the structure and dynamics of key regions are embedded in the primary sequence. These findings reveal the structure-dynamics basis for the embodiment of both enzymatic precision and substrate promiscuity in the RNA component of the RNase P. Mapping the conformational space of the RNase P RNA demonstrates a new general approach to studying RNA structure and dynamics.

Structural determinants of inverted Alu-mediated backsplicing revealed by -MaP and -JuMP

Biogenesis of circular RNA usually involves a backsplicing reaction where the downstream donor site is ligated to the upstream acceptor site by the spliceosome. For this reaction to occur, it is hypothesized that these sites must be in proximity. Inverted repeat sequences, such as Alu elements, in the upstream and downstream introns are predicted to base-pair and represent one mechanism for inducing proximity. Here, we investigate the pre-mRNA structure of the human HIPK3 gene at exon 2, which forms a circular RNA via backsplicing. We leverage multiple chemical probing techniques, including the recently developed SHAPE- JuMP strategy, to characterize secondary and tertiary interactions in the pre- mRNA that govern backsplicing. Our data confirm that the antisense Alu elements, AluSz(-) and AluSq2(+) in the upstream and downstream introns, form a highly- paired interaction. Circularization requires formation of long-range Alu-mediated base pairs but does not require the full-length AluSq2(+). In addition to confirming long-range base pairs, our SHAPE-JuMP data identified multiple long-range interactions between non-pairing nucleotides. Genome-wide analysis of inverted repeats flanking circular RNAs confirm that their presence favors circularization, but the overall effect is modest. Together these results suggest that secondary structure considerations alone cannot fully explain backsplicing and additional interactions are key.

Cryo-EM Structure of raiA ncRNA From Clostridium Reveals a New RNA 3D Fold

Advancements in genome-wide sequence analysis have led to the discovery of numerous novel bacterial non-coding RNAs (ncRNAs). These ncRNAs have been categorized into various RNA families and classes based on their size, structure, function, and evolutionary relationships. One such ncRNA family, raiA, is notably abundant in the bacterial phyla Firmicutes and Actinobacteria and is remarkably well-conserved across many Gram-positive bacteria. In this study, we integrated cryo-electron microscopy single-particle analysis with computational modeling and biochemical techniques to elucidate the structural characteristics of raiA from Clostridium sp. CAG 138. Our findings reveal the globular 3D fold of raiA, providing valuable structural insights. This analysis paves the way for future investigations into the functional properties of raiA, potentially uncovering new regulatory mechanisms in bacterial ncRNAs.

Targeting RNA with small molecules, from RNA structures to precision medicines: IUPHAR review: 40

RNA plays important roles in regulating both health and disease biology in all kingdoms of life. Notably, RNA can form intricate three-dimensional structures, and their biological functions are dependent on these structures. Targeting the structured regions of RNA with small molecules has gained increasing attention over the past decade, because it provides both chemical probes to study fundamental biology processes and lead medicines for diseases with unmet medical needs. Recent advances in RNA structure prediction and determination and RNA biology have accelerated the rational design and development of RNA-targeted small molecules to modulate disease pathology. However, challenges remain in advancing RNA-targeted small molecules towards clinical applications. This review summarizes strategies to study RNA structures, to identify small molecules recognizing these structures, and to augment the functionality of RNA-binding small molecules. We focus on recent advances in developing RNA-targeted small molecules as potential therapeutics in a variety of diseases, encompassing different modes of actions and targeting strategies. Furthermore, we present the current gaps between early-stage discovery of RNA-binding small molecules and their clinical applications, as well as a roadmap to overcome these challenges in the near future.

Beyond ligand binding: Single molecule observation reveals how riboswitches integrate multiple signals to balance bacterial gene regulation

Riboswitches are specialized RNA structures that orchestrate gene expression in response to sensing specific metabolite or ion ligands, mostly in bacteria. Upon ligand binding, these conformationally dynamic RNA motifs undergo structural changes that control critical gene expression processes such as transcription termination and translation initiation, thereby enabling cellular homeostasis and adaptation. Because RNA folds rapidly and co-transcriptionally, riboswitches make use of the low complexity of RNA sequences to adopt alternative, transient conformations on the heels of the transcribing RNA polymerase (RNAP), resulting in kinetic partitioning that defines the regulatory outcome. This review summarizes single molecule microscopy evidence that has begun to unveil a sophisticated network of dynamic, kinetically balanced interactions between riboswitch architecture and the gene expression machinery that, together, integrate diverse cellular signals.

Identification of RNA structures and their roles in RNA functions

The development of high-throughput RNA structure profiling methods in the past decade has greatly facilitated our ability to map and characterize different aspects of RNA structures transcriptome-wide in cell populations, single cells and single molecules. The resulting high-resolution data have provided insights into the static and dynamic nature of RNA structures, revealing their complexity as they perform their respective functions in the cell. In this Review, we discuss recent technical advances in the determination of RNA structures, and the roles of RNA structures in RNA biogenesis and functions, including in transcription, processing, translation, degradation, localization and RNA structure-dependent condensates. We also discuss the current understanding of how RNA structures could guide drug design for treating genetic diseases and battling pathogenic viruses, and highlight existing challenges and future directions in RNA structure research.

Structural and biophysical dissection of RNA conformational ensembles

RNA's ability to form and interconvert between multiple secondary and tertiary structures is critical to its functional versatility and the traditional view of RNA structures as static entities has shifted towards understanding them as dynamic conformational ensembles. In this review we discuss RNA structural ensembles and their dynamics, highlighting the concept of conformational energy landscapes as a unifying framework for understanding RNA processes such as folding, misfolding, conformational changes, and complex formation. Ongoing advancements in cryo-electron microscopy and chemical probing techniques are significantly enhancing our ability to investigate multiple structures adopted by conformationally dynamic RNAs, while traditional methods such as nuclear magnetic resonance spectroscopy continue to play a crucial role in providing high-resolution, quantitative spatial and temporal information. We discuss how these methods, when used synergistically, can provide a comprehensive understanding of RNA conformational ensembles, offering new insights into their regulatory functions.

Intronic RNA secondary structural information captured for the human MYC pre-mRNA

To address the lack of intronic reads in secondary structure probing data for the human MYC pre-mRNA, we developed a method that combines spliceosomal inhibition with RNA probing and sequencing. Here, the SIRP-seq method was applied to study the secondary structure of human MYC RNAs by chemically probing HeLa cells with dimethyl sulfate in the presence of the small molecule spliceosome inhibitor pladienolide B. Pladienolide B binds to the SF3B complex of the spliceosome to inhibit intron removal during splicing, resulting in retained intronic sequences. This method was used to increase the read coverage over intronic regions of MYC. The purpose for increasing coverage across introns was to generate complete reactivity profiles for intronic sequences via the DMS-MaPseq approach. Notably, depth was sufficient for analysis by the program DRACO, which was able to deduce distinct reactivity profiles and predict multiple secondary structural conformations as well as their suggested stoichiometric abundances. The results presented here provide a new method for intronic RNA secondary structural analyses, as well as specific structural insights relevant to MYC RNA splicing regulation and therapeutic targeting.

Hfq mediates transcriptome-wide RNA structurome reprogramming under virulence-inducing conditions in a phytopathogen

Although RNA structures play important roles in regulating gene expression, the mechanism and function of mRNA folding in plant bacterial pathogens remain elusive. Therefore, we perform dimethyl sulfate sequencing (DMS-seq) on the Pseudomonas syringae under nutrition-rich and -deficient conditions, revealing that the mRNA structure changes substantially in the minimal medium (MM) that tunes global translation efficiency (TE), thereby inducing virulence. This process is led by the increased expression of hfq, which is directly activated by transcription regulators RpoS and CysB. The co-occurrence of Hfq and RpoS in diverse bacteria and the deep conservation of Hfq Y25 is critical for RNA-mediated regulation and implicates the wider biological importance of mRNA structure and feedback loops in the control of global gene expression.

SHAPEwarp-web: sequence-agnostic search for structurally homologous RNA regions across databases of chemical probing data

RNA molecules perform a variety of functions in cells, many of which rely on their secondary and tertiary structures. Chemical probing methods coupled with high-throughput sequencing have significantly accelerated the mapping of RNA structures, and increasingly large datasets of transcriptome-wide RNA chemical probing data are becoming available. Analogously to what has been done for decades in the protein world, this RNA structural information can be leveraged to aid the discovery of structural similarity to a known RNA (or RNA family), which, in turn, can inform about the function of transcripts. We have previously developed SHAPEwarp, a sequence-agnostic method for the search of structurally homologous RNA segments in a database of reactivity profiles derived from chemical probing experiments. In its original implementation, however, SHAPEwarp required substantial computational resources, even for moderately sized databases, as well as significant Linux command line know-how. To address these limitations, we introduce here SHAPEwarp-web, a user-friendly web interface to rapidly query large databases of RNA chemical probing data for structurally similar RNAs. Aside from featuring a completely rewritten core, which speeds up by orders of magnitude the search inside large databases, the web server hosts several high-quality chemical probing databases across multiple species. SHAPEwarp-web is available from https://shapewarp.incarnatolab.com.

The RNA secondary structure of androgen receptor-FL and V7 transcripts reveals novel regulatory regions

The androgen receptor (AR) is a ligand-dependent nuclear transcription factor belonging to the steroid hormone nuclear receptor family. Due to its roles in regulating cell proliferation and differentiation, AR is tightly regulated to maintain proper levels of itself and the many genes it controls. AR dysregulation is a driver of many human diseases including prostate cancer. Though this dysregulation often occurs at the RNA level, there are many unknowns surrounding post-transcriptional regulation of AR mRNA, particularly the role that RNA secondary structure plays. Thus, a comprehensive analysis of AR transcript secondary structure is needed. We address this through the computational and experimental analyses of two key isoforms, full length (AR-FL) and truncated (AR-V7). Here, a combination of in-cell RNA secondary structure probing experiments (targeted DMS-MaPseq) and computational predictions were used to characterize the static structural landscape and conformational dynamics of both isoforms. Additionally, in-cell assays were used to identify functionally relevant structures in the 5' and 3' UTRs of AR-FL. A notable example is a conserved stem loop structure in the 5'UTR of AR-FL that can bind to Poly(RC) Binding Protein 2 (PCBP2). Taken together, our results reveal novel features that regulate AR expression.

Ribonanza: deep learning of RNA structure through dual crowdsourcing

Prediction of RNA structure from sequence remains an unsolved problem, and progress has been slowed by a paucity of experimental data. Here, we present Ribonanza, a dataset of chemical mapping measurements on two million diverse RNA sequences collected through Eterna and other crowdsourced initiatives. Ribonanza measurements enabled solicitation, training, and prospective evaluation of diverse deep neural networks through a Kaggle challenge, followed by distillation into a single, self-contained model called RibonanzaNet. When fine tuned on auxiliary datasets, RibonanzaNet achieves state-of-the-art performance in modeling experimental sequence dropout, RNA hydrolytic degradation, and RNA secondary structure, with implications for modeling RNA tertiary structure.

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.

Phylogenetic and Chemical Probing Information as Soft Constraints in RNA Secondary Structure Prediction

Extrinsic, experimental information can be incorporated into thermodynamics-based RNA folding algorithms in the form of pseudo-energies. Evolutionary conservation of RNA secondary structure elements is detectable in alignments of phylogenetically related sequences and provides evidence for the presence of certain base pairs that can also be converted into pseudo-energy contributions. We show that the centroid base pairs computed from a consensus folding model such as RNAalifold result in a substantial improvement of the prediction accuracy for single sequences. Evidence for specific base pairs turns out to be more informative than a position-wise profile for the conservation of the pairing status. A comparison with chemical probing data, furthermore, strongly suggests that phylogenetic base pairing data are more informative than position-specific data on (un)pairedness as obtained from chemical probing experiments. In this context we demonstrate, in addition, that the conversion of signal from probing data into pseudo-energies is possible using thermodynamic structure predictions as a reference instead of known RNA structures.

Heterobifunctional small molecules to modulate RNA function

RNA has diverse cellular functionality, including regulating gene expression, protein translation, and cellular response to stimuli, due to its intricate structures. Over the past decade, small molecules have been discovered that target functional structures within cellular RNAs and modulate their function. Simple binding, however, is often insufficient, resulting in low or even no biological activity. To overcome this challenge, heterobifunctional compounds have been developed that can covalently bind to the RNA target, alter RNA sequence, or induce its cleavage. Herein, we review the recent progress in the field of RNA-targeted heterobifunctional compounds using representative case studies. We identify critical gaps and limitations and propose a strategic pathway for future developments of RNA-targeted molecules with augmented functionalities.

RNA: The Unsuspected Conductor in the Orchestra of Macromolecular Crowding

This comprehensive Review delves into the chemical principles governing RNA-mediated crowding events, commonly referred to as granules or biological condensates. We explore the pivotal role played by RNA sequence, structure, and chemical modifications in these processes, uncovering their correlation with crowding phenomena under physiological conditions. Additionally, we investigate instances where crowding deviates from its intended function, leading to pathological consequences. By deepening our understanding of the delicate balance that governs molecular crowding driven by RNA and its implications for cellular homeostasis, we aim to shed light on this intriguing area of research. Our exploration extends to the methodologies employed to decipher the composition and structural intricacies of RNA granules, offering a comprehensive overview of the techniques used to characterize them, including relevant computational approaches. Through two detailed examples highlighting the significance of noncoding RNAs, NEAT1 and XIST, in the formation of phase-separated assemblies and their influence on the cellular landscape, we emphasize their crucial role in cellular organization and function. By elucidating the chemical underpinnings of RNA-mediated molecular crowding, investigating the role of modifications, structures, and composition of RNA granules, and exploring both physiological and aberrant phase separation phenomena, this Review provides a multifaceted understanding of the intriguing world of RNA-mediated biological condensates.

Effects of mRNA conformational switching on translational noise in gene circuits

Intragenic translational heterogeneity describes the variation in translation at the level of transcripts for an individual gene. A factor that contributes to this source of variation is the mRNA structure. Both the composition of the thermodynamic ensemble, i.e., the stationary distribution of mRNA structures, and the switching dynamics between those play a role. The effect of the switching dynamics on intragenic translational heterogeneity remains poorly understood. We present a stochastic translation model that accounts for mRNA structure switching and is derived from a Markov model via approximate stochastic filtering. We assess the approximation on various timescales and provide a method to quantify how mRNA structure dynamics contributes to translational heterogeneity. With our approach, we allow quantitative information on mRNA switching from biophysical experiments or coarse-grain molecular dynamics simulations of mRNA structures to be included in gene regulatory chemical reaction network models without an increase in the number of species. Thereby, our model bridges a gap between mRNA structure kinetics and gene expression models, which we hope will further improve our understanding of gene regulatory networks and facilitate genetic circuit design.

RNAvigate: efficient exploration of RNA chemical probing datasets

Chemical probing technologies enable high-throughput examination of diverse structural features of RNA, including local nucleotide flexibility, RNA secondary structure, protein and ligand binding, through-space interaction networks, and multistate structural ensembles. Deep understanding of RNA structure-function relationships typically requires evaluating a system under structure- and function-altering conditions, linking these data with additional information, and visualizing multilayered relationships. Current platforms lack the broad accessibility, flexibility and efficiency needed to iterate on integrative analyses of these diverse, complex data. Here, we share the RNA visualization and graphical analysis toolset RNAvigate, a straightforward and flexible Python library that automatically parses 21 standard file formats (primary sequence annotations, per- and internucleotide data, and secondary and tertiary structures) and outputs 18 plot types. RNAvigate enables efficient exploration of nuanced relationships between multiple layers of RNA structure information and across multiple experimental conditions. Compatibility with Jupyter notebooks enables nonburdensome, reproducible, transparent and organized sharing of multistep analyses and data visualization strategies. RNAvigate simplifies and accelerates discovery and characterization of RNA-centric functions in biology.

Integrated analysis of lactate-related genes identifies POLRMT as a novel marker promoting the proliferation, migration and energy metabolism of hepatocellular carcinoma via Wnt/β-Catenin signaling

Hepatocellular carcinoma (HCC) is a prevalent and deadly form of cancer globally with typically unfavorable outcomes. Increasing research suggests that lactate serves as an important carbon contributor to cellular metabolism and holds a crucial part in the progression, sustenance, and treatment response of tumors. However, the contribution of lactate-related genes (LRGs) in HCC is still unclear. In this study, we analyzed TCGA datasets and screened 21 differentially expressed LRGs related to long-term survivals in HCC patients. Pan-cancer assays revealed that 21 LRGs expression exhibited a dysregulated level in man types of tumors and associated with clinical prognosis of tumor patients. The analysis of 21 LRGs successfully classified HCC samples into two molecular subtypes, and these two subtypes showed significant differences in clinical information, gene expression, and immune characteristics. Subsequently, based on the aforementioned 21 LRGs, a novel prognostic signature (DTYMK, IRAK1, POLRMT, MPV17, UQCRH, PDSS1, SLC16A3, SPP1 and LDHD) was generated by LASSO-Cox regression analysis. Survival assays demonstrated that the signature performed well in predicting the overall survival of patients with HCC. The results of Gene Set Variation Analysis indicated that the high GSVA scores were associated with poor prognosis. Moreover, we also investigated the correlation between GSVA scores and various signaling pathways in HCC. Among the nine prognostic genes, our attention focused on POLRMT which was highly expressed in HCC specimens based on TCGA datasets and several HCC cell lines. In addition, functional assays indicated that POLRMT distinctly promoted the proliferation, migration and energy metabolism of HCC cells via regulating Wnt/β-Catenin signaling. Overall, through the establishment of a novel prognostic signature, we have provided potential clinical value for assessing the prognosis of HCC patients. Furthermore, our study has identified the high expression of POLRMT in HCC and demonstrated its crucial role in HCC cell proliferation. These findings hold great importance in advancing our understanding of the pathophysiology of HCC, identifying new therapeutic targets, and improving patient survival rates.

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

RNA Secondary Structure Modeling Following the IPANEMAP Workflow

The structure of RNA molecules and their complexes are crucial for understanding biology at the molecular level. Resolving these structures holds the key to understanding their manifold structure-mediated functions ranging from regulating gene expression to catalyzing biochemical processes. Predicting RNA secondary structure is a prerequisite and a key step to accurately model their three dimensional structure. Although dedicated modelling software are making fast and significant progresses, predicting an accurate secondary structure from the sequence remains a challenge. Their performance can be significantly improved by the incorporation of experimental RNA structure probing data. Many different chemical and enzymatic probes have been developed; however, only one set of quantitative data can be incorporated as constraints for computer-assisted modelling. IPANEMAP is a recent workflow based on RNAfold that can take into account several quantitative or qualitative data sets to model RNA secondary structure. This chapter details the methods for popular chemical probing (DMS, CMCT, SHAPE-CE, and SHAPE-Map) and the subsequent analysis and structure prediction using IPANEMAP.

Engineering U1-Based Tetracycline-Inducible Riboswitches to Control Gene Expression in Mammals

Synthetic riboswitches are promising regulatory devices due to their small size, lack of immunogenicity, and ability to fine-tune gene expression in the absence of exogenous trans-acting factors. Based on a gene inhibitory system developed at our lab, termed U1snRNP interference (U1i), we developed tetracycline (TC)-inducible riboswitches that modulate mRNA polyadenylation through selective U1 snRNP recruitment. First, we engineered different TC-U1i riboswitches, which repress gene expression unless TC is added, leading to inductions of gene expression of 3-to-4-fold. Second, we developed a technique called Systematic Evolution of Riboswitches by Exponential Enrichment (SEREX), to isolate riboswitches with enhanced U1 snRNP binding capacity and activity, achieving inducibilities of up to 8-fold. Interestingly, by multiplexing riboswitches we increased inductions up to 37-fold. Finally, we demonstrated that U1i-based riboswitches are dose-dependent and reversible and can regulate the expression of reporter and endogenous genes in culture cells and mouse models, resulting in attractive systems for gene therapy applications. Our work probes SEREX as a much-needed technology for the in vitro identification of riboswitches capable of regulating gene expression in vivo.

Assessment of three-dimensional RNA structure prediction in CASP15

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty-two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and x-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as noncanonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

RNAvigate: Efficient exploration of RNA chemical probing datasets

Chemical probing technologies enable high-throughput examination of diverse structural features of RNA including local nucleotide flexibility, RNA secondary structure, protein- and ligand-binding, through-space interaction networks, and multi-state structural ensembles. Performing these experiments, by themselves, does not directly lead to biological insight. Instead, deep understanding of RNA structure-function relationships typically requires evaluating a system under structure- and function-altering conditions, linking these data with additional information, and visualizing multi-layered relationships. Current platforms lack the broad accessibility, flexibility, and efficiency needed to iterate on integrative analyses of these diverse, complex data. Here, we share the RNA visualization and graphical analysis toolset RNAvigate, a straightforward and flexible Python library. RNAvigate currently automatically parses twenty-one standard file formats (primary sequence annotations, per- and internucleotide data, and secondary and tertiary structures) and outputs eighteen plot types. These features enable efficient exploration of nuanced relationships between chemical probing data, RNA structure, and motif annotations across multiple experimental samples. Compatibility with Jupyter Notebooks enables non-burdensome, reproducible, transparent and organized sharing of multi-step analyses and data visualization strategies. RNAvigate simplifies examination of multi-layered RNA structure information and accelerates discovery and characterization of RNA-centric functions in biology.

Observation of coordinated RNA folding events by systematic cotranscriptional RNA structure probing

RNA begins to fold as it is transcribed by an RNA polymerase. Consequently, RNA folding is constrained by the direction and rate of transcription. Understanding how RNA folds into secondary and tertiary structures therefore requires methods for determining the structure of cotranscriptional folding intermediates. Cotranscriptional RNA chemical probing methods accomplish this by systematically probing the structure of nascent RNA that is displayed from an RNA polymerase. Here, we describe a concise, high-resolution cotranscriptional RNA chemical probing procedure called variable length Transcription Elongation Complex RNA structure probing (TECprobe-VL). We demonstrate the accuracy and resolution of TECprobe-VL by replicating and extending previous analyses of ZTP and fluoride riboswitch folding and mapping the folding pathway of a ppGpp-sensing riboswitch. In each system, we show that TECprobe-VL identifies coordinated cotranscriptional folding events that mediate transcription antitermination. Our findings establish TECprobe-VL as an accessible method for mapping cotranscriptional RNA folding pathways.

Integrative analysis reveals that SLC38A1 promotes hepatocellular carcinoma development via PI3K/AKT/mTOR signaling via glutamine mediated energy metabolism

Although hepatocellular carcinoma (HCC) is rather frequent, little is known about the molecular pathways underlying its development, progression, and prognosis. In the current study, we comprehensively analyzed the deferentially expressed metabolism-related genes (MRGs) in HCC based on TCGA datasets attempting to discover the potentially prognostic genes in HCC. The up-regulated MRGs were further subjected to analyze their prognostic values and protein expressions. Twenty-seven genes were identified because their high expressions were significant in OS, PFS, DFS, DSS, and HCC tumor samples. They were then used for GO, KEGG, methylation, genetics changes, immune infiltration analyses. Moreover, we established a prognostic model in HCC using univariate assays and LASSO regression based on these MRGs. Additionally, we also found that SLC38A1, an amino acid metabolism closely related transporter, was a potential prognostic gene in HCC, and its function in HCC was further studied using experiments. We found that the knockdown of SLC38A1 notably suppressed the growth and migration of HCC cells. Further studies revealed that SLC38A1 modulated the development of HCC cells by regulating PI3K/AKT/mTOR signaling via glutamine mediated energy metabolism. In conclusion, this study identified the potentially prognostic MRGs in HCC and uncovered that SLC38A1 regulated HCC development and progression by regulating PI3K/AKT/mTOR signaling via glutamine mediated energy metabolism, which might provide a novel marker and potential therapeutic target in HCC.

TransUFold: Unlocking the structural complexity of short and long RNA with pseudoknots

The RNA secondary structure is like a blueprint that holds the key to unlocking the mysteries of RNA function and 3D structure. It serves as a crucial foundation for investigating the complex world of RNA, making it an indispensable component of research in this exciting field. However, pseudoknots cannot be accurately predicted by conventional prediction methods based on free energy minimization, which results in a performance bottleneck. To this end, we propose a deep learning-based method called TransUFold to train directly on RNA data annotated with structure information. It employs an encoder-decoder network architecture, named Vision Transformer, to extract long-range interactions in RNA sequences and utilizes convolutions with lateral connections to supplement short-range interactions. Then, a post-processing program is designed to constrain the model's output to produce realistic and effective RNA secondary structures, including pseudoknots. After training TransUFold on benchmark datasets, we outperform other methods in test data on the same family. Additionally, we achieve better results on longer sequences up to 1600 nt, demonstrating the outstanding performance of Vision Transformer in extracting long-range interactions in RNA sequences. Finally, our analysis indicates that TransUFold produces effective pseudoknot structures in long sequences. As more high-quality RNA structures become available, deep learning-based prediction methods like Vision Transformer can exhibit better performance.

Assessment of three-dimensional RNA structure prediction in CASP15

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and X-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as non-canonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

Dissecting the Molecular Mechanisms Driving Electropathology in Atrial Fibrillation: Deployment of RNA Sequencing and Transcriptomic Analyses

Despite many efforts to treat atrial fibrillation (AF), the most common progressive and age-related cardiac tachyarrhythmia in the Western world, the efficacy is still suboptimal. A plausible reason for this is that current treatments are not directed at underlying molecular root causes that drive electrical conduction disorders and AF (i.e., electropathology). Insights into AF-induced transcriptomic alterations may aid in a deeper understanding of electropathology. Specifically, RNA sequencing (RNA-seq) facilitates transcriptomic analyses and discovery of differences in gene expression profiles between patient groups. In the last decade, various RNA-seq studies have been conducted in atrial tissue samples of patients with AF versus controls in sinus rhythm. Identified differentially expressed molecular pathways so far include pathways related to mechanotransduction, ECM remodeling, ion channel signaling, and structural tissue organization through developmental and inflammatory signaling pathways. In this review, we provide an overview of the available human AF RNA-seq studies and highlight the molecular pathways identified. Additionally, a comparison is made between human RNA-seq findings with findings from experimental AF model systems and we discuss contrasting findings. Finally, we elaborate on new exciting RNA-seq approaches, including single-nucleotide variants, spatial transcriptomics and profiling of different populations of total RNA, small RNA and long non-coding RNA.

GCNfold: A novel lightweight model with valid extractors for RNA secondary structure prediction

RNA secondary structure is essential for predicting the tertiary structure and understanding RNA function. Recent research tends to stack numerous modules to design large deep-learning models. This can increase the accuracy to more than 70%, as well as significant training costs and prediction efficiency. We proposed a model with three feature extractors called GCNfold. Structure Extractor utilizes a three-layer Graph Convolutional Network (GCN) to mine the structural information of RNA, such as stems, hairpin, and internal loops. Structure and Sequence Fusion embeds structural information into sequences with Transformer Encoders. Long-distance Dependency Extractor captures long-range pairwise relationships by UNet. The experiments indicate that GCNfold has a small number of parameters, a fast inference speed, and a high accuracy among all models with over 80% accuracy. Additionally, GCNfold-Small takes only 90ms to infer an RNA secondary structure and can achieve close to 90% accuracy on average. The GCNfold code is available on Github https://github.com/EnbinYang/GCNfold.

Led-Seq: ligation-enhanced double-end sequence-based structure analysis of RNA

Structural analysis of RNA is an important and versatile tool to investigate the function of this type of molecules in the cell as well as in vitro. Several robust and reliable procedures are available, relying on chemical modification inducing RT stops or nucleotide misincorporations during reverse transcription. Others are based on cleavage reactions and RT stop signals. However, these methods address only one side of the RT stop or misincorporation position. Here, we describe Led-Seq, a new approach based on lead-induced cleavage of unpaired RNA positions, where both resulting cleavage products are investigated. The RNA fragments carrying 2', 3'-cyclic phosphate or 5'-OH ends are selectively ligated to oligonucleotide adapters by specific RNA ligases. In a deep sequencing analysis, the cleavage sites are identified as ligation positions, avoiding possible false positive signals based on premature RT stops. With a benchmark set of transcripts in Escherichia coli, we show that Led-Seq is an improved and reliable approach based on metal ion-induced phosphodiester hydrolysis to investigate RNA structures in vivo.

A systematic survey of LU domain-containing proteins reveals a novel human gene, LY6A, which encodes the candidate ortholog of mouse Ly-6A/Sca-1 and is aberrantly expressed in pituitary tumors

The Ly-6 and uPAR (LU) domain-containing proteins represent a large family of cell-surface markers. In particular, mouse Ly-6A/Sca-1 is a widely used marker for various stem cells; however, its human ortholog is missing. In this study, based on a systematic survey and comparative genomic study of mouse and human LU domain-containing proteins, we identified a previously unannotated human gene encoding the candidate ortholog of mouse Ly-6A/Sca-1. This gene, hereby named LY6A, reversely overlaps with a lncRNA gene in the majority of exonic sequences. We found that LY6A is aberrantly expressed in pituitary tumors, but not in normal pituitary tissues, and may contribute to tumorigenesis. Similar to mouse Ly-6A/Sca-1, human LY6A is also upregulated by interferon, suggesting a conserved transcriptional regulatory mechanism between humans and mice. We cloned the full-length LY6A cDNA, whose encoded protein sequence, domain architecture, and exon-intron structures are all well conserved with mouse Ly-6A/Sca-1. Ectopic expression of the LY6A protein in cells demonstrates that it acts the same as mouse Ly-6A/Sca-1 in their processing and glycosylphosphatidylinositol anchoring to the cell membrane. Collectively, these studies unveil a novel human gene encoding a candidate biomarker and provide an interesting model gene for studying gene regulatory and evolutionary mechanisms.

Plant organellar RNA maturation

Plant organellar RNA metabolism is run by a multitude of nucleus-encoded RNA-binding proteins (RBPs) that control RNA stability, processing, and degradation. In chloroplasts and mitochondria, these post-transcriptional processes are vital for the production of a small number of essential components of the photosynthetic and respiratory machinery-and consequently for organellar biogenesis and plant survival. Many organellar RBPs have been functionally assigned to individual steps in RNA maturation, often specific to selected transcripts. While the catalog of factors identified is ever-growing, our knowledge of how they achieve their functions mechanistically is far from complete. This review summarizes the current knowledge of plant organellar RNA metabolism taking an RBP-centric approach and focusing on mechanistic aspects of RBP functions and the kinetics of the processes they are involved in.

Targeting miRNAs and Other Non-Coding RNAs as a Therapeutic Approach: An Update

Since the discovery of the first microRNAs (miRNAs, miRs), the understanding of miRNA biology has expanded substantially. miRNAs are involved and described as master regulators of the major hallmarks of cancer, including cell differentiation, proliferation, survival, the cell cycle, invasion, and metastasis. Experimental data indicate that cancer phenotypes can be modified by targeting miRNA expression, and because miRNAs act as tumor suppressors or oncogenes (oncomiRs), they have emerged as attractive tools and, more importantly, as a new class of targets for drug development in cancer therapeutics. With the use of miRNA mimics or molecules targeting miRNAs (i.e., small-molecule inhibitors such as anti-miRS), these therapeutics have shown promise in preclinical settings. Some miRNA-targeted therapeutics have been extended to clinical development, such as the mimic of miRNA-34 for treating cancer. Here, we discuss insights into the role of miRNAs and other non-coding RNAs in tumorigenesis and resistance and summarize some recent successful systemic delivery approaches and recent developments in miRNAs as targets for anticancer drug development. Furthermore, we provide a comprehensive overview of mimics and inhibitors that are in clinical trials and finally a list of clinical trials based on miRNAs.

Single-Molecule Correlated Chemical Probing: A Revolution in RNA Structure Analysis

RNA molecules convey biological information both in their linear sequence and in their base-paired secondary and tertiary structures. Chemical probing experiments, which involve treating an RNA with a reagent that modifies conformationally dynamic nucleotides, have broadly enabled examination of short- and long-range RNA structure in diverse contexts, including in living cells. For decades, chemical probing experiments have been interpreted in a per-nucleotide way, such that the reactivity measured at each nucleotide reports the average structure at a position over all RNA molecules within a sample. However, there are numerous important cases where per-nucleotide chemical probing falls short, including for RNAs that are bound by proteins, RNAs that form complex higher order structures, and RNAs that sample multiple conformations.Recent experimental and computational innovations have started a revolution in RNA structure analysis by transforming chemical probing into a massively parallel, single-molecule experiment. Enabled by a specialized reverse transcription strategy called mutational profiling (MaP), multiple chemical modification events can be measured within individual RNA molecules. Nucleotides that communicate structurally through direct base pairing or large-scale folding-unfolding transitions will react with chemical probes in a correlated manner, thereby revealing structural complexity hidden to conventional approaches. These single-molecule correlated chemical probing (smCCP) experiments can be interpreted to directly identify nucleotides that base pair (the PAIR-MaP strategy) and to reveal long-range, through-space structural communication (RING-MaP). Correlated probing can also define the thermodynamic populations of complex RNA ensembles (DANCE-MaP). Complex RNA-protein networks can be interrogated by cross-linking proteins to RNA and measuring correlations between cross-linked positions (RNP-MaP).smCCP thus visualizes RNA secondary and higher-order structure with unprecedented accuracy, defining novel structures, RNA-protein interaction networks, time-resolved dynamics, and allosteric structural switches. These strategies are not mutually exclusive; in favorable cases, multiple levels of RNA structure ─ base pairing, through-space structural communication, and equilibrium ensembles ─ can be resolved concurrently. The physical experimentation required for smCCP is profoundly simple, and experiments are readily performed in cells on RNAs of any size, including large noncoding RNAs and mRNAs. Single-molecule correlated chemical probing is paving the way for a new generation of biophysical studies on RNA in living systems.

RNA:RNA interaction in ternary complexes resolved by chemical probing

RNA regulation can be performed by a second targeting RNA molecule, such as in the microRNA regulation mechanism. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) probes the structure of RNA molecules and can resolve RNA:protein interactions, but RNA:RNA interactions have not yet been addressed with this technique. Here, we apply SHAPE to investigate RNA-mediated binding processes in RNA:RNA and RNA:RNA-RBP complexes. We use RNA:RNA binding by SHAPE (RABS) to investigate microRNA-34a (miR-34a) binding its mRNA target, the silent information regulator 1 (mSIRT1), both with and without the Argonaute protein, constituting the RNA-induced silencing complex (RISC). We show that the seed of the mRNA target must be bound to the microRNA loaded into RISC to enable further binding of the compensatory region by RISC, while the naked miR-34a is able to bind the compensatory region without seed interaction. The method presented here provides complementary structural evidence for the commonly performed luciferase-assay-based evaluation of microRNA binding-site efficiency and specificity on the mRNA target site and could therefore be used in conjunction with it. The method can be applied to any nucleic acid-mediated RNA- or RBP-binding process, such as splicing, antisense RNA binding, or regulation by RISC, providing important insight into the targeted RNA structure.

Probing the dynamic RNA structurome and its functions

RNA is a key regulator of almost every cellular process, and the structures adopted by RNA molecules are thought to be central to their functions. The recent fast-paced evolution of high-throughput sequencing-based RNA structure mapping methods has enabled the rapid in vivo structural interrogation of entire cellular transcriptomes. Collectively, these studies are shedding new light on the long underestimated complexity of the structural organization of the transcriptome - the RNA structurome. Moreover, recent analyses are challenging the view that the RNA structurome is a static entity by revealing how RNA molecules establish intricate networks of alternative intramolecular and intermolecular interactions and that these ensembles of RNA structures are dynamically regulated to finely tune RNA functions in living cells. This new understanding of how RNA can shape cell phenotypes has important implications for the development of RNA-targeted therapeutic strategies.

Observation of coordinated cotranscriptional RNA folding events

RNA begins to fold as it is transcribed by an RNA polymerase. Consequently, RNA folding is constrained by the direction and rate of transcription. Understanding how RNA folds into secondary and tertiary structures therefore requires methods for determining the structure of cotranscriptional folding intermediates. Cotranscriptional RNA chemical probing methods accomplish this by systematically probing the structure of nascent RNA that is displayed from RNA polymerase. Here, we have developed a concise, high-resolution cotranscriptional RNA chemical probing procedure called Transcription Elongation Complex RNA structure probing-Multilength (TECprobe-ML). We validated TECprobe-ML by replicating and extending previous analyses of ZTP and fluoride riboswitch folding, and mapped the folding pathway of a ppGpp-sensing riboswitch. In each system, TECprobe-ML identified coordinated cotranscriptional folding events that mediate transcription antitermination. Our findings establish TECprobe-ML as an accessible method for mapping cotranscriptional RNA folding pathways.

Biophysical characterization of the structure of a SARS-CoV-2 self-amplifying RNA (saRNA) vaccine

The current SARS-Covid-2 (SARS-CoV-2) pandemic has led to an acceleration of messenger ribonucleic acid (mRNA) vaccine technology. The development of production processes for these large mRNA molecules, especially self-amplifying mRNA (saRNA), has required concomitant development of analytical characterization techniques. Characterizing the purity, shape and structure of these biomolecules is key to their successful performance as drug products. This article describes the biophysical characterization of the Imperial College London Self-amplifying viral RNA vaccine (IMP-1) developed for SARS-CoV-2. A variety of analytical techniques have been used to characterize the IMP-1 RNA molecule. In this article, we use ultraviolet spectroscopy, dynamic light scattering, size-exclusion chromatography small-angle X-ray scattering and circular dichroism to determine key biophysical attributes of IMP-1. Each technique provides important information about the concentration, size, shape, structure and purity of the molecule.