Probing the dynamic RNA structurome and its functions

Mechanistic analysis of Riboswitch Ligand interactions provides insights into pharmacological control over gene expression

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

RNA dynamics from experimental and computational approaches

Conformational dynamics is crucial for the biological function of RNA molecules and for their potential as therapeutic targets. This meeting report outlines key "take-home" messages that emerged from the presentations and discussions during the CECAM workshop "RNA dynamics from experimental and computational approaches" in Paris, June 26-28, 2023.

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.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

Pulsed Dipolar EPR for Self-Limited Complexes of Oligonucleotides Studies

Pulsed electron-electron double resonance (PELDOR) spectroscopy is a powerful method for determining nucleic acid (NA) structure and conformational dynamics. PELDOR with molecular dynamics (MD) simulations opens up unique possibilities for defining the conformational ensembles of flexible, three-dimensional, self-assembled complexes of NA. Understanding the diversity and structure of these complexes is vital for uncovering matrix and regulative biological processes in the human body and artificially influencing them for therapeutic purposes. To explore the reliability of PELDOR and MD simulations, we site-specifically attached nitroxide spin labels to oligonucleotides, which form self-assembled complexes between NA chains and exhibit significant conformational flexibility. The DNA complexes assembled from a pair of oligonucleotides with different linker sizes showed excellent agreement between the distance distributions obtained from PELDOR and calculated from MD simulations, both for the mean inter-spin distance and the distance distribution width. These results prove that PELDOR with MD simulations has significant potential for studying the structure and dynamics of conformational flexible complexes of NA.

RNA structure in alternative splicing regulation: from mechanism to therapy

Alternative splicing is a highly intricate process that plays a crucial role in post-transcriptional regulation and significantly expands the functional proteome of a limited number of coding genes in eukaryotes. Its regulation is multifactorial, with RNA structure exerting a significant impact. Aberrant RNA conformations lead to dysregulation of splicing patterns, which directly affects the manifestation of disease symptoms. In this review, the molecular mechanisms of RNA secondary structure-mediated splicing regulation are summarized, with a focus on the complex interplay between aberrant RNA conformations and disease phenotypes resulted from splicing defects. This study also explores additional factors that reshape structural conformations, enriching our understanding of the mechanistic network underlying structure-mediated splicing regulation. In addition, an emphasis has been placed on the clinical role of targeting aberrant splicing corrections in human diseases. The principal mechanisms of action behind this phenomenon are described, followed by a discussion of prospective development strategies and pertinent challenges.

Unraveling the complexity: Advanced methods in analyzing DNA, RNA, and protein interactions

Exploring the intricate interplay within and between nucleic acids, as well as their interactions with proteins, holds pivotal significance in unraveling the molecular complexities steering cancer initiation and progression. To investigate these interactions, a diverse array of highly specific and sensitive molecular techniques has been developed. The selection of a particular technique depends on the specific nature of the interactions. Typically, researchers employ an amalgamation of these different techniques to obtain a comprehensive and holistic understanding of inter- and intramolecular interactions involving DNA-DNA, RNA-RNA, DNA-RNA, or protein-DNA/RNA. Examining nucleic acid conformation reveals alternative secondary structures beyond conventional ones that have implications for cancer pathways. Mutational hotspots in cancer often lie within sequences prone to adopting these alternative structures, highlighting the importance of investigating intra-genomic and intra-transcriptomic interactions, especially in the context of mutations, to deepen our understanding of oncology. Beyond these intramolecular interactions, the interplay between DNA and RNA leads to formations like DNA:RNA hybrids (known as R-loops) or even DNA:DNA:RNA triplex structures, both influencing biological processes that ultimately impact cancer. Protein-nucleic acid interactions are intrinsic cellular phenomena crucial in both normal and pathological conditions. In particular, genetic mutations or single amino acid variations can alter a protein's structure, function, and binding affinity, thus influencing cancer progression. It is thus, imperative to understand the differences between wild-type (WT) and mutated (MT) genes, transcripts, and proteins. The review aims to summarize the frequently employed methods and techniques for investigating interactions involving nucleic acids and proteins, highlighting recent advancements and diverse adaptations of each technique.

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.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

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.

RNATACs: Multispecific small molecules targeting RNA by induced proximity

