Predicting novel RNA-RNA interactions

The purpose of this article is to give a brief, yet concise overview of the current computational methods for predicting novel RNA-RNA interactions, that is interactions whose characteristic features we do not yet know. We start by briefly reviewing experimentally confirmed examples of RNA-RNA interactions before introducing computational methods for predicting RNA-RNA interactions. We will focus primarily on the interactions between different RNA molecules, that is trans RNA-RNA interactions, and will only discuss methods for predicting RNA structure, that is cis-only RNA-RNA interactions, where this helps to gain a better understanding. We conclude by discussing the merits of the different approaches and provide an outlook on probably and desirable future developments in the field.     

Long-range RNA-RNA interactions circularize the dengue virus genome

Secondary and tertiary RNA structures present in viral RNA genomes play essential regulatory roles during translation, RNA replication, and assembly of new viral particles. In the case of flaviviruses, RNA-RNA interactions between the 5' and 3' ends of the genome have been proposed to be required for RNA replication. We found that two RNA elements present at the ends of the dengue virus genome interact in vitro with high affinity. Visualization of individual molecules by atomic force microscopy revealed that physical interaction between these RNA elements results in cyclization of the viral RNA. Using RNA binding assays, we found that the putative cyclization sequences, known as 5' and 3' CS, present in all mosquito-borne flaviviruses, were necessary but not sufficient for RNA-RNA interaction. Additional sequences present at the 5' and 3' untranslated regions of the viral RNA were also required for RNA-RNA complex formation. We named these sequences 5' and 3' UAR (upstream AUG region). In order to investigate the functional role of 5'-3' UAR complementarity, these sequences were mutated either separately, to destroy base pairing, or simultaneously, to restore complementarity in the context of full-length dengue virus RNA. Nonviable viruses were recovered after transfection of dengue virus RNA carrying mutations either at the 5' or 3' UAR, while the RNA containing the compensatory mutations was able to replicate. Since sequence complementarity between the ends of the genome is required for dengue virus viability, we propose that cyclization of the RNA is a required conformation for viral replication.     

Short-range RNA-RNA crosslinking methods to determine rRNA structure and interactions

We describe details of procedures to analyze RNA-RNA crosslinks made by far-UV irradiation (< 300 nm) or made by irradiation with near-UV light (320-365 nm) on RNA containing photosensitive nucleotides, in the present case containing 4-thiouridine. Zero-length crosslinks of these types must occur because of the close proximity of the participants through either specific interactions or transient contacts in the folded RNA structure, so they are valuable monitors of the conformation of the RNA. Procedures to produce crosslinks in the 16S ribosomal RNA and between the 16S rRNA and mRNA or tRNA are described. Gel electrophoresis conditions are described that separate the products according to their structure to allow the determination of the number and frequency of the crosslinking products. Gel electrophoresis together with an ultracentrifugation procedure for the efficient recovery of RNA from the polyacrylamide gels allows the purification of molecules containing different crosslinks. These separation techniques allow the analysis of the sites of crosslinking by primer extension and RNA sequencing techniques. The procedures are applicable to other types of RNA molecules with some differences to control levels of crosslinking and separation conditions.     

Determinants of apical loop-internal loop RNA-RNA interactions involving the HCV IRES

Domain II of the hepatitis C virus internal ribosome entry site is a major RNA structure involved in the viral mRNA translation. It comprises four different structural domains. We performed in vitro selection against the apical loop of the domain II and we identified RNA aptamers folding as an imperfect hairpin with an internal loop of interacting with the apical loop of the domain II. This RNA-RNA interaction creates apical loop-internal loop complex. The aptamer binds the target with an apparent K(d) of 35nM. In this study, the main structural elements of the target and the aptamer involved in the formation of the complex are characterized by mutation, deletion, and RNase probing analysis. We demonstrate that a complementary loop flanked by G,C rich upper and lower stems are crucial for such RNA-RNA interactions.     

In-gel probing of individual RNA conformers within a mixed population reveals a dimerization structural switch in the HIV-1 leader

Definitive secondary structural mapping of RNAs in vitro can be complicated by the presence of more than one structural conformer or multimerization of some of the molecules. Until now, probing a single structure of conformationally flexible RNA molecules has typically relied on introducing stabilizing mutations or adjusting buffer conditions or RNA concentration. Here, we present an in-gel SHAPE (selective 2′OH acylation analysed by primer extension) approach, where a mixed structural population of RNA molecules is separated by non-denaturing gel electrophoresis and the conformers are individually probed within the gel matrix. Validation of the technique using a well-characterized RNA stem-loop structure, the HIV-1 trans-activation response element, showed that authentic structure was maintained and that the method was accurate and highly reproducible. To further demonstrate the utility of in-gel SHAPE, we separated and examined monomeric and dimeric species of the HIV-1 packaging signal RNA. Extensive differences in acylation sensitivity were seen between monomer and dimer. The results support a recently proposed structural switch model of RNA genomic dimerization and packaging, and demonstrate the discriminatory power of in-gel SHAPE.     

Hierarchical folding of multiple sequence alignments for the prediction of structures and RNA-RNA interactions

Background  

Many regulatory non-coding RNAs (ncRNAs) function through complementary binding with mRNAs or other ncRNAs, e.g., microRNAs, snoRNAs and bacterial sRNAs. Predicting these RNA interactions is essential for functional studies of putative ncRNAs or for the design of artificial RNAs. Many ncRNAs show clear signs of undergoing compensating base changes over evolutionary time. Here, we postulate that a non-negligible part of the existing RNA-RNA interactions contain preserved but covarying patterns of interactions.     

