Probing RNA structure and RNA-ligand complexes with chemical probes

Native Top-Down Mass Spectrometry Uncovers Two Distinct Binding Motifs of a Functional Neomycin-Sensing Riboswitch Aptamer

Understanding how ligands bind to ribonucleic acids (RNA) is important for understanding RNA recognition in biological processes and drug development. Here, we have studied neomycin B binding to neomycin-sensing riboswitch aptamer constructs by native top-down mass spectrometry (MS) using electrospray ionization (ESI) and collisionally activated dissociation (CAD). Our MS data for a 27 nt aptamer construct reveal the binding site and ligand interactions, in excellent agreement with the structure derived from nuclear magnetic resonance (NMR) studies. Strikingly, for an extended 40 nt aptamer construct, which represents the sequence with the highest regulatory factor for riboswitch function, we identified two binding motifs for neomycin B binding, one corresponding to the bulge-loop motif of the 27 nt construct and the other one in the minor groove of the lower stem, which according to the MS data are equally populated. By replacing a noncanonical with a canonical base pair in the lower stem of the 40 nt aptamer, we can reduce binding to the minor groove motif from ∼50 to ∼30%. Conversely, the introduction of a CUG/CUG motif in the lower stem shifts the binding equilibrium in favor of minor groove binding. The MS data reveal site-specific and stoichiometry-resolved information on aminoglycoside binding to RNA that is not directly accessible by other methods and underscore the role of noncanonical base pairs in RNA recognition by aminoglycosides.

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.

RNA Structure Analysis by Chemical Probing with DMS and CMCT

RNA structure is important for understanding RNA function and stability within a cell. Chemical probing is a well-established and convenient method to evaluate the structure of an RNA. Several structure-sensitive chemicals can differentiate paired and unpaired nucleotides. This chapter specifically addresses the use of DMS and CMCT. Although exhibiting different affinities, the combination of these two chemical reagents enables screening of all four nucleobases. DMS and CMCT are only reactive with exposed unpaired nucleotides. We have used this method to analyze the effect of the RNA chaperone Hfq on the conformation of the 16S rRNA. The strategy here described may be applied for the study of many other RNA-binding proteins and RNAs.

RNA Remodeling by RNA Chaperones Monitored by RNA Structure Probing

RNA structure probing enables the characterization of RNA secondary structures by established procedures such as the enzyme- or chemical-based detection of single- or double-stranded regions. A specific type of application involves the detection of changes of RNA structures and conformations that are induced by proteins with RNA chaperone activity. This chapter outlines a protocol to analyze RNA structures in vitro in the presence of an RNA-binding protein with RNA chaperone activity. For this purpose, we make use of the methylating agents dimethyl sulfate (DMS) and 1-cyclohexyl-3-(2-morpholinoethyl) carbodiimide metho-p-toluenesulfonate (CMCT). DMS and CMCT specifically modify nucleotides that are not involved in base-pairing or tertiary structure hydrogen bonding and that are not protected by a ligand such as a protein. Modified bases are identified by primer extension. As an example, we describe how the RNA chaperone activity of an isoform of the RNA-binding protein AUF1 induces the flaviviral RNA switch required for viral genome cyclization and viral replication.This chapter includes comprehensive protocols for in vitro synthesis of RNA, ³²P-5'-end labeling of DNA primers, primer extension, as well as the preparation and running of analytical gels. The described methodology should be applicable to any other RNA and protein of interest to identify protein-directed RNA remodeling.

Advances in RNA 3D Structure Modeling Using Experimental Data

RNA is a unique bio-macromolecule that can both record genetic information and perform biological functions in a variety of molecular processes, including transcription, splicing, translation, and even regulating protein function. RNAs adopt specific three-dimensional conformations to enable their functions. Experimental determination of high-resolution RNA structures using x-ray crystallography is both laborious and demands expertise, thus, hindering our comprehension of RNA structural biology. The computational modeling of RNA structure was a milestone in the birth of bioinformatics. Although computational modeling has been greatly improved over the last decade showing many successful cases, the accuracy of such computational modeling is not only length-dependent but also varies according to the complexity of the structure. To increase credibility, various experimental data were integrated into computational modeling. In this review, we summarize the experiments that can be integrated into RNA structure modeling as well as the computational methods based on these experimental data. We also demonstrate how computational modeling can help the experimental determination of RNA structure. We highlight the recent advances in computational modeling which can offer reliable structure models using high-throughput experimental data.

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.

High-throughput determination of RNA structures

RNA performs and regulates a diverse range of cellular processes, with new functional roles being uncovered at a rapid pace. Interest is growing in how these functions are linked to RNA structures that form in the complex cellular environment. A growing suite of technologies that use advances in RNA structural probes, high-throughput sequencing and new computational approaches to interrogate RNA structure at unprecedented throughput are beginning to provide insights into RNA structures at new spatial, temporal and cellular scales.

Obstacles to Scanning by RNase E Govern Bacterial mRNA Lifetimes by Hindering Access to Distal Cleavage Sites

The diversity of mRNA lifetimes in bacterial cells is difficult to reconcile with the relaxed cleavage site specificity of RNase E, the endonuclease most important for governing mRNA degradation. This enzyme has generally been thought to locate cleavage sites by searching freely in three dimensions. However, our results now show that its access to such sites in 5'-monophosphorylated RNA is hindered by obstacles-such as bound proteins or ribosomes or coaxial small RNA (sRNA) base pairing-that disrupt the path from the 5' end to those sites and prolong mRNA lifetimes. These findings suggest that RNase E searches for cleavage sites by scanning linearly from the 5'-terminal monophosphate along single-stranded regions of RNA and that its progress is impeded by structural discontinuities encountered along the way. This discovery has major implications for gene regulation in bacteria and suggests a general mechanism by which other prokaryotic and eukaryotic regulatory proteins can be controlled.