RNA-targeting small molecules (rSMs) have become an attractive modality to tackle traditionally undruggable proteins and expand the druggable space. Among many innovative concepts, RNA-targeting chimeras (RNATACs) represent a new class of multispecific, induced proximity small molecules that act by chemically bringing RNA targets into proximity with an endogenous RNA effector, such as a ribonuclease (RNase). Depending on the RNA effector, RNATACs can alter the stability, localization, translation, or splicing of the target RNA. Although still in its infancy, this new modality has the potential for broad applications in the future to treat diseases with high unmet need. In this review, we discuss potential advantages of RNATACs, recent progress in the field, and challenges to this cutting-edge technology.

Selective Arylation of RNA 2'-OH Groups via SNAr Reaction with Trialkylammonium Heterocycles

Small-molecule reactions at the 2'-OH groups of RNA enable useful applications for transcriptome technology and biology. To date, all reactions have involved carbonyl acylation and mechanistically related sulfonylation, limiting the types of modifications and properties that can be achieved. Here we report that electron-deficient heteroaryl species selectively react with 2'-OH groups of RNA in water via SNAr chemistry. In particular, trialkyl-ammonium (TAA)-activated aromatic heterocycles, prepared in one step from aryl chloride precursors, give high conversions to aryl ether adducts with RNAs in aqueous buffer in ~2-3 h. Remarkably, a TAA triazine previously used only for reaction with carboxylic acids, shows unprecedented selectivity for RNA over water, reacting rapidly with 2'-OH groups while exhibiting a half-life in water of >10 days. We further show that a triazine aryl species can be used as a probe at trace-level yields to map RNA structure in vitro. Finally, we prepare a number of functionalized trialkylammonium triazine reagents and show that they can be used to covalently label RNA efficiently for use in vitro and in living cells. This direct arylation chemistry offers a simple and distinct structural scaffold for post-synthetic RNA modification, with potential utility in multiple applications in transcriptome research.

CD3 aptamers promote expansion and persistence of tumor-reactive T cells for adoptive T cell therapy in cancer

The CD3/T cell receptor (TCR) complex is responsible for antigen-specific pathogen recognition by T cells, and initiates the signaling cascade necessary for activation of effector functions. CD3 agonistic antibodies are commonly used to expand T lymphocytes in a wide range of clinical applications, including in adoptive T cell therapy for cancer patients. A major drawback of expanding T cell populations ex vivo using CD3 agonistic antibodies is that they expand and activate T cells independent of their TCR antigen specificity. Therapeutic agents that facilitate expansion of T cells in an antigen-specific manner and reduce their threshold of T cell activation are therefore of great interest for adoptive T cell therapy protocols. To identify CD3-specific T cell agonists, several RNA aptamers were selected against CD3 using Systematic Evolution of Ligands by EXponential enrichment combined with high-throughput sequencing. The extent and specificity of aptamer binding to target CD3 were assessed through surface plasma resonance, P32 double-filter assays, and flow cytometry. Aptamer-mediated modulation of the threshold of T cell activation was observed in vitro and in preclinical transgenic TCR mouse models. The aptamers improved efficacy and persistence of adoptive T cell therapy by low-affinity TCR-reactive T lymphocytes in melanoma-bearing mice. Thus, CD3-specific aptamers can be applied as therapeutic agents which facilitate the expansion of tumor-reactive T lymphocytes while conserving their tumor specificity. Furthermore, selected CD3 aptamers also exhibit cross-reactivity to human CD3, expanding their potential for clinical translation and application in the future.

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.

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.

Isoform-specific RNA structure determination using Nano-DMS-MaP

RNA structure determination is essential to understand how RNA carries out its diverse biological functions. In cells, RNA isoforms are readily expressed with partial variations within their sequences due, for example, to alternative splicing, heterogeneity in the transcription start site, RNA processing or differential termination/polyadenylation. Nanopore dimethyl sulfate mutational profiling (Nano-DMS-MaP) is a method for in situ isoform-specific RNA structure determination. Unlike similar methods that rely on short sequencing reads, Nano-DMS-MaP employs nanopore sequencing to resolve the structures of long and highly similar RNA molecules to reveal their previously hidden structural differences. This Protocol describes the development and applications of Nano-DMS-MaP and outlines the main considerations for designing and implementing a successful experiment: from bench to data analysis. In cell probing experiments can be carried out by an experienced molecular biologist in 3-4 d. Data analysis requires good knowledge of command line tools and Python scripts and requires a further 3-5 d.

