 Research
 Open Access
 Published:
A model of hematopoietic bone marrow apoptosis during growth factor deprivation in combination with a cytokine
Theoretical Biology and Medical Modelling volume 15, Article number: 8 (2018)
Abstract
Background
The process by which blood cells are formed is referred to as hematopoiesis. This process involves a complex sequence of phases that blood cells must complete. During hematopoiesis, a small fraction of cells undergo cell death. Causes of cell death are dependent upon various factors; one such factor being growth factor deprivation.
Methods
In this paper, a mathematical model of hematopoiesis during growth factor deprivation is presented. The model consists of a set of three coupled differential delay equations. Phase plane and linear stability analysis are performed in order to locate and determine stability of fixed points. Numerical simulations of the governing equations are run and provide a visual display of the behavior of the stem cell population undergoing growth factor deprivation. In addition, the effect of cytokine administration is incorporated in the model in an effort to understand how cytokine administration can offset the negative effects of apoptosis caused by growth factor deprivation.
Conclusions
The model produces qualitatively similar results to that observed during serum deprivation. The model captures apoptosis levels of cells at different time points. Additionally, it is shown that cytokine administration stabilizes the stem cell count.
Background
Hematopoiesis is the biological process by which hematopoietic stem cells (HSC) located in the bone marrow generate blood cells [1]. HSCs are responsible for the daily production of blood cells estimated to be approximately 10^{11}–10^{12} in a healthy individual [2]. Understanding the exact process that occurs during the cell cycle is a significant problem because there are several diseases associated with inadequate production and maintenance of blood cells. Multicellular organism development relies on a delicate balance between cell proliferation and cell death (apoptosis) [3]. While cell proliferation results in the development of various tissues and organs, cellular death leads to refinement of the various tissues and ultimately produces functional organs. Imbalances in cell proliferation and cell death result in a variety of diseases ranging from neurodegeneration to tumorigenesis [4, 5]. Cells undergoing serum (blood plasma not including fibrinogens) deprivation exhibit signs of apoptosis [6]. In some cases, cells subject to cytokines may show signs of proliferation [7]. In other cases, cytokines may exert a negative feedback on hematopoietic mechanisms in conditions of high cell counts in circulation [8]. As such, cytokines act as regulatory proteins that are released by cells of the immune system and control the selfrenewing process of stem cells [9]. There are numerous types of cytokines some of which are erythropoietin (EPO), granulocyte colonystimulating factor (GCSF), and Interleukin1 through Interleukin12. The cytokine EPO, for example, has been shown to intensify the production of red blood cells in response to oxygen levels in the body [10]. Interleukin3 has been shown to enhance the survival of HSCs while GCSF has been shown to contribute to mobilization of HSCs [11,12,13]. Based on previously published work, we develop a mathematical model that captures the effect of growth factor deprivation and cytokine administration [14, 15]. Our goal is to model the phenomenon observed in the presence of cytokines offsetting the negative effects of growth factor deprivation.
The model presented in this paper captures the behavior observed during serum deprivation. Based on flow cytometric studies conducted by others, the set of coupled differential delay equations captures apoptosis measured as PS plasma translocation levels of cells at different time points [16,17,18].
Additionally, we incorporate the effect of cytokines in our model. Cytokines have experimentally been observed to trigger an increase in proliferation of quiescent human stem cells [19]. Therefore, it is important to include a mathematical term to account for the effect of cytokines in the model equations in an effort to determine if cytokines can stabilize cell count. In order to better understand the dynamics of the model, phase plane analysis is carried out. Linear stability analysis is also performed in order to locate and determine stability of fixed points and simulations of the model are run in XPPAUT [20].
Materials and methods
The model
Using the model previously constructed and published in [15], we expand to include the effect of growth factor deprivation on a population of proliferating and nonproliferating cells. This model describes the number of proliferating and nonproliferating stem cells by a set of coupled delay differential equations. As opposed to our previous work, we have assumed a fixed oxygen concentration in the model. The modeling set of equations is:
In the above equations, P represents the number of proliferating stem cells and Nrepresents the number of nonproliferating cells. β(N) (eq. (4)) measures the rate of cell reentry into proliferation, β_{0} is the maximal rate of cell transit from resting phase to S phase, n measures the sensitivity of the rate of cell transit from G0 phase to S phase, and θ is the G0 stem cell population at which the rate of cell movement from G0 into proliferation is onehalf of its maximal value. δ is the rate of random cell loss by escape to the periphery and γ is the rate of cell loss due to apoptosis. τ is the time required for a cell to complete one cycle of the proliferation phase. The notation N_{ τ }, for example, represents N(t − τ), thus introducing a time delay into the equations. The values of parameters stated above are based on estimates given by in [21,22,23,24,25].
Growth factor deprivation contributes to cell death and, therefore, we subtract \( \overline{r}(s) rP \) in eq. (1). We define \( \overline{r}(s) \) as a maximal rate of cell death as a function of growth factor deprivation, s, and we define r as a uniformly distributed random variable from the interval (0, 0.0005) so that the variability of the effect of growth factor deprivation is within a .05% interval. We keep track of the number of apoptotic cells by introducing the variable A whose rate of change is given by eq. (3). Note that cell death can be attributed to both growth factor deprivation and apoptosis and determined through the parameter γ.
The term β_{ c }(P)N in eqs. (1) and (2) models the effect of cytokines on the proliferating and nonproliferating stem cells and is a Hill function. This term models the effect of cytokines by creating a cell loss of the nonproliferating population and a cell gain of the proliferating population when the proliferating cell population is low. This represents the transition of nonproliferating to proliferating cells as a result of cytokine administration. Since the effects of cytokines decay with time, we multiply by g_{ c }(t) which causes the effect of cytokine administration to exponentially decay. The parameter s_{1}determines the rate of decay of the cytokines and τ_{1}is a parameter that sets the time at which the effect of cytokines begins to decay. The parameter values in the term β_{ c }(P) were chosen to be consistent with realistic time frames of cytokines as noted in [26].
Steadystate solutions of the model and stability analysis
Steadystate solutions for the proliferating and nonproliferating cells have been calculated and linear stability analysis of eqs. (1) and (2) has been carried out as in [14, 15]. Two sets of steadystate solutions were found. The first set is P^{∗} = N^{∗} = 0. The second and more interesting steadystate solution is:
The nontrivial steady state exists only if
Also, \( \gamma =\overline{r}(s)r \) and \( \frac{\delta +{\beta}_{0,c}}{\gamma \overline{r}(s)r}>0 \).
Stable solutions of eqs. (1) and (2) exist only under the condition:
\( \omega \tau <{\cos}^{1}\left(\frac{A}{B}\right) \), where A = δ + β_{0}F + β_{0, c}, B = − 2β_{0}exp^{−γτ}F and \( F={\theta}^n\left[\frac{\theta +{N}^{\ast n}\left(1n\right)}{{\left({\theta}^n+{N}^{\ast n}\right)}^2}\right] \).
With the restriction that \( \left\frac{A}{B}\right<1. \).
Biologically, γ captures the apoptosis rate which has to be within a certain range for the system to remain active. The condition in eq. (9) reflects the fact that if apoptosis exceeds certain values, hematopoiesis will cease leading to death. Lastly, for homeostasis to prevail, the collection of parameter values must yield dynamic equilibrium (as given by eqs. (7) and (8)). Parameter values used in our model are given in Table 1 and were found to satisfy the above conditions.
Results
The system of eqs. (1)–(6) were solved numerically using the software XPPAUT under various conditions. See Figs. 1 and 2 below. Simulations were run in the baseline case meaning that growth factor deprivation and cytokine administration were not included in the model. The initial conditions were chosen so that P=1000, N=10,000 and A = 0. The ratio of P vs N is consistent with HSC measurements conducted using young adult bone marrow [27]. In this case, the number of proliferating stem cells exhibited damped oscillations and then reached a steady state value of approximately 775 cells. Next, simulations were run with growth factor deprivation included in the model. With growth factor deprivation, the number of proliferating stem cells continued to initially exhibit damped oscillations but at smaller values of the proliferating cell count. As was expected, the number of proliferating cells reached a smaller equilibrium value (approximately 674 cells) than in the baseline case. A comparison of the number of apoptotic cells in the simulations run in the baseline case versus simulations run with growth factor deprivation show that cell death occurred at a higher rate when growth factor deprivation was incorporated in the modeling set of equations. Finally, simulations were run with growth factor deprivation and cytokine administration both included in the model. In this case, cytokine administration was able to counter the adverse effect of growth factor deprivation. Cytokine administration boosts the proliferating cell count. With an appropriate strength (dosage) of the cytokine administration, the number of proliferating stem cells underdoing growth factor deprivation can be brought back to a value that is in close proximity to the baseline case. Fig. 1 shows a comparison of the number of proliferating stem cells in the three cases: baseline, with growth factor deprivation but not cytokine administration, and with both growth factor deprivation and cytokine administration. The number of proliferating cells is initially the same in the simulations. Serum deprivation causes the number of proliferating cells to decrease in comparison to the baseline case and cytokine administration causes the number of proliferating cells to shoot up in comparison to the baseline and serum deprived cases. As the number of proliferating cells settles towards the equilibrium (at approximately t=150 hrs), the baseline case and case with joint growth factor deprivation plus cytokine administration are relatively consistent with one another, while the case solely with growth factor deprivation shows lower numbers of proliferating cells. However, the effect of the cytokine begins to wear off (due to g_{ c }(t)) causing the number of proliferating cells to decay towards the same steady state value as is reached solely when serum deprivation is present.
In Fig. 2, it can be seen that with growth factor deprivation, the number of apoptotic cells increases as expected. When cytokine in administered in conjunction with growth factor deprivation, the number of proliferating cells greatly increases. This increase in the population of proliferating cells also means that there are more cells available to undergo random cell loss (as determined by γ). Thus, in addition to the number of proliferating cells increasing, the number of apoptotic cells also increases in the presence of cytokines.
Discussion
Based on the simulation results seen in Fig. 1, one can deduce that both cytokine administration and nocytokine administration will ultimately converge to the same steady levels provided the time period is long enough. This is consistent with homeostatic principles and the overall tendency of a hematopoietic system to restore itself to a normal state. Nonetheless, cytokine administration mobilizes stem cell proliferation in the early stages and results in higher number of proliferating cells in the early stages. Again, this observation is consistent to what is seen in patients who are treated with hematopoietic inducing agents. In similar manner, growth factor deprivation, a scenario typically seen in patients subjected to chemotherapy, there is an increase in the number of proliferating cells. In the first case, this is a result of the presence of proliferating inducing agent and in the latter a result of increased apoptosis due to a chemotherapeutic agent. However, in a scenario of increased apoptosis, the number of proliferating cells is decreased in the long term and unable to reach normal levels. As expected, a low number of proliferating cells at the beginning of the simulation yields a low number of proliferating cells at the end of the simulation.
The number of apoptotic cells is higher in both growth factor deprivation and growth factor deprivation with a cytokine compared with the baseline case, as seen in Fig. 2. The overall effect of the cytokine is to induce a higher number of proliferating cells that compensates for the loss of apoptotic cells. As such, the benefit of using a cytokine is not through the prevention of apoptosis under a stressful event but more to an increase in the number of proliferating cells.
Conclusion
In conclusion, the model presented exhibits an endpoint where the levels of proliferating stem cells in the presence of cytokine stabilize to a level that is comparable to levels seen in the absence of growth factor deprivation. The theoretical outcome of this model is that cytokines or any inducing hematopoietic agent allow for an increased value of proliferating cells but do not prevent apoptosis. Administration of cytokines during a stressful event, such as chemotherapy, is a practical way of maintaining equilibrium between apoptotic and proliferating.
Abbreviations
 EPO:

Erythropoietin
 FACS:

Flow Cytometry
 GCSF:

GranulocyteColony Stimulating Factor
 HSC:

Hematopoietic stem cells
References
 1.
Birbrair A, Frenette PS. Niche heterogeneity in the bone marrow. Ann N Y Acad Sci. 2016;1370(1):82–96.
 2.
Parslow TG, Stites DP, Terr AI, Imboden JB, Medical Immunology (1st ed.), McGraw Hill Professional; 2001.
 3.
Alenzi FQ, Alenazi BQ, Ahmad SY, Salem ML, AlJabri AA, Wyse RK. The haemopoietic stem cell: between apoptosis and self renewal. Yale J Biol Med. 2009;82(1):7–18.
 4.
Das M, Chatterjee S, Basak P, Das P, Pereira JA, Dutta RK, Chaklader M, Chaudhuri S, Law S. The bone marrow stem stromal imbalancea key feature of disease progression in case of myelodysplastic mouse model. J Stem Cells. 2010;5(2):49–64.
 5.
Lindvall O, Kokaia Z. Stem cells in human neurodegenerative disorders — time for clinical translation? J Clin Invest. 2010;120(1):29–40.
 6.
Zhu W, Chen J, Cong X, Hu S, Chen X. Hypoxia and growth factor deprivationinduced apoptosis in mesenchymal stem cells. Stem Cells. 2006;24(2):416–25.
 7.
