CRH
—
pg/mL
ACTH
—
pg/mL
Cortisol
—
µg/dL
Dopamine
—
ng/mL
Norepineph.
—
ng/mL
Oxytocin
—
pg/mL
Hormone dynamics over time
ODE system — bidirectional PK coupling
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