#whillier4



#reactions
v_1:
 {1e-15}glne = {1.0}glni
  (glne*Vmaxtransglne)/((Kmtransglne + glne/volet)*volet)
v_10:
 {1.0}glni = {1.0}ammi + {1.0}glui
  (Kdgln*glni)/voliwt
v_11:
 {1e-15}alae = {1.0}alai
  (Kdalae*alae)/volet + (alae*Vmaxtransalae)/((Kmtransalae + alae/volet)*volet)
v_12:
 {1.0}alai = {1e-15}alae
  (Kdalai*alai)/voliwt + (alai*Vmaxtransalai)/((Kmtransalai + alai/voliwt)*voliwt)
v_13:
 {1.0}aketoe = {1.0}aketoi
  (aketoe*Vmaxtransaketoe)/((Kmtransaketoe + aketoe/volet)*volet)
v_14:
 {1.0}aketoi = {1.0}aketoe
  (aketoi*Vmaxtransaketoi)/((Kmtransaketoi + aketoi/voliwt)*voliwt)
v_15:
 {1.0}alai + {1.0}aketoi = {1.0}glui + {1.0}pyrui
  (ALATRASE*((k13r1*k13r2*k13r3*k13r4*aketoi*alai)/voliwt**2 - (k13f1*k13f2*k13f3*k13f4*glui*pyrui)/voliwt**2))/(((k13f4*k13r1*k13r2 + k13r1*k13r2*k13r3)*aketoi)/voliwt + ((k13f2*k13f3*k13f4 + k13f3*k13f4*k13r1)*pyrui)/voliwt + (glui*(k13f1*k13f2*k13f4 + k13f1*k13f2*k13r3 + ((k13f1*k13f4*k13r2 + k13f1*k13r2*k13r3)*aketoi)/voliwt + ((k13f1*k13f2*k13f3 + k13f1*k13f3*k13f4)*pyrui)/voliwt))/voliwt + (alai*(k13f2*k13r3*k13r4 + k13r1*k13r3*k13r4 + ((k13r1*k13r2*k13r4 + k13r2*k13r3*k13r4)*aketoi)/voliwt + ((k13f2*k13f3*k13r4 + k13f3*k13r1*k13r4)*pyrui)/voliwt))/voliwt)
v_16:
 {1.0}glui + {1.0}cysi + {1e-15}atpi = {1e-15}adpi + {1.0}phosi + {1.0}ggci
  (GGCSASE*k2f1*k2f2*k2f3*k2f4*k2f5*k2f6*k2r7*atpi*cysi*glui)/((k2f4*k2f5*k2f6*k2r1*k2r2*k2r7 + k2f5*k2f6*k2r1*k2r2*k2r3*k2r7 + (k2f2*k2f3*k2f4*k2f5*k2f6*k2r7*atpi*cysi)/voliwt**2 + (k2f3*k2f4*k2f5*k2f6*k2r1*k2r7*atpi)/voliwt + (k2f5*k2f6*k2f7*gshi*(k2r1*k2r2*k2r3 + k2f4*(k2r1*k2r2 + (k2f3*atpi*(k2r1 + (k2f2*cysi)/voliwt))/voliwt)))/voliwt + (k2r6*k2r7*ggci*(k2r1*k2r2*k2r3 + k2f4*(k2r1*k2r2 + (k2f3*atpi*(k2r1 + (k2f2*cysi)/voliwt))/voliwt))*(k2f5 + (k2r5*adpi)/voliwt))/voliwt + (k2f1*k2r7*glui*(k2f4*k2f5*k2f6*k2r2 + k2f5*k2f6*k2r3*(k2r2 + (k2f2*cysi)/voliwt) + (k2f2*k2f3*k2f4*k2r5*adpi*atpi*cysi)/voliwt**3 + (k2f2*k2f3*k2f4*k2f5*atpi*cysi)/voliwt**2 + (k2f2*k2f3*k2f4*k2f6*atpi*cysi)/voliwt**2 + (k2f2*k2f3*k2f5*k2f6*atpi*cysi)/voliwt**2 + (k2f3*k2f4*k2f5*k2f6*atpi)/voliwt + (k2f2*k2f4*k2f5*k2f6*cysi)/voliwt))/voliwt)*voliwt**3)
