Dimensionality reduction via path integration for computing mRNA distributions

内容摘要

Inherent stochasticity in gene expression leads to distributions of mRNA copy numbers in a population of identical cells. These distributions are determined primarily by the multitude of states of a gene promoter, each driving transcription at a different rate. In an era where single-cell mRNA copy number data are more and more available, there is an increasing need for fast computations of mRNA distributions. In this paper, we present a method for computing separate distributions for each species of mRNA molecules, i. e. mRNAs that have been either partially or fully processed post-transcription. The method involves the integration over all possible realizations of promoter states, which we cast into a set of linear ordinary differential equations of dimension $M\times n_j$, where $M$ is the number of available promoter states and $n_j$ is the mRNA copy number of species $j$ up to which one wishes to compute the probability distribution. This approach is superior to solving the Master equation (ME) directly in two ways: a) the number of coupled differential equations in the ME approach is $M\timesΛ_1\timesΛ_2\times ...\timesΛ_L$, where $Λ_j$ is the cutoff for the probability of th