Facile synthesis and evaluation of a dual-functioning furoyl probe for in-cell SHAPE

Recent analysis of transcriptomes has revealed that RNA molecules perform a myriad of functions beyond coding for proteins. RNA molecules can fold into complex secondary and tertiary structures, which are critical for regulating their function. Selective Hydroxyl Acylation analyzed by Primer Extension, or SHAPE is a common method for probing RNA structure in and outside of cells. Recent developments in SHAPE include the design of acyl imidazole acylating electrophiles with alkyl azides to enrich the sites of SHAPE adduct formation. Enrichment is key for next-generation sequencing experiments as it dramatically improves the signal. In a recent comparison of different structures of such reagents, we realized that furoyl acylating reagents form hyper-stable ester adducts with hydroxyls. This prompted us to design, synthesize and test a novel dual-functioning SHAPE probe (FAI-N₃), which has the stable furoyl scaffold and the alkyl azide for enrichment. Herein we present the results that show FAI-N₃ is a suitable probe for RNA structure analysis by SHAPE and that it can be used for enrichment of SHAPE adducts. These results strongly demonstrate that FAI-N₃ is an ideal probe for structure probing in cells and will be very useful for sequencing-based analysis of SHAPE.

Mapping RNA Structure In Vitro with SHAPE Chemistry and Next-Generation Sequencing (SHAPE-Seq)

Mapping RNA structure with selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry has proven to be a versatile method for characterizing RNA structure in a variety of contexts. SHAPE reagents covalently modify RNAs in a structure-dependent manner to create adducts at the 2'-OH group of the ribose backbone at nucleotides that are structurally flexible. The positions of these adducts are detected using reverse transcriptase (RT) primer extension, which stops one nucleotide before the modification, to create a pool of cDNAs whose lengths reflect the location of SHAPE modification. Quantification of the cDNA pools is used to estimate the "reactivity" of each nucleotide in an RNA molecule to the SHAPE reagent. High reactivities indicate nucleotides that are structurally flexible, while low reactivities indicate nucleotides that are inflexible. These SHAPE reactivities can then be used to infer RNA structures by restraining RNA structure prediction algorithms. Here, we provide a state-of-the-art protocol describing how to perform in vitro RNA structure probing with SHAPE chemistry using next-generation sequencing to quantify cDNA pools and estimate reactivities (SHAPE-Seq). The use of next-generation sequencing allows for higher throughput, more consistent data analysis, and multiplexing capabilities. The technique described herein, SHAPE-Seq v2.0, uses a universal reverse transcription priming site that is ligated to the RNA after SHAPE modification. The introduced priming site allows for the structural analysis of an RNA independent of its sequence.

Comparative Analysis Reveals Furoyl in Vivo Selective Hydroxyl Acylation Analyzed by Primer Extension Reagents Form Stable Ribosyl Ester Adducts

RNA molecules depend on structural elements that are critical for cellular function. Chemical methods for probing RNA structure have emerged as a necessary component of characterizing RNA function. As such, understanding the limitations and idiosyncrasies of these methods is essential for their utility. Selective hydroxyl acylation has emerged as a common method for analyzing RNA structure. Ester products as a result of 2'-hydroxyl acylation can then be identified through reverse transcription or mutational enzyme profiling. The central aspect of selective hydroxyl acylation analyzed by primer extension (SHAPE) experiments is the fact that stable ester adducts are formed on the 2'-hydroxyl. Despite its importance, there has not been a direct comparison of SHAPE electrophiles for their ability to make stable RNA adducts. Herein, we conduct a systematic analysis of hydrolysis stability experiments to demonstrate that furoyl imidazole SHAPE reagents form stable ester adducts even at elevated temperatures. We also demonstrate that the acylation reaction with the furoyl acylimidaole SHAPE reagent can be controlled with dithiothreitol quenching, even in live cells. These results are important for our understanding of the biochemical details of the SHAPE experiment.

Accumulation of Stable Full-Length Circular Group I Intron RNAs during Heat-Shock

Group I introns in nuclear ribosomal RNA of eukaryotic microorganisms are processed by splicing or circularization. The latter results in formation of full-length circular introns without ligation of the exons and has been proposed to be active in intron mobility. We applied qRT-PCR to estimate the copy number of circular intron RNA from the myxomycete Didymium iridis. In exponentially growing amoebae, the circular introns are nuclear and found in 70 copies per cell. During heat-shock, the circular form is up-regulated to more than 500 copies per cell. The intron harbours two ribozymes that have the potential to linearize the circle. To understand the structural features that maintain circle integrity, we performed chemical and enzymatic probing of the splicing ribozyme combined with molecular modeling to arrive at models of the inactive circular form and its active linear counterpart. We show that the two forms have the same overall structure but differ in key parts, including the catalytic core element P7 and the junctions at which reactions take place. These differences explain the relative stability of the circular species, demonstrate how it is prone to react with a target molecule for circle integration and thus supports the notion that the circular form is a biologically significant molecule possibly with a role in intron mobility.

Mechanistic study of base-pairing small regulatory RNAs in bacteria

