#whillier3



#reactions
v_1:
 $pool = {1.0}cysi
  (cyse*Kdcyse)/volet + (cyse*Vmaxtranscyse)/((Kmtranscyse + cyse/volet)*volet)
v_10:
 {1.0}glni = $pool
  (Vmaxtransglni*glni)/((Kmtransglni + glni/voliwt)*voliwt)
v_11:
 $pool = {1.0}glni
  (glne*Vmaxtransglne)/((Kmtransglne + glne/volet)*volet)
v_12:
 {1.0}cysi = $pool
  (Kdcysi*cysi)/voliwt + (Vmaxtranscysi*cysi)/((Kmtranscysi + cysi/voliwt)*voliwt)
v_13:
 {1.0}glyi = $pool
  (Kdglyi*glyi)/voliwt + (Vmaxatransglyi*glyi)/((Kmatransglyi + glyi/voliwt)*voliwt) + (Vmaxgtransglyi*glyi)/((Kmgtransglyi + glyi/voliwt)*voliwt)
v_14:
 {1.0}gssgi = {1.0}gssge + {1.0}phosi
  (atpi*GSSGTRANS*k51*k52*k53*k54*k55*k56*gssgi)/((((k52*k53*k54*k55*k56 + k53*k54*k55*k56*k5min1)*gssgi)/voliwt + (adpi*(k54*k55*k56*k5min1*k5min2 + k55*k56*k5min1*k5min2*k5min3 + k56*k5min1*k5min2*k5min3*k5min4 + (k5min1*k5min2*k5min3*k5min4*k5min5*gssge)/volet))/voliwt + (atpi*(k51*k52*k54*k55*k56 + k51*k52*k55*k56*k5min3 + k51*k52*k56*k5min3*k5min4 + (k51*k52*k5min3*k5min4*k5min5*gssge)/volet + (adpi*(k51*k54*k55*k56*k5min2 + k51*k55*k56*k5min2*k5min3 + k51*k56*k5min2*k5min3*k5min4 + (k51*k5min2*k5min3*k5min4*k5min5*gssge)/volet))/voliwt + (gssgi*(k51*k52*k53*k54*k55 + k51*k52*k53*k54*k56 + k51*k52*k53*k55*k56 + k51*k53*k54*k55*k56 + k51*k52*k53*k56*k5min4 + ((k51*k52*k53*k54*k5min5 + k51*k52*k53*k5min4*k5min5)*gssge)/volet))/voliwt))/voliwt)*voliwt**2)
v_15:
 $pool = {1.0}glyi
  (glye*Kdglye)/volet + (glye*Vmaxatransglye)/((Kmatransglye + glye/volet)*volet) + (glye*Vmaxgtransglye)/((Kmgtransglye + glye/volet)*volet)
v_16:
 {1.0}alai = $pool
  (Kdalai*alai)/voliwt + (Vmaxtransalai*alai)/((Kmtransalai + alai/voliwt)*voliwt)
v_17:
 $pool = {1.0}alai
  (alae*Kdalae)/volet + (alae*Vmaxtransalae)/((Kmtransalae + alae/volet)*volet)
v_18:
 {1.0}aketoi = $pool
  (Vmaxtransaketoi*aketoi)/((Kmtransaketoi + aketoi/voliwt)*voliwt)
v_19:
 $pool = {1.0}aketoi
  (aketoe*Vmaxtransaketoe)/((Kmtransaketoe + aketoe/volet)*volet)
v_2:
 {1.0}glni = {1.0}ammi + {1.0}glui
  (Kdgln*glni)/voliwt
v_20:
 {1.0}naci = {1.0}nace
  (Kdnaci*naci)/voliwt
v_21:
 {1.0}nace = {1.0}naci
  (Kdnace*nace)/volet
v_3:
 {1.0}naci = {1.0}acetate + {1.0}cysi
  (k111*k112*k113*NACASE*naci)/((k112*k113 + k113*k11min1 + ((k112*k11min3 + k11min1*k11min3)*cysi)/voliwt + ((k111*k112 + k111*k113)*naci)/voliwt)*voliwt)
