Probing the structure of RNAs in solution

SEISMICgraph: a web-based tool for RNA structure data visualization

In recent years, RNA has been increasingly recognized for its essential roles in biology, functioning not only as a carrier of genetic information but also as a dynamic regulator of gene expression through its interactions with other RNAs, proteins, and itself. Advances in chemical probing techniques have significantly enhanced our ability to identify RNA secondary structures and understand their regulatory roles. These developments, alongside improvements in experimental design and data processing, have greatly increased the resolution and throughput of structural analyses. Here, we introduce SEISMICgraph, a web-based tool designed to support RNA structure research by offering data visualization and analysis capabilities for a variety of chemical probing modalities. SEISMICgraph enables simultaneous comparison of data across different sequences and experimental conditions through a user-friendly interface that requires no programming expertise. We demonstrate its utility by investigating known and putative riboswitches and exploring how RNA modifications influence their structure and binding. SEISMICgraph's ability to rapidly visualize adenine-dependent structural changes and assess the impact of pseudouridylation on these transitions provides novel insights and establishes a roadmap for numerous future applications.

Targeting sgRNA N secondary structure as a way of inhibiting SARS-CoV-2 replication

SARS-CoV-2 is a betacoronavirus that causes COVID-19, a global pandemic that has resulted in many infections, deaths, and socio-economic challenges. The virus has a large positive-sense, single-stranded RNA genome of ∼30 kb, which produces subgenomic RNAs (sgRNAs) through discontinuous transcription. The most abundant sgRNA is sgRNA N, which encodes the nucleocapsid (N) protein. In this study, we probed the secondary structure of sgRNA N and a shorter model without a 3' UTR in vitro, using the SHAPE (selective 2'-hydroxyl acylation analyzed by a primer extension) method and chemical mapping with dimethyl sulfate and 1-cyclohexyl-(2-morpholinoethyl) carbodiimide metho-p-toluene sulfonate. We revealed the secondary structure of sgRNA N and its shorter variant for the first time and compared them with the genomic RNA N structure. Based on the structural information, we designed gapmers, siRNAs and antisense oligonucleotides (ASOs) to target the N protein coding region of sgRNA N. We also generated eukaryotic expression vectors containing the complete sequence of sgRNA N and used them to screen for new SARS-CoV-2 gene N expression inhibitors. Our study provides novel insights into the structure and function of sgRNA N and potential therapeutic tools against SARS-CoV-2.

A new reagent for in vivo structure probing of RNA G and U residues that improves RNA structure prediction alone and combined with DMS

A key to understanding the roles of RNA in regulating gene expression is knowing their structures in vivo. One way to obtain this information is through probing the structures of RNA with chemicals. To probe RNA structure directly in cells, membrane-permeable reagents that modify the Watson-Crick (WC) face of unpaired nucleotides can be used. Although dimethyl sulfate (DMS) has led to substantial insight into RNA structure, it has limited nucleotide specificity in vivo, with WC face reactivity only at adenine (A) and cytosine (C) at neutral pH. The reagent 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide (EDC) was recently shown to modify the WC face of guanine (G) and uracil (U). Although useful at lower concentrations in experiments that measure chemical modifications by reverse transcription (RT) stops, at higher concentrations necessary for detection by mutational profiling (MaP), EDC treatment leads to degradation of RNA. Here, we demonstrate EDC-stimulated degradation of RNA in Gram-negative and Gram-positive bacteria. In an attempt to overcome these limitations, we developed a new carbodiimide reagent, 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide methiodide (ETC), which we show specifically modifies unpaired Gs and Us in vivo without substantial degradation of RNA. We establish ETC as a probe for MaP and optimize the RT conditions and computational analysis in Escherichia coli Importantly, we demonstrate the utility of ETC as a probe for improving RNA structure prediction both alone and with DMS.

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

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

RNA-Based Sensor Systems for Affordable Diagnostics in the Age of Pandemics

In the era of the COVID-19 pandemic, the significance of point-of-care (POC) diagnostic tools has become increasingly vital, driven by the need for quick and precise virus identification. RNA-based sensors, particularly toehold sensors, have emerged as promising candidates for POC detection systems due to their selectivity and sensitivity. Toehold sensors operate by employing an RNA switch that changes the conformation when it binds to a target RNA molecule, resulting in a detectable signal. This review focuses on the development and deployment of RNA-based sensors for POC viral RNA detection with a particular emphasis on toehold sensors. The benefits and limits of toehold sensors are explored, and obstacles and future directions for improving their performance within POC detection systems are presented. The use of RNA-based sensors as a technology for rapid and sensitive detection of viral RNA holds great potential for effectively managing (dealing/coping) with present and future pandemics in resource-constrained settings.

