Abstract:Abstract RNA molecules can adopt stable secondary and tertiary structures, which is essential in mediating physical interactions with other partners such as RNA binding proteins (RBPs) and in carrying out their cellular functions. In vivo and in vitro experiments such as RNAcompete and eCLIP have revealed in vitro binding preferences of RBPs to RNA oligomers and in vivo binding sites in cells. Analysis of these binding data showed that the structure properties of the RNAs in these binding sites are important determinants of the binding events; however, it has been a challenge to incorporate the structure information into an interpretable model. Here we describe a new approach, RNANetMotif, which takes predicted secondary structure of thousands of RNA sequences bound by an RBP as input and uses a graph theory approach to recognize enriched subgraphs. These enriched subgraphs are in essence shared sequence-structure elements that are important in RBP-RNA binding. To validate our approach, we performed RNA structure modeling via discrete molecular dynamics folding simulations for selected 4 RBPs, and RNA-protein docking for LIN28. The simulation results, e.g., solvent accessibility and energetics, further support the biological relevance of the discovered network subgraphs. Author Summary RNA binding proteins (RBPs) regulate every aspect of RNA biology, including splicing, translation, transportation, and degradation. High-throughput technologies such as eCLIP have identified thousands of binding sites for a given RBP throughout the genome. It has been shown by earlier studies that, in addition to nucleotide sequences, the structure and conformation of RNAs also play important role in RBP-RNA interactions. Analogous to protein-protein interactions or protein-DNA interactions, it is likely that there exist intrinsic sequence-structure motifs common to these RNAs that underlie their binding specificity to specific RBPs. It is known that RNAs form energetically favorable secondary structures, which can be represented as a graph, with nucleotides being nodes and backbone covalent bonds and base-pairing hydrogen bonds representing edges. We hypothesize that these graphs can be mined by graph theory approaches to identify sequence-structure motifs as enriched sub-graphs. In this article, we described the details of this approach, termed RNANetMotif and associated new concepts, namely EKS (Extended K-mer Subgraphs) and GraphK graph search algorithm. To test the utility of our approach, we conducted 3D structure modeling of selected RNA sequences through molecular dynamics (MD) folding simulation and evaluated the significance of the discovered RNA motifs by comparing their spatial exposure with other regions on the RNA. We believe that this approach has the novelty of treating the RNA sequence as a graph and RBP binding sites as enriched subgraph, which has broader applications beyond RBP-RNA interactions.

PERFUMES: pipeline to extract RNA functional motifs and exposed structures

Discovering putative peptides encoded from non-coding RNAs in ribosome profiling data of Arabidopsis thaliana.

Web-based platform for analysis of RNA folding from high throughput chemical probing data

ssHMM: Extracting intuitive sequence-structure motifs from high-throughput RNA-binding protein data

Concurrent prediction of RNA secondary structures with pseudoknots and local 3D motifs in an integer programming framework

RNANetMotif: identifying sequence-structure RNA network motifs in RNA-protein binding sites

Dissecting the energetic architecture within an RNA tertiary structural motif via high-throughput thermodynamic measurements

AnnapuRNA: A scoring function for predicting RNA-small molecule binding poses

A computational approach to finding RNA tertiary motifs in genomic sequences: a case study

DecoyFinder: Identification of Contaminants in Sets of Homologous RNA Sequences

RNA tertiary structure prediction in CASP15 by the GeneSilico group: Folding simulations based on statistical potentials and spatial restraints

RNAscape: geometric mapping and customizable visualization of RNA structure

Genome-Wide Identification of Stable RNA Secondary Structures Across Multiple Organisms Using Chemical Probing Data: Insights into Short Structural Motifs and RNA-Targeting Therapeutics

RNAvigate: efficient exploration of RNA chemical probing datasets

A Computational Approach to Finding RNA Tertiary Motifs in Genomic Sequences

Rsite: a Computational Method to Identify the Functional Sites of Noncoding RNAs

The Multiscale Ernwin/SPQR RNA Structure Prediction Pipeline

RNAelem: an algorithm for discovering sequence-structure motifs in RNA bound by RNA-binding proteins

RNAex: an RNA Secondary Structure Prediction Server Enhanced by High-Throughput Structure-Probing Data.

High-throughput biochemistry in RNA sequence space: predicting structure and function

Scalable and Interpretable Identification of Minimal Undesignable RNA Structure Motifs with Rotational Invariance