ANNEX PUBLISHERS

Journal of Computational Systems Biology

ISSN: 2455-7625

Open Access
Case report
Max Screen >>

Development of the Analytical Model for the Assessment of the Efficiencies of Different Therapeutic Modalities in Leukemia

Received Date: December 29, 2014 Accepted Date: August 03, 2015 Published Date: August 07, 2015

Copyright: © 2015 Majumder D. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Related article at Pubmed, Google Scholar

Abstract

Different therapeutic regimes like chemotherapy, cytokine based immunotherapy and/or stem cell transplantation are generally be the suggested options for treatment of leukemia. However, the selection of the best therapeutic strategy for an individual patient remains uncertain till now. Several analytical models are proposed for each of the therapeutic strategies; however, no single analytical model is available which can make a comparative assessment regarding the long-term therapeutic efficacy among these strategies. This, in turn, may limit the clinical application of such analytical models. To address this issue we have developed an analytical model that can make a comparative assessment among these therapeutic regimes along with their different combinations. The model also has the flexibility to test their efficacies in different types of leukemia. Contrary to the previous models here we have incorporated several other important constraints like chemotherapeutic toxicity and chronic effects of HLA incompatibility that may limit therapeutic success. The simulation studies with the developed model indicate towards the immune potentiation as the determining factor to bring the leukemia free survival. It also reveals that the time optimization for immune boosting is another important crucial factor for providing leukemia free survival.

Keywords: Leukemia; Difference delay equation; Cachexia; Chemotherapy; HLA matching

Introduction

Several mathematical models have been proposed to study the complex dynamical behavior of the hematopoiesis. These models have attempted to reveal the pathogenesis of different hematological diseases (HD) including hematological malignancies (HM). Amongst, the first work was targeted with acute myelogenous leukemias (AML) where systems of delay ordinary difference equations (DODE) are considered for modeling the hematopoietic dynamics in different compartments of hematopoietic lineages and the effect of chemotherapeutic intervention [1-3]. This modeling strategy has also been used in explaining the behaviour of periodic chronic myelogeneous leukemia and cyclical neutropenia diseases [4,5]. Similar modeling strategy was further extended for designing of long term targeted chemotherapeutic strategies (imatinib) in chronic myelogenous leukemias (CML) [6] and the effect of autologous vaccination procedure by T-cells [7]. Using a system of nonlinear delay differential equations, further work suggests that stability of the system is sensitive to the T cell related parameters like T cell stimulation and/or supply of T cell into the system [8].

However, models without delay have also been proposed by several workers [9-11]. These sorts of modeling approaches are helpful for the control and stability analysis of the system. Presently available analytical models for hematopoiesis consisting with multicompartmental models and the incorporation of ordinary differential delay equations makes it feasible to study the replication potential of the hematopoietic cells in several compartments as well as the chemotherapeutic efficiency, whereas other models without delay consideration are formulated for the simplicity of stability and control analysis.

For achieving better tolerance of the hematopoietic cell transplantation, hematopoietic stem cell (HSC) transplantation is suggested. Stem cell transplantation (SCT), an emerging therapeutic procedure in hematology involves a complex procedure,however, little has been done with the mathematical modeling, in spite of rich medical and clinical literatures are available in this field. In this regard works of DeConde et al, 2005 [12] can be mentioned where delay differential equations are used to explain the time evolution of six populations after allo-SCT in chronic myeloid leukemia and concluded that remission depends on the allo-reactive T cell immune function [7]. Previous modeling approaches were further modified to incorporate the drug resistant mutation in CML and suggest for combination therapy as a means of long-term remission in CML undergoing imatinib treatment [13]. Further works suggests that relapse is due to presence of dormant stem cells and/or symmetric division of stem cell acquire for the development of drug resistant phenotype [14,15].

However, allo-SCT transplantation is a complex therapeutic procedure as it may require intermittent intervention at different time points due to non-availability of ideal HLA matching in clinical scenarios. Hence, a flexible model with different input variables would be the ideal for the assessment of the potentiality of different therapeutic modalities for varying circumstances. In several instances of leukemia, partial remission is achieved by the application of myeloablative chemotherapy (MAC) and/or radiotherapy.Allogenic hematopoietic stem cell transplantation (allo-HSCT) is therefore suggested for achieving long-term leukemia free survival (LFS). Donor lymphocyte alloreactivity against the host hematopoietic system translates powerful graft-versus-leukemia (GVL) reactions [16]. However if the transplanted cells have the HLA incompatibility with the host then a severe graft-versushost diseases (GVHD) occurs which favors the opportunistic infection within the host. This has the significance of transplanted related mortality. Allo-HSC transplantation is regarded as the double edge sword as it protects the host from leukemia relapse and simultaneously the death of the host.

Though it is difficult to get, however complete remission through transplantation is possible with the ideal HLA matching with the syngenic (monogygotic twin of the patient) donor cells. Allo-HSC transplantation has two folds benefits – on one hand the conditioning regimen consisting of MAC destroys the hematopoietic process of patient which in turn reduce the tumor load in the patient, and following procedure of allo-HSC transplantation exerts a long lasting immunologic GVL by the donor derived T-lymphocytes which helps to remove the residual diseases of the patient after the MAC [17].

Now-a-days, allo-HSC transplantation (HSCT) is done with three sources – bone marrow (BM), peripheral blood (PBL) and umbilical cord blood (UCB). In recent time majority of cases HSCT is performed with PBL as the sources of HSC. For harvesting of HSC from the PBL it is needed to mobilize the HSC to the peripheral blood of the donor by the application of GMCSF or low dose cyclophosphamide. Allo-HSC of UCB origin is the better option compared to the adult HSC because of its greater tolerance within the host body due to low expression of the HLA on the stem cell surface [18].

Attempt has been made to study the post transplant immune responses in CML [2]; however, the management of the post transplant complications is still a challenge for different types of leukemias and yet to be determined. Though cancer vaccination strategy has been designed [7], however, in designing this model several parametric estimations like kinetic coefficient, probability of cancer cell death is not a very clinically feasible approach while one considers the individual leukemia cases. Therefore it is needed to develop a generalized model for different leukemias that could be more clinically compatible in the assessment of different therapeutic procedures including transplantation and its complications. Here we propose a simplified but clinically compatible analytical model for the assessment of the outcome of chemotherapeutic strategies, immune potentiation, HSCT and post-HSCT complications.

Methods

In present model we have considered a difference delay equation based three compartmental system model [19] having 7 cell types - the stem cell compartment having hematopoietic stem cell (S), progenitor compartment having three progenitor cells of three lineages namely erythroid (P1), leukocytoid (P2) and megakaryocytoid (P3) and matured cell compartment having mature cells of three lineages namely RBC (B1), WBC (B2) and platelet (THB) (B3). In each compartment three different category cells – normal cells (g), drug sensitive cells (s) and drug resistive cells (r) are considered. Hence in first compartment three cell types (Sg, Ss and Sr), in second compartment total nine cell types (P1g, P2g, P3g, P1s, P2s, P3s, P1r, P2r and P3r) and in third compartment total nine cell types (B1g, B2g, B3g, B1s, B2s, B3s, B1r, B2r and B3r) are considered. Therefore in present system model we have considered total 21 types of cells (Figure 1).

The equation (1) is used as a general form to represent the cellular dynamics of each cell type of any three compartment of the hematopoietic system at discrete time interval (k).

In the above equations (1) – (3) Nx(k) represents the number of cell on k-th day. N represents the cell type and suffix x represents the category of that concerned cell type. Nxp(k-dk) is the precursor cell type from which the present cell type was differentiated dk time ago. Hence dk is the delay time required for cell maturation by the action of different cytokines, cellular signaling molecules and transcription factors. Each cell type grows exponentially with its own multiplication rate (mNx) and decays with its own apoptosis rate (aNx). Thus, for normal stem cell these variables are denoted with (mSg) and (aSg). Similarly for drug sensitive P1 type cell, these variables are represented by (mP1s) and (aP1s). These notational schemes are followed for other cell types and are represented in Table 1A-C. It is assumed that multiplication rate of each cell type at the matured cells compartment is zero.

For the cells which are differentiable such as stem cell (S), erythroid (P1), leukocytoid (P2) and megakaryocytoid (P3), they will have c1 = 1 in their system equation. For the system equation of RBC, WBC and platelet, c1 = 0 as they do not differentiate. When the effect of differentiation is absent, that is, when k is less than the delay time (dk) required for the formation of the cell type to the down-stream cell type of the same lineage, c2 = 0. Otherwise c2 will be a positive quantity. In the present model, during the stage of differentiation from stem cell to progenitor cell stage, value of c2 = 1/3 (as it is assumed that stem cell is distributed equally to three lineages after differentiation). However, in the system equation for the cells of the subsequent compartment, c2 =1 (as each progenitor cell type is differentiated into a single lineage only).

Again Nxdr is the differentiation rate of any cell type and Nxpdr is the differentiation rate of its concerned precursor cell type. For example, in the system equation of stem cells has only Nxdr and for sensitive category stem cell, Nxdr is represented as Ssdr. Similarly, for drug resistive category of P1 lineage cell type, Nxdr = P1rdr and Nxpdr = Sxdr; for P3 lineage cell type of drug sensitive category, Nxdr = P3sdr and Nxpdr = Ssdr and differentiation rates for other cell types are listed in Table 1A-C.

Now, if k<dk (delay time required for the formation of a progenitor cell from its previous compartment cell of the same lineage) the system equation of concerned differentiated cell types can be represented through the following generalized equation:

In equations (3) - (5), the value of f is different for different cell types. The above mentioned generalized equations may be used to develop the system model in form of a difference delay equation [19] which is mentioned below.

While considering the conversion from concerned cell type (i.e., Nx) to other two cell categories of same type concerned conversion rates represented by CN xy and CN xz are subtracted from the diagonal elements of A matrix. Again a(i-2,j-1), a(i-2,j), a(i-1,j-2), a(i-1,j), a(i,j-2) and a(i,j-1) elements in matrix A where i = j = 3, 6, 9, 12, 15, 18, 21 can be modified with the corresponding conversion rates to incorporate conversion from other two categories to concerned cell type. For example, putting i = j = 3 in (i,j), we get six matrix elements of A namely a(1,2), a(1,3), a(2,1), a(2,3), a(3,1), a(3,2); where a(1,2) and a(1,3) are the conversion rates from Sr and Ss to Sg i.e., CSsg and CSrg respectively (Table 1A). For other category cells of stem cell type the other elements become operative similarly; for P1, i = j = 6 and likewise.

Physiological feed-back control system

The number of different normal cells is maintained through physiological homeostatic mechanism. In the mature cell compartment if the number is decreased below a certain level then positive feedback system causes the respective progenitor cells to increase its number. Similarly, if the number increases, then it produces negative feedback. Functionally within the physiological system, if RBC counts decreases, then reticulocytes (autocrine feedback) and erythroblasts are stimulated by erythropoietin (EPO) to form mature RBC. Similarly, when platelet counts decreases, the normal count (number) is maintained by megakaryocytes and megakaryoblast. Likewise when WBC cell count (say T-lymphocyte) decreases, its' number is increased by IL-2 or other cytokines by stimulating the corresponding precursor cells or lymphoblasts. The negative feedback occurs within the leukocyte lineage when the number increases. For lymphoid lineage a feedback system is also considered from mature cell level upto the HSC level. Autocrine regulation is considered for blast cell level and HSC by the effect of cytokine; for example, regulation by EPO, IL-2, thrombopoietin and stem cell factor at erythroblast, lymphoblast, megakaryoblast and HSC level respectively. In the model, both autocrine and inter-compartmental feedback mechanisms are considered. Inter-compartmental feedback mechanism is being effective by modulating the differentiation rate and apoptosis rate (Figure 2).

Host immune system

It is considered that under normal condition, host immune system is capable of removing leukemic cells [20-22]. This has been incorporated by modifying equation (7a) into equation (7b):

