The dynamics of simple gene-network motifs subject to extrinsic fluctuations

Cellular processes do not follow deterministic rules; even in identical environments genetically identical cells can make random choices leading to different phenotypes. This randomness originates from fluctuations present in the biomolecular interaction networks. Most previous work has been focused on the intrinsic noise (IN) of these networks. Yet, especially for high-copy-number biomolecules, extrinsic or environmental noise (EN) has been experimentally shown to dominate the variation. Here, we develop an analytical formalism that allows for calculation of the effect of EN on gene-expression motifs. We introduce a method for modeling bounded EN as an auxiliary species in the master equation. The method is fully generic and is not limited to systems with small EN magnitudes. We focus our study on motifs that can be viewed as the building blocks of genetic switches: a nonregulated gene, a self-inhibiting gene, and a self-promoting gene. The role of the EN properties (magnitude, correlation time, and distribution) on the statistics of interest are systematically investigated, and the effect of fluctuations in different reaction rates is compared. Due to its analytical nature, our formalism can be used to quantify the effect of EN on the dynamics of biochemical networks and can also be used to improve the interpretation of data from single-cell gene-expression experiments.

💡 Research Summary

This paper presents a comprehensive analytical framework for quantifying the impact of extrinsic (environmental) noise on simple gene‑regulatory motifs. While most prior studies have focused on intrinsic noise arising from the stochasticity of biochemical reactions, experimental evidence shows that for high‑copy‑number molecules extrinsic fluctuations often dominate. To address this, the authors model bounded extrinsic noise as an auxiliary species within the master equation, allowing arbitrary noise magnitude, correlation time, and distribution (including non‑Gaussian statistics).

The starting point is a minimal gene‑expression model with protein production at rate F(n) and degradation at rate ν n. The intrinsic‑noise‑only master equation is solved analytically, and a WKB (large‑N) approximation yields a Hamilton‑Jacobi formulation with Hamiltonian H(x,pₓ)=f(x)(e^{pₓ}−1)+x(e^{−pₓ}−1).

Extrinsic noise is introduced by letting a kinetic parameter (e.g., the degradation rate) fluctuate as ξ(t) with ⟨ξ⟩=1, variance σ²_ex, and exponential autocorrelation τ_c. ξ is generated by a fast mRNA‑protein auxiliary circuit whose stationary distribution is negative‑binomial (or, for large copy numbers, gamma). The auxiliary protein count k=K ξ defines a rescaled noise variable \tilde ξ=ρ ξ (ρ=K/N). After adiabatically eliminating the short‑lived auxiliary mRNA, the combined system reduces to a two‑dimensional Hamiltonian (Eq. 8) that captures both intrinsic and extrinsic stochasticity.

Two analytically tractable limits are examined:

1. White‑noise limit (τ_c → 0) – Fast extrinsic fluctuations are treated as an effective white noise. Solving the Hamilton‑Jacobi equation yields an explicit expression for the momentum pₓ (Eq. 9). Integrating pₓ gives the action S(x) and, via the Gaussian approximation, the stationary probability density. The variance of the protein distribution becomes
    σ²_obs = N x* (1 + x* V τ_c) /