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.

New RoxS sRNA Targets Identified in Bacillus subtilis by Pulsed SILAC.

Non-coding RNAs (sRNA) play a key role in controlling gene expression in bacteria, typically by base-pairing with ribosome binding sites to block translation. The modification of ribosome traffic along the mRNA generally affects its stability. However, a few cases have been described in bacteria where sRNAs can affect translation without a major impact on mRNA stability. To identify new sRNA targets in Bacillus subtilis potentially belonging to this class of mRNAs, we used pulsed-SILAC (stable isotope labeling by amino acids in cell culture) to label newly synthesized proteins after short expression of the RoxS sRNA, the best characterized sRNA in this bacterium. RoxS sRNA was previously shown to interfere with the expression of genes involved in central metabolism, permitting control of the NAD+/NADH ratio in B. subtilis. In this study, we confirmed most of the known targets of RoxS, showing the efficiency of the method. We further expanded the number of mRNA targets encoding enzymes of the TCA cycle and identified new targets. One of these is YcsA, a tartrate dehydrogenase that uses NAD+ as co-factor, in excellent agreement with the proposed role of RoxS in management of NAD+/NADH ratio in Firmicutes. IMPORTANCE Non-coding RNAs (sRNA) play an important role in bacterial adaptation and virulence. The identification of the most complete set of targets for these regulatory RNAs is key to fully identifying the perimeter of its function(s). Most sRNAs modify both the translation (directly) and mRNA stability (indirectly) of their targets. However, sRNAs can also influence the translation efficiency of the target primarily, with little or no impact on mRNA stability. The characterization of these targets is challenging. We describe here the application of the pulsed SILAC method to identify such targets and obtain the most complete list of targets for a defined sRNA.

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.

Identification of an RNA sponge that controls the RoxS riboregulator of central metabolism in Bacillus subtilis.

sRNAs are a taxonomically-restricted but transcriptomically-abundant class of post-transcriptional regulators. While of major importance for adaption to the environment, we currently lack global-scale methodology enabling target identification, especially in species without known RNA hub proteins (e.g. Hfq). Using psoralen RNA cross-linking and Illumina-sequencing we identify RNA-RNA interacting pairs in vivo in Bacillus subtilis, resolving previously well-described interactants. Although sRNA-sRNA pairings are rare (compared with sRNA-mRNA), we identify a robust example involving the conserved sRNA RoxS and an unstudied sRNA RosA (Regulator of sRNA A). We show RosA to be the first confirmed RNA sponge described in a Gram-positive bacterium. RosA interacts with at least two sRNAs, RoxS and FsrA. The RosA/RoxS interaction not only affects the levels of RoxS but also its processing and regulatory activity. We also found that the transcription of RosA is repressed by CcpA, the key regulator of carbon-metabolism in B. subtilis. Since RoxS is already known to be transcriptionally controlled by malate via the transcriptional repressor Rex, its post-transcriptional regulation by CcpA via RosA places RoxS in a key position to control central metabolism in response to varying carbon sources.

Assay of Bacillus subtilis Ribonuclease Activity In Vitro.

Ribonucleases can cleave RNAs internally in endoribonucleolytic mode or remove one nucleotide at a time from either the 5' or 3' end through exoribonuclease action. To show direct implication of an RNase in a specific pathway of RNA maturation or decay requires the setting up of in vitro assays with purified enzymes and substrates. This chapter complements Chapter 24 on assays of ribonuclease action in vivo by providing detailed protocols for the assay of B. subtilis RNases with prepared substrates in vitro.

Analysis of Bacillus subtilis Ribonuclease Activity In Vivo.

Ribonucleases remodel RNAs to render them functional or to send them on their way toward degradation. In our laboratory, we study these pathways in detail using a plethora of different techniques. These can range from the isolation of RNAs in various RNase mutants to determine their implication in maturation or decay pathways by Northern blot, to proving their direct roles in RNA cleavage reactions using purified enzymes and transcribed substrates in vitro. In this chapter, we provide in-depth protocols for the techniques we use daily in the laboratory to assay RNase activity in vivo, with detailed notes on how to get these methods to work optimally. This chapter complements Chapter 25 on assays of ribonuclease action in vitro.

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.

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.

Non-canonical translation initiation of the spliced mRNA encoding the human T-cell leukemia virus type 1 basic leucine zipper protein.

Human T-cell leukemia virus type 1 (HTLV-1) is the etiological agent of adult T-cell leukemia (ATL). The HTLV-1 basic leucine zipper protein (HBZ) is expressed in all cases of ATL and is directly associated with virus pathogenicity. The two isoforms of the HBZ protein are synthesized from antisense messenger RNAs (mRNAs) that are either spliced (sHBZ) or unspliced (usHBZ) versions of the HBZ transcript. The sHBZ and usHBZ mRNAs have entirely different 5'untranslated regions (5'UTR) and are differentially expressed in cells, with the sHBZ protein being more abundant. Here, we show that differential expression of the HBZ isoforms is regulated at the translational level. Translation initiation of the usHBZ mRNA relies on a cap-dependent mechanism, while the sHBZ mRNA uses internal initiation. Based on the structural data for the sHBZ 5'UTR generated by SHAPE in combination with 5' and 3' deletion mutants, the minimal region harboring IRES activity was mapped to the 5'end of the sHBZ mRNA. In addition, the sHBZ IRES recruited the 40S ribosomal subunit upstream of the initiation codon, and IRES activity was found to be dependent on the ribosomal protein eS25 and eIF5A.

Distinct roles for the IIId2 sub-domain in pestivirus and picornavirus internal ribosome entry sites.

Viral internal ribosomes entry site (IRES) elements coordinate the recruitment of the host translation machinery to direct the initiation of viral protein synthesis. Within hepatitis C virus (HCV)-like IRES elements, the sub-domain IIId(1) is crucial for recruiting the 40S ribosomal subunit. However, some HCV-like IRES elements possess an additional sub-domain, termed IIId2, whose function remains unclear. Herein, we show that IIId2 sub-domains from divergent viruses have different functions. The IIId2 sub-domain present in Seneca valley virus (SVV), a picornavirus, is dispensable for IRES activity, while the IIId2 sub-domains of two pestiviruses, classical swine fever virus (CSFV) and border disease virus (BDV), are required for 80S ribosomes assembly and IRES activity. Unlike in SVV, the deletion of IIId2 from the CSFV and BDV IRES elements impairs initiation of translation by inhibiting the assembly of 80S ribosomes. Consequently, this negatively affects the replication of CSFV and BDV. Finally, we show that the SVV IIId2 sub-domain is required for efficient viral RNA synthesis and growth of SVV, but not for IRES function. This study sheds light on the molecular evolution of viruses by clearly demonstrating that conserved RNA structures, within distantly related RNA viruses, have acquired different roles in the virus life cycles.