where the functional effective (immunity related) killing of the malignant cells (killeff) by B2g becomes effective depending upon the minimum cell number of B2g (Min B2g) and is operated by killim and CF. The assumed efficiency of immune related killing (killeff) by the host immune cells are in the order of B2s>B2r=P2s>P2r>Ss>Sr. It is assumed that initially there is total (functional) immune suppression i.e., killim = 0. Application of any immune-modulatory agent is required to make the immune related killing effective. It is also assumed that CF is a time varying variable that saturates up to 60% immunity level of the normal population after application of immune-modulatory agent for a certain period of time and reaching this level will be dependent upon the applied dose of that immune-modulatory agent.

Cancer cachexia

About 50% of the cancer patients suffer from a progressive weight loss due to atropy of adipose tissue and skeletal muscle. This phenomenon is called cachexia. This is more often associated with preterminal cancer patients. Cachexia may be present in the early stages of tumor growth before any signs or symptoms of malignancy [23]. Recent time several tumor related soluble factors are identified those are responsible for cancer cachexia [24,25]. Though cancer cachexia is not a term frequently used in leukemia; however in recent time by ultrastructural analysis it has been shown that there is some morphological de-arrangement in the red cell membrane (phenotype effect) in different leukemias including chronic and myelodysplastic syndrome [26]. Such phenomenon is known as cancer cachexia and is considered in the model.

Leukmic blast cells may produce much more cachexia effect. It is considered that resistant variety of malignant cells will produce the most detrimental effect compared to other cell types. Two types of toxicity are considered here – phenotypic and genotypic. The gradations of phenotype effect are implemented by three factors MTox1, MTox2 and MTox3 and genotype effect implemented by factor Tox. Considered gradation of toxicity intensities are in the following order: MTox1 > MTox2 > MTox3 > Tox (Figure 3). Hence a(i,j) of equation (7b) will be modified as:

Again these three factors are multiplied with three different factors to implement the severity of killing intensity in different compartments. It has been assumed that the order of killing intensity is decreasing from matured cell compartment to stem cell level. Hence MTox1, MTox2 and MTox3 are again multiplied with another three factors ch_MS (= 0.02) for stem cell, ch_MP (= 1.5) for progenitor cell compartment and ch_MB (= 5) for matured cell compartment. It is also assumed that ch_MP (= 0) when P2r (k) + P2s (k) <10,000.

Myeloablative chemotherapy model

Conventionally, in clinical practice patients are treated with a myeloablative (MYL) chemotherapeutic regime. This has been incorporated in the model (Figure 4). However, before applying the drug strategy, cells are allowed to grow for certain period of time (say, 30 days). On the consecutive day myeloablative chemotherapy is applied. MYL is introduced to the system model as a specific subtractive term with a fixed percentage to the existing number of a specific cell population at discrete interval of time. The percentage of cell killing is dependent on the drug sensitivity of a particular cell type.

In the system model there is a provision of changing the drug dose, the number of cycles of drug application and interval between two successive drug applications. So in the system equations a(i,j) of equation (7c) have been further modified as:

where, dmyl(k) is the amount of drug present within the system denoted either the amount of applied drug (drugmyldose) on the day of drug application or drug retention (drugmylret) on the subsequent days as a fraction of the previous day's applied drug (it is assumed that in each corresponding day of drug application a certain amount of drug is cleared from the system);

Drug application is started on the day kmylsd and nmyl is the number of cycles and tmyld is the interval of two subsequent drug applications in days and Cmyl is the ON/OFF switch for drug application. MylNx is the myeloablative drug sensitivity of the concerned cell Nx . The drug sensitivity MylNx can be varied accordingly from different lineages and/or cell types (Table 2A).

Model for the supportive therapy to MYL: High dose of MYL therapy produces reduction in RBC and platelet level in the patient. In clinical practice, this has been tackled through intermittent transfusion of exogenous RBC and platelet. These sorts of therapeutic management are also incorporated in the model. So equation (7) is modified as

On the day of MYL application, minimal level (a set value) of RBC (RBCllm) and platelets (Plateletllm) are checked and if found lower than that level, then RBC and/or platelet are transfused. If both the values are sufficient then drug will be delivered. If either value is insufficient, then, as a precautionary measure, an upper lower level of the other (that is, if RBCllm is less then Plateletulm or if Plateletllm is less then RBCulm) will be checked and if both are found below, then transfusion decision of both cell types will be made accordingly. Drug is not delivered on the day of any type of transfusion is applied. And in that case MYL drug application will be delayed by Sft days from kmylsd and after which subsequent days will be applied. This strategy prevents the unnecessary delay in subsequent drug application. The small delay period after transfusion is considered due to restoration of patient's physiological condition for tolerating the subsequent chemotherapy load; however, model has the flexibility to change or avoid this gap period by making Sft=0. Again, transfused cells are represented through following relations -

TRBC and TPLATELET represent amount of transfused RBC and platelets respectively. E1 and E2 may be fraction or whole number. MylTPlatelet and MylTRBC are the MYL drug sensitivities of TRBC and TPLATELET respectively. Application strategy of blood transfusion is demonstrated in Figure 5.

Hematopoietic Stem Cell (HSC) transplantation model

With a mismatched HLA between donor and host, the transplanted donor cells will also be recognized as foreign cell by the host immune system and vice versa. This produces Graft versus host disease (GVHD) (Figure 6). It is to be noted here that generally MYL therapy (as mentioned earlier) is applied as the conditioning regime before stem cell transplantation.

Transplanted cells' behavior within recipient's system: Transplanted cells' behavior can be represented by modifying the equation (8). In equation (8) each of the matrixes will be updated as follows. Hence,

where a(i,j), b1(i,j), b2(i,j), b3(i,j), b4(i,j) representing the elements of [28×28] matrixes A, B1, B2, B3, B4 respectively and RBCtransfuse(i,j), PLATELETtransfuse(i,j) representing the elements of [28×1] matrixes RBCtransfuse and PLATELETtransfuse respectively. In equation (8a) the effect of cachexia on the donor cells are also included. In equation (8a), the effect of MYL drug on the donor cell is ignored as at the time of transplantation the quantity of residual drug in the system has been considered to be negligible. In HLA mismatched case, transplanted donor cells are killed due to GVHD reaction is incorporated into the system model by a factor killb2g→D. Transplantation with 100% HLA matching G = 0; otherwise G>0, that means with the increase in HLA mismatch level the value of G (Graft versus host diseases, GVHD) will increase. The killing intensity killmul is operative for GVHD effect that can be different for different cell types of donor.

Host cells' behavior towards transplanted cells: Model has been formulated with an assumption that donor stem cell (Se) transplantation is made when chemotherapeutic drug level within the system is minimum say 10% of the last administered drug dose. However, this can be changed according to the choice of the investigator. Alloreactive lymphocytes (B2e) are developed [along with other lineages i.e., erythrocytes (B1e) and platelets (B3e)] from transplanted stem cells lead to an immune response causing Graft-versus-tumor (GVT) (or GVL) [killb2e→hm] and/or GVHD [killb2e→hg] effect. So a(i,j) represented in equation (8a) is modified further as:

In equation (8a-1) hm stands for host leukemic cells (Ss, Sr, P1s, P1r, P2s, P2r, P3s, P3r, B1s, B1r, B2s, B2r, B3s, B3r) and hg stands host normal cells (i.e., Sg, P1g, P2g, P3g, B1g, B2g and B3g). The GVHD reaction is operated into the system by modifying the kill factors killb2g→D and killb2e→hg. The killing intensity killmul is operative for GVHD and GVL effect that can be different for different cell types of host.

Nonlinear effect of GVHD: It is assumed that with the increase in the degree of HLA mismatch, effect of GVHD also increases simultaneously. The effect of GVHD (killb2e→hg) works on all the normal cells of the recipient. GVHD produces several aberrant unknown but cascade of biochemical reactions within the recipient that gradually increases with the progress of time.

It is assumed that the effect of GVHD will be operative when minimum cell numbers of B2e is developed from transplanted HSC. As a result with the progress of time GVHD may be operative nonlinearly depending on the cell numbers of B2e. Killing efficiency of B2e to other hematopoietic cells of donor will be governed with a multiplying factor effectprocell which is governed though the following conditions.

In the above relation, effect_p1g, effect_p2g and effect_p3g have some constant values (Table 4). The equation (8a-1) will be modified further as follows:

In HLA mismatch the GVHD related killing intensity by host cells (killmu1) can be different for different cell types of donor cells. In the model it is assumed that host lymphocytes are unable to kill the malignant cells and at the time of transplantation the effects of cancer cachexia on donor as well as host stem cells are nil, as donor cells are transplanted in a condition when the tumor load is minimum. Different kill factors that have been considered in GVHD are represented in Table 4.

In our model we have considered the effect of chronic GVHD (cGVHD). Clinically the effect of cGVHD is assessed by different qualitative based clinical observations of changes in liver, gut, lungs [27]. Mathematically this can be represented through a marker (MR) that has been considered in the model as a time varying squared function of the number of leukocytes of donor origin. This function will vary according to the effect of different grades of mismatch (G) or GVHD. This is represented by the following relation, where CalGVHD is a time varying calibration factor.

Cytokine therapy model

Due to immune suppression leukemic patients are susceptible to infections [28]. Such susceptibility becomes more pronounced after MYL therapy. Clinically this could be circumvented by the application of cytokines (CYT) like interferon, GM-CSF and/ or IL-2 [29-33]. Application of cytokines can also augment the number of immunocytes. This, in turn, may also kill the residual malignant cells. Such therapeutic strategy has incorporated in the model equation. This cytokine application scheme can also be operated single or in combination of other therapeutic strategy. Application strategy of MYL/Chemotherapy with transfusion followed by cytokine therapy is shown in Figure 7.

It is assumed that these cytokines will affect the multiplication rate of the progenitor cells of the leukocytic lineages, which in turn, also affect the differentiation process for the formation of mature immunocytes. To introduce this sorts of therapeutic procedure, the a(i,j) of equation (8a-2) will be modified as follows:

kCYTsd, nCYT, tCYTd, mCYT and CCYT indicate the starting day of cytokine (drug) application, the number of cycles and the interval of two subsequent cytokine drug applications (in days), cytokine induced effect (additive) on the multiplication rate of recipient's target cells and the ON/OFF switch for cytokine (drug) application. This modeling scheme can also be operative after HSC transplantation, if needed and then application of cytokine will affect the transplanted cells of donor origin through mDCYT.

Maintenance of transplantation by immunosuppressive drug

To keep the chronic GVHD under control immunosuppressive drugs (dim) are generally suggested. Present model has the provision in the assessment of such therapeutic management. For modeling it has been considered that immunosuppressive drugs affect different leukocytic subpopulations of both host and transplanted cells i.e., P2g, P2r, P2s, B2g, B2r, B2s, P2e and B2e. Depending upon the severity of GVHD /matching, clinical decision regarding the immunosuppressant drug (dim) strategy i.e., drug dose, interval of consecutive drug application and duration becomes important. Such issues are addressed in the developed model. Therefore a(i,j) of equation (8a-3) will be modified as follows.

dim(k) is the amount of immuno-suppressive drug present within the system denote either the amount of applied drug (immunosuppressive) on the day of drug application (drugimdose)or drug retention (drugimret) on the subsequent days as a fraction of the previous day's applied drug (it is assumed that in each corresponding day of drug application a certain amount of drug is cleared from the system). Drug application is started on the day kimsd, nim is the number of cycles, timd is the interval of two subsequent drug applications in days and Cim is the ON/OFF switch for drug application is introduced through above scheme. And dimsensh and dimsensD are two positive fractional numbers representing the immunosuppressive drug sensitivities of the host cells and donor cells respectively (Table 2C). Further the controlling efficacy of immunosuppressive drug will depend on degree of GVHD (G), minimum drug effective set day (GVHDset), drug dose and drug type by affecting CalGVHD with a fractional multiplier (GVHDres). The operational relationship is represented as follows:

Simulation Results

A rigorous simulation exercise with equations (1-8) has been carried out in MATLAB 6.5. Simulation studies may help in revealing the rationale of different therapeutic strategies for HM and HD.

Freely growing tumor

