// 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
// vgp: Ph GLY => G1P 
// vpglm: G1P => G6P 
// vpgi: G6P => F6P 
// vpfk: F6P atp => FDP adp 
// vald: FDP => DHAP GAP 
// vtpi: GAP => DHAP 
// vgapdh: GAP NAD Ph => DPG NADH 
// vpgk: DPG adp => P3G atp 
// vpgm: P3G => P2G 
// ven: P2G => PEP 
// vpk: PEP adp => PYR atp 
// vldh: PYR NADH => LAC NAD 

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

  // variable definitions
  real compartment = 1 L;
  real Vfgly = .15;
  real KgpAigly = .002;
  real KgpAipi = .0047;
  real KgpApi = .004;
  real KgpAglyf = .0017;
  real KgpAglyb = 1.5E-4;
  real KgpAg1p = .0027;
  real KgpAig1p = .0101;
  real KgpBipi = .0046;
  real KgpBiglyf = .015;
  real KgpBpi = 2E-4;
  real KgpBg1p = .0015;
  real KgpBiglyb = .0044;
  real KgpBig1p = .0074;
  real Kgpamp = 9.7E-5;
  real fracA = .4;
  real fracB = .6;
  real nH = 1.75;
  real Vfpglm = .48;
  real Kpglmg1p = 6.3E-5;
  real Kpglmg6p = 3E-5;
  real Vbpgi = .88;
  real Kpgig6p = 4.8E-4;
  real Kpgif6p = 1.19E-4;
  real Vfpfk = .056;
  real Kpfkfdp = .0042;
  real KpfkfdpT = .0042;
  real Kpfkadp = .00271;
  real KpfkadpT = .00271;
  real Kpfkf6p = 1.8E-4;
  real Kpfkf6pT = .02;
  real Kpfkatp = 8E-5;
  real KpfkatpT = 2.5E-4;
  real Kpfkamp = 6E-5;
  real Kpfkiatp = 8.7E-4;
  real Vfald = .104;
  real Kaldfdp = 5E-4;
  real Kaldgap = .001;
  real Kalddhap = .002;
  real Vftpi = 12;
  real Ktpigap = 3.2E-4;
  real Ktpidhap = 6.1E-4;
  real Vfgad = 1.65;
  real Kgapdhgap = 2.5E-6;
  real Kgapdhnad = 9E-5;
  real Kgapdhpi = 2.9E-4;
  real Kgapdhnadh = 3.3E-6;
  real Kgapdh13dpg = 8E-7;
  real Vbpgk = 1.12;
  real Kpgk13dpg = 1.9E-6;
  real Kpgkadp = 8.3E-5;
  real Kpgk3pg = .0012;
  real Kpgkatp = 3.5E-4;
  real Vfpgm = 1.12;
  real Kpgm3pg = 2E-4;
  real Kpgm2pg = 1.4E-5;
  real Vfen = .5;
  real Ken2pg = 1E-4;
  real Kenpep = 3.7E-4;
  real Vfpk = 3;
  real Kpkpep = 8E-5;
  real Kpkpyr = .00705;
  real Kpkadp = 3E-4;
  real Kpkatp = .00113;
  real Vfldh = 1.92;
  real Kldhpyr = 3.35E-4;
  real Kldhnadh = 2.167E-6;
  real Kldhlac = .01717;
  real Kldhnad = 8.49E-4;
  real VrevCK = .5;
  real KckPCr = .00111;
  real Kckiatp = .0035;
  real Kckcr = .0038;
  real Kckiadp = 1.35E-4;
  real KckiPCr = .0039;
  real KeqCK = 233;
  real Vfadk = .88;
  real Kadkadp = 3.5E-4;
  real Kadkamp = 3.2E-4;
  real Kadkatp = 2.7E-4;
  real k = .075;
  real kout = .2;
  real dn = .01;
  real en = .01;
  real Lo = 13;
  real amp = 9.5209831E-6;
  real Ph = .0307416 M;
  real GLY = .112 M;
  real G1P(time) M;
  real G6P(time) M;
  real F6P(time) M;
  real atp = .00779842 M;
  real FDP(time) M;
  real adp = 4.0507956E-4 M;
  real DHAP(time) M;
  real GAP(time) M;
  real NAD(time) M;
  real DPG(time) M;
  real NADH(time) M;
  real P3G(time) M;
  real P2G(time) M;
  real PEP(time) M;
  real PYR(time) M;
  real LAC = .0020513 M;
  real vgp(time) katal;
  real vpglm(time) katal;
  real vpgi(time) katal;
  real vpfk(time) katal;
  real vald(time) katal;
  real vtpi(time) katal;
  real vgapdh(time) katal;
  real vpgk(time) katal;
  real vpgm(time) katal;
  real ven(time) katal;
  real vpk(time) katal;
  real vldh(time) katal;

  // equations
  when (time=time.min) G1P = 5E-5;
  G1P:time = (vgp + -1*vpglm)/compartment;
  when (time=time.min) G6P = 7.5E-4;
  G6P:time = (vpglm + -1*vpgi)/compartment;
  when (time=time.min) F6P = 2E-4;
  F6P:time = (vpgi + -1*vpfk)/compartment;
  when (time=time.min) FDP = 7E-5;
  FDP:time = (vpfk + -1*vald)/compartment;
  when (time=time.min) DHAP = 7E-5;
  DHAP:time = (vald + vtpi)/compartment;
  when (time=time.min) GAP = 3E-5;
  GAP:time = (vald + -1*vtpi + -1*vgapdh)/compartment;
  when (time=time.min) NAD = 5E-4;
  NAD:time = (-1*vgapdh + vldh)/compartment;
  when (time=time.min) DPG = 6.5E-5;
  DPG:time = (vgapdh + -1*vpgk)/compartment;
  when (time=time.min) NADH = 5E-7;
  NADH:time = (vgapdh + -1*vldh)/compartment;
  when (time=time.min) P3G = 5E-5;
  P3G:time = (vpgk + -1*vpgm)/compartment;
  when (time=time.min) P2G = 5E-6;
  P2G:time = (vpgm + -1*ven)/compartment;
  when (time=time.min) PEP = 1E-5;
  PEP:time = (ven + -1*vpk)/compartment;
  when (time=time.min) PYR = 1.44E-4;
  PYR:time = (vpk + -1*vldh)/compartment;
  vgp = fracA*(Vfgly*Ph*GLY/(KgpAigly*KgpApi)/(1+GLY/KgpAglyf+Ph/KgpApi+GLY*Ph/(KgpAglyf*KgpAipi)+GLY/KgpAglyb+G1P/KgpAg1p+GLY*G1P/(KgpAig1p*KgpAglyb))-Vfgly*KgpAglyb*KgpAig1p/(KgpAigly*KgpApi*.31)*G1P*GLY/(KgpAglyb*KgpAig1p)/(1+GLY/KgpAglyf+Ph/KgpApi+GLY*Ph/(KgpAglyf*KgpAipi)+GLY/KgpAglyb+G1P/KgpAg1p+GLY*G1P/(KgpAig1p*KgpAglyb)))+fracB*(amp^nH/(Kgpamp^nH*.02)/(1+amp^nH/(Kgpamp^nH*.02))*(Vfgly*Ph*GLY/(KgpBiglyf*KgpBpi))/(1+GLY/KgpBipi+Ph/KgpBiglyf+GLY/KgpBiglyb+G1P/KgpBig1p+GLY*Ph/(KgpBiglyf*KgpBpi)+GLY*G1P/(KgpBg1p*KgpBiglyb))-amp^nH/(Kgpamp^nH*.02)/(1+amp^nH/(Kgpamp^nH*.02))*(Vfgly*KgpBg1p*KgpBiglyb/(KgpBiglyf*KgpBpi*.31)*G1P*GLY/(KgpBg1p*KgpBiglyb))/(1+GLY/KgpBipi+Ph/KgpBiglyf+GLY/KgpBiglyb+G1P/KgpBig1p+GLY*Ph/(KgpBiglyf*KgpBpi)+GLY*G1P/(KgpBg1p*KgpBiglyb)));
  vpglm = (Vfpglm*G1P/Kpglmg1p-Vfpglm*Kpglmg6p/(Kpglmg1p*16.62)*G6P/Kpglmg6p)/(1+G1P/Kpglmg1p+G6P/Kpglmg6p);
  vpgi = (Vbpgi*Kpgig6p/Kpgif6p*.45*G6P/Kpgig6p-Vbpgi*F6P/Kpgif6p)/(1+F6P/Kpgif6p+G6P/Kpgig6p);
  vpfk = Vfpfk*atp*F6P*(1+Kpfkatp*Kpfkf6p*Lo*(1+en*amp/Kpfkamp)^4*(1+atp/Kpfkiatp)^4*((1+atp/KpfkatpT)*(1+F6P/Kpfkf6pT)+FDP/KpfkfdpT+adp*(1+FDP/KpfkfdpT)/KpfkadpT)^3/(KpfkatpT*Kpfkf6pT*(1+amp/Kpfkamp)^4*(1+dn*atp/Kpfkiatp)^4*((1+atp/Kpfkatp)*(1+F6P/Kpfkf6p)+FDP/Kpfkfdp+adp*(1+FDP/Kpfkfdp)/Kpfkadp)^3))/(Kpfkatp*Kpfkf6p*((1+atp/Kpfkatp)*(1+F6P/Kpfkf6p)+FDP/Kpfkfdp+adp*(1+FDP/Kpfkfdp)/Kpfkadp)*(1+Lo*(1+en*amp/Kpfkamp)^4*(1+atp/Kpfkiatp)^4*((1+atp/KpfkatpT)*(1+F6P/Kpfkf6pT)+FDP/KpfkfdpT+adp*(1+FDP/KpfkfdpT)/KpfkadpT)^4/((1+amp/Kpfkamp)^4*(1+dn*atp/Kpfkiatp)^4*((1+atp/Kpfkatp)*(1+F6P/Kpfkf6p)+FDP/Kpfkfdp+adp*(1+FDP/Kpfkfdp)/Kpfkadp)^4)))-.00411742907728414*Vfpfk*adp*FDP*(1+Kpfkatp*Kpfkf6p*Lo*(1+en*amp/Kpfkamp)^4*(1+atp/Kpfkiatp)^4*((1+atp/KpfkatpT)*(1+F6P/Kpfkf6pT)+FDP/KpfkfdpT+adp*(1+FDP/KpfkfdpT)/KpfkadpT)^3/(KpfkatpT*Kpfkf6pT*(1+amp/Kpfkamp)^4*(1+dn*atp/Kpfkiatp)^4*((1+atp/Kpfkatp)*(1+F6P/Kpfkf6p)+FDP/Kpfkfdp+adp*(1+FDP/Kpfkfdp)/Kpfkadp)^3))/(Kpfkatp*Kpfkf6p*((1+atp/Kpfkatp)*(1+F6P/Kpfkf6p)+FDP/Kpfkfdp+adp*(1+FDP/Kpfkfdp)/Kpfkadp)*(1+Lo*(1+en*amp/Kpfkamp)^4*(1+atp/Kpfkiatp)^4*((1+atp/KpfkatpT)*(1+F6P/Kpfkf6pT)+FDP/KpfkfdpT+adp*(1+FDP/KpfkfdpT)/KpfkadpT)^4/((1+amp/Kpfkamp)^4*(1+dn*atp/Kpfkiatp)^4*((1+atp/Kpfkatp)*(1+F6P/Kpfkf6p)+FDP/Kpfkfdp+adp*(1+FDP/Kpfkfdp)/Kpfkadp)^4)));
  vald = (Vfald*FDP/Kaldfdp-Vfald*Kaldgap*Kalddhap/(Kaldfdp*9.499999999999999E-5)*DHAP*GAP/(Kaldgap*Kalddhap))/(1+FDP/Kaldfdp+GAP/Kaldgap+DHAP/Kalddhap);
  vtpi = (Vftpi*GAP/Ktpigap-Vftpi*Ktpidhap/(Ktpigap*19.2)*DHAP/Ktpidhap)/(1+GAP/Ktpigap+DHAP/Ktpidhap);
  vgapdh = ((-11.2359550561798)*Vfgad*DPG*NADH/(Kgapdhgap*Kgapdhnad*Kgapdhpi)+Vfgad*GAP*NAD*Ph/(Kgapdhgap*Kgapdhnad*Kgapdhpi))/(1+DPG/Kgapdh13dpg+GAP/Kgapdhgap+NAD/Kgapdhnad+GAP*NAD/(Kgapdhgap*Kgapdhnad)+NADH/Kgapdhnadh+DPG*NADH/(Kgapdh13dpg*Kgapdhnadh)+Ph/Kgapdhpi+GAP*NAD*Ph/(Kgapdhgap*Kgapdhnad*Kgapdhpi));
  vpgk = (57109*Vbpgk*adp*DPG/(Kpgk3pg*Kpgkatp)-Vbpgk*atp*P3G/(Kpgk3pg*Kpgkatp))/(1+adp/Kpgkadp+atp/Kpgkatp+DPG/Kpgk13dpg+adp*DPG/(Kpgk13dpg*Kpgkadp)+P3G/Kpgk3pg+atp*P3G/(Kpgk3pg*Kpgkatp));
  vpgm = (Vfpgm*P3G/Kpgm3pg-Vfpgm*Kpgm2pg/(Kpgm3pg*.18)*P2G/Kpgm2pg)/(1+P3G/Kpgm3pg+P2G/Kpgm2pg);
  ven = (Vfen*P2G/Ken2pg-Vfen*Kenpep/(Ken2pg*.49)*PEP/Kenpep)/(1+P2G/Ken2pg+PEP/Kenpep);
  vpk = (Vfpk*PEP*adp/(Kpkpep*Kpkadp)-Vfpk*Kpkpyr*Kpkatp/(Kpkpep*Kpkadp*10304)*PYR*atp/(Kpkpyr*Kpkatp))/(1+PEP/Kpkpep+adp/Kpkadp+PEP*adp/(Kpkpep*Kpkadp)+atp/Kpkatp+PYR/Kpkpyr+PYR*atp/(Kpkpyr*Kpkatp));
  vldh = (Vfldh*PYR*NADH/(Kldhpyr*Kldhnadh)-Vfldh*Kldhlac*Kldhnad/(Kldhpyr*Kldhnadh*16198)*LAC*NAD/(Kldhlac*Kldhnad))/(1+PYR/Kldhpyr+NADH/Kldhnadh+PYR*NADH/(Kldhpyr*Kldhnadh)+LAC/Kldhlac+NAD/Kldhnad+LAC*NAD/(Kldhlac*Kldhnad));

  // variable properties
  compartment.sbmlRole="compartment";
  Vfgly.sbmlRole="parameter";
  KgpAigly.sbmlRole="parameter";
  KgpAipi.sbmlRole="parameter";
  KgpApi.sbmlRole="parameter";
  KgpAglyf.sbmlRole="parameter";
  KgpAglyb.sbmlRole="parameter";
  KgpAg1p.sbmlRole="parameter";
  KgpAig1p.sbmlRole="parameter";
  KgpBipi.sbmlRole="parameter";
  KgpBiglyf.sbmlRole="parameter";
  KgpBpi.sbmlRole="parameter";
  KgpBg1p.sbmlRole="parameter";
  KgpBiglyb.sbmlRole="parameter";
  KgpBig1p.sbmlRole="parameter";
  Kgpamp.sbmlRole="parameter";
  fracA.sbmlRole="parameter";
  fracB.sbmlRole="parameter";
  nH.sbmlRole="parameter";
  Vfpglm.sbmlRole="parameter";
  Kpglmg1p.sbmlRole="parameter";
  Kpglmg6p.sbmlRole="parameter";
  Vbpgi.sbmlRole="parameter";
  Kpgig6p.sbmlRole="parameter";
  Kpgif6p.sbmlRole="parameter";
  Vfpfk.sbmlRole="parameter";
  Kpfkfdp.sbmlRole="parameter";
  KpfkfdpT.sbmlRole="parameter";
  Kpfkadp.sbmlRole="parameter";
  KpfkadpT.sbmlRole="parameter";
  Kpfkf6p.sbmlRole="parameter";
  Kpfkf6pT.sbmlRole="parameter";
  Kpfkatp.sbmlRole="parameter";
  KpfkatpT.sbmlRole="parameter";
  Kpfkamp.sbmlRole="parameter";
  Kpfkiatp.sbmlRole="parameter";
  Vfald.sbmlRole="parameter";
  Kaldfdp.sbmlRole="parameter";
  Kaldgap.sbmlRole="parameter";
  Kalddhap.sbmlRole="parameter";
  Vftpi.sbmlRole="parameter";
  Ktpigap.sbmlRole="parameter";
  Ktpidhap.sbmlRole="parameter";
  Vfgad.sbmlRole="parameter";
  Kgapdhgap.sbmlRole="parameter";
  Kgapdhnad.sbmlRole="parameter";
  Kgapdhpi.sbmlRole="parameter";
  Kgapdhnadh.sbmlRole="parameter";
  Kgapdh13dpg.sbmlRole="parameter";
  Vbpgk.sbmlRole="parameter";
  Kpgk13dpg.sbmlRole="parameter";
  Kpgkadp.sbmlRole="parameter";
  Kpgk3pg.sbmlRole="parameter";
  Kpgkatp.sbmlRole="parameter";
  Vfpgm.sbmlRole="parameter";
  Kpgm3pg.sbmlRole="parameter";
  Kpgm2pg.sbmlRole="parameter";
  Vfen.sbmlRole="parameter";
  Ken2pg.sbmlRole="parameter";
  Kenpep.sbmlRole="parameter";
  Vfpk.sbmlRole="parameter";
  Kpkpep.sbmlRole="parameter";
  Kpkpyr.sbmlRole="parameter";
  Kpkadp.sbmlRole="parameter";
  Kpkatp.sbmlRole="parameter";
  Vfldh.sbmlRole="parameter";
  Kldhpyr.sbmlRole="parameter";
  Kldhnadh.sbmlRole="parameter";
  Kldhlac.sbmlRole="parameter";
  Kldhnad.sbmlRole="parameter";
  VrevCK.sbmlRole="parameter";
  KckPCr.sbmlRole="parameter";
  Kckiatp.sbmlRole="parameter";
  Kckcr.sbmlRole="parameter";
  Kckiadp.sbmlRole="parameter";
  KckiPCr.sbmlRole="parameter";
  KeqCK.sbmlRole="parameter";
  Vfadk.sbmlRole="parameter";
  Kadkadp.sbmlRole="parameter";
  Kadkamp.sbmlRole="parameter";
  Kadkatp.sbmlRole="parameter";
  k.sbmlRole="parameter";
  kout.sbmlRole="parameter";
  dn.sbmlRole="parameter";
  en.sbmlRole="parameter";
  Lo.sbmlRole="parameter";
  amp.sbmlRole="parameter";
  Ph.sbmlRole="species";
  Ph.sbmlCompartment="compartment";
  GLY.sbmlRole="species";
  GLY.sbmlCompartment="compartment";
  G1P.sbmlRole="species";
  G1P.sbmlCompartment="compartment";
  G6P.sbmlRole="species";
  G6P.sbmlCompartment="compartment";
  F6P.sbmlRole="species";
  F6P.sbmlCompartment="compartment";
  atp.sbmlRole="species";
  atp.sbmlCompartment="compartment";
  FDP.sbmlRole="species";
  FDP.sbmlCompartment="compartment";
  adp.sbmlRole="species";
  adp.sbmlCompartment="compartment";
  DHAP.sbmlRole="species";
  DHAP.sbmlCompartment="compartment";
  GAP.sbmlRole="species";
  GAP.sbmlCompartment="compartment";
  NAD.sbmlRole="species";
  NAD.sbmlCompartment="compartment";
  DPG.sbmlRole="species";
  DPG.sbmlCompartment="compartment";
  NADH.sbmlRole="species";
  NADH.sbmlCompartment="compartment";
  P3G.sbmlRole="species";
  P3G.sbmlCompartment="compartment";
  P2G.sbmlRole="species";
  P2G.sbmlCompartment="compartment";
  PEP.sbmlRole="species";
  PEP.sbmlCompartment="compartment";
  PYR.sbmlRole="species";
  PYR.sbmlCompartment="compartment";
  LAC.sbmlRole="species";
  LAC.sbmlCompartment="compartment";
  vgp.sbmlRole="rate";
  vpglm.sbmlRole="rate";
  vpgi.sbmlRole="rate";
  vpfk.sbmlRole="rate";
  vald.sbmlRole="rate";
  vtpi.sbmlRole="rate";
  vgapdh.sbmlRole="rate";
  vpgk.sbmlRole="rate";
  vpgm.sbmlRole="rate";
  ven.sbmlRole="rate";
  vpk.sbmlRole="rate";
  vldh.sbmlRole="rate";
}