Gene expression noise is affected differentially by feedback in burst frequency and burst size

A stochastic model of autoregulated bursty gene expression by Kumar et al. [Phys. Rev. Lett. 113, 268105 (2014)] has been exactly solved in steady-state conditions under the implicit assumption that protein numbers are sufficiently large such that fluctuations in protein numbers due to reversible protein-promoter binding can be ignored. Here we derive an alternative model that takes into account these fluctuations and hence can be used to study low protein number effects. The...

Mathematical models of gene regulatory networks are widely used to study cell fate changes and transcriptional regulation. When designing such models, it is important to accurately account for sources of stochasticity. However, doing so can be computationally expensive and analytically untractable, posing limits on the extent of our explorations and on parameter inference. Here, we explore this challenge using the example of a simple auto-negative feedback motif, in which we ...

In this article we demonstrate that the so-called bursting production of molecular species during gene expression may be an artifact caused by low time resolution in experimental data collection and not an actual burst in production. We reach this conclusion through an analysis of a two-stage and binary model for gene expression, and demonstrate that in the limit when mRNA degradation is much faster than protein degradation they are equivalent. The negative binomial distribut...

Noise-induced oscillations in individual cells are usually characterized by a non-monotonic power spectrum with an oscillatory autocorrelation function. Here we develop an analytical approach of stochastic oscillations in a minimal stochastic gene expression model including promoter state switching, protein synthesis and degradation, as well as a genetic feedback loop. The autocorrelated function, power spectrum, and probability distribution of protein concentration fluctuati...

Signal-processing molecules inside cells are often present at low copy number, which necessitates probabilistic models to account for intrinsic noise. Probability distributions have traditionally been found using simulation-based approaches which then require estimating the distributions from many samples. Here we present in detail an alternative method for directly calculating a probability distribution by expanding in the natural eigenfunctions of the governing equation, wh...

Expression of many genes varies as a cell transitions through different cell-cycle stages. How coupling between stochastic expression and cell cycle impacts cell-to-cell variability (noise) in the level of protein is not well understood. We analyze a model, where a stable protein is synthesized in random bursts, and the frequency with which bursts occur varies within the cell cycle. Formulas quantifying the extent of fluctuations in the protein copy number are derived and dec...

We show how one may analytically compute the stationary density of the distribution of molecular constituents in populations of cells in the presence of noise arising from either bursting transcription or translation, or noise in degradation rates arising from low numbers of molecules. We have compared our results with an analysis of the same model systems (either inducible or repressible operons) in the absence of any stochastic effects, and shown the correspondence between ...

The canonical model of mRNA expression is the telegraph model, describing a gene that switches on and off, subject to transcription and decay. It describes steady-state mRNA distributions that subscribe to transcription in bursts with first-order decay, referred to as super-Poissonian expression. Using a telegraph-like model, I propose an answer to the question of why gene expression is bursty in the first place, and what benefits it confers. Using analytics for the entropy p...

We study a stochastic model of gene expression, in which protein production has a form of random bursts whose size distribution is arbitrary, whereas protein decay is a first-order reaction. We find exact analytical expressions for the time evolution of the cumulant-generating function for the most general case when both the burst size probability distribution and the model parameters depend on time in an arbitrary (e.g. oscillatory) manner, and for arbitrary initial conditio...

A key goal of systems biology is the predictive mathematical description of gene regulatory circuits. Different approaches are used such as deterministic and stochastic models, models that describe cell growth and division explicitly or implicitly etc. Here we consider simple systems of unregulated (constitutive) gene expression and compare different mathematical descriptions systematically to obtain insight into the errors that are introduced by various common approximations...