Boolean network model predicts cell cycle sequence of fission yeast

A Boolean network model of the cell-cycle regulatory network of fission yeast (Schizosaccharomyces Pombe) is constructed solely on the basis of the known biochemical interaction topology. Simulating the model in the computer, faithfully reproduces the known sequence of regulatory activity patterns along the cell cycle of the living cell. Contrary to existing differential equation models, no parameters enter the model except the structure of the regulatory circuitry. The dynamical properties of the model indicate that the biological dynamical sequence is robustly implemented in the regulatory network, with the biological stationary state G1 corresponding to the dominant attractor in state space, and with the biological regulatory sequence being a strongly attractive trajectory. Comparing the fission yeast cell-cycle model to a similar model of the corresponding network in S. cerevisiae, a remarkable difference in circuitry, as well as dynamics is observed. While the latter operates in a strongly damped mode, driven by external excitation, the S. pombe network represents an auto-excited system with external damping.

💡 Research Summary

The paper presents a Boolean network model of the cell‑cycle regulatory circuitry in the fission yeast Schizosaccharomyces pombe, constructed solely from the known topology of biochemical interactions. Each node in the network represents a key regulatory protein or gene (such as Cdc2, Cdc13, Rum1, Ste9, etc.) and can occupy one of two discrete states: active (1) or inactive (0). Logical update rules are derived directly from experimentally established activation and inhibition relationships, eliminating the need for kinetic parameters that are required in ordinary differential equation (ODE) models.

Simulation of the model proceeds by iterating the Boolean update functions from a randomly chosen initial configuration. The state space, comprising 2^N possible configurations (where N is the number of nodes), is explored exhaustively. The results show that the overwhelming majority of initial states converge to a single fixed point that corresponds to the biologically observed G1 phase. This fixed point functions as the dominant attractor of the system, indicating that the network architecture itself encodes a robust “resting” state.

Beyond the attractor, the model reproduces the precise temporal order of regulatory activity observed in living cells: G1 → S → G2 → M → G1. The transition from one phase to the next follows a strongly attractive trajectory in state space. At each step, a specific subset of nodes switches on, triggering the next subset in a cascade that proceeds without external cues. In this sense the network behaves as an auto‑excited system: internal feedback loops generate the necessary pulses to drive the cycle, while external damping mechanisms prevent runaway activation.

A comparative analysis with a Boolean model of the budding yeast Saccharomyces cerevisiae reveals striking differences. Although both networks share many homologous components, the S. pombe circuit exhibits self‑sustained oscillatory dynamics, whereas the S. cerevisiae circuit is strongly damped and requires periodic external excitation (e.g., nutrient signals) to maintain cycling. This contrast suggests divergent evolutionary strategies: S. pombe relies on an internally driven clock, while S. cerevisiae integrates environmental information more tightly into its cell‑cycle control.

The authors also test the robustness of the S. pombe model by imposing virtual mutations: permanently silencing individual nodes or altering logical rules. In most cases the system still reaches the G1 attractor or follows a modified but still cyclic trajectory. This resilience indicates that the network possesses multiple redundant pathways, allowing it to tolerate perturbations without catastrophic failure—a hallmark of biological robustness.

Overall, the study demonstrates that a parameter‑free Boolean representation can capture the essential dynamical behavior of a complex biological process. By showing that the circuit topology alone dictates both the stability of the G1 state and the directed progression through the cell‑cycle phases, the work provides a powerful framework for analyzing large regulatory networks where kinetic data are scarce. Potential applications include synthetic biology design (where logical circuits must be predictable), disease modeling (e.g., cancers with cell‑cycle dysregulation), and comparative evolutionary studies of cell‑cycle control across species.