Merging concepts - coupling an agent-based model of hematopoietic stem cells with an ODE model of granulopoiesis

BackgroundHematopoiesis is a complex process involving different cell types and feedback mechanisms mediated by cytokines. This complexity stimulated various models with different scopes and applications. A combination of complementary models promises to provide their mutual confirmation and to explain a broader range of scenarios. Here we propose a combination of an ordinary differential equation (ODE) model of human granulopoiesis and an agent-based model (ABM) of hematopoietic stem cell (HSC) organization. The first describes the dynamics of bone marrow cell stages and circulating cells under various perturbations such as G-CSF treatment or chemotherapy. In contrast to the ODE model describing cell numbers, our ABM focuses on the organization of individual cells in the stem population.ResultsWe combined the two models by replacing the HSC compartment of the ODE model by a difference equation formulation of the ABM. In this hybrid model, regulatory mechanisms and parameters of the original models were kept unchanged except for a few specific improvements: (i) Effect of chemotherapy was restricted to proliferating HSC and (ii) HSC regulation in the ODE model was replaced by the intrinsic regulation of the ABM. Model simulations of bleeding, chronic irradiation and stem cell transplantation revealed that the dynamics of hybrid and ODE model differ markedly in scenarios with stem cell damage. Despite these differences in response to stem cell damage, both models explain clinical data of leukocyte dynamics under four chemotherapy regimens.ConclusionsABM and ODE model proved to be compatible and were combined without altering the structure of both models. The new hybrid model introduces model improvements by considering the proliferative state of stem cells and enabling a cell cycle-dependent effect of chemotherapy. We demonstrated that it is able to explain and predict granulopoietic dynamics for a large variety of scenarios such as irradiation, bone marrow transplantation, chemotherapy and growth factor applications. Therefore, it promises to serve as a valuable tool for studies in a broader range of clinical applications, in particular where stem cell activation and proliferation are involved.

[1]  M. Goodell,et al.  Hematopoietic potential of stem cells isolated from murine skeletal muscle. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[2]  Ingo Roeder,et al.  Characterization and quantification of clonal heterogeneity among hematopoietic stem cells: a model-based approach. , 2008, Blood.

[3]  Fabien Crauste,et al.  A Mathematical Study of the Hematopoiesis Process with Applications to Chronic Myelogenous Leukemia , 2009, SIAM J. Appl. Math..

[4]  Sebastian Bonhoeffer,et al.  Dynamic variation in cycling of hematopoietic stem cells in steady state and inflammation , 2011, The Journal of experimental medicine.

[5]  Doron Levy,et al.  Modeling Imatinib-Treated Chronic Myelogenous Leukemia: Reducing the Complexity of Agent-Based Models , 2008, Bulletin of mathematical biology.

[6]  J. Till,et al.  A direct measurement of the radiation sensitivity of normal mouse bone marrow cells. , 1961, Radiation research.

[7]  M. Entman,et al.  Stem Cell Plasticity in Muscle and Bone Marrow , 2001, Annals of the New York Academy of Sciences.

[8]  Jinzhi Lei,et al.  Neutrophil dynamics in response to chemotherapy and G-CSF. , 2012, Journal of theoretical biology.

[9]  J. Chute,et al.  Transplantation of vascular endothelial cells mediates the hematopoietic recovery and survival of lethally irradiated mice. , 2007, Blood.

[10]  A. Vescovi,et al.  Turning brain into blood: a hematopoietic fate adopted by adult neural stem cells in vivo. , 1999, Science.

[11]  Andreas Trumpp,et al.  Hematopoietic Stem Cells Reversibly Switch from Dormancy to Self-Renewal during Homeostasis and Repair , 2008, Cell.

[12]  L. Lajtha,et al.  Kinetic Model of a Bone‐Marrow Stem‐Cell Population , 1962, British journal of haematology.

[13]  R. Schofield The relationship between the spleen colony-forming cell and the haemopoietic stem cell. , 1978, Blood cells.

