The mechanisms of RNA SHAPE chemistry

Two riboswitch classes that share a common ligand-binding fold show major differences in the ability to accommodate mutations

Riboswitches are structured RNAs that sense small molecules to control expression. Prequeuosine1 (preQ1)-sensing riboswitches comprise three classes (I, II and III) that adopt distinct folds. Despite this difference, class II and III riboswitches each use 10 identical nucleotides to bind the preQ1 metabolite. Previous class II studies showed high sensitivity to binding-pocket mutations, which reduced preQ1 affinity and impaired function. Here, we introduced four equivalent mutations into a class III riboswitch, which maintained remarkably tight preQ1 binding. Co-crystal structures of each class III mutant showed compensatory interactions that preserve the fold. Chemical modification analysis revealed localized RNA flexibility changes for each mutant, but molecular dynamics (MD) simulations suggested that each mutation was not overtly destabilizing. Although impaired, class III mutants retained tangible gene-regulatory activity in bacteria compared to equivalent preQ1-II variants; mutations in the preQ1-pocket floor were tolerated better than wall mutations. Principal component analysis of MD trajectories suggested that the most functionally deleterious wall mutation samples different motions compared to wildtype. Overall, the results reveal that formation of compensatory interactions depends on the context of mutations within the overall fold and that functionally deleterious mutations can alter long-range correlated motions that link the riboswitch binding pocket with distal gene-regulatory sequences.

RNA infrastructure profiling illuminates transcriptome structure in crowded spaces

RNAs fold into compact structures and undergo protein interactions in cells. These occluded environments can block reagents that probe the underlying RNAs. Probes that can analyze structure in crowded settings can shed light on RNA biology. Here, we employ 2'-OH-reactive probes that are small enough to access folded RNA structure underlying close molecular contacts within cells, providing considerably broader coverage for intracellular RNA structural analysis. The data are analyzed first with well-characterized human ribosomal RNAs and then applied transcriptome-wide to polyadenylated transcripts. The smallest probe acetylimidazole (AcIm) yields 80% greater structural coverage than larger conventional reagent NAIN3, providing enhanced structural information in hundreds of transcripts. The acetyl probe also provides superior signals for identifying m6A modification sites in transcripts, particularly in sites that are inaccessible to a standard probe. Our strategy enables profiling RNA infrastructure, enhancing analysis of transcriptome structure, modification, and intracellular interactions, especially in spatially crowded settings.

Functional Validation of SAM Riboswitch Element A from Listeria monocytogenes

SreA is one of seven candidate S-adenosyl methionine (SAM) class I riboswitches identified in Listeria monocytogenes, a saprophyte and opportunistic foodborne pathogen. SreA precedes genes encoding a methionine ATP-binding cassette (ABC) transporter, which imports methionine and is presumed to regulate transcription of its downstream genes in a SAM-dependent manner. The proposed role of SreA in controlling the transcription of genes encoding an ABC transporter complex may have important implications for how the bacteria senses and responds to the availability of the metabolite SAM in the diverse environments in which L. monocytogenes persists. Here we validate SreA as a functional SAM-I riboswitch through ligand binding studies, structure characterization, and transcription termination assays. We determined that SreA has both a structure and SAM binding properties similar to those of other well-characterized SAM-I riboswitches. Despite the apparent structural similarities to previously described SAM-I riboswitches, SreA induces transcription termination in response to comparatively lower (nanomolar) ligand concentrations. Furthermore, SreA is a leaky riboswitch that permits some transcription of the downstream gene even in the presence of millimolar SAM, suggesting that L. monocytogenes may "dampen" the expression of genes for methionine import but likely does not turn them "OFF".

On the aqueous origins of the condensation polymers of life

Water is essential for life as we know it, but it has paradoxically been considered inimical to the emergence of life. Proteins and nucleic acids have sustained evolution and life for billions of years, but both are condensation polymers, suggesting that their formation requires the elimination of water. This presents intrinsic challenges at the origins of life, including how condensation polymer synthesis can overcome the thermodynamic pressure of hydrolysis in water and how nucleophiles can kinetically outcompete water to yield condensation products. The answers to these questions lie in balancing thermodynamic activation and kinetic stability. For peptides, an effective strategy is to directly harness the energy trapped in prebiotic molecules, such as nitriles, and avoid the formation of fully hydrolysed monomers. In this Review, we discuss how chemical energy can be built into precursors, retained, and released selectively for polymer synthesis. Looking to the future, the outstanding goals include how nucleic acids can be synthesized, avoiding the formation of fully hydrolysed monomers and what caused information to flow from nucleic acids to proteins.

Chemical diversity of reagents that modify RNA 2'-OH in water: a review

Electrophilic water-soluble compounds have proven versatile in reacting selectively with 2'-OH groups in RNA, enabling structure mapping, probing, caging, labeling, crosslinking, and conjugation of RNAs in vitro and in living cells. While early work focused on one or two types of reagents with limited properties, recent studies have greatly diversified the structure, properties, and applications of these reagents. Here we review the scope of documented RNA hydroxyl-reactive species reported to date, with an eye to the effects of chemical structure on reactivity with RNA and other useful properties. Multiple forms of carbonyl electrophiles are now known to react at the 2'-OH, and recently, sulfonyl and aryl electrophiles have also been documented to form bonds there in high yields as well. In addition to electrophilicity, data also point to significant effects of reagent stability, steric bulk, and chirality on reaction yields and selectivity. Finally, we outline reagent properties and principles that define utility in applications with RNA, with an eye to the design of future reagents.

Chemical and Enzyme-Mediated Chemical Reactions for Studying Nucleic Acids and Their Modifications

Nucleic acids are genetic information-carrying molecules inside cells. Apart from basic nucleotide building blocks, there exist various naturally occurring chemical modifications on nucleobase and ribose moieties, which greatly increase the encoding complexity of nuclei acids, contribute to the alteration of nucleic acid structures, and play versatile regulation roles in gene expression. To study the functions of certain nucleic acids in various biological contexts, robust tools to specifically label and identify these macromolecules and their modifications, and to illuminate their structures are highly necessary. In this review, we summarize recent technique advances of using chemical and enzyme-mediated chemical reactions to study nucleic acids and their modifications and structures. By highlighting the chemical principles of these techniques, we aim to present a perspective on the advancement of the field as well as to offer insights into developing specific chemical reactions and precise enzyme catalysis utilized for nucleic acids and their modifications.

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.

