1
|
Mackey MC. Periodic hematological disorders: Quintessential examples of dynamical diseases. CHAOS (WOODBURY, N.Y.) 2020; 30:063123. [PMID: 32611100 DOI: 10.1063/5.0006517] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/04/2020] [Accepted: 05/21/2020] [Indexed: 06/11/2023]
Abstract
This paper summarizes the evidence supporting the classification of cyclic neutropenia as a dynamical disease and periodic chronic myelogenous leukemia is also considered. The unsatisfactory state of knowledge concerning the genesis of cyclic thrombocytopenia and periodic autoimmune hemolytic anemia is detailed.
Collapse
Affiliation(s)
- Michael C Mackey
- Department of Physiology, Department of Physics, and Department of Mathematics McGill University, Montreal, Quebec H4X 2C1, Canada
| |
Collapse
|
2
|
Andersen M, Hasselbalch HC, Kjær L, Skov V, Ottesen JT. Global dynamics of healthy and cancer cells competing in the hematopoietic system. Math Biosci 2020; 326:108372. [PMID: 32442449 DOI: 10.1016/j.mbs.2020.108372] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2019] [Revised: 05/06/2020] [Accepted: 05/06/2020] [Indexed: 01/08/2023]
Abstract
Stem cells in the bone marrow differentiate to ultimately become mature, functioning blood cells through a tightly regulated process (hematopoiesis) including a stem cell niche interaction and feedback through the immune system. Mutations in a hematopoietic stem cell can create a cancer stem cell leading to a less controlled production of malfunctioning cells in the hematopoietic system. This was mathematically modelled by Andersen et al. (2017) including the dynamic variables: healthy and cancer stem cells and mature cells, dead cells and an immune system response. Here, we apply a quasi steady state approximation to this model to construct a two dimensional model with four algebraic equations denoted the simple cancitis model. The two dynamic variables are the clinically available quantities JAK2V617F allele burden and the number of white blood cells. The simple cancitis model represents the original model very well. Complete phase space analysis of the simple cancitis model is performed, including proving the existence and location of globally attracting steady states. Hence, parameter values from compartments of stem cells, mature cells and immune cells are directly linked to disease and treatment prognosis, showing the crucial importance of early intervention. The simple cancitis model allows for a complete analysis of the long term evolution of trajectories. In particular, the value of the self renewal of the hematopoietic stem cells divided by the self renewal of the cancer stem cells is found to be an important diagnostic marker and perturbing this parameter value at intervention allows the model to reproduce clinical data. Treatment at low cancer cell numbers allows returning to healthy blood production while the same intervention at a later disease stage can lead to eradication of healthy blood producing cells. Assuming the total number of white blood cells is constant in the early cancer phase while the allele burden increases, a one dimensional model is suggested and explicitly solved, including parameters from all original compartments. The solution explicitly shows that exogenous inflammation promotes blood cancer when cancer stem cells reproduce more efficiently than hematopoietic stem cells.
Collapse
Affiliation(s)
- Morten Andersen
- IMFUFA, Department of Science and Environment, Roskilde University, Denmark.
| | - Hans C Hasselbalch
- Department of Haematology, Zealand University Hospital, Roskilde, Denmark
| | - Lasse Kjær
- Department of Haematology, Zealand University Hospital, Roskilde, Denmark
| | - Vibe Skov
- Department of Haematology, Zealand University Hospital, Roskilde, Denmark
| | - Johnny T Ottesen
- IMFUFA, Department of Science and Environment, Roskilde University, Denmark
| |
Collapse
|
3
|
Lei J. A general mathematical framework for understanding the behavior of heterogeneous stem cell regeneration. J Theor Biol 2020; 492:110196. [PMID: 32067937 DOI: 10.1016/j.jtbi.2020.110196] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2019] [Revised: 12/28/2019] [Accepted: 02/11/2020] [Indexed: 11/21/2022]
Abstract
Stem cell heterogeneity is essential for homeostasis in tissue development. This paper establishes a general mathematical framework to model the dynamics of stem cell regeneration with cell heterogeneity and random transitions of epigenetic states. The framework generalizes the classical G0 cell cycle model and incorporates the epigenetic states of individual cells represented by a continuous multidimensional variable. In the model, the kinetic rates of cell behaviors, including proliferation, differentiation, and apoptosis, are dependent on their epigenetic states, and the random transitions of epigenetic states between cell cycles are represented by an inheritance probability function that describes the conditional probability of cell state changes. Moreover, the model can be extended to include genotypic changes and describe the process of gene mutation-induced tumor development. The proposed mathematical framework provides a generalized formula that helps us to understand various dynamic processes of stem cell regeneration, including tissue development, degeneration, and abnormal growth.
Collapse
Affiliation(s)
- Jinzhi Lei
- Zhou Pei-Yuan Center for Applied Mathematics, MOE Key Laboratory of Bioinformatics, Tsinghua University, Beijing 100084, China.
| |
Collapse
|
4
|
Mir P, Klimiankou M, Findik B, Hähnel K, Mellor-Heineke S, Zeidler C, Skokowa J, Welte K. New insights into the pathomechanism of cyclic neutropenia. Ann N Y Acad Sci 2020; 1466:83-92. [PMID: 32083314 DOI: 10.1111/nyas.14309] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2019] [Revised: 12/27/2019] [Accepted: 01/14/2020] [Indexed: 01/10/2023]
Abstract
Cyclic neutropenia (CyN) is a hematologic disorder in which peripheral blood absolute neutrophil counts (ANCs) show cycles of approximately 21-day intervals. The majority of CyN patients harbor ELANE mutations, but the mechanism of ANC cycling is unclear. We performed analysis of bone marrow (BM) subpopulations in CyN patients at the peak and the nadir of the ANC cycle and detected high proportions of BM hematopoietic stem cells (HSCs) and hematopoietic stem and progenitor cells (HSPCs) at the nadir of the ANC cycle, as compared with the peak. BM HSPCs produced fewer granulocyte colony-forming unit colonies at the ANC peak. To investigate the mechanism of cycling, we found that mRNA expression levels of ELANE and unfolded protein response (UPR)-related genes (ATF6, BiP (HSPA5), CHOP (DDIT3), and PERK (EIF2AK3)) were elevated, but antiapoptotic genes (Bcl-2 (BCL2) and bcl-xL (BCL2L1)) were reduced in CD34+ cells tested at the ANC nadir. Moreover, HSPCs revealed increased levels of reactive oxygen species and gH2AX at the ANC nadir. We suggest that in CyN patients, some HSPCs escape the UPR-induced endoplasmic reticulum (ER) stress and proliferate in response to granulocyte colony-stimulating factor (G-CSF) to a certain threshold at which UPR again affects the majority of HSPCs. There is a cyclic balance between ER stress-induced apoptosis of HSPCs and compensatory G-CSF-stimulated HSPC proliferation followed by granulocytic differentiation.
Collapse
Affiliation(s)
- Perihan Mir
- Department of Oncology, Hematology, Immunology, Rheumatology and Pulmonology, University Hospital Tübingen, Tübingen, Germany
| | - Maksim Klimiankou
- Department of Oncology, Hematology, Immunology, Rheumatology and Pulmonology, University Hospital Tübingen, Tübingen, Germany
| | - Betuel Findik
- Department of Oncology, Hematology, Immunology, Rheumatology and Pulmonology, University Hospital Tübingen, Tübingen, Germany
| | - Karin Hähnel
- Department of Oncology, Hematology, Immunology, Rheumatology and Pulmonology, University Hospital Tübingen, Tübingen, Germany
| | - Sabine Mellor-Heineke
- Department of Oncology, Hematology and Bone Marrow Transplantation, Hannover Medical School, Hannover, Germany
| | - Cornelia Zeidler
- Department of Oncology, Hematology and Bone Marrow Transplantation, Hannover Medical School, Hannover, Germany
| | - Julia Skokowa
- Department of Oncology, Hematology, Immunology, Rheumatology and Pulmonology, University Hospital Tübingen, Tübingen, Germany
| | - Karl Welte
- Department of Pediatric Hematology, Oncology and Bone Marrow Transplantation, University Children's Hospital Tübingen, Tübingen, Germany
| |
Collapse
|
5
|
Knauer F, Stiehl T, Marciniak-Czochra A. Oscillations in a white blood cell production model with multiple differentiation stages. J Math Biol 2019; 80:575-600. [DOI: 10.1007/s00285-019-01432-6] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2019] [Revised: 07/02/2019] [Indexed: 12/15/2022]
|
6
|
Zhuge C, Mackey MC, Lei J. Origins of oscillation patterns in cyclical thrombocytopenia. J Theor Biol 2019; 462:432-445. [DOI: 10.1016/j.jtbi.2018.11.024] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Revised: 11/22/2018] [Accepted: 11/26/2018] [Indexed: 10/27/2022]
|
7
|
Situ Q, Lei J. A mathematical model of stem cell regeneration with epigenetic state transitions. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2018; 14:1379-1397. [PMID: 29161866 DOI: 10.3934/mbe.2017071] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
In this paper, we study a mathematical model of stem cell regeneration with epigenetic state transitions. In the model, the heterogeneity of stem cells is considered through the epigenetic state of each cell, and each epigenetic state defines a subpopulation of stem cells. The dynamics of the subpopulations are modeled by a set of ordinary differential equations in which epigenetic state transition in cell division is given by the transition probability. We present analysis for the existence and linear stability of the equilibrium state. As an example, we apply the model to study the dynamics of state transition in breast cancer stem cells.
Collapse
Affiliation(s)
- Qiaojun Situ
- Zhou Pei-Yuan Center for Applied Mathematics, Tsinghua University, Beijing 100084, China
| | - Jinzhi Lei
- Zhou Pei-Yuan Center for Applied Mathematics, MOE Key Laboratory of Bioinformatics, Tsinghua University, Beijing 100084, China
| |
Collapse
|
8
|
Abstract
Severe congenital neutropenias are a heterogeneous group of rare haematological diseases characterized by impaired maturation of neutrophil granulocytes. Patients with severe congenital neutropenia are prone to recurrent, often life-threatening infections beginning in their first months of life. The most frequent pathogenic defects are autosomal dominant mutations in ELANE, which encodes neutrophil elastase, and autosomal recessive mutations in HAX1, whose product contributes to the activation of the granulocyte colony-stimulating factor (G-CSF) signalling pathway. The pathophysiological mechanisms of these conditions are the object of extensive research and are not fully understood. Furthermore, severe congenital neutropenias may predispose to myelodysplastic syndromes or acute myeloid leukaemia. Molecular events in the malignant progression include acquired mutations in CSF3R (encoding G-CSF receptor) and subsequently in other leukaemia-associated genes (such as RUNX1) in a majority of patients. Diagnosis is based on clinical manifestations, blood neutrophil count, bone marrow examination and genetic and immunological analyses. Daily subcutaneous G-CSF administration is the treatment of choice and leads to a substantial increase in blood neutrophil count, reduction of infections and drastic improvement of quality of life. Haematopoietic stem cell transplantation is the alternative treatment. Regular clinical assessments (including yearly bone marrow examinations) to monitor treatment course and detect chromosomal abnormalities (for example, monosomy 7 and trisomy 21) as well as somatic pre-leukaemic mutations are recommended.
Collapse
Affiliation(s)
- Julia Skokowa
- Department of Hematology, Oncology, Clinical Immunology, University of Tübingen, Tübingen, Germany
| | - David C Dale
- Department of Medicine, University of Washington, Seattle, Washington, USA
| | - Ivo P Touw
- Department of Hematology, Erasmus University Medical Center, Rotterdam, The Netherlands
| | - Cornelia Zeidler
- Department of Hematology and Oncology, Medical School Hannover, Hannover, Germany
| | - Karl Welte
- University Children's Hospital, Department of General Pediatrics and Pediatric Hematology and Oncology, Hoppe-Seyler-Str. 1, Tübingen 72076, Germany
| |
Collapse
|
9
|
Craig M, Humphries AR, Mackey MC. A Mathematical Model of Granulopoiesis Incorporating the Negative Feedback Dynamics and Kinetics of G-CSF/Neutrophil Binding and Internalization. Bull Math Biol 2016; 78:2304-2357. [PMID: 27324993 DOI: 10.1007/s11538-016-0179-8] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2015] [Accepted: 05/19/2016] [Indexed: 11/24/2022]
Abstract
We develop a physiological model of granulopoiesis which includes explicit modelling of the kinetics of the cytokine granulocyte colony-stimulating factor (G-CSF) incorporating both the freely circulating concentration and the concentration of the cytokine bound to mature neutrophils. G-CSF concentrations are used to directly regulate neutrophil production, with the rate of differentiation of stem cells to neutrophil precursors, the effective proliferation rate in mitosis, the maturation time, and the release rate from the mature marrow reservoir into circulation all dependent on the level of G-CSF in the system. The dependence of the maturation time on the cytokine concentration introduces a state-dependent delay into our differential equation model, and we show how this is derived from an age-structured partial differential equation model of the mitosis and maturation and also detail the derivation of the rest of our model. The model and its estimated parameters are shown to successfully predict the neutrophil and G-CSF responses to a variety of treatment scenarios, including the combined administration of chemotherapy and exogenous G-CSF. This concomitant treatment was reproduced without any additional fitting to characterize drug-drug interactions.
Collapse
Affiliation(s)
- M Craig
- Faculté de Pharmacie, Université de Montréal, Montréal, QC, H3T 1J4, Canada.
- Program for Evolutionary Dynamics, Harvard University, Cambridge, MA, 02138, USA.
| | - A R Humphries
- Department of Mathematics and Statistics, McGill University, Montréal, QC, H3A 0B9, Canada
| | - M C Mackey
- Departments of Mathematics, Physics and Physiology, McGill University, Montréal, QC, H3G 1Y6, Canada
| |
Collapse
|
10
|
Craig M, Humphries AR, Nekka F, Bélair J, Li J, Mackey MC. Neutrophil dynamics during concurrent chemotherapy and G-CSF administration: Mathematical modelling guides dose optimisation to minimise neutropenia. J Theor Biol 2015; 385:77-89. [PMID: 26343861 DOI: 10.1016/j.jtbi.2015.08.015] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2014] [Revised: 06/10/2015] [Accepted: 08/20/2015] [Indexed: 11/18/2022]
Abstract
The choice of chemotherapy regimens is often constrained by the patient's tolerance to the side effects of chemotherapeutic agents. This dose-limiting issue is a major concern in dose regimen design, which is typically focused on maximising drug benefits. Chemotherapy-induced neutropenia is one of the most prevalent toxic effects patients experience and frequently threatens the efficient use of chemotherapy. In response, granulocyte colony-stimulating factor (G-CSF) is co-administered during chemotherapy to stimulate neutrophil production, increase neutrophil counts, and hopefully avoid neutropenia. Its clinical use is, however, largely dictated by trial and error processes. Based on up-to-date knowledge and rational considerations, we develop a physiologically realistic model to mathematically characterise the neutrophil production in the bone marrow which we then integrate with pharmacokinetic and pharmacodynamic (PKPD) models of a chemotherapeutic agent and an exogenous form of G-CSF (recombinant human G-CSF, or rhG-CSF). In this work, model parameters represent the average values for a general patient and are extracted from the literature or estimated from available data. The dose effect predicted by the model is confirmed through previously published data. Using our model, we were able to determine clinically relevant dosing regimens that advantageously reduce the number of rhG-CSF administrations compared to original studies while significantly improving the neutropenia status. More particularly, we determine that it could be beneficial to delay the first administration of rhG-CSF to day seven post-chemotherapy and reduce the number of administrations from ten to three or four for a patient undergoing 14-day periodic chemotherapy.
Collapse
Affiliation(s)
- Morgan Craig
- Faculté de Pharmacie, Université de Montréal, Montréal, QC, Canada H3C 3J7; Centre for Applied Mathematics in Bioscience and Medicine (CAMBAM), McGill University, Montreal, QC, Canada H3G 1Y6.
| | - Antony R Humphries
- Department of Mathematics and Statistics, McGill University, Montreal, QC, Canada H3A 0B9; Centre for Applied Mathematics in Bioscience and Medicine (CAMBAM), McGill University, Montreal, QC, Canada H3G 1Y6; Centre de recherches mathématiques, Université de Montréal, Montréal, QC, Canada H3C 3J7.
| | - Fahima Nekka
- Centre for Applied Mathematics in Bioscience and Medicine (CAMBAM), McGill University, Montreal, QC, Canada H3G 1Y6; Centre de recherches mathématiques, Université de Montréal, Montréal, QC, Canada H3C 3J7.
| | - Jacques Bélair
- Département de mathématiques et de statistique, Université de Montréal, Montréal, QC, Canada H3C 3J7; Centre for Applied Mathematics in Bioscience and Medicine (CAMBAM), McGill University, Montreal, QC, Canada H3G 1Y6; Centre de recherches mathématiques, Université de Montréal, Montréal, QC, Canada H3C 3J7.
| | - Jun Li
- Faculté de Pharmacie, Université de Montréal, Montréal, QC, Canada H3C 3J7; Centre for Applied Mathematics in Bioscience and Medicine (CAMBAM), McGill University, Montreal, QC, Canada H3G 1Y6; Centre de recherches mathématiques, Université de Montréal, Montréal, QC, Canada H3C 3J7.
| | - Michael C Mackey
- Department of Mathematics and Statistics, McGill University, Montreal, QC, Canada H3A 0B9; Centre for Applied Mathematics in Bioscience and Medicine (CAMBAM), McGill University, Montreal, QC, Canada H3G 1Y6; Departments of Physiology and Physics, McGill University, Montreal, QC, Canada H3G 1Y6.
| |
Collapse
|
11
|
Bernard S, Crauste F. Optimal linear stability condition for scalar differential equations with distributed delay. ACTA ACUST UNITED AC 2015. [DOI: 10.3934/dcdsb.2015.20.1855] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
12
|
Abstract
Recently, there has been significant activity in the mathematical community, aimed at developing quantitative tools for studying leukemia and lymphoma. Mathematical models have been applied to evaluate existing therapies and to suggest novel therapies. This article reviews the recent contributions of mathematical modeling to leukemia and lymphoma research. These developments suggest that mathematical modeling has great potential in this field. Collaboration between mathematicians, clinicians, and experimentalists can significantly improve leukemia and lymphoma therapy.
Collapse
Affiliation(s)
- Geoffrey Clapp
- Department of Mathematics, University of Maryland, College Park, MD 20742
| | - Doron Levy
- Department of Mathematics, University of Maryland, College Park, MD 20742; Center for Scientific Computation and Mathematical Modeling (CSCAMM), University of Maryland, College Park, MD 20742, USA
| |
Collapse
|
13
|
Dale DC, Mackey MC. Understanding, treating and avoiding hematological disease: better medicine through mathematics? Bull Math Biol 2014; 77:739-57. [PMID: 25213154 DOI: 10.1007/s11538-014-9995-x] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2014] [Accepted: 07/08/2014] [Indexed: 10/24/2022]
Abstract
This paper traces the experimental, clinical and mathematical modeling efforts to understand a periodic hematological disease-cyclical neutropenia. It is primarily a highly personal account by two scientists from quite different backgrounds of their interactions over almost 40 years and their attempts to understand this intriguing disease. It's also a story of their efforts to offer effective treatments for the patients who suffer from cyclic neutropenia and other conditions causing neutropenia and infections.
Collapse
Affiliation(s)
- David C Dale
- Department of Medicine, University of Washington, Seattle, WA, 98195, USA,
| | | |
Collapse
|
14
|
Sánchez-Taltavull D, Alarcón T. Robustness of differentiation cascades with symmetric stem cell division. J R Soc Interface 2014; 11:20140264. [PMID: 24718457 DOI: 10.1098/rsif.2014.0264] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Stem cells (SCs) perform the task of maintaining tissue homeostasis by both self-renewal and differentiation. While it has been argued that SCs divide asymmetrically, there is also evidence that SCs undergo symmetric division. Symmetric SC division has been speculated to be key for expanding cell numbers in development and regeneration after injury. However, it might lead to uncontrolled growth and malignancies such as cancer. In order to explore the role of symmetric SC division, we propose a mathematical model of the effect of symmetric SC division on the robustness of a population regulated by a serial differentiation cascade and we show that this may lead to extinction of such population. We examine how the extinction likelihood depends on defining characteristics of the population such as the number of intermediate cell compartments. We show that longer differentiation cascades are more prone to extinction than systems with less intermediate compartments. Furthermore, we have analysed the possibility of mixed symmetric and asymmetric cell division against invasions by mutant invaders in order to find optimal architecture. Our results show that more robust populations are those with unfrequent symmetric behaviour.
Collapse
Affiliation(s)
- Daniel Sánchez-Taltavull
- Centre de Recerca Matemàtica, , Edifici C, Campus de Bellaterra, 08193 Bellaterra (Barcelona), Spain
| | | |
Collapse
|
15
|
Understanding and Treating Cytopenia Through Mathematical Modeling. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2014; 844:279-302. [DOI: 10.1007/978-1-4939-2095-2_14] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
16
|
Adimy M, Angulo O, Marquet C, Sebaa L. A mathematical model of multistage hematopoietic cell lineages. ACTA ACUST UNITED AC 2014. [DOI: 10.3934/dcdsb.2014.19.1] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
17
|
Brooks G, Provencher G, Lei J, Mackey MC. Neutrophil dynamics after chemotherapy and G-CSF: the role of pharmacokinetics in shaping the response. J Theor Biol 2012; 315:97-109. [PMID: 22981924 DOI: 10.1016/j.jtbi.2012.08.028] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2012] [Revised: 07/21/2012] [Accepted: 08/23/2012] [Indexed: 11/17/2022]
Abstract
Chemotherapy has profound effects on the hematopoietic system, most notably leading to neutropenia. Granulocyte colony stimulating factor (G-CSF) is often used to deal with this neutropenia, but the response is highly variable. In this paper we examine the role of pharmacokinetics and delivery protocols in shaping the neutrophil responses to chemotherapy and G-CSF. Neutrophil responses to different protocols of chemotherapy administration with varying dosages, infusion times, and schedules are studied through a mathematical model. We find that a single dose of chemotherapy produces a damped oscillation in neutrophil levels, and short-term applications of chemotherapy can induce permanent oscillations in neutrophil level if there is a bistability in the system. In addition, we confirm previous findings [Zhuge et al., J. Theor. Biol., 293(2012), 111-120] that when periodic chemotherapy is given, there is a significant period of delivery that induces resonance in the system and exacerbates the corresponding neutropenia. The width of this resonant period peak increases with the recovery rate after a single chemotherapy, which is given by the real part of the dominant eigenvalue pair at the steady state, and both are determined by a single cooperativity coefficient in the feedback function for the neutrophils. Our numerical studies show that the neutropenia caused by chemotherapy can be overcome if G-CSF is given early after chemotherapy but can actually be worsened if G-CSF is given later, consistent with results reported in Zhuge et al. (2012). The nadir in neutrophil level is found to be more sensitive to the dosage of chemotherapy than that of the G-CSF. Furthermore, dependence of our results with changes in key pharmacokinetic parameters as well as initial functions are studied. Thus, this study illuminates the potential for destructive resonance leading to neutropenia in response to periodic chemotherapy, and explores and explains why the timing of G-CSF is so crucial for successful reversal of chemotherapy induced neutropenia.
Collapse
Affiliation(s)
- Grace Brooks
- Department of Physiology, Centre for Applied Mathematics in Bioscience and Medicine, McGill University, Montreal, QC, Canada H4X 2C1
| | | | | | | |
Collapse
|
18
|
Zhuge C, Lei J, Mackey MC. Neutrophil dynamics in response to chemotherapy and G-CSF. J Theor Biol 2012; 293:111-20. [DOI: 10.1016/j.jtbi.2011.10.017] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2011] [Revised: 10/13/2011] [Accepted: 10/14/2011] [Indexed: 11/25/2022]
|
19
|
Nakaoka S, Aihara K. Mathematical study on kinetics of hematopoietic stem cells--theoretical conditions for successful transplantation. JOURNAL OF BIOLOGICAL DYNAMICS 2011; 6:836-854. [PMID: 22873618 DOI: 10.1080/17513758.2011.588343] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
Numerous haematological diseases occur due to dysfunctions during homeostasis processes of blood cell production. Haematopoietic stem cell transplantation (HSCT) is a therapeutic option for the treatment of haematological malignancy and congenital immunodeficiency. Today, HSCT is widely applied as an alternative method to bone marrow transplantation; however, HSCT can be a risky procedure because of potential side effects and complications after transplantations. Although an optimal regimen to achieve successful HSCT while maintaining quality of life is to be developed, even theoretical considerations such as the evaluations of successful engraftments and proposals of clinical management strategies have not been fully discussed yet. In this paper, we construct and investigate mathematical models that describe the kinetics of hematopoietic stem cell self-renewal and granulopoiesis under the influence of growth factors. Moreover, we derive theoretical conditions for successful HSCT, primarily on the basis of the idea that the basic reproduction number R (0) represents a threshold condition for a population to successfully grow in a given steady-state environment. Successful engraftment of transplanted haematopoietic stem cells (HSCs) is subsequently ensured by employing a concept of dynamical systems theory known as 'persistence'. On the basis of the implications from the modelling study, we discuss how the conditions derived for a successful HSCT are used to link to experimental studies.
Collapse
Affiliation(s)
- Shinji Nakaoka
- FIRST, Aihara Innovative Mathematical Modelling Project, Japan Science and Technology Agency, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8505, Japan.
| | | |
Collapse
|