With the initial parametric values as mentioned in Table 1, extensive simulation exercises have been carried out. Simulation studies show that under freely growing condition, there is an exponential growth for the malignant cells while a gradual decay of normal cells (nonmalignant) in the long run (Figure 9) (Table 5). With the increase in malignant cell population, normal cells of the other lineages i.e., erythroid and megakaryocytoid lineages are gradually decreased. Simulation study with the increase in multiplication rate of Sg shows that nonmalignant stem cell population falls earlier than the previous condition due to generation of more number of malignant cells due to the conversion. Disappearance of normal stem has been found earlier with increase in multiplication rate of drug resistive leukoblast cell. These earlier fall in normal stem cell population is due to development of extra toxicity burden to the nonmalignant cell population. Free growth pattern of this simulation resembles the acute leukemia case.

Application of myeloablative/chemotherapy drug regime

Simulation is carried with an assumption that myeloablative (MYL) chemotherapeutic regime is applied on the 30th day after the diagnosis of leukemia (day 1). Drug regime is considered as the amount of total drug = number of cycles × drug dose. Drug dose, cycle and interval between two successive drug application and MYL drug sensitivity for every cell are presented in Table 2. Efficiency of drug also depends on drug retention, in simulation it is considered as 80% of the previous day. In each case drug clearance rate per day has been also considered as 20% of the previous day (Table 2). With MYL drug application there is a gradual decay of different cell populations. After stoppage of drug regimen the malignant cells again start growing. In the peripheral blood, the effect of regaining of cells compared to the precursor level is delayed, however, the effect is more prominent in the stem and progenitor cells present in the bone marrow. This indicates that MYL chemotherapy is unable to bring the leukemic patients in remission (Figure 10) (Table 5). It is to be noted here that though RBC and platelet count reduced due to MYL therapy and after maintaining a reduced level for a longer period of time, both reduces further due to effect of cancer cachexia that has been incorporated into systems equations. Simulation is also carried out with varying number of dose, cycle and interval between two successive drug applications. With the increase in dose, though there is more chance of removal of malignant stem cells however, that may bring down the RBC and platelet population very low and with this level patient may not survive in reality [35-40] i.e., when there is no exogenous RBC and/or platelet transfusion, Tex=0.

High dose Myeloablative chemotherapy application: System was simulated with different combinations of drug doses (Drugmyldose), cycle numbers (nmyl), intervals (tmyld) with the same initial conditions as depicted in Table 1 and followings are noted.

Combination I: Drugmyldose=0.7, nmyl=8, tmyld=15, Tex=0 (Figure 10A)
Combination II: Drugmyldose=0.8, nmyl=8, tmyld=15, Tex=0
Combination III: Drugmyldose=0.7, nmyl=9, tmyld=15, Tex=0
Combination IV: Drugmyldose=0.7, nmyl=8, tmyld=20, Tex=0

Observation with Combination II: When the system was simulated with increase in MYL dose, keeping other conditions unchanged, it was observed that nonmalignant stem cell survived for longer duration (15 days more) than the Combination I. Again after completion with Combination II the nonmalignant stem cell count took more time (136 days more) to reach normal level compared to Combination I. The system was found to continue with lower RBC level than the Combination I.

Observation with Combination III: With Combination III it was found that after completion of MYL therapy, the nonmalignant stem cell count took more time (34 days more) to reach normal level compared to Combination I however quicker by 102 days compared to Combination II. Again simulation with Combination III (keeping other conditions unchanged) it was found that nonmalignant stem cell survived for shorter duration (3 days less) compared to the simulation with Combination II, but 12 days more than with Combination I. However with this combination the system continued to stay with lower RBC cell count compared to Combination I though higher than Combination II.

Observation with Combination IV: Simulation with the increased intervals (Combination IV) it was found that nonmalignant stem cell survived for only 10 days more compared to Combination I; however, 5 days and 2 days lesser than with Combination II and Combination III respectively. Completion of MYL therapy with Combination IV nonmalignant stem cell count took 52 days more than Combination III to reach normal cell count level, however earlier by 54 days than Combination II but longer by 82 days than Combination I. And the system continued to stay with lower RBC cell count compared to Combination I though higher than Combination II or Combination III.

Low dose chemotherapeutic drug application: The system was simulated with the same initial conditions as depicted in Table 1 but with following combinations of drug doses (Drugmyldose), cycle numbers (nmyl) and intervals (tmyld).
Combination V: Drugmyldose=0.35, nmyl=8, tmyld=15, Tex=0 (Figure 10B)
Combination VI: Drugmyldose=0.40, nmyl=8, tmyld=15, Tex=0
Combination VII: Drugmyldose=0.35, nmyl=9, tmyld=15, Tex=0
Combination VIII: Drugmyldose=0.35, nmyl=8, tmyld=20, Tex=0

Observation with Combination VI: When the system was simulated with slightly increased lower dose of chemotherapeutic drug (Combination VI) then it was found that after completion of drug application, the nonmalignant stem cell count took almost same time as with Combination V (only 2 days more in later case) to reach normal cell count level. Again simulation with Combination VI shows that nonmalignant stem cell survived for longer duration (11 days more) than Combination V. It was found that with this combination, the system continued to stay with lower RBC cell count compared to Combination V.

Observation with Combination VII: When the system was simulated with increased cycle number (Combination VII) it was found that nonmalignant stem cell survived for 9 days more compared to Combination V, same as with combination VI; but after completion of Chemotherapy with Combination VII the nonmalignant stem cell count took same time as Combination V to reach normal cell count level though quicker than combination VI (only 2 days earlier). When simulated with Combination VII the system continued to stay with lower RBC cell count though higher than with combination VI.

Observation with Combination VIII: After completion of chemotherapeutic drug application with Combination VIII, the nonmalignant stem cell count took 5 days lesser than both the Combination V and VII to reach normal cell count level; however, 7 days earlier than Combination VI. Again simulation with Combination VIII shows that nonmalignant stem cell survived only 4 days more compared to Combination V, 7 days lesser compared to Combination VI and 5 days lesser compared to Combination VII. It was observed that the system continued to stay with lower RBC cell count compared to Combination V though higher than Combination VII.

In all the above cases though P2r and P2s were reducing initially with the application of drug but with the stoppage of drug they started growing exponentially. Similar observations have been found for malignant cells in mature cell compartment.

Supportive therapy to conventional MYL/chemotherapy – RBC and platelet transfusion

RBC and platelet transfusions are occasionally needed during high dose chemotherapy, as though the number of both RBC and platelet is regained after stoppage of MYL/chemo therapeutic drug, but it may reach so minimal level that a patient would not survive (as shown in Figure 10b). Hence during the application of MYL/chemotherapy drug exogenous RBC and platelet transfusions are needed as supportive therapy (i.e., Tex=1). The effect of such supportive treatment procedure is included in model and studied through simulation. It is assumed that before application of each MYL/chemotherapy dose, the need for supportive therapy is being checked. On the day of MYL drug application if RBCs and/or platelets levels were found below RBCllm and Plateletllm respectively, packed cells of 1,00,000 of RBC and 50,000 of platelets are assumed to be transfused; however, for lower dose of chemotherapeutic drug application, packed cells of 8,00,000 of RBC and 8,000 of platelets are assumed to be transfused (Figure 11). Transfused RBC and platelets have the apoptosis rate which are denoted by aT_RBC (= 0.07747) and aT_Platelet (=0.017) respectively.

System was simulated in presence of RBC and platelet transfusions with Drugmyldose=0.7, nmyl=8 and tmyld=15. It has been observed that though the recovery of the nonmalignant stem cell population into the normal cell count level was delayed by 109 days but they reduced to zero by the toxicity generated by the malignant cells almost at the same time (only 4 days delayed in presence of transfusion compared to Combination I).

System was simulated in presence and absence of RBC and platelet transfusion with Combination V and observed that the recovery of the nonmalignant stem cell population into the normal level was same as the Combination V, but in presence of transfusion with Combination V the nonmalignant stem cell population remained in normal level 5 days more. However, RBC level after the completion of Combination V with transfusion showed slightly lesser than the Combination V alone. Though transfusion helped in maintaining the normal RBC and platelet level in the system but the gap period for transfusion allows the malignancy to increase that, in turn, generates toxicity burden and thereby decrease the nonmalignant cell population.

Cytokine therapy after conventional MYL/chemotherapy

For the treatment of leukemia through immunotherapy, particularly, cytokine based therapies are suggested after MYL therapy. It is assumed that applied cytokine affects only P2g cell types (through differentiation process) and increases B2g cell population which in turn, help in removing leukemic cells from the system. Simulation results are shown in Table 5. With present parametric values simulations suggest that cytokine drug application after MYL are unsuccessful in removing leukemic cells from the system (Figure 12); however, cytokine is applied after low dose chemotherapeutic drug application then malignant cells can successfully be removed. This observation may corroborate the findings of O' Brien et al 1995 [29], where authors showed that residual leukemia after chemotherapy was successfully removed with interferon treatment.

Myeloablative (high) or low dose chemotherapeutic drug applications followed by cytokine drug applications were observed when the system was simulated with following combinations of cytokine dug application after High dose MYL (Drugmyldose=0.7, nmyl=8, tmyld=15, Tex=1) and low dose chemotherapy (Drugmyldose=0.35, nmyl=8, tmyld=15, Tex=1). The initial parametric values including the assumed effects of cytokine drug application are presented in Table 2. It was observed that cytokine drug dose (DrugCYTdose), drug application cycle (nCYT), interval between two successive drug applications (tCYTd), minimum number of WBC (MinB2g) required to start immunity related killing (i.e., to turn killim =1 ) in presence of cytokine drug and malignancy killing efficiency (killeff) play important role in terms of removal of malignancy and maintenance of normal cell population. Rigorous simulation exercises have been carried out to find out the optimum values for cytokine drug dose /application cycles/intervals/minimum WBC level to start immunity related killing (MinB2g) and malignancy killing efficiency (killeff) by cytokine mediated killing.

Application of low dose chemotherapy (Drugmyldose=0.35, nmyl=8, tmyld=15, Tex=1) followed by cytokine drug application: With parametric variations the followings observations are noted.

Observation I: The system was simulated with different drug (DrugCYTdose) level when nCYT =15, tCYTd=25, MinB2g =224 and killeff =80. The cytokine drug was started from day 150 i.e., 5 days after the application of last dose of chemotherapy. It was observed that drug resistive stem cell disappeared from the system on day 568, 661 and 706 when DrugCYTdose were 0.8, 0.76 and 0.755 respectively. i.e., with application of increased drug level the removal of malignancy becomes faster. And in all the cases non-malignant cell population in all compartments were found either maintaining normal cell count level or growing to reach normal cell count level.

Observation II: The system was simulated with different killeff when DrugCYTdose=1.2, nCYT =15, tCYTd=25 and MinB2g = 224. The cytokine drug was started from day 160 i.e., 15 days after the application of last dose of chemotherapy. It was observed that drug resistive stem cell disappeared from the system on day 320, 341 and 387 when killeff were 100, 90 and 80 respectively. i.e., with application of increased killing efficiency the removal of malignancy becomes faster. And in all the cases, non-malignant cell populations in all level were found either maintaining normal cell count level or growing to reach normal cell count level.

