A mathematical model of HPA axis dynamics and impacts of alcohol consumption

Gergley, Matthew; Chellamuthu, Vinodh

doi:10.2478/ijmce-2026-0006

Figures & Tables

The dotted curves represent the solutions to u(c) in the full model, Eq.(21). The solid curve represents the solutions to u(c) in the simplified model, Eq.(23). Parameters from Table 2.

The grey region represents bistability in the (y1*,y2*)(y_1^*,y_2^*) plane for κ=13×104\kappa = {1 \over 3} \times {10^4}. The three solid blue lines are defined as y2*=3y1*,y2*=y1*y_2^* = 3y_1^*,y_2^* = y_1^*, and y2*=13y1*y_2^* = {1 \over 3}y_1^*. The red dotted and dashed curves are defined by Eq.(35) and Eq.(36). Parameters used are from Table 2, and the red dot indicates the intersection of Eq.(35) and Eq.(36) when h = 1 for reference.

Dotted red - Fast c nullcline, full model; Solid red - Fast y nullcline, simplified model; Dashed blue - Slow cs nullcline, full model; Solid blue - Slow x nullcline, simplified model.

h2(t, 33, μ): These plots represent the variation in the rate h2, using g2, when fixing an age to x = 33 and applying different drinking schedules and alcohol types/concentrations (See Appendix 7). Recall t = 0 corresponds to 8:00 AM. (g2 and x = 33 were chosen for simulation purposes. The plots utilizing g1, g3 can be found in the GitHub within Appendix 7).

Plots for consuming ten drinks following the everyday schedule for beer, wine, and liquor. The time of drinking is 9:00 PM, 13 hours after peak cortisol levels (8:00 AM). The top subplot for each plot shows the linear decay of the BAC level for each occurrence of drinking. In the bottom plots, the solid purple line represents the individual’s circadian rhythm while the black dotted line represents the circadian rhythm if no alcohol was consumed. The vertical red dotted lines (t˜=24,48,⋯,168 hours)(\tilde t = 24,48, \cdots,168{\rm{ hours}}) represent a 24 hour block of time (8:00AM - 8:00AM (next day)).

Plots for consuming a random number of drinks each day (Nj = randint{1, 10}), following the everyday schedule for beer, wine, and liquor and varying W as needed. The value of Nj for each day can be found in Appendix 7.

Plots for consuming a random number of drinks each day, following the everyday schedule for beer, wine, and liquor. The value of Nj for each day can be found in Appendix 7.

Dimensional parameter values_

Parameter	Value	Description	Source(s)
c¯∞{{\bar c}_\infty }	0.2	Minimal stored baseline CRH	[28–30]
b	0.6	Stored CRH decay rate as a function of cortisol	[28–30]
t_c	69.3	CRH biosynthesis timescale	[28–30]
q₀	28.0	Maximum release rate of CRH in basal state	[28–30]
I₀	1.0	Basal level of the external stimuli	[28–30]
k	2.83	Relates stored CRH to CRH release rate	[28–30]
gc, max	42.0	Maximum auto/paracrine effect of CRH in the pituitary	[28–30]
n	5	Hill coefficient describing the self-up-regulation of CRH	[28–30]
q1−1q_1^{ - 1}	25.0	Circulating CRH conc. at half-maximum self-up-regulation	[28–30]
q₂	1.8	Ratio of CRH and cortisol decay rates	[28–30]
p2−1p_2^{ - 1}	0.067	or-complex conc. for half-maximum negative feedback	[28–30]
P₃	7.2	Ratio of ACTH and cortisol decay rates	[28–30]
p₄	0.05	(or-complex conc.)² at half-maximum positive feedback on r production	[28–30]
P₅	0.11	Basal GR production rate by pituitary	[28–30]
P₆	2.9	Ratio of GR and cortisol decay rates	[28–30]
ω	0.045	Frequency of 24hr circadian rhythm	[28–30]
z_B	0.6	Fluid ounces of alcohol for 12oz. beer	[25]
z_W	0.75	Fluid ounces of alcohol for 6.25oz. wine	[25]
z_L	0.8	Fluid ounces of alcohol for 2oz. liquor	[25]
β	0.00015	The alcohol elimination rate in kg/L/hr	[24, 31]
σ	0.6	Volume of distribution (a constant relating the distribution of water in the body in L/kg)	[24]
δ	0.8	The density of ethanol (0.8 oz. per fluid ounce)	[24]
W_f	2732.8	Average weight of a female in oz.	[32]
W_m	3196.8	Average weight of a male in oz.	[32]
A	30	Number of days in the period of concern
BAC_max	0.0008	Legal limit of BAC in the United States in kg/L	[33]
λ	(varies) with λ ∈ [0,7]	Number of days per week (7 days) where alcohol was consumed of any quantity and type.