Lotem J, Sachs L. Cytokine control of developmental programs in normal hematopoiesis and leukemia. Oncogene. 2002;21(21):3284–94.
 8.
Zhang CC, Lodish HF. Cytokines regulating hematopoietic stem cell function. Curr Opin Hematol. 2008;(4):307–11.
 9.
Lacy P. Secretion of cytokines and chemokines by innate immune cells. Front Immunol. 2015;6:190.
 10.
Ogunshola OO, Bogdanova AY. Epo and nonhematopoietic cells: what do we know? Methods Mol Biol. 2013;982:13–41.
 11.
Bryder D, Jacobsen SE. Interleukin3 supports expansion of longterm multilineage repopulating activity after multiple stem cell divisions in vitro. Blood. 2000;96(5):1748–55.
 12.
Bendall LJ, Bradstock KF. GCSF: from granulopoietic stimulant to bone marrow stem cell mobilizing agent. Cytokine Growth Factor Rev. 2014;25(4):355–67.
 13.
U. Duhrsen U, Villeval JL, Boyd J, Kannourakis G, Morstyn G, Metcalf D, Effects of recombinant human granulocyte colonystimulating factor on hematopoietic progenitor cells in cancer patients. Blood 1988;72: 2074–2081.
 14.
Mouser CL, Antoniou ES, Tadros J, Vassiliou EK. A model of hematopoietic stem cell proliferation under the influence of a chemotherapeutic agent in combination with a hematopoietic inducing agent. Theor Biol Med Model. 2014;11:4.
 15.
