Probing RNA structure with chemical reagents and enzymes

The RNA interactome in the Hallmarks of Cancer

Ribonucleic acid (RNA) molecules are indispensable for cellular homeostasis in healthy and malignant cells. However, the functions of RNA extend well beyond that of a protein-coding template. Rather, both coding and non-coding RNA molecules function through critical interactions with a plethora of cellular molecules, including other RNAs, DNA, and proteins. Deconvoluting this RNA interactome, including the interacting partners, the nature of the interaction, and dynamic changes of these interactions in malignancies has yielded fundamental advances in knowledge and are emerging as a novel therapeutic strategy in cancer. Here, we present an RNA-centric review of recent advances in the field of RNA-RNA, RNA-protein, and RNA-DNA interactomic network analysis and their impact across the Hallmarks of Cancer. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

RNA structural probing of guanine and uracil nucleotides in yeast

RNA structure can be essential for its cellular function. Therefore, methods to investigate the structure of RNA in vivo are of great importance for understanding the role of cellular RNAs. RNA structure probing is an indirect method to asess the three-dimensional structure of RNA by analyzing the reactivity of different nucleotides to chemical modifications. Dimethyl sulfate (DMS) is a well-established compound that reports on base pairing context of adenine (A) and cytidine (C) in-vitro and in-vivo, but is not reactive to guanine (G) or uracil (U). Recently, new compounds were used to modify Gs and Us in plant, bacteria, and human cells. To complement the scope of RNA structural probing by chemical modifications in the model organism yeast, we analyze the effectiveness of guanine modification by the glyoxal family in Saccharomyces cerevisiae and Candida albicans. We show that within glyoxal family of compounds, phenylglyoxal (PGO) is the best guanine probe for structural probing in S. cerevisiae and C. albicans. Further, we show that PGO treatment does not affect the processing of different RNA species in the cell and is not toxic for the cells under the conditions we have established for RNA structural probing. We also explore the effectiveness of uracil modification by Cyclohexyl-3-(2-Morpholinoethyl) Carbodiimide metho-p-Toluenesulfonate (CMCT) in vivo and demonstrate that uracils can be modified by CMCT in S. cerevisiae in vivo. Our results provide the conditions for in vivo probing the reactivity of guanine and uracil nucleotides in RNA structures in yeast and offer a valuable tool for studying RNA structure and function in two widely used yeast model systems.

Using 5' 32 P-labeled Primer and Reverse Transcription to Probe RNA Structure

RNA molecules perform numerous cellular functions necessary for cell viability, some of which can depend on the RNA's structure. Therefore, it is important to study RNA structure and RNA folding to better understand the molecular basis of these functions. RNA chemical mapping strategies to elucidate RNA structural changes involve using chemical reagents that form adducts or cleave RNA. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) measures RNA flexibility by modification of the 2' hydroxyl groups of flexible nucleotides. These RNA adducts can be detected using ³² P-labeled primers and reverse transcription (RT) followed by PAGE analysis. This strategy can reveal the base-paired regions of the RNA and provide insight into tertiary structure and solvent accessibility. This protocol provides a method to interrogate RNA structure using furoyl acylimidazole (FAI). © 2023 Wiley Periodicals LLC. Basic Protocol 1: Reverse transcription (RT) primer labeling with ³² P radionuclide Basic Protocol 2: Characterization of RNA structure with radiolabeled primer and reverse transcription (RT).

Overview of Chemical Methods to Probe RNA Structure with Radionucleotides

Structural features of RNA play an important role in its capability to perform various functions in biological systems. To probe structural features, chemical probes are used to conjugate or cleave RNA at solvent-accessible sites, differentiating between flexible and constrained regions. These conjugates or cleaved products are then detected using reverse transcription (RT), where enzymatic RNA-dependent DNA primer extension is abruptly halted at the conjugation site or cleavage site. Here, we provide an overview of methods to probe RNA structure in vitro using radioactively labeled DNA primers, which provide a highly sensitive method to visualize RT stop sites with gel electrophoresis. © 2023 Wiley Periodicals LLC.

G-quadruplexes from non-coding RNAs

Non-coding RNAs (ncRNAs) are significant regulators of gene expression in a wide range of biological processes, such as transcription, RNA maturation, or translation. ncRNAs interplay with proteins or other RNAs through not only classical sequence-based mechanisms but also unique higher-order structures such as RNA G-quadruplexes (rG4s). rG4s are predictably formed in guanine-rich sequences and are closely related to various human diseases, such as tumors, neurodegenerative diseases, and infections. This review focuses on the vital role of rG4s in ncRNAs, particularly lncRNAs and miRNAs. We outline the dynamic balance between rG4s and RNA stem-loop/hairpin structures and the interplay between ncRNAs and interactors, thereby modulating gene expression and disease progression. A complete understanding of the biological regulatory role and mechanism of rG4s in ncRNAs affirms the critical importance of folding into the appropriate three-dimensional structure in maintaining or modulating the functions of ncRNAs. It makes them novel therapeutic targets for adjusting potential-G4-containing-ncRNAs-associated diseases.

Secondary Structure of Subgenomic RNA M of SARS-CoV-2

SARS-CoV-2 belongs to the Coronavirinae family. Like other coronaviruses, SARS-CoV-2 is enveloped and possesses a positive-sense, single-stranded RNA genome of ~30 kb. Genomic RNA is used as the template for replication and transcription. During these processes, positive-sense genomic RNA (gRNA) and subgenomic RNAs (sgRNAs) are created. Several studies presented the importance of the genomic RNA secondary structure in SARS-CoV-2 replication. However, the structure of sgRNAs has remained largely unsolved so far. In this study, we probed the sgRNA M model of SARS-CoV-2 in vitro. The presented model molecule includes 5'UTR and a coding sequence of gene M. This is the first experimentally informed secondary structure model of sgRNA M, which presents features likely to be important in sgRNA M function. The knowledge of sgRNA M structure provides insights to better understand virus biology and could be used for designing new therapeutics.

