The Molecular Clockwork
At the heart of every cell in your body ticks a molecular clock — a self-sustaining oscillator built from interlocking feedback loops of gene transcription and protein degradation. The core loop centers on four gene families: CLOCK and BMAL1 (the activators) and PER and CRY (the repressors). CLOCK/BMAL1 heterodimers bind E-box promoter elements to drive PER and CRY transcription. The resulting proteins accumulate over hours, complex together, translocate to the nucleus, and directly inhibit CLOCK/BMAL1 — silencing their own genes. As PER/CRY are degraded by proteasomes, the brake lifts and the cycle begins anew.
The Goodwin Oscillator
Brian Goodwin's 1965 model distilled this feedback into its mathematical essence: a molecule that represses its own production through a nonlinear (Hill-type) function, with time delays introduced by intermediate steps. The critical requirement for sustained oscillation is sufficient nonlinearity — the Hill coefficient must exceed a threshold that depends on the number of intermediate steps. This simulation implements a refined Goodwin model calibrated to mammalian circadian parameters, showing how transcription rate, degradation rate, and cooperativity jointly determine period and amplitude.
Period-Determining Mechanisms
The approximately 24-hour period is not hardwired into any single rate constant — it emerges from the concert of many. However, certain steps are particularly influential. Casein kinase 1ε (CK1ε) phosphorylation of PER proteins marks them for ubiquitination and proteasomal degradation. The tau mutation in hamster CK1ε accelerates PER degradation, shortening the period to 20 hours. In humans, a PER2 mutation that removes the CK1ε phosphorylation site causes Familial Advanced Sleep Phase Syndrome — people who fall asleep at 7 PM and wake at 3 AM, running on a ~22-hour clock.
Robustness and Fragility
The circadian clock is remarkably robust to perturbations in most parameters — a property called structural stability. Temperature compensation keeps the period nearly constant from 25°C to 37°C (unlike most biochemical reactions that accelerate with heat). Yet the clock is exquisitely sensitive to light input, which can shift its phase by hours — the basis for entrainment to the solar cycle. This simulation reveals the parameter space where oscillations are robust versus the boundaries where rhythmicity collapses.