[14]  V. Vainsteina,et al.  The complex effect of granulocyte colony-stimulating factor on human granulopoiesis analyzed by a new physiologically-based mathematical model , 2005 .

[15]  M. Tavassoli,et al.  Alterations of bone marrow sinus endothelium induced by ionizing irradiation: implications in the homing of intravenously transplanted marrow cells. , 1992, Blood cells.

[16]  O. Ringdén,et al.  A Phase 3, Randomized, Placebo-controlled Trial of Filgrastim in Patients with Haematological Malignancies Undergoing Matched-related Allogeneic Bone Marrow Transplantation , 2008, Archives of drug information.

[17]  S. Lenhoff,et al.  Granulocyte interactions with GM-CSF and G-CSF secretion by endothelial cells and monocytes. , 1999, European cytokine network.

[18]  Olivier Gandrillon,et al.  Mathematical study of feedback control roles and relevance in stress erythropoiesis. , 2010, Journal of theoretical biology.

[19]  C. D. Dunn Current concepts in erythropoiesis , 1983 .

[20]  Dirk Hasenclever,et al.  Two-weekly or 3-weekly CHOP chemotherapy with or without etoposide for the treatment of elderly patients with aggressive lymphomas: results of the NHL-B2 trial of the DSHNHL. , 2004, Blood.

[21]  W. Fleming,et al.  Adult blood vessels restore host hematopoiesis following lethal irradiation. , 2002, Experimental hematology.

[22]  Markus Scholz,et al.  A biomathematical model of human thrombopoiesis under chemotherapy. , 2010, Journal of theoretical biology.

[23]  K. Moore,et al.  Stem Cells and Their Niches , 2006, Science.

[24]  M. Loeffler,et al.  Stem Cell Fate Analysis Revisited: Interpretation of Individual Clone Dynamics in the Light of a New Paradigm of Stem Cell Organization , 2007, Journal of biomedicine & biotechnology.

[25]  Elisa Laurenti,et al.  Hematopoiesis: a human perspective. , 2012, Cell stem cell.

[26]  L. Sachs,et al.  In vitro control of the development of macrophage and granulocyte colonies. , 1966, Proceedings of the National Academy of Sciences of the United States of America.

[27]  A Howell,et al.  The kinetics of human granulopoiesis following treatment with granulocyte colony-stimulating factor in vivo. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[28]  Markus Scholz,et al.  Pharmacokinetic and -dynamic modelling of G-CSF derivatives in humans , 2012, Theoretical Biology and Medical Modelling.

[29]  Ivar Østby,et al.  Mathematical modeling of granulocyte reconstitution after high-dose chemotherapy with stem cell support: effect of post-transplant G-CSF treatment. , 2004, Journal of theoretical biology.

[30]  M. Scholz,et al.  Modelling Human Granulopoiesis under Poly-chemotherapy with G-CSF Support , 2005, Journal of mathematical biology.

[31]  E. Cronkite,et al.  Structural, cytologic and autoradiographic (H3-thymidine) changes in the bone marrow following total body irradiation. , 1961, The American journal of pathology.

[32]  M C Mackey,et al.  Unified hypothesis for the origin of aplastic anemia and periodic hematopoiesis. , 1978, Blood.

[33]  M. Loeffler,et al.  Somatic mutation, monoclonality and stochastic models of stem cell organization in the intestinal crypt. , 1993, Journal of theoretical biology.

[34]  C. Schmoor,et al.  Influence of recombinant human granulocyte colony-stimulating factor (filgrastim) on hematopoietic recovery and outcome following allogeneic bone marrow transplantation (BMT) from volunteer unrelated donors , 1999, Bone Marrow Transplantation.

[35]  Jinzhi Lei,et al.  Neutrophil dynamics after chemotherapy and G-CSF: the role of pharmacokinetics in shaping the response. , 2012, Journal of theoretical biology.