RNA G-quadruplexes (rG4s): genomics and biological functions

G-quadruplexes (G4s) are non-classical DNA or RNA secondary structures that have been first observed decades ago. Over the years, these four-stranded structural motifs have been demonstrated to have significant regulatory roles in diverse biological processes, but challenges remain in detecting them globally and reliably. Compared to DNA G4s (dG4s), the study of RNA G4s (rG4s) has received less attention until recently. In this review, we will summarize the innovative high-throughput methods recently developed to detect rG4s on a transcriptome-wide scale, highlight the many novel and important functions of rG4 being discovered in vivo across the tree of life, and discuss the key biological questions to be addressed in the near future.

Decoding LncRNAs

Non-coding RNAs (ncRNAs) have been considered as unimportant additions to the transcriptome. Yet, in light of numerous studies, it has become clear that ncRNAs play important roles in development, health and disease. Long-ignored, long non-coding RNAs (lncRNAs), ncRNAs made of more than 200 nucleotides have gained attention due to their involvement as drivers or suppressors of a myriad of tumours. The detailed understanding of some of their functions, structures and interactomes has been the result of interdisciplinary efforts, as in many cases, new methods need to be created or adapted to characterise these molecules. Unlike most reviews on lncRNAs, we summarize the achievements on lncRNA studies by taking into consideration the approaches for identification of lncRNA functions, interactomes, and structural arrangements. We also provide information about the recent data on the involvement of lncRNAs in diseases and present applications of these molecules, especially in medicine.

From structure to function: Route to understanding lncRNA mechanism

RNAs have emerged as a major target for diagnostics and therapeutics approaches. Regulatory nonprotein-coding RNAs (ncRNAs) in particular display remarkable versatility. They can fold into complex structures and interact with proteins, DNA, and other RNAs, thus modulating activity, localization, or interactome of multi-protein complexes. Thus, ncRNAs confer regulatory plasticity and represent a new layer of regulatory control. Interestingly, long noncoding RNAs (lncRNAs) tend to acquire complex secondary and tertiary structures and their function-in many cases-is dependent on structural conservation rather than primary sequence conservation. Whereas for many proteins, structure and its associated function are closely connected, for lncRNAs, the structural domains that determine functionality and its interactome are still not well understood. Numerous approaches for analyzing the structural configuration of lncRNAs have been developed recently. Here, will provide an overview of major experimental approaches used in the field, and discuss the potential benefit of using combinatorial strategies to analyze lncRNA modes of action based on structural information.

An evolutionarily conserved RNA structure in the functional core of the lincRNA Cyrano

The wide prevalence and regulated expression of long noncoding RNAs (lncRNAs) highlight their functional roles, but the molecular basis for their activities and structure-function relationships remains to be investigated, with few exceptions. Among the relatively few lncRNAs conserved over significant evolutionary distances is the long intergenic noncoding RNA (lincRNA) Cyrano (orthologous to human OIP5-AS1), which contains a region of 300 highly conserved nucleotides within tetrapods, which in turn contains a functional stretch of 26 nt of deep conservation. This region binds to and facilitates the degradation of the microRNA miR-7, a short ncRNA with multiple cellular functions, including modulation of oncogenic expression. We probed the secondary structure of Cyrano in vitro and in cells using chemical and enzymatic probing, and validated the results using comparative sequence analysis. At the center of the functional core of Cyrano is a cloverleaf structure maintained over the >400 million years of divergent evolution that separates fish and primates. This strikingly conserved motif provides interaction sites for several RNA-binding proteins and masks a conserved recognition site for miR-7. Conservation in this region strongly suggests that the function of Cyrano depends on the formation of this RNA structure, which could modulate the rate and efficiency of degradation of miR-7.

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.

Probing RNA structure in vivo

RNA structure underpins many essential functions in biology. New chemical reagents and techniques for probing RNA structure in living cells have emerged in recent years. High-throughput, genome-wide techniques such as Structure-seq2 and DMS-MaPseq exploit nucleobase modification by dimethylsulfate (DMS) to obtain complete structuromes, and are applicable to multiple domains of life and conditions. New reagents such as 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC), glyoxal, and nicotinoyl azide (NAz) greatly expand the capabilities of nucleobase probing in cells. Additionally, ribose-targeting reagents in selective 2'-hydroxyl acylation and primer extension (SHAPE) detect RNA flexibility in vivo. These techniques, coupled with crosslinking nucleobases in psoralen analysis of RNA interactions and structures (PARIS), provide new and diverse ways to elucidate RNA secondary and tertiary structure in vivo and genome-wide.

Probing In Vivo Structure of Individual mRNA 3' Isoforms Using Dimethyl Sulfate

The DMS region extraction and deep sequencing (DREADS) procedure was designed to probe RNA structure in vivo and to link this structural information to specific 3' isoforms. Growing cells are treated with the alkylating agent dimethyl sulfate (DMS), which enters easily into cells and modifies RNA molecules at solvent-exposed A and C residues. RNA is isolated, and sequencing libraries are constructed in a manner that preserves the identities of individual mRNA isoforms arising from alternative cleavage/polyadenylation sites. During the cDNA synthesis step of library construction, the progress of reverse transcriptase (RT) is blocked when it encounters a DMS modification on the RNA, leading to disproportionate cDNA termination adjacent to DMS-modified positions. After paired-end deep sequencing, the downstream end of each sequenced fragment is mapped to a specific cleavage/poly(A) site representing an individual mRNA 3' isoform. The upstream mapped end of the sequenced fragment defines where the RT reaction stopped. Over the population of all sequenced fragments derived from a particular isoform, A and C positions that are overrepresented next to the upstream endpoints in the DMS sample (relative to a parallel untreated control) are inferred to have been DMS modified, and hence solvent exposed. This method thus allows in vivo structural information obtained using DMS to be linked to individual mRNA 3' isoforms. © 2019 by John Wiley & Sons, Inc.

