Stochastic epigenetic dynamics of gene switching

Epigenetic modifications of histones crucially affect eukaryotic gene activity, while the epigenetic histone state is largely determined by the binding of specific factors such as the transcription factors (TFs) to DNA. Here, the way in which the TFs and the histone state are dynamically correlated is not obvious when the TF synthesis is regulated by the histone state. This type of feedback regulatory relation is ubiquitous in gene networks to determine cell fate in differentiation and other cell transformations. To gain insights into such dynamical feedback regulations, we theoretically analyze a model of epigenetic gene switching by extending the Doi-Peliti operator formalism of reaction kinetics to the problem of coupled molecular processes. Spin-1 and spin-1/2 coherent-state representations are introduced to describe stochastic reactions of histones and binding or unbinding of TFs in a unified way, which provides a concise view of the effects of timescale difference among these molecular processes; even in the case that binding or unbinding of TFs to or from DNA is adiabatically fast, the slow nonadiabatic histone dynamics gives rise to a distinct circular flow of the probability flux around basins in the landscape of the gene state distribution, which leads to hysteresis in gene switching. In contrast to the general belief that the change in the amount of TF precedes the histone state change, flux drives histones to be modified prior to the change in the amount of TF in self-regulating circuits. Flux-landscape analyses shed light on the nonlinear nonadiabatic mechanism of epigenetic cell fate decision making.