Zhu2015 - Combined gemcitabine and birinapant in pancreatic cancer cells - basic PD model

  public model
Model Identifier
BIOMD0000000668
Short description
Zhu2015 - Combined gemcitabine and birinapant in pancreatic cancer cells - basic PD model
Mathematical model to illustrate the effectiveness of combination chemotherapy involving gemcitabine and birinapant against pancreatic cancer.

This model is described in the article:

Zhu X, Straubinger RM, Jusko WJ.
J Pharmacokinet Pharmacodyn 2015 Oct; 42(5): 477-496

Abstract:

Combination chemotherapy is standard treatment for pancreatic cancer. However, current drugs lack efficacy for most patients, and selection and evaluation of new combination regimens is empirical and time-consuming. The efficacy of gemcitabine, a standard-of-care agent, combined with birinapant, a pro-apoptotic antagonist of Inhibitor of Apoptosis Proteins (IAPs), was investigated in pancreatic cancer cells. PANC-1 cells were treated with vehicle, gemcitabine (6, 10, 20 nM), birinapant (50, 200, 500 nM), and combinations of the two drugs. Temporal changes in cell numbers, cell cycle distribution, and apoptosis were measured. A basic pharmacodynamic (PD) model based on cell numbers, and a mechanism-based PD model integrating all measurements, were developed. The basic PD model indicated that synergistic effects occurred in both cell proliferation and death processes. The mechanism-based model captured key features of drug action: temporary cell cycle arrest in S phase induced by gemcitabine alone, apoptosis induced by birinapant alone, and prolonged cell cycle arrest and enhanced apoptosis induced by the combination. A drug interaction term Ψ was employed in the models to signify interactions of the combination when data were limited. When more experimental information was utilized, Ψ values approaching 1 indicated that specific mechanisms of interactions were captured better. PD modeling identified the potential benefit of combining gemcitabine and birinapant, and characterized the key interaction pathways. An optimal treatment schedule of pretreatment with gemcitabine for 24-48 h was suggested based on model predictions and was verified experimentally. This approach provides a generalizable modeling platform for exploring combinations of cytostatic and cytotoxic agents in cancer cell culture studies.

This model is hosted on BioModels Database and identified by: BIOMD0000000668.

To cite BioModels Database, please use: Chelliah V et al. BioModels: ten-year anniversary. Nucl. Acids Res. 2015, 43(Database issue):D542-8.

To the extent possible under law, all copyright and related or neighbouring rights to this encoded model have been dedicated to the public domain worldwide. Please refer to CC0 Public Domain Dedication for more information.

Format
SBML (L2V4)
Related Publication
  • Mechanism-based mathematical modeling of combined gemcitabine and birinapant in pancreatic cancer cells.
  • Zhu X, Straubinger RM, Jusko WJ
  • Journal of pharmacokinetics and pharmacodynamics , 10/ 2015 , Volume 42 , Issue 5 , pages: 477-496 , PubMed ID: 26252969
  • Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, State University of New York at Buffalo, Buffalo, NY, 14214, USA.
  • Combination chemotherapy is standard treatment for pancreatic cancer. However, current drugs lack efficacy for most patients, and selection and evaluation of new combination regimens is empirical and time-consuming. The efficacy of gemcitabine, a standard-of-care agent, combined with birinapant, a pro-apoptotic antagonist of Inhibitor of Apoptosis Proteins (IAPs), was investigated in pancreatic cancer cells. PANC-1 cells were treated with vehicle, gemcitabine (6, 10, 20 nM), birinapant (50, 200, 500 nM), and combinations of the two drugs. Temporal changes in cell numbers, cell cycle distribution, and apoptosis were measured. A basic pharmacodynamic (PD) model based on cell numbers, and a mechanism-based PD model integrating all measurements, were developed. The basic PD model indicated that synergistic effects occurred in both cell proliferation and death processes. The mechanism-based model captured key features of drug action: temporary cell cycle arrest in S phase induced by gemcitabine alone, apoptosis induced by birinapant alone, and prolonged cell cycle arrest and enhanced apoptosis induced by the combination. A drug interaction term Ψ was employed in the models to signify interactions of the combination when data were limited. When more experimental information was utilized, Ψ values approaching 1 indicated that specific mechanisms of interactions were captured better. PD modeling identified the potential benefit of combining gemcitabine and birinapant, and characterized the key interaction pathways. An optimal treatment schedule of pretreatment with gemcitabine for 24-48 h was suggested based on model predictions and was verified experimentally. This approach provides a generalizable modeling platform for exploring combinations of cytostatic and cytotoxic agents in cancer cell culture studies.
Contributors
Submitter of the first revision: Vijayalakshmi Chelliah
Submitter of this revision: administrator
Modellers: administrator, Vijayalakshmi Chelliah

Metadata information

is (2 statements)
BioModels Database MODEL1604270000
BioModels Database BIOMD0000000668

isDescribedBy (3 statements)
isInstanceOf (5 statements)
KEGG Drug Gemcitabine (USAN/INN)
NCIt Combination Chemotherapy
Brenda Tissue Ontology pancreas
BioModels Database MODEL1604270000
KEGG Drug Birinapant (USAN/INN)

hasPart (2 statements)
occursIn (1 statement)
Brenda Tissue Ontology pancreas


Curation status
Curated

Tags

Connected external resources

SBGN view in Newt Editor

Name Description Size Actions

Model files

BIOMD0000000668_url.xml SBML L2V4 representation of Zhu2015 - Combined gemcitabine and birinapant in pancreatic cancer cells - basic PD model 86.00 KB Preview | Download

Additional files

