Why do detailed circadian models (like the Leloup-Goldbeter model) include explicit phosphorylation and nuclear transport steps, rather than simply using a time delay?
Think about your answer, then reveal below.
Model answer: Detailed models include these biochemical steps because they serve multiple functions beyond simply introducing delay. Phosphorylation by kinases like CK1 (casein kinase 1) determines protein stability — mutations in CK1 binding sites cause familial advanced sleep phase syndrome by accelerating PER degradation and shortening the period. Nuclear transport gates when repressor complexes reach their transcriptional targets. Dimerization (PER-TIM, PER-CRY) creates stoichiometric relationships that affect amplitude and robustness. Each step adds nonlinearity that collectively generates the ultrasensitivity needed for oscillation without requiring an unrealistically steep Hill coefficient at any single step. Delay-based models capture the period correctly but cannot predict the effects of specific mutations, drug perturbations, or the amplitude and waveform of the oscillation.
The distinction between delay models and mechanistic ODE models parallels the broader systems biology tradeoff between parsimony and predictive detail. DDEs are mathematically elegant and analytically tractable, but the Leloup-Goldbeter approach — modeling each phosphorylation state, complex, and compartment — enables direct connection to experimental perturbations (kinase inhibitors, nuclear export block, mutation of specific phosphorylation sites), making it the standard for circadian modeling in practice.