Visualization and characterization of RNA–protein interactions in living cells

RNA–protein interactions are the structural and functional basis of significant numbers of RNA molecules. RNA–protein interaction assays though, still mainly depend on biochemical tests in vitro. Here, we establish a convenient and reliable RNA fluorescent three-hybrid (rF3H) method to detect/interrogate the interactions between RNAs and proteins in cells. A GFP tagged highly specific RNA trap is constructed to anchor the RNA of interest to an artificial or natural subcellular structure, and RNA–protein interactions can be detected and visualized by the enrichment of RNA binding proteins (RBPs) at these structures. Different RNA trapping systems are developed and detection of RNA–protein complexes at multiple subcellular structures are assayed. With this new toolset, interactions between proteins and mRNA or noncoding RNAs are characterized, including the interaction between a long noncoding RNA and an epigenetic modulator. Our approach provides a flexible and reliable method for the characterization of RNA–protein interactions in living cells.     

The role of disorder in RNA binding affinity and specificity

Most RNA-binding modules are small and bind few nucleotides. RNA-binding proteins typically attain the physiological specificity and affinity for their RNA targets by combining several RNA-binding modules. Here, we review how disordered linkers connecting RNA-binding modules govern the specificity and affinity of RNA–protein interactions by regulating the effective concentration of these modules and their relative orientation. RNA-binding proteins also often contain extended intrinsically disordered regions that mediate protein–protein and RNA–protein interactions with multiple partners. We discuss how these regions can connect proteins and RNA resulting in heterogeneous higher-order assemblies such as membrane-less compartments and amyloid-like structures that have the characteristics of multi-modular entities. The assembled state generates additional RNA-binding specificity and affinity properties that contribute to further the function of RNA-binding proteins within the cellular environment.     

Evidence that avian reovirus σNS is an RNA chaperone: implications for genome segment assortment

Reoviruses are important human, animal and plant pathogens having 10–12 segments of double-stranded genomic RNA. The mechanisms controlling the assortment and packaging of genomic segments in these viruses, remain poorly understood. RNA–protein and RNA–RNA interactions between viral genomic segment precursors have been implicated in the process. While non-structural viral RNA-binding proteins, such as avian reovirus σNS, are essential for virus replication, the mechanism by which they assist packaging is unclear. Here we demonstrate that σNS assembles into stable elongated hexamers in vitro, which bind single-stranded nucleic acids with high affinity, but little sequence specificity. Using ensemble and single molecule fluorescence spectroscopy, we show that σNS also binds to a partially double-stranded RNA, resulting in gradual helix unwinding. The hexamer can bind multiple RNA molecules and exhibits strand-annealing activity, thus mediating conversion of metastable, intramolecular stem-loops into more stable heteroduplexes. We demonstrate that the ARV σNS acts as an RNA chaperone facilitating specific RNA–RNA interactions between genomic precursors during segment assortment and packaging.     

Genome-Wide Double-Stranded RNA Sequencing Reveals the Functional Significance of Base-Paired RNAs in Arabidopsis

The functional structure of all biologically active molecules is dependent on intra- and inter-molecular interactions. This is especially evident for RNA molecules whose functionality, maturation, and regulation require formation of correct secondary structure through encoded base-pairing interactions. Unfortunately, intra- and inter-molecular base-pairing information is lacking for most RNAs. Here, we marry classical nuclease-based structure mapping techniques with high-throughput sequencing technology to interrogate all base-paired RNA in Arabidopsis thaliana and identify ∼200 new small (sm)RNA–producing substrates of RNA–DEPENDENT RNA POLYMERASE6. Our comprehensive analysis of paired RNAs reveals conserved functionality within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs, as well as a novel population of functional RNAs, many of which are the precursors of smRNAs. Finally, we identify intra-molecular base-pairing interactions to produce a genome-wide collection of RNA secondary structure models. Although our methodology reveals the pairing status of RNA molecules in the absence of cellular proteins, previous studies have demonstrated that structural information obtained for RNAs in solution accurately reflects their structure in ribonucleoprotein complexes. Furthermore, our identification of RNA–DEPENDENT RNA POLYMERASE6 substrates and conserved functional RNA domains within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs using this approach strongly suggests that RNA molecules are correctly folded into their secondary structure in solution. Overall, our findings highlight the importance of base-paired RNAs in eukaryotes and present an approach that should be widely applicable for the analysis of this key structural feature of RNA.     

