Connecting protein and mRNA burst distributions for stochastic models of gene expression

Transcriptional pulsing has been observed in both prokaryotes and eukaryotes and plays a crucial role in cell to cell variability of protein and mRNA numbers. The issue is how the time constants associated with episodes of transcriptional bursting impact cellular mRNA and protein distributions and reciprocally, to what extent experimentally observed distributions can be attributed to transcriptional pulsing. We address these questions by investigating the exact time-dependent...

The intrinsic stochasticity of gene expression can lead to large variations in protein levels across a population of cells. To explain this variability, different sources of mRNA fluctuations ('Poisson' and 'Telegraph' processes) have been proposed in stochastic models of gene expression. Both Poisson and Telegraph scenario models explain experimental observations of noise in protein levels in terms of 'bursts' of protein expression. Correspondingly, there is considerable int...

Inside individual cells, expression of genes is stochastic across organisms ranging from bacterial to human cells. A ubiquitous feature of stochastic expression is burst-like synthesis of gene products, which drives considerable intercellular variability in protein levels across an isogenic cell population. One common mechanism by which cells control such stochasticity is negative feedback regulation, where a protein inhibits its own synthesis. For a single gene that is expre...

Due to the stochastic nature of biochemical processes, the copy number of any given type of molecule inside a living cell often exhibits large temporal fluctuations. Here, we develop analytic methods to investigate how the noise arising from a bursting input is reshaped by a transport reaction which is either linear or of the Michaelis-Menten type. A slow transport rate smoothes out fluctuations at the output end and minimizes the impact of bursting on the downstream cellular...

The processes, resulting in the transcription of RNA, are intrinsically noisy. It was observed experimentally that the synthesis of mRNA molecules is driven by short, burst-like, events. An accurate prediction of the protein level often requires one to take these fluctuations into account. Here, we consider the stochastic model of gene expression regulated by small RNAs. Small RNA post-transcriptional regulation is achieved by base-pairing with mRNA. We show that in a strong ...

Many of the existing stochastic models of gene expression contain the first-order decay reaction term that may describe active protein degradation or dilution. If the model variable is interpreted as the molecule number, and not concentration, the decay term may also approximate the loss of protein molecules due to cell division as a continuous degradation process. The seminal model of that kind leads to gamma distributions of protein levels, whose parameters are defined by t...

The intrinsic stochasticity of gene expression can lead to large variability of protein levels across a population of cells. Variability (or noise) in protein distributions can be modulated by cellular mechanisms of gene regulation; in particular, there is considerable interest in understanding the role of post-transcriptional regulation. To address this issue, we propose and analyze a stochastic model for post-transcriptional regulation of gene expression. The analytical sol...

This paper considers adiabatic reduction in both discrete and continuous models of stochastic gene expression. In gene expression models, the concept of bursting is a production of several molecules simultaneously and is generally represented as a compound Poisson process of random size. In a general two-dimensional birth and death discrete model, we prove that under specific assumptions and scaling (that are characteristics of the mRNA-protein system) an adiabatic reduction ...

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...

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...