The role of multiple marks in epigenetic silencing and the emergence of a stable bivalent chromatin state

We introduce and analyze a minimal model of epigenetic silencing in budding yeast, built upon known biomolecular interactions in the system. Doing so, we identify the epigenetic marks essential for the bistability of epigenetic states. The model explicitly incorporates two key chromatin marks, namely H4K16 acetylation and H3K79 methylation, and explores whether the presence of multiple marks lead to a qualitatively different systems behavior. We find that having both modifications is important for the robustness of epigenetic silencing. Besides the silenced and transcriptionally active fate of chromatin, our model leads to a novel state with bivalent (i.e., both active and silencing) marks under certain perturbations (knock-out mutations, inhibition or enhancement of enzymatic activity). The bivalent state appears under several perturbations and is shown to result in patchy silencing. We also show that the titration effect, owing to a limited supply of silencing proteins, can result in counter-intuitive responses. The design principles of the silencing system is systematically investigated and disparate experimental observations are assessed within a single theoretical framework. Specifically, we discuss the behavior of Sir protein recruitment, spreading and stability of silenced regions in commonly-studied mutants (e.g., sas2, dot1) illuminating the controversial role of Dot1 in the systems biology of yeast silencing.

💡 Research Summary

The paper presents a minimal, mathematically tractable model of epigenetic silencing in budding yeast that explicitly incorporates two well‑characterized histone modifications: H4K16 acetylation (Ac) and H3K79 methylation (Me). Building on known biochemical interactions, the authors formulate a stochastic Markov‑chain framework in which each nucleosomal unit can occupy one of three states—silenced, active, or bivalent—depending on the local balance of Ac and Me and the availability of Sir proteins. Sir complexes are assumed to be present in limited total quantity, introducing a titration effect: excessive Sir binding at one locus reduces the pool available for spreading to neighboring loci.

The model demonstrates that the coexistence of both marks is essential for robust bistability. Ac acts as a negative feedback, preventing Sir binding, whereas Me provides positive feedback that promotes Sir recruitment and spreading. When both feedbacks operate together, the system exhibits two deep attractors (silenced and active) that are resistant to perturbations. Parameter sweeps reveal that reducing Ac (sas2 deletion) leads to over‑binding of Sir and uncontrolled expansion of silenced domains, while reducing Me (dot1 deletion) impairs Sir propagation, rendering silencing unstable and allowing active chromatin to invade.

A striking prediction is the emergence of a third, “bivalent” state under specific perturbations. In this state, Ac and Me coexist on the same nucleosome, allowing only partial Sir occupancy. The result is a patchy, heterogeneous silencing pattern that matches experimental observations of mixed chromatin marks in certain mutant backgrounds. The bivalent state also appears when enzymatic activities are artificially suppressed or enhanced, illustrating that the system’s qualitative behavior is not solely determined by the presence of a single mark but by the dynamic balance between them.

The titration effect further explains counter‑intuitive experimental findings. Because Sir proteins are limiting, increasing their total concentration can paradoxically reduce the overall spread of silencing: excess Sir saturates high‑affinity sites, preventing redistribution to distal loci. Conversely, lowering Sir levels can free the pool for more uniform spreading, expanding silenced regions. This mechanism reconciles reports of Sir over‑expression diminishing silencing with classic models that predict the opposite.

The authors systematically test the model against a range of yeast mutants (sas2Δ, dot1Δ, double mutants, Sir over‑expression, and pharmacological inhibition of Dot1). In each case, the simulated outcomes reproduce key phenotypes: loss of bistability, excessive silencing, or the appearance of bivalent chromatin. The analysis clarifies the controversial role of Dot1, showing that its methyltransferase activity is not merely a silencing antagonist but a necessary component for stabilizing Sir‑mediated repression under conditions where Ac is present.

Overall, the study provides a unified theoretical framework that integrates multiple epigenetic marks, protein scarcity, and feedback loops to explain the robustness and flexibility of yeast silencing. It highlights design principles—dual‑mark regulation and limited silencing factor supply—that likely extend to more complex eukaryotic systems and offers a basis for engineering synthetic epigenetic circuits with predictable bistable or multistable behavior.