Remodeling of ribonucleoprotein complexes with DExH/D RNA helicases

The DExH/D protein family is the largest group of enzymes in eukaryotic RNA metabolism. DExH/D proteins are mainly known for their ability to unwind RNA duplexes in an ATP-dependent fashion. However, it has become clear in recent years that these DExH/D RNA helicases are also involved in the ATP-dependent remodeling of RNA–protein complexes. Here we review recent studies that highlight physiological roles of DExH/D proteins in the displacement of proteins from RNA. We further discuss work with simple RNA–protein complexes in vitro, which illuminates mechanisms by which DExH/D proteins remove proteins from RNA. Although we are only beginning to understand how DExH/D proteins remodel RNA–protein complexes, these studies have shown that an ‘RNA helicase’ does not per se require cofactors to displace proteins from RNA, that protein displacement does not necessarily involve RNA duplex unwinding, and that not all DExH/D proteins are able to disassemble the same range of ribonucleoproteins.     

RNAnue: efficient data analysis for RNA–RNA interactomics

RNA–RNA inter- and intramolecular interactions are fundamental for numerous biological processes. While there are reasonable approaches to map RNA secondary structures genome-wide, understanding how different RNAs interact to carry out their regulatory functions requires mapping of intermolecular base pairs. Recently, different strategies to detect RNA–RNA duplexes in living cells, so called direct duplex detection (DDD) methods, have been developed. Common to all is the Psoralen-mediated in vivo RNA crosslinking followed by RNA Proximity Ligation to join the two interacting RNA strands. Sequencing of the RNA via classical RNA-seq and subsequent specialised bioinformatic analyses the result in the prediction of inter- and intramolecular RNA–RNA interactions. Existing approaches adapt standard RNA-seq analysis pipelines, but often neglect inherent features of RNA–RNA interactions that are useful for filtering and statistical assessment. Here we present RNAnue, a general pipeline for the inference of RNA–RNA interactions from DDD experiments that takes into account hybridisation potential and statistical significance to improve prediction accuracy. We applied RNAnue to data from different DDD studies and compared our results to those of the original methods. This showed that RNAnue performs better in terms of quantity and quality of predictions.     

End-to-end crosstalk within the hepatitis C virus genome mediates the conformational switch of the 3′X-tail region

The hepatitis C virus (HCV) RNA genome contains multiple structurally conserved domains that make long-distance RNA–RNA contacts important in the establishment of viral infection. Microarray antisense oligonucelotide assays, improved dimethyl sulfate probing methods and 2′ acylation chemistry (selective 2’-hydroxyl acylation and primer extension, SHAPE) showed the folding of the genomic RNA 3′ end to be regulated by the internal ribosome entry site (IRES) element via direct RNA–RNA interactions. The essential cis-acting replicating element (CRE) and the 3′X-tail region adopted different 3D conformations in the presence and absence of the genomic RNA 5′ terminus. Further, the structural transition in the 3′X-tail from the replication-competent conformer (consisting of three stem-loops) to the dimerizable form (with two stem-loops), was found to depend on the presence of both the IRES and the CRE elements. Complex interplay between the IRES, the CRE and the 3′X-tail region would therefore appear to occur. The preservation of this RNA–RNA interacting network, and the maintenance of the proper balance between different contacts, may play a crucial role in the switch between different steps of the HCV cycle.     

Insights into structure, dynamics and hydration of locked nucleic acid (LNA) strand-based duplexes from molecular dynamics simulations

