Received date: September 28, 2015; Accepted date: November 26, 2015; Published date: December 02, 2015
Citation: Dhar PK, Naskar TK, Majumder D (2015) Analytical Model for the Assessment of Efficiency of Stem Cell Transplantation with Suicidal Gene Construct for the Treatment of Leukemia. Oncol Trans Res 1:103. doi:10.4172/jotr.1000103
Copyright: © 2015 Dhar PK, et al. 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.
Visit for more related articles at Journal of Oncology Translational Research
Hematopoietic stem cell transplantation is now being the emergent methodology for leukemia treatment. Due to HLA mis-match, there is chance of transplantation related mortality. However, with the increase in HLA mis-match between donor cells and recipient, there is more chance of complete removal of leukemic cells in host (patient). To tackle this, recently stem cells transplantation with suicidal tk-gene construct is being suggested. Due to unavailability of suitable analytical methods this option has limited applications in clinical cases. Present work provides an analytical platform to test the efficacy of this therapeutic procedure.
Hematopoietic stem cell transplantation (HSCT) is considered to be the emerging therapeutic procedure for leukemia treatment. And hence a lot of clinical trials have been done across the globe. Generally allogenic transplantation is being made and the success rate depends on HLA (Human Leukocytic Antigen) incompatibility between the donor and receipient. The degree of HLA incompatibility enhances the killing efficiency of leukemic cell(s); contrarily it decreases the quality of life of the concerned patient due to enhanced GVHD (Graft Versus Host Disease). Therefore stem cell transplantation is being regarded as the ‘double edge sword’ [1,2]. Hence HSCT requires a judicious choice of donor cell selection. Though there is availability of rich medical and clinical literature on HSCT; however, most of the literatures are based on empirical observations and physicians personal experiences. This is due to unavailability of proper analytical method for such therapeutic procedure.
Though several analytical models are available for leukemia pathogenesis and dynamics, imatinib treatment, myeloablative chemotherapy, immune potentiation effect for different types of leukemia (CML, Chronic Myelogenous Leukemia; AML, Acute Myeloblastic Leukemia); but no single model was available to test the efficacy of different combinations of therapy for any particular leukemia type for any particular individual patient [3-15]. Recent time different therapeutic strategies along with their combinations can be assessed for different types of leukemia with a single analytical model . Through this model the efficacy of myeloablative and low dose chemotherapy along with supportive therapy like RBC and platelet transfusion, cytokine therapy, imatinib therapy, allogenic HSCT with different degree of HLA matching, effect of GVHD and immunosuppressive drug after HSCT can be tested.
In order to optimize the situation i.e., to enhance malignancy removal together with improved safety of allogenic transplantation, recently allogenic HSCT having HLA incompatibility with a tk (thymidine kinase) gene construct has been suggested. This therapeutic procedure has two way advantages – one, complete removal of leukemic cells(Graft Versus Leukemia reaction, GVL) due to high HLA mismatch and the other, tk gene construct provide the advantage of controlling GVHD by selective killing of transplanted cells with the application of pro-drug Gangcyclovir [17,18]. Though it seems to be safer, however such therapeutic procedure is a risky one as with the delay in application of pro-drug gangcyclovir, there is a high chance of transplanted related mortality. A suitable analytical method is required that can predict the minimal residual disease (MRD) and the optimal timing for application of pro-drug.
In the previous work , different therapies namely myeloablative chemotherapy (myl), cytokine (CYT), RBC (TRBC) and platelet (TPLATELET) transfusion (transfuse in subscript or tf), allogenic HSCT (Se), immuo-suppressive drug (im) in case of HLA incompatibility are being operative through following equation.
x(k)= Ax(k-1) +B1x(k-dk1) +B2x(k-dk2) +B3x(k-dk3) +B4x(k-dk4) + RBCtransfuseTRBC(k) + PLATELETtransfuseTPLATELET(k)
or, [x1(k); x2(k); …; x28(k)]=[a1,1, …, a1,28; …; a28,1, …, a28,28] ×[x1(k-1); x2(k-1); …; x28(k-1)] +[b11,1, …, b11,28; …; b128,1, …, b128,28] ×[x1(k-dk1); x1(k-dk1); … ; x28(k-dk1)] +[b21,1, …, b21,28; …; b128,1, …, b128,28] ×[x1(k-dk2); x1(k-dk2); … ; x28(k-dk2)] +[b31,1, …, b31,28; …; b328,1, …, b328,28] ×[x1(k-dk3); x2(k-dk3); …; x28(k-dk3)] +[b41,1, …, b41,28; …; b428,1, …, b428,28] ×[x1(kdk4); x2(k-dk4); …; x28(k-dk4)] +[RBCtf1,1; RBCtf2,1; …; RBCtf28,1] ×TRBC(k) +[PLATELETtf1,1; PLATELETtf2,1; …; PLATELETtf28,1] ×TPLATELET(k),
where the elements of matrix A are a1,1(k) = (1+mSg-aSg-Sgdr-CSgs-CSgr -Cmyl×MylSg×dmyl(k) -Cim×imsensSg×dim(k)), a2,2(k) = (1+mSs-aSs-Ssdr-CSsg- CSsr-Cmyl×MylSs×dmyl(k)-Cim ×imsensSs ×dim(k)), a3,3(k)=(1+mSr-aSr-Srdr- CSrs-CSrg -Cmyl×MylSr×dmyl(k) -Cim×imsensSr×dim(k)), a4,4(k)=(1+mP1g- aP1g-P1gdr-CP1gs-CP1gr -Cmyl×MylP1g×dmyl(k)-Cim×imsensP1g×dim(k)), a5,5(k)= (1 +mP1s -aP1s -P1sdr -CP1sg -CP1sr -Cmyl×MylP1s×dmyl(k) -Cim ×imsensP1s×dim(k)), a6,6(k)=(1 +mP1r -aP1r-P1rdr -CP1rg -CP1rs -Cmyl×MylP1r×dmyl(k)-Cim×imsensP1r×dim(k)), a7,7(k)=(1+mP2g-aP2g- P2gdr-CP2gs-CP2gr-Cmyl×MylP2g×dmyl(k)-Cim×imsensP2g ×dim(k) +CCYT×dCYT(k)×mCYT), a8,8(k)=(1+mP2s-aP2s -P2sdr -CP2sg -CP2sr -Cmyl×MylP2s×dmyl(k) -Cim×imsensP2s ×dim(k)), a9,9(k)=(1+mP2r-aP2r- P2rdr-CP2rg-CP2rs-Cmyl×MylP2r ×dmyl(k)-Cim×imsensP2r×dim(k)), a10,10(k)= (1+mP3g-aP3g-P3gdr-CP3gs -CP3gr-Cmyl×MylP3g×dmyl(k)- Cim×imsensP3g×dim(k)), a11,11(k)=(1 +mP3s -aP3s -P3sdr -CP3sg -CP3sr -Cmyl×MylP3s×dmyl(k) -Cim ×imsensP3s×dim(k)), a12,12(k) = (1 +mP3r -aP3r -P3rdr -CP3rg -CP3rs-Cmyl×MylP3r×dmyl(k)-Cim×imsensP3r×dim(k)), a13,13(k)= (1- aB1g -CB1gs -CB1gr -Cmyl×MylB1g×dmyl(k) -Cim×imsensB1g ×dim(k)), a14,14(k) = (1-aB1s-CB1sg-CB1sr -Cmyl×MylB1s×dmyl(k) -Cim×imsensB1s×dim(k)), a15,15(k) = (1 -aB1r -CB1rg -CB1rs -Cmyl ×MylB1r×dmyl(k)-Cim×imsensB1r×dim(k)), a16,16(k) =(1-aB2g-CB2gs -CB2gr -Cmyl×MylB2g×dmyl(k) -Cim×imsensB2g×dim(k)), a17,17(k) = (1 -aB2s -CB2sg -CB2sr -Cmyl×MylB2s×dmyl(k) -Cim×imsensB2s ×dim(k)), a18,18(k) = (1 -aB2r -CB2rg -CB2rs -Cmyl ×MylB2r×dmyl(k) -Cim×imsensB2r×dim(k)), a19,19(k) = (1 -aB3g -CB3gs -CB3gr -Cmyl×MylB3g×dmyl(k) -Cim×imsensB3g×dim(k)), a20,20(k) = (1 -aB3s -CB3sg -CB3sr -Cmyl×MylB3s×dmyl(k) -Cim ×imsensB3s×dim(k)), a21,21(k) = (1 -aB3r -CB3rg -CB3rs -Cmyl×MylB3r×dmyl(k) -Cim ×imsensB3r×dim(k)), a22,22(k)=(1 +mSe -aSe –Sedr -Cim ×imsensSe ×dim(k)), a23,23(k)=(1+mP1e-aP1e-P1edr -Cim ×imsensP1e×dim(k)), a24,24(k) = (1 +mP2e -aP2e -P2edr -Cim×imsensP2e×dim(k) +CCYT×dCYT(k)×mDCYT), a25,25(k)=(1 +mP3e-aP3e-P3edr-Cim×imsensP3e ×dim(k)), a26,26(k) = (1 -aB1e -Cim×imsensB1e×dim(k)), a27,27(k)=(1-aB2e -Cim×imsensB2e× dim(k)), a28,28(k)=(1-aB3e-Cim×imsensB3e×dim(k)). Multiplication rate, apoptosis rate and differentiation rate of concerned cell type are denoted by m, a and dr respectively.
Due to mutability one category cell is converted to another one, say normal (g) to resistant (r) or sensitive (s) category of different types of cells [stem cell (S), erythroblast (P1), lymphoblast (P2), megakaryoblast (P3), RBC (B1), WBC (B2), platelet (B3)]. The conversion rates (C) from one category (represented by first subscript) to another (represented by second subscript) are represented by the following marix elements.
a1,2 =CSsg, a1,3 =CSrg, a2,1 =CSgs, a2,3 =CSrs, a3,1 =CSgr, a3,2=CSsr, a4,5 =CP1sg, a4,6 = CP1rg, a5,4 = CP1gs, a5,6 =CP1rs, a6,4 =CP1gr, a6,5 = CP1sr, a7,8 = CP2sg, a7,9 = CP2rg, a8,7 =CP2gs, a8,9 =CP2rs, a9,7 =CP2gr, a9,8 =CP2sr, a10,11 = CP3sg, a10,12 = CP3rg, a11,10 = CP3gs, a11,12 =CP3rs, a12,10 =CP3gr, a12,11 = CP3sr, a13,14= CB1sg, a13,15 =CB1rg, a14,13 =CB1gs, a14,15 =CB1rs, a15,13 =CB1gr, a15,14 = CB1sr, a16,17 = CB2sg, a16,18 = CB2rg, a17,16 =CB2gs, a17,18 =CB2rs, a18,16 =CB2gr, a18,17 =CB2sr, a19,20 =CB3sg, a19,21 =CB3rg, a20,19 = CB3gs, a20,11 =CB3rs, a21,19 =CB3gr, a21,20 =CB3sr. In case of simulation study conversion rates were assumed to be zero.
Following matrix elements are representing the immunity related killing of malignant cells and are time varying functions depending on host lymphocyte (B2g) cell count.
a2,16(k)=-B2g_Ss(k), a3,16(k)=-B2g_Sr(k), a8,16(k)=-B2g_P2s(k), a9,16(k) = -B2g_P2r(k), a17,16(k) = -B2g_B2s(k)+ CB2gs, a18,16(k) = -B2g_ B2r(k) + CB2gr.
Following matrix elements are changed due to malignancy related toxicity development .
a1,2= -0.001+ CSsg, a1,3= -0.001+ CSrg, a4,5= -1+ CP1sg, a4,6= -1 +CP1rg, a10,11 = -1+ CP3sg, a10,12 = -1+ CP3rg, a22,2 = -0.001, a22,3 = -0.001, a23,5= -1, a23,6= -1, a25,11= -1, a25,12= -1.
Following matrix elements are changed due to cachexia related killing of host cells and exogenous transplanted (T) cells (T=1, after transplantation, otherwise T=0) by P2r development.
a1,9 = -0.000011, a4,9 = -0.000825, a7,9=-0.000825+CP2rg, a10,9= -0.000825, a13,9 = -0.00275, a16,9 = -0.000275, a19,9 = -0.000275, a22,9=T×0.000011, a23,9=-T×0.000825, a24,9=-T×0.000825, a25,9 = -T×0.000825, a26,9 = -T×0.000011, a27,9 = -T×0.000011, a28,9 =-T×0.000011.
Similarly, following matrix elements are changed due to cachexia related killing by P2s development.
a1,8= -0.000006, a4,8 = -0.00045, a7,8= -0.00045+CP2sg, a10,8= -0.00045, a13,8= -0.0015, a16,8= -0.0015, a19,8= -0.0015, a22,8 = T ×0.000006, a23,8=-T×0.00045, a24,8=-T×0.00045, a25,8=-T×0.00045, a26,8=-T×0.0015, a27,8=-T×0.0015, a28,8=-T×0.0015.
Similarly, following matrix elements are changed due to cachexia related killing by B2r development.
a1,18= -0.000003, a4,18= -0.0003375, a7,18= -0.0003375, a10,18= -0.0003375, a13,18= -0.00075, a16,18= -0.00075 +CB2rg, a19,18 = -0.00075, a22,18 =-T×0.000003, a23,18= -T×0.0003375, a24,18= -T×0.0003375, a25,18=-T×0.0003375, a26,18=-T×0.00075, a27,18= -T×0.00075, a28,18 = -T× 0.00075.
Likewise, following matrix elements are changed due to cachexia related killing by B2s development.
a1,17= -0.000003, a4,17 = -0.0003375, a7,17 = -0.0003375, a10,17= -0.0003375, a13,17= -0.00075, a16,17 = -0.00075+CB2sg, a19,17= - 0.00075, a22,17= -T×0.000003, a23,17 = -T×0.0003375, a24,17= -T×0.0003375, a25,17 = -T×0.0003375, a26,17= -T×0.00075, a27,17= -T×0.00075, a28,17= -T×0.00075.
Following matrix elements are changed due to GVHD related killing of donor cell by host lymphocytes.
a22,16 = -T ×G ×st_mul ×killb2g→D, a23,16= -T ×G ×pro_mul × effect_ p1e×killb2g→D, a24,16= -T ×G ×pro_mul ×effect_p2e × killb2g→D, a25,16= -T ×G ×pro_mul×effect_p3e ×killb2g→D, a26,16 =-T×G×mat_mul×killb2g→RBC, a27,16 =-T×G×mat_mul×killb2g→D, a28,16= -T×G×mat_mul×killb2g→D. The degree of HLA mis-match is represented through G.
Following matrix elements are changed due to GVHD related killing of host normal cells by donor lymphocytes. a1,27(k)= -T×G×st_mul×killb2e→g, a4,27(k) = -T×G×pro_mul× killb2e→g × GVHD_kill_p1g(k-1), a7,27(k) = -T×G×pro_mul× killb2e→g×GVHD_kill_ p2g(k-1), a10,27(k) = -T×G×pro_mul× killb2e→g × GVHD_kill_p3g(k-1), a13,27 = - T×G×mat_mul× killb2e→RBC, a16,27 = -T×G×mat_mul×killb2e→g, a19,27 = -T× G×mat_mul×killb2e→g.
Following matrix elements are changed due to GVL related killing of malignant cells by donor lymphocytes.
a2,27=-T×st_mul×killb2e→hm, a3,27=-T×st_mul×killb2e→hm, a5,27=- T×pro_mul×killb2e→hm, a6,27=-T×pro_mul×killb2e→hm, a8,27=-T× pro_mul × killb2e→hm, a9,27 = - T×pro_mul×killb2e→hm, a11,27=-T× pro_mul×killb2e→hm, a12,27 =-T ×pro_mul ×killb2e→hm, a14,27 =-T× mat_mul×killb2e→hm, a15,27 = -T×mat_mul×killb2e→hm, a17,27 = -T ×mat_mul×killb2e→hm, a18,27 =-T×mat_mul×killb2e→hm, a20,27 = -T×mat_mul×killb2e→hm, a21,27 = -T×mat_ mul×killb2e→hm. The parametric values of st_mul, pro_mul, mat_mul, killb2g→D, killb2g→RBC, killb2e→g, killb2e→RBC, killb2e→g, killb2e→hm are indicated in . The rest of the elements of matrix A are zero.
The matrix elements of B1, B2, B3, B4, RBCtransfuse and PLATELETtransfuse matrixes are as follows.
b14,1 = b17,1= b110,1= Sgdr/3, b15,2= b18,2= b111,2 = Ssdr/3, b16,3 = b19,3 = b112,3= Srdr/3, b213,4 = P1gdr, b214,5 = P1sdr, b215,6 = P1rdr, b226,23 = P1edr, b316,7= P2gdr, b317,8 = P2sdr, b318,9 = P2rdr, b327,24 = P2edr, b419,10 = P3gdr, b420,11 = P3sdr, b421,12 = P3rdr, b428,25 = P3edr, RBCtf13,1 =1, PLATELETtf19,1=1. The other matrix elements are equal to zero.
To implement the effect of pro-drug Gangcyclovir (img) the elements of matrix A are modified as follows:
a1,1(k) = (1+mSg-aSg-Sgdr-CSgs-CSgr -Cmyl×MylSg×dmyl(k) -Cg× imgsensSg×dimg(k)), a2,2(k)=(1+mSs-aSs-Ssdr-CSsg-CSsr-Cmyl ×MylSs×dmyl(k) -Cg×imgsensSs×dimg(k)), a3,3(k)=(1+mSr-aSr-Srdr- CSrs- CSrg -Cmyl×MylSr×dmyl(k) -Cg×imgsensSr×dimg(k)), a4,4(k) = (1+mP1g -aP1g -P1gdr -CP1gs -CP1gr -Cmyl×MylP1g×dmyl(k) -Cg ×imgsensP1g×dimg(k)), a5,5(k)=(1+mP1s-aP1s-P1sdr-CP1sg-CP1sr-Cmyl×MylP1s×dmyl(k) -Cg×imgsensP1s×dimg(k)), a6,6(k)=(1+mP1r-aP1r-P1rdr-CP1rg -CP1rs- Cmyl×MylP1r×dmyl(k) -Cg×imgsensP1r ×dimg(k)), a7,7(k)=(1+mP2g-aP2g-P2gdr- CP2gs-CP2gr-Cmyl×MylP2g ×dmyl(k)-Cg×imgsensP2g×dimg(k)+CCYT×dCYT(k) ×MCYT), a8,8(k)= (1+mP2s-aP2s -P2sdr-CP2sg-CP2sr -Cmyl×MylP2s×dmyl(k)- Cg× imgsensP2s×dimg(k)), a9,9(k)=(1+mP2r-aP2r-P2rdr-CP2rg-CP2rs- Cmyl×MylP2r×dmyl(k)-Cg×imgsensP2r×dimg(k)), a10,10(k)= (1+mP3g-aP3g- P3gdr-CP3gs-CP3gr-Cmyl×MylP3g×dmyl(k) -Cg×imgsensP3g ×dimg(k)), a11,11(k)=(1+mP3s-aP3s-P3sdr-CP3sg-CP3sr-Cmyl× MylP3s ×dmyl(k)- Cg×imgsensP3s×dimg(k)), a12,12(k)=(1 +mP3r-aP3r -P3rdr -CP3rg -CP3rs -Cmyl×MylP3r&tim es;dmyl(k) -Cg×imgsensP3r ×dimg(k)), a13,13(k)=(1-aB1g-CB1gs- CB1gr -Cmyl×MylB1g×dmyl(k) -Cg×imgsensB1g×dimg(k)), a14,14(k)=(1-aB1s- CB1sg-CB1sr-Cmyl× MylB1s×dmyl(k)-Cg×imgsensB1s×dimg(k)), a15,15(k) = (1-aB1r-CB1rg- CB1rs-Cmyl×MylB1r×dmyl(k)- Cg×imgsensB1r×dimg(k)), a16,16(k)= (1- aB2g -CB2gs -CB2gr -Cmyl×MylB2g×dmyl(k) -Cg×imgsensB2g× dimg(k)), a17,17(k)=(1-aB2s-CB2sg-CB2sr-Cmyl×MylB2s×dmyl(k) -Cg× imgsensB2s×dimg(k)), a18,18(k)=(1-aB2r-CB2rg -CB2rs -Cmyl×MylB2r ×dmyl(k)-Cg×imgsensB2r×dimg(k)), a19,19(k)=(1-aB3g-CB3gs-CB3gr- Cmyl×MylB3g×dmyl(k)-Cg×imgsensB3g×dimg(k)), a20,20(k)=(1-aB3s-CB3sg- CB3sr -Cmyl×MylB3s×dmyl(k) -Cg×imgsensB3s×dimg(k)), a21,21(k) = (1 -aB3r -CB3rg -CB3rs -Cmyl×MylB3r×dmyl(k) -Cg× imgsensB3r×dimg(k)), a22,22(k)=(1+mSe-aSe- Sedr -Cim×imsensSe× dim(k)), a23,23(k)=(1+mP1e- aP1e-P1edr -Cg×imgsensP1e×dimg(k)), a24,24(k) = (1+ mP2e -aP2e -P2edr -Cg×imgsensP2e×dimg(k) +CCYT× dCYT(k)×mDCYT), a25,25(k) = (1+mP3e-aP3e- P3edr -Cg×imgsensP3e× dimg(k)), a26,26(k)=(1-aB1e -Cg×imgsensB1e×dimg(k)), a27,27(k) = (1-aB2e -Cg×imgsensB2e×dimg(k)), a28,28(k) = (1-aB3e -Cg×imgsensB3e ×dimg(k)). The other elements of all the matrixes remain unchanged.
For initial parametric values we have followed . Free growth condition is depicted in Figures 1 and 2. Simulation observations suggest that with the increase in malignant cell population there is gradual decrease in the normal cell types in all the lineages. Conventionally in clinical practice HSCT is carried out after myeloablative chemotherapy and followed by different supportive therapy (i.e., transfusion of RBC and/or platelet). Here simulation is also carried out with this scheme. Simulation is carried out considering a worst leukemic condition. The differentiation rate of sensitive cell types is same as of resistive cell type in all the lineages and one-third cells become successful to be matured from progenitor cells.
Figure 1: Removal of leukemia after conventional myeloablative chemotherapy followed by stem cell transplantation with suicidal gene construct. Plots showing the behavior of stem cells (i), erythroblasts (ii), leukoblasts (iii), megakaryoblasts (iv) under freely growing condition (a) and with conventional chemotherapy followed by HSCT with tk-gene construct therapy (b). In all the plots, x-axis represents ‘Days’ and y-axis represents ‘Counts’. In all the plots, green, blue, black and magenta lines represent nonmalignant, sensitive type malignant, resistive type malignant cell populations and cells with donor origin respectively.
Figure 2: Removal of leukemia after conventional myeloablative chemotherapy followed by stem cell transplantation with suicidal gene construct. Plots showing the behavior of RBCs (i), WBCs (ii), platelets (iii) under freely growing condition (a) and with conventional chemotherapy followed by HSCT with tkgene construct therapy (b). In all the plots, x-axis represents ‘Days’ and y-axis represents ‘Counts’. In all the plots, green, blue, black and magenta lines represent nonmalignant, sensitive type malignant, resistive type malignant cell populations and cells with donor origin respectively.
Previously it has been shown that the day of transplantation is very crucial for successful outcome of stem cell transplantation. As per previous report  here, HSCT has been also carried out on day 180. Here simulation has been carried out with 100% HLA mismatch (high GVHD). Under this condition it is observed that mature lymphocytes from the transplanted stem cells becomes operative (i.e., lymphocytes of donor origin are developed and starts killing of leukemic cells) after day ~201 and leukemic stem cells are removed from system on day ~217.
The application of doses of pro-drug in two different schemes are considered as 0.01 and 0.02 and drug sensitivities of different cell types are imgsensSe = imgsensP1e = imgsensP2e = imgsensP3e = imgsensB1e = imgsensB2e = imgsensB3e = 0.92 both for the first and the second doses. While the sensitivities of other cell types are assumed to be zero.
The application of pro-drug Gangcyclovir is applied on day 210 with a low dose; however, application of high dose is started on day 240 along the continuation of low dose (Figure 3). Our simulation results suggest that if the low dose application is absent or only low dose is continued then there is chance of GVHD related mortality. Contrarily if high dose is applied early then there is chance of residual leukemia due to suppression of development of lymphocytes of donor origin. So, optimization in the application of drug schedule and drug dose is very important in the successful outcome of HSCT with tk-gene construct.
This therapeutic model is developed from the earlier three compartmental hematopoiesis model where methodology of incorporation of stochastic process at the parametric level has already been shown . To keep the simulation study simple we have avoided this in the present simulation runs. However one can incorporate the necessary stochastic components to test the effect of stochastic behaviors of individual system parameters towards the cancer dynamics as the developed model is flexible in nature.
Simulation runs showed that with allo-HSC transplantation with high HLA mismatch may be helpful in removal of malignancy but it requires a minimum kill factor (minimum killing efficiency of the developed lymphocytes from transplanted cells of donor origin). Though with the high HLA mis-match between donor and recipient there is possibility of high GVHD and GVHD related death of the recipient, but transplanted stem cells having the suicidal (tk) gene construct provides a control measure to GVHD mortality. Our study indicates that this scheme provides the opportunity of successful outcome of HSCT without any GVHD related detrimental effect on recipient. However, our model and simulation study also indicate that optimization of drug schedule (drug dose and time of drug application) is a crucial factor for the successful outcome of this therapy. Since earlier application of higher dose may stop the development of the donor lymphocyte thereby restrict the removal of malignant cell while after malignancy removal if lymphocytes donor origin are not removed from the host system then it may cause GVHD related mortality. Our model also indicates that the threshold of drug sensitivity of the lymphocytes is another important criteria and needed to check before going to such therapy and depending upon this factor, drug dose and drug application time are needed to be adjusted further.
Recent time for leukemia therapy, HSCT with suicidal tk gene construct is suggested by the different experts of the concerned scientific community [17,18]. However, due to unavailability of suitable analytical method, this procedure may be far from clinical practice. Here it is shown that the previously developed model  for leukemia therapy can be fitted for the assessment of such therapeutic procedure with a minor modification of variables and thus, provides the desired platform to bring such excellent therapeutic procedure in reality.
Presently, systems biologists are dependent on reverse engineering and from this stand-point, simulation runs with fitting of different parametric values of different systems component to an analytical model can help to identify the controlling variable(s) and their extent (boundary conditions) as well as its impact on the systems dynamical behavior. After applying any newer treatment procedure for a shorter duration of time, reassessment of the variables is needed and if there is any deviation between the simulation output and real life situation, then readjustment of initial parametric values of different variables are needed. Hence it will follow the predict-observe-correct cycle [16,20,21]. It is needless to point out here that as the developed model is based on difference equation and different feedbacks are incorporated in the form of inequation, so the model has the flexibility to change its parametric values of different variables and hence investigator can change the parametric values of any variables. Moreover, all the considered variables like cell counts, multiplication rate, differentiation time and leukemic cell killing capacity of immunocytes or drug can be captured by different hematological, immunological and cell culture based investigations.
Though HSCT with tk-gene construct has an immense potentiality for leukemia treatment; however, still is in laboratory phase. Presently cancers including leukemia treatment of individual patients are dependent on population based clinical trials. Recent time emphasis is imparted towards individualized treatment procedure. Though several large databases on leukemia are available; however, treatment dynamical data in public domain are unavailable. Moreover, clinically relevant suitable analytical platform is also unavailable; hence, priority in the area of Cancer Systems Biology is given towards the development of analytical model with a strong clinical rationality; so that, in future the efficacy and outcome of any newer cancer treatment procedure can be tested apriori by fitting the individual patients’ initial parametric values of different variables along with different therapeutic schedule of clinicians’ choice [16,20,22]. Towards this goal, this analytical model is developed.
Authors acknowledge the critical comments of the eminent founder members of the Society for Systems Biology & Translational Research.