CRH
—
pg/mL
ACTH
—
pg/mL
Cortisol
—
µg/dL
Dopamine
—
ng/mL
Norepineph.
—
ng/mL
Oxytocin
—
pg/mL
Hormone dynamics over time
⚡ Free 120-second test drive
Run the live HPA-axis engine for 120 seconds — try the scenarios, adjust doses, and watch the feedback loops respond in real time. No signup.
Test drive ended
Subscribe to keep simulating — full API access, longer time horizons, and data export.
See plans →
ODE system — bidirectional PK coupling
click any equation to tune its constants ↓
▸ HPA axis — CRH
dCRH/dt = K_syn + stress·2 + K_NE·NE
− K_deg·CRH − K_fb·f(CORT)·CRH
− K_deg·CRH − K_fb·f(CORT)·CRH
▸ ACTH
dACTH/dt = K_syn·CRH − K_deg·ACTH
▸ Cortisol — oxytocin buffered
dCORT/dt = K_syn·ACTH
− (K_deg + K_OT·OT)·CORT
− (K_deg + K_OT·OT)·CORT
▸ Dopamine — Hill suppression
dDA/dt = K_syn·(1 − 0.6·CORT/(CORT+20))
+ K_OT·OT − K_deg·DA
+ K_OT·OT − K_deg·DA
▸ Norepinephrine
dNE/dt = K_syn + stress·1.5 − K_deg·NE
▸ Oxytocin — social coupling
dOT/dt = K_syn + social·0.5 + K_DA·DA
− (K_deg + K_CORT·CORT)·OT
− (K_deg + K_CORT·CORT)·OT
🧩 Model Builder — extend the system
Add your own signals and couplings on top of the validated base model, then simulate. Exports to SBML.