Observation III: The system was simulated with different minimum WBC level (MinB2g) to start immunity related killing when DrugCYTdose=1, nCYT =15, tCYTd=25 and killeff =80. The cytokine drug was started from day 160 i.e., 15 days after the application of last chemo dose. It was observed that drug resistive stem cell disappeared from the system on day 393, 393 and 412 when (MinB2g) was 223, 222 and 224 respectively; i.e., with the lowering the minimum effective WBC level ((MinB2g) for starting the immunity related killing, malignancy removal becomes faster. And in every the cases and non-malignant cell population in all level were found either maintaining normal cell count level or growing to reach normal cell count level. But when simulated with (MinB2g) =225 with other conditions unchanged malignancy found to increase and all the non-malignant cell populations reduced to zero (stem cell disappeared from the system reduced to zero on day 711). In this case the minimum MinB2g has been found as 224. In this present parametric setting above this value the system will fail to control malignancy. In case of cytokine application (DrugCYTdose = 1) after chemotherapeutic strategy with transfusion though there is a persistence of B1r and B1s cells having a number >2000 on day 1500 however with a decaying tendency.

Application of high dose MYL (Drugmyldose=0.7, nmyl=8, tmyld=15, Tex=1) followed by cytokine drug application: With parametric variations the followings observations are noted.

The system was simulated with different drug (DrugCYTdose) levels when nCYT =15, tCYTd=25, MinB2g =224 and killeff =80. The cytokine drug was started from day 170 i.e., 15 days after the application of last MYL dose. It was observed that keeping all initial conditions as depicted in Table 1 unchanged with the increase in DrugCYTdose level, nonmalignant cells remained in the system for longer time and recovery of nonmalignant stem cell to normal cell count level were faster (nonmalignant stem cell reached to normal cell count after completion of MYL on day 338th, 304th and 291st when DrugCYTdose level were assumed as 1, 4 and 6 respectively). For DrugCYTdose =6 the malignancy were found to disappear from all compartment (resistive stem cell disappeared from the system on day 363) and normal cell populations were found either in normal cell count level or growing to reach normal cell count level. Though malignancy were found to be disappeared from all compartmental level with DrugCYTdose = 6 but it may be unrealistic (as DrugCYTdose×mCYT =6×0.07=0.42 i.e., 42% extra growth in Leukoblast per day).

The failure to control malignancy with MYL followed by cytokine is due to the fact that high dose MYL bring down normal leukoblastic cells into such a lower level that may not be possible to boost to such high level so quickly to control malignancy. In present system model it has been assumed that with high dose MYL application when number of P2g<30, the Sg>dt and P2gdt are increased to 30 days and 20 days respectively than their normal value; whereas Sgdr and P2gdr are decreased to 0.0165×{3-(mSg + aSg)} cells/day and 0.765×{1-(mP2g + aP2g)} cells/day respectively) with respect to their normal value and aSg gets increased (to 0.48 cells/day) with respect to its normal value.

Allogenic hematopoietic stem cell (HSC) transplantation with 100% HLA matching

For simulation it is assumed that HSC transplantation has been done 25 days (on 180th day after the day of diagnosis) after the completion of last dose of myeloablative (MYL) drug regime (cycles 8, drug dose 0.7 and interval 15 days). The day for HSC transplantation is chosen when drug level within the system reaches below a desired level (say < 1% of the applied drug). It is to be noted here that in simulation MYL drug was applied on day 30th, 45th, 60th, 75th, 95th, 115th, 135th, 155th. The unequal gap between successive drug application were due to transfusion of RBC and Platelet on 90th, 110th, 130th and 150th days for low RBC and Platelet count. The assumed characteristics parameters of the transplanted HSC are presented in Table 3. As there is, no GVHD, so the maturated donor lymphocytes that are produced from the transplanted HSC will kill only the residual malignant cells with a factor kill b2e→hm = 7.142 (Table 4). With the present simulation condition, it has been observed that around ~324 days after HSC transplantation all the malignant cells were totally disappeared from the recipient's system, no relapse was noticed even after 2000 days (Figure 13) (Table 5). Chimera (Mixture of donor and host cells) in different cellular level is observed. Simulation has been carried out with different killing efficiency of lymphocytes originated form transplanted cells.

Observation 1: Similar type of observations (as in Figure 13) were noted for killb2e→hm = 7.14, and 7.13 with other conditions unchanged.

Observation 2: With the other conditions unchanged simulation exercises with varying killing efficiencies show that killb2e→hm = 7.13 is the minimum required value for the complete removal of malignancy from the system. It signifies that the exogenous cell should have minimum killing efficiency for complete malignancy removal. Simulation with the value below 7.13 there is a decreasing trend of normal cell level due to toxicity generated by the increasing level of malignant cell. However, when the system was simulated with increased number of transplanted HSC (Se = 12) having killb2e→hm = 7.125, the malignancy is removed from the system and normal cell population is sustained. It signifies that donor stem cells with lower killing efficiency can be compensated with increasing its number.

Observation 3: With other conditions unchanged when simulation is carried out with Se = 12 and killb2e→hm = 7.125; however, transplantation is assumed to be on 190th day after the day of diagnosis, malignancy exist, in all level.

Observation 4: With other conditions unchanged donor stem cells with killb2e→hm = 7.125 and transplantation is assumed to be on 190th day after the day of diagnosis with more number of HSC (Se = 20) there is a complete removal of malignancy from the system. It signifies that the delay (in days) in transplantation may compensated (for removing leukemic cells) by increased donor cell number by multiple source of cord blood engraftment.

Therefore delay in transplantation can be compensated either with increased donor cell number or donor cells having property of generating lymphocytes of higher killing efficiency. Moreover, in the model, the application of MYL drug strategy (drug dose and cycle number and interval) can be changed/adjusted. This allows the investigator to check the condition of attaining a minimal tumor load suitable for transplantation.

Allogenic hematopoietic stem cell (HSC) transplantation with HLA mis-matching

HLA mismatch is manifested as graft versus host diseases (GVHD). Simulation is carried out considering chronic GVHD with two grades: high and low (Figure 14). The condition of HLA mismatch can also be studied through this model. The effects of HLA mismatch are incorporated in the model using several killing factors: exogenous lymphocytes developed from transplanted HSC reacted against the recipient (host) malignant cells (killb2e→hm) (same as for 100% HLA matching case) and host normal cells (killb2e→g). The effect of GVHD is manifested through development of GVHD against normal tissue host like liver, gut and lungs. And the intensity of such GVHD related tissue damage can be measured using a marker (MR), where MR= f(B2e, G).

Simulation has been carried out in different conditions and results are presented as below.

Observation 1: Simulation study shows that in case of transplantation with HLA mismatch, removal of malignancy is found to be much faster than 100% HLA match case (Table 5). For example Sr reduced to zero on day 504th when transplanted with 100% HLA matched donor cell; however, when transplanted with low and high mismatched donor, removal of Sr is on day 224th and 217th respectively.

Observation 2: As killb2e→hm increases while other factors unchanged normal cell population also reduces faster (Table 5) along with the faster removal of malignancy.

Observation 3: Keeping other conditions unchanged with delay in transplantation though normal cell survived for longer period but malignancy removal was found delayed relative to the previous cases.

Observation 4: Transplantation with high mis-match shows faster removal of malignancy faster than transplantation with lower mis-matched donor and it also reduces normal cell population earlier even more than free growth condition.

With the present parametric settings (Table 4) it is observed that with the increase in HLA mis-matching, GVHD appears much earlier (Figure 14). The effect is much more prominent on the mature cell level than the precursor cell level.

Application of immunosuppressive drug scheduling for HLA mismatch cases

When allogenic hematopoietic stem cells (allo-HSC) with a certain degree of HLA mis-match are transplanted, patients need to be treated with an immunosuppressive drug for the management of GVHD. It is worthwhile to mention here that depending upon the severity of mismatch, determination of immunosuppressive (IMS) drug dose and drug schedule could be the important aspects for clinical cases. The developed model has the provision to test the efficacy of the IMS drug for the management of GVHD. This model has the flexibility of choosing IMS drug scheduling strategy.

In simulation different IMS drug strategies have been applied to different HLA mis-match cases as mentioned previously. In all the cases, allo-HSC transplantation has been performed on day 180 after completion of conventional MYL drug regime. Simulations have been carried out considering two different types of IMS drugs (i.e., IMS-1 and IMS-2). Both types of drugs have the same mechanism of action; however, IMS-2, signifies higher dose (double) than IMS-1 and once one drug is applied it continued throughout the simulation period. Strategy started 30 days after HSCT and IMS-2 drug strategy started 60 days after HSCT (Figure 15).

IMS drug dose, drug application interval, drug retention and the considered sensitivities of these drugs to different cells are depicted in Table 2C and results of such simulations are presented in Table 5. Simulation study indicates that IMS drug application in controlled level can keep GVHD under control to remove the malignancy from the system and thereby increase the sustaining period of normal cell population.

However, with the application of excessive IMS drug dose (keeping other conditions unchanged) malignancy may relapse in system.

Remission by chemotherapy /chemo-responsive leukemia

So far our simulation exercises have shown the management strategies for chemo-unresponsive leukemia (worst) case. However, several leukemia cases respond to the conventional chemotherapy regimes. Simulation study with the increase in drug sensitivity ~10 times for the malignant cells only as mentioned in Table 2 and with the application of 8 cycles of high dose MYL in presence of transfusion can bring the considered worst case (as in Table 1) into remission; however, it does not bring to remission when 6 cycles of the same therapeutic regime is applied.

Simulation study with the developed model also carried out to show the features of chemo-responsive leukemia (remission) cases with change in some initial parametric values (Case II). Free growth condition of such case is initialized by setting mSs = 0.07, mSr = 0.075, aSs = 0.032, aSr =0.032, mP2s = 0.0345, mP2r =0.0345, aP2s = 0.00359, aP2r =0.00359, 1.3 times, 1.3 times and 1.5 times of Srdr, P2sdr, and P2rdr of worst case respectively while rest of the values are kept unchanged as mentioned in Table 1 (Figure 16). Remission is possible with the application of 6 cycles of high dose of MYL with 5.5 times MYL drug sensitivity towards malignant cells as mentioned in Table 2. Simulation study with this case depicts that the number of malignant stem cells have slower growth rate than the worst case and the leukemic/malignant cells have greater drug sensitivity

Keeping all the parametric values and the MYL drug regime same as of the Case II remission case while kept back the parametric value of mSr and aSr as of the worst case fails to bring the system into long-term remission and simulation study depicts the presence of residual disease (Figure 17).

Chronic leukemia cases

Previously simulation studies resemble acute leukemia cases. However, change in initial parametric values in the developed model can show the features of chronic leukemia cases. Chronic myeloid leukemia (CML) is considered as the differentiation defect with an arrest of cells at different stages of morphological maturation [40]. Free growth condition is initialized by changing Srdr, P2sdr and P2rdr to 1.35 times, 1.55 times and 1.95 times respectively of worst case while rest of the values are kept unchanged as mentioned in Table 1 (Figure 18). In this case drug application is considered daily with drug dose 0.037 (6/1000 th of total high dose MYL drug) for 1500 days. It is to be noted here that in CML (chronic myeloid leukemia) cases are generally treated with the targeted drug imatinib and this drug is being applied daily [41-43]. Simulation is carried out with two different drug sensitivities for malignant leukemic cells – one simulation is done with the same sensitivity as mentioned in Table 2 (Case I) and another simulation is carried out with the 5.5 times of drug sensitivity for malignant cells (Case II). However in both the simulations it is also considered that such drug strategy has minimal effect on the normal cells of the hematopoietic system. Simulations studies depict that with the increase in drug sensitivity of malignant cells can able to bring the leukemic free survival in the long-run.

Discussion

In reality medical treatments are dependent on two approaches – population based randomized clinical trials and qualitative judgments of the subjective experiences of clinicians. Hence treatment dealing procedure of an individual patient becomes empirical. However, physicalist approaches demand for mechanistic explanation. Mathematical models with dynamical nature may provide the way-out to an understanding of a disease process in a much more quantitative manner. It is also true that till now the qualitative judgments and subjective experiences cannot be denied while dealing with the individual clinical cases [44,45]. This fact is mostly denied as mere information and/or unable to comprehend its importance with the cultural training of physicalist school of thought [46]. Ironically majority of developed mathematical models are done with abstract assumptions of variables that are unable to hold any applicative significance to clinical practice for dealing of individual patients though the model behavior synchronizes some real data. Biological discipline has a long history of training program to comprehend the importance of biological variability and analyzing the natural phenomena. Though physiology and medicine both deal life processes, but due to technical and instrumental dependency, both follow reductionism; contrarily, biological analysis is viewed as and confined within the subjective experiences and as a result there is absence of time varying dynamical data (of individual cases) in medical literature. Systems Biology and its translation towards Systems Medicine seems to be essential to meet this gap. Towards this, presently, development of analytical model with rationality can be done [47,48]. Simulation studies present here may reflect the subjective experiences of clinicians and the developed model may be helpful in providing the mathematical rationality towards their domain knowledge.