Structural polymorphism of the nucleic acids in pentanucleotide repeats associated with the neurological disorder CANVAS

Short tandem repeats are inherently unstable during DNA replication depending on repeat length, and the expansion of the repeat length in the human genome is responsible for repeat expansion disorders. Pentanucleotide AAGGG and ACAGG repeat expansions in intron 2 of the gene encoding replication factor C subunit 1 (RFC1) cause cerebellar ataxia, neuropathy, vestibular areflexia syndrome (CANVAS) and other phenotypes of late-onset cerebellar ataxia. Herein, we reveal the structural polymorphism of the RFC1 repeats associated with CANVAS in vitro. Single-stranded AAGGG repeat DNA formed a hybrid-type G-quadruplex, whereas its RNA formed a parallel-type G-quadruplex with three layers. The RNA of the ACAGG repeat formed hairpin structure comprising C-G and G-C base pairs with A:A and GA:AG mismatched repeats. Furthermore, both pathogenic repeat RNAs formed more rigid structures than those of the nonpathogenic repeat RNAs. These findings provide novel insights into the structural polymorphism of the RFC1 repeats, which may be closely related to the disease mechanism of CANVAS.

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

RNA Secondary Structure Modeling Following the IPANEMAP Workflow

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

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

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

RNA structure determination: From 2D to 3D

RNA molecules serve a wide range of functions that are closely linked to their structures. The basic structural units of RNA consist of single- and double-stranded regions. In order to carry out advanced functions such as catalysis and ligand binding, certain types of RNAs can adopt higher-order structures. The analysis of RNA structures has progressed alongside advancements in structural biology techniques, but it comes with its own set of challenges and corresponding solutions. In this review, we will discuss recent advances in RNA structure analysis techniques, including structural probing methods, X-ray crystallography, nuclear magnetic resonance, cryo-electron microscopy, and small-angle X-ray scattering. Often, a combination of multiple techniques is employed for the integrated analysis of RNA structures. We also survey important RNA structures that have been recently determined using various techniques.

RNA Secondary Structure Analysis Using RNAstructure

RNAstructure is a user-friendly program for the prediction and analysis of RNA secondary structure. It is available as a web server, a program with a graphical user interface, or a set of command line tools. The programs are available for Microsoft Windows, macOS, or Linux. This article provides protocols for prediction of RNA secondary structure (using the web server, the graphical user interface, or the command line) and high-affinity oligonucleotide binding sites to a structured RNA target (using the graphical user interface). © 2023 Wiley Periodicals LLC. Basic Protocol 1: Predicting RNA secondary structure using the RNAstructure web server Alternate Protocol 1: Predicting secondary structure and base pair probabilities using the RNAstructure graphical user interface Alternate Protocol 2: Predicting secondary structure and base pair probabilities using the RNAstructure command line interface Basic Protocol 2: Predicting binding affinities of oligonucleotides complementary to an RNA target using OligoWalk.

Nano-DMS-MaP allows isoform-specific RNA structure determination

Genome-wide measurements of RNA structure can be obtained using reagents that react with unpaired bases, leading to adducts that can be identified by mutational profiling on next-generation sequencing machines. One drawback of these experiments is that short sequencing reads can rarely be mapped to specific transcript isoforms. Consequently, information is acquired as a population average in regions that are shared between transcripts, thus blurring the underlying structural landscape. Here, we present nanopore dimethylsulfate mutational profiling (Nano-DMS-MaP)-a method that exploits long-read sequencing to provide isoform-resolved structural information of highly similar RNA molecules. We demonstrate the value of Nano-DMS-MaP by resolving the complex structural landscape of human immunodeficiency virus-1 transcripts in infected cells. We show that unspliced and spliced transcripts have distinct structures at the packaging site within the common 5' untranslated region, likely explaining why spliced viral RNAs are excluded from viral particles. Thus, Nano-DMS-MaP is a straightforward method to resolve biologically important transcript-specific RNA structures that were previously hidden in short-read ensemble analyses.

Rock, scissors, paper: How RNA structure informs function

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

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

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

Strategies for targeting RNA with small molecule drugs

INTRODUCTION

Historically, therapeutic treatment of disease has been restricted to targeting proteins. Of the approximately 20,000 translated human proteins, approximately 1600 are associated with diseases. Strikingly, less than 15% of disease-associated proteins are predicted or known to be 'druggable.' While the concept and narrative of protein druggability continue to evolve with the development of novel technological and pharmacological advances, most of the human proteome remains undrugged. Recent genomic studies indicate that less than 2% of the human genome encodes for proteins, and while as much as 75% of the genome is transcribed, RNA has largely been ignored as a druggable target for therapeutic interventions.