Differentiation and classification of RNA motifs using small molecule-based pattern recognition

Understanding how to design small molecules that target coding and non-coding RNA has the potential to exponentially increase the number of therapeutically-relevant druggable targets, which are currently mostly proteins. However, there is limited information on the principles at the basis of RNA recognition. In this chapter, we describe a pattern-based technique that can be used for the simultaneous elucidation of RNA motifs and small molecule features for RNA selective recognition, termed Pattern Recognition of RNA by Small Molecules (PRRSM). We provide protocols for the computational design and synthetic preparation of an RNA training set as well as how to perform the assay in plate reader format. Furthermore, we provide details on how to perform and interpret the statistical analysis and indicate possible future extensions of the technique. By combining insights into characteristics of the small molecules and of the RNA that leads to differentiation, PRRSM promises to accelerate the elucidation of the determinants at the basis of RNA recognition.

Unveiling the druggable RNA targets and small molecule therapeutics

The increasing appreciation for the crucial roles of RNAs in infectious and non-infectious human diseases makes them attractive therapeutic targets. Coding and non-coding RNAs frequently fold into complex conformations which, if effectively targeted, offer opportunities to therapeutically modulate numerous cellular processes, including those linked to undruggable protein targets. Despite the considerable skepticism as to whether RNAs can be targeted with small molecule therapeutics, overwhelming evidence suggests the challenges we are currently facing are not outside the realm of possibility. In this review, we highlight the most recent advances in molecular techniques that have sparked a revolution in understanding the RNA structure-to-function relationship. We bring attention to the application of these modern techniques to identify druggable RNA targets and to assess small molecule binding specificity. Finally, we discuss novel screening methodologies that support RNA drug discovery and present examples of therapeutically valuable RNA targets.

Secondary structure of the segment 5 genomic RNA of influenza A virus and its application for designing antisense oligonucleotides

Influenza virus causes seasonal epidemics and dangerous pandemic outbreaks. It is a single stranded (-)RNA virus with a segmented genome. Eight segments of genomic viral RNA (vRNA) form the virion, which are then transcribed and replicated in host cells. The secondary structure of vRNA is an important regulator of virus biology and can be a target for finding new therapeutics. In this paper, the secondary structure of segment 5 vRNA is determined based on chemical mapping data, free energy minimization and structure-sequence conservation analysis for type A influenza. The revealed secondary structure has circular folding with a previously reported panhandle motif and distinct novel domains. Conservations of base pairs is 87% on average with many structural motifs that are highly conserved. Isoenergetic microarray mapping was used to additionally validate secondary structure and to discover regions that easy bind short oligonucleotides. Antisense oligonucleotides, which were designed based on modeled secondary structure and microarray mapping, inhibit influenza A virus proliferation in MDCK cells. The most potent oligonucleotides lowered virus titer by ~90%. These results define universal for type A structured regions that could be important for virus function, as well as new targets for antisense therapeutics.

UGA stop codon readthrough to translate intergenic region of Plautia stali intestine virus does not require RNA structures forming internal ribosomal entry site

The translation of capsid proteins of Plautia stali intestine virus (PSIV), encoded in its second open reading frame (ORF2), is directed by an internal ribosomal entry site (IRES) located in the intergenic region (IGR). Owing to the specific properties of PSIV IGR in terms of nucleotide length and frame organization, capsid proteins are also translated via stop codon readthrough in mammalian cultured cells as an extension of translation from the first ORF (ORF1) and IGR. To delineate stop codon readthrough in PSIV, we determined requirements of cis-acting elements through a molecular genetics approach applied in both cell-free translation systems and cultured cells. Mutants with deletions from the 3' end of IGR revealed that almost none of the sequence of IGR is necessary for readthrough, apart from the 5'-terminal codon CUA. Nucleotide replacement of this CUA trinucleotide or change of the termination codon from UGA severely impaired readthrough. Chemical mapping of the IGR region of the most active 3' deletion mutant indicated that this defined minimal element UGACUA, together with its downstream sequence, adopts a single-stranded conformation. Stimulatory activities of downstream RNA structures identified to date in gammaretrovirus, coltivirus, and alphavirus were not detected in the context of PSIV IGR, despite the presence of structures for IRES. To our knowledge, PSIV IGR is the first example of stop codon readthrough that is solely defined by the local hexamer sequence, even though the sequence is adjacent to an established region of RNA secondary/tertiary structures.

In vivo RNA structural probing of uracil and guanine base-pairing by 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC)

Many biological functions performed by RNAs arise from their in vivo structures. The structure of the same RNA can differ in vitro and in vivo owing in part to the influence of molecules ranging from protons to secondary metabolites to proteins. Chemical reagents that modify the Watson-Crick (WC) face of unprotected RNA bases report on the absence of base-pairing and so are of value to determining structures adopted by RNAs. Reagents have thus been sought that can report on the native RNA structures that prevail in living cells. Dimethyl sulfate (DMS) and glyoxal penetrate cell membranes and inform on RNA secondary structure in vivo through modification of adenine (A), cytosine (C), and guanine (G) bases. Uracil (U) bases, however, have thus far eluded characterization in vivo. Herein, we show that the water-soluble carbodiimide 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC) is capable of modifying the WC face of U and G in vivo, favoring the former nucleobase by a factor of ∼1.5, and doing so in the eukaryote rice, as well as in the Gram-negative bacterium Escherichia coli While both EDC and glyoxal target Gs, EDC reacts with Gs in their typical neutral state, while glyoxal requires Gs to populate the rare anionic state. EDC may thus be more generally useful; however, comparison of the reactivity of EDC and glyoxal may allow the identification of Gs with perturbed pK_as in vivo and genome-wide. Overall, use of EDC with DMS allows in vivo probing of the base-pairing status of all four RNA bases.

