// 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: Glu ATP => G6P ADP 
// v2: G6P => F6P 
// v3: F6P ATP => FBP ADP 
// v4: FBP => DHAP GAP 
// v5: GAP => DHAP 
// v6: GAP NAD P => G13DP NADH H 
// v7: G13DP ADP => P3G ATP 
// v8: P3G => P2G 
// v9: P2G => PEP 
// v10: PEP ADP => Pyr ATP P 
// v11: Pyr NADH => Lac NAD 
// v12: G13DP => G23DP 
// v13: G23DP => P3G P 

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

  // variable definitions
  real compartment = 1 L;
  real KmHK = 50;
  real Km1HK = 1E3;
  real Km2HK = 50;
  real Km3HK = 1E3;
  real KeHK = 7.75E6;
  real VmaxHK = 186;
  real KmGPI = 200;
  real Km1GPI = 80;
  real KeGPI = .329;
  real VmaxGPI = 9E3;
  real KmPFK = 100;
  real Km1PFK = 100;
  real Km2PFK = 600;
  real Km3PFK = 600;
  real Km4PFK = 4E3;
  real Km5PFK = 3E4;
  real KePFK = 1.44E7;
  real nPFK = 5;
  real VmaxPFK = 5190;
  real KmALD = 8;
  real Km1ALD = 10;
  real Km2ALD = 40;
  real Km3ALD = 40;
  real Km4ALD = 15;
  real KeALD = 140;
  real VmaxALD = 1380;
  real KmTPI = 650;
  real Km1TPI = 1300;
  real KeTPI = 9;
  real VmaxTPI = 525000;
  real KmGAPDH = 45;
  real Km1GAPDH = 3200;
  real Km2GAPDH = 1.6E-14;
  real Km3GAPDH = 1.6E-14;
  real Km4GAPDH = 21;
  real Km5GAPDH = 31;
  real Km6GAPDH = .67;
  real Km7GAPDH = 100;
  real Km8GAPDH = 45;
  real Km9GAPDH = 1;
  real Km10GAPDH = 64;
  real KeGAPDH = 1.34;
  real VmaxGAPDH = 39000;
  real KmPGK = 80;
  real Km1PGK = 2;
  real Km2PGK = 200;
  real Km3PGK = 1100;
  real Km4PGK = 1600;
  real Km5PGK = 205;
  real KePGK = 2.5;
  real VmaxPGK = 87000;
  real KmPGM = 168;
  real Km1PGM = 26;
  real KePGM = .17;
  real VmaxPGM = 6600;
  real KmENOLASE = 46;
  real Km1ENOLASE = 140;
  real Km2ENOLASE = 111;
  real Km3ENOLASE = 46;
  real Km4ENOLASE = 140;
  real KeENOLASE = 4;
  real VmaxENOLASE = 3E3;
  real KmPK = 225;
  real Km1PK = 474;
  real Km2PK = 2E3;
  real Km3PK = 3E3;
  real Km4PK = 3400;
  real Km5PK = 255;
  real Km6PK = 5;
  real KePK = 425;
  real nPK = 1;
  real VmaxPK = 4200;
  real KmLDH = 2.45;
  real Km1LDH = 137;
  real Km2LDH = 503;
  real Km3LDH = 1070;
  real Km4LDH = 8.44;
  real Km5LDH = 107;
  real Km6LDH = 101;
  real Km7LDH = 7330;
  real Km8LDH = 228;
  real KeLDH = 4;
  real VmaxLDH = 55500;
  real KmBPGM = 40;
  real Km1BPGM = .013;
  real VmaxBPGM = 51000;
  real KmBPGP = 20;
  real Km1BPGP = .006;
  real VmaxBPGP = 8.1;
  real Mg = 300;
  real Glu = 5E3 M;
  real ATP = 3E3 M;
  real G6P(time) M;
  real ADP = 300 M;
  real F6P(time) M;
  real FBP(time) M;
  real DHAP(time) M;
  real GAP(time) M;
  real NAD(time) M;
  real P = 1E3 M;
  real G13DP(time) M;
  real NADH(time) M;
  real H = .01 M;
  real GAPNAD = 0 M;
  real GAPNADH = 0 M;
  real P3G(time) M;
  real P2G(time) M;
  real PEP(time) M;
  real Pyr = 60 M;
  real Lac = 1200 M;
  real G23DP = 4500 M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;
  real v5(time) katal;
  real v6(time) katal;
  real v7(time) katal;
  real v8(time) katal;
  real v9(time) katal;
  real v10(time) katal;
  real v11(time) katal;
  real v12(time) katal;
  real v13(time) katal;

  // equations
  when (time=time.min) G6P = 0;
  G6P:time = (v1 + -1*v2)/compartment;
  when (time=time.min) F6P = 0;
  F6P:time = (v2 + -1*v3)/compartment;
  when (time=time.min) FBP = 0;
  FBP:time = (v3 + -1*v4)/compartment;
  when (time=time.min) DHAP = 0;
  DHAP:time = (v4 + v5)/compartment;
  when (time=time.min) GAP = 0;
  GAP:time = (v4 + -1*v5 + -1*v6)/compartment;
  when (time=time.min) NAD = 60;
  NAD:time = (-1*v6 + v11)/compartment;
  when (time=time.min) G13DP = 0;
  G13DP:time = (v6 + -1*v7 + -1*v12)/compartment;
  when (time=time.min) NADH = 1;
  NADH:time = (v6 + -1*v11)/compartment;
  when (time=time.min) P3G = 0;
  P3G:time = (v7 + -1*v8 + v13)/compartment;
  when (time=time.min) P2G = 0;
  P2G:time = (v8 + -1*v9)/compartment;
  when (time=time.min) PEP = 0;
  PEP:time = (v9 + -1*v10)/compartment;
  v1 = VmaxHK*(Glu*ATP/(KmHK*Km1HK)-G6P*ADP/KeHK)/(1+Glu/KmHK+ATP/Km1HK+G6P/Km2HK+ADP/Km3HK+Glu*ATP/(KmHK*Km1HK)+G6P*ADP/(Km2HK*Km3HK));
  v2 = VmaxGPI/KmGPI*(G6P-F6P/KeGPI)/(1+G6P/KmGPI+F6P/Km1GPI);
  v3 = 1/(1+.1122^nPFK*(1+ATP/Km1PFK)^4*(1+Mg/Km4PFK)^4/((1+F6P/KmPFK+FBP/Km2PFK)^4*(1+P/Km5PFK)^4))*VmaxPFK*(F6P*ATP/(KmPFK*Km1PFK)-FBP*ADP/KePFK)/(1+F6P/KmPFK+ATP/Km1PFK+FBP/Km2PFK+ADP/Km3PFK+F6P*ATP/(KmPFK*Km1PFK)+FBP*ADP/(Km2PFK*Km3PFK));
  v4 = VmaxALD*(FBP/KmALD-DHAP*GAP/KeALD)/(1+FBP/KmALD+DHAP/Km1ALD+(Km3ALD*GAP+DHAP*GAP)/(Km1ALD*Km2ALD)+Km3ALD*FBP*GAP/(Km1ALD*Km2ALD*Km4ALD));
  v5 = VmaxTPI/KmTPI*(GAP-DHAP/KeTPI)/(1+GAP/KmTPI+DHAP/Km1TPI);
  v6 = VmaxGAPDH*(GAP*NAD*P/(KmGAPDH*Km1GAPDH*Km2GAPDH)-G13DP*NADH*H/(KeGAPDH*Km3GAPDH*Km4GAPDH))/((GAP/Km2GAPDH+G13DP/Km3GAPDH+GAP*P/(Km1GAPDH*Km2GAPDH))*(1+GAP/Km5GAPDH)+Km6GAPDH*NADH*H/(Km3GAPDH*Km4GAPDH)+Km7GAPDH*NAD*P/(KmGAPDH*Km1GAPDH*Km2GAPDH)+NAD*GAP/(Km8GAPDH*Km2GAPDH)+NAD*G13DP/(Km3GAPDH*Km8GAPDH)+GAPNAD*P/(KmGAPDH*Km1GAPDH*Km2GAPDH)+Km7GAPDH*NAD*P*G13DP/(KmGAPDH*Km2GAPDH*Km1GAPDH*Km9GAPDH)+NADH*G13DP*H/(Km3GAPDH*Km4GAPDH)+GAPNADH*H/(Km2GAPDH*Km10GAPDH)+Km6GAPDH*NADH*P*H/(Km1GAPDH*Km3GAPDH*Km4GAPDH)+GAP*NADH*H*P/(Km1GAPDH*Km2GAPDH*Km10GAPDH)+Km6GAPDH*NADH*P*H*G13DP/(Km9GAPDH*Km1GAPDH*Km3GAPDH*Km4GAPDH));
  v7 = VmaxPGK*((G13DP*ADP/(KmPGK*Km1PGK)-P3G*ATP/(KePGK*Km2PGK*Km3PGK))/(1+G13DP/Km4PGK+ADP/KmPGK+P3G/Km5PGK+ATP/Km2PGK+G13DP*ADP/(KmPGK*Km1PGK)+P3G*ATP/(Km2PGK*Km3PGK)));
  v8 = VmaxPGM/KmPGM*(P3G-P2G/KePGM)/(1+P3G/KmPGM+P2G/Km1PGM);
  v9 = VmaxENOLASE*(P2G*Mg/(KmENOLASE*Km1ENOLASE)-PEP*Mg/(KeENOLASE*Km2ENOLASE*Km3ENOLASE))/(1+P2G/Km4ENOLASE+Mg/KmENOLASE+PEP/Km2ENOLASE+P2G*Mg/(KmENOLASE*Km1ENOLASE)+PEP*Mg/(Km2ENOLASE*Km3ENOLASE));
  v10 = 1/(1+.0631^nPK*(1+ATP/Km4PK)^4/((1+PEP/Km5PK+Pyr/Km2PK)^4*(1+FBP/Km6PK)^4))*VmaxPK*(PEP*ADP/(KmPK*Km1PK)-Pyr*ATP/(KePK*Km2PK*Km3PK))/(1+PEP/KmPK+ADP/Km1PK+Pyr/Km2PK+ATP/Km3PK+PEP*ADP/(KmPK*Km1PK)+Pyr*ATP/(Km2PK*Km3PK));
  v11 = VmaxLDH*((NADH*Pyr/(KmLDH*Km1LDH)-NAD*Lac/(KeLDH*Km2LDH*Km3LDH))/(NADH/KmLDH+NAD/Km3LDH+NADH*Pyr/(KmLDH*Km1LDH)+NAD*Lac/(Km2LDH*Km3LDH)+(1+Km4LDH*Pyr/(KmLDH*Km1LDH)+Km5LDH*Lac/(Km2LDH*Km3LDH))*(1+Pyr/Km6LDH)+NADH*Lac*Km5LDH/(KmLDH*Km2LDH*Km3LDH)+NAD*Pyr*Km4LDH/(KmLDH*Km1LDH*Km3LDH)+NADH*Lac*Pyr/(KmLDH*Km1LDH*Km7LDH)+NAD*Pyr*Lac/(Km8LDH*Km2LDH*Km3LDH)));
  v12 = VmaxBPGM*G13DP/(KmBPGM+Km1BPGM*G13DP+G23DP);
  v13 = VmaxBPGP*G23DP/(KmBPGP*(1+(P2G+P3G)/Km1BPGP)+G23DP);

  // variable properties
  compartment.sbmlRole="compartment";
  KmHK.sbmlRole="parameter";
  Km1HK.sbmlRole="parameter";
  Km2HK.sbmlRole="parameter";
  Km3HK.sbmlRole="parameter";
  KeHK.sbmlRole="parameter";
  VmaxHK.sbmlRole="parameter";
  KmGPI.sbmlRole="parameter";
  Km1GPI.sbmlRole="parameter";
  KeGPI.sbmlRole="parameter";
  VmaxGPI.sbmlRole="parameter";
  KmPFK.sbmlRole="parameter";
  Km1PFK.sbmlRole="parameter";
  Km2PFK.sbmlRole="parameter";
  Km3PFK.sbmlRole="parameter";
  Km4PFK.sbmlRole="parameter";
  Km5PFK.sbmlRole="parameter";
  KePFK.sbmlRole="parameter";
  nPFK.sbmlRole="parameter";
  VmaxPFK.sbmlRole="parameter";
  KmALD.sbmlRole="parameter";
  Km1ALD.sbmlRole="parameter";
  Km2ALD.sbmlRole="parameter";
  Km3ALD.sbmlRole="parameter";
  Km4ALD.sbmlRole="parameter";
  KeALD.sbmlRole="parameter";
  VmaxALD.sbmlRole="parameter";
  KmTPI.sbmlRole="parameter";
  Km1TPI.sbmlRole="parameter";
  KeTPI.sbmlRole="parameter";
  VmaxTPI.sbmlRole="parameter";
  KmGAPDH.sbmlRole="parameter";
  Km1GAPDH.sbmlRole="parameter";
  Km2GAPDH.sbmlRole="parameter";
  Km3GAPDH.sbmlRole="parameter";
  Km4GAPDH.sbmlRole="parameter";
  Km5GAPDH.sbmlRole="parameter";
  Km6GAPDH.sbmlRole="parameter";
  Km7GAPDH.sbmlRole="parameter";
  Km8GAPDH.sbmlRole="parameter";
  Km9GAPDH.sbmlRole="parameter";
  Km10GAPDH.sbmlRole="parameter";
  KeGAPDH.sbmlRole="parameter";
  VmaxGAPDH.sbmlRole="parameter";
  KmPGK.sbmlRole="parameter";
  Km1PGK.sbmlRole="parameter";
  Km2PGK.sbmlRole="parameter";
  Km3PGK.sbmlRole="parameter";
  Km4PGK.sbmlRole="parameter";
  Km5PGK.sbmlRole="parameter";
  KePGK.sbmlRole="parameter";
  VmaxPGK.sbmlRole="parameter";
  KmPGM.sbmlRole="parameter";
  Km1PGM.sbmlRole="parameter";
  KePGM.sbmlRole="parameter";
  VmaxPGM.sbmlRole="parameter";
  KmENOLASE.sbmlRole="parameter";
  Km1ENOLASE.sbmlRole="parameter";
  Km2ENOLASE.sbmlRole="parameter";
  Km3ENOLASE.sbmlRole="parameter";
  Km4ENOLASE.sbmlRole="parameter";
  KeENOLASE.sbmlRole="parameter";
  VmaxENOLASE.sbmlRole="parameter";
  KmPK.sbmlRole="parameter";
  Km1PK.sbmlRole="parameter";
  Km2PK.sbmlRole="parameter";
  Km3PK.sbmlRole="parameter";
  Km4PK.sbmlRole="parameter";
  Km5PK.sbmlRole="parameter";
  Km6PK.sbmlRole="parameter";
  KePK.sbmlRole="parameter";
  nPK.sbmlRole="parameter";
  VmaxPK.sbmlRole="parameter";
  KmLDH.sbmlRole="parameter";
  Km1LDH.sbmlRole="parameter";
  Km2LDH.sbmlRole="parameter";
  Km3LDH.sbmlRole="parameter";
  Km4LDH.sbmlRole="parameter";
  Km5LDH.sbmlRole="parameter";
  Km6LDH.sbmlRole="parameter";
  Km7LDH.sbmlRole="parameter";
  Km8LDH.sbmlRole="parameter";
  KeLDH.sbmlRole="parameter";
  VmaxLDH.sbmlRole="parameter";
  KmBPGM.sbmlRole="parameter";
  Km1BPGM.sbmlRole="parameter";
  VmaxBPGM.sbmlRole="parameter";
  KmBPGP.sbmlRole="parameter";
  Km1BPGP.sbmlRole="parameter";
  VmaxBPGP.sbmlRole="parameter";
  Mg.sbmlRole="parameter";
  Glu.sbmlRole="species";
  Glu.sbmlCompartment="compartment";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  G6P.sbmlRole="species";
  G6P.sbmlCompartment="compartment";
  ADP.sbmlRole="species";
  ADP.sbmlCompartment="compartment";
  F6P.sbmlRole="species";
  F6P.sbmlCompartment="compartment";
  FBP.sbmlRole="species";
  FBP.sbmlCompartment="compartment";
  DHAP.sbmlRole="species";
  DHAP.sbmlCompartment="compartment";
  GAP.sbmlRole="species";
  GAP.sbmlCompartment="compartment";
  NAD.sbmlRole="species";
  NAD.sbmlCompartment="compartment";
  P.sbmlRole="species";
  P.sbmlCompartment="compartment";
  G13DP.sbmlRole="species";
  G13DP.sbmlCompartment="compartment";
  NADH.sbmlRole="species";
  NADH.sbmlCompartment="compartment";
  H.sbmlRole="species";
  H.sbmlCompartment="compartment";
  GAPNAD.sbmlRole="species";
  GAPNAD.sbmlCompartment="compartment";
  GAPNADH.sbmlRole="species";
  GAPNADH.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";
  G23DP.sbmlRole="species";
  G23DP.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  v5.sbmlRole="rate";
  v6.sbmlRole="rate";
  v7.sbmlRole="rate";
  v8.sbmlRole="rate";
  v9.sbmlRole="rate";
  v10.sbmlRole="rate";
  v11.sbmlRole="rate";
  v12.sbmlRole="rate";
  v13.sbmlRole="rate";
}