// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// v1: PEP EI => PyrPI 
// v2: PyrPI => EIP Pyr 
// v3: HPr EIP => EIPHPr 
// v4: EIPHPr => HPrP EI 
// v5: HPrP EIIA => HPrPIIA 
// v6: HPrPIIA => EIIAP HPr 
// v7: EIICB EIIAP => EIIAPIICB 
// v8: EIIAPIICB => EIICBP EIIA 
// v9: EIICBP Glc => EIICBPGlc 
// v10: EIICBPGlc => EIICB GlcP 

math main {
  realDomain time second;
  time.min=0;
  extern time.max;
  extern time.delta;

  // variable definitions
  real compartment = 1 L;
  real EI(time) uM;
  real PyrPI(time) uM;
  real EIP(time) uM;
  real HPr(time) uM;
  real EIPHPr(time) uM;
  real HPrP(time) uM;
  real EIIA(time) uM;
  real HPrPIIA(time) uM;
  real EIIAP(time) uM;
  real EIICB(time) uM;
  real EIIAPIICB(time) uM;
  real EIICBP(time) uM;
  real EIICBPGlc(time) uM;
  real PEP = 2800 uM;
  real Pyr = 900 uM;
  real GlcP = 50 uM;
  real Glc = 500 uM;
  real v1(time) umol/min;
  real k1f = 1960;
  real k1r = 4.8E5;
  real v2(time) umol/min;
  real k2f = 108000;
  real k2r = 294;
  real v3(time) umol/min;
  real k3f = 14000;
  real k3r = 14000;
  real v4(time) umol/min;
  real k4f = 84000;
  real k4r = 3360;
  real v5(time) umol/min;
  real k5f = 21960;
  real k5r = 21960;
  real v6(time) umol/min;
  real k6f = 4392;
  real k6r = 3384;
  real v7(time) umol/min;
  real k7f = 880;
  real k7r = 880;
  real v8(time) umol/min;
  real k8f = 2640;
  real k8r = 960;
  real v9(time) umol/min;
  real k9f = 260;
  real k9r = 389;
  real v10(time) umol/min;
  real k10f = 4800;
  real k10r = .0054;

  // equations
  when (time=time.min) EI = 3;
  EI:time = (-1*v1 + v4)/compartment;
  when (time=time.min) PyrPI = 0;
  PyrPI:time = (v1 + -1*v2)/compartment;
  when (time=time.min) EIP = 2;
  EIP:time = (v2 + -1*v3)/compartment;
  when (time=time.min) HPr = 25;
  HPr:time = (-1*v3 + v6)/compartment;
  when (time=time.min) EIPHPr = 0;
  EIPHPr:time = (v3 + -1*v4)/compartment;
  when (time=time.min) HPrP = 25;
  HPrP:time = (v4 + -1*v5)/compartment;
  when (time=time.min) EIIA = 20;
  EIIA:time = (-1*v5 + v8)/compartment;
  when (time=time.min) HPrPIIA = 0;
  HPrPIIA:time = (v5 + -1*v6)/compartment;
  when (time=time.min) EIIAP = 20;
  EIIAP:time = (v6 + -1*v7)/compartment;
  when (time=time.min) EIICB = 5;
  EIICB:time = (-1*v7 + v10)/compartment;
  when (time=time.min) EIIAPIICB = 0;
  EIIAPIICB:time = (v7 + -1*v8)/compartment;
  when (time=time.min) EIICBP = 5;
  EIICBP:time = (v8 + -1*v9)/compartment;
  when (time=time.min) EIICBPGlc = 0;
  EIICBPGlc:time = (v9 + -1*v10)/compartment;
  v1 = compartment*(k1f*PEP*EI-k1r*PyrPI);
  v2 = compartment*(k2f*PyrPI-k2r*Pyr*EIP);
  v3 = compartment*(k3f*EIP*HPr-k3r*EIPHPr);
  v4 = compartment*(k4f*EIPHPr-k4r*EI*HPrP);
  v5 = compartment*(k5f*HPrP*EIIA-k5r*HPrPIIA);
  v6 = compartment*(k6f*HPrPIIA-k6r*HPr*EIIAP);
  v7 = compartment*(k7f*EIIAP*EIICB-k7r*EIIAPIICB);
  v8 = compartment*(k8f*EIIAPIICB-k8r*EIIA*EIICBP);
  v9 = compartment*(k9f*EIICBP*Glc-k9r*EIICBPGlc);
  v10 = compartment*(k10f*EIICBPGlc-k10r*EIICB*GlcP);

  // variable properties
  compartment.sbmlRole="compartment";
  compartment.sbmlName="cell";
  EI.sbmlRole="species";
  EI.sbmlCompartment="compartment";
  PyrPI.sbmlRole="species";
  PyrPI.sbmlCompartment="compartment";
  EIP.sbmlRole="species";
  EIP.sbmlCompartment="compartment";
  HPr.sbmlRole="species";
  HPr.sbmlCompartment="compartment";
  EIPHPr.sbmlRole="species";
  EIPHPr.sbmlCompartment="compartment";
  HPrP.sbmlRole="species";
  HPrP.sbmlCompartment="compartment";
  EIIA.sbmlRole="species";
  EIIA.sbmlCompartment="compartment";
  HPrPIIA.sbmlRole="species";
  HPrPIIA.sbmlCompartment="compartment";
  EIIAP.sbmlRole="species";
  EIIAP.sbmlCompartment="compartment";
  EIICB.sbmlRole="species";
  EIICB.sbmlCompartment="compartment";
  EIIAPIICB.sbmlRole="species";
  EIIAPIICB.sbmlCompartment="compartment";
  EIICBP.sbmlRole="species";
  EIICBP.sbmlCompartment="compartment";
  EIICBPGlc.sbmlRole="species";
  EIICBPGlc.sbmlCompartment="compartment";
  PEP.sbmlRole="species";
  PEP.sbmlCompartment="compartment";
  Pyr.sbmlRole="species";
  Pyr.sbmlCompartment="compartment";
  GlcP.sbmlRole="species";
  GlcP.sbmlCompartment="compartment";
  Glc.sbmlRole="species";
  Glc.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  k1f.sbmlRole="parameter";
  k1r.sbmlRole="parameter";
  v2.sbmlRole="rate";
  k2f.sbmlRole="parameter";
  k2r.sbmlRole="parameter";
  v3.sbmlRole="rate";
  k3f.sbmlRole="parameter";
  k3r.sbmlRole="parameter";
  v4.sbmlRole="rate";
  k4f.sbmlRole="parameter";
  k4r.sbmlRole="parameter";
  v5.sbmlRole="rate";
  k5f.sbmlRole="parameter";
  k5r.sbmlRole="parameter";
  v6.sbmlRole="rate";
  k6f.sbmlRole="parameter";
  k6r.sbmlRole="parameter";
  v7.sbmlRole="rate";
  k7f.sbmlRole="parameter";
  k7r.sbmlRole="parameter";
  v8.sbmlRole="rate";
  k8f.sbmlRole="parameter";
  k8r.sbmlRole="parameter";
  v9.sbmlRole="rate";
  k9f.sbmlRole="parameter";
  k9r.sbmlRole="parameter";
  v10.sbmlRole="rate";
  k10f.sbmlRole="parameter";
  k10r.sbmlRole="parameter";
}