Extensive Structural Differences of Closely Related 3' mRNA Isoforms: Links to Pab1 Binding and mRNA Stability

Alternative polyadenylation generates numerous 3' mRNA isoforms that can vary in biological properties, such as stability and localization. We developed methods to obtain transcriptome-scale structural information and protein binding on individual 3' mRNA isoforms in vivo. Strikingly, near-identical mRNA isoforms can possess dramatically different structures throughout the 3' UTR. Analyses of identical mRNAs in different species or refolded in vitro indicate that structural differences in vivo are often due to trans-acting factors. The level of Pab1 binding to poly(A)-containing isoforms is surprisingly variable, and differences in Pab1 binding correlate with the extent of structural variation for closely spaced isoforms. A pattern encompassing single-strandedness near the 3' terminus, double-strandedness of the poly(A) tail, and low Pab1 binding is associated with mRNA stability. Thus, individual 3' mRNA isoforms can be remarkably different physical entities in vivo. Sequences responsible for isoform-specific structures, differential Pab1 binding, and mRNA stability are evolutionarily conserved, indicating biological function.

Time-Resolved Hydroxyl Radical Footprinting of RNA with X-Rays

RNA footprinting by hydroxyl radical cleavage provides 'snapshots' of RNA tertiary structure or protein interactions that bury the RNA backbone. Generation of hydroxyl radicals with a high-flux synchrotron X-ray beam provides analysis on a short timescale (5-100 msec), which enables the structures of folding intermediates or other transient conformational states to be determined in biochemical solutions or cells. This article provides protocols for using synchrotron beamlines for hydroxyl radical footprinting. © 2018 by John Wiley & Sons, Inc.

Reactivity and Specificity of RNase T1, RNase A, and RNase H toward Oligonucleotides of RNA Containing 8-Oxo-7,8-dihydroguanosine

Understanding how oxidatively damaged RNA interacts with ribonucleases is important because of its proposed role in the development and progression of disease. Thus, understanding structural aspects of RNA containing lesions generated under oxidative stress, as well as its interactions with other biopolymers, is fundamental. We explored the reactivity of RNase A, RNase T₁, and RNase H toward oligonucleotides of RNA containing 8-oxo-7,8-dihydroguanosine (8oxoG). This is the first example that addresses this relationship and will be useful for understanding (1) how these RNases can be used to characterize the structural impact that this lesion has on RNA and (2) how oxidatively modified RNA may be handled intracellularly. 8-OxoG was incorporated into 10-16-mers of RNA, and its reactivity with each ribonuclease was assessed via electrophoretic analyses, circular dichroism, and the use of other C8-purine-modified analogues (8-bromoguanosine, 8-methoxyguanosine, and 8-oxoadenosine). RNase T₁ does not recognize sites containing 8-oxoG, while RNase A recognizes and cleaves RNA at positions containing this lesion while differentiating if it is involved in H-bonding. The selectivity of RNase A followed the order C > 8-oxoG ≈ U. In addition, isothermal titration calorimetry showed that an 8-oxoG-C3'-methylphosphate derivative can inhibit RNase A activity. Cleavage patterns obtained from RNase H displayed changes in reactivity in a sequence- and concentration-dependent manner and displayed recognition at sites containing the modification in some cases. These data will aid in understanding how this modification affects reactivity with ribonucleases and will enable the characterization of global and local structural changes in oxidatively damaged RNA.

Influenza virus segment 5 (+)RNA - secondary structure and new targets for antiviral strategies

Influenza A virus is a threat for humans due to seasonal epidemics and occasional pandemics. This virus can generate new strains that are dangerous through nucleotide/amino acid changes or through segmental recombination of the viral RNA genome. It is important to gain wider knowledge about influenza virus RNA to create new strategies for drugs that will inhibit its spread. Here, we present the experimentally determined secondary structure of the influenza segment 5 (+)RNA. Two RNAs were studied: the full-length segment 5 (+)RNA and a shorter construct containing only the coding region. Chemical mapping data combined with thermodynamic energy minimization were used in secondary structure prediction. Sequence/structure analysis showed that the determined secondary structure of segment 5 (+)RNA is mostly conserved between influenza virus type A strains. Microarray mapping and RNase H cleavage identified accessible sites for oligonucleotides in the revealed secondary structure of segment 5 (+)RNA. Antisense oligonucleotides were designed based on the secondary structure model and tested against influenza virus in cell culture. Inhibition of influenza virus proliferation was noticed, identifying good targets for antisense strategies. Effective target sites fall within two domains, which are conserved in sequence/structure indicating their importance to the virus.

Interpreting Reverse Transcriptase Termination and Mutation Events for Greater Insight into the Chemical Probing of RNA

Chemical probing has the power to provide insight into RNA conformation in vivo and in vitro, but interpreting the results depends on methods to detect the chemically modified nucleotides. Traditionally, the presence of modified bases was inferred from their ability to halt reverse transcriptase during primer extension and the locations of termination sites observed by electrophoresis or sequencing. More recently, modification-induced mutations have been used as a readout for chemical probing data. Given the variable propensity for mismatch incorporation and read-through with different reverse transcriptases, we examined how termination and mutation events compare to each other in the same chemical probing experiments. We found that mutations and terminations induced by dimethyl sulfate probing are both specific for methylated bases, but these two measures have surprisingly little correlation and represent largely nonoverlapping indicators of chemical modification data. We also show that specific biases for modified bases depend partly on local sequence context and that different reverse transcriptases show different biases toward reading a modification as a stop or a mutation. These results support approaches that incorporate analysis of both termination and mutation events into RNA probing experiments.

