Stochastic and deterministic simulations of a delayed genetic oscillation model: Investigating the validity of reductions

Abstract

Quasi-stationary approximations are commonly used in order to simplify and reduce the number of equations of genetic circuit models. Protein/protein and protein/DNA binding reactions are considered to occur on much shorter time scale than protein production and degradation processes and often tacitly assumed at a quasi-equilibrium. Taking a biologically inspired, typical, small, abstract, negative feedback, genetic circuit model as study case, we investigate in this paper how different quasi-stationary approximations change the system behaviour both in deterministic and stochastic frameworks. We investigate the consistence between the deterministic and stochastic behaviours of our time-delayed negative feedback genetic circuit model with different implementations of quasi-stationary approximations. Quantitative and qualitative differences are observed among the various reduction schemes and with the underlying microscopic model, for biologically reasonable ranges and combinations of the microscopic model kinetic rates. The different reductions do not behave in the same way: correlations and amplitudes of the stochastic oscillations are not equally captured and the population behaviour is not always in consistence with the deterministic curves.