[36]  Willi Jäger,et al.  Modeling of asymmetric cell division in hematopoietic stem cells--regulation of self-renewal is essential for efficient repopulation. , 2009, Stem cells and development.

[37]  K. Akashi,et al.  Hematopoietic developmental pathways: on cellular basis , 2007, Oncogene.

[38]  Dirk Hasenclever,et al.  Two-weekly or 3-weekly CHOP chemotherapy with or without etoposide for the treatment of young patients with good-prognosis (normal LDH) aggressive lymphomas: results of the NHL-B1 trial of the DSHNHL. , 2004, Blood.

[39]  A. Fauci,et al.  Alternate-day prednisone. Leukocyte kinetics and susceptibility to infections. , 1974, The New England journal of medicine.

[40]  I. Glauche,et al.  Cellular aging leads to functional heterogeneity of hematopoietic stem cells: a modeling perspective , 2011, Aging cell.

[41]  A. Trumpp,et al.  Therapy of chronic myeloid leukaemia can benefit from the activation of stem cells: simulation studies of different treatment combinations , 2012, British Journal of Cancer.

[42]  A. Alexander DER LYMPHOZYT ALS GEMEINSAME STAMMZELLE DER VERSCHIEDENEN BLUTELEMENTE IN DER EMBRYONALEN ENTWICKLUNG UND IM POSTFETALEN LEBEN DER SäUGETIERE , 2009 .

[43]  Markus Scholz,et al.  A computational model of human granulopoiesis to simulate the hematotoxic effects of multicycle polychemotherapy. , 2004, Blood.

[44]  Ivar Østby,et al.  A mathematical model for reconstitution of granulopoiesis after high dose chemotherapy with autologous stem cell transplantation , 2003, Journal of mathematical biology.

[45]  Michael C Mackey,et al.  A mathematical model of hematopoiesis: II. Cyclical neutropenia. , 2005, Journal of theoretical biology.

[46]  Laurent Pujo-Menjouet,et al.  Adding self-renewal in committed erythroid progenitors improves the biological relevance of a mathematical model of erythropoiesis. , 2008, Journal of theoretical biology.

[47]  Michael Cross,et al.  Lineage Specification of Hematopoietic Stem Cells: Mathematical Modeling and Biological Implications , 2007, Stem cells.

[48]  M. Loeffler,et al.  Stem cells: attributes, cycles, spirals, pitfalls and uncertainties. Lessons for and from the crypt. , 1990, Development.

[49]  H E Wichmann,et al.  Quantification of the cell kinetic effects of G-CSF using a model of human granulopoiesis. , 1993, Experimental hematology.

[50]  M M Wintrobe,et al.  Leukokinetic studies. 13. A non-steady-state kinetic evaluation of the mechanism of cortisone-induced granulocytosis. , 1968, The Journal of clinical investigation.

[51]  M. Loeffler,et al.  Mathematical Modeling of Genesis and Treatment of Chronic Myeloid Leukemia , 2008, Cells Tissues Organs.

[52]  Priv.-Doz. Dr. Hans-Peter Lohrmann,et al.  Cytotoxic Drugs and the Granulopoietic System , 1981, Recent Results in Cancer Research.

[53]  Michael C Mackey,et al.  A mathematical model of hematopoiesis--I. Periodic chronic myelogenous leukemia. , 2005, Journal of theoretical biology.

[54]  C. Peterson,et al.  Dynamical modelling of haematopoiesis: an integrated view over the system in homeostasis and under perturbation , 2013, Journal of The Royal Society Interface.

[55]  Michael C Mackey,et al.  G-CSF treatment of canine cyclical neutropenia: a comprehensive mathematical model. , 2007, Experimental hematology.

[56]  Ingo Roeder,et al.  Dynamic modeling of imatinib-treated chronic myeloid leukemia: functional insights and clinical implications , 2006, Nature Medicine.

[57]  Ingo Roeder,et al.  A novel dynamic model of hematopoietic stem cell organization based on the concept of within-tissue plasticity. , 2002, Experimental hematology.