AREAS COVERED

This review delineates the theory and techniques involved in the development of small molecule inhibitors of RNAs from brute force, high-throughput screening technologies to de novo molecular design using computational machine and deep learning. We will also highlight the potential pitfalls and limitations of targeting RNA with small molecules.

EXPERT OPINION

Although significant advances have recently been made in developing systems to identify small molecule inhibitors of RNAs, many challenges remain. Focusing on RNA structure and ligand binding sites may help bring drugging RNA in line with traditional protein drug targeting.

RNA Chemical Labeling with Site-Specific, Relative Quantification by Mass Spectrometry for the Structural Study of a Neomycin-Sensing Riboswitch Aptamer Domain

High-resolution mass spectrometry was used for the label-free, direct localization and relative quantification of CMC+ -modifications of a neomycin-sensing riboswitch aptamer domain in the absence and presence of the aminoglycoside ligands neomycin B, ribostamycin, and paromomycin. The chemical probing and MS data for the free riboswitch show high exposure to solvent of the uridine nucleobases U7, U8, U13, U14, U18 as part of the proposed internal and apical loops, but those of U10 and U21 as part of the proposed internal loop were found to be far less exposed than expected. Thus, our data are in better agreement with the proposed secondary structure of the riboswitch in complexes with aminoglycosides than with that of free RNA. For the riboswitch in complexes with neomycin B, ribostamycin, and paromomycin, we found highly similar CMC+ -modification patterns and excellent agreement with previous NMR studies. Differences between the chemical probing and MS data in the absence and presence of the aminoglycoside ligands were quantitative rather than qualitative (i. e., the same nucleobases were labeled, but to different extents) and can be rationalized by stabilization of both the proposed bulge and the apical loop by aminoglycoside binding. Our study shows that chemical probing and mass spectrometry can provide important structural information and complement other techniques such as NMR spectroscopy.

Genome-wide analysis of the in vivo tRNA structurome reveals RNA structural and modification dynamics under heat stress

RNA structure plays roles in myriad cellular events including transcription, translation, and RNA processing. Genome-wide analyses of RNA secondary structure in vivo by chemical probing have revealed critical structural features of mRNAs and long ncRNAs. Here, we examine the in vivo secondary structure of a small RNA class, tRNAs. Study of tRNA structure is challenging because tRNAs are heavily modified and strongly structured. We introduce "tRNA structure-seq," a new workflow that accurately determines in vivo secondary structures of tRNA. The workflow combines dimethyl sulfate (DMS) probing, ultra-processive RT, and mutational profiling (MaP), which provides mutations opposite DMS and natural modifications thereby allowing multiple modifications to be identified in a single read. We applied tRNA structure-seq to E. coli under control and stress conditions. A leading folding algorithm predicts E. coli tRNA structures with only ∼80% average accuracy from sequence alone. Strikingly, tRNA structure-seq, by providing experimental restraints, improves structure prediction under in vivo conditions to ∼95% accuracy, with more than 14 tRNAs predicted completely correctly. tRNA structure-seq also quantifies the relative levels of tRNAs and their natural modifications at single nucleotide resolution, as validated by LC-MS/MS. Our application of tRNA structure-seq yields insights into tRNA structure in living cells, revealing that it is not immutable but has dynamics, with partial unfolding of secondary and tertiary tRNA structure under heat stress that is correlated with a loss of tRNA abundance. This method is applicable to other small RNAs, including those with natural modifications and highly structured regions.

UFold: fast and accurate RNA secondary structure prediction with deep learning

For many RNA molecules, the secondary structure is essential for the correct function of the RNA. Predicting RNA secondary structure from nucleotide sequences is a long-standing problem in genomics, but the prediction performance has reached a plateau over time. Traditional RNA secondary structure prediction algorithms are primarily based on thermodynamic models through free energy minimization, which imposes strong prior assumptions and is slow to run. Here, we propose a deep learning-based method, called UFold, for RNA secondary structure prediction, trained directly on annotated data and base-pairing rules. UFold proposes a novel image-like representation of RNA sequences, which can be efficiently processed by Fully Convolutional Networks (FCNs). We benchmark the performance of UFold on both within- and cross-family RNA datasets. It significantly outperforms previous methods on within-family datasets, while achieving a similar performance as the traditional methods when trained and tested on distinct RNA families. UFold is also able to predict pseudoknots accurately. Its prediction is fast with an inference time of about 160 ms per sequence up to 1500 bp in length. An online web server running UFold is available at https://ufold.ics.uci.edu. Code is available at https://github.com/uci-cbcl/UFold.

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.