BIOMD0000000668-biopax2.owl Auto-generated BioPAX (Level 2) 16.29 KB Preview | Download
BIOMD0000000668-biopax3.owl Auto-generated BioPAX (Level 3) 25.78 KB Preview | Download
BIOMD0000000668.m Auto-generated Octave file 9.02 KB Preview | Download
BIOMD0000000668.pdf Auto-generated PDF file 168.77 KB Preview | Download
BIOMD0000000668.png Auto-generated Reaction graph (PNG) 31.03 KB Preview | Download
BIOMD0000000668.sci Auto-generated Scilab file 154.00 Bytes Preview | Download
BIOMD0000000668.svg Auto-generated Reaction graph (SVG) 21.75 KB Preview | Download
BIOMD0000000668.vcml Auto-generated VCML file 900.00 Bytes Preview | Download
BIOMD0000000668.xpp Auto-generated XPP file 5.97 KB Preview | Download
BIOMD0000000668_urn.xml Auto-generated SBML file with URNs 85.93 KB Preview | Download
MODEL1604270000_basic.cps Curated and annotated COPASI file. 95.72 KB Preview | Download
MODEL1604270000_basic.sedml SED-ML file to produce a similar figure to figure 4a of the reference publication. Concentration of birinapant is set to zero while a parameter scan varies the concentration of gemcitabine from 0 to 20 in increments of 5 nM. 4.13 KB Preview | Download

  • Model originally submitted by : Vijayalakshmi Chelliah
  • Submitted: Apr 27, 2016 3:53:09 PM
  • Last Modified: Feb 7, 2018 1:16:38 PM
Revisions
  • Version: 2 public model Download this version
    • Submitted on: Feb 7, 2018 1:16:38 PM
    • Submitted by: administrator
    • With comment: Current curated version of Zhu2015_basic_ADAPT5
  • Version: 1 public model Download this version
    • Submitted on: Apr 27, 2016 3:53:09 PM
    • Submitted by: Vijayalakshmi Chelliah
    • With comment: Original import of Zhu2015_basic_ADAPT5.txt

(*) You might be seeing discontinuous revisions as only public revisions are displayed here. Any private revisions unpublished model revision of this model will only be shown to the submitter and their collaborators.

Legends
: Variable used inside SBML models


Species
Species Initial Concentration/Amount
Sti b1

delay
0.0 mmol
Inh b

Chemotherapy ; Inhibition of Cancer Cell Growth ; Chemotherapy
0.0 mmol
Sti g2

delay
0.0 mmol
Sti g3

delay
0.0 mmol
Ra

PANC-1 cell ; Adhesion
307000.0 mmol
Sti b4

delay
0.0 mmol
Rd

PANC-1 cell ; Detached ; cell death
1940.0 mmol
Sti b

Chemotherapy ; Positive Regulation of Cell Death ; Chemotherapy
0.0 mmol
Sti g1

delay
0.0 mmol
Reactions
Reactions Rate Parameters
Sti_b1 => Sti_b2 Pancreas*ktau_b*Sti_b1 ktau_b = 0.611
Inh_b = Imax_b*C_b^Hi_b/((Psi_i*IC50_b)^Hi_b+C_b^Hi_b) [] Hi_b = 1.06; C_b = 0.0; IC50_b = 145.0; Imax_b = 0.375; Psi_i = 1.0
Sti_b => Sti_b1 Pancreas*ktau_b*Sti_b ktau_b = 0.611
Sti_g2 => Sti_g3 Pancreas*ktau_g*Sti_g2 ktau_g = 0.086
Ra = (1-Inh_g)*(1-Inh_b)*ModelValue_3*Ra*(1-Ra/ModelValue_2)-(1+Sti_g4)*(1+Sti_b4)*ModelValue_4*Ra (1-Inh_g)*(1-Inh_b)*ModelValue_3*Ra*(1-Ra/ModelValue_2)-(1+Sti_g4)*(1+Sti_b4)*ModelValue_4*Ra ModelValue_4 = 3.85E-4; ModelValue_3 = 0.0209; ModelValue_2 = 5490000.0
Sti_b3 => Sti_b4 Pancreas*ktau_b*Sti_b3 ktau_b = 0.611
Rd = (1+Sti_g4)*(1+Sti_b4)*ModelValue_4*Ra-ModelValue_4*Rd (1+Sti_g4)*(1+Sti_b4)*ModelValue_4*Ra-ModelValue_4*Rd ModelValue_4 = 3.85E-4
Sti_b = Smax_b*C_b^Hs_b/((Psi_s*SC50_b)^Hs_b+C_b^Hs_b) [] Psi_s = 1.0; C_b = 0.0; SC50_b = 168.0; Smax_b = 17.5; Hs_b = 0.984
Sti_g => Sti_g1 Pancreas*ktau_g*Sti_g ktau_g = 0.086
Curator's comment:
(added: 07 Feb 2018, 13:10:49, updated: 07 Feb 2018, 13:10:49)
Similar figures of figure 4 of the reference publication have been produced with Psi_i=1 and Psi_s=1 as opposed to Psi_i=0.582 and Psi_s=0.583 for figures 4E and 4F as listed in table 1 of the reference publication. The figures illustrate the number of proliferating (left) and dead (right) cells after 100 hours of gemcitabine (top), birinapant (middle) and combination (bottom) treatment. Black, blue, green and red curves correspond to increasing drug concentration, as used in the reference publication. The number of proliferating cells decreased after a combination of 20nM gemcitabine + 500 nM birinapant (figure 4E, bottom left, red curve). The simulations were performed in COPASI V4.22 (Build 170) and figures were generated in MATLAB R2014.