Second-Generation Chiral Amino Acid Derivatives Afford High Stereoselectivity and Stability in Aqueous RNA Acylation

Activated acyl species have proven versatile in the esterification of 2'-OH groups in RNA, enabling structure mapping, caging, profiling, and labeling of the biopolymer. Nearly all reagents developed for this reaction have been achiral; however, a recent study reported that simple chiral amino acid acylimidazole derivatives could yield diastereoselective reactions at RNA 2'-OH in water, enabling up to 4:1 selectivity in screening. Here, we investigated the effect of steric bulk on the stereoselectivity of RNA reaction and on the stability of adducts with a library of 36 chiral acylimidazole scaffolds with increasing steric demand. The results document the highest stereoselectivity yet achieved in RNA acylation reactions, with as high as >99:1 diastereoselectivity at >70% conversion. Also notably, the bulky adducts were found to have markedly improved stability on RNA.

Dynamics of Protein-RNA Interfaces Using All-Atom Molecular Dynamics Simulations

Facing the current challenges posed by human health diseases requires the understanding of cell machinery at a molecular level. The interplay between proteins and RNA is key for any physiological phenomenon, as well protein-RNA interactions. To understand these interactions, many experimental techniques have been developed, spanning a very wide range of spatial and temporal resolutions. In particular, the knowledge of tridimensional structures of protein-RNA complexes provides structural, mechanical, and dynamical pieces of information essential to understand their functions. To get insights into the dynamics of protein-RNA complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on nine different protein-RNA complexes with different functions and interface size by taking into account the bound and unbound forms. First, we characterized structural changes upon binding and, for the RNA part, the change in the puckering. Second, we extensively analyzed the interfaces, their dynamics and structural properties, and the structural waters involved in the binding, as well as the contacts mediated by them. Based on our analysis, the interfaces rearranged during the simulation time showing alternative and stable residue-residue contacts with respect to the experimental structure.

Bioorthogonal masked acylating agents for proximity-dependent RNA labelling

RNA localization is highly regulated, with subcellular organization driving context-dependent cell physiology. Although proximity-based labelling technologies that use highly reactive radicals or carbenes provide a powerful method for unbiased mapping of protein organization within a cell, methods for unbiased RNA mapping are scarce and comparably less robust. Here we develop α-alkoxy thioenol and chloroenol esters that function as potent acylating agents upon controlled ester unmasking. We pair these probes with subcellular-localized expression of a bioorthogonal esterase to establish a platform for spatial analysis of RNA: bioorthogonal acylating agents for proximity labelling and sequencing (BAP-seq). We demonstrate that, by selectively unmasking the enol probe in a locale of interest, we can map RNA distribution in membrane-bound and membrane-less organelles. The controlled-release acylating agent chemistry and corresponding BAP-seq method expand the scope of proximity labelling technologies and provide a powerful approach to interrogate the cellular organization of RNAs.

Analysis of Telomerase RNA Structure in Physcomitrium patens Indicates Functionally Relevant Transitions Between OPEN and CLOSED Conformations

Telomerase RNA (TR) conformation determines its function as a template for telomere synthesis and as a scaffold for the assembly of the telomerase nucleoprotein complex. Experimental analyses of TR secondary structure using DMS-Map Seq and SHAPE-Map Seq techniques show its CLOSED conformation as the consensus structure where the template region cannot perform its function. Our data show that the apparent discrepancy between experimental results and predicted TR functional conformation, mostly ignored in published studies, can be explained using data analysis based on single-molecule structure prediction from individual sequencing reads by the recently established DaVinci method. This method results in several clusters of secondary structures reflecting the structural dynamics of TR, possibly related to its multiple functional states. Interestingly, the presumed active (OPEN) conformation of TR corresponds to a minor fraction of TR under in vivo conditions. Therefore, structural polymorphism and dynamic TR transitions between CLOSED and OPEN conformations may be involved in telomerase activity regulation as a switch that functions independently of total TR transcript levels.

2'-OH as a universal handle for studying intracellular RNAs

RNA plays pivotal roles in most cellular processes, serving as both the traditional carrier of genetic information and as a key regulator of cellular functions. The advent of chemical technologies has contributed critically to the analysis of cellular RNA structures, functions, and interactions. Many of these methods and molecules involve the utilization of chemically reactive handles in RNAs, either introduced externally or inherent within the polymer itself. Among these handles, the 2'-hydroxyl (2'-OH) group has emerged as an exceptionally well-suited and general chemical moiety for the modification and profiling of RNAs in intracellular studies. In this review, we provide an overview of the recent advancements in intracellular applications of acylation at the 2'-OH group of RNA. We outline progress made in probing RNA structure and interactomes, controlling RNA function, RNA imaging, and analyzing RNA-small molecule interactions, all achieved in living cells through this simple chemical handle on the biopolymer.

RNA Infrastructure Profiling Illuminates Transcriptome Structure in Crowded Spaces

RNAs can fold into compact three-dimensional structures, and most RNAs undergo protein interactions in the cell. These compact and occluded environments can block the ability of structure-probing agents to provide useful data about the folding and modification of the underlying RNA. The development of probes that can analyze structure in crowded settings, and differentiate the proximity of interactions, can shed new light on RNA biology. To this end, here we employ 2'-OH-reactive probes that are small enough to access folded RNA structure underlying many close molecular contacts within cells, providing considerably broader coverage for intracellular RNA structural analysis. We compare reverse transcriptase stops in RNA-Seq data from probes of small and standard size to assess RNA-protein proximity and evaluate solvent-exposed tunnels adjacent to RNA. The data are analyzed first with structurally characterized complexes (human 18S and 28S RNA), and then applied transcriptome-wide to polyadenylated transcripts in HEK293 cells. In our transcriptome profile, the smallest probe acetylimidazole (AcIm) yields 80% greater structural coverage than larger conventional reagent NAIN3, providing enhanced structural information in hundreds of transcripts. We further show that acetyl probes provide superior signals for identifying m6A modification sites in transcripts, and provide information regarding methylation sites that are inaccessible to a larger standard probe. RNA infrastructure profiling (RISP) enables enhanced analysis of transcriptome structure, modification, and interactions in living cells, especially in spatially crowded settings.