Identification of N6-methyladenosine reader proteins

The reversible N⁶-methyladenosine (m⁶A) modification of eukaryotic messenger RNAs (mRNAs) is a widespread regulatory mechanism that impacts every step in the mRNA life cycle. The effect of m⁶A on mRNA fate depends on the binding of "m⁶A reader" proteins - RNA binding proteins that specifically bind to RNAs containing m⁶A. Here, we describe an RNA pull-down method that can be used to identify novel m⁶A reader proteins starting from a known m⁶A-modified site in cellular or viral RNA. We further describe how a combination of immunoprecipitation-based sequencing methods can be used to identify m⁶A-modified sites bound by an m⁶A reader protein on a transcriptome-wide level. The discovery of new m⁶A reader proteins and their m⁶A-modified targets would provide further insight into the mechanisms and functions of m⁶A in the cell.

Understanding the Functions of Long Non-Coding RNAs through Their Higher-Order Structures

Although thousands of long non-coding RNAs (lncRNAs) have been discovered in eukaryotes, very few molecular mechanisms have been characterized due to an insufficient understanding of lncRNA structure. Therefore, investigations of lncRNA structure and subsequent elucidation of the regulatory mechanisms are urgently needed. However, since lncRNA are high molecular weight molecules, which makes their crystallization difficult, obtaining information about their structure is extremely challenging, and the structures of only several lncRNAs have been determined so far. Here, we review the structure-function relationships of the widely studied lncRNAs found in the animal and plant kingdoms, focusing on the principles and applications of both in vitro and in vivo technologies for the study of RNA structures, including dimethyl sulfate-sequencing (DMS-seq), selective 2'-hydroxyl acylation analyzed by primer extension-sequencing (SHAPE-seq), parallel analysis of RNA structure (PARS), and fragmentation sequencing (FragSeq). The aim of this review is to provide a better understanding of lncRNA biological functions by studying them at the structural level.

Transcriptome-wide interrogation of RNA secondary structure in living cells with icSHAPE

icSHAPE (in vivo click selective 2-hydroxyl acylation and profiling experiment) captures RNA secondary structure at a transcriptome-wide level by measuring nucleotide flexibility at base resolution. Living cells are treated with the icSHAPE chemical NAI-N3 followed by selective chemical enrichment of NAI-N3-modified RNA, which provides an improved signal-to-noise ratio compared with similar methods leveraging deep sequencing. Purified RNA is then reverse-transcribed to produce cDNA, with SHAPE-modified bases leading to truncated cDNA. After deep sequencing of cDNA, computational analysis yields flexibility scores for every base across the starting RNA population. The entire experimental procedure can be completed in ∼5 d, and the sequencing and bioinformatics data analysis take an additional 4-5 d with no extensive computational skills required. Comparing in vivo and in vitro icSHAPE measurements can reveal in vivo RNA-binding protein imprints or facilitate the dissection of RNA post-transcriptional modifications. icSHAPE reactivities can additionally be used to constrain and improve RNA secondary structure prediction models.

Probing-directed identification of novel structured RNAs

Transcripts often harbor RNA elements, which regulate cell processes co- or post-transcriptionally. The functions of many regulatory RNA elements depend on their structure, thus it is important to determine the structure as well as to scan genomes for structured elements. State of the art ab initio approaches to predict structured RNAs rely on DNA sequence analysis. They use 2 major types of information inferred from a sequence: thermodynamic stability of an RNA structure and evolutionary footprints of base-pair interactions. In recent years, chemical probing of RNA has arisen as an alternative source of structural information. RNA probing experiments detect positions accessible to specific types of chemicals or enzymes indicating their propensity to be in a paired or unpaired state. There exist several strategies to integrate probing data into RNA secondary structure prediction algorithms that substantially improve the prediction quality. However, whether and how probing data could contribute to detection of structured RNAs remains an open question. We previously developed the energy-based approach RNASurface to detect locally optimal structured RNA elements. Here, we integrate probing data into the RNASurface energy model using a general framework. We show that the use of experimental data allows for better discrimination of ncRNAs from other transcripts. Application of RNASurface to genome-wide analysis of the human transcriptome with PARS data identifies previously undetectable segments, with evidence of functionality for some of them.

Experiment-Assisted Secondary Structure Prediction with RNAstructure

Experimental probing data can be used to improve the accuracy of RNA secondary structure prediction. The software package RNAstructure can take advantage of enzymatic cleavage data, FMN cleavage data, traditional chemical modification reactivity data, and SHAPE reactivity data for secondary structure modeling. This chapter provides protocols for using experimental probing data with RNAstructure to restrain or constrain RNA secondary structure prediction.

Progress and challenges for chemical probing of RNA structure inside living cells

Proper gene expression is essential for the survival of every cell. Once thought to be a passive transporter of genetic information, RNA has recently emerged as a key player in nearly every pathway in the cell. A full description of its structure is critical to understanding RNA function. Decades of research have focused on utilizing chemical tools to interrogate the structures of RNAs, with recent focus shifting to performing experiments inside living cells. This Review will detail the design and utility of chemical reagents used in RNA structure probing. We also outline how these reagents have been used to gain a deeper understanding of RNA structure in vivo. We review the recent merger of chemical probing with deep sequencing. Finally, we outline some of the hurdles that remain in fully characterizing the structure of RNA inside living cells, and how chemical biology can uniquely tackle such challenges.

Current perspectives on RNA secondary structure probing

The range of roles played by structured RNAs in biological systems is vast. At the same time as we are learning more about the importance of RNA structure, recent advances in reagents, methods and technology mean that RNA secondary structural probing has become faster and more accurate. As a result, the capabilities of laboratories that already perform this type of structural analysis have increased greatly, and it has also become more widely accessible. The present review summarizes established and recently developed techniques. The information we can derive from secondary structural analysis is assessed, together with the areas in which we are likely to see exciting developments in the near future.

