Research Article  Open Access
M. Ivan Ariful Fathoni, Fajar AdiKusumo, Gunardi Gunardi, Susanna Hilda Hutajulu, "Dynamics of a Breast Cancer Model for Neutropenia Case due to Chemotherapy Effects", International Journal of Differential Equations, vol. 2021, Article ID 3401639, 8 pages, 2021. https://doi.org/10.1155/2021/3401639
Dynamics of a Breast Cancer Model for Neutropenia Case due to Chemotherapy Effects
Abstract
Breast cancer is a type of carcinoma with a high prevalence. The treatment of breast cancer through chemotherapy can cause a risk to healthy cells throughout the body. The neutrophil is one of the cells that is influenced by chemotherapy drugs. Chemotherapyinduced neutropenia is one of the most common toxic effects experienced by patients and often threatens chemotherapy to use efficiency. In this paper, we introduce an interaction model between blood components, i.e., neutrophil, lymphocytes, and albumin, with chemotherapy drugs. The model is important to understand the neutropenia effect due to chemotherapy in mathematical perspective and to calculate breast cancer patients’ survival level. Our model is a fourdimensional system of the firstorder ODE with 13dimensional parameter space. We focus our study for understanding the steadystate conditions and the bifurcations when the parameter values are varied. Here, we also study the role of albumin for reducing the neutropenia effects for breast cancer patients mathematically, where the results can be used as an alternative solution for treating neutropenia in a breast cancer case.
1. Introduction
Breast cancer is a type of a high prevalence carcinoma. It can occur in both men and women, where the higher prevalence is in women. Data from WHO in 2018 showed that the number of new breast cancer cases was 2,088,849 or 11.6% of the world’s total number of cancer cases. This number is only different from 5,027 cases of lung cancer, which occupies the first position. Based on this number, 626,679 of them were the death cases. Asia is in the top rank of the number of breast cancer incidence and mortality in the world [1]. The disease is treatable if it is detected earlier and is still in its earler stages. The higher the stage implies, the smaller the life expectancy is. There is no cancer that has a 100% mortality rate; it means that the treatment can be done to get a cure, especially if treatment is started earlier [2].
There are several studies being carried out in cancer modeling and treatment. Schättler et al. [3] developed a mathematical model with minimal parameters for lowdose metronomic chemotherapy. The model takes into account the angiogenic signaling between the tumor and the vasculature and the tumor inhibitory effects of tumorimmune system interactions. Liu and Yang [4] developed a cancer treatment model with radiotherapy followed by chemotherapy. Jordão and Tavares [5] studied the dynamics of healthy and diseased cells. The model for healthy cells is based on the p53 gene network, cell cycle control mechanism, the MAP kinase pathway, and the PI3K or Akt pathway, a cancer cell model embodied by deregulation in a healthy cell model. SolisPerez et al. in [6] created a mathematical model that considers the population dynamics among cancer stem cells, tumor cells, healthy cells, the effects of excess estrogen, and the body’s natural immune response to cell populations.
Several mathematical models that are more specific to breast cancer have been introduced before. Enderling et al. in [7] modeled the breast cancer development, the local treatment, and the recurrence. Furthermore, Prameswari and AdiKusumo [8] constructed a model of the phase transition regulation on cell cycle to characterize breast cancer cells’ growth in the earlier stage, which is indicated by the anomalies of cell proliferation. The model knows the behavior changes in cell concentration used as an early indicator of cancer growth. In the same year, Kusuma and AdiKusumo [9] introduced a mathematical model of DNA damage as the estrogen response involving the G1/S transition phase in the cell cycle. The model is a foundation to understand the initial behavior of breast cancer.
Treatment of breast cancer through chemotherapy can pose a risk to healthy cells throughout the body. The neutrophil is white blood cells that have a short life span (6–10 hours). They provide the first line of defense part of the body’s immune system. Neutrophils are the most common types of leukocytes found in human peripheral blood. The neutrophil concentration is about 60% of white blood cells in human blood. They are major cellular contributors to inflammation, mediating the initial phases of the inflammatory response. They can eliminate microbes through various mechanisms and are essential for the defense of the human body. Acquired neutrophil disorders are usually caused by cytotoxic chemotherapy or immunosuppressive agents. Primary neutrophil disorders usually result from the deactivation of individual gene mutations that result in impaired neutrophil number or function [10].
Neutropenia is defined as a decreasing number of neutrophils in the circulating blood. Chemotherapyinduced neutropenia is one of the most common toxic effects experienced by patients and often threatens the efficiency of chemotherapy use. Neutropenia usually occurs 7 to 10 days after chemotherapy. It can increase the risk of a complication of infection and fever known as febrile neutropenia (FN) [11]. In several papers, they consider the prediction of the occurrence of neutropenia in breast cancer patients. A mathematical model for optimizing neutropenia minimization was introduced in the study by Craig et al. [12]. A model that was developed is based on the characterization of neutrophil production in the bone marrow. The model was integrated with a pharmacokinetic and pharmacodynamic (PKPD) model of a chemotherapy agent. In another case with the model proposed by Craig et al. [12], we will also look at a model of neutrophil interaction with variables that, according to the data, affect the survival of breast cancer patients, namely, lymphocytes, chemotherapy drugs, and albumin.
Lymphocytes are white blood cells that play a role in maintaining the immune system. Lymphocytes are one of the main components of the immune system and are indispensable in the immune response to cancer. Lymphocytes are divided into two types, namely, T lymphocytes and B lymphocytes. T lymphocytes kill bacteria and viruses, while B lymphocytes function to make antibodies that fight antigens. The neutrophillymphocyte ratio (NLR) has been evaluated in various types of cancer. It is proposed as a marker of inflammation to predict outcomes. Increased NLR in the peripheral blood of patients with cancer could be a sign of a manifestation of systemic inflammation induced by cancer cells [14]. Hartono’s study [15] discussed the assessment of the neutrophillymphocyte and albumin ratio before and after chemotherapy. The results of the study demonstrated a reduction in NLR and an increase in albumin after chemotherapy.
Albumin is the most substantial protein in blood plasma. It is maintaining oncotic pressure, transporting fatty acids to the liver, transporting drugs, shortening the halflife of the drug, transporting bilirubin, binding + ions a buffer solution, transporting thyroid hormones, and transporting other hormones, especially those that can dissolve in fat. Albumin also functions as a transport for various substances, including bilirubin, fatty acids, metals, ions, hormones, and drugs. It plays a vital role in the regeneration and repair process of cells, as an agent that forms cell bonds. Albumin also helps to give signal for the immune system about signs of damage to body cells.
Chemotherapy causes hematopoietic progenitor cells in the bone marrow to experience depletion and acute myelosuppression [16]. In conditions of severe hemopoiesis stress, stem cells can repair themselves, proliferate, and differentiate into various blood cell series in the periphery. However, suppose the selfrenewal ability of hemopoiesis stem cells is impaired. In that case, there will be permanent damage to the hematopoietic system resulting in bone marrow failure. It can occur when the bone marrow is exposed to high doses of chemotherapy agents or radiation [17]. The effect of neutropenia appears to be higher than that of cytopenia of the erythroid and platelet series because the composition of neutrophil precursors occupies the most significant number in the bone marrow, which is more than 50% [18].
In this paper, we consider a new mathematical model on breast cancer therapy that includes the neutropenia effects due to the chemotherapy. Our model is a 4dimensional system of firstorder ODE based on the interaction between the concentrations of neutrophils (N), chemotherapy drugs (C), lymphocytes (L), and albumin (A). We focus our analysis to understand the role of parameters to the steadystate conditions of the system. In this case, we employ bifurcation analysis to study the existence and stability of the equilibria. It is important for calculating the survival level of the breast cancer patients.
2. Model Formulation
The interaction between the neutrophils (N), chemotherapy drugs (C), lymphocytes (L), and albumin (A) is shown in Figure 1. The interaction between those substances is created by the results of some medical papers (see [10, 15–18]), and the correlation analysis of the data from National Health and Nutrition Examination Survey (see [14]).
Our model is formulated as follows:
By the transfer diagram shown in Figure 1, we construct a mathematical model that shows the interactions between the neutrophil, lymphocyte, and albumin concentrations due to the appearance of chemotherapy drugs concentration (see system (1)). System (1) is a four dimensional system of the firstorder ODE in indimensional parameter space. The parameters used in the model are assumed to be positive, with the information in Table 1.

