Liu2017 - chemotherapy targeted model of tumor immune system

  public model
Model Identifier
Short description
Its a mathematical model reflecting chemotherapy response in tumor immune interaction system. Model encoded by Sarubini Kananathan and annotated+submitted to Biomodels by Krishna Tiwari
Related Publication
  • Dynamics of a tumor-immune model considering targeted chemotherapy
  • Peng Liu, Xijun Liu
  • Chaos, Solitons and Fractals , 1/ 2020 , pages: 7-13 , DOI: 10.1016/j.chaos.2017.03.002
  • Department of Mechanics, Tianjin University, Tianjin, 300354, PR China Tianjin Key Laboratory of Nonlinear Dynamics and Chaos Control, Tianjin, 300072, PR China
  • Considering the targeted chemotherapy, a mathematical model of tumor-immune system was constructed on the basis of de Pillis’s model. In this paper, we conducted qualitative analysis on the mathematical model, including the positivity and boundedness of solutions, local stability and global stability of equi- librium solutions. Some numerical simulations were given to illustrate the analytic results. Comparing the targeted chemotherapy model with regular chemotherapy model, we found that the targeted chemother- apy was benefit to kill tumor cells.
Submitter of the first revision: Krishna Kumar Tiwari
Submitter of this revision: Krishna Kumar Tiwari
Modellers: Krishna Kumar Tiwari

Metadata information

is (2 statements)
BioModels Database MODEL2004060002
BioModels Database BIOMD0000000930

hasTaxon (1 statement)
Taxonomy Homo sapiens

hasProperty (2 statements)
Gene Ontology immune response to tumor cell
Mathematical Modelling Ontology Ordinary differential equation model

isDescribedBy (1 statement)
occursIn (1 statement)

Curation status


Connected external resources

SBGN view in Newt Editor

Name Description Size Actions

Model files

Liu2017.xml SBML V2.4 file for the model 54.01 KB Preview | Download

Additional files

Liu2017.cps COPASI 4.27(built217) file 92.36 KB Preview | Download
Liu2017.sedml SEDML file for the model 4.90 KB Preview | Download

  • Model originally submitted by : Krishna Kumar Tiwari
  • Submitted: Apr 6, 2020 7:10:47 PM
  • Last Modified: Apr 6, 2020 7:10:47 PM
  • Version: 2 public model Download this version
    • Submitted on: Apr 6, 2020 7:10:47 PM
    • Submitted by: Krishna Kumar Tiwari
    • With comment: Automatically added model identifier BIOMD0000000930
: Variable used inside SBML models

Species Initial Concentration/Amount
Chemotherapeutic drug concentration M

0.45 mmol
Tumor Cell Population T

Neoplastic Cell
1.0E7 mmol
Circulating lymphocyte population C

BTO:0000775 ; OBI 0000181
6.25E10 mmol
Effector immune cell population N

OBI 1110016
300000.0 mmol
Reactions Rate Parameters
Chemotherapeutic_drug_concentration__M => Tumor_Microenvironment*gamma*Chemotherapeutic_drug_concentration__M gamma = 0.9
=> Tumor_Cell_Population__T Tumor_Microenvironment*a*Tumor_Cell_Population__T*(1-b*Tumor_Cell_Population__T) a = 0.431; b = 1.02E-14
Circulating_lymphocyte_population__C => ; Chemotherapeutic_drug_concentration__M Tumor_Microenvironment*Kc*(1-eta)*Chemotherapeutic_drug_concentration__M*Circulating_lymphocyte_population__C Kc = 0.6; eta = 0.0
=> Effector_immune_cell_population__N; Tumor_Cell_Population__T Tumor_Microenvironment*g*Tumor_Cell_Population__T/(s+Tumor_Cell_Population__T)*Effector_immune_cell_population__N s = 20.2; g = 0.015
Effector_immune_cell_population__N => ; Chemotherapeutic_drug_concentration__M Tumor_Microenvironment*Kn*(1-eta)*Chemotherapeutic_drug_concentration__M*Effector_immune_cell_population__N Kn = 0.6; eta = 0.0
=> Effector_immune_cell_population__N Tumor_Microenvironment*alpha1 alpha1 = 12000.0
Effector_immune_cell_population__N => Tumor_Microenvironment*mu*Effector_immune_cell_population__N mu = 0.0412
=> Chemotherapeutic_drug_concentration__M Tumor_Microenvironment*Vm Vm = 0.45
Effector_immune_cell_population__N => ; Tumor_Cell_Population__T Tumor_Microenvironment*p*Tumor_Cell_Population__T*Effector_immune_cell_population__N p = 2.0E-11
=> Circulating_lymphocyte_population__C Tumor_Microenvironment*alpha2 alpha2 = 7.5E8
Curator's comment:
(added: 06 Apr 2020, 19:10:23, updated: 06 Apr 2020, 19:10:23)
Model encoded and plot generated using COPASI 4.27 (Built417). Figure 2A of the literature is reproduced. Run time is 50 days with interval 100.