v_4:
 {1.0}aketoi + {1.0}alai = {1.0}glui + {1.0}pyrui
  (ALATRASE*((k13min1*k13min2*k13min3*k13min4*aketoi*alai)/voliwt**2 - (k131*k132*k133*k134*glui*pyrui)/voliwt**2))/(((k134*k13min1*k13min2 + k13min1*k13min2*k13min3)*aketoi)/voliwt + ((k132*k133*k134 + k133*k134*k13min1)*pyrui)/voliwt + (glui*(k131*k132*k134 + k131*k132*k13min3 + ((k131*k134*k13min2 + k131*k13min2*k13min3)*aketoi)/voliwt + ((k131*k132*k133 + k131*k133*k134)*pyrui)/voliwt))/voliwt + (alai*(k132*k13min3*k13min4 + k13min1*k13min3*k13min4 + ((k13min1*k13min2*k13min4 + k13min2*k13min3*k13min4)*aketoi)/voliwt + ((k132*k133*k13min4 + k133*k13min1*k13min4)*pyrui)/voliwt))/voliwt)
v_5:
 {1.0}cysi + {1.0}glui = {1.0}ggci + {1.0}phosi
  (atpi*GGCSASE*k21*k22*k23*k24*k25*k26*k2min7*cysi*glui)/((k24*k25*k26*k2min1*k2min2*k2min7 + k25*k26*k2min1*k2min2*k2min3*k2min7 + (atpi*k22*k23*k24*k25*k26*k2min7*cysi)/voliwt**2 + (atpi*k23*k24*k25*k26*k2min1*k2min7)/voliwt + (k25*k26*k27*gshi*(k2min1*k2min2*k2min3 + k24*(k2min1*k2min2 + (atpi*k23*(k2min1 + (k22*cysi)/voliwt))/voliwt)))/voliwt + (k2min6*k2min7*ggci*(k2min1*k2min2*k2min3 + k24*(k2min1*k2min2 + (atpi*k23*(k2min1 + (k22*cysi)/voliwt))/voliwt))*(k25 + (adpi*k2min5)/voliwt))/voliwt + (k21*k2min7*glui*(k24*k25*k26*k2min2 + k25*k26*k2min3*(k2min2 + (k22*cysi)/voliwt) + (adpi*atpi*k22*k23*k24*k2min5*cysi)/voliwt**3 + (atpi*k22*k23*k24*k25*cysi)/voliwt**2 + (atpi*k22*k23*k24*k26*cysi)/voliwt**2 + (atpi*k22*k23*k25*k26*cysi)/voliwt**2 + (atpi*k23*k24*k25*k26)/voliwt + (k22*k24*k25*k26*cysi)/voliwt))/voliwt)*voliwt**3)
v_6:
 {1.0}ggci + {1.0}glyi = {1.0}gshi + {1.0}phosi
  (atpi*GSHSASE*k31*k32*k33*k34*k35*k36*ggci*glyi)/((k34*k35*k36*k3min1*k3min2 + k35*k36*k3min1*k3min2*k3min3 + (atpi*k32*k33*k34*k35*k36*ggci)/voliwt**2 + (atpi*k33*k34*k35*k36*k3min1)/voliwt + (gshi*(k34*k35*k3min1*k3min2*k3min6 + k35*k3min1*k3min2*k3min3*k3min6 + (atpi*ggci*(k32*k33*k34*k35*k3min6 + (adpi*k32*k33*k34*k3min5*k3min6)/voliwt))/voliwt**2 + (adpi*(k34*k3min1*k3min2*k3min5*k3min6 + k3min1*k3min2*k3min3*k3min5*k3min6))/voliwt + (atpi*(k33*k34*k35*k3min1*k3min6 + (adpi*k33*k34*k3min1*k3min5*k3min6)/voliwt))/voliwt))/voliwt + (glyi*(k31*k34*k35*k36*k3min2 + k31*k35*k36*k3min2*k3min3 + (atpi*k31*k33*k34*k35*k36)/voliwt + (ggci*(k31*k32*k34*k35*k36 + k31*k32*k35*k36*k3min3 + (atpi*(k31*k32*k33*k34*k35 + k31*k32*k33*k34*k36 + k31*k32*k33*k35*k36 + (adpi*k31*k32*k33*k34*k3min5)/voliwt))/voliwt))/voliwt))/voliwt)*voliwt**3)
v_7:
 {1.0}ggci = {1.0}cysi + {1.0}oxopi
  (GGCTASE*k41*k42*k43*ggci)/((k42*k43 + k43*k4min1 + ((k42*k4min3 + k4min1*k4min3)*cysi)/voliwt + ((k41*k42 + k41*k43)*ggci)/voliwt)*voliwt)
v_8:
 {1.0}gshi = {0.5}gssgi
  (k61*gshi**2)/voliwt**2
