Mathematical modeling reveals how the speed of endocrine regulation should affect baseline and stress-induced glucocorticoid levels

Unpredictable environmental changes displace individuals from homeostasis and elicit a stress response. In vertebrates, the stress response is mediated mainly by glucocorticoids (GCs) which initiate physiological changes while minimizing allostatic overload. Individuals and species vary consistently in baseline and stress-induced GC levels and the speed with which GC levels can be upregulated or downregulated, but the extent to which variation in hormone regulation influences baseline and stress-induced GC levels is unclear. Using mathematical modeling, we tested how GC regulation rate, frequencies and durations of acute stressors, fitness functions, and allostatic overload affect GC levels during control and acute stress periods. As GC regulation rate slows, baseline and acute stress-induced GC levels become more similar. When the speed of up- and downregulation decreased, hormone levels became more linked to anticipated future conditions to avoid fitness costs of mismatching a new environmental state. More frequent acute stressors caused baseline and acute stress-induced GC levels to converge. When fitness was more tightly linked to hormone levels during acute stress periods than during control states, the speed of upregulation influenced optimal hormone levels more than the downregulation rate. With allostatic overload costs included, predicted GC levels were lower and more dependent on the frequency of past acute stressors. Our results show the value of optimality modeling to study the hormonal response to stressors and suggest GC levels depend on past and anticipated future environmental states as well as individual differences in hormone regulation. •Hormone regulation rates affect optimal baseline and stress-induced hormone levels.•Frequent stressors cause baseline and stress-induced hormone levels to converge.•Stronger selection on stress-induced hormone levels impacts upregulation rates most.•Allostatic overload costs cause optimal hormone levels to depend on past stressors.•Optimality modeling identifies assumptions and exposes their effects.