Dataset Information

Predicting sequence and structural specificities of RNA binding regions recognized by splicing factor SRSF1.

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.

SUBMITTER: Wang X

PROVIDER: S-EPMC3287504 | biostudies-literature | 2011 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Predicting sequence and structural specificities of RNA binding regions recognized by splicing factor SRSF1.

Wang Xin X Juan Liran L Lv Junjie J Wang Kejun K Sanford Jeremy R JR Liu Yunlong Y

BMC genomics 20111223

<h4>Background</h4>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 prot ...[more]

PMID: 22369183

Similar Datasets

Project description:A combination of cDNA sequencing of the complete coding region, protein comparisons, binding site mapping, and electron microscopic imaging has permitted the formulation of a structural model of cytotactin. Cytotactin is a large extracellular matrix glycoprotein that displays a restricted tissue distribution during development. Although there appears to be a single cytotactin gene, multiple cytotactin polypeptides and mRNAs are detected in a variety of tissues. We report here the sequences and relationships of cDNAs that encode the complete amino acid sequences of two cytotactin polypeptides in chicken brain. The translated cDNA sequences agree with those obtained by direct analysis of cytotactin and fragments of the molecule. All regions of the polypeptides appear to be identical except for a 273 amino acid segment found in the larger but not in the smaller. At their amino termini, both polypeptides contain a cysteine-rich segment that probably includes those residues that link monomers into hexamers. This segment is followed by 13 epidermal growth factor-like (EGFL) repeats and then 8 consecutive segments that each resemble the type III repeats found in fibronectin. At their carboxyl termini, the polypeptides are similar to the beta and gamma chains of fibrinogen, including a calcium-binding segment. The additional sequence in the large polypeptide is inserted after the fifth type III repeat and includes three additional type III repeats. On RNA transfer blot analyses, cytotactin cDNA probes detected a 6.4-kilobase (kb) component in both brain and gizzard and larger mRNAs in both tissues, but those in gizzard were larger by about 1 kb than those in brain. A probe specific to the insert did not hybridize to the 6.4-kb mRNA in either tissue but detected the larger mRNAs in both tissues. At least a portion of the insert is thus present in both tissues, but there may be additional inserts in the gizzard mRNAs. The proposed model of cytotactin specifies the orientation of the polypeptides, the localization of interchain disulfide bonds, the structural elements constituting the thin and thick segments (EGFL repeats and type III repeats, respectively), the terminal fibrinogen-like nodular region, and the relative location of the cell-binding region.

Project description:BackgroundMotivated by the increased amount of data on protein-RNA interactions and the availability of complete genome sequences of several organisms, many computational methods have been proposed to predict binding sites in protein-RNA interactions. However, most computational methods are limited to finding RNA-binding sites in proteins instead of protein-binding sites in RNAs. Predicting protein-binding sites in RNA is more challenging than predicting RNA-binding sites in proteins. Recent computational methods for finding protein-binding sites in RNAs have several drawbacks for practical use.ResultsWe developed a new support vector machine (SVM) model for predicting protein-binding regions in mRNA sequences. The model uses sequence profiles constructed from log-odds scores of mono- and di-nucleotides and nucleotide compositions. The model was evaluated by standard 10-fold cross validation, leave-one-protein-out (LOPO) cross validation and independent testing. Since actual mRNA sequences have more non-binding regions than protein-binding regions, we tested the model on several datasets with different ratios of protein-binding regions to non-binding regions. The best performance of the model was obtained in a balanced dataset of positive and negative instances. 10-fold cross validation with a balanced dataset achieved a sensitivity of 91.6%, a specificity of 92.4%, an accuracy of 92.0%, a positive predictive value (PPV) of 91.7%, a negative predictive value (NPV) of 92.3% and a Matthews correlation coefficient (MCC) of 0.840. LOPO cross validation showed a lower performance than the 10-fold cross validation, but the performance remains high (87.6% accuracy and 0.752 MCC). In testing the model on independent datasets, it achieved an accuracy of 82.2% and an MCC of 0.656. Testing of our model and other state-of-the-art methods on a same dataset showed that our model is better than the others.ConclusionsSequence profiles of log-odds scores of mono- and di-nucleotides were much more powerful features than nucleotide compositions in finding protein-binding regions in RNA sequences. But, a slight performance gain was obtained when using the sequence profiles along with nucleotide compositions. These are preliminary results of ongoing research, but demonstrate the potential of our approach as a powerful predictor of protein-binding regions in RNA. The program and supporting data are available at http://bclab.inha.ac.kr/RBPbinding .

Dataset Information

Predicting sequence and structural specificities of RNA binding regions recognized by splicing factor SRSF1.

Background

Results

Conclusion

Publications

Predicting sequence and structural specificities of RNA binding regions recognized by splicing factor SRSF1.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets