Insights into the mechanisms of thymus involution and regeneration by modeling the glucocorticoid-induced perturbation of thymocyte populations dynamics

T-cells develop in the thymus and based on CD4 and CD8 expressions there are four main thymocyte populations in a normal mouse thymus. Currently, there are several mathematical models that describe the dynamics of thymocyte populations in a normal thymus, but only a few of them model the transient perturbation of their homeostasis. Our aim is to model the perturbation in the dynamics of each thymocyte population which is induced by the administration of a glucocorticoid, i.e. dexamethasone. The proposed approach relies on extending a four compartment thymus model based on differential equations by adding perturbation terms either globally (at the level of each equation) or locally (at the level of proliferation, death, and transfer rates). By fitting the perturbed model with experimental data on mice thymi collected before and after the administration of dexamethasone, it was possible to estimate the relevant parameters using a population-based stochastic search method. The fitted model is further used to conduct a quantitative analysis on the differentiated impact of dexamethasone on each T-cell population and on proliferation, death, and transfer processes. The obtained quantitative information on the perturbation could be used to explore and modify the flow of thymocytes between thymus compartments in order to elucidate the mechanisms of thymus involution and its subsequent regeneration. Since glucocorticoids are raised in many pathological situations, such a model could be useful in evaluating the impact of diseases on thymocyte dynamics in the thymus.