Mutation signature filtering enables high-fidelity RNA structure probing at all four nucleobases with DMS

Chemical probing experiments have transformed RNA structure analysis, enabling high-throughput measurement of base-pairing in living cells. Dimethyl sulfate (DMS) is one of the most widely used structure probing reagents and has played a pivotal role in enabling next-generation single-molecule probing analyses. However, DMS has traditionally only been able to probe adenine and cytosine nucleobases. We previously showed that, using appropriate conditions, DMS can also be used to interrogate base-pairing of uracil and guanines in vitro at reduced accuracy. However, DMS remained unable to informatively probe guanines in cells. Here, we develop an improved DMS mutational profiling (MaP) strategy that leverages the unique mutational signature of N1-methylguanine DMS modifications to enable high-fidelity structure probing at all four nucleotides, including in cells. Using information theory, we show that four-base DMS reactivities convey greater structural information than current two-base DMS and SHAPE probing strategies. Four-base DMS experiments further enable improved direct base-pair detection by single-molecule PAIR analysis, and ultimately support RNA structure modeling at superior accuracy. Four-base DMS probing experiments are straightforward to perform and will broadly facilitate improved RNA structural analysis in living cells.

Mutation signature filtering enables high-fidelity RNA structure probing at all four nucleobases with DMS

Chemical probing experiments have transformed RNA structure analysis, enabling high-throughput measurement of base-pairing in living cells. Dimethyl sulfate (DMS) is one of the most widely used structure probing reagents and has played a prominent role in enabling next-generation single-molecule probing analyses. However, DMS has traditionally only been able to probe adenine and cytosine nucleobases. We previously showed that, using appropriate conditions, DMS can also be used to interrogate base-pairing of uracil and guanines in vitro at reduced accuracy. However, DMS remained unable to informatively probe guanines in cells. Here, we develop an improved DMS mutational profiling (MaP) strategy that leverages the unique mutational signature of N 1 -methylguanine DMS modifications to enable robust, high-fidelity structure probing at all four nucleotides, including in cells. Using information theory, we show that four-base DMS reactivities convey greater structural information than comparable two-base DMS and SHAPE probing strategies. Four-base DMS experiments further enable improved direct base-pair detection by single-molecule PAIR analysis, and ultimately support RNA structure modeling at superior accuracy. Four-base DMS probing experiments are easily performed and will broadly facilitate improved RNA structural analysis in living cells.

A riboswitch separated from its ribosome-binding site still regulates translation

Riboswitches regulate downstream gene expression by binding cellular metabolites. Regulation of translation initiation by riboswitches is posited to occur by metabolite-mediated sequestration of the Shine-Dalgarno sequence (SDS), causing bypass by the ribosome. Recently, we solved a co-crystal structure of a prequeuosine1-sensing riboswitch from Carnobacterium antarcticum that binds two metabolites in a single pocket. The structure revealed that the second nucleotide within the gene-regulatory SDS, G34, engages in a crystal contact, obscuring the molecular basis of gene regulation. Here, we report a co-crystal structure wherein C10 pairs with G34. However, molecular dynamics simulations reveal quick dissolution of the pair, which fails to reform. Functional and chemical probing assays inside live bacterial cells corroborate the dispensability of the C10-G34 pair in gene regulation, leading to the hypothesis that the compact pseudoknot fold is sufficient for translation attenuation. Remarkably, the C. antarcticum aptamer retained significant gene-regulatory activity when uncoupled from the SDS using unstructured spacers up to 10 nucleotides away from the riboswitch-akin to steric-blocking employed by sRNAs. Accordingly, our work reveals that the RNA fold regulates translation without SDS sequestration, expanding known riboswitch-mediated gene-regulatory mechanisms. The results infer that riboswitches exist wherein the SDS is not embedded inside a stable fold.

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.

Diverse Reagent Scaffolds Provide Differential Selectivity of 2'-OH Acylation in RNA

RNA 2'-OH acylation is widely used both for mapping structure and for conjugating RNA, generally relying on selective reactions with unpaired nucleotides over paired ones. Common reagents for this acylation have been chiefly restricted to two similar aryl scaffolds, leaving open the question of how more broadly varied reagent structure might affect selectivity. Here, we prepared a set of 10 structurally diverse acylimidazole reagents and employed deep sequencing to profile their reactivity and selectivity in an RNA library of systematically varied structure. We show that structure-directed reactivity profiles vary significantly with the reagent scaffold, and we document new acylating agents that have altered selectivity profiles, including reagents that show elevated selectivity within loops, as well as compounds with reduced off-target reactivity in loop closing base pairs. Interestingly, we also show that the simplest reagent (acetylimidazole) is cell permeable and is small enough to map RNA structure in the presence of protein contacts that block other reagents. Finally, we describe reagents that show elevated selectivity within small loops, with applications in site-selective labeling. The results provide new tools for improved conjugation and mapping of RNA.

Chemical Probing of RNA Structure In Vivo Using SHAPE-MaP and DMS-MaP

The functional roles of RNAs are often regulated by their structure. Selective 2' hydroxyl acylation analyzed by primer extension (SHAPE) and dimethyl sulfate (DMS) base reactivity can be employed to investigate the flexibility of nucleotides and correlate it to the constraints imparted by base-pairing and/or protein-binding. In vivo, a multitude of proteins could bind an RNA molecule, regulating its structure and function. Hence, to obtain a more comprehensive view of the RNA structure-function relationship in vivo, it may be required to characterize both the RNA structure and the RNA-protein interaction network. In this chapter, we describe methods for characterizing the in vivo nucleotide flexibility of RNA in cells by SHAPE-MaP (SHAPE by Mutational Profiling) and DMS-MaP. In another chapter, we will discuss the characterization of RNA-protein interaction network by RNP-MaP.

Acylation probing of "generic" RNA libraries reveals critical influence of loop constraints on reactivity

