// 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: Glc ATP => F6P ADP 
// v2: F6P ATP => FBP ADP 
// v3: FBP => 2GAP 
// v4: GAP NAD => BPG NADH 
// v5: BPG ADP => PEP ATP 
// v6: PEP ADP => Pyr ATP 
// v7: Pyr => ACA 
// v8: ACA NADH => EtOH NAD 
// v9: AMP ATP => 2ADP 
// v10: F6P => P 
// v11: ATP => P 
// v12: ADP => P 
// v13: AMP => P 
// v14: NADH => P 
// v15: NAD => P 
// v16: Glc => P 
// v17: F6P => P 
// v18: FBP => P 
// v19: PEP => P 
// v20: GAP => P 
// v21: BPG => P 
// v22: Pyr => P 
// v23: ACA => P 
// v24: EtOH => P 
// v25: ATPo => ATP 
// v26: ADPo => ADP 
// v27: NADHo => NADH 
// v28: NADo => NAD 
// v29: Glco => Glc 

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

  // variable definitions
  real compartment = 1 L;
  real k0 = .011;
  real V1m = .5;
  real K1Glc = .1;
  real K1ATP = .063;
  real V2m = 1.5;
  real K2 = .0016;
  real kappa2 = .017;
  real K2mATP = .01;
  real k3f = 1;
  real k3b = 50;
  real V4m = 10;
  real K4G = 1;
  real K4N = 1;
  real k5f = 1;
  real k5r = .5;
  real V6m = 10;
  real K6PEP = .2;
  real K6ADP = .3;
  real V7m = 2;
  real K7 = .3;
  real k8f = 1;
  real k8r = 1.43E-4;
  real kf9 = 10;
  real kr9 = 10;
  real kf10 = .05;
  real Glc(time) M;
  real ATP(time) M;
  real F6P(time) M;
  real ADP(time) M;
  real FBP(time) M;
  real GAP(time) M;
  real NAD(time) M;
  real BPG(time) M;
  real NADH(time) M;
  real PEP(time) M;
  real Pyr(time) M;
  real ACA(time) M;
  real EtOH(time) M;
  real AMP(time) M;
  real P = 0 M;
  real ATPo = 3.5 M;
  real ADPo = 1.1 M;
  real NADHo = .24 M;
  real NADo = 4 M;
  real Glco = 50 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;
  real v14(time) katal;
  real v15(time) katal;
  real v16(time) katal;
  real v17(time) katal;
  real v18(time) katal;
  real v19(time) katal;
  real v20(time) katal;
  real v21(time) katal;
  real v22(time) katal;
  real v23(time) katal;
  real v24(time) katal;
  real v25(time) katal;
  real v26(time) katal;
  real v27(time) katal;
  real v28(time) katal;
  real v29(time) katal;

  // equations
  when (time=time.min) Glc = .0112817;
  Glc:time = (-1*v1 + -1*v16 + v29)/compartment;
  when (time=time.min) ATP = 4.49064;
  ATP:time = (-1*v1 + -1*v2 + v5 + v6 + -1*v9 + -1*v11 + v25)/compartment;
  when (time=time.min) F6P = .65939;
  F6P:time = (v1 + -1*v2 + -1*v10 + -1*v17)/compartment;
  when (time=time.min) ADP = .108367;
  ADP:time = (v1 + v2 + -1*v5 + -1*v6 + 2*v9 + -1*v12 + v26)/compartment;
  when (time=time.min) FBP = .00770135;
  FBP:time = (v2 + -1*v3 + -1*v18)/compartment;
  when (time=time.min) GAP = .00190919;
  GAP:time = (2*v3 + -1*v4 + -1*v20)/compartment;
  when (time=time.min) NAD = 3.62057;
  NAD:time = (-1*v4 + v8 + -1*v15 + v28)/compartment;
  when (time=time.min) BPG = .299109;
  BPG:time = (v4 + -1*v5 + -1*v21)/compartment;
  when (time=time.min) NADH = .616118;
  NADH:time = (v4 + -1*v8 + -1*v14 + v27)/compartment;
  when (time=time.min) PEP = .00211259;
  PEP:time = (v5 + -1*v6 + -1*v19)/compartment;
  when (time=time.min) Pyr = .00422702;
  Pyr:time = (v6 + -1*v7 + -1*v22)/compartment;
  when (time=time.min) ACA = .0738334;
  ACA:time = (v7 + -1*v8 + -1*v23)/compartment;
  when (time=time.min) EtOH = .339813;
  EtOH:time = (v8 + -1*v24)/compartment;
  when (time=time.min) AMP = .00261149;
  AMP:time = (-1*v9 + -1*v13)/compartment;
  v1 = V1m*ATP*Glc/(K1Glc+Glc)/(K1ATP+ATP);
  v2 = V2m*F6P^2/(K2+K2*kappa2*ATP^2/AMP^2+F6P^2)*(ATP/(K2mATP+ATP));
  v3 = k3f*FBP-k3b*GAP^2;
  v4 = V4m*GAP*NAD/(K4G+GAP)/(K4N+NAD);
  v5 = k5f*BPG*ADP-k5r*PEP*ATP;
  v6 = V6m*ADP*PEP/(K6PEP+PEP)/(K6ADP+ADP);
  v7 = V7m*Pyr/(K7+Pyr);
  v8 = k8f*ACA*NADH-k8r*EtOH*NAD;
  v9 = kf9*AMP*ATP-kr9*ADP^2;
  v10 = kf10*F6P;
  v11 = k0*ATP;
  v12 = k0*ADP;
  v13 = k0*AMP;
  v14 = k0*NADH;
  v15 = k0*NAD;
  v16 = k0*Glc;
  v17 = k0*F6P;
  v18 = k0*FBP;
  v19 = k0*PEP;
  v20 = k0*GAP;
  v21 = k0*BPG;
  v22 = k0*Pyr;
  v23 = k0*ACA;
  v24 = k0*EtOH;
  v25 = k0*ATPo;
  v26 = k0*ADPo;
  v27 = k0*NADHo;
  v28 = k0*NADo;
  v29 = k0*Glco;

  // variable properties
  compartment.sbmlRole="compartment";
  k0.sbmlRole="parameter";
  V1m.sbmlRole="parameter";
  K1Glc.sbmlRole="parameter";
  K1ATP.sbmlRole="parameter";
  V2m.sbmlRole="parameter";
  K2.sbmlRole="parameter";
  kappa2.sbmlRole="parameter";
  K2mATP.sbmlRole="parameter";
  k3f.sbmlRole="parameter";
  k3b.sbmlRole="parameter";
  V4m.sbmlRole="parameter";
  K4G.sbmlRole="parameter";
  K4N.sbmlRole="parameter";
  k5f.sbmlRole="parameter";
  k5r.sbmlRole="parameter";
  V6m.sbmlRole="parameter";
  K6PEP.sbmlRole="parameter";
  K6ADP.sbmlRole="parameter";
  V7m.sbmlRole="parameter";
  K7.sbmlRole="parameter";
  k8f.sbmlRole="parameter";
  k8r.sbmlRole="parameter";
  kf9.sbmlRole="parameter";
  kr9.sbmlRole="parameter";
  kf10.sbmlRole="parameter";
  Glc.sbmlRole="species";
  Glc.sbmlCompartment="compartment";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  F6P.sbmlRole="species";
  F6P.sbmlCompartment="compartment";
  ADP.sbmlRole="species";
  ADP.sbmlCompartment="compartment";
  FBP.sbmlRole="species";
  FBP.sbmlCompartment="compartment";
  GAP.sbmlRole="species";
  GAP.sbmlCompartment="compartment";
  NAD.sbmlRole="species";
  NAD.sbmlCompartment="compartment";
  BPG.sbmlRole="species";
  BPG.sbmlCompartment="compartment";
  NADH.sbmlRole="species";
  NADH.sbmlCompartment="compartment";
  PEP.sbmlRole="species";
  PEP.sbmlCompartment="compartment";
  Pyr.sbmlRole="species";
  Pyr.sbmlCompartment="compartment";
  ACA.sbmlRole="species";
  ACA.sbmlCompartment="compartment";
  EtOH.sbmlRole="species";
  EtOH.sbmlCompartment="compartment";
  AMP.sbmlRole="species";
  AMP.sbmlCompartment="compartment";
  P.sbmlRole="species";
  P.sbmlCompartment="compartment";
  ATPo.sbmlRole="species";
  ATPo.sbmlCompartment="compartment";
  ADPo.sbmlRole="species";
  ADPo.sbmlCompartment="compartment";
  NADHo.sbmlRole="species";
  NADHo.sbmlCompartment="compartment";
  NADo.sbmlRole="species";
  NADo.sbmlCompartment="compartment";
  Glco.sbmlRole="species";
  Glco.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";
  v14.sbmlRole="rate";
  v15.sbmlRole="rate";
  v16.sbmlRole="rate";
  v17.sbmlRole="rate";
  v18.sbmlRole="rate";
  v19.sbmlRole="rate";
  v20.sbmlRole="rate";
  v21.sbmlRole="rate";
  v22.sbmlRole="rate";
  v23.sbmlRole="rate";
  v24.sbmlRole="rate";
  v25.sbmlRole="rate";
  v26.sbmlRole="rate";
  v27.sbmlRole="rate";
  v28.sbmlRole="rate";
  v29.sbmlRole="rate";
}