In all three kingdoms of life, RNA is not only involved in the expression of genetic information, but also carries out extremely diverse cellular functions. This versatility is essentially due to the fact that RNA molecules can exploit the power of base pairing to allow them to fold into a wide variety of structures through which they can perform diverse roles, but also to selectively target and bind to other nucleic acids. This is true in particular for bacterial small regulatory RNAs that act by imperfect base-pairing with target mRNAs, and thereby control their expression through different mechanisms. Here we outline an overview of in vivo and in vitro approaches that are currently used to gain mechanistic insights into how these sRNAs control gene expression in bacteria.

Tertiary structure mapping of the pri-miRNA miR-17~92

The understanding of RNA in regulating gene expression has exploded over the past 15 years. MicroRNAs (miRNAs) have vastly expanded the role of RNA in gene regulation beyond spliceosomal, ribosomal, and messenger RNAs. Approximately one half of miRNAs are polycistronic, where two or more miRNAs are encoded on a single pri-miRNA transcript, termed a miRNA cluster. The six miRNAs of the miR-17~92 cluster are contained within a ~800 nucleotide region within intron 3 of the cl13orf25 ~7 kb pri-miRNA transcript. We recently reported on the tertiary structured domain of miR-17~92 and its role in modulating miRNA biogenesis. The key finding was that the cluster structure explained the differential processing of the miRNA hairpins by Drosha. This work demonstrated the need to consider pri-miRNA tertiary structure in miRNA biogenesis. Since biochemical structure probing is typically performed on relatively short RNAs (≤200 nucleotides), we had to adapt these methodologies for application on large RNAs (~800 nucleotide miR-17~92 pri-miRNA). We present here our adaptation of a protection footprinting method using ribonucleases to probe the structure of the ~800 nucleotide miR-17~92 pri-miRNA. We outline the technical difficulties involved in probing large RNAs and data visualization using denaturing polyacrylamide gel electrophoresis and how we adapted the existing approaches to probe large RNAs. The methodology outlined here is generally applicable to large RNAs including long noncoding RNAs (lncRNA).

Mod-seq: A High-Throughput Method for Probing RNA Secondary Structure

It has become increasingly clear that large RNA molecules, especially long noncoding RNAs, function in almost all gene regulatory processes (Cech & Steitz, 2014). Many large RNAs appear to be structural scaffolds for assembly of important RNA/protein complexes. However, the structures of most large cellular RNA molecules are currently unknown (Hennelly & Sanbonmatsu, 2012). While chemical probing can reveal single-stranded regions of RNA, traditional approaches to identify sites of chemical modification are time consuming. Mod-seq is a high-throughput method used to map chemical modification sites on RNAs of any size, including complex mixtures of RNA. In this protocol, we describe preparation of Mod-seq high-throughput sequencing libraries from chemically modified RNA. We also describe a software package "Mod-seeker," which is a compilation of scripts written in Python, for the analysis of Mod-seq data. Mod-seeker returns statistically significant modification sites, which can then be used to aid in secondary structure prediction.

Structure-specific ribonucleases for MS-based elucidation of higher-order RNA structure

Supported by high-throughput sequencing technologies, structure-specific nucleases are experiencing a renaissance as biochemical probes for genome-wide mapping of nucleic acid structure. This report explores the benefits and pitfalls of the application of Mung bean (Mb) and V1 nuclease, which attack specifically single- and double-stranded regions of nucleic acids, as possible structural probes to be employed in combination with MS detection. Both enzymes were found capable of operating in ammonium-based solutions that are preferred for high-resolution analysis by direct infusion electrospray ionization (ESI). Sequence analysis by tandem mass spectrometry (MS/MS) was performed to confirm mapping assignments and to resolve possible ambiguities arising from the concomitant formation of isobaric products with identical base composition and different sequences. The observed products grouped together into ladder-type series that facilitated their assignment to unique regions of the substrate, but revealed also a certain level of uncertainty in identifying the boundaries between paired and unpaired regions. Various experimental factors that are known to stabilize nucleic acid structure, such as higher ionic strength, presence of Mg(II), etc., increased the accuracy of cleavage information, but did not completely eliminate deviations from expected results. These observations suggest extreme caution in interpreting the results afforded by these types of reagents. Regardless of the analytical platform of choice, the results highlighted the need to repeat probing experiments under the most diverse possible conditions to recognize potential artifacts and to increase the level of confidence in the observed structural information.

Probing RNA folding by hydroxyl radical footprinting

In recent years RNA molecules have emerged as central players in the regulation of gene expression. Many of these noncoding RNAs possess well-defined, complex, three-dimensional structures which are essential for their biological function. In this context, much effort has been devoted to develop computational and experimental techniques for RNA structure determination. Among available experimental tools to investigate the higher-order folding of structured RNAs, hydroxyl radical probing stands as one of the most informative and reliable ones. Hydroxyl radicals are oxidative species that cleave the nucleic acid backbone solely according to the solvent accessibility of individual phosphodiester bonds, with no sequence or secondary structure specificity. Therefore, the cleavage pattern obtained directly reflects the degree of protection/exposure to the solvent of each section of the molecule under inspection, providing valuable information about how these different sections interact together to form the final three-dimensional architecture. In this chapter we describe a robust, accurate and very sensitive hydroxyl radical probing method that can be applied to any structured RNA molecule and is suitable to investigate RNA folding and RNA conformational changes induced by binding of a ligand.

Mapping RNA structure in vitro using nucleobase-specific probes