Computational analysis of RNA structures with chemical probing data

RNAs play various roles, not only as the genetic codes to synthesize proteins, but also as the direct participants of biological functions determined by their underlying high-order structures. Although many computational methods have been proposed for analyzing RNA structures, their accuracy and efficiency are limited, especially when applied to the large RNAs and the genome-wide data sets. Recently, advances in parallel sequencing and high-throughput chemical probing technologies have prompted the development of numerous new algorithms, which can incorporate the auxiliary structural information obtained from those experiments. Their potential has been revealed by the secondary structure prediction of ribosomal RNAs and the genome-wide ncRNA function annotation. In this review, the existing probing-directed computational methods for RNA secondary and tertiary structure analysis are discussed.

High-throughput characterization of protein-RNA interactions

RNA-binding proteins (RBPs) are important regulators of eukaryotic gene expression. Genomes typically encode dozens to hundreds of proteins containing RNA-binding domains, which collectively recognize diverse RNA sequences and structures. Recent advances in high-throughput methods for assaying the targets of RBPs in vitro and in vivo allow large-scale derivation of RNA-binding motifs as well as determination of RNA–protein interactions in living cells. In parallel, many computational methods have been developed to analyze and interpret these data. The interplay between RNA secondary structure and RBP binding has also been a growing theme. Integrating RNA–protein interaction data with observations of post-transcriptional regulation will enhance our understanding of the roles of these important proteins.

Secondary structure of a conserved domain in an intron of influenza A M1 mRNA

Influenza A virus utilizes RNA throughout infection. Little is known, however, about the roles of RNA structures. A previous bioinformatics survey predicted multiple regions of influenza A virus that are likely to generate evolutionarily conserved and stable RNA structures. One predicted conserved structure is in the pre-mRNA coding for essential proteins, M1 and M2. This structure starts 79 nucleotides downstream of the M2 mRNA 5' splice site. Here, a combination of biochemical structural mapping, mutagenesis, and NMR confirms the predicted three-way multibranch structure of this RNA. Imino proton NMR spectra reveal no change in secondary structure when 80 mM KCl is supplemented with 4 mM MgCl2. Optical melting curves in 1 M NaCl and in 100 mM KCl with 10 mM MgCl2 are very similar, with melting temperatures ∼14 °C higher than that for 100 mM KCl alone. These results provide a firm basis for designing experiments and potential therapeutics to test for function in cell culture.

Insights into RNA structure and function from genome-wide studies

A comprehensive understanding of RNA structure will provide fundamental insights into the cellular function of both coding and non-coding RNAs. Although many RNA structures have been analysed by traditional biophysical and biochemical methods, the low-throughput nature of these approaches has prevented investigation of the vast majority of cellular transcripts. Triggered by advances in sequencing technology, genome-wide approaches for probing the transcriptome are beginning to reveal how RNA structure affects each step of protein expression and RNA stability. In this Review, we discuss the emerging relationships between RNA structure and the regulation of gene expression.

Massively parallel RNA chemical mapping with a reduced bias MAP-seq protocol

Chemical mapping methods probe RNA structure by revealing and leveraging correlations of a nucleotide's structural accessibility or flexibility with its reactivity to various chemical probes. Pioneering work by Lucks and colleagues has expanded this method to probe hundreds of molecules at once on an Illumina sequencing platform, obviating the use of slab gels or capillary electrophoresis on one molecule at a time. Here, we describe optimizations to this method from our lab, resulting in the MAP-seq protocol (Multiplexed Accessibility Probing read out through sequencing), version 1.0. The protocol permits the quantitative probing of thousands of RNAs at once, by several chemical modification reagents, on the time scale of a day using a tabletop Illumina machine. This method and a software package MAPseeker ( http://simtk.org/home/map_seeker ) address several potential sources of bias, by eliminating PCR steps, improving ligation efficiencies of ssDNA adapters, and avoiding problematic heuristics in prior algorithms. We hope that the step-by-step description of MAP-seq 1.0 will help other RNA mapping laboratories to transition from electrophoretic to next-generation sequencing methods and to further reduce the turnaround time and any remaining biases of the protocol.

Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo

RNA plays a dual role as an informational molecule and a direct effector of biological tasks. The latter function is enabled by RNA’s ability to adopt complex secondary and tertiary folds and thus has motivated extensive computational^1–2 and experimental^3–8 efforts for determining RNA structures. Existing approaches for evaluating RNA structure have been largely limited to in vitro systems, yet the thermodynamic forces which drive RNA folding in vitro may not be sufficient to predict stable RNA structures in vivo⁵. Indeed, the presence of RNA binding proteins and ATP-dependent helicases can influence which structures are present inside cells. Here we present an approach for globally monitoring RNA structure in native conditions in vivo with single nucleotide precision. This method is based on in vivo modification with dimethyl sulfate (DMS), which reacts with unpaired adenine and cytosine residues⁹, followed by deep sequencing to monitor modifications. Our data from yeast and mammalian cells are in excellent agreement with known mRNA structures and with the high-resolution crystal structure of the Saccharomyces cerevisiae ribosome¹⁰. Comparison between in vivo and in vitro data reveals that in rapidly dividing cells there are vastly fewer structured mRNA regions in vivo than in vitro. Even thermostable RNA structures are often denatured in cells, highlighting the importance of cellular processes in regulating RNA structure. Indeed, analysis of mRNA structure under ATP-depleted conditions in yeast reveals that energy-dependent processes strongly contribute to the predominantly unfolded state of mRNAs inside cells. Our studies broadly enable the functional analysis of physiological RNA structures and reveal that, in contrast to the Anfinsen view of protein folding, thermodynamics play an incomplete role in determining mRNA structure in vivo.