Universal and strain specific structure features of segment 8 genomic RNA of influenza A virus-application of 4-thiouridine photocrosslinking

RNA structure in the influenza A virus (IAV) has been the focus of several studies that have shown connections between conserved secondary structure motifs and their biological function in the virus replication cycle. Questions have arisen on how to best recognize and understand the pandemic properties of IAV strains from an RNA perspective, but determination of the RNA secondary structure has been challenging. Herein, we used chemical mapping to determine the secondary structure of segment 8 viral RNA (vRNA) of the pandemic A/California/04/2009 (H1N1) strain of IAV. Additionally, this long, naturally occurring RNA served as a model to evaluate RNA mapping with 4-thiouridine (4sU) crosslinking. We explored 4-thiouridine as a probe of nucleotides in close proximity, through its incorporation into newly transcribed RNA and subsequent photoactivation. RNA secondary structural features both universal to type A strains and unique to the A/California/04/2009 (H1N1) strain were recognized. 4sU mapping confirmed and facilitated RNA structure prediction, according to several rules: 4sU photocross-linking forms efficiently in the double-stranded region of RNA with some flexibility, in the ends of helices, and across bulges and loops when their structural mobility is permitted. This method highlighted three-dimensional properties of segment 8 vRNA secondary structure motifs and allowed to propose several long-range three-dimensional interactions. 4sU mapping combined with chemical mapping and bioinformatic analysis could be used to enhance the RNA structure determination as well as recognition of target regions for antisense strategies or viral RNA detection.

Chemical and Enzymatic Probing of Viral RNAs: From Infancy to Maturity and Beyond

RNA molecules are key players in a variety of biological events, and this is particularly true for viral RNAs. To better understand the replication of those pathogens and try to block them, special attention has been paid to the structure of their RNAs. Methods to probe RNA structures have been developed since the 1960s; even if they have evolved over the years, they are still in use today and provide useful information on the folding of RNA molecules, including viral RNAs. The aim of this review is to offer a historical perspective on the structural probing methods used to decipher RNA structures before the development of the selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) methodology and to show how they have influenced the current probing techniques. Actually, these technological breakthroughs, which involved advanced detection methods, were made possible thanks to the development of next-generation sequencing (NGS) but also to the previous works accumulated in the field of structural RNA biology. Finally, we will also discuss how high-throughput SHAPE (hSHAPE) paved the way for the development of sophisticated RNA structural techniques.

General Principles for the Detection of Modified Nucleotides in RNA by Specific Reagents

Epitranscriptomics heavily rely on chemical reagents for the detection, quantification, and localization of modified nucleotides in transcriptomes. Recent years have seen a surge in mapping methods that use innovative and rediscovered organic chemistry in high throughput approaches. While this has brought about a leap of progress in this young field, it has also become clear that the different chemistries feature variegated specificity and selectivity. The associated error rates, e.g., in terms of false positives and false negatives, are in large part inherent to the chemistry employed. This means that even assuming technically perfect execution, the interpretation of mapping results issuing from the application of such chemistries are limited by intrinsic features of chemical reactivity. An important but often ignored fact is that the huge stochiometric excess of unmodified over-modified nucleotides is not inert to any of the reagents employed. Consequently, any reaction aimed at chemical discrimination of modified versus unmodified nucleotides has optimal conditions for selectivity that are ultimately anchored in relative reaction rates, whose ratio imposes intrinsic limits to selectivity. Here chemical reactivities of canonical and modified ribonucleosides are revisited as a basis for an understanding of the limits of selectivity achievable with chemical methods.

RNA Footprinting Using Small Chemical Reagents

RNA is a pivotal element of the cell which is most of the time found in complex with protein(s) in a cellular environment. RNA can adopt three-dimensional structures that may form specific binding sites not only for proteins but for all sorts of molecules. Since the early days of molecular biology, strategies to probe RNA structure have been developed. Such probes are small molecules or RNases that most of the time specifically react with single strand nucleotides. The precise reaction or cleavage site can be mapped by reverse transcription. It appears that nucleotides in close contact or in proximity of a ligand are no longer reactive to these probes. Carrying the RNA probing experiment in parallel in presence and absence of a ligand yield differences that are known as the ligand "footprint." Such footprints allow for the identification of the precise site of the ligand interaction, but also reveals RNA structural rearrangement upon ligand binding. Here we provide an experimental and analytical workflow to carry RNA footprinting experiments.

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.

SHAPE Directed Discovery of New Functions in Large RNAs

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