Ranges of φ¯(t˜)\bar \varphi (\tilde t) by setting t˜=ti=0\tilde t = {t_i} = 0 in Eq_(5) with Nj = randint{ 1,10} (See Appendix 7) for each dj and changing z and W as needed_ In the everyday schedule, the value of φ¯(t˜)\bar \varphi (\tilde t) varies daily, since Nj is not fixed_ See Figure 6_

	Beer	Wine	Liquor
Male	0.3128≤φ¯(t˜)≤2.5025 (hours)0.3128 \le \bar \varphi (\tilde t) \le 2.5025{\rm{ (hours)}}	0.3910≤φ¯(t˜)≤3.5191 (hours)0.3910 \le \bar \varphi (\tilde t) \le 3.5191{\rm{ (hours)}}	0.4171≤φ¯(t˜)≤3.7538 (hours)0.4171 \le \bar \varphi (\tilde t) \le 3.7538{\rm{ (hours)}}
Female	0.3659≤φ¯(t˜)≤2.9274 (hours)0.3659 \le \bar \varphi (\tilde t) \le 2.9274{\rm{ (hours)}}	0.9148≤φ¯(t˜)≤4.1167 (hours)0.9148 \le \bar \varphi (\tilde t) \le 4.1167{\rm{ (hours)}}	0.4879≤φ¯(t˜)≤3.4153 (hours)0.4879 \le \bar \varphi (\tilde t) \le 3.4153{\rm{ (hours)}}

Maximum values of φ¯(t˜)\bar \varphi \left( {\tilde t} \right) which occur when t˜=ti=0\tilde t = {t_i} = 0 in Eq_(5) with Nj = 10 for each dj and changing z and W as needed_ In the everyday schedule, the value of φ¯(t˜)\bar \varphi \left( {\tilde t} \right) is uniform across each day since Nj is fixed_ This is not the case in Figure 6, however_ See Table 4_

	Beer	Wine	Liquor
Male	max(φ¯(t˜))=3.182 hours\max (\bar \varphi (\tilde t)) = 3.182\;{\rm{hours}}	max(φ¯(t˜))=3.910 hours\max (\bar \varphi (\tilde t)) = 3.182\;{\rm{hours}}	max(φ¯(t˜))=4.171 hours\max (\bar \varphi (\tilde t)) = 4.171\;{\rm{hours}}
Female	max(φ¯(t˜))=3.659 hours\max (\bar \varphi (\tilde t)) = 3.659\;{\rm{hours}}	max(φ¯(t˜))=4.574 hours\max (\bar \varphi (\tilde t)) = 4.574\;{\rm{hours}}	max(φ¯(t˜))=4.879 hours\max (\bar \varphi (\tilde t)) = 4.879\;{\rm{hours}}

j_ijmce-2026-0006_tab_000

College	Everyday	Weekdays	Weekend	Sporadic
0 \| 0 \| 0 \| 1 \| 1 \| 1 \| 1	1 \| 1 \| 1 \| 1 \| 1 \| 1 \| 1	1 \| 1 \| 1 \| 1 \| 1 \| 0 \| 0	0 \| 0 \| 0 \| 0 \| 0 \| 1 \| 1	α \| α \| α \| α \| α \| α \| α