RNAs have to adopt specific three-dimensional structures to fulfill their biological functions. Therefore exploring RNA structure is of interest to understand RNA-dependent processes. Chemical probing in vitro is a very powerful tool to investigate RNA molecules under a variety of conditions. Among the most frequently used chemical reagents are the nucleobase-specific probes dimethyl sulfate (DMS), 1-cyclohexyl-3-(2-morpholinoethyl) carbodiimide metho-p-toluenesulfonate (CMCT) and β-ethoxy-α-ketobutyraldehyde (kethoxal). These chemical reagents modify nucleotides which are not involved in hydrogen bonding or protected by a ligand, such as proteins or metabolites. Upon performing modification reactions with all three chemicals the accessibility of all four nucleobases can be determined. With this fast and inexpensive method local changes in RNA secondary and tertiary structure, as well as the formation of contacts between RNA and its ligands can be detected independent of the RNA's length.

Chemical probing of RNA in living cells

RNAs need to adopt a specific architecture to exert their task in cells. While significant progress has been made in describing RNA folding landscapes in vitro, understanding intracellular RNA structure formation is still in its infancy. This is in part due to the complex nature of the cellular environment but also to the limited availability of suitable methodologies. To assess the intracellular structure of large RNAs, we recently applied a chemical probing technique and a metal-induced cleavage assay in vivo. These methods are based on the fact that small molecules, like dimethyl sulfate (DMS), or metal ions, such as Pb(2+), penetrate and spread throughout the cell very fast. Hence, these chemicals are able to modify accessible RNA residues or to induce cleavage of the RNA strand in the vicinity of a metal ion in living cells. Mapping of these incidents allows inferring information on the intracellular conformation, metal ion binding sites or ligand-induced structural changes of the respective RNA molecule. Importantly, in vivo chemical probing can be easily adapted to study RNAs in different cell types.

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 centrality of RNA for engineering gene expression

Synthetic biology holds promise as both a framework for rationally engineering biological systems and a way to revolutionize how we fundamentally understand them. Essential to realizing this promise is the development of strategies and tools to reliably and predictably control and characterize sophisticated patterns of gene expression. Here we review the role that RNA can play towards this goal and make a case for why this versatile, designable, and increasingly characterizable molecule is one of the most powerful substrates for engineering gene expression at our disposal. We discuss current natural and synthetic RNA regulators of gene expression acting at key points of control – transcription, mRNA degradation, and translation. We also consider RNA structural probing and computational RNA structure predication tools as a way to study RNA structure and ultimately function. Finally, we discuss how next-generation sequencing methods are being applied to the study of RNA and to the characterization of RNA's many properties throughout the cell.

Multilevel functional and structural defects induced by two pathogenic mitochondrial tRNA mutations

Point mutations in hmtRNAs (human mitochondrial tRNAs) can cause various disorders, such as CPEO (chronic progressive external ophthalmoplegia) and MM (mitochondrial myopathy). Mitochondrial tRNALeu, especially the UUR codon isoacceptor, is recognized as a hot spot for pathogenic mtDNA point mutations. Thus far, 40 mutations have been reported in hmtRNAsLeu. In the present paper, we describe the wide range of effects of two substitutions found in the TΨC arms of two hmtRNAsLeu isoacceptors. The G52A substitution, corresponding to the pathogenic G12315A mutation in tRNALeu(CUN), and G3283A in tRNALeu(UUR) exhibited structural changes in the outer corner of the tRNA shape as shown by RNase probing. These mutations also induced reductions in aminoacylation, 3'-end processing and base modification processes. The main effects of the A57G substitution, corresponding to mutations A12320G in tRNALeu(CUN) and A3288G in tRNALeu(UUR), were observed on the aminoacylation activity and binding to hmEF-Tu (human mitochondrial elongation factor Tu). These observations suggest that the wide range of effects may amplify the deleterious impact on mitochondrial protein synthesis in vivo. The findings also emphasize that an exact understanding of tRNA dysfunction is critical for the future development of therapies for mitochondrial diseases.

Small RNA modules confer different stabilities and interact differently with multiple targets

Bacterial Hfq-associated small regulatory RNAs (sRNAs) parallel animal microRNAs in their ability to control multiple target mRNAs. The small non-coding MicA RNA represses the expression of several genes, including major outer membrane proteins such as ompA, tsx and ecnB. In this study, we have characterised the RNA determinants involved in the stability of MicA and analysed how they influence the expression of its targets. Site-directed mutagenesis was used to construct MicA mutated forms. The 5′linear domain, the structured region with two stem-loops, the A/U-rich sequence or the 3′ poly(U) tail were altered without affecting the overall secondary structure of MicA. The stability and the target regulation abilities of the wild-type and the different mutated forms of MicA were then compared. The 5′ domain impacted MicA stability through an RNase III-mediated pathway. The two stem-loops showed different roles and disruption of stem-loop 2 was the one that mostly affected MicA stability and abundance. Moreover, STEM2 was found to be more important for the in vivo repression of both ompA and ecnB mRNAs while STEM1 was critical for regulation of tsx mRNA levels. The A/U-rich linear sequence is not the only Hfq-binding site present in MicA and the 3′ poly(U) sequence was critical for sRNA stability. PNPase was shown to be an important exoribonuclease involved in sRNA degradation. In addition to the 5′ domain of MicA, the stem-loops and the 3′ poly(U) tail are also shown to affect target-binding. Disruption of the 3′U-rich sequence greatly affects all targets analysed. In conclusion, our results have shown that it is important to understand the “sRNA anatomy” in order to modulate its stability. Furthermore, we have demonstrated that MicA RNA can use different modules to regulate its targets. This knowledge can allow for the engineering of non-coding RNAs that interact differently with multiple targets.

Protein-RNA footprinting: an evolving tool