Thermodynamics of RNA-RNA binding

BACKGROUND: Reliable prediction of RNA-RNA binding energies is crucial, e.g. for the understanding on RNAi, microRNA-mRNA binding and antisense interactions. The thermodynamics of such RNA-RNA interactions can be understood as the sum of two energy contributions: (1) the energy necessary to 'open' the binding site and (2) the energy gained from hybridization. METHODS: We present an extension of the standard partition function approach to RNA secondary structures that computes the probabilities Pu[i, j] that a sequence interval [i, j] is unpaired. RESULTS: Comparison with experimental data shows that Pu[i, j] can be applied as a significant determinant of local target site accessibility for RNA interference (RNAi). Furthermore, these quantities can be used to rigorously determine binding free energies of short oligomers to large mRNA targets. The resource consumption is comparable with a single partition function computation for the large target molecule. We can show that RNAi efficiency correlates well with the binding energies of siRNAs to their respective mRNA target. AVAILABILITY: RNAup will be distributed as part of the Vienna RNA Package, www.tbi.univie.ac.at/~ivo/RNA/     

Single molecule mechanical probing of the SNARE protein interactions

Exocytotic release of neurotransmitters is mediated by the ternary soluble N-ethyl maleimide-sensitive fusion protein attachment protein receptors (SNAREs) complex, comprised of syntaxin (Sx), synaptosome-associated protein of 25 kDa (SNAP25), and synaptobrevin 2 (Sb2). Since exocytosis involves the nonequilibrium process of association and dissociation of bonds between molecules of the SNARE complex, dynamic measurements at the single molecule level are necessary for a detailed understanding of these interactions. To address this issue, we used the atomic force microscope in force spectroscopy mode to show from single molecule investigations of the SNARE complex, that Sx1A and Sb2 are zippered throughout their entire SNARE domains without the involvement of SNAP25. When SNAP25B is present in the complex, it creates a local interaction at the 0 (ionic) layer by cuffing Sx1A and Sb2. Force loading rate studies indicate that the ternary complex interaction is more stable than the Sx1A-Sb2 interaction.     

Long range RNA-RNA interactions in the 30 S ribosomal subunit of E. coli.

We have attempted to identify long-range interactions in the tertiary structure of RNA in the E. coli 30 S ribosome. Native subunits were cleaved with ribonuclease and separated into nucleoprotein fragments which were deproteinized and fractionated into multi-oligonucleotide complexes under conditions intended to preserve RNA-RNA interactions. The final products were denatured by urea and heat and their constituent oligonucleotides resolved and sequenced. Many complexes contained complementary sequences known to be bound together in the RNA secondary structure, attesting to the validity of the technique. Other co-migrating oligonucleotides, not joined in the secondary structure, contained mutually complementary sequences in locations that allow base-pairing interaction without disrupting pre-existing secondary structure. In seven instances the complementary relationship was found to have been preserved during phylogenetic diversification.     

Multiplexed imaging for probing RAS-RAF interactions in living cells

GTP-bound RAS interacts with its protein effectors in response to extracellular stimuli, leading to chemical inputs for downstream pathways. Significant progress has been made in measuring these reversible protein-protein interactions (PPIs) in various cell-free environments. Yet, acquiring high sensitivity in heterogeneous solutions remains challenging. Here, using an intermolecular fluorescence resonance energy transfer (FRET) biosensing approach, we develop a method to visualize and localize HRAS-CRAF interactions in living cells. We demonstrate that the EGFR activation and the HRAS-CRAF complex formation can be concurrently probed in a single cell. This biosensing strategy discriminates EGF-stimulated HRAS-CRAF interactions at the cell and organelle membranes. In addition, we provide quantitative FRET measurements for assessing these transient PPIs in a cell-free environment. Finally, we prove the utility of this approach by showing that an EGFR-binding compound is a potent inhibitor of HRAS-CRAF interactions. The outcomes of this work form a fundamental basis for further explorations of the spatiotemporal dynamics of various signaling networks.     

Long-distance RNA-RNA interactions and conserved sequence elements affect potato virus X plus-strand RNA accumulation.

Conserved octanucleotide sequences located upstream of two major potato virus X (PVX) subgenomic RNAs (sgRNAs), as well as elements in the 5' end of the genome, affect accumulation of sgRNA. To determine if complementarity between these sequences is important for PVX RNA accumulation, we analyzed the effects of mutations within these elements and compensatory mutations in a tobacco protoplast system and in plants. Mutations in the 5' nontranslated region (NTR mutants) that reduced complementarity resulted in lower genomic RNA (gRNA) and sgRNA levels, whereas mutations to the octanucleotide elements affected only the corresponding sgRNA levels. However, for both the NTR and octanucleotide mutants, the extent of reductions in RNA levels did not directly correlate with the degree of complementarity, suggesting that the sequences of these elements are also important. Mutants containing changes in the NTR and compensatory changes in one of the octanucleotide elements restored levels of gRNA and the other sgRNA species with an unaltered octanucleotide element to those of wild-type. Although compensatory changes significantly increased levels of the sgRNA species with the modified octanucleotide element, levels were not restored to those of wild-type. Our data indicate that long distance RNA-RNA interactions and the sequences of the interacting elements are required for PVX plus-strand RNA accumulation.