3. Existence of Equilibrium Solutions
System (1) has eight equilibrium points, but we consider the analysis of the first seven equilibrium points. Calculating the last equilibrium point is very complicated, so we left it as an open problem in this paper. The existences of those equilibrium points are presented in Theorem 1 and Theorem 2.
Theorem 1 (the existences of equilibrium points to of System (1)). (i)The equilibrium point exists if and (ii)The equilibrium points and exist for all parameter values(iii)The equilibrium point exists if it satisfies the conditions in Table 2(iv)The equilibrium point exists if (v)The equilibrium point exists if

Theorem 2 (the existences of equilibrium points of System (1)). The equilibrium point with and where , , , and has two possible existences due to value . If value is negative, then condition of existence is . If value is positive, then condition of existence is .
Then, the results of Theorem 1 are summarized in Table 3 in the appendix. Based on the results of the theorem, we show the area in parameter spaces that two or more equilibria coexist in Table 4.

 
means that the condition for the existence of is number 20 in Table 2. 
4. Stability Analysis
We apply the linearization method to determine the stability of the equilibrium point of System (1), and we have the Jacobian matrix as follows:where , , , and
The eigenvalues can be obtained from the characteristic equation evaluated at each point of equilibrium. The stability analysis discussed in this study is only at the equilibrium point to . The stability analysis of the equilibrium point will be carried out by considering the numerical simulation in special cases. The stability analysis of the equilibrium point T4 will be carried out by considering the numerical simulation in special cases.
4.1. Equilibrium Point
Linearization of System (1) near the equilibrium produces the characteristic equation as follows:where and
From (3), it can be seen that there is a positive eigenvalue, i.e., . It can be concluded that the equilibrium point is unstable.
4.2. Equilibrium Point
The characteristic equation of the linearization of System (1) near the equilibrium point is , so that its eigenvalues are and In this case, we have a positive definite eigenvalue, so the equilibrium point of is unstable.
4.3. Equilibrium Point
The characteristic equation of the linearization of System (1) near the equilibrium point is where the eigenvalues are , and The positive eigenvalues indicate that the equilibrium point of is unstable.
4.4. Equilibrium Point
The Jacobian matrix obtained from the linearization of System (1) near the equilibrium point iswhere , . It is assumed that the values of , and are positive. The characteristic equation formed iswhere and
which has the eigenvalues where and . Based on the equations of , and , the stability conditions for the equilibrium point of are obtained, namely, , and .
4.5. Equilibrium Point
Linearization of System (1) near the equilibrium point obtains the characteristic equation as follows:where , all of which are assumed to be positive. The terms of the stability of this point are analyzed using the Routh–Hurwitz criteria, and we obtain
4.6. Equilibrium Point
Linearization of System (1) near the equilibrium produces the characteristic equation as follows:where , all of which are assumed to be positive. The terms of the stability of this point are analyzed using the Routh–Hurwitz criteria with the result as follows:
The stability of the equilibrium points and is very important for neutropenia. If those equilibria are stable, the neutrophils will decrease and go to zero. This situation triggers the neutropenia case due to the chemotherapy. The other important equilibrium point is . The stability of this equilibrium shows the reduction of lymphocyte which is one of the main component in the immune system against the cancer.
5. Numerical Simulations
Before treatment, all breast cancer patients will perform a complete blood laboratory examination, including examining the concentration of neutrophils, lymphocytes, and albumin. If the patient is indicated that there is breast cancer, then the patient will be treated medically through chemotherapy. However, chemotherapy drugs also consider the side effects of neutropenia. Mathematical models are needed to understand and to anticipate the side effects of neutropenia. Simulations of the changes in neutrophils, chemotherapy drugs, lymphocytes, and albumin concentration during chemotherapy are performed to predict neutropenia.
The simulation was carried out based on the assumption that the concentration of neutrophils was 0.8 (80%) and that of lymphocytes was 0.2 (20%) from white blood cells, and that of albumin was 0.7 (70%) from blood plasma. Chemotherapy concentration is assumed to be 0.4 (40%) from the maximum limit capacity of chemotherapy drugs in the body. Medically, neutrophils, albumin, and neutrophils in human blood circulation cannot be zero, so the first case used in this simulation is directed towards the equilibrium point of . By using the parameter values , we obtained dynamic changes in the model as shown in the curves of Figures 2 and 3. Numerical simulations are also used to verify the results of the dynamic analysis that has been carried out. Based on the predetermined parameter values, using four different initial conditions, the solution will go to the equilibrium point , as shown in Figure 4.
Figures 2 and 3 show the dynamic changes of the four observed variables, where neutrophils at the beginning of chemotherapy increased until the second day and then decreased significantly to equilibrium at 0.725. It is consistent with the medical facts of the side effects of neutropenia after chemotherapy. The chemotherapy concentration at the beginning of chemotherapy is 0.4, gradually decreasing until it reaches zero. Albumin plays a role in forming white blood cells; in this case, the albumin has increased to 0.96 before dropping to 0.8 on day 25. The increase in albumin correlates with an increase in lymphocytes. Lymphocytes that play a role against cancer cells and chemotherapy have increased to 0.456. Based on these dynamic changes, it is known that the neutrophillymphocyte ratio has decreased, where at the beginning of chemotherapy, it has a ratio of 4. After chemotherapy, it drops to 1.6.
The conditions in Figures 2 and 3 show the occurrence of neutropenia, where the neutrophil concentration is significantly reduced. One alternative to prevent neutropenia is to reduce albumin. In the second case, we take the same parameter value as first case for , and , by changing value. Albumin can be reduced by reducing its maximum capacity () to 0.68 so that the albumin concentrates towards the steadystate at this value. In this case, the neutrophils raised at the beginning of chemotherapy will not decrease significantly, and it reaches the steadystate, which is at 0.8. The lymphocytes continued to increase towards the final state at 0.31. The neutrophillymphocyte ratio in this second case was reduced to 2.58 after chemotherapy. This second case simulation is shown in the graph in Figure 5.
In this section, we can find dynamic simulations of the interaction of blood components, represented by neutrophils, lymphocytes, and albumin, with the effects of chemotherapy drugs. Based on a mathematical model simulation, it can be seen that the regulation of albumin concentration in blood has a role in preventing and predicting the side effects of neutropenia due to chemotherapy drugs. So, if the patient’s condition requires chemotherapy treatment, medical experts can reduce the level of albumin concentration in the blood. It is an alternative solution to prevent neutropenia based on a mathematical approach.
6. Concluding Remarks
Breast cancer is a dangerous disease with the side effects of treatment. A side effect of chemotherapy for breast cancer that needs to be watched out for is neutropenia, where the neutrophils in the white blood cells are reduced to being abnormal. The dynamics of the three components of blood cells with chemotherapy concentrations have been constructed using a differential equation system (1) based on data and medical facts described previously.
Simulations are carried out to determine changes in the number of neutrophils, chemotherapy, lymphocytes, and albumin concentrations from time to time if given the initial conditions. By knowing the dynamics that will occur, we can control changes in the concentration of the three components of blood cells to return to normal after chemotherapy. This study formed a dynamic model of blood components and chemotherapy based on hospital data and medical literature.
The mathematical model constructed in this study can describe the dynamics of changes in the concentration of neutrophils, chemotherapy, lymphocytes, and albumin in the patient’s body to predict the possibility of neutropenia. One alternative in preventing neutropenia based on the model in this study is to reduce the albumin concentration. The results of this study can be used as an alternative solution for treating neutropenia in breast cancer patients, in mathematical point of view.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors thank the Deputy for Research and Strengthening Development of the Ministry of Research and Technology Indonesia (the National Research and Innovation Agency) which has provided research funding through Penelitian Dasar Unggulan Perguruan Tinggi (PDUPT) with contract agreement number 1683/UN1/DITLIT/DITLIT/PT/2021. The authors also thank the Department of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Gadjah Mada for the support in research facilities.
References
 Globocan, 2018, Cancer today: international agency for research on cancer, https://gco.iarc.fr/today/data/factsheets/cancers/20Breastfactsheet.pdf.
 M. C. Oz and M. F. Roizen, You: The Owner’s Manual: An Insider’s Guide to the Body that Will Make You Healthier and Younger, Harper Collins, New York, NY, USA, 2013.
 H. Schättler, U. Ledzewicz, and B. Amini, “Dynamical properties of A minimally parameterized mathematical model for metronomic chemotherapy,” Journal of Mathematical Biology, vol. 72, no. 5, pp. 1255–1280, 2016. View at: Publisher Site  Google Scholar
 Z. Liu and C. Yang, “A mathematical model of cancer treatment by radiotherapy followed by chemotherapy,” Mathematics and Computers in Simulation, vol. 124, pp. 1–15, 2016. View at: Publisher Site  Google Scholar
 G. Jordão and J. N. Tavares, “Mathematical models in cancer therapy,” Biosystems, vol. 162, pp. 12–23, 2017. View at: Publisher Site  Google Scholar
 J. E. SolísPérez, J. F. GómezAguilar, and A. Atangana, “A fractional mathematical model of breast cancer competition model,” Chaos, Solitons & Fractals, vol. 127, pp. 38–54, 2019. View at: Publisher Site  Google Scholar
 H. Enderling, M. A. J. Chaplain, A. R. A. Anderson, and J. S. Vaidya, “A mathematical model of breast cancer development, local treatment and recurrence,” Journal of Theoretical Biology, vol. 246, no. 2, pp. 245–259, 2007. View at: Publisher Site  Google Scholar
 K. Prameswari and F. AdiKusumo, “A mathematical model for phase transition regulation on the cell cycle in early stage of breast cancer,” in AIP Conference Proceedings, vol. 2192, no. 1, AIP Publishing LLC, 2019. View at: Google Scholar
 M. F. Kusuma and F. AdiKusumo, “A mathematical modelling for estradiol influence on DNA damage response and G1/S transition phase regulations in early stage of breast cancer,” in AIP Conference Proceedings, vol. 2192, no. 1, AIP Publishing LLC, 2019. View at: Google Scholar
 R. Suzuki, X. Wei, P. K. Allen, J. D. Cox, R. Komaki, and S. H. Lin, “Prognostic significance of total lymphocyte count, neutrophiltolymphocyte ratio, and platelettolymphocyte ratio in limitedstage smallcell lung cancer,” Clinical Lung Cancer, vol. 20, no. 2, pp. 117–123, 2018. View at: Google Scholar
 E. Å. Lundqvist, “Principles of chemotherapy,” International Journal of Gynecology & Obstetrics, vol. 119, Article ID 119S2, pp. S151–S154, 2012. View at: Publisher Site  Google Scholar
 M. Craig, A. R. Humphries, F. Nekka, J. Bélair, J. Li, and M. C. Mackey, “Neutrophil dynamics during concurrent chemotherapy and GCSF administration: mathematical modelling guides dose optimisation to minimise neutropenia,” Journal of Theoretical Biology, vol. 385, pp. 77–89, 2015. View at: Publisher Site  Google Scholar
 H. K. Lehman and B. H. Segal, “The role of neutrophils in host defense and disease,” The Journal of Allergy and Clinical Immunology, vol. 145, no. 6, pp. 1535–1544, 2020. View at: Publisher Site  Google Scholar
 National Health and Nutrition Examination Survey, Complete Blood Count with 5part Differential—Whole Blood (CBC_G), National Bureau of Economic Research, Cambridge, MA, USA, 2014, http://data.nber.org/nhanes/20112012/CBC_G.htm.
 B. Hartono, V. S. Pontoh, and M. A. Merung, “Penilaian jumlah neutrofil, limfosit dan trombosit, kadar protein reaktif C, kadar albumin, rasio neutrofil limfosit, serta rasio trombosit limfosit sebelum dan setelah terapi pada penderita karsinoma payudara,” Jurnal Biomedik: JBM, vol. 7, 3 pages, 2015. View at: Publisher Site  Google Scholar
 V. T. Devita Jr, S. Hellman, and S. A. Rosenberg, Principles and Practise of Oncology, Lippincot Williams and Wilkins, Philadelphia, PA, USA, 5th edition, 2008.
 K. R. Georgiou, B. K. Foster, and C. J. Xian, “Damage and recovery of the bone marrow microenvironment induced by cancer chemotherapy  potential regulatory role of chemokine CXCL12/receptor CXCR4 signalling,” Current Molecular Medicine, vol. 10, no. 5, pp. 440–453, 2010. View at: Publisher Site  Google Scholar
 F. Crea, E. Giovannetti, P. L. Zinzani, and R. Danesi, “Pharmacologic rationale for early GCSF prophylaxis in cancer patients and role of pharmacogenetics in treatment optimization,” Critical Reviews in Oncology, vol. 72, no. 1, pp. 21–44, 2009. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2021 M. Ivan Ariful Fathoni et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.