Predicting Human mRNA Isoform Levels from Site-Specific Splicing Kinetics in silico

Splicing of pre-mRNA can result in multiple possible mRNA isoforms per gene due to alternative splicing. The frequency at which individual isoforms occur depends on the intrinsic splicing kinetics of the pre-mRNA as well as intracellular chemical conditions. Computational modeling can potentially provide a platform to rapidly assess how variations in intracellular and environmental conditions, for example differential levels of regulatory splicing proteins, affect kinetics and resulting mRNA isoforms. Overcoming the vast combinatoric possibilities of splicing, however, has remained a significant challenge in modeling its kinetics. Here we report the development of a stochastic kinetic model of splicing that is extensible to most protein-coding genes in the human genome. Our model allows for variations in site-specific reaction rates as well as the ability to introduce additional splicing factors. We experimentally validate the predictive capability of our computational model by exploring the spliced isoform ratio of a target gene (SRSF6) under normoxia and hypoxia. This work provides a resource for quantitative, computational analysis of pre-mRNA splicing, allowing for a rapid computational-experimental approach to assess biological hypotheses. Author summarymRNA splicing is a key step in human gene expression with deep impact on the molecular processes determining cell physiology and affecting development and disease. Computational models of splicing are highly attractive to understand life processes but so far have been limited in directly accounting for the chemistry of splicing and are not extensible to most of the human genome. We report here a model that overcomes both of these challenges, providing a computational benchtop to probe splicing kinetics for most genes in the human genome. This model has the capability to rapidly pose biological hypotheses for experimental validation. As a targeted demonstration, we explore the effects of the pathologically-relevant chemical condition hypoxia on the pre-mRNA of a single gene.