// 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
// vak: asp atp => aspp adp 
// vasd: aspp nadph => asa nadp pi 
// vhdh: asa nadph => hs nadp 
// vtsy: hsp => thr pi 
// vhk: hs atp => hsp adp 

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

  // variable definitions
  real compartment = 1 L;
  real lys = .46;
  real pH = 7.5;
  real F1 = 1;
  real prot = 1;
  real vm11 = .68987;
  real k11 = .97;
  real k1thr = .167;
  real alpha = 2.47;
  real nak1 = 4.09;
  real k1aspp = .017;
  real k1atp = .98;
  real k1adp = .25;
  real keqak = 6.39E-4;
  real vm13 = .33488;
  real k13 = .323;
  real k1lys = .391;
  real nak3 = 2.78;
  real k13atp = .225;
  real k13aspp = .017;
  real k13adp = .25;
  real vm2f = .920322;
  real k2eq = 56.41503346;
  real k2aspp = .022;
  real k2nadph = .0287;
  real k2asa = .11;
  real k2nadp = .144;
  real k2p = 10;
  real vm3f = 3.6707;
  real k3eq = 3.1622777E3;
  real k3asa = .245;
  real k3asp = 10;
  real k3nadph = .037;
  real k3hs = 3.39;
  real k3nadp = .067;
  real k3thr = .094;
  real nhdh1 = 1.41;
  real alpha3 = 3.93;
  real vm4f = .54579;
  real k4atp = .072;
  real k4hs = .11;
  real k4thr = 1.09;
  real k4lys = 9.45;
  real k4iatp = 4.35;
  real k4ihs = 4.7;
  real vm5 = .2392;
  real k5hsp = .307;
  real asp = 1.34 M;
  real atp = 1.31 M;
  real aspp(time) M;
  real adp = .17 M;
  real nadph = .56 M;
  real asa(time) M;
  real nadp = .63 M;
  real pi = 5 M;
  real hs(time) M;
  real hsp(time) M;
  real thr = 3.49 M;
  real vak(time) katal;
  real vasd(time) katal;
  real vhdh(time) katal;
  real vtsy(time) katal;
  real vhk(time) katal;

  // equations
  when (time=time.min) aspp = 0;
  aspp:time = (vak + -1*vasd)/compartment;
  when (time=time.min) asa = 0;
  asa:time = (vasd + -1*vhdh)/compartment;
  when (time=time.min) hs = 0;
  hs:time = (vhdh + -1*vhk)/compartment;
  when (time=time.min) hsp = 0;
  hsp:time = (-1*vtsy + vhk)/compartment;
  vak = F1*(vm11*prot*(asp*atp-aspp*adp/keqak)/((k11*(1+(thr/k1thr)^nak1)/(1+(thr/(alpha*k1thr))^nak1)+k11*aspp/k1aspp+asp)*(k1atp*(1+adp/k1adp)+atp))+vm13*prot*(asp*atp-aspp*adp/keqak)/((1+(lys/k1lys)^nak3)*(k13*(1+aspp/k13aspp)+asp)*(k13atp*(1+adp/k13adp)+atp)));
  vasd = vm2f*prot*(aspp*nadph-asa*nadp*3.141592653589793/k2eq)/((k2aspp*(1+asa/k2asa)*(1+3.141592653589793/k2p)+aspp)*(k2nadph*(1+nadp/k2nadp)+nadph));
  vhdh = vm3f*prot*(asa*nadph-hs*nadp/k3eq)/((1+asp/k3asp)*(1+(thr/k3thr)^nhdh1)/(1+(thr/(alpha3*k3thr))^nhdh1)*(k3asa+asa+hs*k3asa/k3hs)*(k3nadph*(1+nadp/k3nadp)+nadph));
  vtsy = vm5*prot*hsp/(hsp+k5hsp);
  vhk = vm4f*prot*hs*atp/((1+lys/k4lys)*(atp+k4atp*(1+hs/k4ihs))*(hs+k4hs*(1+thr/k4thr)*(1+atp/k4iatp)));

  // variable properties
  compartment.sbmlRole="compartment";
  lys.sbmlRole="parameter";
  pH.sbmlRole="parameter";
  F1.sbmlRole="parameter";
  prot.sbmlRole="parameter";
  vm11.sbmlRole="parameter";
  k11.sbmlRole="parameter";
  k1thr.sbmlRole="parameter";
  alpha.sbmlRole="parameter";
  nak1.sbmlRole="parameter";
  k1aspp.sbmlRole="parameter";
  k1atp.sbmlRole="parameter";
  k1adp.sbmlRole="parameter";
  keqak.sbmlRole="parameter";
  vm13.sbmlRole="parameter";
  k13.sbmlRole="parameter";
  k1lys.sbmlRole="parameter";
  nak3.sbmlRole="parameter";
  k13atp.sbmlRole="parameter";
  k13aspp.sbmlRole="parameter";
  k13adp.sbmlRole="parameter";
  vm2f.sbmlRole="parameter";
  k2eq.sbmlRole="parameter";
  k2aspp.sbmlRole="parameter";
  k2nadph.sbmlRole="parameter";
  k2asa.sbmlRole="parameter";
  k2nadp.sbmlRole="parameter";
  k2p.sbmlRole="parameter";
  vm3f.sbmlRole="parameter";
  k3eq.sbmlRole="parameter";
  k3asa.sbmlRole="parameter";
  k3asp.sbmlRole="parameter";
  k3nadph.sbmlRole="parameter";
  k3hs.sbmlRole="parameter";
  k3nadp.sbmlRole="parameter";
  k3thr.sbmlRole="parameter";
  nhdh1.sbmlRole="parameter";
  alpha3.sbmlRole="parameter";
  vm4f.sbmlRole="parameter";
  k4atp.sbmlRole="parameter";
  k4hs.sbmlRole="parameter";
  k4thr.sbmlRole="parameter";
  k4lys.sbmlRole="parameter";
  k4iatp.sbmlRole="parameter";
  k4ihs.sbmlRole="parameter";
  vm5.sbmlRole="parameter";
  k5hsp.sbmlRole="parameter";
  asp.sbmlRole="species";
  asp.sbmlCompartment="compartment";
  atp.sbmlRole="species";
  atp.sbmlCompartment="compartment";
  aspp.sbmlRole="species";
  aspp.sbmlCompartment="compartment";
  adp.sbmlRole="species";
  adp.sbmlCompartment="compartment";
  nadph.sbmlRole="species";
  nadph.sbmlCompartment="compartment";
  asa.sbmlRole="species";
  asa.sbmlCompartment="compartment";
  nadp.sbmlRole="species";
  nadp.sbmlCompartment="compartment";
  pi.sbmlRole="species";
  pi.sbmlCompartment="compartment";
  hs.sbmlRole="species";
  hs.sbmlCompartment="compartment";
  hsp.sbmlRole="species";
  hsp.sbmlCompartment="compartment";
  thr.sbmlRole="species";
  thr.sbmlCompartment="compartment";
  vak.sbmlRole="rate";
  vasd.sbmlRole="rate";
  vhdh.sbmlRole="rate";
  vtsy.sbmlRole="rate";
  vhk.sbmlRole="rate";
}