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.

Transcriptome-wide Identification of Single-Stranded RNA Binding Proteins.

Transcriptome-wide Discovery of Coding and Noncoding RNA-binding Proteins.

RNA-Binding Proteome Analysis and Functional Explorations

Identification of whole-cell dsRNA-binding proteins by phase separation

Large-scale Enrichment of RNA-binding Proteins and Binding Site Identification Using the Psoralen Cycloaddition Reaction

Transcriptome-Wide Mapping of Protein-RNA Interactions.

Multiobjective Genome-Wide RNA-Binding Event Identification From CLIP-Seq Data

An RNA tagging approach for system-wide RNA-binding proteome profiling and dynamics investigation upon transcription inhibition

Profiling of RNA-binding protein binding sites by in situ reverse transcription-based sequencing

A large-scale binding and functional map of human RNA-binding proteins

Systematic identification of post-transcriptional regulatory modules

RBPBind: Quantitative Prediction of Protein-RNA Interactions

Large-scale map of RNA-binding protein interactomes across the mRNA life cycle

Capture and Identification of RNA-binding Proteins by Using Click Chemistry-assisted RNA-interactome Capture (CARIC) Strategy

Identification of high-confidence RNA regulatory elements by combinatorial classification of RNA–protein binding sites

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

Capture RIC-seq reveals positional rules of PTBP1-associated RNA loops in splicing regulation

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

lhCLIP reveals the in vivo RNA–RNA interactions recognized by hnRNPK

Global profiling of RNA-binding protein target sites by LACE-seq

Analysis of Sequencing Data for Probing RNA Secondary Structures and Protein-Rna Binding in Studying Posttranscriptional Regulations.