v_17:
 {1.0}ggci + {1.0}glyi + {1e-15}atpi = {1e-15}adpi + {1.0}phosi + {1.0}gshi
  (GSHSASE*k3f1*k3f2*k3f3*k3f4*k3f5*k3f6*atpi*ggci*glyi)/((k3f4*k3f5*k3f6*k3r1*k3r2 + k3f5*k3f6*k3r1*k3r2*k3r3 + (k3f2*k3f3*k3f4*k3f5*k3f6*atpi*ggci)/voliwt**2 + (k3f3*k3f4*k3f5*k3f6*k3r1*atpi)/voliwt + (gshi*(k3f4*k3f5*k3r1*k3r2*k3r6 + k3f5*k3r1*k3r2*k3r3*k3r6 + (atpi*ggci*(k3f2*k3f3*k3f4*k3f5*k3r6 + (k3f2*k3f3*k3f4*k3r5*k3r6*adpi)/voliwt))/voliwt**2 + ((k3f4*k3r1*k3r2*k3r5*k3r6 + k3r1*k3r2*k3r3*k3r5*k3r6)*adpi)/voliwt + (atpi*(k3f3*k3f4*k3f5*k3r1*k3r6 + (k3f3*k3f4*k3r1*k3r5*k3r6*adpi)/voliwt))/voliwt))/voliwt + (glyi*(k3f1*k3f4*k3f5*k3f6*k3r2 + k3f1*k3f5*k3f6*k3r2*k3r3 + (k3f1*k3f3*k3f4*k3f5*k3f6*atpi)/voliwt + (ggci*(k3f1*k3f2*k3f4*k3f5*k3f6 + k3f1*k3f2*k3f5*k3f6*k3r3 + (atpi*(k3f1*k3f2*k3f3*k3f4*k3f5 + k3f1*k3f2*k3f3*k3f4*k3f6 + k3f1*k3f2*k3f3*k3f5*k3f6 + (k3f1*k3f2*k3f3*k3f4*k3r5*adpi)/voliwt))/voliwt))/voliwt))/voliwt)*voliwt**3)
v_18:
 {1.0}ggci = {1.0}oxopi + {1.0}cysi
  (GGCTASE*k4f1*k4f2*k4f3*ggci)/((k4f2*k4f3 + k4f3*k4r1 + ((k4f2*k4r3 + k4r1*k4r3)*cysi)/voliwt + ((k4f1*k4f2 + k4f1*k4f3)*ggci)/voliwt)*voliwt)
v_19:
 {1.0}gssgi + {1e-15}atpi = {1e-15}adpi + {1.0}phosi + {1.0}gssge
  (GSSGTRANS*k5f1*k5f2*k5f3*k5f4*k5f5*k5f6*atpi*gssgi)/((((k5f2*k5f3*k5f4*k5f5*k5f6 + k5f3*k5f4*k5f5*k5f6*k5r1)*gssgi)/voliwt + (adpi*(k5f4*k5f5*k5f6*k5r1*k5r2 + k5f5*k5f6*k5r1*k5r2*k5r3 + k5f6*k5r1*k5r2*k5r3*k5r4 + (k5r1*k5r2*k5r3*k5r4*k5r5*gssge)/volet))/voliwt + (atpi*(k5f1*k5f2*k5f4*k5f5*k5f6 + k5f1*k5f2*k5f5*k5f6*k5r3 + k5f1*k5f2*k5f6*k5r3*k5r4 + (k5f1*k5f2*k5r3*k5r4*k5r5*gssge)/volet + (adpi*(k5f1*k5f4*k5f5*k5f6*k5r2 + k5f1*k5f5*k5f6*k5r2*k5r3 + k5f1*k5f6*k5r2*k5r3*k5r4 + (k5f1*k5r2*k5r3*k5r4*k5r5*gssge)/volet))/voliwt + (gssgi*(k5f1*k5f2*k5f3*k5f4*k5f5 + k5f1*k5f2*k5f3*k5f4*k5f6 + k5f1*k5f2*k5f3*k5f5*k5f6 + k5f1*k5f3*k5f4*k5f5*k5f6 + k5f1*k5f2*k5f3*k5f6*k5r4 + ((k5f1*k5f2*k5f3*k5f4*k5r5 + k5f1*k5f2*k5f3*k5r4*k5r5)*gssge)/volet))/voliwt))/voliwt)*voliwt**2)