Previously DODE is being applied by several workers for mathematical modeling of both chronic and acute leukemia. Hence it is imperative to state that the same mathematical strategy can be applicable for both chronic and acute leukemia with varying in parametric values for different variables. Such mathematical strategies are also used for the understanding of the efficacy of different therapeutic strategies like myeloablative chemotherapy, tyrosine kinase inhibitor and/or T-cell based therapy, but different workers have used different variables in the development of the model with a focus on different therapeutic strategies. Hence each model is being apart from other. However, in clinic there is a need to evaluate and/ or predict the response of a particular individual patient to a particular therapeutic strategy among the different available therapeutic options. Hence there is an immense need for the development of an analytical model that can make a comparative assessment among the different available therapeutic options. To address this issue we have developed our model using DODE with an incorporation of feedback at different systems level. Our model can be aptly be fitted with both acute and chronic leukemia cases by setting the initial parametric values of the variables. Similarly setting of different drug parameters can be used for the assessment of the efficacy of different drugs that are being used currently in conventional clinical practice. The model has the flexibility in chemotherapeutic drug application (dose, duration and cycles), so the model can also be suited for the assessment and prediction regarding the dynamical features of different patients with different patho-physiological states.

The present model includes the complex nonlinear behavior of the different time varying system elements of the hematopoietic system and their inter-relationships. The complex inter-relationships have been reflected when observed with minor parametric variations. In freely growing condition it has been found that minor increase in multiplication rate of any normal cell type failed to increase its population rather it decreases. It is due to the consideration of mutability of normal cells to the malignant phenotypes. That in turn generates toxicity cachexia to bring down the normal cell population. This may be major cause of infection and bleeding in leukemia.

With the inclusion of such complex nature in the system model, the management of MYL drug application becomes more difficult for the maintenance of normal cell count of the normal cell population in one hand and removal of malignancy on other hand. With increase in MYL dose amount the malignant cell count can be found to decrease. However, simulation shows that excessive drug application for removal of malignancy may collapse the normal cell population due to the drug related toxicity which in turn, enhances the proneness to infection and bleeding. This reflects real clinical scenario. Contrarily, residual disease may exist due to application of lower dose of chemotherapy.

So application of MYL drug dose should be controlled judiciously with a tradeoff between maintenance of normal cell count for longer period of time versus elimination of malignant cell populations from system. This can only possible if there is an increase in drug clearance rate and/or drug sensitivity for malignant calls.

To circumvent the MYL drug related toxicity, RBC and/or platelet transfusion is the regular practice in clinic. This has been included in our system model. Simulation study shows that application of transfusion strategy may make a delay in the consecutive MYL drug application; which in turn, may allow the malignant cells to grow. This again becomes detrimental for the patient due to toxicity generation by the malignant cells.

Our simulation study shows failure of MYL therapy with or without transfusion. This signifies that we have considered a worst clinical case. Simulation has been carried out further to reveal the expected management procedure for drug resistive worst leukemia case. Considering this worst case, further simulation shows that the application of MYL drug followed by cytokine is also unable to bring the desired success under the same leukemic condition; however the application of lower doses of chemotherapeutic drug (preparatory regime) followed by cytokine strategy may bring the desired therapeutic success. This is due to prior application of MYL before cytokine application is unable to potentiate the immune-competent lymphocyte due to removal of cells of the existing leukocytic lineage. This observation corroborates the earlier observation which indicates that the leukemia free survival is dependent on the existing functional T-cell [7,12].

The considered cell count can be routinely assessed by hematological investigation, different immunological parameters like killing efficiency of lymphocytes can be assessed with different immunological investigations, and the considered chronic GVHD can be assessed by different qualitative clinical investigations.

The present model has the provision to test the efficacy of transplantation under different conditions in terms of HLA matched condition, transplantation day, delay in transplantation, malignancy killing efficiency by the lymphocyte of donor origin and number of transplanted cells. The considered cell number for transplantation with respect to volume can be equated in terms of body weight. With 100% HLA matching, simulation with donor cell having lower killing efficiency no improvement in terms of malignancy removal was noticed though there was no GVHD effect but GVL is failed. However increase in transplanted donor stem cell number (in reality this can be done with multiple donor engraftment) keeping other conditions unchanged, are able to remove the malignancy by developing GVL to threshold level. Again when simulated with delay in transplantation day there is a failure in controlling leukemia. Since during delay time leukemic cells reached to such level that develop toxicity which in turn, kills donor cell population. So determination of minimal tumor load on the day of transplantation is another major consideration. With this model, determination of such criteria can also be assessed. Hence, simulations suggest that this failure may be compensated by increasing initial transplanted cell number or donor cell having higher killing efficiency. This can be achieved by increasing the lymphocyte population of donor origin by applying cytokines. Model has such provision to incorporate such effect by mDCYT.

With increase in HLA mis-match keeping other factors unchanged removal of malignancy also becomes faster. Though transplantation with high GVHD are able to remove malignancy faster than transplantation with lower mis-matched donor but it also reduces normal cell population faster even more than free growth condition. To tackle this, immunosuppressive drug application is suggested. Model has the provision in the assessment of such therapeutic strategy. Even the model flexibility is also made in a manner to test the second choice of immunosuppressive drug application, if first one fails. Our simulation again suggests that judicious choice of immunosuppressive drug application may eliminate both malignancy and GVHD from the system, as malignancy may predominate in the system due to application of excessive amount of immunosuppressive drug. Excessive application of immunosuppressive drug may enhance the susceptibility of the patient to infection and bleeding. This is due to removal of functional donor lymphocytes by the excessive amount of immunosuppressive drug. However such efficacy of immunosuppressive drug will also depend further on the degree of HLA mis-match.

It is worthwhile to mention here that with the availability of individual patient's dynamical data, fine adjustment of simulation studies with this model could be possible. It is to be mentioned here that our developed model is flexible in nature. In our model all the parametric values of every variable can be changed and can be adjusted accordingly to the individual clinical cases as well as choice of the investigator. The major objectivity of the present work is to provide a mathematical/analytical platform to make a comparative assessment among the different available therapeutic options [47,48]; so that, a clinician can make a prior assessment regarding the efficacy of a particular therapeutic strategy before applying it to an individual patient. This model can also be aptly fitted to test the efficacy of different therapeutic strategies on different individual cases with different types of leukemia. Moreover, model can be applied for prognostic assessment in individual cases. This quantitative assessment tool will help to transcendent the empirically based clinical practice in future.

Conclusion

Recent time several therapeutic modalities for the treatment of leukemia are being suggested by the scientific community; however, due to individual variation in response to a particular therapy it is difficult to rationalize which therapeutic modality would be the choice of the treatment for a particular individual patient. The developed model can help in the apriori assessment regarding the efficacies of different therapeutic modalities. Simulation studies with this developed model indicate that under the same leukemic condition, when conventional chemotherapeutic regime is failed to provide long-term leukemia free survival immune-boosting mechanism by any means like cytokine application or allo-HSCT could be the determining factor for long-term leukemia free survival. Our simulation study also provides the insight that the time of immune stimulation is another determining factor for the success of leukemia therapy. The developed model is flexible in nature. Clinicians can change any parametric value according to the need of the individual clinical cases and make an in-silico assessment regarding the suitability of a particular therapeutic regime depending on the individual leukemia cases before applying it in reality.

Acknowledgement

Authors acknowledge the eminent founder members of Society for Systems Biology & Translational Research for their valuable comments to enrich this work.

5 Colijn C, Mackey MC (2005b) A mathematical model of hematopoiesis -II. Cyclical neutropenia. J Theo Biol 237: 133-46.
6 Michor F, Hughes TP, Iwasa Y, Branford S, Shah NP, et al. (2005) Dynamics of chronic myeloid leukemia Nature 435: 1267-70.
16 Horowitz MM, Gale RP, Sondel PM, Goldman JM, Kersey J, et al. (1990) Graft-versus-leukemia reactions after bone marrow transplantation. Blood 75: 555-62.
20 Majumder D, Bandyopadhyay D, Chandra S, Mukhopadhayay A, Mukherjee N (2005) Analysis of HLA class Ia transcripts in human leukemia. Immunogenetics 57: 579-89.
25 Argiles JM, Busquets S, Stemmler B, Lopez-Soriano FJ (2014) Cancer Cachexia: Understanding the Molecular Basis. Nature Rev Cancer 14: 754-62.
27 Dignan FL, Amrolia, Clark A, Cornish J, Jackson G, et al. (2012) Diagnosis and management of chronic graft-versus-host diseases. Br J Haematol 158: 30-45.
33 Talpaz M, Hehlmann R, Quintas-Cardama A, Mercer J, Cortes J (2013) Re-emergence of interferon-α in the treatment of chronic myeloid leukemia. Leukemia 27: 803-12.
34 Ross JF, Crockett CL Jr, Finch SC, Emerson CP (1954) The mechanism of anemia in leukemia and malignant lymphoma. BMQ 5: 8-10.
35 Teuffel O, Stanulla M, Cario G, Ludwig WD, Rottgers S, et al (2008) Anemia and survival in childhood acute lymphoblastic leukemia. Haematologica 93: 1652-7.
37 Cancer Connect.com, Anemia: A manageable side effect leukemia, Assessed on 14 Feb, 2015.
38 Rickles FR, Falanga A, Montesinos P, Sanz MA, Brenner B, et al. (2007) Bleeding and thrombosis in acute leukemia: what does the future look like? Thromb Res120: S99-106.
39 Cancer.gov, Chronic myelogenous leukemia treatment. Assessed on 08 March, 2015.

Journal of Computational Systems Biology

