Abstract
Motivation: RNA binding proteins (RBPs) play important roles in post-transcriptional control of gene expression, including splicing, transport, polyadenylation and RNA stability. To model protein–RNA interactions by considering all available sources of information, it is necessary to integrate the rapidly growing RBP experimental data with the latest genome annotation, gene function, RNA sequence and structure. Such integration is possible by matrix factorization, where current approaches have an undesired tendency to identify only a small number of the strongest patterns with overlapping features. Because protein–RNA interactions are orchestrated by multiple factors, methods that identify discriminative patterns of varying strengths are needed.Results: We have developed an integrative orthogonality-regularized nonnegative matrix factorization (iONMF) to integrate multiple data sources and discover non-overlapping, class-specific RNA binding patterns of varying strengths. The orthogonality constraint halves the effective size of the factor model and outperforms other NMF models in predicting RBP interaction sites on RNA. We have integrated the largest data compendium to date, which includes 31 CLIP experiments on 19 RBPs involved in splicing (such as hnRNPs, U2AF2, ELAVL1, TDP-43 and FUS) and processing of 3’UTR (Ago, IGF2BP). We show that the integration of multiple data sources improves the predictive accuracy of retrieval of RNA binding sites. In our study the key predictive factors of protein–RNA interactions were the position of RNA structure and sequence motifs, RBP co-binding and gene region type. We report on a number of protein-specific patterns, many of which are consistent with experimentally determined properties of RBPs.Availability and implementation: The iONMF implementation and example datasets are available at https://github.com/mstrazar/ionmf.Contact: tomaz.curk@fri.uni-lj.siSupplementary information: Supplementary data are available at Bioinformatics online.
Highlights
Reference and details about all experiments used in the study are listed
Depending on the experimental protocol used (PARCLIP, CLIPSEQ, iCLIP, HITSCLIP) we report number of cross-linking clusters and number of individual sites for each experiment used
>HSU14570 Human Alu-Sb2 subfamily consensus sequence. (rev. complement) TTTTTTTTGAGACGGAGTCTCGCTCTGTCGCCCAGGCCGGACTGCGGACTGCAGTGGCGCAATCTCGGCTCACTGCAAGCT TCCGCTTCCCGGGTTCACGCCATTCTCCTGCCTCAGCCTCCCCAGTAGCTGGGACTACAGGCGCCCGCCACCGCGCCCGG CTAATTTTTTGTATTTTTAGTAGAGACGGGGTTTCACCTTGTTAGCCAGGATGGTCTCGATCTCCTGACCTCATGATCCAC CCCGCCTCGGCCTCCCAAAGTGCTGGGATTACAGGCGTGAGCCACCGCGCCCGGCC
Summary
The general matrix factorization problem can be solved with different optimization approaches; these include:. Non-convexity of the problem follows from observing that there exist equivalent solutions X = WUT UH, where U is any unitary matrix of appropriate size. To learn a factor model with iONMF we propose the following optimization problem with respect to W and Hi for i = 1, ..., N : N minW,H − 2Hi) − λi λi = −XTi W + HiWT W + α(2HiHTi Hi − Hi) To satisfy the Karush-Kuhn-Tucker optimality conditions at a stationary point we must have: which leads to the following update rules: Hi ◦ λi = 0 H2i ◦ (λ+i − λ−i ) = 0. This is exactly the update rule in Equation 3 (see main text)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.