The reactivity of RNA 2'-OH acylation is broadly useful both in probing structure and in preparing conjugates. To date, this reactivity has been analyzed in limited sets of biological RNA sequences, leaving open questions of how reactivity varies inherently without regard to sequence in structured contexts. We constructed and probed "generic" structured RNA libraries using homogeneous loop sequences, employing deep sequencing to carry out a systematic survey of reactivity. We find a wide range of RNA reactivities among single-stranded sequences, with nearest neighbors playing substantial roles. Remarkably, certain small loops are found to be far more reactive on average (up to 4,000-fold) than single-stranded RNAs, due to conformational constraints that enhance reactivity. Among loops, we observe large variations in reactivity based on size, type, and position. The results lend insights into RNA designs for achieving high-efficiency local conjugation and provide new opportunities to refine structure analysis.

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

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

Lessons Learned and Yet-to-Be Learned on the Importance of RNA Structure in SARS-CoV-2 Replication

SARS-CoV-2, the etiological agent responsible for the COVID-19 pandemic, is a member of the virus family Coronaviridae, known for relatively extensive (~30-kb) RNA genomes that not only encode for numerous proteins but are also capable of forming elaborate structures. As highlighted in this review, these structures perform critical functions in various steps of the viral life cycle, ultimately impacting pathogenesis and transmissibility. We examine these elements in the context of coronavirus evolutionary history and future directions for curbing the spread of SARS-CoV-2 and other potential human coronaviruses. While we focus on structures supported by a variety of biochemical, biophysical, and/or computational methods, we also touch here on recent evidence for novel structures in both protein-coding and noncoding regions of the genome, including an assessment of the potential role for RNA structure in the controversial finding of SARS-CoV-2 integration in “long COVID” patients. This review aims to serve as a consolidation of previous works on coronavirus and more recent investigation of SARS-CoV-2, emphasizing the need for improved understanding of the role of RNA structure in the evolution and adaptation of these human viruses.

ShapeSorter: a fully probabilistic method for detecting conserved RNA structure features supported by SHAPE evidence

There is an increased interest in the determination of RNA structures in vivo as it is now possible to probe them in a high-throughput manner, e.g. using SHAPE protocols. By now, there exist a range of computational methods that integrate experimental SHAPE-probing evidence into computational RNA secondary structure prediction. The state-of-the-art in this field is currently provided by computational methods that employ the minimum-free energy strategy for prediction RNA secondary structures with SHAPE-probing evidence. These methods, however, rely on the assumption that transcripts in vivo fold into the thermodynamically most stable configuration and ignore evolutionary evidence for conserved RNA structure features. We here present a new computational method, ShapeSorter, that predicts RNA structure features without employing the thermodynamic strategy. Instead, ShapeSorter employs a fully probabilistic framework to identify RNA structure features that are supported by evolutionary and SHAPE-probing evidence. Our method can capture RNA structure heterogeneity, pseudo-knotted RNA structures as well as transient and mutually exclusive RNA structure features. Moreover, it estimates P-values for the predicted RNA structure features which allows for easy filtering and ranking. We investigate the merits of our method in a comprehensive performance benchmarking and conclude that ShapeSorter has a significantly superior performance for predicting base-pairs than the existing state-of-the-art methods.

Insights into the secondary and tertiary structure of the Bovine Viral Diarrhea Virus Internal Ribosome Entry Site

The internal ribosome entry site (IRES) RNA of bovine viral diarrhoea virus (BVDV), an economically significant Pestivirus, is required for the cap-independent translation of viral genomic RNA. Thus, it is essential for viral replication and pathogenesis. We applied a combination of high-throughput biochemical RNA structure probing (SHAPE-MaP) and in silico modelling approaches to gain insight into the secondary and tertiary structures of BVDV IRES RNA. Our study demonstrated that BVDV IRES RNA in solution forms a modular architecture composed of three distinct structural domains (I-III). Two regions within domain III are represented in tertiary interactions to form an H-type pseudoknot. Computational modelling of the pseudoknot motif provided a fine-grained picture of the tertiary structure and local arrangement of helices in the BVDV IRES. Furthermore, comparative genomics and consensus structure predictions revealed that the pseudoknot is evolutionarily conserved among many Pestivirus species. These studies provide detailed insight into the structural arrangement of BVDV IRES RNA H-type pseudoknot and encompassing motifs that likely contribute to the optimal functionality of viral cap-independent translation element.

Conjugation of RNA via 2'-OH acylation: Mechanisms determining nucleotide reactivity

The acylation reactivity of RNA 2'-OH groups has proven broadly useful for labeling and mapping RNA. Here we perform kinetics studies to test the mechanisms governing this reaction, and we find strong steric and inductive effects modulating reactivity. The results shed light on new strategies for improved conjugation and mapping.

Structural characterization of a new subclass of panicum mosaic virus-like 3' cap-independent translation enhancer

Canonical eukaryotic mRNA translation requires 5'cap recognition by initiation factor 4E (eIF4E). In contrast, many positive-strand RNA virus genomes lack a 5'cap and promote translation by non-canonical mechanisms. Among plant viruses, PTEs are a major class of cap-independent translation enhancers located in/near the 3'UTR that recruit eIF4E to greatly enhance viral translation. Previous work proposed a single form of PTE characterized by a Y-shaped secondary structure with two terminal stem-loops (SL1 and SL2) atop a supporting stem containing a large, G-rich asymmetric loop that forms an essential pseudoknot (PK) involving C/U residues located between SL1 and SL2. We found that PTEs with less than three consecutive cytidylates available for PK formation have an upstream stem-loop that forms a kissing loop interaction with the apical loop of SL2, important for formation/stabilization of PK. PKs found in both subclasses of PTE assume a specific conformation with a hyperreactive guanylate (G*) in SHAPE structure probing, previously found critical for binding eIF4E. While PTE PKs were proposed to be formed by Watson-Crick base-pairing, alternative chemical probing and 3D modeling indicate that the Watson-Crick faces of G* and an adjacent guanylate have high solvent accessibilities. Thus, PTE PKs are likely composed primarily of non-canonical interactions.

In Vitro Methods to Decipher the Structure of Viral RNA Genomes

