Abstract:A toggle switch consists of two genes that mutually repress each other. This regulatory motif is active during cell differentiation and is thought to act as a memory device, being able to choose and maintain cell fate decisions. Commonly, this switch has been modeled in a deterministic framework where transcription and translation are lumped together. In this description, bistability occurs for transcription factor cooperativity, whereas autoactivation leads to a tristable system with an additional undecided state. In this contribution, we study the stability and dynamics of a two-stage gene expression switch within a probabilistic framework inspired by the properties of the Pu/Gata toggle switch in myeloid progenitor cells. We focus on low mRNA numbers, high protein abundance, and monomeric transcription-factor binding. Contrary to the expectation from a deterministic description, this switch shows complex multiattractor dynamics without autoactivation and cooperativity. Most importantly, the four attractors of the system, which only emerge in a probabilistic two-stage description, can be identified with committed and primed states in cell differentiation. To begin, we study the dynamics of the system and infer the mechanisms that move the system between attractors using both the quasipotential and the probability flux of the system. Next, we show that the residence times of the system in one of the committed attractors are geometrically distributed. We derive an analytical expression for the parameter of the geometric distribution, therefore completely describing the statistics of the switching process and elucidate the influence of the system parameters on the residence time. Moreover, we find that the mean residence time increases linearly with the mean protein level. This scaling also holds for a one-stage scenario and for autoactivation. Finally, we study the implications of this distribution for the stability of a switch and discuss the influence of the stability on a specific cell differentiation mechanism. Our model explains lineage priming and proposes the need of either high protein numbers or long-term modifications such as chromatin remodeling to achieve stable cell fate decisions. Notably, we present a system with high protein abundance that nevertheless requires a probabilistic description to exhibit multistability, complex switching dynamics, and lineage priming.

Inter-state switching in stochastic gene expression: Exact solution, an adiabatic limit and oscillations in molecular distributions

Stochastic phenotype transition of a single cell in an intermediate region of gene-state switching

Stability and Multiattractor Dynamics of a Toggle Switch Based on a Two-Stage Model of Stochastic Gene Expression

Approximation scheme based on effective interactions for stochastic gene regulation

Stochastic gene expression as a many body problem

Bistable switching asymptotics for the self regulating gene

Exactly solvable models of stochastic gene expression

Statistical physics of a model binary genetic switch with linear feedback

Metastable behavior in Markov processes with internal states

Stochastic Switching in Gene Networks Can Occur by a Single-Molecule Event or Many Molecular Steps.

Stochastic switching of delayed feedback suppresses oscillations in genetic regulatory systems

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

Intrinsic fluctuations in stochastic delay systems: theoretical description and application to a simple model of gene regulation

Efficient analysis of stochastic gene dynamics in the non-adiabatic regime using piecewise deterministic Markov processes

Stochastic gene expression with delay

Stochastic Bloch-Redfield theory: Quantum jumps in a solid-state environment

Exact solution of stochastic gene expression models with bursting, cell cycle and replication dynamics

Fluctuating-rate model with multiple gene states

Eliminating fast reactions in stochastic simulations of biochemical networks: a bistable genetic switch

Distinguishing the Rates of Gene Activation from Phenotypic Variations

Transient hysteresis and inherent stochasticity in gene regulatory networks