Tables at a glance
Table 1
Table 2
Table 3
Table 4
Table 5
Figures at a glance
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 10
Figure 11
Figure 12
Figure 13
Figure 14
Figure 15
Figure 16
Figure 17
Figure 18
Figure 1: Different cellular components of model
Figure 2: Considered feedback system for maintenance of normal hematopoietic homeostasis
Figure 3: Cachexia effect by the developed malignant cells
Figure 4: Scheme for application of Myeloablative (MYL) therapy
Figure 5: Scheme for application of supportive (RBC/platelet transfusion) therapy in Myloablative (MYL) chemotherapy
Figure 6: Considered GVHD effects. Tox (0) represents negligible GVHD effect on RBC of host/donor origin
Figure 7: Scheme for application of Myloablative (Myl)/ chemotherapy along with supportive (RBC/platelet transfusion) therapy followed by cytokine therapy
Figure 8: Algorithmic scheme for MYL therapy followed by HLA matched/mismatched HSC transplantation and immunosuppressive application
Figure 9: Plots showing the growth of different cell types of the hematopoietic system under freely growing condition. In (a), the dynamical nature of different types of stem cells [in a (i)], erythroblast cells [in a (ii)], leukoblast cells [in a (iii)] and megakaryoblast cells [in a (iv)]; and in (b), the different types of RBCs [in b (i)], WBCs [in b (ii)] and platelets [in b (iii)] are shown. In all the plots of (a) and (b), x-axis represents 'Days' and y-axis represents 'Counts'
Figure 10: Plots showing the effect of myeloablative drug regime (cycle (nmyl) 8, drug dose (Drugmyldose) 0.7 (MYL) /0.35 (chemo), interval (tmyld) 15 days) (Combination I and Combination V) on the growth of the cells of hematopoietic system. In (a) high dose Myeloablative and low dose chemotherapeutic drug profile, in (b) the dynamical nature of different types of stem cells [in b (i)], erythroblast cells [in b (ii)], leukoblast cells [in b (iii)] and megakaryoblast cells [in b (iv)]; in (c), the different types of RBCs [in c (i)], WBCs [in c (ii)] and platelets [in c (iii)] are shown. In all the plots of (b) and (c), x-axis represents 'Days' and y-axis represents 'Counts'
Figure 11: Plots showing the effects of intermittent RBC and platelet transfusion with myeloablative drug (in A) (Combination I) and Chemotherapeutic drug (in B) (Combination V) regime (cycle (nmyl) 8, drug dose (Drugmyldose) 0.7 (MYL) /0.35 (Chemo), interval (tmyld) 15 days) on the growth of the cells of hematopoietic system. In (a) high dose MYL and low dose chemotherapeutic drug profile; in (b) the transfused RBCs and transfused platelets during MYL and chemotherapeutic drug application; in (c), the dynamical nature of different types of stem cells [in c (i)], erythroblast cells [in c (ii)], leukoblast cells [in c (iii)] and megakaryoblast cells [in c (iv)] in both cases and in (d), the different types of RBCs [in d (i)], WBCs [in d (ii)] and platelets [in d (iii)] in both cases are shown. In all the plots of (c) and (d), x-axis represents 'Days' and y-axis represents 'Counts'
Figure 12: Effects of cytokine therapy after myeloablative/chemotherapeutic drug regime (cycle 8, drug dose 0.7 (MYL)/0.35 (chemo), interval 15 days with RBC and platelet transfusion). In (a) plots represent cytokine drug application profile and the change in cell mediated immunity (CMI) profile due to application of cytokine; in (b), plots showing the dynamical behavior of stem cells [in b(i)], erythroblast [in b(ii)], leukoblast [in b(iii)] and megakaryoblast [in b(iv)] in both cases; and in (c), plots represents RBCs [in c(i)], WBCs [in c(ii)] and platelets [in c(iii)] in both cases. In all the plots of (b) and (c), x-axis represents 'Days' and y-axis represents 'Counts'
Figure 13: Effect of allogenic HSC transplantation with 100% HLA matching with number of donor stem cell (Sesub>) = 11, kill b2e→ hm = 7.142, day of transplantation 180 from the day of diagnosis. In (a), plots showing the dynamical behavior of stem cells (i), erythroblasts (ii), leukoblasts (iii) and megakaryoblasts (iv); in (b), plots showing the dynamical behavior of RBCs (i), WBCs (ii), and platelets (iii). In all the plots of (a) and (b), x-axis represents 'Days' and y-axis represents 'Counts'
Figure 14: Effect of allo-HSC transplantation with two different grades [low in (A) and high in (B)] of HLA mis-matching. In (a), plots showing MYL drug profile (i) and transfused RBC and platelets profile (ii); In (b) plots showing the behavior of stem cells (i), erythroblasts (ii), leukoblasts (iii),megakaryoblasts (iv), RBCs (v), WBCs (vi), platelets (vii), damaging effect on the concerned tissue due to GVHD (viii). In all the plots from (i) to (vii) of (b), x-axis represents 'Days' and y-axis represents 'Counts' and in plot (viii) of (b), x-axis represents 'Days' and y-axis represents 'Tissue damaging intensity'
Figure 15a: Immuno-suppressive drugs profiles in HLA mis-matching condition
Figure 15b: Effect of immuno-suppressive drugs on the dynamics of progenitor cells in different grades (Low and high) of HLA mis-matching condition. In all the plots, x-axis represents 'Days' and y-axis represents 'Counts'
Figure 15c: Effect of immuno-suppressive drugs on the dynamics of mature cells in different grades (Low and high) of HLA mis-matching condition. In all the plots, x-axis represents 'Days' and y-axis represents 'Counts'
Figure 15d: GVHD management (damaging effect on the concerned tissue) with immuno-suppressive drugs (IMS1 and IMS2) in different grades [Low in (A) and high in (B)] of HLA mis-matching condition. In all the plots, x-axis represents 'Days' and y-axis represents 'Tissue Damaging Intensity'
Figure 16: Effect of conventional chemotherapy in remission leukemia case. In (A), plots showing MYL drug profile (i) and transfused RBC and platelets profile (ii); In (B) plots showing the behavior of stem cells (i), erythroblasts (ii), leukoblasts (iii), megakaryoblasts (iv), RBCs (v), WBCs (vi), platelets (vii) under freely growing condition (a) and with conventional chemotherapy (b). In all the plots of (b), x-axis represents 'Days' and y-axis represents 'Counts
Figure 17: Presence of residual disease /leukemia after conventional chemotherapy. In (A), plots showing MYL drug profile (i) and transfused RBC and platelets profile (ii); In (B) plots showing the behavior of stem cells (i), erythroblasts (ii), leukoblasts (iii), megakaryoblasts (iv), RBCs (v), WBCs (vi), platelets (vii) under freely growing condition (a) and with conventional chemotherapy (b). In all the plots of (B), x-axis represents 'Days' and y-axis represents 'Counts'
Figure 18: Chronic leukemia case with low daily doses of chemotherapy (Imatinib). In (A), plots showing chemotherapy drug profile; In (B) plots showing the behavior of stem cells (i), erythroblasts (ii), leukoblasts (iii), megakaryoblasts (iv), RBCs (v), WBCs (vi), platelets (vii) under freely growing condition (a), no remission with blast crisis after chemotherapeutic treatment (b) and long-term remission after treatment with a drug with higher sensitivity for malignant cells (c). In all the plots of (B), x-axis represents 'Days' and y-axis represents 'Counts'
Variables Stem cell
Normal type Sensitive type Resistive type
Symbol Value Symbol Value Symbol Value
Initial number Sg 11 Ss 2 Sr 1
Multiplication rate mSg 0.065 mSs 0.075 mSr 0.085
Apoptosis rate aSg 0.032 aSs> 0.0315 aSr 0.031
Differentiation rate Sgdr cf2*(cf1-( mSg + aSg )) Ssdr cf2*(1-( mSs + aSs )) Srdr cf2*(1-( mSr + aSr ))
Differentiation delay
time
Sgdt 14 Ssdt 14 Srdt 14
Conversion rate CSgs 0.0052 CSsg 0.0012 CSrg 0.0022
CSgr 0.0068 CSsr 0.0053 CSrs 0.0032
A. Stem cell compartment
Variables Erythroblast cell
Normal type Sensitive type Resistive type
Symbol Value Symbol Value Symbol Value
Initial number P1g 380 P1s 15 P1r 5
Multiplication rate mP1g 0.013 mP1s 0.0131 mP1r 0.0132
Apoptosis rate aP1g 0.01339 aP1s 0.01338 aP1r 0.013375
Differentiation rate P1gdr cf4*(cf1-( mP1g + aP1g )) P1sdr cf3*(1-( mP1s + aP1s )) P1rdr cf3*(1-( mP1r + aP1r ))
Differentiation delay
time
P1gdt 2 P1sdt 2 P1rdt 2
Conversion rate CP1gs 0.0023 CP1sg 0.0012 CP1rg 0.0012
CP1gr 0.0012 CP1sr 0.0023 CP1rs 0.0023
Leukoblast cell
Initial number P2g 40 P2s 50 P2r 10
Multiplication rate mP2g 0.0345 mP2s 0.033 mP2r 0.043
Apoptosis rate aP2g 0.0359 aP2s 0.00358 aP2r 0.003575
Differentiation rate P2gdt cf4*(cf1-(mP2g + aP2g )) P2sdr cf3*(1-( mP2s + a P2s )) P2rdr cf3*(1-( mP2r + aP2r ))
Differentiation delay time P2gdt 7 P2sdt 7 P2rdt 7
Conversion rate CP2gs 0.0012 CP2sg 0.0012 CP2rg 0.0023
CP2gr 0.0023 CP2sr 0.0023 CP2rs 0.001
Megakaryoblast cell
Initial number P3g 90 P3s 5 P3r 5
Multiplication rate mP3g 0.075 mP3s 0.015 mP3r 0.015
Apoptosis rate aP3g 0.0802 aP3s 0.00805 aP3r 0.00805
Differentiation rate P3gdr cf4*(cf1-(mP3g + aP3g )) P3sdr cf3*(1-(mP3s + aP3s )) P3rdr cf3*(1-( mP3r + aP3r ))
Differentiation delay time P3gdt 3 P3sdt 3 P3rdt 3
Conversion rate CP3gs 0.0012 CP3sg 0.0010 CP3rg 0.0010
CP3gr 0.0023 CP3sr 0.0012 CP3rs 0.0012
B. Progenitor cell compartment
Variables Erythroblast cell
Normal type Sensitive type Resistive type
Symbol Value Symbol Value Symbol Value
Initial number B1g 4750000 B1s 1875000 B1r 62500
Apoptosis rate aB1g 0.00007747 aB1s 0.000005746 aB1r 0.000004745
Conversion rate CB1gs 0.0023 CB1sg 0.0012 CB1rg 0.00121
CB1gr 0.0012 CB1sr 0.0001 CB1rs 0.0012
Leukocytes cell
Initial number B2g 4000 B2s 5000 B2r 1000
Apoptosis rate aB2g 0.012999996 aB2s 0.012999996 aB2r 0.012999996
Conversion rate CB2gs 0.0012 CB2sg 0.002 CB2rg 0.0012
CB2gr 0.0021 CB2sr 0.0012 CB2rs 0.001
Platelet cell
Initial number B3g 450000 B3s 25000 B3r 25000
Apoptosis rate aB3g 0.00017 aB3s 0.000016 aB3r 0.000151
Conversion rate CB3s 0.0012 CB3sg 0.0013 CB3rg 0.0012
CB3gr 0.0023 CB3sr 0.0022 CB3rs 0.0023
† cf1 (=3), cf2 (=0.0365), cf3 (=0.0165) and cf4 (=1) are calibration factors.
Unit of delay is day and unit of rate is cells day-1. The parametric values of normal cells are taken from Dhar et al, 2012.
C. Terminally mature cell compartment
Table 1: Initial parametric values of the host hematopoietic system
Variable Symbol Value Unit Comments
MYL/Chemotherapeutic drug dose drugmyldose 0.7 (MYL)
0.35 (Chemo)
Fraction of cell killing
per day
Assumed value
Duration between two successive MYL/ Chemotherapeutic drug application tmyld 15 (MYL),
15 (Chemo)
Days Drug ( MYL/Chemo) applied once in
15 days
Number of cycle nmyl 8 Number Number of repetition of drug application
Amount of MYL/ Chemotherapeutic drug retention on the subsequent day of drug application drugmylret 80% % of drug of the
previous day
Retention of 80% drug of the previous
day drug amount
MYL drug sensitivity of sensitive type stem cells MylSs 0.4 (MYL) % of cells

40% cells are sensitive to MYL/Chemo
drug.

Values are assumed.

0.4(Chemo)
MYL drug sensitivity of resistive type stem cells MylSr 0.3 (MYL) % of cells

30% cells are sensitive to MYL/Chemo
drug.

Values are assumed.

0.3 (Chemo)
MYL drug sensitivity of sensitive type P2s,B2s cells MylP2s MylB2s 0.4 (MYL) % of cells

40% cells are sensitive to MYL/Chemo
drug.

Values are assumed.

0.4 (Chemo)
MYL drug sensitivity of resistive type P2r, B2r cells MylP2r MylB2r 0.3 (MYL) % of cells

30% cells are sensitive to MYL/Chemo
drug.

Values are assumed.

0.3 (Chemo)
MYL drug sensitivity of all good cells & P1r, B1r, P3r, B3r, P1s, B1s, P3s, B3s cells

MylSg MylPlg MylP2g

MylP3g MylBlg MylB2g

MylB3g MyPlr MylBlr

MylP3r MylB3r MylPls

MylBls MylP3s MylB3s

0.2 (MYL) % of cells

20% cells are sensitive to MYL/Chemo
drug.

Values are assumed.

0.2 (Chemo)
A. MYL/Chemo drug related parameters
Variable Symbol Value Unit Comments
Cytokine drug dose drugCYTdose 1 Fraction of cell killing per day Assumed value.
Duration between two successive Cytokine drug application tCYTd 15 Days Drug applied once in 15 days
Number of cycle tCYT   Number Number of repetition of drug application
Amount of Cytokine drug retention on the subsequent day of drug application drugCYTret 10% % of drug of the previous day Retention of 10% drug of the previous day drug amount
Cytokine induced effect (additive) on the multiplication rate of recipient’s progenitor lukocytic lineage cell mCYT 0.07 Cells/day Assumed value

B. Cytokine drug related parameters