RNA viruses encode essential information in their genomes as conserved structural elements that are involved in efficient viral protein synthesis, replication, and encapsidation. These elements can also establish complex networks of RNA-RNA interactions, the so-called RNA interactome, to shape the viral genome and control different events during intracellular infection. In recent years, targeting these conserved structural elements has become a promising strategy for the development of new antiviral tools due to their sequence and structural conservation. In this context, RNA-based specific therapeutic strategies, such as the use of siRNAs have been extensively pursued to target the genome of different viruses. Importantly, siRNA-mediated targeting is not a straightforward approach and its efficiency is highly dependent on the structure of the target region. Therefore, the knowledge of the viral structure is critical for the identification of potentially good target sites. Here, we describe detailed protocols used in our laboratory for the in vitro study of the structure of viral RNA genomes. These protocols include DMS (dimethylsulfate) probing, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) analysis, and HMX (2'-hydroxyl molecular interference). These methodologies involve the use of high-throughput analysis techniques that provide extensive information about the 3D folding of the RNA under study and the structural tuning derived from the interactome activity. They are therefore a good tool for the development of new RNA-based antiviral compounds.

Progress toward SHAPE Constrained Computational Prediction of Tertiary Interactions in RNA Structure

As more sequencing data accumulate and novel puzzling genetic regulations are discovered, the need for accurate automated modeling of RNA structure increases. RNA structure modeling from chemical probing experiments has made tremendous progress, however accurately predicting large RNA structures is still challenging for several reasons: RNA are inherently flexible and often adopt many energetically similar structures, which are not reliably distinguished by the available, incomplete thermodynamic model. Moreover, computationally, the problem is aggravated by the relevance of pseudoknots and non-canonical base pairs, which are hardly predicted efficiently. To identify nucleotides involved in pseudoknots and non-canonical interactions, we scrutinized the SHAPE reactivity of each nucleotide of the 188 nt long lariat-capping ribozyme under multiple conditions. Reactivities analyzed in the light of the X-ray structure were shown to report accurately the nucleotide status. Those that seemed paradoxical were rationalized by the nucleotide behavior along molecular dynamic simulations. We show that valuable information on intricate interactions can be deduced from probing with different reagents, and in the presence or absence of Mg2+. Furthermore, probing at increasing temperature was remarkably efficient at pointing to non-canonical interactions and pseudoknot pairings. The possibilities of following such strategies to inform structure modeling software are discussed.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

SHAPER: A Web Server for Fast and Accurate SHAPE Reactivity Prediction

Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) chemical probing serves as a convenient and efficient experiment technique for providing information about RNA local flexibility. The local structural information contained in SHAPE reactivity data can be used as constraints in 2D/3D structure predictions. Here, we present SHAPE predictoR (SHAPER), a web server for fast and accurate SHAPE reactivity prediction. The main purpose of the SHAPER web server is to provide a portal that uses experimental SHAPE data to refine 2D/3D RNA structure selection. Input structures for the SHAPER server can be obtained through experimental or computational modeling. The SHAPER server can accept RNA structures with single or multiple conformations, and the predicted SHAPE profile and correlation with experimental SHAPE data (if provided) for each conformation can be freely downloaded through the web portal. The SHAPER web server is available at http://rna.physics.missouri.edu/shaper/.

Direct Mapping of Higher-Order RNA Interactions by SHAPE-JuMP

Higher-order structure governs function for many RNAs. However, discerning this structure for large RNA molecules in solution is an unresolved challenge. Here, we present SHAPE-JuMP (selective 2'-hydroxyl acylation analyzed by primer extension and juxtaposed merged pairs) to interrogate through-space RNA tertiary interactions. A bifunctional small molecule is used to chemically link proximal nucleotides in an RNA structure. The RNA cross-link site is then encoded into complementary DNA (cDNA) in a single, direct step using an engineered reverse transcriptase that "jumps" across cross-linked nucleotides. The resulting cDNAs contain a deletion relative to the native RNA sequence, which can be detected by sequencing, that indicates the sites of cross-linked nucleotides. SHAPE-JuMP measures RNA tertiary structure proximity concisely across large RNA molecules at nanometer resolution. SHAPE-JuMP is especially effective at measuring interactions in multihelix junctions and loop-to-helix packing, enables modeling of the global fold for RNAs up to several hundred nucleotides in length, facilitates ranking of structural models by consistency with through-space restraints, and is poised to enable solution-phase structural interrogation and modeling of complex RNAs.

In Vivo RNA Chemical Footprinting Analysis in Archaea

RNA structural conformation and dynamics govern the functional properties of all RNA/RNP. Accordingly, defining changes of RNA structure and dynamics in various conditions may provide detailed insight into how RNA structural properties regulate the function of RNA/RNP. Traditional chemical footprinting analysis using chemical modifiers allows to sample the dynamics and conformation landscape of diverse RNA/RNP. However, many chemical modifiers are limited in their capacity to provide unbiased information reflecting the in vivo RNA/RNP structural landscape. In the recent years, the development of selective-2'-hydroxyl acylation analyzed by primer extension (SHAPE) methodology that uses powerful new chemical modifiers has significantly improved in vitro and in vivo structural probing of secondary and tertiary interactions of diverse RNA species at the single nucleotide level.Although the original discovery of Archaea as an independent domain of life is intimately linked to the technological development of RNA analysis, our understanding of in vivo RNA structural conformation and dynamics in this domain of life remains scarce.This protocol describes the in vivo use of SHAPE chemistry in two evolutionary divergent model Archaea, Sulfolobus acidocaldarius and Haloferax volcanii.

Mapping the RNA Chaperone Activity of the T. brucei Editosome Using SHAPE Chemical Probing

Mitochondrial pre-mRNAs in African trypanosomes adopt intricately folded, highly stable 2D and 3D structures. The RNA molecules are substrates of a U-nucleotide-specific insertion/deletion-type RNA editing reaction, which is catalyzed by a 0.8 MDa protein complex known as the editosome. RNA binding to the editosome is followed by a chaperone-mediated RNA remodeling reaction. The reaction increases the dynamic of specifically U-nucleotides to lower their base-pairing probability and as a consequence generates a simplified RNA folding landscape that is critical for the progression of the editing reaction cycle. Here we describe a chemical mapping method to quantitatively monitor the chaperone-driven structural changes of pre-edited mRNAs upon editosome binding. The method is known as selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE). SHAPE is based on the differential electrophilic modification of ribose 2'-hydroxyl groups in structurally constraint (double-stranded) versus structurally unconstrained (single-stranded) nucleotides. Electrophilic anhydrides such as 1-methyl-7-nitroisatoic anhydride are used as probing reagents, and the ribose 2'-modified nucleotides are mapped as abortive cDNA synthesis products. As a result, SHAPE allows the identification of all single-stranded and base-paired regions in a given RNA, and the data are used to compute experimentally derived RNA 2D structures. A side-by-side comparison of the RNA 2D folds in the pre- and post-chaperone states finally maps the chaperone-induced dynamic of the different pre-mRNAs with single-nucleotide resolution.