v_2:
 {1.0}glni = {1e-15}glne
  (glni*Vmaxtransglni)/((Kmtransglni + glni/voliwt)*voliwt)
v_20:
 {1.0}gssgi + {1e-15}nadphi = {2.0}gshi + {1e-15}nadpi
  (GSSGRTASE*(-((N2*gshi**2*nadpi)/voliwt**3) + (N1*gssgi*nadphi)/voliwt**2))/(D1 + (D17*gshi**2*gssgi*nadphi)/voliwt**4 + (D18*gshi**2*gssgi*nadpi)/voliwt**4 + (D14*gshi**2*nadphi)/voliwt**3 + ((D12 + D13)*gshi*gssgi*nadphi)/voliwt**3 + (D16*gshi**2*nadpi)/voliwt**3 + (D15*gshi*gssgi*nadpi)/voliwt**3 + (D9*gshi**2)/voliwt**2 + (D7*gshi*nadphi)/voliwt**2 + (D6*gssgi*nadphi)/voliwt**2 + ((D10 + D11)*gshi*nadpi)/voliwt**2 + (D8*gssgi*nadpi)/voliwt**2 + (D4*gshi)/voliwt + (D3*gssgi)/voliwt + (D2*nadphi)/voliwt + (D5*nadpi)/voliwt)
v_21:
 {1.0}gshi = {0.5}gssgi
  k6f1*(gshi/voliwt)**2
v_3:
 {1e-15}cyse = {1.0}cysi
  (Kdcyse*cyse)/volet + (cyse*Vmaxtranscyse)/((Kmtranscyse + cyse/volet)*volet)
v_4:
 {1.0}cysi = {1e-15}cyse
  (Kdcysi*cysi)/voliwt + (cysi*Vmaxtranscysi)/((Kmtranscysi + cysi/voliw0)*voliwt)
v_5:
 {1.0}nace = {1.0}naci
  (Kdnace*nace)/volet
v_6:
 {1.0}naci = {1.0}nace
  (Kdnaci*naci)/voliwt
v_7:
 {1.0}naci = {1.0}cysi + {1.0}acetate
  (k11f1*k11f2*k11f3*NACASE*naci)/((k11f2*k11f3 + k11f3*k11r1 + ((k11f2*k11r3 + k11r1*k11r3)*cysi)/voliwt + ((k11f1*k11f2 + k11f1*k11f3)*naci)/voliwt)*voliwt)
v_8:
 {1e-15}glye = {1.0}glyi
  (Kdglye*glye)/volet + (glye*Vmaxatransglye)/((Kmatransglye + glye/volet)*volet) + (glye*Vmaxgtransglye)/((Kmgtransglye + glye/volet)*volet)
v_9:
 {1.0}glyi = {1e-15}glye
  (Kdglyi*glyi)/voliwt + (glyi*Vmaxatransglyi)/((Kmatransglyi + glyi/voliwt)*voliwt) + (glyi*Vmaxgtransglyi)/((Kmgtransglyi + glyi/voliwt)*voliwt)


