×

The minimal model of the hypothalamic-pituitary-adrenal axis. (English) Zbl 1311.92075

Summary: This paper concerns ODE modeling of the hypothalamic-pituitary- adrenal axis (HPA axis) using an analytical and numerical approach, combined with biological knowledge regarding physiological mechanisms and parameters. The three hormones, CRH, ACTH, and cortisol, which interact in the HPA axis are modeled as a system of three coupled, nonlinear differential equations. Experimental data shows the circadian as well as the ultradian rhythm. This paper focuses on the ultradian rhythm. The ultradian rhythm can mathematically be explained by oscillating solutions. Oscillating solutions to an ODE emerges from an unstable fixed point with complex eigenvalues with a positive real parts and a non-zero imaginary parts. The first part of the paper describes the general considerations to be obeyed for a mathematical model of the HPA axis. In this paper we only include the most widely accepted mechanisms that influence the dynamics of the HPA axis, i.e. a negative feedback from cortisol on CRH and ACTH. Therefore we term our model the minimal model. The minimal model, encompasses a wide class of different realizations, obeying only a few physiologically reasonable demands. The results include the existence of a trapping region guaranteeing that concentrations do not become negative or tend to infinity. Furthermore, this treatment guarantees the existence of a unique fixed point. A change in local stability of the fixed point, from stable to unstable, implies a Hopf bifurcation; thereby, oscillating solutions may emerge from the model. Sufficient criteria for local stability of the fixed point, and an easily applicable sufficient criteria guaranteeing global stability of the fixed point, is formulated. If the latter is fulfilled, ultradian rhythm is an impossible outcome of the minimal model and all realizations thereof. The second part of the paper concerns a specific realization of the minimal model in which feedback functions are built explicitly using receptor dynamics. Using physiologically reasonable parameter values, along with the results of the general case, it is demonstrated that un-physiological values of the parameters are needed in order to achieve local instability of the fixed point. Small changes in physiologically relevant parameters cause the system to be globally stable using the analytical criteria. All simulations show a globally stable fixed point, ruling out periodic solutions even when an investigation of the ‘worst case parameters’ is performed.

MSC:

92C40 Biochemistry, molecular biology
92C45 Kinetics in biochemical problems (pharmacokinetics, enzyme kinetics, etc.)
34C60 Qualitative investigation and simulation of ordinary differential equation models
34D23 Global stability of solutions to ordinary differential equations
Full Text: DOI

References:

