// 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
// vADH: NADH acald => NAD etoh 
// vAK: 2ADP => AMP ATP 
// vALDO: FBP => G3P DHAP 
// vATPases: ATP => ADP 
// vDHases: NAD => NADH 
// vENO: n2PG => PEP 
// vGAPDH: G3P NAD => NADH n13BPG 
// vHKrev: Glu ATP => ADP G6P 
// vHPhi: G6P => F6P 
// vPFORALDH: NADH pyr => NAD acald 
// vPGAM: n3PG => n2PG 
// vPGK: ADP n13BPG => n3PG ATP 
// vPPDK: PEP AMP PPhi => pyr ATP Phi 
// vPPhiPFK: PPhi F6P => FBP Phi 
// vPPhisynthesis: ATP => AMP PPhi 
// vTPhi: DHAP => G3P 
// vglycogensynthesis: ATP G6P => ADP PPhi glycogen 
// vglycogendegradation: glycogen Phi => G6P 
// vn3PGDH: NAD n3PG => NADH n3POHpyr 

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

  // variable definitions
  real compartment = 1 L;
  real HKrevparVmf = .095;
  real HKrevparalfa = 1;
  real HKrevparKa = .025;
  real HKrevparKb = .121;
  real HKrevparKeq = 656;
  real HKrevparKi = .036;
  real HKrevparbeta = 1;
  real HKrevpargama = 1;
  real HKrevparKp = 1;
  real HKrevparKq = .235;
  real HPhiparKms = .61;
  real HPhiparKmp = .46;
  real HPhiparVf = .233;
  real HPhiparVr = .206;
  real PPhiPFKparVmf = .213;
  real PPhiPFKparalfa = .5;
  real PPhiPFKparKa = .695;
  real PPhiPFKparKb = .38;
  real PPhiPFKparVmr = .346;
  real PPhiPFKparbeta = 2;
  real PPhiPFKparKp = .109;
  real PPhiPFKparKq = 2.3;
  real ALDOparVmf = .16;
  real ALDOparKs = .028;
  real ALDOparVmr = .284;
  real ALDOparKp = .21;
  real ALDOparKq = .26;
  real TPhiparKms = .445;
  real TPhiparKmp = .32;
  real TPhiparVf = 4.366;
  real TPhiparVr = 6.098;
  real GAPDHparVmf = .405;
  real GAPDHparalfa = 1;
  real GAPDHparKa = .083;
  real GAPDHparKb = .043;
  real GAPDHparVmr = 3.68;
  real GAPDHparbeta = 1;
  real GAPDHparKp = .016;
  real GAPDHparKq = .06;
  real PGKparVmf = 3.182;
  real PGKparalfa = 1;
  real PGKparKa = .125;
  real PGKparKb = .04;
  real PGKparVmr = 1.742;
  real PGKparbeta = 1;
  real PGKparKp = .505;
  real PGKparKq = .061;
  real PGAMparKms = .83;
  real PGAMparKmp = .106;
  real PGAMparVf = .116;
  real PGAMparVr = .104;
  real ENOparKms = .06;
  real ENOparKmp = .102;
  real ENOparVf = .508;
  real ENOparVr = .103;
  real PPDKparVmf = .304;
  real PPDKparalfa = .01;
  real PPDKparKa = .03;
  real PPDKparKb = .002;
  real PPDKparKc = .091;
  real PPDKparVmr = .019;
  real PPDKparbeta = 2;
  real PPDKparKp = .3;
  real PPDKparKq = .28;
  real PPDKparKr = 1;
  real PFORALDHparVmf = .3;
  real PFORALDHparalfa = 1;
  real PFORALDHparKa = .18;
  real PFORALDHparKb = .14;
  real PFORALDHparVmr = .01;
  real PFORALDHparbeta = 1;
  real PFORALDHparKp = .1;
  real PFORALDHparKq = .2;
  real ADHparVm = .4;
  real ADHparKa = .05;
  real ADHparKb = .15;
  real glycogensynthesispark = .0015;
  real glycogendegradationpark = 1E-4;
  real n3PGDHparVm = .01;
  real n3PGDHparKa = .087;
  real n3PGDHparKb = .212;
  real ATPasesparK = .002;
  real AKpark1 = 2.26;
  real AKpark2 = 1;
  real PPhisynthesispark = .0065;
  real DHasespark1 = .1;
  real DHasespark2 = 1.5;
  real NADH(time) M;
  real acald(time) M;
  real NAD(time) M;
  real etoh = 10 M;
  real ADP(time) M;
  real AMP(time) M;
  real ATP(time) M;
  real FBP(time) M;
  real G3P(time) M;
  real DHAP(time) M;
  real n2PG(time) M;
  real PEP(time) M;
  real n13BPG(time) M;
  real Glu = 5 M;
  real G6P(time) M;
  real F6P(time) M;
  real pyr(time) M;
  real n3PG(time) M;
  real PPhi = .45 M;
  real Phi = 5 M;
  real glycogen = 1 M;
  real n3POHpyr = 1 M;
  real vADH(time) katal;
  real vAK(time) katal;
  real vALDO(time) katal;
  real vATPases(time) katal;
  real vDHases(time) katal;
  real vENO(time) katal;
  real vGAPDH(time) katal;
  real vHKrev(time) katal;
  real vHPhi(time) katal;
  real vPFORALDH(time) katal;
  real vPGAM(time) katal;
  real vPGK(time) katal;
  real vPPDK(time) katal;
  real vPPhiPFK(time) katal;
  real vPPhisynthesis(time) katal;
  real vTPhi(time) katal;
  real vglycogensynthesis(time) katal;
  real vglycogendegradation(time) katal;
  real vn3PGDH(time) katal;

  // equations
  when (time=time.min) NADH = .05;
  NADH:time = (-1*vADH + vDHases + vGAPDH + -1*vPFORALDH + vn3PGDH)/compartment;
  when (time=time.min) acald = .001;
  acald:time = (-1*vADH + vPFORALDH)/compartment;
  when (time=time.min) NAD = 1.5;
  NAD:time = (vADH + -1*vDHases + -1*vGAPDH + vPFORALDH + -1*vn3PGDH)/compartment;
  when (time=time.min) ADP = 3.3;
  ADP:time = (-2*vAK + vATPases + vHKrev + -1*vPGK + vglycogensynthesis)/compartment;
  when (time=time.min) AMP = 1.6;
  AMP:time = (vAK + -1*vPPDK + vPPhisynthesis)/compartment;
  when (time=time.min) ATP = 5;
  ATP:time = (vAK + -1*vATPases + -1*vHKrev + vPGK + vPPDK + -1*vPPhisynthesis + -1*vglycogensynthesis)/compartment;
  when (time=time.min) FBP = .43;
  FBP:time = (-1*vALDO + vPPhiPFK)/compartment;
  when (time=time.min) G3P = .36;
  G3P:time = (vALDO + -1*vGAPDH + vTPhi)/compartment;
  when (time=time.min) DHAP = 1.15;
  DHAP:time = (vALDO + -1*vTPhi)/compartment;
  when (time=time.min) n2PG = 1;
  n2PG:time = (-1*vENO + vPGAM)/compartment;
  when (time=time.min) PEP = 1;
  PEP:time = (vENO + -1*vPPDK)/compartment;
  when (time=time.min) n13BPG = 1;
  n13BPG:time = (vGAPDH + -1*vPGK)/compartment;
  when (time=time.min) G6P = 6.2;
  G6P:time = (vHKrev + -1*vHPhi + -1*vglycogensynthesis + vglycogendegradation)/compartment;
  when (time=time.min) F6P = 1.1;
  F6P:time = (vHPhi + -1*vPPhiPFK)/compartment;
  when (time=time.min) pyr = .92;
  pyr:time = (-1*vPFORALDH + vPPDK)/compartment;
  when (time=time.min) n3PG = 1;
  n3PG:time = (-1*vPGAM + vPGK + -1*vn3PGDH)/compartment;
  vADH = ADHparVm*NADH*acald/(ADHparKa*ADHparKb+ADHparKb*NADH+NADH*acald);
  vAK = AKpark1*ADP^2-AKpark2*ATP*AMP;
  vALDO = (ALDOparVmf*(FBP/ALDOparKs)-ALDOparVmr*(DHAP*G3P/(ALDOparKp+ALDOparKq)))/(1+FBP/ALDOparKs+DHAP/ALDOparKp+G3P/ALDOparKq+DHAP*G3P/(ALDOparKp*ALDOparKq));
  vATPases = ATPasesparK*ATP;
  vDHases = DHasespark1*NAD-DHasespark2*NADH;
  vENO = (ENOparVf*n2PG/ENOparKms-ENOparVr*PEP/ENOparKmp)/(1+n2PG/ENOparKms+PEP/ENOparKmp);
  vGAPDH = (GAPDHparVmf*(NAD*G3P/(GAPDHparalfa*GAPDHparKa*GAPDHparKb))-GAPDHparVmr*(n13BPG*NADH/(GAPDHparbeta*GAPDHparKp*GAPDHparKq)))/(1+NAD/GAPDHparKa+G3P/GAPDHparKb+NAD*G3P/(GAPDHparalfa*GAPDHparKa*GAPDHparKb)+n13BPG*NADH/(GAPDHparbeta*GAPDHparKp*GAPDHparKq)+n13BPG/GAPDHparKp+NADH/GAPDHparKq);
  vHKrev = HKrevparVmf/(HKrevparalfa*HKrevparKa*HKrevparKb)*(Glu*ATP-G6P*ADP/HKrevparKeq)/(1+Glu/HKrevparKa+ATP/HKrevparKb+Glu*ATP/(HKrevparalfa*HKrevparKa*HKrevparKb)+AMP/HKrevparKi+Glu*AMP/(HKrevparbeta*HKrevparKa*HKrevparKi)+AMP*ATP/(HKrevpargama*HKrevparKi*HKrevparKb)+G6P/HKrevparKp+ADP/HKrevparKq+G6P*ADP/(HKrevparalfa*HKrevparKp*HKrevparKq));
  vHPhi = (HPhiparVf*G6P/HPhiparKms-HPhiparVr*F6P/HPhiparKmp)/(1+G6P/HPhiparKms+F6P/HPhiparKmp);
  vPFORALDH = (PFORALDHparVmf*(NADH*pyr/(PFORALDHparalfa*PFORALDHparKa*PFORALDHparKb))-PFORALDHparVmr*(acald*NAD/(PFORALDHparbeta*PFORALDHparKp*PFORALDHparKq)))/(1+NADH/PFORALDHparKa+pyr/PFORALDHparKb+NADH*pyr/(PFORALDHparalfa*PFORALDHparKa*PFORALDHparKb)+acald*NAD/(PFORALDHparbeta*PFORALDHparKp*PFORALDHparKq)+acald/PFORALDHparKp+NAD/PFORALDHparKq);
  vPGAM = (PGAMparVf*n3PG/PGAMparKms-PGAMparVr*n2PG/PGAMparKmp)/(1+n3PG/PGAMparKms+n2PG/PGAMparKmp);
  vPGK = (PGKparVmf*(n13BPG*ADP/(PGKparalfa*PGKparKa*PGKparKb))-PGKparVmr*(n3PG*ATP/(PGKparbeta*PGKparKp*PGKparKq)))/(1+n13BPG/PGKparKa+ADP/PGKparKb+n13BPG*ADP/(PGKparalfa*PGKparKa*PGKparKb)+n3PG*ATP/(PGKparbeta*PGKparKp*PGKparKq)+n3PG/PGKparKp+ATP/PGKparKq);
  vPPDK = (PPDKparVmf*(PEP*AMP*PPhi/(PPDKparalfa*PPDKparKa*PPDKparKb*PPDKparKc))-PPDKparVmr*(pyr*ATP*Phi/(PPDKparbeta*PPDKparKp*PPDKparKq*PPDKparKr)))/(1+PEP/PPDKparKa+PEP*AMP/(PPDKparalfa*PPDKparKa*PPDKparKb)+PEP*AMP*PPhi/(PPDKparalfa*PPDKparbeta*PPDKparKa*PPDKparKb*PPDKparKc)+pyr*ATP*Phi/(PPDKparalfa*PPDKparbeta*PPDKparKp*PPDKparKq*PPDKparKr)+ATP*Phi/(PPDKparalfa*PPDKparKq*PPDKparKr)+Phi/PPDKparKr);
  vPPhiPFK = (PPhiPFKparVmf*(F6P*PPhi/(PPhiPFKparalfa*PPhiPFKparKa*PPhiPFKparKb))-PPhiPFKparVmr*(FBP*Phi/(PPhiPFKparbeta*PPhiPFKparKp*PPhiPFKparKq)))/(1+F6P/PPhiPFKparKa+PPhi/PPhiPFKparKb+F6P*PPhi/(PPhiPFKparalfa*PPhiPFKparKa*PPhiPFKparKb)+FBP*Phi/(PPhiPFKparbeta*PPhiPFKparKp*PPhiPFKparKq)+FBP/PPhiPFKparKp+Phi/PPhiPFKparKq);
  vPPhisynthesis = PPhisynthesispark*ATP;
  vTPhi = (TPhiparVf*DHAP/TPhiparKms-TPhiparVr*G3P/TPhiparKmp)/(1+DHAP/TPhiparKms+G3P/TPhiparKmp);
  vglycogensynthesis = glycogensynthesispark*G6P*ATP;
  vglycogendegradation = glycogendegradationpark*glycogen*Phi;
  vn3PGDH = n3PGDHparVm*NAD*n3PG/(n3PGDHparKa*n3PGDHparKb+n3PGDHparKb*NAD+NAD*n3PG);

  // variable properties
  compartment.sbmlRole="compartment";
  HKrevparVmf.sbmlRole="parameter";
  HKrevparalfa.sbmlRole="parameter";
  HKrevparKa.sbmlRole="parameter";
  HKrevparKb.sbmlRole="parameter";
  HKrevparKeq.sbmlRole="parameter";
  HKrevparKi.sbmlRole="parameter";
  HKrevparbeta.sbmlRole="parameter";
  HKrevpargama.sbmlRole="parameter";
  HKrevparKp.sbmlRole="parameter";
  HKrevparKq.sbmlRole="parameter";
  HPhiparKms.sbmlRole="parameter";
  HPhiparKmp.sbmlRole="parameter";
  HPhiparVf.sbmlRole="parameter";
  HPhiparVr.sbmlRole="parameter";
  PPhiPFKparVmf.sbmlRole="parameter";
  PPhiPFKparalfa.sbmlRole="parameter";
  PPhiPFKparKa.sbmlRole="parameter";
  PPhiPFKparKb.sbmlRole="parameter";
  PPhiPFKparVmr.sbmlRole="parameter";
  PPhiPFKparbeta.sbmlRole="parameter";
  PPhiPFKparKp.sbmlRole="parameter";
  PPhiPFKparKq.sbmlRole="parameter";
  ALDOparVmf.sbmlRole="parameter";
  ALDOparKs.sbmlRole="parameter";
  ALDOparVmr.sbmlRole="parameter";
  ALDOparKp.sbmlRole="parameter";
  ALDOparKq.sbmlRole="parameter";
  TPhiparKms.sbmlRole="parameter";
  TPhiparKmp.sbmlRole="parameter";
  TPhiparVf.sbmlRole="parameter";
  TPhiparVr.sbmlRole="parameter";
  GAPDHparVmf.sbmlRole="parameter";
  GAPDHparalfa.sbmlRole="parameter";
  GAPDHparKa.sbmlRole="parameter";
  GAPDHparKb.sbmlRole="parameter";
  GAPDHparVmr.sbmlRole="parameter";
  GAPDHparbeta.sbmlRole="parameter";
  GAPDHparKp.sbmlRole="parameter";
  GAPDHparKq.sbmlRole="parameter";
  PGKparVmf.sbmlRole="parameter";
  PGKparalfa.sbmlRole="parameter";
  PGKparKa.sbmlRole="parameter";
  PGKparKb.sbmlRole="parameter";
  PGKparVmr.sbmlRole="parameter";
  PGKparbeta.sbmlRole="parameter";
  PGKparKp.sbmlRole="parameter";
  PGKparKq.sbmlRole="parameter";
  PGAMparKms.sbmlRole="parameter";
  PGAMparKmp.sbmlRole="parameter";
  PGAMparVf.sbmlRole="parameter";
  PGAMparVr.sbmlRole="parameter";
  ENOparKms.sbmlRole="parameter";
  ENOparKmp.sbmlRole="parameter";
  ENOparVf.sbmlRole="parameter";
  ENOparVr.sbmlRole="parameter";
  PPDKparVmf.sbmlRole="parameter";
  PPDKparalfa.sbmlRole="parameter";
  PPDKparKa.sbmlRole="parameter";
  PPDKparKb.sbmlRole="parameter";
  PPDKparKc.sbmlRole="parameter";
  PPDKparVmr.sbmlRole="parameter";
  PPDKparbeta.sbmlRole="parameter";
  PPDKparKp.sbmlRole="parameter";
  PPDKparKq.sbmlRole="parameter";
  PPDKparKr.sbmlRole="parameter";
  PFORALDHparVmf.sbmlRole="parameter";
  PFORALDHparalfa.sbmlRole="parameter";
  PFORALDHparKa.sbmlRole="parameter";
  PFORALDHparKb.sbmlRole="parameter";
  PFORALDHparVmr.sbmlRole="parameter";
  PFORALDHparbeta.sbmlRole="parameter";
  PFORALDHparKp.sbmlRole="parameter";
  PFORALDHparKq.sbmlRole="parameter";
  ADHparVm.sbmlRole="parameter";
  ADHparKa.sbmlRole="parameter";
  ADHparKb.sbmlRole="parameter";
  glycogensynthesispark.sbmlRole="parameter";
  glycogendegradationpark.sbmlRole="parameter";
  n3PGDHparVm.sbmlRole="parameter";
  n3PGDHparKa.sbmlRole="parameter";
  n3PGDHparKb.sbmlRole="parameter";
  ATPasesparK.sbmlRole="parameter";
  AKpark1.sbmlRole="parameter";
  AKpark2.sbmlRole="parameter";
  PPhisynthesispark.sbmlRole="parameter";
  DHasespark1.sbmlRole="parameter";
  DHasespark2.sbmlRole="parameter";
  NADH.sbmlRole="species";
  NADH.sbmlCompartment="compartment";
  acald.sbmlRole="species";
  acald.sbmlCompartment="compartment";
  NAD.sbmlRole="species";
  NAD.sbmlCompartment="compartment";
  etoh.sbmlRole="species";
  etoh.sbmlCompartment="compartment";
  ADP.sbmlRole="species";
  ADP.sbmlCompartment="compartment";
  AMP.sbmlRole="species";
  AMP.sbmlCompartment="compartment";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  FBP.sbmlRole="species";
  FBP.sbmlCompartment="compartment";
  G3P.sbmlRole="species";
  G3P.sbmlCompartment="compartment";
  DHAP.sbmlRole="species";
  DHAP.sbmlCompartment="compartment";
  n2PG.sbmlRole="species";
  n2PG.sbmlCompartment="compartment";
  PEP.sbmlRole="species";
  PEP.sbmlCompartment="compartment";
  n13BPG.sbmlRole="species";
  n13BPG.sbmlCompartment="compartment";
  Glu.sbmlRole="species";
  Glu.sbmlCompartment="compartment";
  G6P.sbmlRole="species";
  G6P.sbmlCompartment="compartment";
  F6P.sbmlRole="species";
  F6P.sbmlCompartment="compartment";
  pyr.sbmlRole="species";
  pyr.sbmlCompartment="compartment";
  n3PG.sbmlRole="species";
  n3PG.sbmlCompartment="compartment";
  PPhi.sbmlRole="species";
  PPhi.sbmlCompartment="compartment";
  Phi.sbmlRole="species";
  Phi.sbmlCompartment="compartment";
  glycogen.sbmlRole="species";
  glycogen.sbmlCompartment="compartment";
  n3POHpyr.sbmlRole="species";
  n3POHpyr.sbmlCompartment="compartment";
  vADH.sbmlRole="rate";
  vAK.sbmlRole="rate";
  vALDO.sbmlRole="rate";
  vATPases.sbmlRole="rate";
  vDHases.sbmlRole="rate";
  vENO.sbmlRole="rate";
  vGAPDH.sbmlRole="rate";
  vHKrev.sbmlRole="rate";
  vHPhi.sbmlRole="rate";
  vPFORALDH.sbmlRole="rate";
  vPGAM.sbmlRole="rate";
  vPGK.sbmlRole="rate";
  vPPDK.sbmlRole="rate";
  vPPhiPFK.sbmlRole="rate";
  vPPhisynthesis.sbmlRole="rate";
  vTPhi.sbmlRole="rate";
  vglycogensynthesis.sbmlRole="rate";
  vglycogendegradation.sbmlRole="rate";
  vn3PGDH.sbmlRole="rate";
}