Mitochondrial RNA Helicases: Key Players in the Regulation of Plant Organellar RNA Splicing and Gene Expression

Mitochondrial genomes of land plants are large and exhibit a complex mode of gene organization and expression, particularly at the post-transcriptional level. The primary organellar transcripts in plants undergo extensive maturation steps, including endo- and/or exo-nucleolytic cleavage, RNA-base modifications (mostly C-to-U deaminations) and both 'cis'- and 'trans'-splicing events. These essential processing steps rely on the activities of a large set of nuclear-encoded factors. RNA helicases serve as key players in RNA metabolism, participating in the regulation of transcription, mRNA processing and translation. They unwind RNA secondary structures and facilitate the formation of ribonucleoprotein complexes crucial for various stages of gene expression. Furthermore, RNA helicases are involved in RNA metabolism by modulating pre-mRNA maturation, transport and degradation processes. These enzymes are, therefore, pivotal in RNA quality-control mechanisms, ensuring the fidelity and efficiency of RNA processing and turnover in plant mitochondria. This review summarizes the significant roles played by helicases in regulating the highly dynamic processes of mitochondrial transcription, RNA processing and translation in plants. We further discuss recent advancements in understanding how dysregulation of mitochondrial RNA helicases affects the splicing of organellar genes, leading to respiratory dysfunctions, and consequently, altered growth, development and physiology of land plants.

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.

The application of single-molecule optical tweezers to study disease-related structural dynamics in RNA

RNA, a dynamic and flexible molecule with intricate three-dimensional structures, has myriad functions in disease development. Traditional methods, such as X-ray crystallography and nuclear magnetic resonance, face limitations in capturing real-time, single-molecule dynamics crucial for understanding RNA function. This review explores the transformative potential of single-molecule force spectroscopy using optical tweezers, showcasing its capability to directly probe time-dependent structural rearrangements of individual RNA molecules. Optical tweezers offer versatility in exploring diverse conditions, with the potential to provide insights into how environmental changes, ligands and RNA-binding proteins impact RNA behaviour. By enabling real-time observations of large-scale structural dynamics, optical tweezers emerge as an invaluable tool for advancing our comprehension of RNA structure and function. Here, we showcase their application in elucidating the dynamics of RNA elements in virology, such as the pseudoknot governing ribosomal frameshifting in SARS-CoV-2.

RNA three-dimensional structure drives the sequence organization of potato spindle tuber viroid quasispecies

RNA viruses and viroids exist and evolve as quasispecies due to error-prone replication. Quasispecies consist of a few dominant master sequences alongside numerous variants that contribute to genetic diversity. Upon environmental changes, certain variants within quasispecies have the potential to become the dominant sequences, leading to the emergence of novel infectious strains. However, the emergence of new infectious variants remains unpredictable. Using mutant pools prepared by saturation mutagenesis of selected stem and loop regions, our study of potato spindle tuber viroid (PSTVd) demonstrates that mutants forming local three-dimensional (3D) structures similar to the wild type (WT) are more likely to accumulate in PSTVd quasispecies. The selection mechanisms underlying this biased accumulation are likely associated with cell-to-cell movement and long-distance trafficking. Moreover, certain trafficking-defective PSTVd mutants can be spread by functional sister genomes in the quasispecies. Our study reveals that the RNA 3D structure of stems and loops constrains the evolution of viroid quasispecies. Mutants with a structure similar to WT have a higher likelihood of being maintained within the quasispecies and can potentially give rise to novel infectious variants. These findings emphasize the potential of targeting RNA 3D structure as a more robust approach to defend against viroid infections.

rnaCrosslinkOO: an object-oriented R package for the analysis of RNA structural data generated by RNA crosslinking experiments

SUMMARY

RNA (ribonucleic acid) molecules have secondary and tertiary structures in vivo which play a crucial role in cellular processes such as the regulation of gene expression, RNA processing and localization. The ability to investigate these structures will enhance our understanding of their function and contribute to the diagnosis and treatment of diseases caused by RNA dysregulation. However, there are no mature pipelines or packages for processing and analyzing complex in vivo RNA structural data. Here, we present rnaCrosslinkOO (RNA Crosslink Object-Oriented), a novel software package for the comprehensive analysis of data derived from the COMRADES (Crosslinking of Matched RNA and Deep Sequencing) method. rnaCrosslinkOO offers a comprehensive pipeline from raw sequencing reads to the identification and comparison of RNA structural features. It includes read processing and alignment, clustering of duplexes, data exploration, folding and comparisons of RNA structures. rnaCrosslinkOO also enables comparisons between conditions, the identification of inter-RNA interactions, and the incorporation of reactivity data to improve structure prediction.