RNA architecture influences plant biology

The majority of the genome is transcribed to RNA in living organisms. RNA transcripts can form astonishing arrays of secondary and tertiary structures via Watson-Crick, Hoogsteen, or wobble base pairing. In vivo, RNA folding is not a simple thermodynamic event of minimizing free energy. Instead, the process is constrained by transcription, RNA-binding proteins, steric factors, and the microenvironment. RNA secondary structure (RSS) plays myriad roles in numerous biological processes, such as RNA processing, stability, transportation, and translation in prokaryotes and eukaryotes. Emerging evidence has also implicated RSS in RNA trafficking, liquid-liquid phase separation, and plant responses to environmental variations such as temperature and salinity. At molecular level, RSS is correlated with splicing, polyadenylation, protein synthesis, and miRNA biogenesis and functions. In this review, we summarize newly reported methods for probing RSS in vivo and functions and mechanisms of RSS in plant physiology.

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.

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.

Translation inhibition from a distance: The small RNA SgrS silences a ribosomal protein S1-dependent enhancer

Many bacterial small RNAs (sRNAs) efficiently inhibit translation of target mRNAs by forming a duplex that sequesters the Shine-Dalgarno (SD) sequence or start codon and prevents formation of the translation initiation complex. There are a growing number of examples of sRNA-mRNA binding interactions distant from the SD region, but how these mediate translational regulation remains unclear. Our previous work in Escherichia coli and Salmonella identified a mechanism of translational repression of manY mRNA by the sRNA SgrS through a binding interaction upstream of the manY SD. Here, we report that SgrS forms a duplex with a uridine-rich translation-enhancing element in the manY 5' untranslated region. Notably, we show that the enhancer is ribosome-dependent and that the small ribosomal subunit protein S1 interacts with the enhancer to promote translation of manY. In collaboration with the chaperone protein Hfq, SgrS interferes with the interaction between the translation enhancer and ribosomal protein S1 to repress translation of manY mRNA. Since bacterial translation is often modulated by enhancer-like elements upstream of the SD, sRNA-mediated enhancer silencing could be a common mode of gene regulation.

Viral RNA structure analysis using DMS-MaPseq

RNA structure is critically important to RNA viruses in every part of the replication cycle. RNA structure is also utilized by DNA viruses in order to regulate gene expression and interact with host factors. Advances in next-generation sequencing have greatly enhanced the utility of chemical probing in order to analyze RNA structure. This review will cover some recent viral RNA structural studies using chemical probing and next-generation sequencing as well as the advantages of dimethyl sulfate (DMS)-mutational profiling and sequencing (MaPseq). DMS-MaPseq is a robust assay that can easily modify RNA in vitro, in cell and in virion. A detailed protocol for whole-genome DMS-MaPseq from cells transfected with HIV-1 and the structure of TAR as determined by DMS-MaPseq is presented. DMS-MaPseq has the ability to answer a variety of integral questions about viral RNA, including how they change in different environments and when interacting with different host factors.

An RNA pseudoknot stimulates HTLV-1 pro-pol programmed -1 ribosomal frameshifting

Programmed -1 ribosomal frameshifts (-1 PRFs) are commonly used by viruses to regulate their enzymatic and structural protein levels. Human T-cell leukemia virus type 1 (HTLV-1) is a carcinogenic retrovirus that uses two independent -1 PRFs to express viral enzymes critical to establishing new HTLV-1 infections. How the cis-acting RNA elements in this viral transcript function to induce frameshifting is unknown. The objective of this work was to conclusively define the 3' boundary of and the RNA elements within the HTLV-1 pro-pol frameshift site. We hypothesized that the frameshift site structure was a pseudoknot and that its 3' boundary would be defined by the pseudoknot's 3' end. To test these hypotheses, the in vitro frameshift efficiencies of three HTLV-1 pro-pol frameshift sites with different 3' boundaries were quantified. The results indicated that nucleotides included in the longest construct were essential to highly efficient frameshift stimulation. Interestingly, only this construct could form the putative frameshift site pseudoknot. Next, the secondary structure of this frameshift site was determined. The dominant structure was an H-type pseudoknot which, together with the slippery sequence, stimulated frameshifting to 19.4(±0.3)%. The pseudoknot's critical role in frameshift stimulation was directly revealed by examining the impact of structural changes on HTLV-1 pro-pol -1 PRF. As predicted, mutations that occluded pseudoknot formation drastically reduced the frameshift efficiency. These results are significant because they demonstrate that a pseudoknot is important to HTLV-1 pro-pol -1 PRF and define the frameshift site's 3' boundary.

