Mathematical analysis of the limiting behaviors of a chromatin modification circuit

In the last decade, the interactions among histone modifications and DNA methylation and their effect on the DNA structure, i.e., chromatin state, have been identified as key mediators for the maintenance of cell identity, defined as epigenetic cell memory. In this paper, we determine how the positive feedback loops generated by the auto- and cross-catalysis among repressive modifications affect the temporal duration of the cell identity. To this end, we conduct a stochastic analysis of a recently published chromatin modification circuit considering two limiting behaviors: fast erasure rate of repressive histone modifications or fast erasure rate of DNA methylation. In order to perform this mathematical analysis, we first show that the deterministic model of the system is a singular singularly perturbed (SSP) system and use a model reduction approach for SSP systems to obtain a reduced one-dimensional model. We thus analytically evaluate the reduced system’s stationary probability distribution and the mean switching time between active and repressed chromatin states. We then add a computational study of the original reaction model to validate and extend the analytical findings. Our results show that the absence of DNA methylation reduces the bias of the system’s stationary probability distribution toward the repressed chromatin state and the temporal duration of this state’s memory. In the absence of repressive histone modifications, we also observe that the time needed to reactivate a repressed gene with an activating input is less stochastic, suggesting that repressive histone modifications specifically contribute to the highly variable latency of state reactivation.