AVAILABILITY AND IMPLEMENTATION

rnaCrosslinkOO is freely available to noncommercial users and implemented in R, with the source code and documentation accessible at https://CRAN.R-project.org/package=rnaCrosslinkOO. The software is supported on Linux, macOS, and Windows platforms.

Transient states during the annealing of mismatched and bulged oligonucleotides

Oligonucleotide hybridization is crucial in various biological, prebiotic and nanotechnological processes, including gene regulation, non-enzymatic primer extension and DNA nanodevice assembly. Although extensive research has focused on the thermodynamics and kinetics of nucleic acid hybridization, the behavior of complex mixtures and the outcome of competition for target binding remain less well understood. In this study, we investigate the impact of mismatches and bulges in a 12 bp DNA or RNA duplex on its association (kon) and dissociation (koff) kinetics. We find that such defects have relatively small effects on the association kinetics, while the dissociation kinetics vary in a position-dependent manner by up to 6 orders of magnitude. Building upon this observation, we explored a competition scenario involving multiple oligonucleotides, and observed a transient low specificity of probe hybridization to fully versus partially complementary targets in solution. We characterize these long-lived metastable states and their evolution toward equilibrium, and show that sufficiently long-lived mis-paired duplexes can serve as substrates for prebiotically relevant chemical copying reactions. Our results suggest that transient low accuracy states may spontaneously emerge within all complex nucleic acid systems comprising a large enough number of competing strands, with potential repercussions for gene regulation in the realm of modern biology and the prebiotic preservation of genetic information.

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.

Accurate prediction of RNA secondary structure including pseudoknots through solving minimum-cost flow with learned potentials

Pseudoknots are key structure motifs of RNA and pseudoknotted RNAs play important roles in a variety of biological processes. Here, we present KnotFold, an accurate approach to the prediction of RNA secondary structure including pseudoknots. The key elements of KnotFold include a learned potential function and a minimum-cost flow algorithm to find the secondary structure with the lowest potential. KnotFold learns the potential from the RNAs with known structures using an attention-based neural network, thus avoiding the inaccuracy of hand-crafted energy functions. The specially designed minimum-cost flow algorithm used by KnotFold considers all possible combinations of base pairs and selects from them the optimal combination. The algorithm breaks the restriction of nested base pairs required by the widely used dynamic programming algorithms, thus enabling the identification of pseudoknots. Using 1,009 pseudoknotted RNAs as representatives, we demonstrate the successful application of KnotFold in predicting RNA secondary structures including pseudoknots with accuracy higher than the state-of-the-art approaches. We anticipate that KnotFold, with its superior accuracy, will greatly facilitate the understanding of RNA structures and functionalities.

Small RNA structural biochemistry in a post-sequencing era

High-throughput sequencing has had an enormous impact on small RNA research during the past decade. However, sequencing only offers a one-dimensional view of the transcriptome and is often highly biased. Additionally, the 'sequence, map and annotate' approach, used widely in small RNA research, can lead to flawed interpretations of the data, lacking biological plausibility, due in part to database issues. Even in the absence of technical biases, the loss of three-dimensional information is a major limitation to understanding RNA stability, turnover and function. For example, noncoding RNA-derived fragments seem to exist mainly as dimers, tetramers or as nicked forms of their parental RNAs, contrary to widespread assumptions. In this perspective, we will discuss main sources of bias during small RNA-sequencing, present several useful bias-reducing strategies and provide guidance on the interpretation of small RNA-sequencing results, with emphasis on RNA fragmentomics. As sequencing offers a one-dimensional projection of a four-dimensional reality, prior structure-level knowledge is often needed to make sense of the data. Consequently, while less-biased sequencing methods are welcomed, integration of orthologous experimental techniques is also strongly recommended.

Mechanistic Analysis of Riboswitch Ligand Interactions Provides Insights into Pharmacological Control over Gene Expression

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Collision-Induced Unfolding Reveals Disease-Associated Stability Shifts in Mitochondrial Transfer Ribonucleic Acids