Transcription attenuation-derived small RNA rnTrpL regulates tryptophan biosynthesis gene expression in trans

Ribosome-mediated transcription attenuation is a basic posttranscriptional regulation mechanism in bacteria. Liberated attenuator RNAs arising in this process are generally considered nonfunctional. In Sinorhizobium meliloti, the tryptophan (Trp) biosynthesis genes are organized into three operons, trpE(G), ppiD-trpDC-moaC-moeA, and trpFBA-accD-folC, of which only the first one, trpE(G), contains a short ORF (trpL) in the 5′-UTR and is regulated by transcription attenuation. Under conditions of Trp sufficiency, transcription is terminated between trpL and trpE(G), and a small attenuator RNA, rnTrpL, is produced. Here, we show that rnTrpL base-pairs with trpD and destabilizes the polycistronic trpDC mRNA, indicating rnTrpL-mediated downregulation of the trpDC operon in trans. Although all three trp operons are regulated in response to Trp availability, only in the two operons trpE(G) and trpDC the Trp-mediated regulation is controlled by rnTrpL. Together, our data show that the trp attenuator coordinates trpE(G) and trpDC expression posttranscriptionally by two fundamentally different mechanisms: ribosome-mediated transcription attenuation in cis and base-pairing in trans. Also, we present evidence that rnTrpL-mediated regulation of trpDC genes expression in trans is conserved in Agrobacterium and Bradyrhizobium, suggesting that the small attenuator RNAs may have additional conserved functions in the control of bacterial gene expression.

RNA base-pairing complexity in living cells visualized by correlated chemical probing

RNA structure and dynamics are critical to biological function. However, strategies for determining RNA structure in vivo are limited, with established chemical probing and newer duplex detection methods each having deficiencies. Here we convert the common reagent dimethyl sulfate into a useful probe of all 4 RNA nucleotides. Building on this advance, we introduce PAIR-MaP, which uses single-molecule correlated chemical probing to directly detect base-pairing interactions in cells. PAIR-MaP has superior resolution compared to alternative experiments, can resolve multiple sets of pairing interactions for structurally dynamic RNAs, and enables highly accurate structure modeling, including of RNAs containing multiple pseudoknots and extensively bound by proteins. Application of PAIR-MaP to human RNase MRP and 2 bacterial messenger RNA 5' untranslated regions reveals functionally important and complex structures undetected by prior analyses. PAIR-MaP is a powerful, experimentally concise, and broadly applicable strategy for directly visualizing RNA base pairs and dynamics in cells.

RNAdt: An online tutorial and data portal for the RNA structurome era

RNA is not only a passive transporter of genetic information, but also a pivotal player in all domains of life. RNA can regulate gene expression because of its involvement in transcription, mRNA modification and processing, and translation. RNA also possesses other intricate functions such as catalysis, ligand sensing, interaction with biomolecules, response to environment stresses, and information storage. The primary structure of RNA is single stranded, but it always folds into complex secondary and tertiary structures owing to base pairing and effects from the cellular environment. The importance of structure has been increasingly recognized in understanding the myriad functions of RNA. After decades of development, there is a wide range of RNA structure probing techniques. The marriage between structure probing and high-throughput sequencing (HTS) especially enables the measurement of RNA structure on a transcriptomic scale, advancing the advent of the RNA structurome era. Dozens of HTS-associated RNA structure probing methods have been published, so it is urgent to provide a user-friendly and easy-to-use resource for users who are perplexed by selecting the most suitable method for their experiments. Motivated by this demand, we collected currently available HTS-associated RNA structure probing methods and then developed RNAdt (freely accessible at http://www.zhounan.org/rnadt). RNAdt can be used as a web-based tutorial to learn fundamental knowledge of HTS-associated RNA structure probing methods. RNAdt can also be used as a data portal to access HTS data sets from previous RNA structurome studies. At the end of this work, we also provided perspectives on future development of RNA structure probing methods. Our study is expected to facilitate RNA structure probing and ultimately elucidate the connection between RNA structure and biological functions.

Physical and Functional Analysis of Viral RNA Genomes by SHAPE

RNA viruses encode the information required to usurp cellular metabolism and gene regulation and to enable their own replication in two ways: in the linear sequence of their RNA genomes and in higher-order structures that form when the genomic RNA strand folds back on itself. Application of high-resolution SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) structure probing to viral RNA genomes has identified numerous new regulatory elements, defined new principles by which viral RNAs interact with the cellular host and evade host immune responses, and revealed relationships between virus evolution and RNA structure. This review summarizes our current understanding of genome structure-function interrelationships for RNA viruses, as informed by SHAPE structure probing, and outlines opportunities for future studies.

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.

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.