Tackling structures of long noncoding RNAs

RNAs are important catalytic machines and regulators at every level of gene expression. A new class of RNAs has emerged called long non-coding RNAs, providing new insights into evolution, development and disease. Long non-coding RNAs (lncRNAs) predominantly found in higher eukaryotes, have been implicated in the regulation of transcription factors, chromatin-remodeling, hormone receptors and many other processes. The structural versatility of RNA allows it to perform various functions, ranging from precise protein recognition to catalysis and metabolite sensing. While major housekeeping RNA molecules have long been the focus of structural studies, lncRNAs remain the least characterized class, both structurally and functionally. Here, we review common methodologies used to tackle RNA structure, emphasizing their potential application to lncRNAs. When considering the complexity of lncRNAs and lack of knowledge of their structure, chemical probing appears to be an indispensable tool, with few restrictions in terms of size, quantity and heterogeneity of the RNA molecule. Probing is not constrained to in vitro analysis and can be adapted to high-throughput sequencing platforms. Significant efforts have been applied to develop new in vivo chemical probing reagents, new library construction protocols for sequencing platforms and improved RNA prediction software based on the experimental evidence.

The 3' splice site of influenza A segment 7 mRNA can exist in two conformations: a pseudoknot and a hairpin

The 3′ splice site of influenza A segment 7 is used to produce mRNA for the M2 ion-channel protein, which is critical to the formation of viable influenza virions. Native gel analysis, enzymatic/chemical structure probing, and oligonucleotide binding studies of a 63 nt fragment, containing the 3′ splice site, key residues of an SF2/ASF splicing factor binding site, and a polypyrimidine tract, provide evidence for an equilibrium between pseudoknot and hairpin structures. This equilibrium is sensitive to multivalent cations, and can be forced towards the pseudoknot by addition of 5 mM cobalt hexammine. In the two conformations, the splice site and other functional elements exist in very different structural environments. In particular, the splice site is sequestered in the middle of a double helix in the pseudoknot conformation, while in the hairpin it resides in a two-by-two nucleotide internal loop. The results suggest that segment 7 mRNA splicing can be controlled by a conformational switch that exposes or hides the splice site.

Structural architecture of the human long non-coding RNA, steroid receptor RNA activator

While functional roles of several long non-coding RNAs (lncRNAs) have been determined, the molecular mechanisms are not well understood. Here, we report the first experimentally derived secondary structure of a human lncRNA, the steroid receptor RNA activator (SRA), 0.87 kB in size. The SRA RNA is a non-coding RNA that coactivates several human sex hormone receptors and is strongly associated with breast cancer. Coding isoforms of SRA are also expressed to produce proteins, making the SRA gene a unique bifunctional system. Our experimental findings (SHAPE, in-line, DMS and RNase V1 probing) reveal that this lncRNA has a complex structural organization, consisting of four domains, with a variety of secondary structure elements. We examine the coevolution of the SRA gene at the RNA structure and protein structure levels using comparative sequence analysis across vertebrates. Rapid evolutionary stabilization of RNA structure, combined with frame-disrupting mutations in conserved regions, suggests that evolutionary pressure preserves the RNA structural core rather than its translational product. We perform similar experiments on alternatively spliced SRA isoforms to assess their structural features.

Isoenergetic microarrays to study the structure and interactions of DsrA and OxyS RNAs in two- and three-component complexes

Information on the secondary structure and interactions of RNA is important to understand the biological function of RNA as well as in applying RNA as a tool for therapeutic purposes. Recently, the isoenergetic microarray mapping method was developed to improve the prediction of RNA secondary structure. Herein, for the first time, isoenergetic microarrays were used to study the binding of RNA to protein or other RNAs as well as the interactions of two different RNAs and protein in a three-component complex. The RNAs used as models were the regulatory DsrA and OxyS RNAs from Escherichia coli, the fragments of their target mRNAs (fhlA and rpoS), and their complexes with Hfq protein. The collected results showed the advantages and some limitations of microarray mapping.

Binding and cleavage of unstructured RNA by nuclear RNase P

Ribonuclease P (RNase P) is an essential endoribonuclease for which the best-characterized function is processing the 5' leader of pre-tRNAs. Compared to bacterial RNase P, which contains a single small protein subunit and a large catalytic RNA subunit, eukaryotic nuclear RNase P is more complex, containing nine proteins and an RNA subunit in Saccharomyces cerevisiae. Consistent with this, nuclear RNase P has been shown to possess unique RNA binding capabilities. To understand the unique molecular recognition of nuclear RNase P, the interaction of S. cerevisiae RNase P with single-stranded RNA was characterized. Unstructured, single-stranded RNA inhibits RNase P in a size-dependent manner, suggesting that multiple interactions are required for high affinity binding. Mixed-sequence RNAs from protein-coding regions also bind strongly to the RNase P holoenzyme. However, in contrast to poly(U) homopolymer RNA that is not cleaved, a variety of mixed-sequence RNAs have multiple preferential cleavage sites that do not correspond to identifiable consensus structures or sequences. In addition, pre-tRNA(Tyr), poly(U)(50) RNA, and mixed-sequence RNA cross-link with purified RNase P in the RNA subunit Rpr1 near the active site in "Conserved Region I," although the exact positions vary. Additional contacts between poly(U)(50) and the RNase P proteins Rpr2p and Pop4p were identified. We conclude that unstructured RNAs interact with multiple protein and RNA contacts near the RNase P RNA active site, but that cleavage depends on the nature of interaction with the active site.

Genome-wide measurement of RNA secondary structure in yeast

