// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;
unit item=scalar;
unit substance = 1E-6 mole;
unit time = 60 second;
unit volume = litre;
unit area = metre^2;
unit length = metre;

// 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 volume;
  private real EI.amt(time) substance;
  real EI(time) substance/volume;
  real EI.init substance/volume;
  private real PyrPI.amt(time) substance;
  real PyrPI(time) substance/volume;
  real PyrPI.init substance/volume;
  private real EIP.amt(time) substance;
  real EIP(time) substance/volume;
  real EIP.init substance/volume;
  private real HPr.amt(time) substance;
  real HPr(time) substance/volume;
  real HPr.init substance/volume;
  private real EIPHPr.amt(time) substance;
  real EIPHPr(time) substance/volume;
  real EIPHPr.init substance/volume;
  private real HPrP.amt(time) substance;
  real HPrP(time) substance/volume;
  real HPrP.init substance/volume;
  private real EIIA.amt(time) substance;
  real EIIA(time) substance/volume;
  real EIIA.init substance/volume;
  private real HPrPIIA.amt(time) substance;
  real HPrPIIA(time) substance/volume;
  real HPrPIIA.init substance/volume;
  private real EIIAP.amt(time) substance;
  real EIIAP(time) substance/volume;
  real EIIAP.init substance/volume;
  private real EIICB.amt(time) substance;
  real EIICB(time) substance/volume;
  real EIICB.init substance/volume;
  private real EIIAPIICB.amt(time) substance;
  real EIIAPIICB(time) substance/volume;
  real EIIAPIICB.init substance/volume;
  private real EIICBP.amt(time) substance;
  real EIICBP(time) substance/volume;
  real EIICBP.init substance/volume;
  private real EIICBPGlc.amt(time) substance;
  real EIICBPGlc(time) substance/volume;
  real EIICBPGlc.init substance/volume;
  private real PEP.amt(time) substance;
  real PEP substance/volume;
  real PEP.init substance/volume;
  private real Pyr.amt(time) substance;
  real Pyr substance/volume;
  real Pyr.init substance/volume;
  private real GlcP.amt(time) substance;
  real GlcP substance/volume;
  real GlcP.init substance/volume;
  private real Glc.amt(time) substance;
  real Glc substance/volume;
  real Glc.init substance/volume;
  real v1.rate(time) substance/time;
  real k1f = 1960;
  real k1r = 4.8E5;
  real v2.rate(time) substance/time;
  real k2f = 108000;
  real k2r = 294;
  real v3.rate(time) substance/time;
  real k3f = 14000;
  real k3r = 14000;
  real v4.rate(time) substance/time;
  real k4f = 84000;
  real k4r = 3360;
  real v5.rate(time) substance/time;
  real k5f = 21960;
  real k5r = 21960;
  real v6.rate(time) substance/time;
  real k6f = 4392;
  real k6r = 3384;
  real v7.rate(time) substance/time;
  real k7f = 880;
  real k7r = 880;
  real v8.rate(time) substance/time;
  real k8f = 2640;
  real k8r = 960;
  real v9.rate(time) substance/time;
  real k9f = 260;
  real k9r = 389;
  real v10.rate(time) substance/time;
  real k10f = 4800;
  real k10r = .0054;

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

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