// This model generated automatically from SBML
// WARNING: Renamed variable PI to PI_ due to JSim namespace restrictions

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// vut1: PII UTP => PIIUMP PPi 
// vur1: PIIUMP => PII UMP 
// vut2: PIIUMP UTP => PIIUMP2 PPi 
// vur2: PIIUMP2 => PIIUMP UMP 
// vut3: PIIUMP2 UTP => PIIUMP3 PPi 
// vur3: PIIUMP3 => PIIUMP2 UMP 
// vad: GS => GSAMP 
// vdead: GSAMP => GS 
// vgdh: NH KG NADPH => GLU NADP 
// vgog: GLN NADPH KG => 2GLU NADP 
// vgs: GLU ATP NH => GLN ADP 
// vgludem: GLU => AZGLU 
// vazgludem: AZGLU => AZglu 
// vglndem: GLN => AZGLN 
// vazglndem: AZGLN => AZgln 
// vatpase: ADP => ATP 

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

  // variable definitions
  real compartment = 1 L;
  real UT = 6E-4;
  real kcatut = 137;
  real Kglnut = .07;
  real Kutipii = .0018;
  real Kutpii = .003;
  real Kutpiiump = .0035;
  real Kututp = .04;
  real Kutippi = .1135;
  real UR = 6E-4;
  real kcatur = 5.5;
  real Kurpiiump = .0023;
  real Kurump = 8.4;
  real Kglnur = .07;
  real a1 = 1E-22;
  real b1 = .5166;
  real c1 = .5974;
  real d1 = .0387;
  real Vad = .5;
  real Kadpiikg = 1.052E-5;
  real Kadgln = .9714;
  real Kadgs = .001703;
  real e1 = 1E-22;
  real f1 = 2.766;
  real g1 = 3.323;
  real h1 = .2148;
  real i1 = 1E-22;
  real j1 = 1E-22;
  real k1 = 1E-22;
  real l1 = .02316;
  real m1 = .8821;
  real n1 = 8.491;
  real o1 = .8791;
  real Vdead = .5;
  real Kdeadpiikg = 2.274E-6;
  real Kdeadgln = .04444;
  real Kdeadpiiu = 1.805E-5;
  real Kdeadgsa = 2.015E-4;
  real Vgdh = 360;
  real Kgdhkg = .32;
  real Kgdhnh = 1.1;
  real Kgdhglu = 10;
  real Kgdhnadph = .04;
  real Kgdhnadp = .042;
  real Keqgdh = 1290;
  real Kgdhazglu = 2.5;
  real Vgog = 85;
  real Kgoggln = .175;
  real Kgogkg = .007;
  real Kgognadph = .0015;
  real Kgogglu = 11;
  real Kgognadp = .0037;
  real Kgogaz = .65;
  real Vgs = 600;
  real aamp = 10;
  real bamp = 2.3667;
  real camp = .1012;
  real damp = 10.8688;
  real n1amp = 1.1456;
  real n2amp = 19.2166;
  real Kgseq = 460;
  real Kgsatp = .35;
  real Kgsglu = 4.1;
  real Kgsnh = .1;
  real Kgsadp = .0585;
  real Kgspi = 3.7;
  real Kgsgln = 5.65;
  real Keq = 460;
  real Vgludem = 120;
  real Kgludemglu = 8;
  real Kgludemeq = 1E10;
  real Kgludemazglu = .5;
  real Vglndem = 70;
  real Kglndemgln = 2;
  real Kglndemeq = 1E10;
  real Kglndemazgln = .25;
  real Vazglndem = 20;
  real Kazglndemazgln = 1;
  real Kazglndemeq = 1E10;
  real Kazglndemazinter = .5;
  real Vazgludem = 30;
  real Kazgludemazglu = .3;
  real Kazgludemeq = 1E10;
  real Kazgludemazinter = .5;
  real Vadp = 100;
  real Kadp = .5;
  real ATPtot = 5.37;
  real GStot = .014;
  real PIItot = .003;
  real Kd1 = .005;
  real Kd2 = .15;
  real Kd3 = .15;
  real Kd1piiump = .025;
  real Kd2piiump = .15;
  real Kd3piiump = .15;
  real azinter = .1;
  real PI_ = 10;
  real PII(time) M;
  real UTP = .5 M;
  real PIIUMP(time) M;
  real PPi = .05 M;
  real GLN(time) M;
  real PIIUMP2(time) M;
  real PIIUMP3(time) M;
  real UMP = .01 M;
  real GS(time) M;
  real GSAMP(time) M;
  real NH = .05 M;
  real KG = .2 M;
  real NADPH = .15 M;
  real GLU(time) M;
  real NADP = .05 M;
  real AZGLU(time) M;
  real ATP(time) M;
  real ADP(time) M;
  real AZglu = 1 M;
  real AZGLN(time) M;
  real AZgln = 1 M;
  real vut1(time) katal;
  real vur1(time) katal;
  real vut2(time) katal;
  real vur2(time) katal;
  real vut3(time) katal;
  real vur3(time) katal;
  real vad(time) katal;
  real vdead(time) katal;
  real vgdh(time) katal;
  real vgog(time) katal;
  real vgs(time) katal;
  real vgludem(time) katal;
  real vazgludem(time) katal;
  real vglndem(time) katal;
  real vazglndem(time) katal;
  real vatpase(time) katal;

  // equations
  when (time=time.min) PII = .003;
  PII:time = (-1*vut1 + vur1)/compartment;
  when (time=time.min) PIIUMP = 0;
  PIIUMP:time = (vut1 + -1*vur1 + -1*vut2 + vur2)/compartment;
  when (time=time.min) GLN = 1;
  GLN:time = (-1*vgog + vgs + -1*vglndem)/compartment;
  when (time=time.min) PIIUMP2 = 0;
  PIIUMP2:time = (vut2 + -1*vur2 + -1*vut3 + vur3)/compartment;
  when (time=time.min) PIIUMP3 = 0;
  PIIUMP3:time = (vut3 + -1*vur3)/compartment;
  when (time=time.min) GS = .014;
  GS:time = (-1*vad + vdead)/compartment;
  when (time=time.min) GSAMP = 0;
  GSAMP:time = (vad + -1*vdead)/compartment;
  when (time=time.min) GLU = 1;
  GLU:time = (vgdh + 2*vgog + -1*vgs + -1*vgludem)/compartment;
  when (time=time.min) AZGLU = 1;
  AZGLU:time = (vgludem + -1*vazgludem)/compartment;
  when (time=time.min) ATP = 2.685;
  ATP:time = (-1*vgs + vatpase)/compartment;
  when (time=time.min) ADP = 2.685;
  ADP:time = (vgs + -1*vatpase)/compartment;
  when (time=time.min) AZGLN = 1;
  AZGLN:time = (vglndem + -1*vazglndem)/compartment;
  vut1 = kcatut*UT*UTP*PII/(Kutipii*Kututp*(1+GLN/Kglnut)*(1+UTP/Kututp+(PII+PIIUMP+PIIUMP2)/Kutipii+UTP*(PII+PIIUMP+PIIUMP2)/(Kutipii*Kututp)+PPi*UTP*(PII+PIIUMP+PIIUMP2)/(Kutipii*Kutippi*Kututp)+Kutpii*(PIIUMP+PIIUMP2+PIIUMP3)/(Kutipii*Kutpiiump)));
  vur1 = kcatur*UR*PIIUMP/(Kurpiiump*(1+Kglnur/GLN)*(1+(1+UMP/Kurump)*(PIIUMP+PIIUMP2+PIIUMP3)/Kurpiiump));
  vut2 = kcatut*UT*UTP*PIIUMP/(Kutipii*Kututp*(1+GLN/Kglnut)*(1+UTP/Kututp+(PII+PIIUMP+PIIUMP2)/Kutipii+UTP*(PII+PIIUMP+PIIUMP2)/(Kutipii*Kututp)+PPi*UTP*(PII+PIIUMP+PIIUMP2)/(Kutipii*Kutippi*Kututp)+Kutpii*(PIIUMP+PIIUMP2+PIIUMP3)/(Kutipii*Kutpiiump)));
  vur2 = kcatur*UR*PIIUMP2/(Kurpiiump*(1+Kglnur/GLN)*(1+(1+UMP/Kurump)*(PIIUMP+PIIUMP2+PIIUMP3)/Kurpiiump));
  vut3 = kcatut*UT*UTP*PIIUMP2/(Kutipii*Kututp*(1+GLN/Kglnut)*(1+UTP/Kututp+(PII+PIIUMP+PIIUMP2)/Kutipii+UTP*(PII+PIIUMP+PIIUMP2)/(Kutipii*Kututp)+PPi*UTP*(PII+PIIUMP+PIIUMP2)/(Kutipii*Kutippi*Kututp)+Kutpii*(PIIUMP+PIIUMP2+PIIUMP3)/(Kutipii*Kutpiiump)));
  vur3 = kcatur*UR*PIIUMP3/(Kurpiiump*(1+Kglnur/GLN)*(1+(1+UMP/Kurump)*(PIIUMP+PIIUMP2+PIIUMP3)/Kurpiiump));
  vad = Vad*GS*(b1*GLN/Kadgln+3*a1*KG*PII/(Kadpiikg*Kd1*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3)))+3*c1*KG*GLN*PII/(Kadgln*Kadpiikg*Kd1*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3))))/((Kadgs+GS)*(1+GLN/Kadgln+3*KG*PII/(Kadpiikg*Kd1*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3)))+3*KG*GLN*PII/(d1*Kadgln*Kadpiikg*Kd1*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3)))));
  vdead = Vdead*GSAMP*(f1*GLN/Kdeadgln+3*e1*KG*PII/(Kd1*Kdeadpiikg*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3)))+3*h1*KG*GLN*PII/(Kd1*Kdeadgln*Kdeadpiikg*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3)))+g1*KG^3*PIIUMP3/(Kd1piiump*Kd2piiump*Kd3piiump*Kdeadpiiu*(1+3*KG/Kd1piiump+3*KG^2/(Kd1piiump*Kd2piiump)+KG^3/(Kd1piiump*Kd2piiump*Kd3piiump)))+j1*KG^3*GLN*PIIUMP3/(Kd1piiump*Kd2piiump*Kd3piiump*Kdeadgln*Kdeadpiiu*(1+3*KG/Kd1piiump+3*KG^2/(Kd1piiump*Kd2piiump)+KG^3/(Kd1piiump*Kd2piiump*Kd3piiump)))+3*i1*KG^4*PII*PIIUMP3/(Kd1*Kd1piiump*Kd2piiump*Kd3piiump*Kdeadpiikg*Kdeadpiiu*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3))*(1+3*KG/Kd1piiump+3*KG^2/(Kd1piiump*Kd2piiump)+KG^3/(Kd1piiump*Kd2piiump*Kd3piiump)))+3*k1*KG^4*GLN*PII*PIIUMP3/(Kd1*Kd1piiump*Kd2piiump*Kd3piiump*Kdeadgln*Kdeadpiikg*Kdeadpiiu*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3))*(1+3*KG/Kd1piiump+3*KG^2/(Kd1piiump*Kd2piiump)+KG^3/(Kd1piiump*Kd2piiump*Kd3piiump))))/((Kdeadgsa+GSAMP)*(1+GLN/Kdeadgln+3*KG*PII/(Kd1*Kdeadpiikg*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3)))+3*KG*GLN*PII/(Kd1*Kdeadgln*Kdeadpiikg*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3))*l1)+KG^3*PIIUMP3/(Kd1piiump*Kd2piiump*Kd3piiump*Kdeadpiiu*(1+3*KG/Kd1piiump+3*KG^2/(Kd1piiump*Kd2piiump)+KG^3/(Kd1piiump*Kd2piiump*Kd3piiump)))+KG^3*GLN*PIIUMP3/(Kd1piiump*Kd2piiump*Kd3piiump*Kdeadgln*Kdeadpiiu*(1+3*KG/Kd1piiump+3*KG^2/(Kd1piiump*Kd2piiump)+KG^3/(Kd1piiump*Kd2piiump*Kd3piiump))*n1)+3*KG^4*PII*PIIUMP3/(Kd1*Kd1piiump*Kd2piiump*Kd3piiump*Kdeadpiikg*Kdeadpiiu*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3))*(1+3*KG/Kd1piiump+3*KG^2/(Kd1piiump*Kd2piiump)+KG^3/(Kd1piiump*Kd2piiump*Kd3piiump))*m1)+3*KG^4*GLN*PII*PIIUMP3/(Kd1*Kd1piiump*Kd2piiump*Kd3piiump*Kdeadgln*Kdeadpiikg*Kdeadpiiu*(1+3*KG/Kd1+3*KG^2/(Kd1*Kd2)+KG^3/(Kd1*Kd2*Kd3))*(1+3*KG/Kd1piiump+3*KG^2/(Kd1piiump*Kd2piiump)+KG^3/(Kd1piiump*Kd2piiump*Kd3piiump))*o1)));
  vgdh = Vgdh*(KG*NADPH*NH-NADP*GLU/Keqgdh)/(Kgdhkg*Kgdhnadph*Kgdhnh*(1+NADP/Kgdhnadp+NADPH/Kgdhnadph)*(1+NH/Kgdhnh)*(1+KG/Kgdhkg+GLU/Kgdhglu));
  vgog = KG*NADPH*Vgog*GLN/(Kgoggln*Kgogkg*Kgognadph*(1+NADP/Kgognadp+NADPH/Kgognadph)*(1+AZGLU/Kgogaz)*(1+KG/Kgogkg+GLU/Kgogglu)*(1+GLN/Kgoggln+GLU/Kgogglu));
  vgs = aamp*camp*Vgs*((-1)*(3.141592653589793*ADP*GLN/Keq)+NH*ATP*GLU)/(Kgsatp*Kgsglu*Kgsnh*(1+3.141592653589793/Kgspi+ADP/Kgsadp+3.141592653589793*ADP/(Kgsadp*Kgspi)+ATP/Kgsatp)*(1+NH/Kgsnh+GLN/Kgsgln+NH*GLN/(Kgsgln*Kgsnh)+GLU/Kgsglu+NH*GLU/(Kgsglu*Kgsnh))*(1+12^n1amp*(GSAMP/(bamp*GStot))^n1amp)*(1+12^n2amp*(GSAMP/(damp*GStot))^n2amp));
  vgludem = Vgludem*((-1)*(AZGLU/Kgludemeq)+GLU)/(Kgludemglu*(1+AZGLU/Kgludemazglu+GLU/Kgludemglu));
  vazgludem = Vazgludem*((-1)*(azinter/Kazgludemeq)+AZGLU)/(Kazgludemazglu*(1+azinter/Kazgludemazinter+AZGLU/Kazgludemazglu));
  vglndem = Vglndem*((-1)*(AZGLN/Kglndemeq)+GLN)/(Kglndemgln*(1+AZGLN/Kglndemazgln+GLN/Kglndemgln));
  vazglndem = Vazglndem*((-1)*(azinter/Kazglndemeq)+AZGLN)/(Kazglndemazgln*(1+azinter/Kazglndemazinter+AZGLN/Kazglndemazgln));
  vatpase = Vadp*ADP/(Kadp+ADP);

  // variable properties
  compartment.sbmlRole="compartment";
  UT.sbmlRole="parameter";
  kcatut.sbmlRole="parameter";
  Kglnut.sbmlRole="parameter";
  Kutipii.sbmlRole="parameter";
  Kutpii.sbmlRole="parameter";
  Kutpiiump.sbmlRole="parameter";
  Kututp.sbmlRole="parameter";
  Kutippi.sbmlRole="parameter";
  UR.sbmlRole="parameter";
  kcatur.sbmlRole="parameter";
  Kurpiiump.sbmlRole="parameter";
  Kurump.sbmlRole="parameter";
  Kglnur.sbmlRole="parameter";
  a1.sbmlRole="parameter";
  b1.sbmlRole="parameter";
  c1.sbmlRole="parameter";
  d1.sbmlRole="parameter";
  Vad.sbmlRole="parameter";
  Kadpiikg.sbmlRole="parameter";
  Kadgln.sbmlRole="parameter";
  Kadgs.sbmlRole="parameter";
  e1.sbmlRole="parameter";
  f1.sbmlRole="parameter";
  g1.sbmlRole="parameter";
  h1.sbmlRole="parameter";
  i1.sbmlRole="parameter";
  j1.sbmlRole="parameter";
  k1.sbmlRole="parameter";
  l1.sbmlRole="parameter";
  m1.sbmlRole="parameter";
  n1.sbmlRole="parameter";
  o1.sbmlRole="parameter";
  Vdead.sbmlRole="parameter";
  Kdeadpiikg.sbmlRole="parameter";
  Kdeadgln.sbmlRole="parameter";
  Kdeadpiiu.sbmlRole="parameter";
  Kdeadgsa.sbmlRole="parameter";
  Vgdh.sbmlRole="parameter";
  Kgdhkg.sbmlRole="parameter";
  Kgdhnh.sbmlRole="parameter";
  Kgdhglu.sbmlRole="parameter";
  Kgdhnadph.sbmlRole="parameter";
  Kgdhnadp.sbmlRole="parameter";
  Keqgdh.sbmlRole="parameter";
  Kgdhazglu.sbmlRole="parameter";
  Vgog.sbmlRole="parameter";
  Kgoggln.sbmlRole="parameter";
  Kgogkg.sbmlRole="parameter";
  Kgognadph.sbmlRole="parameter";
  Kgogglu.sbmlRole="parameter";
  Kgognadp.sbmlRole="parameter";
  Kgogaz.sbmlRole="parameter";
  Vgs.sbmlRole="parameter";
  aamp.sbmlRole="parameter";
  bamp.sbmlRole="parameter";
  camp.sbmlRole="parameter";
  damp.sbmlRole="parameter";
  n1amp.sbmlRole="parameter";
  n2amp.sbmlRole="parameter";
  Kgseq.sbmlRole="parameter";
  Kgsatp.sbmlRole="parameter";
  Kgsglu.sbmlRole="parameter";
  Kgsnh.sbmlRole="parameter";
  Kgsadp.sbmlRole="parameter";
  Kgspi.sbmlRole="parameter";
  Kgsgln.sbmlRole="parameter";
  Keq.sbmlRole="parameter";
  Vgludem.sbmlRole="parameter";
  Kgludemglu.sbmlRole="parameter";
  Kgludemeq.sbmlRole="parameter";
  Kgludemazglu.sbmlRole="parameter";
  Vglndem.sbmlRole="parameter";
  Kglndemgln.sbmlRole="parameter";
  Kglndemeq.sbmlRole="parameter";
  Kglndemazgln.sbmlRole="parameter";
  Vazglndem.sbmlRole="parameter";
  Kazglndemazgln.sbmlRole="parameter";
  Kazglndemeq.sbmlRole="parameter";
  Kazglndemazinter.sbmlRole="parameter";
  Vazgludem.sbmlRole="parameter";
  Kazgludemazglu.sbmlRole="parameter";
  Kazgludemeq.sbmlRole="parameter";
  Kazgludemazinter.sbmlRole="parameter";
  Vadp.sbmlRole="parameter";
  Kadp.sbmlRole="parameter";
  ATPtot.sbmlRole="parameter";
  GStot.sbmlRole="parameter";
  PIItot.sbmlRole="parameter";
  Kd1.sbmlRole="parameter";
  Kd2.sbmlRole="parameter";
  Kd3.sbmlRole="parameter";
  Kd1piiump.sbmlRole="parameter";
  Kd2piiump.sbmlRole="parameter";
  Kd3piiump.sbmlRole="parameter";
  azinter.sbmlRole="parameter";
  PI_.sbmlRole="parameter";
  PII.sbmlRole="species";
  PII.sbmlCompartment="compartment";
  UTP.sbmlRole="species";
  UTP.sbmlCompartment="compartment";
  PIIUMP.sbmlRole="species";
  PIIUMP.sbmlCompartment="compartment";
  PPi.sbmlRole="species";
  PPi.sbmlCompartment="compartment";
  GLN.sbmlRole="species";
  GLN.sbmlCompartment="compartment";
  PIIUMP2.sbmlRole="species";
  PIIUMP2.sbmlCompartment="compartment";
  PIIUMP3.sbmlRole="species";
  PIIUMP3.sbmlCompartment="compartment";
  UMP.sbmlRole="species";
  UMP.sbmlCompartment="compartment";
  GS.sbmlRole="species";
  GS.sbmlCompartment="compartment";
  GSAMP.sbmlRole="species";
  GSAMP.sbmlCompartment="compartment";
  NH.sbmlRole="species";
  NH.sbmlCompartment="compartment";
  KG.sbmlRole="species";
  KG.sbmlCompartment="compartment";
  NADPH.sbmlRole="species";
  NADPH.sbmlCompartment="compartment";
  GLU.sbmlRole="species";
  GLU.sbmlCompartment="compartment";
  NADP.sbmlRole="species";
  NADP.sbmlCompartment="compartment";
  AZGLU.sbmlRole="species";
  AZGLU.sbmlCompartment="compartment";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  ADP.sbmlRole="species";
  ADP.sbmlCompartment="compartment";
  AZglu.sbmlRole="species";
  AZglu.sbmlCompartment="compartment";
  AZGLN.sbmlRole="species";
  AZGLN.sbmlCompartment="compartment";
  AZgln.sbmlRole="species";
  AZgln.sbmlCompartment="compartment";
  vut1.sbmlRole="rate";
  vur1.sbmlRole="rate";
  vut2.sbmlRole="rate";
  vur2.sbmlRole="rate";
  vut3.sbmlRole="rate";
  vur3.sbmlRole="rate";
  vad.sbmlRole="rate";
  vdead.sbmlRole="rate";
  vgdh.sbmlRole="rate";
  vgog.sbmlRole="rate";
  vgs.sbmlRole="rate";
  vgludem.sbmlRole="rate";
  vazgludem.sbmlRole="rate";
  vglndem.sbmlRole="rate";
  vazglndem.sbmlRole="rate";
  vatpase.sbmlRole="rate";
}