As more RNA molecules with important cellular functions are discovered, there is a strong need to characterize their structures, functions, and interactions. Chemical and enzymatic footprinting methods are used to map RNA secondary and tertiary structure, to monitor ligand interactions and conformational changes, and in the study of protein-RNA interactions. These methods provide data at single-nucleotide resolution that nicely complements the structural information available from X-ray diffraction, nuclear magnetic resonance spectroscopy (NMR), or cryo-electron microscopy. Footprinting methods also complement the dynamic information derived from single-molecule Förster resonance energy transfer. RNA footprinting tools have been used for decades, but we have recently seen spectacular advances, for instance, the use in combination with massive parallel sequencing techniques. Large libraries of RNA molecules (small or large in size) can now be probed in high-throughput manner when RNA footprinting methods are combined with fluorescent probe technologies and automation. In this article, after a brief historical overview, we summarize recent advances in RNA-protein footprinting methodologies that now integrate tools for massive parallel analysis.

Selective 2'-hydroxyl acylation analyzed by protection from exoribonuclease (RNase-detected SHAPE) for direct analysis of covalent adducts and of nucleotide flexibility in RNA

RNA SHAPE chemistry yields quantitative, single-nucleotide resolution structural information based on the reaction of the 2'-hydroxyl group of conformationally flexible nucleotides with electrophilic SHAPE reagents. However, SHAPE technology has been limited by the requirement that sites of RNA modification be detected by primer extension. Primer extension results in loss of information at both the 5' and 3' ends of an RNA and requires multiple experimental steps. Here we describe RNase-detected SHAPE that uses a processive, 3'→5' exoribonuclease, RNase R, to detect covalent adducts in 5'-end-labeled RNA in a one-tube experiment. RNase R degrades RNA but stops quantitatively three and four nucleotides 3' of a nucleotide containing a covalent adduct at the ribose 2'-hydroxyl or the pairing face of a nucleobase, respectively. We illustrate this technology by characterizing ligand-induced folding for the aptamer domain of the Escherichia coli thiamine pyrophosphate riboswitch RNA. RNase-detected SHAPE is a facile, two-day approach that can be used to analyze diverse covalent adducts in any RNA molecule, including short RNAs not amenable to analysis by primer extension and RNAs with functionally important structures at their 5' or 3' ends.

A powerful approach for the selection of 2-aminopurine substitution sites to investigate RNA folding

A precise tertiary structure must be adopted to allow the function of many RNAs in cells. Accordingly, increasing resources have been devoted to the elucidation of RNA structures and the folding of RNAs. 2-Aminopurine (2AP), a fluorescent nucleobase analogue, can be substituted in strategic positions of DNA or RNA molecules to act as site-specific probe to monitor folding and folding dynamics of nucleic acids. Recent studies further demonstrated the potential of 2AP modifications in the assessment of folding kinetics during ligand-induced secondary and tertiary RNA structure rearrangements. However, an efficient way to unambiguously identify reliable positions for 2AP sensors is as yet unavailable and would represent a major asset, especially in the absence of crystallographic or NMR structural data for a target molecule. We report evidence of a novel and direct correlation between the 2'-OH flexibility of nucleotides, observed by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) probing and the fluorescence response following nucleotide substitutions by 2AP. This correlation leads to a straightforward method, using SHAPE probing with benzoyl cyanide, to select appropriate nucleotide sites for 2AP substitution. This clear correlation is presented for three model RNAs of biological significance: the SAM-II, adenine (addA), and preQ(1) class II (preQ(1)cII) riboswitches.

Cap-assisted internal initiation of translation of histone H4

In eukaryotes, a crucial step of translation initiation is the binding of the multifactor complex eIF4F to the 5' end of the mRNA, a prerequisite to recruitment of the activated small ribosomal 43S particle. Histone H4 mRNAs have short 5'UTRs, which do not conform to the conventional scanning-initiation model. Here we show that the ORF of histone mRNA contains two structural elements critical for translation initiation. One of the two structures binds eIF4E without the need of the cap. Ribosomal 43S particles become tethered to this site and directly loaded in the vicinity of the AUG. The other structure, 19 nucleotides downstream of the initiation codon, forms a three-way helix junction, which sequesters the m(7)G cap. This element facilitates direct positioning of the ribosome on the cognate start codon. This unusual translation initiation mode might be considered as a hybrid mechanism between the canonical and the IRES-driven translation initiation process.

The use and abuse of -omes

The diverse fields of Omics research share a common logical structure combining a cataloging effort for a particular class of molecules or interactions, the underlying -ome, and a quantitative aspect attempting to record spatiotemporal patterns of concentration, expression, or variation. Consequently, these fields also share a common set of difficulties and limitations. In spite of the great success stories of Omics projects over the last decade, much remains to be understood not only at the technological, but also at the conceptual level. Here, we focus on the dark corners of Omics research, where the problems, limitations, conceptual difficulties, and lack of knowledge are hidden.

RNA interactions

Noncoding RNAs form an indispensible component of the cellular information processing networks, a role that crucially depends on the specificity of their interactions among each other as well as with DNA and protein. Patterns of intramolecular and intermolecular base pairs govern most RNA interactions. Specific base pairs dominate the structure formation of nucleic acids. Only little details distinguish intramolecular secondary structures from those cofolding molecules. RNA-protein interactions, on the other hand, are strongly dependent on the RNA structure as well since the sequence content of helical regions is largely unreadable, so that sequence specificity is mostly restricted to unpaired loop regions. Conservation of both sequence and structure thus this can give indications of the functioning of the diversity of ncRNAs.

