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

Volume 1 | Issue 1 Development of the Analytical Model for the Assessment of the Efficiencies of Different Therapeutic Modalities in Leukemia Dhar PK1,2 and Majumder D*1,3 1Society for Systems Biology and Translational Research, Kolkata, India 2Department of Electronics and Communication Engineering, Bengal College of Engineering and Technology, Sahid Sukumar Banerjee Sarani, Bidhannagar, Durgapur, Burdwan, West Bengal, India 3Department of Physiology, West Bengal State University, Berunanpukuria, Malikapur, Barasat, Kolkata, India *Corresponding author: Majumder D, Department of Physiology, West Bengal State University, Berunanpukuria, Malikapur, Barasat, Kolkata 700 126, India, E-mail: durjoy@ssbtr.net, durjoy@rocketmail.com Research Article Open Access


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][2][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].[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.
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 (S g , S s and S r ), in second compartment total nine cell types (P1 g , P2 g , P3 g , P1 s , P2 s , P3 s , P1 r , P2 r and P3 r ) and in third compartment total nine cell types (B1 g , B2 g , B3 g , B1 s , B2 s , B3 s , B1 r , B2 r and B3 r ) are considered.Therefore in present system model we have considered total 21 types of cells (Figure 1).

Methods
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) N x (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.N xp (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 ( ) x N m and decays with its own apoptosis rate ( ) x N a .Thus, for normal stem cell these variables are denoted with ( ) g S m and ( ) g S a .Similarly for drug sensitive P1 type cell, these variables are represented by ( ) s P m 1 and ( ) s P a 1 .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.
Again N x dr is the differentiation rate of any cell type and N xp dr is the differentiation rate of its concerned precursor cell type.For example, in the system equation of stem cells has only N x dr and for sensitive category stem cell, N x dr is represented as S s dr.
Similarly, for drug resistive category of P1 lineage cell type, N x dr = P1 r dr and N xp dr = S r dr; for P3 lineage cell type of drug sensitive category, N x dr = P3 s dr and N xp dr = S s dr 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: -------------------------------------------------------------------- (4) and if k≥ dk 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.

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][21][22].This has been incorporated by modifying equation (7a) into equation (7b): where the functional effective (immunity related) killing of the malignant cells (kill eff ) by B2 g becomes effective depending upon the minimum cell number of B2 g ( Min 2 ) and is operated by kill im and CF.The assumed efficiency of immune related killing (kill eff ) by the host immune cells are in the order of B2 s >B2 r =P2 s >P2 r >S s >S r .It is assumed that initially there is total (functional) immune suppression i.e., kill im = 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.Where i = j = 1, 2, 3, ....21 ( ) 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: In equation (7c) it is assumed that Cancer cachexia 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.

Myeloablative chemotherapy model
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: , , 2 , , 1 where, d myl (k) is the amount of drug present within the system denoted either the amount of applied drug (drug myldose ) on the day of drug application or drug retention (drug mylret ) 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);  ---------------(7d) .
Drug application is started on the day k mylsd and n myl is the number of cycles and t myld is the interval of two subsequent drug applications in days and C myl is the ON/OFF switch for drug application.Myl 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 ----- (8) where where 13 and 1 0 else 1 when is the day of MYL application and 1 1 1 when is the day of MYL application, 3 1 1 when is the day of MYL application and 1 1 1 when is the day of MYL application, 3 1 1 when is the day of MYL application and 1 1 1 when is the day of MYL application, 3 1 1 when is the day of MYL application and 1 1 1 when is the day of MYL application, 3 1 1 when is the day of MYL application and 1 1 1 when is the day of MYL application, 3 1 when is the day of MYL application and 3 1 T RBC and T PLATELET represent amount of transfused RBC and platelets respectively.E1 and E2 may be fraction or whole number.

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, ---------(8a) ( )  for , 2, 1 ,