Variable Symbol Value Unit Comments
Immuno-supressive drug dose drugimdose 0.01 & 0.02 Fraction of cell killing
per day
Drug doses (0.01/0.02) are chosen as per
requirement.
Duration between two successive Immuno-supressive drug application timd 1 Days Daily application of drug.
Number of cycle nim 2000 Number Number of cycle is applied as per requirement.
Amount of Immuno-supressive drug retention on the subsequent day of drug application drugimret 80% % of drug of the
preceding day
Retention of 20% drug of the previous
day drug amount
Immuno-supressive drug sensitivity of normal stem cells (Sg) dimsensSg 0 (Dose1)
0 (Dose2)
% of cells Assumed value
Immuno-supressive drug sensitivity of normal progenitor cells (P1g, P2g, P3g) dimsensP1g 0 (Dose1)
0 (Dose2)
% of cells Assumed value
dimsensP2g 0 (Dose1)
0 (Dose2)
dimsensP3g 0 (Dose1)
0 (Dose2)
Immuno-supressive drug sensitivity of normal matured cells (B1g, B2g, B3g) dimsensB1g 0 (Dose1)
0 (Dose2)
% of cells Assumed value
dimsensB2g 0.06 (Dose1)
0.06 (Dose2)
dimsensB3g 0 (Dose1)
0 (Dose2)
Immuno-supressive drug sensitivity of drug sensitive stem cells (Ss) dimsensSs 0 (Dose1)
0 (Dose2)
% of cells Assumed value
Immuno-supressive drug sensitivity of drug sensitive progenitor cells (P1s, P2s, P3s) dimsensP1s 0 (Dose1)
0 (Dose2)
dimsensP2s 0.006 (Dose1)
0.006 (Dose2)
dimsensP3s 0 (Dose1)
0 (Dose2)
% of cells Assumed value
Immuno-supressive drug sensitivity of drug sensitive cells (B1s, B2s,B3s) dimsensB1s 0 (Dose1)
0 (Dose2)
dimsensB2s 0.06 (Dose1)
0.06 (Dose2)
dimsensB3s 0 (Dose1)
0 (Dose2)
% of cells Assumed value
Immuno-supressive drug sensitivity of drug resistive stem cells (Sr) dimsensSr 0 (Dose1)
0 (Dose2)
Immuno-supressive drug sensitivity of drug resistive progenitor cells (P1r, P2r, P3r) dimsensP1r 0 (Dose1)
0 (Dose2)
dimsensP2r 0.0006 (Dose1)
0.0006 (Dose2)
% of cells Assumed value
dimsensP3r 0 (Dose1)
0 (Dose2)
Immuno-supressive drug sensitivity of drug resistive cells (B1r, B2r, B3r) dimsensB1r 0 (Dose1)
0 (Dose2)
dimsensB2r 0.006 (Dose1)
0.006 (Dose2)
% of cells Assumed value
dimsensB3r 0 (Dose1)
0 (Dose2)
Immuno-supressive drug sensitivity of doner stem cells (Se) dimsensSe 0 (Dose1)
0 (Dose2)
Immuno-supressive drug sensitivity of doner cells (P1e, P2e, P3e) dimsensP1e 0 (Dose1)
0 (Dose2)
% of cells Assumed value
dimsensP2e 0 (Dose1)
0 (Dose2)
dimsensP3e 0 (Dose1)
0 (Dose2)
Immuno-supressive drug sensitivity of doner cells (B1e, B2e, B3e) dimsensB1e 0 (Dose1)
0 (Dose2)
% of cells Assumed value
dimsensB2e 0.06 (Dose1)
0.06 (Dose2)
dimsensB3e 0 (Dose1)
0 (Dose2)
C. Immunosuppressive drug related parameters
Table 2: Initial parametric values for applied drugs
Variable Symbol Value Unit
Initial number Se 14 Number of cells
Multiplication rate mSe 0.065 Cells/day
Apoptosis rate aSe 0.032 Cells/day
Differentiation rate Sedr cf2*(cf1-(mSe + aSe )) Cells/day
Differentiation delay time Sedt 14 Days
A. Stem cell compartment
Variable Symbol Value Unit
Erythroid
Initial number P1e 0 Number of cells
Multiplication rate mPle 0.0546 Cells/day
Apoptosis rate aPle 0.011 Cells/day
Differentiation rate P1edr cf4*(cf1-(mP1e + aP1e)) Cells/day
Differentiation delay time P1edt 2 Days
Leukocytic
Initial number P2e 0 Number of cells
Multiplication rate mP2e 0.0545 Cells/day
Apoptosis rate aP2e 0.0312 Cells/day
Differentiation rate P2edr cf4*(cf1-(mP2e +aP2e )) Cells/day
Differentiation delay time P2edt 7 Days
Megakaryocytic
Initial number P3e 0 Number of cells
Multiplication rate mP3e 0.0535 Cells/day
Apoptosis rate aP3e 0.01 Cells/day
Differentiation rate P3edr cf4*(cf1-(mP3e + aP3e)) Cells/day
Differentiation delay time P3edt 3 Days
B. Progenitor cell compartment
Variable Symbol Value Unit
Erythrocytes
Initial number B1e 0 Cell numbers
Apoptosis rate aB1e 0.00007747 Cells/day
Leukocytes
Initial number B2e 0 Cell numbers
Apoptosis rate aB2e 0.012999996 Cells/day
Platelet
Initial number B3e 0 Cell numbers
Apoptosis rate aB3e 0.00017 Cells/day
C. Terminally mature cell compartment
Table 3: Initial (assumed) parametric values of the transplanted cell (Donor). Cell numbers are considered as count/μl of blood
Kill factor Target cells Value [% of cells killing/
day] in 100% HLA
matching
Value [% of cells killing/ day] in GVHD Comments
Low GVHD High GVHD
killb2e→hm Ss, P1s, P2s, P3s, B1s, B2s, B3s,
Sr, P1r, P2r, P3r, B1r, B2r, B3r
.
7.142 7.142×10 7.142×10 Three types of multiplying factors have
been considered
(i) ‘G’ has been considered to control the overall killing intensity. And
G = 0 for 100% matching
G = 18 for low GVHD
G = 19 for high GVHD
(ii) Compartment wise killing intensity are controlled by three more multiplying factors
st_mul = 0.02, pro_mul = 0.04 and mat_
mul = 0 1 for all cases considered.
(iii) In progenitor cell compartment for three types of cells killing intensity are
being controlled by three more multiplying factors –
effect_p1 = 10, effect_p2 = 0.01, effect_p3
= 0.3 for all cases considered.
(iv) Killb2g→RBC and Killb2e→RBC are kept zero
as GVHD may occur in liver, gut and
lungs without any hemolysis.
killb2e→g Sg, P1g, P2g, P3g, B2g, B3g. 0 0.02×5 0.02×10
killb2g→D Se, P1e, P2e, P3e, B1e, B2e, B3e. 0 0.03×5 0.03×10
killb2e→RBC B1g 0 0 0
killb2g→RBC B1e 0 0 0
hm = malignant cells within host (recipient), D = donor cells. All values are calibrated data
Table 4: Different kill factors description for three considered cases (i.e., 100% HLA matching, low GVHD and high GVHD)
Variables Day of Observation Strategy I Free
Growth
Strategy II MYL Strategy III MYL+ Transfusion Strategy IV Chemo Strategy V
Chemo+
Transfusion
Strategy VI
MYL+
Transfusion+
cytokine
Strategy VII Chemo + Transfusion +Cytokine Strategy VIII
100%
HLA
Strategy IX
Low
GVHD
Strategy X
Low
GVHD+
IMS1
Strategy XI
Low
GVHD+
IMS(1&2)
Strategy XII
High
GVHD
Strategy XIII
High
GVHD+
IMS1
Strategy XIV
High
GVHD+
IMS(1&2)
Sg 100 35.0167 10.3242 10.6605 35.3867 35.3867 10.6605 35.3867 10.6605 10.6605 10.6605 10.6605 10.6605 10.6605 10.6605
300 34.9587 10.7536 10.5485 35.6752 35.5680 16.5596 35.0343 10.5778 3.1079 3.2348 3.4121 0 on d282 0 on d282 0 on d283
500 0 on d 410 34.3735 34.7197 33.2508 33.0376 35.4694 35.4189 12.8343 0 on d319 0 on d320 0 on d322  -------  -------  -------
1000 - 0 on d 570 0 on d 574 0 on d 533 0 on d 533 0 on d 574 35.1821 12.1834  -------  -------  -------  -------  -------  -------
Sr 100 8.1256 1.8143 1.9023 3.8662 3.8662 1.9023 3.8662 1.9023 1.9023 1.9023 1.9023 1.9023 1.9023 1.9023
300 559.6793 49.0347 47.6234 167.5038 166.6943 47.6234 56.8266 20.2441 0 on d224 0 on d224 0 on d224 0 on d217 0 on d217 0 on d217
500 3.8550×104 3.3775×103 3.2802×103 1.1537×104 1.1482×104 3.2802×103 0 on d-412 3.9821  --------  --------  --------  --------  --------  --------
1000 1.5179×109 1.3299×108 1.2916×108 4.5428×108 4.5209×108 1.2916×108  ------- 0 on d504  --------  --------  --------  --------  --------  --------
Ss 100 5.8425 0.7649 0.8155 2.1411 2.1411 0.8155 2.1411 0.8155 0.8155 0.8155 0.8155 0.8155 0.8155 0.8155
300 50.9499 1.8765 1.8041 9.9795 9.9146 1.8041 0 on d-175 0 on d227 0 on d209 0 on d209 0 on d209 0 on d207 0 on d207 0 on d207
500 444.3107 16.3642 15.7327 87.0271 86.4609 15.7327  -------  -------  -------  -------  -------  -------  -------  -------
1000 9.9781×104 3.6750×103 3.5332×103 1.9544×104 1.9417×104 3.5332×103  --------  --------  --------  --------  --------  --------  --------  --------
P1g 100 485.8812 33.5282 37.1630 189.3582 189.3582 37.1630 189.3582 37.1630 37.1630 37.1630 37.1630 37.1630 37.1630 37.1630
300 339.9653 49.0031 44.3867 494.6528 488.1867 44.2251 488.1947 44.3874 0 on d246 0 on d246 0 on d246 0 on d236 0 on d236 0 on d236
500 0 on d 321 56.5098 81.1015 0 on d 452 0 on d-453 92.4633 503.7698 379.0082  --------  --------  --------  --------  --------  --------
1000 .--------- 0 on d-504 0 on d506  ---------  --------- 0 on d-506 499.6429 490.9414  --------  --------  --------  --------  --------  --------
P1r 100 5.5453 0.4792 0.4904 1.5004 1.5004 0.4904 1.5004 0.4904 0.4904 0.4904 0.4904 0.4904 0.4904  
300 180.4301 15.6180 15.1194 53.4973 53.2062 15.1194 32.3269 0on d 245 0 on d211 0 on d211 0 on d211 0 on d208 0 on d208 0 on d208
500 1.2344×104 1.0814×103 1.0503×103 3.6943×103 3.6764×103 1.0503×103 45.8008  ---------  ---------  ---------  ---------  ---------  ---------  ---------
1000 4.8603×104 4.2583×107 4.1357×107 1.4546×108 1.4476×108 4.1357×107 5.0563  ---------  ---------  ---------  ---------  ---------  ---------  ---------
P1s 100 11.9047 0.6700 0.7031 2.6453 2.6453 0.7031 2.6453 0.7031 0.7031 0.7031 0.7031 0.7031c 0.7031 0.7031
300 182.2828 15.5669 15.0704 53.4493 53.1580 15.0704 32.3285 0 on d 245 0 on d211 0 on d211 0 on d211 0 on d208 0 on d208 0 on d208
500 1.2297×104 1.0772×103 1.0462×103 3.6798×103 3.6620×103 1.0462×103 45.0584  --------  --------  --------  --------  --------  --------  --------
1000 4.8412×108 4.2415×107 4.1194×107 1.4489×108 1.4419×108 4.1194×107 4.7256  --------  --------  --------  --------  --------  --------  --------
P2g 100 95.6186 5.6304 6.1012 20.4474 20.4474 6.1012 20.4474 6.1012 6.1012 6.1012 6.1012 6.1012 6.1012 6.1012
300 68.8196 29.5808 25.5780 85.3616 83.6121 42.8545 159.0237 25.5787 24.7253 24.7461 24.7731 17.3686 17.4062 17.7578
500 0 on d 305 0 on d-488 0 on d491 0 on d-432 0 on d-433 94.0281 298.9938 151.5945 0 on d484 0 on d493 0 on d514 0 on d410 0 on d410 0 on d422
1000  ---------  ---------  ---------  ---------  --------- 0 on d-506 297.0921 296.0064  ---------  ---------  ---------  ---------  ---------  ---------
P2r 100 107.8105 0.9720 1.0723 9.5575 9.5575 1.0723 9.5575 1.0723 1.0723 1.0723 1.0723  -------- 1.0723  
300 1.2801×104 92.4047 80.2331 529.7196 518.1488 80.2331 0 on d-177 11.8146 0 on d216 0 on d216 0 on d216 0 on d211 0 on d211 0 on d211
500 1.4629×106 1.6885×104 1.5370×104 8.0828×104 7.9463×104 1.5370×104  --------- 0 on d343  -------  -------  -------  ------- ------- -------
1000 1.8926×1011 3.0477×109 2.8349×109 1.3226×1010 1.3044×1010 2.8349×109  ---------  --------  --------  --------  --------  --------  --------  --------
P2s 100 193.1684 0.3969 0.4125 6.9023 6.9023 0.4125 6.9023 0.4125 0.4125 0.4125 0.4125 0.4125 0.4125 0.4125
300 3.2984×103 38.1150 34.9882 144.2417 141.5614 34.9882 0 on d-174 0 on d266 0 on d213 0 on d213 0 on d213 0 on d209 0 on d209 0 on d209
500 8.0273×104 3.3499×103 3.2236×103 1.1649×104 1.1564×104 3.2236×103  --------  -------  -------  -------  -------  -------  -------  -------
1000 1.6266×109 1.3947×108 1.3543×108 4.7659×108 4.7426×108 1.3543×108  --------  -------  -------  -------  -------  -------  -------  -------
P3g 100 398.6151 328.9891 307.4990 406.2133 406.2133 307.4990 406.2133 307.4990 307.4990 307.4990 307.4990 307.4990 307.4990 307.4990
300 386.8417 406.2747 413.3843 407.0047 403.8850 413.2424 403.8922 391.7407 401.8555 402.2933 402.9248 225.8939 226.5821 246.2507
500 0 on d 343 392.7604 386.5695 0 on d-470 0 on d-472 388.9702 408.7169 381.8802 0 on d389 0 on d391 0 on d395 0 on d352 0 on d353 0 on d357
1000      ----- 0 on d535 0 on d537  ------------  ------------ 0 on d-537 409.8584 349.6556  --------  --------  --------  --------  --------  --------
P3r 100 2.8890 0.3525 0.3515 0.9369 0.9369 0.3515 0.9369 0.3515 0.3515 0.3515 0.3515 0.3515 0.3515 0.3515
300 127.2011 11.1339 10.8074 38.0411 37.8529 10.8074 21.1795 0 on d238 0 on d210 0 on d210 0 on d210 0 on d208 0 on d208 0 on d208
500 8.7579×103 767.3020 745.2169 2.6211×103 2.6085×103 745.2169 8.5363  -------  -------  -------  -------  -------  -------  -------
1000 3.4484×108 3.0212×107 2.9343×107 1.0321×108 1.0271×108 2.9343×107 0.0046  -------  -------  -------  -------  -------  -------  -------
P3s 100 2.5471 0.3334 0.3307 0.8599 0.8599 0.3307 0.8599 0.3307 0.3307 0.3307 0.3307 0.3307 0.3307 0.3307
300 119.1560 10.4342 10.1304 35.647 35.4721 10.1304 19.5580 0 on d 237 0 on d210 0 on d210 0 on d210 0 on d207 0 on d207 0 on d207
500 8.2057×103 718.9148 698.2226 2.4558×103 2.4440×103 698.2226 5.9997  --------  --------  --------  --------  --------  --------  --------
1000 3.2309×108 2.8307×107 2.7492×107 9.6697×107 9.6230×107 2.7492×107 9.2519×10-4  --------  --------  --------  --------  --------  --------  --------