#parameters
Kdalae = 1.419e-07
Kdalai = 1.419e-07
Kdcyse = 2.2575e-07
Kdcysi = 2.2575e-07
Kdglye = 3.4572e-08
Kdglyi = 3.4572e-08
Kdnace = 3.612e-08
Kdnaci = 3.612e-08
Kmatransglye = 0.001022
Kmatransglyi = 0.001022
Kmgtransglye = 6.7e-05
Kmgtransglyi = 6.7e-05
Kmtransaketoe = 0.00034
Kmtransaketoi = 0.00034
Kmtransalae = 0.00034
Kmtransalai = 0.00034
Kmtranscyse = 2e-05
Kmtranscysi = 2e-05
Kmtransglne = 0.00025
Kmtransglni = 0.00025
Vmaxatransglye = 4.9923e-11
Vmaxatransglyi = 2.4123e-11
Vmaxgtransglye = 2.0382e-11
Vmaxgtransglyi = 1.18293e-11
Vmaxtransaketoe = 1.02297e-10
Vmaxtransaketoi = 1.02297e-10
Vmaxtransalae = 1.02297e-10
Vmaxtransalai = 1.02297e-10
Vmaxtranscyse = 4.6698e-11
Vmaxtranscysi = 4.386e-11
Vmaxtransglne = 1.26807e-10
Vmaxtransglni = 1.26807e-10
cli = 8.70105e-05
htsus0 = 0.43
k11f1 = 11390000.0
k11f2 = 7000.0
k11f3 = 2000000.0
k11r1 = 3500.0
k11r3 = 10000000.0
k13f1 = 17410.0
k13f2 = 17600.0
k13f3 = 3827.0
k13f4 = 84.34
k13r1 = 54.95
k13r2 = 68300000.0
k13r3 = 100000.0
k13r4 = 18190.0
k2f1 = 66250.0
k2f2 = 443800.0
k2f3 = 368000.0
k2f4 = 100000.0
k2f5 = 161.9
k2f6 = 1000.0
k2f7 = 100000000.0
k2r1 = 100.0
k2r2 = 10.0
k2r3 = 10000.0
k2r5 = 100000.0
k2r6 = 10000000.0
k2r7 = 230000.0
k3f1 = 365000.0
k3f2 = 615000.0
k3f3 = 3395000.0
k3f4 = 1000.0
k3f5 = 843.4
k3f6 = 200.0
k3r1 = 10.0
k3r2 = 316.0
k3r3 = 10000.0
k3r5 = 100000.0
k3r6 = 1000000.0
k4f1 = 7200000.0
k4f2 = 3846.0
k4f3 = 10000.0
k4r1 = 15000.0
k4r3 = 5000000.0
k5f1 = 1000000.0
k5f2 = 100000.0
k5f3 = 4300000.0
k5f4 = 1000.0
k5f5 = 10000.0
k5f6 = 100.0
k5r1 = 100000.0
k5r2 = 100000.0
k5r3 = 100.0
k5r4 = 3000.0
k5r5 = 1000000.0
k7f1 = 85000000.0
k7f2 = 123000000.0
k7f3 = 810.0
k7f4 = 1000000.0
k7f5 = 1000000.0
k7f6 = 7000000.0
k7r1 = 510.0
k7r2 = 7200.0
k7r3 = 1000.0
k7r4 = 50000000.0
k7r5 = 50000000.0
k7r6 = 100000000.0
kale = 1.71e-05
kali = 0.000137385
nae = 0.0002394
nai = 4.5795e-06
vole0 = 0.00171
volhbn = 0.0003741
voli0 = 0.00129
voliw0 = 0.0009159
voltotal = 0.003
voltotalw0 = 0.0026259
default_compartment = 1.0