Characteristic chemical probing patterns of loop motifs improve prediction accuracy of RNA secondary structures

RNA structures play a fundamental role in nearly every aspect of cellular physiology and pathology. Gaining insights into the functions of RNA molecules requires accurate predictions of RNA secondary structures. However, the existing thermodynamic folding models remain less accurate than desired, even when chemical probing data, such as selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reactivities, are used as restraints. Unlike most SHAPE-directed algorithms that only consider SHAPE restraints for base pairing, we extract two-dimensional structural features encoded in SHAPE data and establish robust relationships between characteristic SHAPE patterns and loop motifs of various types (hairpin, internal, and bulge) and lengths (2-11 nucleotides). Such characteristic SHAPE patterns are closely related to the sugar pucker conformations of loop residues. Based on these patterns, we propose a computational method, SHAPELoop, which refines the predicted results of the existing methods, thereby further improving their prediction accuracy. In addition, SHAPELoop can provide information about local or global structural rearrangements (including pseudoknots) and help researchers to easily test their hypothesized secondary structures.

SHAPE Directed Discovery of New Functions in Large RNAs

RNA lies upstream of nearly all biology and functions as the central conduit of information exchange in all cells. RNA molecules encode information both in their primary sequences and in complex structures that form when an RNA folds back on itself. From the time of discovery of mRNA in the late 1950s until quite recently, we had only a rudimentary understanding of RNA structure across vast regions of most messenger and noncoding RNAs. This deficit is now rapidly being addressed, especially by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry, mutational profiling (MaP), and closely related platform technologies that, collectively, create chemical microscopes for RNA. These technologies make it possible to interrogate RNA structure, quantitatively, at nucleotide resolution, and at large scales, for entire mRNAs, noncoding RNAs, and viral RNA genomes. By applying comprehensive structure probing to diverse problems, we and others are showing that control of biological function mediated by RNA structure is ubiquitous across prokaryotic and eukaryotic organisms.Work over the past decade using SHAPE-based analyses has clarified key principles. First, the method of RNA structure probing matters. SHAPE-MaP, with its direct and one-step readout that probes nearly every nucleotide by reaction at the 2'-hydroxyl, gives a more detailed and accurate readout than alternatives. Second, comprehensive chemical probing is essential. Focusing on fragments of large RNAs or using meta-gene or statistical analyses to compensate for sparse data sets misses critical features and often yields structure models with poor predictive power. Finally, every RNA has its own internal structural personality. There are myriad ways in which RNA structure modulates sequence accessibility, protein binding, translation, splice-site choice, phase separation, and other fundamental biological processes. In essentially every instance where we have applied rigorous and quantitative SHAPE technologies to study RNA structure-function interrelationships, new insights regarding biological regulatory mechanisms have emerged. RNA elements with more complex higher-order structures appear more likely to contain high-information-content clefts and pockets that bind small molecules, broadly informing a vigorous field of RNA-targeted drug discovery.The broad implications of this collective work are twofold. First, it is long past time to abandon depiction of large RNAs as simple noodle-like or gently flowing molecules. Instead, we need to emphasize that nearly all RNAs are punctuated with distinctive internal structures, a subset of which modulate function in profound ways. Second, structure probing should be an integral component of any effort that seeks to understand the functional nexuses and biological roles of large RNAs.

A novel SHAPE reagent enables the analysis of RNA structure in living cells with unprecedented accuracy

Due to the mounting evidence that RNA structure plays a critical role in regulating almost any physiological as well as pathological process, being able to accurately define the folding of RNA molecules within living cells has become a crucial need. We introduce here 2-aminopyridine-3-carboxylic acid imidazolide (2A3), as a general probe for the interrogation of RNA structures in vivo. 2A3 shows moderate improvements with respect to the state-of-the-art selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reagent NAI on naked RNA under in vitro conditions, but it significantly outperforms NAI when probing RNA structure in vivo, particularly in bacteria, underlining its increased ability to permeate biological membranes. When used as a restraint to drive RNA structure prediction, data derived by SHAPE-MaP with 2A3 yields more accurate predictions than NAI-derived data. Due to its extreme efficiency and accuracy, we can anticipate that 2A3 will rapidly take over conventional SHAPE reagents for probing RNA structures both in vitro and in vivo.

SHAPE Profiling to Probe Group II Intron Conformational Dynamics During Splicing

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) is a widely used technique for studying the structure and function of RNA molecules. It characterizes the flexibility of single nucleotides in the context of the local RNA structure. Here we describe the application of SHAPE-MaP (mutational profiling) to study different conformational states of the group II intron during the self-splicing reaction.

Light-Dependent Translation Change of Arabidopsis psbA Correlates with RNA Structure Alterations at the Translation Initiation Region

mRNA secondary structure influences translation. Proteins that modulate the mRNA secondary structure around the translation initiation region may regulate translation in plastids. To test this hypothesis, we exposed Arabidopsis thaliana to high light, which induces translation of psbA mRNA encoding the D1 subunit of photosystem II. We assayed translation by ribosome profiling and applied two complementary methods to analyze in vivo RNA secondary structure: DMS-MaPseq and SHAPE-seq. We detected increased accessibility of the translation initiation region of psbA after high light treatment, likely contributing to the observed increase in translation by facilitating translation initiation. Furthermore, we identified the footprint of a putative regulatory protein in the 5' UTR of psbA at a position where occlusion of the nucleotide sequence would cause the structure of the translation initiation region to open up, thereby facilitating ribosome access. Moreover, we show that other plastid genes with weak Shine-Dalgarno sequences (SD) are likely to exhibit psbA-like regulation, while those with strong SDs do not. This supports the idea that changes in mRNA secondary structure might represent a general mechanism for translational regulation of psbA and other plastid genes.