Ribonucleic acids (RNAs) remain challenging targets for structural biology, creating barriers to understanding their vast functions in cellular biology and fully realizing their applications in biotechnology. The inherent dynamism of RNAs creates numerous obstacles in capturing their biologically relevant higher-order structures (HOSs), and as a result, many RNA functions remain unknown. In this study, we describe the development of native ion mobility-mass spectrometry and collision-induced unfolding (CIU) for the structural characterization of a variety of RNAs. We evaluate the ability of these techniques to preserve native structural features in the gas phase across a wide range of functional RNAs. Finally, we apply these tools to study the elusive mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes-associated A3243G mutation. Our data demonstrate that our experimentally determined conditions preserve some solution-state memory of RNAs via the correlated complexity of CIU fingerprints and RNA HOS, the observation of predicted stability shifts in the control RNA samples, and the retention of predicted magnesium binding events in gas-phase RNA ions. Significant differences in collision cross section and stability are observed as a function of the A3243G mutation across a subset of the mitochondrial tRNA maturation pathway. We conclude by discussing the potential application of CIU for the development of RNA-based biotherapeutics and, more broadly, transcriptomic characterization.

RNA structure: implications in viral infections and neurodegenerative diseases

RNA is an intermediary between DNA and protein, a catalyzer of biochemical reactions, and a regulator of genes and transcripts. RNA structures are essential for complicated functions. Recent years have witnessed rapid advancements in RNA secondary structure probing techniques. These technological strides provided comprehensive insights into RNA structures, which significantly contributed to our understanding of diverse cellular regulatory processes, including gene regulation, epigenetic regulation, and post-transactional regulation. Meanwhile, they have facilitated the creation of therapeutic tools for tackling human diseases. Despite their therapeutic applications, RNA structure probing methods also offer a promising avenue for exploring the mechanisms of human diseases, potentially providing the key to overcoming existing research constraints and obtaining the in-depth information necessary for a deeper understanding of disease mechanisms.

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

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

A TeZla micromixer for interrogating the early and broad folding landscape of G-quadruplex via multistage velocity descending

Biomacromolecular folding kinetics involves fast folding events and broad timescales. Current techniques face limitations in either the required time resolution or the observation window. In this study, we developed the TeZla micromixer, integrating Tesla and Zigzag microstructures with a multistage velocity descending strategy. TeZla achieves a significant short mixing dead time (40 µs) and a wide time window covering four orders of magnitude (up to 300 ms). Using this unique micromixer, we explored the folding landscape of c-Myc G4 and its noncanonical-G4 derivatives with different loop lengths or G-vacancy sites. Our findings revealed that c-Myc can bypass folding intermediates and directly adopt a G4 structure in the cation-deficient buffer. Moreover, we found that the loop length and specific G-vacancy site could affect the folding pathway and significantly slow down the folding rates. These results were also cross-validated with real-time NMR and circular dichroism. In conclusion, TeZla represents a versatile tool for studying biomolecular folding kinetics, and our findings may ultimately contribute to the design of drugs targeting G4 structures.

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.

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

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

The myriad roles of RNA structure in the flavivirus life cycle

As positive-sense RNA viruses, the genomes of flaviviruses serve as the template for all stages of the viral life cycle, including translation, replication, and infectious particle production. Yet, they encode just 10 proteins, suggesting that the structure and dynamics of the viral RNA itself helps shepherd the viral genome through these stages. Herein, we highlight advances in our understanding of flavivirus RNA structural elements through the lens of their impact on the viral life cycle. We highlight how RNA structures impact translation, the switch from translation to replication, negative- and positive-strand RNA synthesis, and virion assembly. Consequently, we describe three major themes regarding the roles of RNA structure in flavivirus infections: 1) providing a layer of specificity; 2) increasing the functional capacity; and 3) providing a mechanism to support genome compaction. While the interactions described herein are specific to flaviviruses, these themes appear to extend more broadly across RNA viruses.

Quantifying RNA structures and interactions with a unified reduced chain representation model

