Before revision: The choice between deterministic and stochastic modeling can lead to dramatically different theoretical conclusions regarding the steady state behavior of a gene regulatory network (GRN). This is particularly interesting when low-molecular counts and slow TF-gene binding /unbinding lead to the emergence of new phenotypes in the stochastic model that are not reflected in the corresponding deterministic model. This work uncovers a mechanism that underlies this emergence of multiple modes under slow slow promoter kinetics, and studies it theoretically.
After revision: Phenotypical variability in the absence of genetic variation often reflects complex energetic landscapes associated with underlying gene regulatory networks (GRNs). In this view, different phenotypes are associated with alternative states of complex nonlinear systems: stable attractors in deterministic models or modes of stationary distributions in stochastic descriptions. We provide theoretical and practical characterizations of these landscapes, specifically focusing on stochastic slow promoter kinetics, a time scale relevant when transcription factor binding and unbinding are affected by epigenetic processes like DNA methylation and chromatin remodeling. In this case, largely unexplored except for numerical simulations, adiabatic approximations of promoter kinetics are not appropriate. In contrast to the existing literature, we provide rigorous analytic characterizations of multiple modes under slow slow promoter kinetics, and studies it theoretically.