Stochastic single-gene autoregulation

A detailed stochastic model of single-gene autoregulation is established and its solutions are explored when mRNA dynamics is fast compared with protein dynamics and in the opposite regime. The model includes all the sources of randomness that are intrinsic to the autoregulation process and it considers both transcriptional and post-transcriptional regulation. The time-scale separation allows the derivation of analytic expressions for the equilibrium distributions of protein and mRNA. These distributions are generally well described in the continuous approximation, which is used to discuss the qualitative features of the protein equilibrium distributions as a function of the biological parameters in the fast mRNA regime. The performance of the time-scale approximation is assessed by comparison with simulations of the full stochastic system, and a good quantitative agreement is found for a wide range of parameter values. We show that either unimodal or bimodal equilibrium protein distributions can arise, and we discuss the autoregulation mechanisms associated with bimodality.