RNA is a pivotal molecule that plays critical roles in various cellular processes. Quantifying RNA structures and interactions is essential to understanding RNA function and developing RNA-based therapeutics. Using a unified five-bead model and a non-redundant database, this paper investigates the structural features and interactions of five commonly occurring RNA motifs, i.e., double-stranded helices, hairpin loops, internal/bulge loops, multi-branched junctions, and single-stranded terminal tails. Analyzing detailed distributions of RNA local structural features and base-base interactions reveals a preference for helical structures in both local backbone structures and base orientations. The interactions between adjacent bases exhibit motif-specific and sequence-dependent characteristics, reflecting the distinct topological constraints imposed by different loop-helix connection modes and the varying pairing and stacking interactions among different sequences. These findings shed light on the stability of RNA helices, emphasizing their significance in providing dominant base pairing and stacking interactions for RNA structures and stability. The four non-helix motifs encompass unpaired nucleotide loops and exhibit diverse base-base interactions, contributing to the structural diversity observed in RNA. Overall, the complexity of RNA structure arises from the intricate interplay of base-base interactions.

A Small-Molecule Approach to Bypass In Vitro Selection of New Aptamers: Designer Pre-Ligands Turn Baby Spinach into Sensors for Reactive Inorganic Targets

Fluorescent light-up aptamer (FLAP) systems are promising biosensing platforms that can be genetically encoded. Here, we describe how a single FLAP that works with specific organic ligands can detect multiple, structurally unique, non-fluorogenic, and reactive inorganic targets. We developed 4-O-functionalized benzylidene imidazolinones as pre-ligands with suppressed fluorescent binding interactions with the RNA aptamer Baby Spinach. Inorganic targets, hydrogen sulfide (H2S) or hydrogen peroxide (H2O2), can specifically convert these pre-ligands into the native benzylidene imidazolinones, and thus be detected with Baby Spinach. Adaptation of this approach to live cells opened a new opportunity for top-down construction of whole-cell sensors: Escherichia coli transformed with a Baby Spinach-encoding plasmid and incubated with pre-ligands generated fluorescence in response to exogenous H2S or H2O2. Our approach eliminates the requirement of in vitro selection of a new aptamer sequence for molecular target detection, allows for the detection of short-lived targets, thereby advancing FLAP systems beyond their current capabilities. Leveraging the functional group reactivity of small molecules can lead to cell-based sensors for inorganic molecular targets, exploiting a new synergism between synthetic organic chemistry and synthetic biology.

Effective Antiviral Application of Antisense in Plants by Exploiting Accessible Sites in the Target RNA

Antisense oligodeoxynucleotides (ASOs) have long been used to selectively inhibit or modulate gene expression at the RNA level, and some ASOs are approved for clinical use. However, the practicability of antisense technologies remains limited by the difficulty of reliably predicting the sites accessible to ASOs in complex folded RNAs. Recently, we applied a plant-based method that reproduces RNA-induced RNA silencing in vitro to reliably identify sites in target RNAs that are accessible to small interfering RNA (siRNA)-guided Argonaute endonucleases. Here, we show that this method is also suitable for identifying ASOs that are effective in DNA-induced RNA silencing by RNases H. We show that ASOs identified in this way that target a viral genome are comparably effective in protecting plants from infection as siRNAs with the corresponding sequence. The antiviral activity of the ASOs could be further enhanced by chemical modification. This led to two important conclusions: siRNAs and ASOs that can effectively knock down complex RNA molecules can be identified using the same approach, and ASOs optimized in this way could find application in crop protection. The technology developed here could be useful not only for effective RNA silencing in plants but also in other organisms.

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.

Versatile toolkit for highly-efficient and scarless overexpression of circular RNAs

Circular RNAs (circRNAs) are a class of single-stranded, covalently closed RNA that contain a unique back-splice junction (bsj) sequence created by the ligation of their 5' and 3' ends via spliceosome-catalyzed back-splicing. A key step in illuminating the cellular roles of specific circRNAs is via increasing their expression. This is frequently done by transfecting cells with plasmid DNA containing cloned exons from which the circRNA is transcribed, flanked by sequences that promote back-splicing. We observed that commonly used plasmids lead to the production of circRNAs with molecular scars at the circRNA bsj. Stepwise redesign of the cloning vector corrected this problem, ensuring bona fide circRNAs are produced with their natural bsj at high efficiency. The fidelity of circRNAs produced from this new construct was validated by RNA sequencing and also functionally validated. To increase the utility of this modified resource for expressing circRNA, we developed an expanded set of vectors incorporating this design that (i) enables selection with a variety of antibiotics and fluorescent proteins, (ii) employs a range of promoters varying in promoter strength and (iii) generated a complementary set of lentiviral plasmids for difficult-to-transfect cells. These resources provide a novel and versatile toolkit for high-efficiency and scarless overexpression of circular RNAs that fulfill a critical need for the investigation of circRNA function.

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