DEAD-box protein facilitated RNA folding in vivo

In yeast mitochondria the DEAD-box helicase Mss116p is essential for respiratory growth by acting as group I and group II intron splicing factor. Here we provide the first structure-based insights into how Mss116p assists RNA folding in vivo. Employing an in vivo chemical probing technique, we mapped the structure of the ai5γ group II intron in different genetic backgrounds to characterize its intracellular fold. While the intron adopts the native conformation in the wt yeast strain, we found that the intron is able to form most of its secondary structure, but lacks its tertiary fold in the absence of Mss116p. This suggests that ai5γ is largely unfolded in the mss116-knockout strain and requires the protein at an early step of folding. Notably, in this unfolded state misfolded substructures have not been observed. As most of the protein-induced conformational changes are located within domain D1, Mss116p appears to facilitate the formation of this largest domain, which is the scaffold for docking of other intron domains. These findings suggest that Mss116p assists the ordered assembly of the ai5γ intron in vivo.

Computational discovery of folded RNA domains in genomes and in vitro selected libraries

Structured functional RNAs are conserved on the level of secondary and tertiary structure, rather than at sequence level, and so traditional sequence-based searches often fail to identify them. Structure-based searches are increasingly used to discover known RNA motifs in sequence databases. We describe the application of the program RNABOB, which performs such searches by allowing the user to define a desired motif's sequence, paired and spacer elements and then scans a sequence file for regions capable of assuming the prescribed fold. Structure descriptors of stem-loops, internal loops, three-way junctions, kissing loops, and the hammerhead and hepatitis delta virus ribozymes are shown as examples of implementation of structure-based searches.

Exploring RNA structure by integrative molecular modelling

RNA molecular modelling is adequate to rapidly tackle the structure of RNA molecules. With new structured RNAs constituting a central class of cellular regulators discovered every year, the need for swift and reliable modelling methods is more crucial than ever. The pragmatic method based on interactive all-atom molecular modelling relies on the observation that specific structural motifs are recurrently found in RNA sequences. Once identified by a combination of comparative sequence analysis and biochemical data, the motifs composing the secondary structure of a given RNA can be extruded in three dimensions (3D) and used as building blocks assembled manually during a bioinformatic interactive process. Comparing the models to the corresponding crystal structures has validated the method as being powerful to predict the RNA topology and architecture while being less accurate regarding the prediction of base-base interactions. These aspects as well as the necessary steps towards automation will be discussed.

The Crc global regulator binds to an unpaired A-rich motif at the Pseudomonas putida alkS mRNA coding sequence and inhibits translation initiation

Crc is a key global translational regulator in Pseudomonads that orchestrates the hierarchy of induction of several catabolic pathways for amino acids, sugars, hydrocarbons or aromatic compounds. In the presence of amino acids, which are preferred carbon sources, Crc inhibits translation of the Pseudomonas putida alkS and benR mRNAs, which code for transcriptional regulators of genes required to assimilate alkanes (hydrocarbons) and benzoate (an aromatic compound), respectively. Crc binds to the 5′-end of these mRNAs, but the sequence and/or structure recognized, and the way in which it inhibits translation, were unknown. We have determined the secondary structure of the alkS mRNA 5′-end through its sensitivity to several ribonucleases and chemical reagents. Footprinting and band-shift assays using variant alkS mRNAs have shown that Crc specifically binds to a short unpaired A-rich sequence located adjacent to the alkS AUG start codon. This interaction is stable enough to prevent formation of the translational initiation complex. A similar Crc-binding site was localized at benR mRNA, upstream of the Shine–Dalgarno sequence. This allowed predicting binding sites at other Crc-regulated genes, deriving a consensus sequence that will help to validate new Crc targets and to discriminate between direct and indirect effects of this regulator.

mRNA secondary structures fold sequentially but exchange rapidly in vivo

RNAs adopt defined structures to perform biological activities, and conformational transitions among alternative structures are critical to virtually all RNA-mediated processes ranging from metabolite-activation of bacterial riboswitches to pre-mRNA splicing and viral replication in eukaryotes. Mechanistic analysis of an RNA folding reaction in a biological context is challenging because many steps usually intervene between assembly of a functional RNA structure and execution of a biological function. We developed a system to probe mechanisms of secondary structure folding and exchange directly in vivo using self-cleavage to monitor competition between mutually exclusive structures that promote or inhibit ribozyme assembly. In previous work, upstream structures were more effective than downstream structures in blocking ribozyme assembly during transcription in vitro, consistent with a sequential folding mechanism. However, upstream and downstream structures blocked ribozyme assembly equally well in vivo, suggesting that intracellular folding outcomes reflect thermodynamic equilibration or that annealing of contiguous sequences is favored kinetically. We have extended these studies to learn when, if ever, thermodynamic stability becomes an impediment to rapid equilibration among alternative RNA structures in vivo. We find that a narrow thermodynamic threshold determines whether kinetics or thermodynamics govern RNA folding outcomes in vivo. mRNA secondary structures fold sequentially in vivo, but exchange between adjacent secondary structures is much faster in vivo than it is in vitro. Previous work showed that simple base-paired RNA helices dissociate at similar rates in vivo and in vitro so exchange between adjacent structures must occur through a different mechanism, one that likely involves facilitation of branch migration by proteins associated with nascent transcripts.

Analysis of small RNA in fission yeast; centromeric siRNAs are potentially generated through a structured RNA