Correspondence between t values and times of day_

t (hours past 8:00 AM)	Time of Day
0	8:00 AM
π2{\pi \over 2}	2:00 PM
π	8:00 PM
3π2 {{3\pi } \over 2}\,	2:00 AM (next day)
2π	8:00 AM (next day)

Maximum values of φ¯(t˜)\bar \varphi (\tilde t) which occur when t˜=ti=0\tilde t = {t_i} = 0 in Eq_(5) with Nj = 10 for each dj and changing z and W as needed_ In the everyday schedule, the value of φ¯(t˜)\bar \varphi (\tilde t) uniform across each day_ This is not the case in Fig_ 8, however_ See Table 6_

	Beer	Wine	Liquor
Male	max(φ¯(t˜))=3.182 hours\max (\bar \varphi (\tilde t)) = 3.182\;{\rm{hours}}	max(φ¯(t˜))=3.910 hours\max (\bar \varphi (\tilde t)) = 3.910\;{\rm{hours}}	max(φ¯(t˜))=4.171 hours\max (\bar \varphi (\tilde t)) = 4.171\;{\rm{hours}}
Female	max(φ¯(t˜))=3.659 hours\max (\bar \varphi (\tilde t)) = 3.659\;{\rm{hours}}	max(φ¯(t˜))=4.574 hours\max (\bar \varphi (\tilde t)) = 4.574\;{\rm{hours}}	max(φ¯(t˜))=4.879 hours\max (\bar \varphi (\tilde t)) = 4.879\;{\rm{hours}}

Ranges of φ¯(t˜)\bar \varphi (\tilde t) by setting t˜=ti=0\tilde t = {t_i} = 0 in Eq_(5) with Nj = randint{1, 10} (See Appendix 7) for each dj and changing z and W as needed_ In the everyday schedule, the value of φ¯(t˜)\bar \varphi (\tilde t) varies daily_ See Figure 8_

	Beer	Wine	Liquor
Male	0.3128≤φ¯(t˜)≤2.8153 (hours)0.3128 \le \bar \varphi (\tilde t) \le 2.8153{\rm{ (hours)}}	0.3910≤φ¯(t˜)≤3.5191 (hours)0.3910 \le \bar \varphi (\tilde t) \le 3.5191{\rm{ (hours)}}	0.8342≤φ¯(t˜)≤2.9196 (hours)0.8342 \le \bar \varphi (\tilde t) \le 2.9196{\rm{ (hours)}}
Female	0.3659≤φ¯(t˜)≤3.2933 (hours)0.3659 \le \bar \varphi (\tilde t) \le 3.2933{\rm{ (hours)}}	1.3722≤φ¯(t˜)≤4.1167 (hours)1.3722 \le \bar \varphi (\tilde t) \le 4.1167{\rm{ (hours)}}	0.9758≤φ¯(t˜)≤2.9274 (hours)0.9758 \le \bar \varphi (\tilde t) \le 2.9274{\rm{ (hours)}}

A mathematical model of HPA axis dynamics and impacts of alcohol consumption

Figures & Tables

Fig. 1

Fig. 2

Fig. 3

Fig. 4

Fig. 5

Fig. 6

Fig. 7

Fig. 8

Dimensional parameter values_

j_ijmce-2026-0006_tab_000

Correspondence between t values and times of day_

Ranges of φ¯(t˜)\bar \varphi (\tilde t) by setting t˜=ti=0\tilde t = {t_i} = 0 in Eq_(5) with Nj = randint{1, 10} (See Appendix 7) for each dj and changing z and W as needed_ In the everyday schedule, the value of φ¯(t˜)\bar \varphi (\tilde t) varies daily_ See Figure 8_

Paradigm

My account