The structures of RNA molecules are often important for their function and regulation, yet there are no experimental techniques for genome-scale measurement of RNA structure. Here we describe a novel strategy termed parallel analysis of RNA structure (PARS), which is based on deep sequencing fragments of RNAs that were treated with structure-specific enzymes, thus providing simultaneous in vitro profiling of the secondary structure of thousands of RNA species at single nucleotide resolution. We apply PARS to profile the secondary structure of the messenger RNAs (mRNAs) of the budding yeast Saccharomyces cerevisiae and obtain structural profiles for over 3,000 distinct transcripts. Analysis of these profiles reveals several RNA structural properties of yeast transcripts, including the existence of more secondary structure over coding regions compared with untranslated regions, a three-nucleotide periodicity of secondary structure across coding regions and an anti-correlation between the efficiency with which an mRNA is translated and the structure over its translation start site. PARS is readily applicable to other organisms and to profiling RNA structure in diverse conditions, thus enabling studies of the dynamics of secondary structure at a genomic scale.

Comparisons between chemical mapping and binding to isoenergetic oligonucleotide microarrays reveal unexpected patterns of binding to the Bacillus subtilis RNase P RNA specificity domain

Microarrays with isoenergetic pentamer and hexamer 2'-O-methyl oligonucleotide probes with LNA (locked nucleic acid) and 2,6-diaminopurine substitutions were used to probe the binding sites on the RNase P RNA specificity domain of Bacillus subtilis. Unexpected binding patterns were revealed. Because of their enhanced binding free energies, isoenergetic probes can break short duplexes, merge adjacent loops, and/or induce refolding. This suggests new approaches to the rational design of short oligonucleotide therapeutics but limits the utility of microarrays for providing constraints for RNA structure determination. The microarray results are compared to results from chemical mapping experiments, which do provide constraints. Results from both types of experiments indicate that the RNase P RNA folds similarly in 1 M Na(+) and 10 mM Mg(2+).

Robust and generic RNA modeling using inferred constraints: a structure for the hepatitis C virus IRES pseudoknot domain

RNA function is dependent on its structure, yet three-dimensional folds for most biologically important RNAs are unknown. We develop a generic discrete molecular dynamics-based modeling system that uses long-range constraints inferred from diverse biochemical or bioinformatic analyses to create statistically significant (p < 0.01) nativelike folds for RNAs of known structure ranging from 45 to 158 nucleotides. We then predict the unknown structure of the hepatitis C virus internal ribosome entry site (IRES) pseudoknot domain. The resulting RNA model rationalizes independent solvent accessibility and cryo-electron microscopy structure information. The pseudoknot domain positions the AUG start codon near the mRNA channel and is tRNA-like, suggesting the IRES employs molecular mimicry as a functional strategy.

Structural basis for the methylation of G1405 in 16S rRNA by aminoglycoside resistance methyltransferase Sgm from an antibiotic producer: a diversity of active sites in m7G methyltransferases

Sgm (Sisomicin-gentamicin methyltransferase) from antibiotic-producing bacterium Micromonospora zionensis is an enzyme that confers resistance to aminoglycosides like gentamicin and sisomicin by specifically methylating G1405 in bacterial 16S rRNA. Sgm belongs to the aminoglycoside resistance methyltransferase (Arm) family of enzymes that have been recently found to spread by horizontal gene transfer among disease-causing bacteria. Structural characterization of Arm enzymes is the key to understand their mechanism of action and to develop inhibitors that would block their activity. Here we report the structure of Sgm in complex with cofactors S-adenosylmethionine (AdoMet) and S-adenosylhomocysteine (AdoHcy) at 2.0 and 2.1 Å resolution, respectively, and results of mutagenesis and rRNA footprinting, and protein-substrate docking. We propose the mechanism of methylation of G1405 by Sgm and compare it with other m⁷G methyltransferases, revealing a surprising diversity of active sites and binding modes for the same basic reaction of RNA modification. This analysis can serve as a stepping stone towards developing drugs that would specifically block the activity of Arm methyltransferases and thereby re-sensitize pathogenic bacteria to aminoglycoside antibiotics.

Crystallization of the glmS ribozyme-riboswitch

Procedures that were critical for crystallization of the glmS ribozyme-riboswitch RNA domain from the thermophilic Gram-positive bacterium Thermoanaerobacter tengcongensis are described. Experimental design based on screening multiple variant RNA sequences and techniques used to identify initial crystallization conditions were similar to those employed for most RNAs. However, serendipitous in-drop digestion of one RNA construct at a specific internucleotide linkage was crucial for the growth of high-quality glmS ribozyme crystals. Biochemical analysis of crystalline RNA identified the site of scission and guided design of an optimized bimolecular RNA construct. Finally, modifications of ionic strength and pH of solutions used for stabilization of the crystals were essential for optimal diffraction and binding of the activator glucosamine-6-phosphate, respectively. Although their details are specific to the glmS ribozyme, these general strategies may be useful for analyzing and improving crystals of other RNAs.

Isoenergetic penta- and hexanucleotide microarray probing and chemical mapping provide a secondary structure model for an RNA element orchestrating R2 retrotransposon protein function

LNA (locked nucleic acids, i.e. oligonucleotides with a methyl bridge between the 2′ oxygen and 4′ carbon of ribose) and 2,6-diaminopurine were incorporated into 2′-O-methyl RNA pentamer and hexamer probes to make a microarray that binds unpaired RNA approximately isoenergetically. That is, binding is roughly independent of target sequence if target is unfolded. The isoenergetic binding and short probe length simplify interpretation of binding to a structured RNA to provide insight into target RNA secondary structure. Microarray binding and chemical mapping were used to probe the secondary structure of a 323 nt segment of the 5′ coding region of the R2 retrotransposon from Bombyx mori (R2Bm 5′ RNA). This R2Bm 5′ RNA orchestrates functioning of the R2 protein responsible for cleaving the second strand of DNA during insertion of the R2 sequence into the genome. The experimental results were used as constraints in a free energy minimization algorithm to provide an initial model for the secondary structure of the R2Bm 5′ RNA.