Abstract:BACKGROUND:Stochastic genetic switching driven by intrinsic noise is an important process in gene expression. When the rates of gene activation/inactivation are relatively slow, fast, or medium compared with the synthesis/degradation rates of mRNAs and proteins, the variability of protein and mRNA levels may exhibit very different dynamical patterns. It is desirable to provide a systematic approach to identify their key dynamical features in different regimes, aiming at distinguishing which regime a considered gene regulatory network is in from their phenotypic variations.RESULTS:We studied a gene expression model with positive feedbacks when genetic switching rates vary over a wide range. With the goal of providing a method to distinguish the regime of the switching rates, we first focus on understanding the essential dynamics of gene expression system in different cases. In the regime of slow switching rates, we found that the effective dynamics can be reduced to independent evolutions on two separate layers corresponding to gene activation and inactivation states, and the transitions between two layers are rare events, after which the system goes mainly along deterministic ODE trajectories on a particular layer to reach new steady states. The energy landscape in this regime can be well approximated by using Gaussian mixture model. In the regime of intermediate switching rates, we analyzed the mean switching time to investigate the stability of the system in different parameter ranges. We also discussed the case of fast switching rates from the viewpoint of transition state theory. Based on the obtained results, we made a proposal to distinguish these three regimes in a simulation experiment. We identified the intermediate regime from the fact that the strength of cellular memory is lower than the other two cases, and the fast and slow regimes can be distinguished by their different perturbation-response behavior with respect to the switching rates perturbations.CONCLUSIONS:We proposed a simulation experiment to distinguish the slow, intermediate and fast regimes, which is the main point of our paper. In order to achieve this goal, we systematically studied the essential dynamics of gene expression system when the switching rates are in different regimes. Our theoretical understanding provides new insights on the gene expression experiments.

A Simulation Study on Gene Expression Regulation Via Stochastic Model

Exactly solvable models of stochastic gene expression

Stochastic modeling of regulation of gene expression by multiple small RNAs

Rule-Based Modeling and Simulation of Gene Expression Process

Dynamics of gene expression and the regulatory inference problem

Stochastic fluctuations can reveal the feedback signs of gene regulatory networks at the single-molecule level

Stochastic Simulation to Visualize Gene Expression and Error Correction in Living Cells

Stochastic modeling of the mRNA life process: A generalized master equation

From Homogeneity to Heterogeneity: Refining Stochastic Simulations of Gene Regulation

Multi-species stochastic models: small RNA regulation of gene expression

Distinguishing the Rates of Gene Activation from Phenotypic Variations

Stochastic equations for a self-regulating gene

Gene expression dynamics with stochastic bursts: exact results for a coarse-grained model

Learning a nonlinear dynamical system model of gene regulation: A perturbed steady-state approach

Validation of an algorithm for delay stochastic simulation of transcription and translation in prokaryotic gene expression

Deterministic and stochastic descriptions of gene expression dynamics

Stochastic single-gene auto-regulation

On a stochastic gene expression with pre-mRNA, mRNA and protein contribution

A Machine Learning Approach to Simulate Gene Expression and Infer Gene Regulatory Networks

Sources of stochasticity in constitutive and autoregulated gene expression

Inferring gene regulatory networks from single-cell data: a mechanistic approach