When will RNA get its AlphaFold moment?

The protein structure prediction problem has been solved for many types of proteins by AlphaFold. Recently, there has been considerable excitement to build off the success of AlphaFold and predict the 3D structures of RNAs. RNA prediction methods use a variety of techniques, from physics-based to machine learning approaches. We believe that there are challenges preventing the successful development of deep learning-based methods like AlphaFold for RNA in the short term. Broadly speaking, the challenges are the limited number of structures and alignments making data-hungry deep learning methods unlikely to succeed. Additionally, there are several issues with the existing structure and sequence data, as they are often of insufficient quality, highly biased and missing key information. Here, we discuss these challenges in detail and suggest some steps to remedy the situation. We believe that it is possible to create an accurate RNA structure prediction method, but it will require solving several data quality and volume issues, usage of data beyond simple sequence alignments, or the development of new less data-hungry machine learning methods.

Photoswitchable molecular glue for RNA: reversible photocontrol of structure and function of the ribozyme

Single-stranded RNA folds into a variety of secondary and higher-order structures. Distributions and dynamics of multiple RNA conformations are responsible for the biological function of RNA. We here developed a photoswitchable molecular glue for RNA, which could reversibly control the association of two unpaired RNA regions in response to light stimuli. The photoswitchable molecular glue, NCTA, is an RNA-binding ligand possessing a photoisomerizable azobenzene moiety. Z-NCTA is an active ligand for the target RNA containing 5'-WGG-3'/5'-WGG-3' (W = U or A) site and stabilizes its hybridized state, while its isomer E-NCTA is not. Photoreversible isomerization of NCTA enabled control of the secondary and tertiary structure of the target RNA. The RNA-cleaving activity of hammerhead ribozyme, where appropriate RNA folding is necessary, could be reversibly regulated by photoirradiation in cells treated with NCTA, demonstrating precise photocontrol of RNA structure and function by the photoswitchable molecular glue.

Emerging functional principles of tRNA-derived small RNAs and other regulatory small RNAs

Recent advancements in small RNA sequencing have unveiled a previously hidden world of regulatory small noncoding RNAs (sncRNAs) that extend beyond the well-studied small interfering RNAs, microRNAs, and piwi-interacting RNAs. This exploration, starting with tRNA-derived small RNAs, has led to the discovery of a diverse universe of sncRNAs derived from various longer structured RNAs such as rRNAs, small nucleolar RNAs, small nuclear RNAs, Y RNAs, and vault RNAs, with exciting uncharted functional possibilities. In this perspective, we discuss the emerging functional principles of sncRNAs beyond the well-known RNAi-like mechanisms, focusing on those that operate independent of linear sequence complementarity but rather function in an aptamer-like fashion. Aptamers use 3D structure for specific interactions with ligands and are modulated by RNA modifications and subcellular environments. Given that aptamer-like sncRNA functions are widespread and present in species lacking RNAi, they may represent an ancient functional principle that predates RNAi. We propose a rethinking of the origin of RNAi and its relationship with these aptamer-like functions in sncRNAs and how these complementary mechanisms shape biological processes. Lastly, the aptamer-like function of sncRNAs highlights the need for caution in using small RNA mimics in research and therapeutics, as their specificity is not restricted solely to linear sequence.

RNA folding studies inside peptide-rich droplets reveal roles of modified nucleosides at the origin of life

Compartmentalization of RNA in biopolymer-rich membraneless organelles is now understood to be pervasive and critical for the function of extant biology and has been proposed as a prebiotically plausible way to accumulate RNA. However, compartment-RNA interactions that drive encapsulation have the potential to influence RNA structure and function in compartment- and RNA sequence-dependent ways. Here, we detail next-generation sequencing (NGS) experiments performed in membraneless compartments called complex coacervates to characterize the fold of many different transfer RNAs (tRNAs) simultaneously under the potentially denaturing conditions of these compartments. Notably, we find that natural modifications favor the native fold of tRNAs in these compartments. This suggests that covalent RNA modifications could have played a critical role in metabolic processes at the origin of life.

Pre-mRNA splicing and its cotranscriptional connections

