A mathematical model of the glucocorticoid effects following preventive administration in a model of neuroinflammation

Abstract

A mathematical model of neuroinflammation is proposed that describes the interplay among the peripheral inflammatory response, microglial activation, neural tissue damage, changes in blood–brain barrier permeability, and the effects of preventive glucocorticoid administration. The model is formulated as a system of ordinary differential equations and is intended to investigate the conditions governing the transition from an acute inflammatory response to a stable chronic state. The analysis includes verification of the models’ well-posedness, identification of equilibria, local stability analysis, bifurcation analysis, and parameter sensitivity analysis. Three dynamic regimes were identified: rapid resolution of inflammation, delayed recovery accompanied by persistent changes in the central nervous system, and transition to chronic neuroinflammation. The development of chronic neuroinflammation was shown to be associated with strengthening of the positive feedback circuit between central inflammatory mediators and activated microglia. The parameters exerting the greatest influence on the duration of the inflammatory response, accumulation of neural tissue damage, and decline in the effectiveness of glucocorticoid-mediated regulation were identified. The model may be used to interpret experimental data, assess the contributions of individual components of the inflammatory response, and investigate the conditions underlying the development of chronic neuroinflammation.

Keywords

mathematical modeling / ordinary differential equations / neuroinflammation / glucocorticoids / bifurcation analysis

References

[1] Müller L, Di Benedetto S (2025) Neuroimmune crosstalk in chronic neuroinflammation: microglial interactions and immune modulation. Front Cell Neurosci 19: 1575022. https://doi.org/10.3389/fncel.2025.1575022

[2] Chen Y, Yin P, Chen Q, et al. (2026) Neurodegenerative diseases and immune system: from pathogenic mechanism to therapy. Neural Regen Res 21: 3387-3410. https://doi.org/10.4103/NRR.NRR-D-25-00274

[3] Reichardt SD, Amouret A, Muzzi C, et al. (2021) The role of glucocorticoids in inflammatory diseases. Cells 10: 2921. https://doi.org/10.3390/cells10112921

[4] Stahn C, Buttgereit F (2008) Genomic and nongenomic effects of glucocorticoids. Nat Clin Pract Rheumatol 4: 525-533. https://doi.org/10.1038/ncprheum0898

[5] Samuel S, Nguyen T, Choi HA (2017) Pharmacologic characteristics of corticosteroids. J Neurocrit Care 10: 53-59. https://doi.org/10.18700/jnc.170035

[6] Valotis A, Högger P (2007) Human receptor kinetics and lung tissue retention of the enhanced-affinity glucocorticoid fluticasone furoate. Resp Res 8: 54. https://doi.org/10.1186/1465-9921-8-54

[7] Hochhaus G (2008) Relative receptor affinity comparisons among inhaled/intranasal corticosteroids: perspectives on clinical relevance. Resp Res 9: 75. https://doi.org/10.1186/1465-9921-9-75

[8] Daley-Yates PT (2015) Inhaled corticosteroids: potency, dose equivalence and therapeutic index. Br J Clin Pharmacol 80: 372-380. https://doi.org/10.1111/bcp.12637

[9] Buttgereit F, Brand MD, Burmester GR (1999) Equivalent doses and relative drug potencies for non-genomic glucocorticoid effects: a novel glucocorticoid hierarchy. Biochem Pharmacol 58: 363-368. https://doi.org/10.1016/S0006-2952(99)00090-8

[10] Meijer OC, de Lange ECM, Breimer DD, et al. (1998) Penetration of dexamethasone into brain glucocorticoid targets is enhanced in mdr1A P-glycoprotein knockout mice. Endocrinology 139: 1789-1793. https://doi.org/10.1210/endo.139.4.5917

[11] Karssen AM, Meijer OC, van der Sandt ICJ, et al. (2001) Multidrug resistance P-glycoprotein hampers the access of cortisol but not of corticosterone to mouse and human brain. Endocrinology 142: 2686-2694. https://doi.org/10.1210/endo.142.6.8213

[12] Mason BL, Pariante CM, Thomas SA (2008) A revised role for P-glycoprotein in the brain distribution of dexamethasone, cortisol, and corticosterone in wild-type and ABCB1A/B-deficient mice. Endocrinology 149: 5244-5253. https://doi.org/10.1210/en.2008-0041

[13] Sorrells SF, Sapolsky RM (2007) An inflammatory review of glucocorticoid actions in the CNS. Brain Behav Immun 21: 259-272. https://doi.org/10.1016/j.bbi.2006.11.006

[14] Frank MG, Miguel ZD, Watkins LR, et al. (2010) Prior exposure to glucocorticoids sensitizes the neuroinflammatory and peripheral inflammatory responses to E. coli lipopolysaccharide. Brain Behav Immun 24: 19-30. https://doi.org/10.1016/j.bbi.2009.07.008

[15] Frank MG, Watkins LR, Maier SF (2015) The permissive role of glucocorticoids in neuroinflammatory priming: mechanisms and insights. Curr Opin Endocrinol Diabetes Obes 22: 300-305. https://doi.org/10.1097/MED.0000000000000168