Host cells' behavior towards transplanted cells: Model has been formulated with an assumption that donor stem cell (S e )
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 (B2 e ) are developed [along with other lineages i.e., erythrocytes (B1 e ) and platelets (B3 e )] from transplanted stem cells lead to an immune response causing Graft-versus-tumor (GVT) (or GVL) [kill b2e→hm ] and/or GVHD [kill b2e→hg ] effect.So a (i,j) represented in equation ( 8a) is modified further as: In equation (8a-1) hm stands for host leukemic cells (S s , S r , P1 s , P1 r , P2 s , P2 r , P3 s , P3 r , B1 s , B1 r , B2 s , B2 r , B3 s , B3 r ) and hg stands host normal cells (i.e., S g , P1 g , P2 g , P3 g , B1 g , B2 g and B3 g ).The GVHD reaction is operated into the system by modifying the kill factors kill b2g→D and kill b2e→hg .The killing intensity kill mul is operative for GVHD and GVL effect that can be different for different cell types of host.
- --------------(8a-1) 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 (kill b2e→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 B2 e is developed from transplanted HSC.As a result with the progress of time GVHD may be operative nonlinearly depending on the cell numbers of B2 e .Killing efficiency of B2 e to other hematopoietic cells of donor will be governed with a multiplying factor effect procell which is governed though the following conditions.for , where 3 6 9 12 15 18 21   G kill kill i j 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 (kill mu1 ) 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 Cal GVHD is a time varying calibration factor.
- --------------(8a-2) 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][30][31][32][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.

Cytokine therapy model
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: where 1 2 3 .21 1 1 where 22, 23, 24, .... 28 for , where 2, 3, 8, 9, 17, 18 for 16 1  C CYT and k CYTsd , n CYT , t CYTd , m CYT and C CYT 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 m DCYT.

Maintenance of transplantation by immunosuppressive drug
To keep the chronic GVHD under control immunosuppressive drugs (d im ) 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., P2 g , P2 r , P2 s , B2 g , B2 r , B2 s , P2 e and B2 e .Depending upon the severity of GVHD /matching, clinical decision regarding the immunosuppressant drug (d im ) 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.

CYTdose
where 1 2 3 .21(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 (drug imdose )or drug retention (drug imret ) 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 k imsd , n im is the number of cycles, t imd is the interval of two subsequent drug applications in days and C im is the ON/OFF switch for drug application is introduced through above scheme.And d imsensh and d imsensD 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 (  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 S g 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 30 th 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][36][37][38][39][40] i.e., when there is no exogenous RBC and/or platelet transfusion, T ex =0.

High dose Myeloablative chemotherapy application:
System was simulated with different combinations of drug doses (Drug myldose ), cycle numbers (n myl ), intervals (t myld ) with the same initial conditions as depicted in Table 1 and followings are noted.Combination I: Drug myldose =0.7, n myl =8, t myld =15, T ex =0 (Figure 10A) Combination II: Drug myldose =0.8, n myl =8, t myld =15, T ex =0 Combination III: Drug myldose =0.7, n myl =9, t myld =15, T ex =0 Combination IV: Drug myldose =0.7, n myl =8, t myld =20, T ex =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 (Drug myldose ), cycle numbers (n myl ) and intervals (t myld ).
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.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., T ex =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 RBC llm and Platelet llm 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 a T_RBC (= 0.07747) and a T_Platelet (=0.017) respectively.
System was simulated in presence of RBC and platelet transfusions with Drug myldose =0.7, n myl =8 and t myld =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.For the treatment of leukemia through immunotherapy, particularly, cytokine based therapies are suggested after MYL therapy.It is assumed that applied cytokine affects only P2 g cell types (through differentiation process) and increases B2 g 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.

Cytokine therapy after conventional MYL/chemotherapy
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 (Drug myldose =0.Observation I: The system was simulated with different drug (Drug CYTdose ) level when n CYT =15, t CYTd =25, g B Min 2 =224 and eff kill =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 Drug CYTdose 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 eff kill when Drug CYTdose =1.2, n CYT =15, t CYTd =25 and g B Min 2 =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 eff kill 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 ( =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 g B Min 2 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 (Drug CYTdose = 1) after chemotherapeutic strategy with transfusion though there is a persistence of B1 r and B1 s cells having a number >2000 on day 1500 however with a decaying tendency.Application of high dose MYL (Drug myldose =0.7, n myl =8, t myld =15, T ex =1) followed by cytokine drug application: With parametric variations the followings observations are noted.
The system was simulated with different drug (Drug CYTdose ) levels when n CYT =15, t CYTd =25, g B Min 2 =224 and eff kill =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 Drug CYTdose 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 338 th, 304 th and 291 st when Drug CYTdose level were assumed as 1, 4 and 6 respectively).For Drug CYTdose =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 Drug CYTdose = 6 but it may be unrealistic (as Drug CYTdose ×m CYT =6×0.07=0.42i.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 P2 g <30, the S g dt and P2 g dt are increased to 30 days and 20 days respectively than their normal value; whereas S g dr and P2 g dr are decreased to 0.0165×{3-( For simulation it is assumed that HSC transplantation has been done 25 days (on 180 th 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 30 th , 45 th , 60 th , 75 th , 95 th , 115 th , 135 th , 155 th .The unequal gap between successive drug application were due to transfusion of RBC and Platelet on 90 th , 110 th , 130 th and 150 th 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 (a) (b) =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 (S e = 12) having 2 b e hm kill → = 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.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 (kill b2e→hm ) (same as for 100% HLA matching case) and host normal cells (kill b2e→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(B2 e , 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 S r reduced to zero on day 504 th when transplanted with 100% HLA matched donor cell; however, when transplanted with low and high mismatched donor, removal of S r is on day 224 th and 217 th respectively.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 2: As
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.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.

Application of immunosuppressive drug scheduling for HLA mismatch cases
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.
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.

Remission by chemotherapy /chemo-responsive leukemia
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 3 times and 1.5 times of S r dr, P2 s dr, and P2 r dr 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.

(B)
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' 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 S r dr, P2 s dr and P2 r dr 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][42][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.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 m DCYT.
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.
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 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.

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.
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.

Figure 1 :
Figure 1: Different cellular components of model

Figure 2 :
Figure 2: Considered feedback system for maintenance of normal hematopoietic homeostasis

Figure 3 :
Figure 3: Cachexia effect by the developed malignant cells 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 levelwith 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 P2 r (k) +P2 s (k) <10,000.

Figure 4 :
Figure 4: Scheme for application of Myeloablative (MYL) therapy myeloablative drug sensitivity of the concerned cell x N .The drug sensitivity x N drug sensitivities of T RBC and T PLATELET respectively.Application strategy of blood transfusion is demonstrated in Figure5.

Figure 7 :
Figure 7: Scheme for application of Myloablative (Myl)/ chemotherapy along with supportive (RBC/platelet transfusion) therapy followed by cytokine therapy set GVHD ), drug dose and drug type by affecting GVHD Cal with a fractional multiplier ( res GVHD ).The operational relationship is represented as follows: Application strategy of immunosuppressive after MYL/Chemo along with transplantation is demonstrated in Figure 8.

Figure 8 :
Figure 8: Algorithmic scheme for MYL therapy followed by HLA matched/mismatched HSC transplantation and immunosuppressive application Simulation Results

Figure 9 :
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' (b) Annex Publishers | www.annexpublishers.comVolume 1 | Issue 1

Figure 11 :
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 (n myl ) 8, drug dose (Drug myldose ) 0.7 (MYL) /0.35 (Chemo), interval (t myld ) 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' related killing (i.e., to turn 1 = im kill ) in presence of cytokine drug and malignancy killing efficiency ( eff kill ) 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 ( g B Min 2 ) and malignancy killing efficiency ( eff kill ) by cytokine mediated killing.

Figure 12 :
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'

2 )
to start immunity related killing when Drug CYTdose =1, n CYT =15, t CYTd =25 and eff kill =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 2 g B Min was 223, 222 and 224 respectively; i.e., with the lowering the minimum effective WBC level ( g B Min 2 ) 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 g B Min 2

2 )
} cells/day respectively) with respect to their normal value and g S a gets increased (to 0.48 cells/day) with respect to its normal value.

Figure 13 :Observation 1 :Observation 2 :
Figure 13: Effect of allogenic HSC transplantation with 100% HLA matching with number of donor stem cell (S e ) = 11, 2 b e hm kill → = 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'

Observation 3 :Observation 4 :
With other conditions unchanged when simulation is carried out with S e = 12 and 2 b e hm kill → = 7.125; however, transplantation is assumed to be on 190 th day after the day of diagnosis, malignancy exist, in all level.With other conditions unchanged donor stem cells with 2 b e hmkill → = 7.125 and transplantation is assumed to be on 190th day after the day of diagnosis with more number of HSC (S e = 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.
kill → increases while other factors unchanged normal cell population also reduces faster (Table

Figure 14 :
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'
drug sensitivity of normal matured cells (B1 g , B2 g , B3 g ) drug sensitivity of drug sensitive cells (B1 s , B2 s , drug sensitivity of drug resistive progenitor cells (P1 r , P2 r , P3 r ) Stem cell compartment Immuno-supressive drug sensitivity of drug resistive cells (B1 r , B2 r , drug sensitivity of doner cells (P1 e , P2 e , P3 e ) Immunosuppressive drug related parameters
RBC transfuse and PLATELET transfuse are two (21×1) matrix defined as and On the day of MYL application, minimal level (a set value) of RBC (RBC llm ) and platelets (Platelet llm ) 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 RBC llm is less then Platelet ulm or if Platelet llm is less then RBC ulm ) 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 k mylsd 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 - 7, n myl =8, t myld =15, T ex =1) and low dose chemotherapy (Drug myldose =0.35, n myl =8, t myld =15, T ex =1).The initial parametric values including the assumed effects of cytokine drug application are presented in Table2.It was observed that cytokine drug dose (Drug CYTdose ), drug application cycle (n CYT ), interval between two successive drug applications (t CYTd ), minimum number of WBC (

Table 1 :
Initial parametric values of the host hematopoietic system

Table 2 :
Initial parametric values for applied drugsContd.

Table 4 :
Different kill factors description for three considered cases (i.e., 100% HLA matching, low GVHD and high GVHD)

Table 3 :
Initial (assumed) parametric values of the transplanted cell (Donor).Cell numbers are considered as count/μl of blood Contd.