Abstract:Background RNA-binding proteins (RBPs) play diverse roles in eukaryotic RNA processing. Despite their pervasive functions in coding and noncoding RNA biogenesis and regulation, elucidating the sequence specificities that define protein-RNA interactions remains a major challenge. Recently, CLIP-seq (Cross-linking immunoprecipitation followed by high-throughput sequencing) has been successfully implemented to study the transcriptome-wide binding patterns of SRSF1, PTBP1, NOVA and fox2 proteins. These studies either adopted traditional methods like Multiple EM for Motif Elicitation (MEME) to discover the sequence consensus of RBP's binding sites or used Z-score statistics to search for the overrepresented nucleotides of a certain size. We argue that most of these methods are not well-suited for RNA motif identification, as they are unable to incorporate the RNA structural context of protein-RNA interactions, which may affect to binding specificity. Here, we describe a novel model-based approach-- RNAMotifModeler to identify the consensus of protein-RNA binding regions by integrating sequence features and RNA secondary structures. Results As an example, we implemented RNAMotifModeler on SRSF1 (SF2/ASF) CLIP-seq data. The sequence-structural consensus we identified is a purine-rich octamer 'AGAAGAAG' in a highly single-stranded RNA context. The unpaired probabilities, the probabilities of not forming pairs, are significantly higher than negative controls and the flanking sequence surrounding the binding site, indicating that SRSF1 proteins tend to bind on single-stranded RNA. Further statistical evaluations revealed that the second and fifth bases of SRSF1octamer motif have much stronger sequence specificities, but weaker single-strandedness, while the third, fourth, sixth and seventh bases are far more likely to be single-stranded, but have more degenerate sequence specificities. Therefore, we hypothesize that nucleotide specificity and secondary structure play complementary roles during binding site recognition by SRSF1. Conclusion In this study, we presented a computational model to predict the sequence consensus and optimal RNA secondary structure for protein-RNA binding regions. The successful implementation on SRSF1 CLIP-seq data demonstrates great potential to improve our understanding on the binding specificity of RNA binding proteins.

SPOT-Seq-RNA: Predicting Protein-Rna Complex Structure and RNA-binding Function by Fold Recognition and Binding Affinity Prediction.

Highly accurate and high-resolution function prediction of RNA binding proteins by fold recognition and binding affinity prediction.

Prediction and Validation of the Unexplored RNA‐binding Protein Atlas of the Human Proteome

Structure-Based Prediction of Rna-Binding Domains and Rna-Binding Sites and Application to Structural Genomics Targets

Seq-RBPPred: Predicting RNA-Binding Proteins from Sequence.

Prediction of RNA binding proteins comes of age from low resolution to high resolution.

RBPBind: Quantitative Prediction of Protein-RNA Interactions

Predicting DNA-binding Proteins and Binding Residues by Complex Structure Prediction and Application to Human Proteome.

RBPsuite: RNA-protein binding sites prediction suite based on deep learning

Predicting Sequence and Structural Specificities of RNA Binding Regions Recognized by Splicing Factor SRSF1

Efficient mapping of RNA‐binding residues in RNA‐binding proteins using local sequence features of binding site residues in protein‐RNA complexes

Identification of RNA-binding Sites in Proteins by Integrating Various Sequence Information.

RBind: computational network method to predict RNA binding sites.

Analysis and Prediction of RNA-binding Residues Using Sequence, Evolutionary Conservation, and Predicted Secondary Structure and Solvent Accessibility.

A Novel Method for Predicting RNA-interacting Residues in Proteins Using a Combination of Feature-Based and Sequence Template-Based Methods

Deep Neural Networks for Inferring Binding Sites of RNA-binding Proteins by Using Distributed Representations of RNA Primary Sequence and Secondary Structure.

RBinds: A User-Friendly Server for RNA Binding Site Prediction

RBRDetector: improved prediction of binding residues on RNA-binding protein structures using complementary feature- and template-based strategies.

A Method for Predicting Rna-Protein Interaction and Interaction Sites

RPI-Pred: predicting ncRNA-protein interaction using sequence and structural information

Predicting Dynamic Cellular Protein–rna Interactions by Deep Learning Using in Vivo RNA Structures