Modeling RNA secondary structure folding ensembles using SHAPE mapping data

RNA secondary structure prediction is widely used for developing hypotheses about the structures of RNA sequences, and structure can provide insight about RNA function. The accuracy of structure prediction is known to be improved using experimental mapping data that provide information about the pairing status of single nucleotides, and these data can now be acquired for whole transcriptomes using high-throughput sequencing. Prior methods for using these experimental data focused on predicting structures for sequences assuming that they populate a single structure. Most RNAs populate multiple structures, however, where the ensemble of strands populates structures with different sets of canonical base pairs. The focus on modeling single structures has been a bottleneck for accurately modeling RNA structure. In this work, we introduce Rsample, an algorithm for using experimental data to predict more than one RNA structure for sequences that populate multiple structures at equilibrium. We demonstrate, using SHAPE mapping data, that we can accurately model RNA sequences that populate multiple structures, including the relative probabilities of those structures. This program is freely available as part of the RNAstructure software package.

Efficiency and precision of microRNA biogenesis modes in plants

Many evolutionarily conserved microRNAs (miRNAs) in plants regulate transcription factors with key functions in development. Hence, mutations in the core components of the miRNA biogenesis machinery cause strong growth defects. An essential aspect of miRNA biogenesis is the precise excision of the small RNA from its precursor. In plants, miRNA precursors are largely variable in size and shape and can be processed by different modes. Here, we optimized an approach to detect processing intermediates during miRNA biogenesis. We characterized a miRNA whose processing is triggered by a terminal branched loop. Plant miRNA processing can be initiated by internal bubbles, small terminal loops or branched loops followed by dsRNA segments of 15–17 bp. Interestingly, precision and efficiency vary with the processing modes. Despite the various potential structural determinants present in a single a miRNA precursor, DCL1 is mostly guided by a predominant structural region in each precursor in wild-type plants. However, our studies in fiery1, hyl1 and se mutants revealed the existence of cleavage signatures consistent with the recognition of alternative processing determinants. The results provide a general view of the mechanisms underlying the specificity of miRNA biogenesis in plants.

SnapShot: RNA Structure Probing Technologies

Chemical probing coupled to high-throughput sequencing offers a flexible approach to uncover many aspects of RNA structure relevant to its cellular function. With a wide variety of chemical probes available that each report on different features of RNA molecules, a broad toolkit exists for investigating in vivo and in vitro RNA structure and interactions with other molecules.

The interplay between molecular flexibility and RNA chemical probing reactivities analyzed at the nucleotide level via an extensive molecular dynamics study