v_9:
 {1.0}gssgi = {2.0}gshi
  (GSSGRTASE*(-((N2*nadpi*gshi**2)/voliwt**3) + (N1*nadphi*gssgi)/voliwt**2))/(D1 + (D17*nadphi*gshi**2*gssgi)/voliwt**4 + (D18*nadpi*gshi**2*gssgi)/voliwt**4 + (D14*nadphi*gshi**2)/voliwt**3 + (D16*nadpi*gshi**2)/voliwt**3 + ((D12 + D13)*nadphi*gshi*gssgi)/voliwt**3 + (D15*nadpi*gshi*gssgi)/voliwt**3 + (D7*nadphi*gshi)/voliwt**2 + ((D10 + D11)*nadpi*gshi)/voliwt**2 + (D9*gshi**2)/voliwt**2 + (D6*nadphi*gssgi)/voliwt**2 + (D8*nadpi*gssgi)/voliwt**2 + (D2*nadphi)/voliwt + (D5*nadpi)/voliwt + (D4*gshi)/voliwt + (D3*gssgi)/voliwt)


#parameters
EXTERNAL = 1.0
adpi = 1.263942e-07
aketoe = 2.1375e-08
alae = 5.13e-07
atpi = 1.786005e-06
cyse = 1.71e-08
glne = 8.6013e-07
glye = 4.5657e-07
htsus0 = 0.43
k111 = 11390000.0
k112 = 7000.0
k113 = 2000000.0
k11min1 = 3500.0
k11min3 = 10000000.0
k131 = 17410.0
k132 = 17600.0
k133 = 3827.0
k134 = 84.34
k13min1 = 54.95
k13min2 = 68300000.0
k13min3 = 100000.0
k13min4 = 18190.0
k21 = 66250.0
k22 = 443800.0
k23 = 368000.0
k24 = 100000.0
k25 = 161.9
k26 = 1000.0
k27 = 100000000.0
k2min1 = 100.0
k2min2 = 10.0
k2min3 = 10000.0
k2min5 = 100000.0
k2min6 = 10000000.0
k2min7 = 230000.0
k31 = 365000.0
k32 = 615000.0
k33 = 3395000.0
k34 = 1000.0
k35 = 843.4
k36 = 200.0
k3min1 = 10.0
k3min2 = 316.0
k3min3 = 10000.0
k3min5 = 100000.0
k3min6 = 1000000.0
k41 = 7200000.0
k42 = 3846.0
k43 = 10000.0
k4min1 = 15000.0
k4min3 = 5000000.0
k51 = 1000000.0
k52 = 100000.0
k53 = 4300000.0
k54 = 1000.0
k55 = 10000.0
k56 = 100.0
k5min1 = 100000.0
k5min2 = 100000.0
k5min3 = 100.0
k5min4 = 3000.0
k5min5 = 1000000.0
k71 = 85000000.0
k72 = 123000000.0
k73 = 810.0
k74 = 1000000.0
k75 = 1000000.0
k76 = 7000000.0
k7min1 = 510.0
k7min2 = 7200.0
k7min3 = 1000.0
k7min4 = 50000000.0
k7min5 = 50000000.0
k7min6 = 100000000.0
nadphi = 9.159e-08
nadpi = 4.5795e-10
vole0 = 1.0
voli0 = 1.0
voliw0 = 1.0
voltotalw0 = 1.0
default_compartment = 1.0


#initial values
acetate = 0.0
aketoi = 4.5795e-09
alai = 2.692746e-07
ammi = 1.263942e-07
cysi = 4.5795e-09
ggci = 4.5795e-09
glni = 5.367174e-07
glui = 3.810144e-07
glyi = 4.194822e-07
gshi = 2.93088e-06
gssge = 1.7955e-08
gssgi = 1.263942e-08
nace = 0.0
naci = 0.0
oxopi = 1.263942e-10
phosi = 1.263942e-07
pyrui = 1.8318e-08