† d = day
A. Progenitor cells

Variables Day of Observation Strategy I Free
Growth
Strategy II MYL Strategy III MYL+ Transfusion Strategy IV Chemo Strategy V
Chemo+
Transfusion
Strategy VI
MYL+
Transfusion+
cytokine
Strategy VII Chemo + Transfusion +Cytokine Strategy VIII
100%
HLA
Strategy IX
Low
GVHD
Strategy X
Low
GVHD+
IMS1
Strategy XI
Low
GVHD+
IMS(1&2)
Strategy XII
High
GVHD
Strategy XIII
High
GVHD+
IMS1
Strategy XIV
High
GVHD+
IMS(1&2)
Sg 100 4.7492×106 1.4524×105 1.6188×105 8.6256×105 8.6256×105 1.6188×105 8.6256×105 1.6188×105 1.6188×105 1.6188×105 1.6188×105 1.6188×105 1.6188×105 1.6188×105
300 4.7467×106 1.6504×105 1.5445×104 2.9753×105 2.9427×105 1.5445×104 2.9434×105 1.5452×104 1.5453×104 1.5453×104 1.5453×104 1.5453×104 1.5453×104 1.5453×104
500 4.5616×106 1.4360×104 1.3491×104 2.8748×105 2.8439×105 1.3492×104 2.9463×105 1.5468×104 1.5453×104 1.5453×104 1.5453×104 1.5453×104 1.5453×104 1.5453×104
1000 0 on d639 0 on d583 0 on d583 0 on d639 0 on d639 0 on d583 2.9537×105 1.6165×104 1.5452×104 1.5452×104 1.5452×104 1.5452×104 1.5452×104 1.5452×104
Blr 100 6.2473×104 1.9104×103 2.1294×103 1.1346×104 1.1346×104 2.1294×103 1.1346×104 2.1294×103 2.1294×103 2.1294×103 2.1294×103 2.1294×103 2.1294×103 2.1294×103
300 6.2460×104 220.7006 206.5868 3.9233×103 3.8804×103 206.5868 3.8776×103 142.9936 0 on d245 0 on d245 0 on d245 0 on d233 0 on d233 0 on d233
500 6.5419×104 484.8125 463.0688 4.8226×103 4.7753×103 463.0688 3.9329×103 0 on d371  --------  --------  --------  --------  --------  --------
1000 1.2063×108 1.0563×107 1.0259×107 3.6088×107 3.5913×107 1.0259×107 3.9739×103  --------  --------  --------  --------  --------  --------  --------
Bls 100 1.8740×105 5.7307×103 6.3873×103 3.4035×104 3.4035×104 6.3873×103 3.4035×104 6.3873×103 6.3873×103 6.3873×103 6.3873×103 6.3873×103 6.3873×103 6.3873×103
300 1.8724×105 654.5210 612.5553 1.1741×104 1.1612×104 612.5553 1.1609×104 548.9708 0 on d272 0 on d272 0 on d272 0 on d252 0 on d252 0 on d252
500 1.9003×105 917.0708 867.5400 1.2627×104 1.2494×104 867.5400 1.1655×104 0 on d479  ---------  -------  -------  -------  -------  -------
1000 1.2028×108 1.0523×107 1.0220×107 3.5955×107 3.5781×107 1.0220×107 1.1670×104  --------  --------  --------  --------  --------  --------  --------
B2g 100 3.4467×103 113.6635 126.5753 658.9195 658.9195 126.5753 658.9195 126.5753 126.5753 126.5753 126.5753 126.5753 126.5753 126.5753
300 880.3321 6.1822 6.3190 183.0205 181.4611 9.4576 260.9238 13.0931 0 on d269 0 on d268 0 on d268 0 on d248 0 on d248 0 on d248
500 0 on d317 0 on d321 0 on d322 0 on d355 0 on d356 0 on d538 516.6665 78.3325  --------  --------  --------  --------  --------  --------
1000   ------       ------       ------      ---------  ---------   ------     1.1877×103 736.3299  --------  --------  --------  --------  --------  --------
B2r 100 905.2402 4.3430 5.1049 68.3978 68.3978 5.1049 68.3978 5.1049 5.1049 5.1049 5.1049 5.1049 5.1049 5.1049
300 2.9970×103 13.3488 11.2598 94.9068 92.7326 11.2598 0 on d174 0 on d212 0 on d205 0 on d205 0 on d205 0 on d204 0 on d204 0 on d204
500 2.6537×105 2.8315×103 2.5618×103 1.3918×104 1.3672×104 2.5618×103  -------  ---------  ---------  ---------  ---------  ---------  ---------  ---------
1000 3.4525×1010 5.4765×108 5.0899×108 2.3861×109 2.3528×109 5.0899×108  --------  ---------  ---------  ---------  ---------  ---------  ---------  ---------
B2s 100 4.4867×103 3.0000 3.7334 137.2267 137.2267 3.7334 137.2267 3.7334 3.7334 3.7334 3.7334 3.7334 3.7334 3.7334
300 4.5007×103 6.3631 5.5790 38.3634 37.3888 5.5790 0 on d162 0 on d209 0 on d204 0 on d204 0 on d204 0 on d204 0 on d204 0 on d204
500 2.4313×104 667.6654 638.6114 2.3586×103 2.3374×103 638.6114  -------  --------  --------  --------  --------  --------  --------  --------
1000 3.4020×108 2.8796×107 2.7959×107 9.8424×107 9.7941×107 2.7959×107  -------  --------  --------  --------  --------  --------  --------  --------
B3g 100 4.4945×105 1.3788×104 1.5363×104 8.1702×104 8.1702×104 1.5363×104 8.1702×104 1.5363×104 1.536×104 1.5363×104 1.5363×104 1.5363×104   1.5363×104
300 4.4676×105 1.6707×103 1.5652×103 2.8217×104 2.7909×104 1.5652×103 2.7976×104 1.5692×103 1.5399×103 1.5405×103 1.5414×103 1.5244×103   1.5296×103
500 2.6183×105 0 on d493 0 on d494 1.8145×104 1.8013×104 0 on d494 2.8065×104 1.6567×103 1.2043×103 1.2274×103 1.2729×103 692.5716   832.5999
1000 0 on d538  ---------  --------- 0 on d542 0 on d542  -------- 2.8285×104 1.8623×103 0 on d787 0 on d827 0 on d827 0 on d595 0 on d611 0 on d648
B3r 100 2.4964×104 763.3894 850.8658 4.5339×103 4.5339×103 850.8658 4.5339×103 850.8658 850.8658 850.8658 850.8658 850.8658   850.8658
300 2.4915×104 88.7642 83.1135 1.5670×103 1.5499×103 83.1135 1.5480×103 20.6596 0 on d231 0 on d231 0 on d231 0 on d223 0 on d223 0 on d223
500 2.6628×104 245.1372 234.9932 2.0974×103 2.0777×103 234.9932 1.5663×103 0 on d315  -------  -------  -------  -------   -------
1000 7.1457×107 6.2584×106 6.0783×106 2.1380×107 2.1277×107 6.0783×106 1.5572×103  --------  --------  --------  --------  --------    --------
B3s 100 2.4962×104 763.3145 850.7824 4.5334×103 4.5334×103 850.7824 4.5334×103 850.7824 850.7824 850.7824 850.7824 850.7824   850.7824
300 2.4907×104 88.6082 82.9666 1.5661×103 1.5490×103 82.9666 1.5472×103 20.6143 0 on d231 0 on d231 0 on d231 0 on d223 0 on d223 0 on d223
500 2.6505×104 235.3585 225.5020 2.0634×103 2.0439×103 225.5020 1.5623×103 0 on d315  -------  -------  -------  -------    -------
1000 6.7075×107 5.8745×106 5.7054×106 2.0069×107 1.9972×107 5.7054×106 1.5515×103  -------  -------  -------  -------  -------    -------
† d = day
B. Mature cells
Table 5: Simulation results with different strategies

Partnered Content Networks

  • Cancer Science
  • Vaccine Studies
  • Gynecology
  • Food Nutrition
  • Nursing Science
  • Public Health
  • The Pharma
  • Infectious Disease
  • Neuro Care
  • Catalysis
  • Neonatal Biology
  • Neonatal Disorders
  • Mutation
  • Nanotechnology
  • Toxicology
  • Dark Biotechnology
  • Pollution Toxicology
  • Cell Biology
  • Bioanalytical Research
  • Renal Disorders
  • The Astrophysics
  • Sleep Physiology
  • Epidemiology
  • Histology