Transcription of eukaryotic genes by RNA polymerase II (Pol II) yields RNA precursors containing introns that must be spliced out and the flanking exons ligated together. Splicing is catalyzed by a dynamic ribonucleoprotein complex called the spliceosome. Recent evidence has shown that a large fraction of splicing occurs cotranscriptionally as the RNA chain is extruded from Pol II at speeds of up to 5 kb/minute. Splicing is more efficient when it is tethered to the transcription elongation complex, and this linkage permits functional coupling of splicing with transcription. We discuss recent progress that has uncovered a network of connections that link splicing to transcript elongation and other cotranscriptional RNA processing events.

Small Molecules Targeting Viral RNA

The majority of antivirals available target viral proteins; however, RNA is emerging as a new and promising antiviral target due to the presence of highly structured RNA in viral genomes fundamental for their replication cycle. Here, we discuss methods for the identification of RNA-targeting compounds, starting from the determination of RNA structures either from purified RNA or in living cells, followed by in silico screening on RNA and phenotypic assays to evaluate viral inhibition. Moreover, we review the small molecules known to target the programmed ribosomal frameshifting element of SARS-CoV-2, the internal ribosomal entry site of different viruses, and RNA elements of HIV.

Analysis of lncRNAs in Lupinus mutabilis (Tarwi) and Their Potential Role in Drought Response

Lupinus mutabilis is a legume with high agronomic potential and available transcriptomic data for which lncRNAs have not been studied. Therefore, our objective was to identify, characterize, and validate the drought-responsive lncRNAs in L. mutabilis. To achieve this, we used a multilevel approach based on lncRNA prediction, annotation, subcellular location, thermodynamic characterization, structural conservation, and validation. Thus, 590 lncRNAs were identified by at least two algorithms of lncRNA identification. Annotation with the PLncDB database showed 571 lncRNAs unique to tarwi and 19 lncRNAs with homology in 28 botanical families including Solanaceae (19), Fabaceae (17), Brassicaceae (17), Rutaceae (17), Rosaceae (16), and Malvaceae (16), among others. In total, 12 lncRNAs had homology in more than 40 species. A total of 67% of lncRNAs were located in the cytoplasm and 33% in exosomes. Thermodynamic characterization of S03 showed a stable secondary structure with -105.67 kcal/mol. This structure included three regions, with a multibranch loop containing a hairpin with a SECIS-like element. Evaluation of the structural conservation by CROSSalign revealed partial similarities between L. mutabilis (S03) and S. lycopersicum (Solyc04r022210.1). RT-PCR validation demonstrated that S03 was upregulated in a drought-tolerant accession of L. mutabilis. Finally, these results highlighted the importance of lncRNAs in tarwi improvement under drought conditions.

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.

Emerging Quantitative Biochemical, Structural, and Biophysical Methods for Studying Ribosome and Protein-RNA Complex Assembly

Ribosome assembly is one of the most fundamental processes of gene expression and has served as a playground for investigating the molecular mechanisms of how protein-RNA complexes (RNPs) assemble. A bacterial ribosome is composed of around 50 ribosomal proteins, several of which are co-transcriptionally assembled on a ~4500-nucleotide-long pre-rRNA transcript that is further processed and modified during transcription, the entire process taking around 2 min in vivo and being assisted by dozens of assembly factors. How this complex molecular process works so efficiently to produce an active ribosome has been investigated over decades, resulting in the development of a plethora of novel approaches that can also be used to study the assembly of other RNPs in prokaryotes and eukaryotes. Here, we review biochemical, structural, and biophysical methods that have been developed and integrated to provide a detailed and quantitative understanding of the complex and intricate molecular process of bacterial ribosome assembly. We also discuss emerging, cutting-edge approaches that could be used in the future to study how transcription, rRNA processing, cellular factors, and the native cellular environment shape ribosome assembly and RNP assembly at large.

Sequencing-based analysis of RNA structures in living cells with 2A3 via SHAPE-MaP

Chemical probing of RNA 2'-hydroxyl groups by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) is a rapid and powerful approach for querying RNA structures in living cells. At reverse transcription, sites of chemical modification can be encoded as mutations in the cDNA, a process called mutational profiling (MaP), enabling their detection via high-throughput sequencing. This chapter describes how to synthesize the SHAPE probe 2-aminopyridine-3-carboxylic acid imidazolide (2A3), how to use it to probe RNA structures in living bacteria, and how to generate Illumina-compatible SHAPE-MaP sequencing libraries. The protocol further describes data analysis using the RNA Framework, from raw sequencing data processing to experimentally-driven RNA secondary structure model generation.