// Created by libAntimony v2.8.0
model *Westermark2003_Pancreatic_GlycOsc_extended()
// Compartments and Species:
compartment cell;
species $GLC in cell, G6P_F6P in cell, $F6P in cell, FBP in cell, $G3P in cell;
species DHAP_G3P in cell, $DHAP in cell;
// Assignment Rules:
F6P := (G6P_F6P*KeqGPI)/(1 + KeqGPI);
G3P := (DHAP_G3P*KeqTPI)/(1 + KeqTPI);
DHAP := DHAP_G3P - G3P;
Vgk := (Vgk_min*dw_per_ml)/min_to_sec;
Vpfk := (Vpfk_min*dw_per_ml)/min_to_sec;
Vfba := (Vfba_min*dw_per_ml)/min_to_sec;
Vgapdh := (Vgapdh_min*dw_per_ml)/min_to_sec;
sigma := F6P/Spfk;
// Reactions:
vgk: $GLC -> G6P_F6P; (cell*Vgk*(GLC/Sgk)^hGK)/(1 + (GLC/Sgk)^hGK);
vpfk: G6P_F6P -> FBP; (cell*Vpfk*(F6P/Spfk)^(hpfk - (hpfk - hact)*(FBP/Sfba/(1 + FBP/Sfba))))/((F6P/Spfk)^(hpfk - (hpfk - hact)*(FBP/Sfba/(1 + FBP/Sfba))) + (1 + (FBP/Xpfk)^hx)/(1 + alpha^(hpfk - (hpfk - hact)*(FBP/Sfba/(1 + FBP/Sfba)))*(FBP/Xpfk)^hx));
vfba: FBP -> 2DHAP_G3P; (cell*Vfba*(FBP/Sfba - (G3P*DHAP)/(Pfba*Qfba*KeqFBA)))/(1 + FBP/Sfba + DHAP/Qfba + (G3P*DHAP)/(Pfba*Qfba));
vgapdh: DHAP_G3P => ; (cell*Vgapdh*G3P)/(Sgapdh + G3P);
// Species initializations:
GLC = 10;
G6P_F6P = 3.71728;
FBP = 0.00063612;
DHAP_G3P = 0.00262966;
// Compartment initializations:
cell = 1;
// Variable initializations:
KeqGPI = 0.3;
KeqGPI has dimensionless;
KeqTPI = 0.045455;
KeqTPI has dimensionless;
Vgk has mM_per_s;
Vgk_min = 10;
Vgk_min has mmole_per_min_g;
dw_per_ml = 0.3333;
dw_per_ml has g_per_ml;
min_to_sec = 60;
min_to_sec has dimensionless;
hGK = 1.7;
hGK has dimensionless;
Vpfk has mM_per_s;
Vpfk_min = 100;
Vpfk_min has mmole_per_min_g;
Vfba has mM_per_s;
Vfba_min = 25;
Vfba_min has mmole_per_min_g;
Vgapdh has mM_per_s;
Vgapdh_min = 250;
Vgapdh_min has mmole_per_min_g;
Sgk = 8;
Sgk has mM;
Spfk = 4;
Spfk has mM;
Sfba = 0.005;
Sfba has mM;
Sgapdh = 0.005;
Sgapdh has mM;
Xpfk = 0.01;
Xpfk has mM;
alpha = 5;
alpha has dimensionless;
hx = 2.5;
hx has dimensionless;
hpfk = 2.5;
hpfk has dimensionless;
hact = 1;
hact has dimensionless;
Pfba = 0.5;
Pfba has mM;
Qfba = 0.275;
Qfba has mM;
KeqFBA = 0.1;
KeqFBA has mM;
// Other declarations:
var Vgk, Vpfk, Vfba, Vgapdh, sigma;
const cell, KeqGPI, KeqTPI, Vgk_min, dw_per_ml, min_to_sec, hGK, Vpfk_min;
const Vfba_min, Vgapdh_min, Sgk, Spfk, Sfba, Sgapdh, Xpfk, alpha, hx, hpfk;
const hact, Pfba, Qfba, KeqFBA;
// Unit definitions:
unit substance = 1e-3 mole;
unit time_unit = second;
unit mM = 1e-3 mole / litre;
unit g_per_ml = 1e-3 kilogram / 1e-3 litre;
unit mM_per_s = 1e-3 mole / (second * litre);
unit mmole_per_min_g = 1e-3 mole / (60 second * kilogram);
// Display Names:
substance is "mmole";
time_unit is "seconds";
g_per_ml is "gramm per ml";
mM_per_s is "mM per sec";
mmole_per_min_g is "mmole per (min kg)";
GLC is "intracellular glucose";
F6P is "fructose-6-phosphate";
FBP is "fructose-1,6-biphosphate";
G3P is "glyceraldehyde--phosphate";
DHAP_G3P is "DHAP-G3P pool";
DHAP is "dihydroxyacetone-phosphate";
end