#initial values
acetate = 0.0
adpi = 1.26394e-07
aketoe = 2.1375e-08
aketoi = 4.5795e-09
alae = 5.13e-07
alai = 2.69275e-07
ammi = 1.26394e-07
atpi = 1.786e-06
cyse = 1.71e-08
cysi = 4.5795e-09
ggci = 4.5795e-09
glne = 1.00206e-06
glni = 4.60698e-07
glui = 3.81014e-07
glye = 4.5657e-07
glyi = 4.19482e-07
gshi = 2.93088e-06
gssge = 1.7955e-08
gssgi = 1.26394e-08
nace = 0.0
naci = 0.0
nadphi = 9.159e-08
nadpi = 4.5795e-10
oxopi = 1.26394e-10
phosi = 1.26394e-07
pyrui = 1.8318e-08


#assignment rules
!F D12 = k7f1*k7f2*k7f6*(k7f3 + k7r3)*k7r4
!F voliwt = ((cli + hbi + kali + nai + acetate + adpi + ammi + atpi + cysi + ggci + glni + glui + glyi + gshi + gssgi + naci + nadphi + nadpi + oxopi + phosi)*voltotalw0)/(cle + cli + hbi + kale + kali + nae + nai + acetate + adpi + ammi + atpi + cyse + cysi + ggci + glne + glni + glui + glye + glyi + gshi + gssge + gssgi + nace + naci + nadphi + nadpi + oxopi + phosi)
!F D9 = k7r1*k7r2*k7r3*k7r4*k7r5
!F volet = voltotalw0 - voliwt
!F D2 = k7f1*k7f5*k7f6*(k7f3*k7f4 + k7f4*k7r2 + k7r2*k7r3)
!F D6 = k7f1*k7f2*(k7f3*k7f4*k7f5 + k7f3*k7f4*k7f6 + k7f3*k7f5*k7f6 + k7f4*k7f5*k7f6 + k7f5*k7f6*k7r3)
!F D8 = k7f2*k7f3*k7f4*k7f5*k7r6
!F ALATRASE = 0.000000007*voliwt
!F GSSGRTASE = 0.000000065*voliwt
!F D18 = k7f2*(k7f3 + k7r3)*k7r4*k7r5*k7r6
!F cle = kale + nae
!F D16 = (k7f3*k7r1 + k7r1*k7r2 + k7r1*k7r3 + k7r2*k7r3)*k7r4*k7r5*k7r6
!F GSSGTRANS = 0.00000001*voliwt
!F N2 = k7r1*k7r2*k7r3*k7r4*k7r5*k7r6
!F D11 = k7r1*(k7f3*k7f4 + k7f4*k7r2 + k7r2*k7r3)*k7r5*k7r6
!F D5 = k7f5*k7r1*(k7f3*k7f4 + k7f4*k7r2 + k7r2*k7r3)*k7r6
!F D15 = k7f2*k7f3*k7f4*k7r5*k7r6
!F GGCTASE = 0.000000000332*voliwt
!F D3 = k7f2*k7f3*k7f4*k7f5*k7f6
!F D14 = k7f1*k7r2*k7r3*k7r4*k7r5
!F NACASE = 0.0000000000072*voliwt
!F D4 = k7f6*k7r1*k7r2*k7r3*k7r4
!F GGCSASE = 0.00000004*voliwt
!F D7 = k7f1*k7f6*k7r2*k7r3*k7r4
!F D17 = k7f1*k7f2*(k7f3 + k7r3)*k7r4*k7r5
!F Kdgln = 0.0000146333*voliwt
!F hbi = 0.0212382*volhbn
!F D13 = k7f1*k7f2*k7f3*k7f4*k7r5
!F GSHSASE = 0.00000005*voliwt
!F N1 = k7f1*k7f2*k7f3*k7f4*k7f5*k7f6
!F D10 = k7r1*k7r2*k7r3*k7r4*k7r6
!F D1 = k7f5*k7f6*k7r1*(k7f3*k7f4 + k7f4*k7r2 + k7r2*k7r3)
!F k6f1 = 0.015565*voliwt