#assignment rules
!F NACASE = 7.2e-12*voliwt
!F Kdglye = 0.000026800000000000004*voli0
!F Kmtranscysi = 1/50000
!F Kmtransaketoi = 0.00034
!F Vmaxtransglne = 9.83e-8*voli0
!F hbi = 0.02123818625889349*volhbn
!F Kdalai = 0.00011000000000000002*voli0
!F Kmatransglyi = 511/500000
!F D17 = k71*k72*(k73 + k7min3)*k7min4*k7min5
!F D9 = k7min1*k7min2*k7min3*k7min4*k7min5
!F D7 = k71*k76*k7min2*k7min3*k7min4
!F D11 = k7min1*(k73*k74 + k74*k7min2 + k7min2*k7min3)*k7min5*k7min6
!F D3 = k72*k73*k74*k75*k76
!F Kdcyse = 0.000175*voli0
!F D10 = k7min1*k7min2*k7min3*k7min4*k7min6
!F kali = 0.15*voliw0
!F Kmtranscyse = 1/50000
!F D18 = k72*(k73 + k7min3)*k7min4*k7min5*k7min6
!F Vmaxtransaketoi = 7.93e-8*voli0
!F Kmtransalai = 0.00034
!F D4 = k76*k7min1*k7min2*k7min3*k7min4
!F nae = 0.14*vole0
!F D1 = k75*k76*k7min1*(k73*k74 + k74*k7min2 + k7min2*k7min3)
!F N2 = k7min1*k7min2*k7min3*k7min4*k7min5*k7min6
!F D2 = k71*k75*k76*(k73*k74 + k74*k7min2 + k7min2*k7min3)
!F D5 = k75*k7min1*(k73*k74 + k74*k7min2 + k7min2*k7min3)*k7min6
!F cle = kale + nae
!F D14 = k71*k7min2*k7min3*k7min4*k7min5
!F D13 = k71*k72*k73*k74*k7min5
!F Kmatransglye = 511/500000
!F Vmaxtranscysi = 3.4e-8*voli0
!F k61 = 0.015415*voliwt
!F Vmaxtransalai = 7.93e-8*voli0
!F Kmtransalae = 0.00034
!F GGCSASE = voliwt/25000000
!F Kdcysi = 0.000175*voli0
!F volet = -voliwt + voltotalw0
!F GSSGRTASE = 6.5e-8*voliwt
!F Vmaxtransalae = 7.93e-8*voli0
!F voltotal = 0.003
!F Kdalae = 0.00011000000000000002*voli0
!F nai = 0.005*voliw0
!F Kdnaci = 0.000028*voli0
!F volhbn = 0.29*voltotal*htsus0
!F Kmgtransglyi = 67/1000000
!F Kmtransglni = 1/4000
!F GSHSASE = voliwt/20000000
!F Vmaxtransaketoe = 7.93e-8*voli0
!F ALATRASE = 7.000000000000001e-9*voliwt
!F Vmaxtransglni = 9.83e-8*voli0
!F Kdglyi = 0.000026800000000000004*voli0
!F Vmaxgtransglyi = 9.17e-9*voli0
!F Kmtransglne = 1/4000
!F D15 = k72*k73*k74*k7min5*k7min6
!F Vmaxatransglye = 3.87e-8*voli0
!F GGCTASE = 3.32e-10*voliwt
!F Kmtransaketoe = 0.00034
!F Vmaxtranscyse = 3.62e-8*voli0
!F Kdnace = 0.000028*voli0
!F D6 = k71*k72*(k73*k74*k75 + k73*k74*k76 + k73*k75*k76 + k74*k75*k76 + k75*k76*k7min3)
!F Vmaxatransglyi = 1.8700000000000002e-8*voli0
!F voliwt = ((adpi + atpi + cli + hbi + kali + nadphi + nadpi + nai + acetate + ammi + cysi + ggci + glni + glui + glyi + gshi + gssgi + naci + oxopi + phosi)*voltotalw0)/(adpi + atpi + cle + cli + cyse + glne + glye + hbi + kale + kali + nadphi + nadpi + nae + nai + acetate + ammi + cysi + ggci + glni + glui + glyi + gshi + gssge + gssgi + nace + naci + oxopi + phosi)
!F Kmgtransglye = 67/1000000
!F kale = 0.01*vole0
!F D8 = k72*k73*k74*k75*k7min6
!F Kdgln = 0.000014633333333333333*voliwt
!F N1 = k71*k72*k73*k74*k75*k76
!F Vmaxgtransglye = 1.5800000000000003e-8*voli0
!F cli = 0.095*voliw0
!F D12 = k71*k72*k76*(k73 + k7min3)*k7min4
!F D16 = (k73*k7min1 + k7min1*k7min2 + k7min1*k7min3 + k7min2*k7min3)*k7min4*k7min5*k7min6
!F GSSGTRANS = 1.0e-8*voliwt