Circadian Rhythm: Clock Gene Networks — Modelling the Cellular Oscillator

The core transcription-translation feedback loop

The mammalian molecular clock is built from a transcription-translation feedback loop (TTFL). The transcription factors CLOCK and BMAL1 heterodimerise and drive expression of the Per (Period) and Cry (Cryptochrome) genes. PER and CRY proteins accumulate in the cytoplasm, dimerise, translocate back into the nucleus, and directly inhibit the CLOCK:BMAL1 complex that made them — a negative-feedback loop built entirely from gene expression rather than fast enzyme kinetics.

CLOCK:BMAL1 → activates Per, Cry transcription
PER, CRY (protein) accumulate → dimerise → re-enter nucleus
PER:CRY → represses CLOCK:BMAL1 activity → transcription falls
PER, CRY degrade (via CK1δ/ε-mediated phosphorylation → proteasome)
repression lifts → cycle restarts     — full cycle period ≈ 24 h in mammals

A second, interlocked loop (ROR/REV-ERB acting on Bmal1 itself) adds robustness, but the PER/CRY-CLOCK/BMAL1 loop is the minimal core that a mathematical model needs to reproduce a 24-hour oscillation.

The Goodwin oscillator: a minimal feedback model

Brian Goodwin's 1965 model reduces any such loop to three variables: mRNA (X), an intermediate protein (Y), and the final repressor protein (Z), with Z inhibiting X's own transcription:

dX/dt = v₁ / (1 + (Z/K)^n) − d₁X    // mRNA (e.g. Per/Cry)
dY/dt = v₂X − d₂Y                    // cytoplasmic protein
dZ/dt = v₃Y − d₃Z                    // nuclear repressor (PER:CRY complex)

The repression term is a Hill function with cooperativity exponent n. Goodwin showed that sustained oscillation in this three-variable ODE system requires n > 8 — an unrealistically high degree of molecular cooperativity, since real PER/CRY dimerisation and repression involve far gentler nonlinearity.

Why delay matters: DDEs vs. ODEs

The real molecular clock achieves robust 24-hour oscillation with much lower cooperativity than Goodwin's ODE requires, because the biology introduces an explicit time delay: transcription, nuclear export, translation, protein folding, dimerisation and re-import of PER:CRY together take several hours — not instantaneous, as an ODE implicitly assumes. Replacing the instantaneous feedback term with a delay differential equation (DDE) captures this directly:

dX/dt = v₁ / (1 + (Z(t − τ)/K)^n) − d₁X
// repression acts on X using Z's value τ hours in the past

With a delay τ of a few hours, sustained oscillations appear even at n = 1 (simple, non-cooperative Michaelis-Menten repression) — matching the modest cooperativity actually observed for CLOCK:BMAL1 repression by PER:CRY. Delay, not ultra-sensitivity, is the feature that makes the biological clock robust: explicit transcriptional and translational delay substitutes for cooperativity in generating a limit cycle.

Limit cycles and the Hopf bifurcation

As the delay τ (or, equivalently, the Hill coefficient n) increases past a critical value, the system's single stable fixed point loses stability in a Hopf bifurcation, and a stable limit cycle is born — a closed, self-sustaining orbit in (X, Y, Z) phase space that the system settles into regardless of small perturbations to its initial state. Below threshold the fixed point is stable and mRNA/protein levels settle to constants with no rhythm; at the bifurcation the fixed point loses stability; beyond it a self-sustained ~24 h oscillation emerges whose amplitude is set by the nonlinearity, not by initial conditions; and daily light input then phase-shifts that limit cycle to align with the 24 h solar day (entrainment).

This independence from initial conditions is the mathematical signature of a robust biological clock: a cell coming from any starting state of PER/CRY levels converges onto the same rhythmic orbit, which is why isolated fibroblasts all settle into comparable ~24 h cycles despite starting from different transcriptional states.

Coupled cells: from single oscillator to tissue rhythm

A single cell's clock, on its own, is a noisy oscillator — its period drifts by roughly ±1–2 hours from cell to cell due to gene expression noise. In the suprachiasmatic nucleus, roughly 20,000 such noisy single-cell oscillators are coupled through neuropeptide signalling (chiefly VIP acting on VPAC2 receptors), which synchronises their phases into a single, precise tissue-level rhythm — the same coupled-oscillator synchronisation mathematics used for coupled pendulums or the Kuramoto model. Weak coupling between single-cell oscillators still allows the population average to keep time even if individual cells are noisy, which is why the whole-body clock (light entrainment, jet lag, melatonin output) is so much more temporally precise than any one of its constituent cells.

Frequently asked questions

Try it live

Everything above runs in your browser — open Circadian Rhythm and watch the delay differential equation settle into a stable ~24-hour-scaled limit cycle. Nothing is installed, nothing is uploaded.