Special Function Methods for Bursty Models of Transcription

Abstract

We explore a Markov model used in the analysis of gene expression, involving the bursty production of pre-mRNA, its conversion to mature mRNA, and its consequent degradation. We demonstrate that the integration used to compute the solution of the stochastic system can be approximated by the evaluation of special functions. Furthermore, the form of the special function solution generalizes to a broader class of burst distributions. In light of the broader goal of biophysical parameter inference from transcriptomics data, we apply the method to simulated data, demonstrating effective control of precision and runtime. Finally, we suggest a non-Bayesian approach to reducing the computational complexity of parameter inference to linear order in state space size and number of candidate parameters.

arXiv commentBody: 15 pages, 2 figures, 2 tables. Supplement: 10 pages, 1 figure

References in corpus (6)

• Analytical distributions for stochastic gene expression
• An introduction to infinite HMMs for single molecule data analysis
• Accurate Chemical Master Equation Solution Using Multi-Finite Buffers
• General Relativity and Quantum Mechanics: Towards a Generalization of the Lambert W Function
• Effect of reaction step-size noise on the switching dynamics of stochastic populations
• Taylor series for generalized Lambert W functions

Cited by in corpus (2)

• Inferring stochastic rates from heterogeneous snapshots of particle positions
• Analytic solution of chemical master equations involving gene switching. I: Representation theory and diagrammatic approach to exact solution