Stochastic Gene Expression in Cells: A Point Process Approach

TL;DR

This paper introduces a new point process approach to model stochastic gene expression, removing the exponential assumption and revealing that non-exponential step durations can significantly increase protein number variance.

Contribution

It develops a marked Poisson point process framework for gene expression, enabling analysis without exponential duration assumptions, and provides new formulas for variance under more realistic conditions.

Findings

Classical models may underestimate protein variance.

Non-exponential step durations can significantly increase variance.

The new approach generalizes existing models and reveals counter-intuitive results.

Abstract

This paper investigates the stochastic fluctuations of the number of copies of a given protein in a cell. This problem has already been addressed in the past and closed-form expressions of the mean and variance have been obtained for a simplified stochastic model of the gene expression. These results have been obtained under the assumption that the duration of all the protein production steps are exponentially distributed. In such a case, a Markovian approach (via Fokker-Planck equations) is used to derive analytic formulas of the mean and the variance of the number of proteins at equilibrium. This assumption is however not totally satisfactory from a modeling point of view since the distribution of the duration of some steps is more likely to be Gaussian, if not almost deterministic. In such a setting, Markovian methods can no longer be used. A finer characterization of the fluctuations of the number of proteins is therefore of primary interest to understand the general economy of the cell. In this paper, we propose a new approach, based on marked Poisson point processes, which allows to remove the exponential assumption. This is applied in the framework of the classical three stages models of the literature: transcription, translation and degradation. The interest of the method is shown by recovering the classical results under the assumptions that all the durations are exponentially distributed but also by deriving new analytic formulas when some of the distributions are not anymore exponential. Our results show in particular that the exponential assumption may, surprisingly, underestimate significantly the variance of the number of proteins when some steps are in fact not exponentially distributed. This counter-intuitive result stresses the importance of the statistical assumptions in the protein production process.