// 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
// vHK: GLU MgATP => G6P MgADP 
// vPGI: G6P => F6P 
// vPFK: F6P MgATP => FBP MgADP 
// vALD: FBP => GAP DHAP 
// vTPI: DHAP => GAP 
// vGPDH: GAP Phi => PG13 NADH 
// vPGK: PG13 MgADP => PG3 MgATP 
// vPGM23: PG13 => PG23 
// vPGP23: PG23 => PG3 
// vPGLM: PG3 => PG2 
// vENO: PG2 => PEP 
// vPK: PEP MgADP => PYR MgATP 
// vLDH: PYR NADH => LAC NAD 
// vAK: MgATP AMP => MgADP ADP 
// vATPase: MgATP => MgADP 
// vMgATPdiss: MgATP => Mg ATP 
// vMgADPdiss: MgADP => Mg ADP 
// vMgAMPdiss: MgAMP => Mg AMP 
// vpentoseshunt: .21G6P => .071F6P .071GAP 

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

  // variable definitions
  real compartment = 1 L;
  real Vmax1HK = 15.8;
  real Vmax2HK = 33.2;
  real KHKMgATP = 1.44;
  real KHKMgATPMg = 1.14;
  real KeqHK = 10433;
  real KHKMg = 1.03;
  real KHKG6P = .0045;
  real KHKPG23 = 2.7;
  real KHKMgPG23 = 2.44;
  real VmaxPGI = 935;
  real KPGIG6P = .182;
  real KPGIF6P = .071;
  real KeqPGI = .444;
  real VmaxPFK = 239;
  real KPFKF6P = .1;
  real KPFKMgATP = .069;
  real KPFKMg = .44;
  real KPFKAMP = .033;
  real KPFKATP = .01;
  real L0PFK = .001072;
  real KeqPFK = 1075;
  real VmaxALD = 89.91;
  real KmALDDHAP = .0364;
  real KmALDGAP = .1906;
  real KmALDFBP = .0071;
  real KeqALD = .114;
  real KiALDGAP = .0572;
  real KiiALDGAP = .176;
  real VmaxTPI = 5456.6;
  real KTPIDHAP = .838;
  real KTPIGAP = .428;
  real KeqTPI = .051;
  real VmaxGPDH = 4300;
  real KeqGPDH = .00109;
  real KGPDHGAP = .005;
  real KGPDHPhi = 3.9;
  real KGPDHNAD = .05;
  real KGPDHNADH = .0083;
  real KGPDHPG13 = .0035;
  real VmaxPGK = 5E3;
  real KeqPGK = 2232;
  real KPGKMgADP = .35;
  real KPGKMgATP = .48;
  real KPGKPG13 = .002;
  real KPGKPG3 = 1.2;
  real kPGM23 = 1.1E5;
  real KPGMPG23 = .04;
  real KeqPGM23 = 1E5;
  real VmaxPGP23 = .52;
  real KeqPGP23 = 1E5;
  real KPGP23 = .2;
  real VmaxPGLM = 2E3;
  real KeqPGLM = .05;
  real KPGLMPG3 = 5;
  real KPGLMPG2 = 1;
  real VmaxENO = 1500;
  real KeqENO = 1.2;
  real KENOPEP = 1;
  real KENOPG2 = 1;
  real VmaxPK = 570;
  real KeqPK = 3349;
  real L0PK = 19;
  real KPKPEP = .225;
  real KPKMgADP = .474;
  real KPKFBP = .005;
  real KPKATP = 3.39;
  real kLDH = 1E6;
  real KeqLDH = 16979;
  real VmaxAK = 1380;
  real KeqAK = .4582;
  real KAKATP = .09;
  real KAKADP = .11;
  real KAKAMP = .08;
  real kATPase = 1.26;
  real kMgATPdiss = 1E6;
  real KeqMgATPdiss = .081;
  real kMgADPdiss = 1E6;
  real KeqMgADPdiss = .81;
  real kMgAMPdiss = 1E6;
  real KeqMgAMPdiss = 22.2;
  real kPS = 1;
  real GLU = 1 M;
  real MgATP(time) M;
  real G6P(time) M;
  real MgADP(time) M;
  real Mg(time) M;
  real PG23(time) M;
  real F6P(time) M;
  real FBP(time) M;
  real ATP(time) M;
  real AMP(time) M;
  real GAP(time) M;
  real DHAP(time) M;
  real Phi = .5 M;
  real PG13(time) M;
  real NADH = 4E-4 M;
  real PG3(time) M;
  real PG2(time) M;
  real PEP(time) M;
  real PYR(time) M;
  real LAC = 1.1 M;
  real NAD = .0886 M;
  real ADP(time) M;
  real MgAMP(time) M;
  real vHK(time) katal;
  real vPGI(time) katal;
  real vPFK(time) katal;
  real vALD(time) katal;
  real vTPI(time) katal;
  real vGPDH(time) katal;
  real vPGK(time) katal;
  real vPGM23(time) katal;
  real vPGP23(time) katal;
  real vPGLM(time) katal;
  real vENO(time) katal;
  real vPK(time) katal;
  real vLDH(time) katal;
  real vAK(time) katal;
  real vATPase(time) katal;
  real vMgATPdiss(time) katal;
  real vMgADPdiss(time) katal;
  real vMgAMPdiss(time) katal;
  real vpentoseshunt(time) katal;

  // equations
  when (time=time.min) MgATP = 1.4254;
  MgATP:time = (-1*vHK + -1*vPFK + vPGK + vPK + -1*vAK + -1*vATPase + -1*vMgATPdiss)/compartment;
  when (time=time.min) G6P = .05102;
  G6P:time = (vHK + -1*vPGI + -.21*vpentoseshunt)/compartment;
  when (time=time.min) MgADP = .18343;
  MgADP:time = (vHK + vPFK + -1*vPGK + -1*vPK + vAK + vATPase + -1*vMgADPdiss)/compartment;
  when (time=time.min) Mg = .754;
  Mg:time = (vMgATPdiss + vMgADPdiss + vMgAMPdiss)/compartment;
  when (time=time.min) PG23 = 4.676;
  PG23:time = (vPGM23 + -1*vPGP23)/compartment;
  when (time=time.min) F6P = .02251;
  F6P:time = (vPGI + -1*vPFK + .071*vpentoseshunt)/compartment;
  when (time=time.min) FBP = .00511;
  FBP:time = (vPFK + -1*vALD)/compartment;
  when (time=time.min) ATP = .15312;
  ATP:time = (vMgATPdiss)/compartment;
  when (time=time.min) AMP = .05534;
  AMP:time = (-1*vAK + vMgAMPdiss)/compartment;
  when (time=time.min) GAP = .00489;
  GAP:time = (vALD + vTPI + -1*vGPDH + .071*vpentoseshunt)/compartment;
  when (time=time.min) DHAP = .09614;
  DHAP:time = (vALD + -1*vTPI)/compartment;
  when (time=time.min) PG13 = 5.81E-4;
  PG13:time = (vGPDH + -1*vPGK + -1*vPGM23)/compartment;
  when (time=time.min) PG3 = .16473;
  PG3:time = (vPGK + vPGP23 + -1*vPGLM)/compartment;
  when (time=time.min) PG2 = .00793;
  PG2:time = (vPGLM + -1*vENO)/compartment;
  when (time=time.min) PEP = .0076;
  PEP:time = (vENO + -1*vPK)/compartment;
  when (time=time.min) PYR = .01435;
  PYR:time = (vPK + -1*vLDH)/compartment;
  when (time=time.min) ADP = .19704;
  ADP:time = (vAK + vMgADPdiss)/compartment;
  when (time=time.min) MgAMP = .00188;
  MgAMP:time = (-1*vMgAMPdiss)/compartment;
  vHK = Vmax1HK/KHKMgATP*(MgATP+Vmax2HK/Vmax1HK*(MgATP*Mg/KHKMgATPMg)-G6P*MgADP/KeqHK)/((1+MgATP/KHKMgATP)*(1+Mg/KHKMgATPMg)+Mg/KHKMg+(G6P/KHKG6P+1.55)*(1+Mg/KHKMg)+PG23/KHKPG23+Mg*PG23/(KHKMg*KHKMgPG23));
  vPGI = VmaxPGI*(G6P-F6P/KeqPGI)/(G6P+KPGIG6P*(1+F6P/KPGIF6P));
  vPFK = VmaxPFK*(F6P*MgATP-FBP*MgADP/KeqPFK)/((F6P+KPFKF6P)*(MgATP+KPFKMgATP)*(1+L0PFK*(1+ATP/KPFKATP)^4*(1+Mg/KPFKMg)^4/((1+AMP/KPFKAMP)^4*(1+F6P/KPFKF6P)^4)));
  vALD = VmaxALD/KmALDFBP*(FBP-GAP*DHAP/KeqALD)/(1+FBP/KmALDFBP+GAP/KiALDGAP+DHAP*(GAP+KmALDGAP)/(KmALDDHAP*KiALDGAP)+FBP*GAP/(KmALDFBP*KiiALDGAP));
  vTPI = VmaxTPI*(DHAP-GAP/KeqTPI)/(DHAP+KTPIDHAP*(1+GAP/KTPIGAP));
  vGPDH = VmaxGPDH/(KGPDHNAD*KGPDHGAP*KGPDHPhi)*(NAD*GAP*Phi-PG13*NADH/KeqGPDH)/((-1)+(1+NAD/KGPDHNAD)*(1+GAP/KGPDHGAP)*(1+Phi/KGPDHPhi)+(1+NADH/KGPDHNADH)*(1+PG13/KGPDHPG13));
  vPGK = VmaxPGK/(KPGKMgADP*KPGKPG13)*(MgADP*PG13-MgATP*PG3/KeqPGK)/((-1)+(1+MgADP/KPGKMgADP)*(1+PG13/KPGKPG13)+(1+MgATP/KPGKMgATP)*(1+PG3/KPGKPG3));
  vPGM23 = kPGM23*(PG13-PG23/KeqPGM23)/(1+PG23/KPGMPG23);
  vPGP23 = VmaxPGP23*(PG23-PG3/KeqPGP23)/(PG23+KPGP23);
  vPGLM = VmaxPGLM*(PG3-PG2/KeqPGLM)/(PG3+KPGLMPG3*(1+PG2/KPGLMPG2));
  vENO = VmaxENO*(PG2-PEP/KeqENO)/(PG2+KENOPG2*(1+PEP/KENOPEP));
  vPK = VmaxPK*(PEP*MgADP-PYR*MgATP/KeqPK)/((PEP+KPKPEP)*(MgADP+KPKMgADP)*(1+L0PK*((1+ATP/KPKATP)^4/((1+PEP/KPKPEP)^4*(1+FBP/KPKFBP)^4))));
  vLDH = kLDH*(KeqLDH*PYR*NADH-LAC*NAD);
  vAK = VmaxAK/(KAKATP*KAKAMP)*(MgATP*AMP-MgADP*ADP/KeqAK)/((1+MgATP/KAKATP)*(1+AMP/KAKAMP)+(MgADP+ADP)/KAKADP+MgADP*ADP/KAKADP^2);
  vATPase = kATPase*MgATP;
  vMgATPdiss = kMgATPdiss*(MgATP*KeqMgATPdiss-Mg*ATP);
  vMgADPdiss = kMgADPdiss*(MgADP*KeqMgADPdiss-Mg*ADP);
  vMgAMPdiss = kMgAMPdiss*(MgAMP*KeqMgAMPdiss-Mg*AMP);
  vpentoseshunt = 1*kPS;

  // variable properties
  compartment.sbmlRole="compartment";
  Vmax1HK.sbmlRole="parameter";
  Vmax2HK.sbmlRole="parameter";
  KHKMgATP.sbmlRole="parameter";
  KHKMgATPMg.sbmlRole="parameter";
  KeqHK.sbmlRole="parameter";
  KHKMg.sbmlRole="parameter";
  KHKG6P.sbmlRole="parameter";
  KHKPG23.sbmlRole="parameter";
  KHKMgPG23.sbmlRole="parameter";
  VmaxPGI.sbmlRole="parameter";
  KPGIG6P.sbmlRole="parameter";
  KPGIF6P.sbmlRole="parameter";
  KeqPGI.sbmlRole="parameter";
  VmaxPFK.sbmlRole="parameter";
  KPFKF6P.sbmlRole="parameter";
  KPFKMgATP.sbmlRole="parameter";
  KPFKMg.sbmlRole="parameter";
  KPFKAMP.sbmlRole="parameter";
  KPFKATP.sbmlRole="parameter";
  L0PFK.sbmlRole="parameter";
  KeqPFK.sbmlRole="parameter";
  VmaxALD.sbmlRole="parameter";
  KmALDDHAP.sbmlRole="parameter";
  KmALDGAP.sbmlRole="parameter";
  KmALDFBP.sbmlRole="parameter";
  KeqALD.sbmlRole="parameter";
  KiALDGAP.sbmlRole="parameter";
  KiiALDGAP.sbmlRole="parameter";
  VmaxTPI.sbmlRole="parameter";
  KTPIDHAP.sbmlRole="parameter";
  KTPIGAP.sbmlRole="parameter";
  KeqTPI.sbmlRole="parameter";
  VmaxGPDH.sbmlRole="parameter";
  KeqGPDH.sbmlRole="parameter";
  KGPDHGAP.sbmlRole="parameter";
  KGPDHPhi.sbmlRole="parameter";
  KGPDHNAD.sbmlRole="parameter";
  KGPDHNADH.sbmlRole="parameter";
  KGPDHPG13.sbmlRole="parameter";
  VmaxPGK.sbmlRole="parameter";
  KeqPGK.sbmlRole="parameter";
  KPGKMgADP.sbmlRole="parameter";
  KPGKMgATP.sbmlRole="parameter";
  KPGKPG13.sbmlRole="parameter";
  KPGKPG3.sbmlRole="parameter";
  kPGM23.sbmlRole="parameter";
  KPGMPG23.sbmlRole="parameter";
  KeqPGM23.sbmlRole="parameter";
  VmaxPGP23.sbmlRole="parameter";
  KeqPGP23.sbmlRole="parameter";
  KPGP23.sbmlRole="parameter";
  VmaxPGLM.sbmlRole="parameter";
  KeqPGLM.sbmlRole="parameter";
  KPGLMPG3.sbmlRole="parameter";
  KPGLMPG2.sbmlRole="parameter";
  VmaxENO.sbmlRole="parameter";
  KeqENO.sbmlRole="parameter";
  KENOPEP.sbmlRole="parameter";
  KENOPG2.sbmlRole="parameter";
  VmaxPK.sbmlRole="parameter";
  KeqPK.sbmlRole="parameter";
  L0PK.sbmlRole="parameter";
  KPKPEP.sbmlRole="parameter";
  KPKMgADP.sbmlRole="parameter";
  KPKFBP.sbmlRole="parameter";
  KPKATP.sbmlRole="parameter";
  kLDH.sbmlRole="parameter";
  KeqLDH.sbmlRole="parameter";
  VmaxAK.sbmlRole="parameter";
  KeqAK.sbmlRole="parameter";
  KAKATP.sbmlRole="parameter";
  KAKADP.sbmlRole="parameter";
  KAKAMP.sbmlRole="parameter";
  kATPase.sbmlRole="parameter";
  kMgATPdiss.sbmlRole="parameter";
  KeqMgATPdiss.sbmlRole="parameter";
  kMgADPdiss.sbmlRole="parameter";
  KeqMgADPdiss.sbmlRole="parameter";
  kMgAMPdiss.sbmlRole="parameter";
  KeqMgAMPdiss.sbmlRole="parameter";
  kPS.sbmlRole="parameter";
  GLU.sbmlRole="species";
  GLU.sbmlCompartment="compartment";
  MgATP.sbmlRole="species";
  MgATP.sbmlCompartment="compartment";
  G6P.sbmlRole="species";
  G6P.sbmlCompartment="compartment";
  MgADP.sbmlRole="species";
  MgADP.sbmlCompartment="compartment";
  Mg.sbmlRole="species";
  Mg.sbmlCompartment="compartment";
  PG23.sbmlRole="species";
  PG23.sbmlCompartment="compartment";
  F6P.sbmlRole="species";
  F6P.sbmlCompartment="compartment";
  FBP.sbmlRole="species";
  FBP.sbmlCompartment="compartment";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  AMP.sbmlRole="species";
  AMP.sbmlCompartment="compartment";
  GAP.sbmlRole="species";
  GAP.sbmlCompartment="compartment";
  DHAP.sbmlRole="species";
  DHAP.sbmlCompartment="compartment";
  Phi.sbmlRole="species";
  Phi.sbmlCompartment="compartment";
  PG13.sbmlRole="species";
  PG13.sbmlCompartment="compartment";
  NADH.sbmlRole="species";
  NADH.sbmlCompartment="compartment";
  PG3.sbmlRole="species";
  PG3.sbmlCompartment="compartment";
  PG2.sbmlRole="species";
  PG2.sbmlCompartment="compartment";
  PEP.sbmlRole="species";
  PEP.sbmlCompartment="compartment";
  PYR.sbmlRole="species";
  PYR.sbmlCompartment="compartment";
  LAC.sbmlRole="species";
  LAC.sbmlCompartment="compartment";
  NAD.sbmlRole="species";
  NAD.sbmlCompartment="compartment";
  ADP.sbmlRole="species";
  ADP.sbmlCompartment="compartment";
  MgAMP.sbmlRole="species";
  MgAMP.sbmlCompartment="compartment";
  vHK.sbmlRole="rate";
  vPGI.sbmlRole="rate";
  vPFK.sbmlRole="rate";
  vALD.sbmlRole="rate";
  vTPI.sbmlRole="rate";
  vGPDH.sbmlRole="rate";
  vPGK.sbmlRole="rate";
  vPGM23.sbmlRole="rate";
  vPGP23.sbmlRole="rate";
  vPGLM.sbmlRole="rate";
  vENO.sbmlRole="rate";
  vPK.sbmlRole="rate";
  vLDH.sbmlRole="rate";
  vAK.sbmlRole="rate";
  vATPase.sbmlRole="rate";
  vMgATPdiss.sbmlRole="rate";
  vMgADPdiss.sbmlRole="rate";
  vMgAMPdiss.sbmlRole="rate";
  vpentoseshunt.sbmlRole="rate";
}