(545d) Modeling the Dynamics of Neuroendocrine-Immune Interactions in Collagen-Induced Arthritis

Modeling the Dynamics of Neuroendocrine-Immune Interactions in Collagen-induced Arthritis

R. Rao¹, D. DuBois^{2, 3}, R. Almon^{2, 3}, W.J. Jusko² and I.P. Androulakis^{1, 4}

¹Department of Chemical and Biochemical Engineering, Rutgers, the State University of New Jersey, Piscataway, NJ, 08854, USA

²Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, State University of New York at Buffalo, SUNY Buffalo, NY 14260, USA

³Department of Biological Sciences, State University of New York at Buffalo, SUNY Buffalo, NY 14260, USA

⁴Department of Biomedical Engineering, Rutgers, the State University of New Jersey, Piscataway, NJ, 08854, USA

Rheumatoid arthritis (RA) is a chronic inflammatory disease that affects about 1% of the population in the developed world, and is primarily characterized by synovial inflammation and cartilage deterioration, typically in the joints of the hands and feet (1). The hypothalamic-pituitary-adrenal (HPA) axis, a critical component of the neuroendocrine system, is thought to be an important regulator of disease activity. It has been hypothesized that adrenal insufficiency; the inability of the HPA axis to mount an adequate anti-inflammatory in response to ongoing chronic inflammation in RA might contribute to disease pathophysiology (2). In addition, circadian variation in clinical symptoms has been well documented in RA patients. Moreover, a disruption in the natural circadian rhythms of major mediators of the HPA axis, such as cortisol has been observed (3). A mechanistic understanding of the interplay between the circadian rhythms of the HPA axis and the immune system might provide useful insights into RA pathologies and optimization of treatment regimens.

Animal models of RA have been extensively used to study disease pathogenesis and in the discovery of anti-arthritic drugs (4). A rodent model of collagen induced arthritis (CIA), which involves immunization of the animal with an emulsion of the complete Freund's adjuvant and type II collagen, is considered to be the gold standard in vivo model for RA (5). In this work we discuss the development of a mathematical model that aims to evaluate the dynamic interactions between the HPA axis and pro-inflammatory cytokines, and their downstream effects on paw edema, a key disease end-point in CIA. In particular, we investigate the changes in circadian secretion of corticosterone (CST) and pro-inflammatory cytokines after induction of CIA, that ultimately result in circadian variability in paw edema. Moreover, we recapitulate the phenomenon of adrenal insufficiency that has been observed in both CIA as well as RA patients. These simulations will aid in the characterization of neuroendocrine-immune response in animal models of RA.

Our model builds on earlier work (6-9) and is comprised of three components; a central oscillatory compartment that simulates the circadian secretion of CST in the HPA axis, a peripheral compartment that captures the downstream effects of secreted CST on cytokine expression, and finally, a compartment that accounts for the disease endpoints, specifically, paw edema. The oscillations in the central compartment are generated due to the negative feedback between CST and corticotropin-releasing hormone (CRH)/adrenocorticotropic hormone (ACTH) using a modified Goodwin oscillator model that incorporates Michaelis-Menten kinetics to express the degradation terms in the hypothalamus, pituitary, and adrenals. The HPA axis is positively modulated by the expression of cytokines in the peripheral compartment, which in turn is balanced by the negative feedback of CST. This results in an interlinked positive and negative feedback circuit that enables the investigation of the dynamics of HPA axis-cytokine interactions upon the onset of inflammation. In this model, we assume that the severity of paw edema is in direct correlation to the expression of pro-inflammatory cytokines. Furthermore, we explicitly model the emergence of adrenal insufficiency in response to chronically elevated cytokines by postulating an indirect response mechanism that effectively results in stimulating the clearance of ACTH in the HPA axis. Disease progression is initiated by a delay in the onset of the increased pro-inflammatory cytokine expression by using a simple transduction model (6).

Model predictions demonstrate that elevated cytokine expression results in circadian variability in paw edema progression, reminiscent of the circadian rhythms of disease symptoms exhibited by RA patients. Notably, our model predicts that an increase in mean secretion of CST after the induction of the disease is accompanied by a decrease in the amplitude of the CST oscillation. This phenomenon persists even after the development of adrenal insufficiency and the consequent drop in mean CST secretion toward levels present prior to the onset of the disease. This suppression of circadian rhythm has been experimentally observed in a number of important neuroendocrine and HPA axis mediators during adjuvant induced arthritis, including; ACTH, CST and adrenocortical ornithine decarboxylase (10). Further analysis of possible alterations in the phase of oscillation of CST and the cytokines during disease progression is currently ongoing. 

In conclusion, our model captures key features of chronic inflammation in animal models of RA. In doing so, it highlights the importance of accounting for circadian rhythms while studying the neuroendocrine-immune system interactions in RA.

  References

1.         McInnes IB, Schett G. The Pathogenesis of Rheumatoid Arthritis. N Engl J Med. 2011 Dec 8;365(23):2205–19.

2.         Harbuz MS, Korendowych E, Jessop DS, Crown AL, Li pdfan S. L., Kirwan JR. Hypothalamo-pituitary-adrenal axis dysregulation in patients with rheumatoid arthritis after the dexamethasone/corticotrophin releasing factor test. J Endocrinol. 2003 Jul;178(1):55–60.

3.         Cutolo M, Masi AT. Circadian Rhythms and Arthritis. Rheum Dis Clin N Am. 2005 Feb;31(1):115–29.

4.         Kannan K, Ortmann RA, Kimpel D. Animal models of rheumatoid arthritis and their relevance to human disease. Pathophysiology. 2005 Oct;12(3):167–81.

5.         Brand DD, Latham KA, Rosloniec EF. Collagen-induced arthritis. Nat Protoc. 2007;2(5):1269–75.

6.         Earp JC, DuBois DC, Molano DS, Pyszczynski NA, Keller CE, Almon RR, et al. Modeling Corticosteroid Effects in a Rat Model of Rheumatoid Arthritis I: Mechanistic Disease Progression Model for the Time Course of Collagen-Induced Arthritis in Lewis Rats. J Pharmacol Exp Ther. 2008 Aug;326(2):532–45.

7.         Mavroudis PD, Corbett SA, Calvano SE, Androulakis IP. Mathematical modeling of light-mediated HPA axis activity and downstream implications on the entrainment of peripheral clock genes. Physiol Genomics. 2014 Oct 15;46(20):766–78.

8.         Sriram K, Rodriguez-Fernandez M, Doyle FJ. Modeling Cortisol Dynamics in the Neuro-endocrine Axis Distinguishes Normal, Depression, and Post-traumatic Stress Disorder (PTSD) in Humans. PLoS Comput Biol [Internet]. 2012 Feb [cited 2015 May 1];8(2). Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3280965/

9.         Ramakrishnan R, DuBois DC, Almon RR, Pyszczynski NA, Jusko WJ. Fifth-generation model for corticosteroid pharmacodynamics: application to steady-state receptor down-regulation and enzyme induction patterns during seven-day continuous infusion of methylprednisolone in rats. J Pharmacokinet Pharmacodyn. 2002 Feb;29(1):1–24.

10.       Cardinali DP, Esquifino AI. Circadian disorganization in experimental arthritis. Neurosignals. 2003 Dec;12(6):267–82.