Sieving RNA 3D Structures with SHAPE and Evaluating Mechanisms Driving Sequence-Dependent Reactivity Bias

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemical probing provides local RNA flexibility information at single-nucleotide resolution. In general, SHAPE is thought of as a secondary structure (2D) technology, but we find evidence that robust tertiary structure (3D) information is contained in SHAPE data. Here, we report a new model that achieves a higher correlation between SHAPE data and native RNA 3D structures than the previous 3D structure-SHAPE relationship model. Furthermore, we demonstrate that the new model improves our ability to discern between SHAPE-compatible and -incompatible structures on model decoys. After identifying sequence-dependent bias in SHAPE experiments, we propose a mechanism driving sequence-dependent bias in SHAPE experiments, using replica-exchange umbrella sampling simulations to confirm that the SHAPE sequence bias is largely explained by the stability of the unreacted SHAPE reagent in the binding pocket. Taken together, this work represents multiple practical advances in our mechanistic and predictive understanding of SHAPE technology.

Computationally reconstructing cotranscriptional RNA folding from experimental data reveals rearrangement of non-native folding intermediates

The series of RNA folding events that occur during transcription can critically influence cellular RNA function. Here, we present reconstructing RNA dynamics from data (R2D2), a method to uncover details of cotranscriptional RNA folding. We model the folding of the Escherichia coli signal recognition particle (SRP) RNA and show that it requires specific local structural fluctuations within a key hairpin to engender efficient cotranscriptional conformational rearrangement into the functional structure. All-atom molecular dynamics simulations suggest that this rearrangement proceeds through an internal toehold-mediated strand-displacement mechanism, which can be disrupted with a point mutation that limits local structural fluctuations and rescued with compensating mutations that restore these fluctuations. Moreover, a cotranscriptional folding intermediate could be cleaved in vitro by recombinant E. coli RNase P, suggesting potential cotranscriptional processing. These results from experiment-guided multi-scale modeling demonstrate that even an RNA with a simple functional structure can undergo complex folding and processing during synthesis.

RNAthor - fast, accurate normalization, visualization and statistical analysis of RNA probing data resolved by capillary electrophoresis

RNAs adopt specific structures to perform their functions, which are critical to fundamental cellular processes. For decades, these structures have been determined and modeled with strong support from computational methods. Still, the accuracy of the latter ones depends on the availability of experimental data, for example, chemical probing information that can define pseudo-energy constraints for RNA folding algorithms. At the same time, diverse computational tools have been developed to facilitate analysis and visualization of data from RNA structure probing experiments followed by capillary electrophoresis or next-generation sequencing. RNAthor, a new software tool for the fully automated normalization of SHAPE and DMS probing data resolved by capillary electrophoresis, has recently joined this collection. RNAthor automatically identifies unreliable probing data. It normalizes the reactivity information to a uniform scale and uses it in the RNA secondary structure prediction. Our web server also provides tools for fast and easy RNA probing data visualization and statistical analysis that facilitates the comparison of multiple data sets. RNAthor is freely available at http://rnathor.cs.put.poznan.pl/.

Allosteric regulation of U1 snRNP by splicing regulatory proteins controls spliceosomal assembly

Alternative splicing is responsible for much of the transcriptomic and proteomic diversity observed in eukaryotes and involves combinatorial regulation by many cis-acting elements and trans-acting factors. SR and hnRNP splicing regulatory proteins often have opposing effects on splicing efficiency depending on where they bind the pre-mRNA relative to the splice site. Position-dependent splicing repression occurs at spliceosomal E-complex, suggesting that U1 snRNP binds but cannot facilitate higher order spliceosomal assembly. To test the hypothesis that the structure of U1 snRNA changes during activation or repression, we developed a method to structure-probe native U1 snRNP in enriched conformations that mimic activated or repressed spliceosomal E-complexes. While the core of U1 snRNA is highly structured, the 5' end of U1 snRNA shows different SHAPE reactivities and psoralen crosslinking efficiencies depending on where splicing regulatory elements are located relative to the 5' splice site. A motif within the 5' splice site binding region of U1 snRNA is more reactive toward SHAPE electrophiles when repressors are bound, suggesting U1 snRNA is bound, but less base-paired. These observations demonstrate that splicing regulators modulate splice site selection allosterically.

IPANEMAP: integrative probing analysis of nucleic acids empowered by multiple accessibility profiles

The manual production of reliable RNA structure models from chemical probing experiments benefits from the integration of information derived from multiple protocols and reagents. However, the interpretation of multiple probing profiles remains a complex task, hindering the quality and reproducibility of modeling efforts. We introduce IPANEMAP, the first automated method for the modeling of RNA structure from multiple probing reactivity profiles. Input profiles can result from experiments based on diverse protocols, reagents, or collection of variants, and are jointly analyzed to predict the dominant conformations of an RNA. IPANEMAP combines sampling, clustering and multi-optimization, to produce secondary structure models that are both stable and well-supported by experimental evidences. The analysis of multiple reactivity profiles, both publicly available and produced in our study, demonstrates the good performances of IPANEMAP, even in a mono probing setting. It confirms the potential of integrating multiple sources of probing data, informing the design of informative probing assays.

Intact RNA structurome reveals mRNA structure-mediated regulation of miRNA cleavage in vivo

MicroRNA (miRNA)-mediated cleavage is involved in numerous essential cellular pathways. miRNAs recognize target RNAs via sequence complementarity. In addition to complementarity, in vitro and in silico studies have suggested that RNA structure may influence the accessibility of mRNAs to miRNA-induced silencing complexes (miRISCs), thereby affecting RNA silencing. However, the regulatory mechanism of mRNA structure in miRNA cleavage remains elusive. We investigated the role of in vivo RNA secondary structure in miRNA cleavage by developing the new CAP-STRUCTURE-seq method to capture the intact mRNA structurome in Arabidopsis thaliana. This approach revealed that miRNA target sites were not structurally accessible for miRISC binding prior to cleavage in vivo. Instead, we found that the unfolding of the target site structure plays a key role in miRISC activity in vivo. We found that the single-strandedness of the two nucleotides immediately downstream of the target site, named Target Adjacent nucleotide Motif, can promote miRNA cleavage but not miRNA binding, thus decoupling target site binding from cleavage. Our findings demonstrate that mRNA structure in vivo can modulate miRNA cleavage, providing evidence of mRNA structure-dependent regulation of biological processes.

