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.

Stochastic equations for a self-regulating gene

Bistable switching asymptotics for the self regulating gene

Self-regulating genes. Exact steady state solution by using Poisson representation

Stochastic single-gene auto-regulation

A Stochastic Analysis of Autoregulation of Gene Expression

Stochastic gene expression with delay

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

Dynamics of gene expression and the regulatory inference problem

Adiabatic and Non-Adiabatic Non-Equilibrium Stochastic Dynamics of Single Regulating Genes

A Simulation Study on Gene Expression Regulation Via Stochastic Model

Exactly solvable models of stochastic gene expression

Exact switching time distributions for autoregulated gene expression models with mRNA and protein descriptions

Stochastic modeling of regulation of gene expression by multiple small RNAs

Stochastic Analysis Of An Incoherent Feedforward Genetic Motif

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

Dynamics of gene expression under feedback

Stochastic mRNA production by a three-state gene

Small protein number effects in stochastic models of autoregulated bursty gene expression

Oscillations in the expression of a self-repressed gene induced by a slow transcriptional dynamics

Single-cell stochastic gene expression kinetics with coupled positive-plus-negative feedback

Stochastic timing in gene expression for simple regulatory strategies