The formation of heterochromatin at the centromeres in fission yeast depends on transcription of the outer repeats. These transcripts are processed into siRNAs that target homologous loci for heterochromatin formation. Here, high throughput sequencing of small RNA provides a comprehensive analysis of centromere-derived small RNAs. We found that the centromeric small RNAs are Dcr1 dependent, carry 5'-monophosphates and are associated with Ago1. The majority of centromeric small RNAs originate from two remarkably well-conserved sequences that are present in all centromeres. The high degree of similarity suggests that this non-coding sequence in itself may be of importance. Consistent with this, secondary structure-probing experiments indicate that this centromeric RNA is partially double-stranded and is processed by Dicer in vitro. We further demonstrate the existence of small centromeric RNA in rdp1Delta cells. Our data suggest a pathway for siRNA generation that is distinct from the well-documented model involving RITS/RDRC. We propose that primary transcripts fold into hairpin-like structures that may be processed by Dcr1 into siRNAs, and that these siRNAs may initiate heterochromatin formation independent of RDRC activity.

A new type of IRES within gag coding region recruits three initiation complexes on HIV-2 genomic RNA

Genomic RNA of primate lentiviruses serves both as an mRNA that encodes Gag and Gag-Pol polyproteins and as a propagated genome. Translation of this RNA is initiated by standard cap dependant mechanism or by internal entry of the ribosome. Two regions of the genomic RNA are able to attract initiation complexes, the 5′ untranslated region and the gag coding region itself. Relying on probing data and a phylogenetic study, we have modelled the secondary structure of HIV-1, HIV-2 and SIV_Mac coding region. This approach brings to light conserved secondary-structure elements that were shown by mutations to be required for internal entry of the ribosome. No structural homologies with other described viral or cellular IRES can be identified and lentiviral IRESes show many peculiar properties. Most notably, the IRES present in HIV-2 gag coding region is endowed with the unique ability to recruit up to three initiation complexes on a single RNA molecule. The structural and functional properties of gag coding sequence define a new type of IRES. Although its precise role is unknown, the conservation of the IRES among fast evolving lentiviruses suggests an important physiological role.

Probing RNA structure within living cells

RNA folding is the most fundamental process underlying RNA function. RNA structure and associated folding paradigms have been intensively studied in vitro. However, in vivo RNA structure formation has only been explored to a limited extent. To determine the influence of the cellular environment, which differs significantly from the in vitro refolding conditions, on RNA architecture, we have applied a chemical probing technique to assess the structure of catalytic RNAs in living cells. This method is based on the fact that chemicals like dimethyl sulfate readily penetrate cells and modify specific atoms of RNA bases (N1-A, N3-C), provided that these positions are solvent accessible. By mapping the modified residues, one gains substantial information on the architecture of the target RNA on the secondary and tertiary structure level. This method also allows exploration of interactions of the target RNA with ligands such as proteins, metabolites, or other RNA molecules and associated conformational changes. In brief, in vivo chemical probing is a powerful tool to investigate RNA structure in its natural environment and can be easily adapted to study RNAs in different cell types.

Tertiary network in mammalian mitochondrial tRNAAsp revealed by solution probing and phylogeny

Primary and secondary structures of mammalian mitochondrial (mt) tRNAs are divergent from canonical tRNA structures due to highly skewed nucleotide content and large size variability of D- and T-loops. The nonconservation of nucleotides involved in the expected network of tertiary interactions calls into question the rules governing a functional L-shaped three-dimensional (3D) structure. Here, we report the solution structure of human mt-tRNA^Asp in its native post-transcriptionally modified form and as an in vitro transcript. Probing performed with nuclease S1, ribonuclease V1, dimethylsulfate, diethylpyrocarbonate and lead, revealed several secondary structures for the in vitro transcribed mt-tRNA^Asp including predominantly the cloverleaf. On the contrary, the native tRNA^Asp folds into a single cloverleaf structure, highlighting the contribution of the four newly identified post-transcriptional modifications to correct folding. Reactivities of nucleotides and phosphodiester bonds in the native tRNA favor existence of a full set of six classical tertiary interactions between the D-domain and the variable region, forming the core of the 3D structure. Reactivities of D- and T-loop nucleotides support an absence of interactions between these domains. According to multiple sequence alignments and search for conservation of Leontis–Westhof interactions, the tertiary network core building rules apply to all tRNA^Asp from mammalian mitochondria.

The CUGBP2 splicing factor regulates an ensemble of branchpoints from perimeter binding sites with implications for autoregulation

Alternative pre-mRNA splicing adjusts the transcriptional output of the genome by generating related mRNAs from a single primary transcript, thereby expanding protein diversity. A fundamental unanswered question is how splicing factors achieve specificity in the selection of target substrates despite the recognition of information-poor sequence motifs. The CUGBP2 splicing regulator plays a key role in the brain region-specific silencing of the NI exon of the NMDA R1 receptor. However, the sequence motifs utilized by this factor for specific target exon selection and its role in splicing silencing are not understood. Here, we use chemical modification footprinting to map the contact sites of CUGBP2 to GU-rich motifs closely positioned at the boundaries of the branch sites of the NI exon, and we demonstrate a mechanistic role for this specific arrangement of motifs for the regulation of branchpoint formation. General support for a branch site-perimeter–binding model is indicated by the identification of a group of novel target exons with a similar configuration of motifs that are silenced by CUGBP2. These results reveal an autoregulatory role for CUGBP2 as indicated by its direct interaction with functionally significant RNA motifs surrounding the branch sites upstream of exon 6 of the CUGBP2 transcript itself. The perimeter-binding model explains how CUGBP2 can effectively embrace the branch site region to achieve the specificity needed for the selection of exon targets and the fine-tuning of alternative splicing patterns.