Analysis of a preQ1-I riboswitch in effector-free and bound states reveals a metabolite-programmed nucleobase-stacking spine that controls gene regulation

Riboswitches are structured RNA motifs that recognize metabolites to alter the conformations of downstream sequences, leading to gene regulation. To investigate this molecular framework, we determined crystal structures of a preQ1-I riboswitch in effector-free and bound states at 2.00 Å and 2.65 Å-resolution. Both pseudoknots exhibited the elusive L2 loop, which displayed distinct conformations. Conversely, the Shine-Dalgarno sequence (SDS) in the S2 helix of each structure remained unbroken. The expectation that the effector-free state should expose the SDS prompted us to conduct solution experiments to delineate environmental changes to specific nucleobases in response to preQ1. We then used nudged elastic band computational methods to derive conformational-change pathways linking the crystallographically-determined effector-free and bound-state structures. Pathways featured: (i) unstacking and unpairing of L2 and S2 nucleobases without preQ1-exposing the SDS for translation and (ii) stacking and pairing L2 and S2 nucleobases with preQ1-sequestering the SDS. Our results reveal how preQ1 binding reorganizes L2 into a nucleobase-stacking spine that sequesters the SDS, linking effector recognition to biological function. The generality of stacking spines as conduits for effector-dependent, interdomain communication is discussed in light of their existence in adenine riboswitches, as well as the turnip yellow mosaic virus ribosome sensor.

Challenges and perspectives for structural biology of lncRNAs-the example of the Xist lncRNA A-repeats

Following the discovery of numerous long non-coding RNA (lncRNA) transcripts in the human genome, their important roles in biology and human disease are emerging. Recent progress in experimental methods has enabled the identification of structural features of lncRNAs. However, determining high-resolution structures is challenging as lncRNAs are expected to be dynamic and adopt multiple conformations, which may be modulated by interaction with protein binding partners. The X-inactive specific transcript (Xist) is necessary for X inactivation during dosage compensation in female placental mammals and one of the best-studied lncRNAs. Recent progress has provided new insights into the domain organization, molecular features, and RNA binding proteins that interact with distinct regions of Xist. The A-repeats located at the 5′ end of the transcript are of particular interest as they are essential for mediating silencing of the inactive X chromosome. Here, we discuss recent progress with elucidating structural features of the Xist lncRNA, focusing on the A-repeats. We discuss the experimental and computational approaches employed that have led to distinct structural models, likely reflecting the intrinsic dynamics of this RNA. The presence of multiple dynamic conformations may also play an important role in the formation of the associated RNPs, thus influencing the molecular mechanism underlying the biological function of the Xist A-repeats. We propose that integrative approaches that combine biochemical experiments and high-resolution structural biology in vitro with chemical probing and functional studies in vivo are required to unravel the molecular mechanisms of lncRNAs.

Lead-seq: transcriptome-wide structure probing in vivo using lead(II) ions

The dynamic conformation of RNA molecules within living cells is key to their function. Recent advances in probing the RNA structurome in vivo, including the use of SHAPE (Selective 2'-Hydroxyl Acylation analyzed by Primer Extension) or kethoxal reagents or DMS (dimethyl sulfate), provided unprecedented insights into the architecture of RNA molecules in the living cell. Here, we report the establishment of lead probing in a global RNA structuromics approach. In order to elucidate the transcriptome-wide RNA landscape in the enteric pathogen Yersinia pseudotuberculosis, we combined lead(II) acetate-mediated cleavage of single-stranded RNA regions with high-throughput sequencing. This new approach, termed 'Lead-seq', provides structural information independent of base identity. We show that the method recapitulates secondary structures of tRNAs, RNase P RNA, tmRNA, 16S rRNA and the rpsT 5'-untranslated region, and that it reveals global structural features of mRNAs. The application of Lead-seq to Y. pseudotuberculosis cells grown at two different temperatures unveiled the first temperature-responsive in vivo RNA structurome of a bacterial pathogen. The translation of candidate genes derived from this approach was confirmed to be temperature regulated. Overall, this study establishes Lead-seq as complementary approach to interrogate intracellular RNA structures on a global scale.

Secondary Structure of Chloroplast mRNAs In Vivo and In Vitro

mRNA secondary structure can influence gene expression, e.g., by influencing translation initiation. The probing of in vivo mRNA secondary structures is therefore necessary to understand what determines the efficiency and regulation of gene expression. Here, in vivo mRNA secondary structure was analyzed using dimethyl sulfate (DMS)-MaPseq and compared to in vitro-folded RNA. We used an approach to analyze specific, full-length transcripts. To test this approach, we chose low, medium, and high abundant mRNAs. We included both monocistronic and multicistronic transcripts. Because of the slightly alkaline pH of the chloroplast stroma, we could probe all four nucleotides with DMS. The structural information gained was evaluated using the known structure of the plastid 16S rRNA. This demonstrated that the results obtained for adenosines and cytidines were more reliable than for guanosines and uridines. The majority of mRNAs analyzed were less structured in vivo than in vitro. The in vivo secondary structure of the translation initiation region of most tested genes appears to be optimized for high translation efficiency.

The chemistry and applications of RNA 2'-OH acylation

RNA is a versatile biomolecule with a broad range of biological functions that go far beyond its initially described role as a simple information carrier. The development of chemical methods to control, manipulate and modify RNA has the potential to yield new insights into its many functions and properties. Traditionally, most of these methods involved the chemical modification of RNA structure using solid-state synthesis or enzymatic transformations. However, over the past 15 years, the direct functionalization of RNA by selective acylation of the 2'-hydroxyl (2'-OH) group has emerged as a powerful alternative that enables the simple modification of both synthetic and transcribed RNAs. In this Review, we discuss the chemical properties and design of effective reagents for RNA 2'-OH acylation, highlighting the unique problem of 2'-OH reactivity in the presence of water. We elaborate on how RNA 2'-OH acylation is being exploited to develop selective chemical probes that enable interrogation of RNA structure and function, and describe new developments and applications in the field.