Locked nucleic acid (LNA) is a chemically modified nucleic acid with its sugar ring locked in an RNA-like (C3′-endo) conformation. LNAs show extraordinary thermal stabilities when hybridized with DNA, RNA or LNA itself. We performed molecular dynamics simulations on five isosequential duplexes (LNA–DNA, LNA–LNA, LNA–RNA, RNA–DNA and RNA–RNA) in order to characterize their structure, dynamics and hydration. Structurally, the LNA–DNA and LNA–RNA duplexes are found to be similar to regular RNA–DNA and RNA–RNA duplexes, whereas the LNA–LNA duplex is found to have its helix partly unwound and does not resemble RNA–RNA duplex in a number of properties. Duplexes with an LNA strand have on average longer interstrand phosphate distances compared to RNA–DNA and RNA–RNA duplexes. Furthermore, intrastrand phosphate distances in LNA strands are found to be shorter than in DNA and slightly shorter than in RNA. In case of induced sugar puckering, LNA is found to tune the sugar puckers in partner DNA strand toward C3′-endo conformations more efficiently than RNA. The LNA–LNA duplex has lesser backbone flexibility compared to the RNA–RNA duplex. Finally, LNA is less hydrated compared to DNA or RNA but is found to have a well-organized water structure.     

Anatomy of noncovalent interactions between the nucleobases or ribose and π-containing amino acids in RNA–protein complexes

A set of >300 nonredundant high-resolution RNA–protein complexes were rigorously searched for π-contacts between an amino acid side chain (W, H, F, Y, R, E and D) and an RNA nucleobase (denoted π–π interaction) or ribose moiety (denoted sugar–π). The resulting dataset of >1500 RNA–protein π-contacts were visually inspected and classified based on the interaction type, and amino acids and RNA components involved. More than 80% of structures searched contained at least one RNA–protein π-interaction, with π–π contacts making up 59% of the identified interactions. RNA–protein π–π and sugar–π contacts exhibit a range in the RNA and protein components involved, relative monomer orientations and quantum mechanically predicted binding energies. Interestingly, π–π and sugar–π interactions occur more frequently with RNA (4.8 contacts/structure) than DNA (2.6). Moreover, the maximum stability is greater for RNA–protein contacts than DNA–protein interactions. In addition to highlighting distinct differences between RNA and DNA–protein binding, this work has generated the largest dataset of RNA–protein π-interactions to date, thereby underscoring that RNA–protein π-contacts are ubiquitous in nature, and key to the stability and function of RNA–protein complexes.     

Ribosomal protein L1 recognizes the same specific structural motif in its target sites on the autoregulatory mRNA and 23S rRNA

The RNA-binding ability of ribosomal protein L1 is of profound interest since the protein has a dual function as a ribosomal protein binding rRNA and as a translational repressor binding its mRNA. Here, we report the crystal structure of ribosomal protein L1 in complex with a specific fragment of its mRNA and compare it with the structure of L1 in complex with a specific fragment of 23S rRNA determined earlier. In both complexes, a strongly conserved RNA structural motif is involved in L1 binding through a conserved network of RNA–protein H-bonds inaccessible to the solvent. These interactions should be responsible for specific recognition between the protein and RNA. A large number of additional non-conserved RNA–protein H-bonds stabilizes both complexes. The added contribution of these non-conserved H-bonds makes the ribosomal complex much more stable than the regulatory one.     

Recognition of viral RNA stem–loops by the tandem double-stranded RNA binding domains of PKR

In humans, the double-stranded RNA (dsRNA)-activated protein kinase (PKR) is expressed in late stages of the innate immune response to viral infection by the interferon pathway. PKR consists of tandem dsRNA binding motifs (dsRBMs) connected via a flexible linker to a Ser/Thr kinase domain. Upon interaction with viral dsRNA, PKR is converted into a catalytically active enzyme capable of phosphorylating a number of target proteins that often results in host cell translational repression. A number of high-resolution structural studies involving individual dsRBMs from proteins other than PKR have highlighted the key features required for interaction with perfectly duplexed RNA substrates. However, viral dsRNA molecules are highly structured and often contain deviations from perfect A-form RNA helices. By use of small-angle X-ray scattering (SAXS), we present solution conformations of the tandem dsRBMs of PKR in complex with two imperfectly base-paired viral dsRNA stem–loops; HIV-1 TAR and adenovirus VA_I-AS. Both individual components and complexes were purified by size exclusion chromatography and characterized by dynamic light scattering at multiple concentrations to ensure monodispersity. SAXS ab initio solution conformations of the individual components and RNA–protein complexes were determined and highlight the potential of PKR to interact with both stem and loop regions of the RNA. Excellent agreement between experimental and model-based hydrodynamic parameter determination heightens our confidence in the obtained models. Taken together, these data support and provide a framework for the existing biochemical data regarding the tolerance of imperfectly base-paired viral dsRNA by PKR.     