[16] Foster-Powell A, Rostami-Hodjegan A, Meno-Tetang G, et al. (2025) Mathematical modeling of neuroinflammation in neurodegenerative diseases. CPT Pharmacometrics Syst Pharmacol 14: 1908-1922. https://doi.org/10.1002/psp4.70064

[17] Han J, Zhang Z, Zhang P, et al. (2025) The roles of microglia and astrocytes in neuroinflammation of Alzheimer’s disease. Front Neurosci 19: 1575453. https://doi.org/10.3389/fnins.2025.1575453

[18] Vaughan LE, Ranganathan PR, Kumar RG, et al. (2018) A mathematical model of neuroinflammation in severe clinical traumatic brain injury. J Neuroinflammation 15: 345. https://doi.org/10.1186/s12974-018-1384-1

[19] Batulin D, Lagzi F, Vezzani A, et al. (2022) A mathematical model of neuroimmune interactions in epileptogenesis for discovering treatment strategies. iScience 25: 104343. https://doi.org/10.1016/j.isci.2022.104343

[20] Torres N, Molina E, Pujo-Menjouet L (2025) An optimal control problem for anti-inflammatory treatments of Alzheimer’s disease. J Math Biol 91: 1. https://doi.org/10.1007/s00285-025-02227-8

[21] Eftimie R, Gillard JJ, Cantrell DA (2016) Mathematical models for immunology: current state of the art and future research directions. Bull Math Biol 78: 2091-2134. https://doi.org/10.1007/s11538-016-0214-9

[22] Chamberland É, Moravveji S, Doyon N, et al. (2024) A computational model of Alzheimer’s disease at the nano, micro, and macroscales. Front Neuroinform 18: 1348113. https://doi.org/10.3389/fninf.2024.1348113

[23] Lebedeva EYa, Kukushkina ES, Tyutunnik TV, et al. (2025) Glucocorticosteroid-induced behavioral changes during endotoxemia in rats. Med Acad J 25: 82-89. https://doi.org/10.17816/MAJ631323

[24] Heneka MT, Carson MJ, Khoury JE, et al. (2015) Neuroinflammation in Alzheimer’s disease. Lancet Neurol 14: 388-405. https://doi.org/10.1016/S1474-4422(15)70016-5

[25] Hindmarsh AC (1983) ODEPACK, a systematized collection of ODE solvers, In: Stepleman, R.S., Scientific Computing, Amsterdam: North-Holland, 55-64.

[26] Petzold L (1983) Automatic selection of methods for solving stiff and nonstiff systems of ordinary differential equations. SIAM J Sci Stat Comput 4: 136-148. https://doi.org/10.1137/0904010

[27] Virtanen P, Gommers R, Oliphant TE, et al. (2020) SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat Methods 17: 261-272. https://doi.org/10.1038/s41592-019-0686-2

[28] Sheiner LB, Stanski DR, Vozeh S, et al. (1979) Simultaneous modeling of pharmacokinetics and pharmacodynamics: application to d-tubocurarine. Clin Pharmacol Ther 25: 358-371. https://doi.org/10.1002/cpt1979253358

[29] Holford NHG, Sheiner LB (1982) Kinetics of pharmacologic response. Pharmacol Ther 16: 143-166. https://doi.org/10.1016/0163-7258(82)90051-1

[30] Brent RP (1973) Algorithms for Minimization without Derivatives, Englewood Cliffs: Prentice-Hall.

[31] Kuznetsov YA (2004) Elements of Applied Bifurcation Theory, 3 Eds., New York: Springer. https://doi.org/10.1007/978-1-4757-3978-7

[32] Jolliffe IT (2002) Principal Component Analysis, 2 Eds., New York: Springer. https://doi.org/10.1007/b98835

[33] Morris MD (1991) Factorial sampling plans for preliminary computational experiments. Technometrics 33: 161-174. https://doi.org/10.1080/00401706.1991.10484804

[34] Marino S, Hogue IB, Ray CJ, et al. (2008) A methodology for performing global uncertainty and sensitivity analysis in systems biology. J Theor Biol 254: 178-196. https://doi.org/10.1016/j.jtbi.2008.04.011

[35] Saltelli A, Annoni P, Azzini I, et al. (2010) Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index. Comput Phys Commun 181: 259-270. https://doi.org/10.1016/j.cpc.2009.09.018

[36] Hartman P (2002) Ordinary Differential Equations, 2 Eds., Philadelphia: SIAM. https://doi.org/10.1137/1.9780898719222

[37] Smith HL (1995) Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperative Systems, Providence: American Mathematical Society. https://doi.org/10.1090/surv/041

[38] Hirsch MW, Smith HL (2006) Monotone dynamical systems, In: Cañada A, Drábek P, Fonda A, Author, Handbook of Differential Equations: Ordinary Differential Equations, Amsterdam: Elsevier, 239-357. https://doi.org/10.1016/S1874-5725(05)80006-9

[39] Khalil HK (2002) Nonlinear Systems, 3 Eds., Upper Saddle River: Prentice Hall.

[40] Hadamard J (1902) Sur les problèmes aux dérivées partielles et leur signification physique. Princeton Univ Bull 13: 49-52.