[1] Allen LJS (2007) An introduction to mathematical biology. Pearson Education, Inc
[2] Andersen M, Vinther F (2010) Mathematical modeling of the hypothalamic-pituitary-adrenal axis, IMFUFA tekst 469. Roskilde University, Universitetsvej 1, 4000 Roskilde, Denmark. http://milne.ruc.dk/imfufatekster/pdf/469.pdf
[3] Bairagi N, Chatterjee S, Chattopadhyay J (2008) Variability in the secretion of corticotropin-releasing hormone adrenocortcotropic hormone and cortisol and understanding of the hypothalamic-pituitary-adrenal axis–a mathematical study based on clinical evidence. Math Med Biol 25: 37–63 · Zbl 1146.92008 · doi:10.1093/imammb/dqn003
[4] Ben-Zvi A, Vernon SD, Broderick G (2009) Model-based therapeutic correction of hypothalamic-pituitary-adrenal axis dysfunction. PLoS Comput Biol 5:e1000273
[5] Bingzhenga L, Zhenye Z, Liansong C (1990) A mathematical model of the regulation system of the secretion of glucocorticoids. J Biol Phys 17(4): 221–233 · doi:10.1007/BF00386598
[6] Carroll BJ, Cassidy F, Naftolowitz D, Tatham NE, Wilson WH, Iranmanesh A, Liu PY, Veldhuis JD (2007) Pathophysiology of hypercortisolism in depression. Acta Psychiatr Scand 115: 90–103 · doi:10.1111/j.1600-0447.2007.00967.x
[7] Christensen KC, Hannesson K, Hansen DF, Jensen LF, Nielsen KHM (2007) Matematisk modellering af hpa-aksen. Master’s thesis, Roskilde Universitet, Universitetsvej 1, 4000 Roskilde, Denmark
[8] Chrousos G (1998) Editorial: ultradian, circadian, and stress-related hypothalamic-pituitary-adrenal axis activity–a dynamic digital-to-analog modulation. Endocrinology 139: 437–440 · doi:10.1210/en.139.2.437
[9] Conrad M, Hubold C, Fischer B, Peters A (2009) Modeling the hypothalamus-pituitary-adrenal system: homeostasis by interacting positive and negative feedback. J Biol Phys 35: 149–162 · doi:10.1007/s10867-009-9134-3
[10] de Kloet ER, Vreugdenhil E, Oitzl MS, Joels M (1998) Brain corticosteroid receptor balance in health and disease. Endocr Rev 19(3): 269–301 · doi:10.1210/er.19.3.269
[11] Felig P, Frohman LA (2001) Endocrinology and metabolism, 4th edn. McGraw-Hill, New York
[12] Griffin JE, Ojeda SR (2004) Textbook of endocrine physiology, 5th edn. Oxford University Press, New York
[13] Gupta S, Aslakson E, Gurbaxani BM, Vernon SD (2007) Inclusion of the glucocorticoid receptor in a hypothalamic pituitary adrenal axis model reveals bistability. Theor Biol Med Model 4: 8 · doi:10.1186/1742-4682-4-8
[14] Hashimoto K, Nishioka T, Numata Y, Ogasa T, Kageyama J, Suemaru S (1993) Plasma levels of corticotropin-releasing hormone in hypothalamic–pituitary–adrenal disorders and chronic renal failure. Acta Endocrinol 128(6): 503–507
[15] Jelic S, Cupic Z, Kolar-Anic L (2005) Mathematical modeling of the hypothalamic-pituitary-adrenal system activity. Math Biosci 197: 173–187 · Zbl 1076.92018 · doi:10.1016/j.mbs.2005.06.006
[16] Keenan DM, Veldhuis JD (2003) Cortisol feedback state governs adrenocorticotropin secretory-burst shape, frequency, and mass in a dual-waveform construct: time of day-dependent regulation. Am J Physiol Regul Integr Comp Physiol 285: 950–961
[17] Keenan DM, Licinio J, Veldhuis JD (2001) A feedback-controlled ensemble model of the stress-responsive hypothalamo-pituitaryadrenal axis. Proc Natl Acad Sci USA 98(7): 4028–4033 · doi:10.1073/pnas.051624198
[18] Kyrylov V, Severyanova LA, Zhiliba A (2004) The ultradian pulsatility and nonlinear effects in the hypothalamic-pituitary-adrenal axis. In: The 2004 international conference on health sciences simulation (HSS 2004), San Diego, California
[19] Kyrylov V, Severyanova LA, Vieira A (2005) Modeling robust oscillatory behaviour of the hypothalamic-pituitary-adrenal axis. IEEE Trans Biomed Eng 52: 1977–1983 · doi:10.1109/TBME.2005.857671
[20] Landau L, Lifshitz E (2008) Course of theoretical physics, vol 5. Statistical physics 3rd edn, part 1. Elsevier, Amsterdam
[21] Liu YW, Hu ZH, Peng JH, Liu BZ (1999) A dynamical model for the pulsatile secretion of the hypothalamo-pituary-adrenal axis. Math Comput Model 29(4): 103–110 · doi:10.1016/S0895-7177(99)00043-6
[22] Murray J (2002) Mathematical biology. Part I: an introduction, 3rd edn. Springer, New York · Zbl 1006.92001
[23] Savic D, Jelic S (2005) A mathematical model of the hypothalamo-pituitary-adrenocortical system and its stability analysis. Chaos Solitons Fractals 26: 427–436 · Zbl 1065.92016 · doi:10.1016/j.chaos.2005.01.013
[24] Savic D, Jelic S (2006) Stability of a general delay differential model of the hypothalamo-pituitary-adrenocortical system. Int J Bifurcation Chaos 16: 3079–3085 · Zbl 1136.92318 · doi:10.1142/S0218127406016665
[25] Tortora GJ, Derrickson B (2006) Principles of anatomy and physiology, 11th edn. Wiley, New York
[26] Veldhuis JD, Iranmanesh A, Lizarralde G, Johnson ML (1989) Amplitude modulation of a burstlike mode of cortisol secretion subserves the circadian glucocorticoid rhythm. Am J Physiol Endocrinol Metab 257: 6–14
[27] Wilson JD, Foster DW (1992) Williams textbook of endocrinology, 8th edn. W. B. Saunders Company, Philadelphia
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.