Alternative splicing is a precisely controlled process that determines whether an exon will be included or skipped in the mature mRNA transcript. Factors that control alternative splicing bind to RNA sequence motifs in the exon or flanking introns and guide tissue and developmental specific splicing events. CUGBP2 is a dual functional regulator of alternative splicing that can cause inclusion or skipping of a target exon, depending on the context of its binding motifs. Previously, the mechanisms of regulation by this protein and the positional significance of its target motifs have not been characterized. In this study, the authors dissected the mechanism of exon skipping by CUGBP2 and demonstrate that a specific configuration of motifs at the perimeters of a functional reference point are intimately involved in this event. Furthermore, this mechanism of regulation is shown to have general significance because novel CUGBP2 target exons contain a similar arrangement of motifs. The most interesting of this group is an exon within the CUGBP2 transcript itself. This study underscores the importance of a functional reference point in the specificity of regulation by an alternative splicing factor and reveals a novel autoregulatory role for CUGBP2.

Noncovalent probes for the investigation of structure and dynamics of protein-nucleic acid assemblies: the case of NC-mediated dimerization of genomic RNA in HIV-1

The nature of specific RNA-RNA and protein-RNA interactions involved in the process of genome dimerization and isomerization in HIV-1, which is mediated in vitro by stemloop 1 (SL1) of the packaging signal and by the nucleocapsid (NC) domain of the viral Gag polyprotein, was investigated by using archetypical nucleic acid ligands as noncovalent probes. Small-molecule ligands make contact with their target substrates through complex combinations of H-bonds, salt bridges, and hydrophobic interactions. Therefore, their binding patterns assessed by electrospray ionization mass spectrometry can provide valuable insights into the factors determining specific recognition between species involved in biopolymer assemblies. In the case of SL1, dimerization and isomerization create unique structural features capable of sustaining stable interactions with classic nucleic acid ligands. The binding modes exhibited by intercalators and minor groove binders were adversely affected by the significant distortion of the duplex formed by palindrome annealing in the kissing-loop (KL) dimer, whereas the modes observed for the corresponding extended duplex (ED) confirmed a more regular helical structure. Consistent with the ability to establish electrostatic interactions with highly negative pockets typical of helix anomalies, polycationic aminoglycosides bound to the stem-bulge motif conserved in all SL1 conformers, to the unpaired nucleotides located at the hinge between kissing hairpins in KL, and to the exposed bases flanking the palindrome duplex in ED. The patterns afforded by intercalators and minor groove binders did not display detectable variations when the corresponding NC-SL1 complexes were submitted to probing. In contrast, aminoglycosides displayed the ability to compete with the protein for overlapping sites, producing opposite effects on the isomerization process. Indeed, displacing NC from the stem-bulges of the KL dimer induced inhibition of stem melting and decreased the efficiency of isomerization. Competition for the hinge region, instead, eliminated the NC stabilization of a grip motif formed by nucleobases of opposite strands, thus facilitating the strand-exchange required for isomerization. These noncovalent probes provided further evidence that the structural context of the actual binding sites has significant influence on the chaperone activities of NC, which should be taken in account when developing potential drug candidates aimed at disrupting genome dimerization and isomerization in HIV-1.

Rescue of internal initiation of translation by RNA complementation provides evidence for a distribution of functions between individual IRES domains

Picornavirus RNAs initiate translation using an internal ribosome entry site (IRES)-dependent mechanism. The IRES element of foot-and-mouth disease virus (FMDV) is organized in domains, being different from each other in RNA structure and RNA-protein interaction. Wild-type transcripts provided in trans rescue defective FMDV IRES mutants. Complementation, however, was partial since translation efficiency of the mutant RNAs was up to 10% of the wild type IRES. We report here that mutations diminishing the RNA-RNA interaction capacity induced a decrease in IRES rescue. On the other hand, IRES transcripts bearing mutations that reorganize the RNA structure of the apical region of central domain, although weakly, complement defective IRES that are unable to interact with the initiation factor eIF4G in a separate domain. Together, these results suggest that IRES rescue may involve RNA-mediated contacts between defective elements, each carrying a defect in a separate domain but having the complementing one with the appropriate structural orientation and/or ribonucleoprotein composition. Our observations further support the essential role of the central domain of the FMDV IRES during protein synthesis and demonstrate that there is a division of functions between the IRES domains.

A novel function for fragile X mental retardation protein in translational activation

Fragile X syndrome, the most frequent form of inherited mental retardation, is due to the absence of Fragile X Mental Retardation Protein (FMRP), an RNA-binding protein involved in several steps of RNA metabolism. To date, two RNA motifs have been found to mediate FMRP/RNA interaction, the G-quartet and the "kissing complex," which both induce translational repression in the presence of FMRP. We show here a new role for FMRP as a positive modulator of translation. FMRP specifically binds Superoxide Dismutase 1 (Sod1) mRNA with high affinity through a novel RNA motif, SoSLIP (Sod1 mRNA Stem Loops Interacting with FMRP), which is folded as three independent stem-loop structures. FMRP induces a structural modification of the SoSLIP motif upon its interaction with it. SoSLIP also behaves as a translational activator whose action is potentiated by the interaction with FMRP. The absence of FMRP results in decreased expression of Sod1. Because it has been observed that brain metabolism of FMR1 null mice is more sensitive to oxidative stress, we propose that the deregulation of Sod1 expression may be at the basis of several traits of the physiopathology of the Fragile X syndrome, such as anxiety, sleep troubles, and autism.