Antoniou ES, Mouser CL, Rosar ME, Tadros J, Vassiliou EK. Hematopoietic stem cell proliferation modeling under the influence of hematopoieticinducing agent. Shock. 2009;32(5):471–7.
 16.
Wlodkowic W, Skommer J, Darzynkiewicz Z. Flow cytometrybased apoptosis detection. Methods Mol Biol. 2009;559:10.
 17.
Van Engeland M, Nieland LJ, Ramaekers FC, Schutte B, Reutelingsperger CP. Annexin Vaffinity assay: a review on an apoptosis detection system based on phosphatidylserine exposure. Cytometry. 1998;31(1):1–9.
 18.
Zhu W, Chen J, Cong X, Hu S, Chen X. Hypoxia and serum deprivationinduced apoptosis in mesenchymal stem cells. Stem Cells. 2006;24:416–25.
 19.
Shah AJ, Smogorzewska EM, Hannum C, Crooks GM. Flt3 ligand induces proliferation of quiescent human bone marrow CD34+CD38 cells and maintains progenitor cells in vitro. Blood. 1996;87(9):3563–70.
 20.
Ermentrout B. Simulating, analyzing, and animating dynamical systems : a guide to XPPAUT for researchers and students. Philadelphia, USA: Society for Industrial and Applied Mathematics; 2002.
 21.
Bernard S, Belair J, Mackey MC. Bifurcation in a whitebloodcell production model. Comptes Rendus Biol. 2004;372:201–10.
 22.
Bernard S, Belair J, Mackey MC. Oscillations in cyclical neutropenia: new evidence based on mathematical modeling. J Theor Biol. 2001;223:283–98.
 23.
Mackey MC. Cell kinetic status of haematopoietic stem cells. Cell Prolif. 2001;32:71–83.
 24.
Berlier JL, Rigutto S, Dalla Valle A, Lechanteur J, Soyfoo MS, Gangji V, Rasschaert J. Adenosine riphosphate prevents serum deprivationinduced apoptosis in human mesenchymal stem cells via activation of the MAPK signaling pathways. Stem Cells. 2015;33(1):211–8.
 25.
Brugger W, Wolfgang M, Heimfeld S, Berenson RJ, Mertelsmann R, Kanz L. Ex vivo expansion of enriched peripheral blood CD34+ progenitor cells by stem cell factor, interleukinlβ (ILlβ), IL6, IL3, interferonγ, and erythropoietin. Blood. 1993;81(10):2579–84.
 26.
Yan X, Lowe PJ, Fink M, Berghout A, Balser S, Krzyzanski W. Population pharmacokinetic and pharmacodynamic modelbased comparability assessment of a recombinant human epoetin alfa and the biosimilar HX575. J Clin Pharmacol. 2012;52(11):1624–44.
 27.
Pang WW, Price EA, Sahoo D, Beerman I, Maloney WJ, Rossi DJ, Schrier SL, Weissman IL. Human bone marrow hematopoietic stem cells are increased in frequency and myeloidbiased with age. Proc Natl Acad Sci U S A. 2011;108(50):20012–7.
Funding
This work was supported by grants from the Center for Research, College of Science and Health, William Paterson University of New Jersey and from the Office of Research and Sponsored Programs of Kean University.
Availability of data and materials
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Author information
Affiliations
Contributions
CM derived the terms to be added to the mathematical model. CM ran the numerical simulations and determined the parameter values for the model. EA calculated the fixed points and stability analysis. EV provided insight and clinical data on biology and did literature search for clinical data. CM and EV wrote the document together. EA verified the equations and proofread the document. All authors read and approved the final manuscript.
Corresponding author
Correspondence to Christina L. Mouser.
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Hematopoiesis
 Growth factor deprivation
 Differential delay equations