He2017  A mathematical model of pancreatic cancer with two kinds of treatments
Model Identifier
BIOMD0000000811
Short description
This is a mathematical model of pancreatic cancer which includes descriptions of regulatory T cell activity and inhibition therapy. Descriptions of cytokine induced killer immunotherapy are also included.
Format
SBML
(L2V4)
Related Publication
 A mathematical model of pancreatic cancer with two kinds of treatments
 He, D.H., Xu, J.X.
 Journal of Biological Systems , 3/ 2017 , Volume 25 , Issue 1 , pages: 83104 , DOI: 10.1142/S021833901750005X
 Department of Mathematics, Zhejiang International Studies University, HangZhou, 310012, China
 In this paper, we investigate a mathematical model of pancreatic cancer, which extends the existing pancreatic cancer models with regulatory T cells (Tregs) and Treg inhibitory therapy. The model consists of tumorimmune interaction and immune suppression from Tregs. In the absence of treatments, we first characterize the system dynamics by locating equilibrium points and determining stability properties. Next, cytokine induced killer (CIK) immunotherapy is incorporated. Numerical simulations of prognostic results illustrate that the median overall survival associated with treatment can be prolonged approximately from 7 to 13 months, which is consistent with the clinical data. Furthermore, we consider cyclophosphamide (CTX) therapy as well as the combined therapy with CIK and CTX. Intensive simulation results suggest that both CTX therapy and the combined CIK/CTX therapy can reduce the number of Tregs and increase the overall survival (OS), but Tregs and tumor cells will gradually rise to equilibrium state as long as therapies are ceased.
Contributors
Submitter of the first revision: Johannes Meyer
Submitter of this revision: Rahuman Sheriff
Modellers: Rahuman Sheriff, Johannes Meyer
Submitter of this revision: Rahuman Sheriff
Modellers: Rahuman Sheriff, Johannes Meyer
Metadata information
is (2 statements)
isDerivedFrom (1 statement)
hasProperty (4 statements)
isDescribedBy (1 statement)
isDerivedFrom (1 statement)
hasProperty (4 statements)
Mathematical Modelling Ontology
Ordinary differential equation model
NCIt CytokineInduced Killer Cells
NCIt Neoplastic Growth or Spread
NCIt Pancreatic Neoplasm
NCIt CytokineInduced Killer Cells
NCIt Neoplastic Growth or Spread
NCIt Pancreatic Neoplasm
isDescribedBy (1 statement)
Curation status
Curated
Tags
Connected external resources
Name  Description  Size  Actions 

Model files 

He2017.xml  SBML L2V4 Representation of He2017  A mathematical model of pancreatic cancer with two kinds of treatments  113.04 KB  Preview  Download 
Additional files 

He2017.cps  COPASI file of He2017  A mathematical model of pancreatic cancer with two kinds of treatments  180.06 KB  Preview  Download 
He2017.sedml  SEDML file of He2017  A mathematical model of pancreatic cancer with two kinds of treatments  4.46 KB  Preview  Download 
 Model originally submitted by : Johannes Meyer
 Submitted: Sep 12, 2019 11:49:27 AM
 Last Modified: Oct 5, 2021 1:31:18 PM
Revisions

Version: 4
 Submitted on: Oct 5, 2021 1:31:18 PM
 Submitted by: Rahuman Sheriff
 With comment: Automatically added model identifier BIOMD0000000811

Version: 2
 Submitted on: Sep 12, 2019 11:49:27 AM
 Submitted by: Johannes Meyer
 With comment: Automatically added model identifier BIOMD0000000811
(*) You might be seeing discontinuous revisions as only public revisions are displayed here. Any private revisions of this model will only be shown to the submitter and their collaborators.
Legends
: Variable used inside SBML models
: Variable used inside SBML models
Species
Species  Initial Concentration/Amount 

N Killer natural killer cell 
4.816E8 item 
H T Helper helper Tlymphocyte 
2.1168E9 item 
R T Regulatory regulatory Tlymphocyte 
1.5876E8 item 
E CD8 CD8Positive TLymphocyte 
8.736E8 item 
P PSC pancreatic stellate cell 
5600000.0 item 
Reactions
Reactions  Rate  Parameters 

=> N_Killer; H_T_Helper  compartment*p_n*H_T_Helper*N_Killer/(g_n*tau_1_alpha_1+H_T_Helper)  g_n = 0.3; tau_1_alpha_1 = 2.2483E11; p_n = 0.125 
=> H_T_Helper  compartment*a_h  a_h = 360000.0 
R_T_Regulatory =>  compartment*delta_r*R_T_Regulatory  delta_r = 0.023 
=> E_CD8; N_Killer, C_PCC  compartment*r_e*N_Killer*C_PCC  r_e = 5.0E12 
=> E_CD8; H_T_Helper  compartment*p_e*H_T_Helper*E_CD8/(g_e*tau_1_alpha_1+H_T_Helper)  tau_1_alpha_1 = 2.2483E11; p_e = 0.125; g_e = 0.3 
E_CD8 =>  compartment*b_e*E_CD8  b_e = 0.02 
=> N_Killer  compartment*a_n  a_n = 130000.0 
=> R_T_Regulatory; H_T_Helper  compartment*b_r*H_T_Helper  b_r = 4.0E4 
P_PSC =>  compartment*lambda_p*P_PSC  lambda_p = 0.015 
H_T_Helper => ; R_T_Regulatory  compartment*delta_h*R_T_Regulatory*H_T_Helper  delta_h = 1.0E10 
Curator's comment:
(added: 12 Sep 2019, 11:49:17, updated: 12 Sep 2019, 11:49:17)
(added: 12 Sep 2019, 11:49:17, updated: 12 Sep 2019, 11:49:17)
Reproduced plot of Figure 1 in the original publication. Initial conditions are as indicated in the figure caption.
Model simulated and plot produced using COPASI 4.24 (Build 197).