RNA aptamers that specifically bind to a 16S ribosomal RNA decoding region construct

RNA–RNA recognition is a critical process in controlling many key biological events, such as translation and ribozyme functions. The recognition process governing RNA–RNA interactions can involve complementary Watson–Crick (WC) base pair binding, or can involve binding through tertiary structural interaction. Hence, it is of interest to determine which of the RNA–RNA binding events might emerge through an in vitro selection process. The A-site of the 16S rRNA decoding region was chosen as the target, both because it possesses several different RNA structural motifs, and because it is the rRNA site where codon/anticodon recognition occurs requiring recognition of both mRNA and tRNA. It is shown here that a single family of RNA molecules can be readily selected from two different sizes of RNA library. The tightest binding aptamer to the A-site 16S rRNA construct, 109.2-3, has its consensus sequences confined to a stem–loop region, which contains three nucleotides complementary to three of the four nucleotides in the stem–loop region of the A-site 16S rRNA. Point mutations on each of the three nucleotides on the stem–loop of the aptamer abolish its binding capacity. These studies suggest that the RNA aptamer 109.2-3 interacts with the simple 27 nt A-site decoding region of 16S rRNA through their respective stem–loops. The most probable mode of interaction is through complementary WC base pairing, commonly referred to as a loop–loop ‘kissing’ motif. High affinity binding to the other structural motifs in the decoding region were not observed.     

A long-range RNA–RNA interaction between the 5′ and 3′ ends of the HCV genome

The RNA genome of the hepatitis C virus (HCV) contains multiple conserved structural cis domains that direct protein synthesis, replication, and infectivity. The untranslatable regions (UTRs) play essential roles in the HCV cycle. Uncapped viral RNAs are translated via an internal ribosome entry site (IRES) located at the 5′ UTR, which acts as a scaffold for recruiting multiple protein factors. Replication of the viral genome is initiated at the 3′ UTR. Bioinformatics methods have identified other structural RNA elements thought to be involved in the HCV cycle. The 5BSL3.2 motif, which is embedded in a cruciform structure at the 3′ end of the NS5B coding sequence, contributes to the three-dimensional folding of the entire 3′ end of the genome. It is essential in the initiation of replication. This paper reports the identification of a novel, strand-specific, long-range RNA–RNA interaction between the 5′ and 3′ ends of the genome, which involves 5BSL3.2 and IRES motifs. Mutants harboring substitutions in the apical loop of domain IIId or in the internal loop of 5BSL3.2 disrupt the complex, indicating these regions are essential in initiating the kissing interaction. No complex was formed when the UTRs of the related foot and mouth disease virus were used in binding assays, suggesting this interaction is specific for HCV sequences. The present data firmly suggest the existence of a higher-order structure that may mediate a protein-independent circularization of the HCV genome. The 5′–3′ end bridge may have a role in viral translation modulation and in the switch from protein synthesis to RNA replication.     

Directed evolution of orthogonal RNA–RBP pairs through library-vs-library in vitro selection

RNA-binding proteins (RBPs) and their RNA ligands play many critical roles in gene regulation and RNA processing in cells. They are also useful for various applications in cell biology and synthetic biology. However, re-engineering novel and orthogonal RNA–RBP pairs from natural components remains challenging while such synthetic RNA–RBP pairs could significantly expand the RNA–RBP toolbox for various applications. Here, we report a novel library-vs-library in vitro selection strategy based on Phage Display coupled with Systematic Evolution of Ligands by EXponential enrichment (PD-SELEX). Starting with pools of 1.1 × 10¹² unique RNA sequences and 4.0 × 10⁸ unique phage-displayed L7Ae-scaffold (LS) proteins, we selected RNA–RBP complexes through a two-step affinity purification process. After six rounds of library-vs-library selection, the selected RNAs and LS proteins were analyzed by next-generation sequencing (NGS). Further deconvolution of the enriched RNA and LS protein sequences revealed two synthetic and orthogonal RNA–RBP pairs that exhibit picomolar affinity and >4000-fold selectivity.