Why the Lysogenic State of Phage λ Is So Stable: A Mathematical Modeling Approach

We develop a mathematical model of the phage λ lysis/lysogeny switch, taking into account recent experimental evidence demonstrating enhanced cooperativity between the left and right operator regions. Model parameters are estimated from available experimental data. The model is shown to have a single stable steady state for these estimated parameter values, and this steady state corresponds to the lysogenic state. When the CI degradation rate (γ_cI) is slightly increased from its normal value (γ_cI 0.0 min⁻¹), two additional steady states appear (through a saddle-node bifurcation) in addition to the lysogenic state. One of these new steady states is stable and corresponds to the lytic state. The other steady state is an (unstable) saddle node. The coexistence these two globally stable steady states (the lytic and lysogenic states) is maintained with further increases of γ_cI until γ_cI 0.35 min⁻¹, when the lysogenic steady state and the saddle node collide and vanish (through a reverse saddle node bifurcation) leaving only the lytic state surviving. These results allow us to understand the high degree of stability of the lysogenic state because, normally, it is the only steady state. Further implications of these results for the stability of the phage λ switch are discussed, as well as possible experimental tests of the model.