Determination of the tridimensional structure of ribonucleic acid molecules is fundamental for understanding their function in the cell. A common method to investigate RNA structures of large molecules is the use of chemical probes such as SHAPE (2^'-hydroxyl acylation analyzed by primer extension) reagents, DMS (dimethyl sulfate) and CMCT (1-cyclohexyl-3-(2-morpholinoethyl) carbodiimide metho-p-toluene sulfate), the reaction of which is dependent on the local structural properties of each nucleotide. In order to understand the interplay between local flexibility, sugar pucker, canonical pairing and chemical reactivity of the probes, we performed all-atom molecular dynamics simulations on a set of RNA molecules for which both tridimensional structure and chemical probing data are available and we analyzed the correlations between geometrical parameters and the chemical reactivity. Our study confirms that SHAPE reactivity is guided by the local flexibility of the different chemical moieties but suggests that a combination of multiple parameters is needed to better understand the implications of the reactivity at the molecular level. This is also the case for DMS and CMCT for which the reactivity appears to be more complex than commonly accepted.

Advances that facilitate the study of large RNA structure and dynamics by nuclear magnetic resonance spectroscopy

The characterization of functional yet nonprotein coding (nc) RNAs has expanded the role of RNA in the cell from a passive player in the central dogma of molecular biology to an active regulator of gene expression. The misregulation of ncRNA function has been linked with a variety of diseases and disorders ranging from cancers to neurodegeneration. However, a detailed molecular understanding of how ncRNAs function has been limited; due, in part, to the difficulties associated with obtaining high-resolution structures of large RNAs. Tertiary structure determination of RNA as a whole is hampered by various technical challenges, all of which are exacerbated as the size of the RNA increases. Namely, RNAs tend to be highly flexible and dynamic molecules, which are difficult to crystallize. Biomolecular nuclear magnetic resonance (NMR) spectroscopy offers a viable alternative to determining the structure of large RNA molecules that do not readily crystallize, but is itself hindered by some technical limitations. Recently, a series of advancements have allowed the biomolecular NMR field to overcome, at least in part, some of these limitations. These advances include improvements in sample preparation strategies as well as methodological improvements. Together, these innovations pave the way for the study of ever larger RNA molecules that have important biological function. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

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.

In Vivo Genome-Wide RNA Structure Probing with Structure-seq

In vivo genome-wide RNA structure probing provides a global view of RNA structure as it occurs in the cell and can assist in elucidating important functional aspects of RNA structure. Structure-seq2 provides high-quality data on transcriptome-wide RNA structure in vivo but contains numerous steps that require technical precision. In this chapter we present the steps needed to produce high-quality structural data with a focus on controls and troubleshooting. Structure-seq2 can be applied to numerous organisms including plants, humans, and bacteria and is amenable to a wide variety of RNA-modifying chemicals including DMS, glyoxal, and SHAPE reagents. Notably, the data generated by the method highlighted here can be readily analyzed using our StructureFold2 computational pipeline.

Partner-specific prediction of RNA-binding residues in proteins: A critical assessment

RNA-protein interactions play essential roles in regulating gene expression. While some RNA-protein interactions are "specific", that is, the RNA-binding proteins preferentially bind to particular RNA sequence or structural motifs, others are "non-RNA specific." Deciphering the protein-RNA recognition code is essential for comprehending the functional implications of these interactions and for developing new therapies for many diseases. Because of the high cost of experimental determination of protein-RNA interfaces, there is a need for computational methods to identify RNA-binding residues in proteins. While most of the existing computational methods for predicting RNA-binding residues in RNA-binding proteins are oblivious to the characteristics of the partner RNA, there is growing interest in methods for partner-specific prediction of RNA binding sites in proteins. In this work, we assess the performance of two recently published partner-specific protein-RNA interface prediction tools, PS-PRIP, and PRIdictor, along with our own new tools. Specifically, we introduce a novel metric, RNA-specificity metric (RSM), for quantifying the RNA-specificity of the RNA binding residues predicted by such tools. Our results show that the RNA-binding residues predicted by previously published methods are oblivious to the characteristics of the putative RNA binding partner. Moreover, when evaluated using partner-agnostic metrics, RNA partner-specific methods are outperformed by the state-of-the-art partner-agnostic methods. We conjecture that either (a) the protein-RNA complexes in PDB are not representative of the protein-RNA interactions in nature, or (b) the current methods for partner-specific prediction of RNA-binding residues in proteins fail to account for the differences in RNA partner-specific versus partner-agnostic protein-RNA interactions, or both.

2'-Fluoro-Pyrimidine-Modified RNA Aptamers Specific for Lipopolysaccharide Binding Protein (LBP)

Lipopolysaccaride binding protein (LBP), a glycosylated acute phase protein, plays an important role in the pathophysiology of sepsis. LBP binds with high affinity to the lipid part of bacterial lipopolysaccaride (LPS). Inhibition of the LPS-LBP interaction or blockage of LBP-mediated transfer of LPS monomers to CD14 may be therapeutical strategies to prevent septic shock. LBP is also of interest as a biomarker to identify septic patients at high risk for death, as LBP levels are elevated during early stages of severe sepsis. As a first step toward such potential applications, we isolated aptamers specific for murine LBP (mLBP) by in vitro selection from a library containing a 60-nucleotide randomized region. Modified RNA pools were transcribed in the presence of 2′-fluoro-modified pyrimidine nucleotides to stabilize transcripts against nuclease degradation. As verified for one aptamer experimentally, the selected aptamers adopt a “three-helix junction” architecture, presenting single-stranded 7-nt (5′-YGCTTCY) or 6-nt (5′-RTTTCY) consensus sequences in their core. The best binder (aptamer A011; K_d of 270 nM for binding to mLBP), characterized in more detail by structure probing and boundary analysis, was demonstrated to bind with high specificity to murine LBP.

The RNA Base-Pairing Problem and Base-Pairing Solutions

RNA molecules are folded into structures and complexes to perform a wide variety of functions. Determination of RNA structures and their interactions is a fundamental problem in RNA biology. Most RNA molecules in living cells are large and dynamic, posing unique challenges to structure analysis. Here we review progress in RNA structure analysis, focusing on methods that use the "cross-link, proximally ligate, and sequence" principle for high-throughput detection of base-pairing interactions in living cells. Beginning with a comparison of commonly used methods in structure determination and a brief historical account of psoralen cross-linking studies, we highlight the important features of cross-linking methods and new biological insights into RNA structures and interactions from recent studies. Further improvement of these cross-linking methods and application to previously intractable problems will shed new light on the mechanisms of the "modern RNA world."