#kouril3

FIX: ga

#reactions
v_1:
 {1.0}atp + {1.0}p3g = {1.0}adp + {1.0}bpg
  protPGK*(-VmrPGK*adp*bpg/(KpgkADP*KpgkbisP) + VmfPGK*(atp*p3g)/(KpgkATP*KpgkP3G))/((1 + adp/KiADP)*(1 + p3g/KpgkP3G*(1 + atp/KpgkATP) + bpg/KpgkbisP*(1 + adp/KpgkADP)))
v_2:
 {1.0}bpg + {1.0}nadph = {1.0}gap + {1.0}nadp
  protGAPDH * (-Vmarev * gap * nadp / (Kgap * Knadp) * (pi / Kpi) + Vmfor * bpg * nadph / (KbisP * Knadph)) * (gap / Kgap * pi / Kpi + bpg / KbisP)**(n - 1) / ((1 + nadp / Knadp + nadph / Knadph) * (1 + (gap / Kgap)**n + (pi / Kpi)**n + (gap / Kgap * (pi / Kpi) + bpg / KbisP)**n))
v_3:
 {1.0}gap = {1.0}dhap
  protTIM*(-VmrTIM*dhap/KtimDHAP + VmfTIM*gap/KtimGAP)/(1 + dhap/KtimDHAP + gap/KtimGAP + pep/KtimPEP + p3g/KtimP3G)
v_4:
 {1.0}dhap + {1.0}gap = {1.0}f6p
  protALDPase*(VmfALD*gap*dhap/(KaldGAP*KaldDHAP))/(1 + gap/KaldGAP + dhap/KaldDHAP + gap*dhap/(KaldGAP*KaldDHAP))
v_5:
 {1.0}gap = $pool
  kdgap * gap
v_6:
 {1.0}dhap = $pool
  kddhap * dhap
v_7:
 {1.0}bpg = $pool
  kdbpg * bpg
v_8:
 {1.0}glc + {1.0}nadp = {1.0}nadph + {1.0}ga
  kGDH*nadp*glc
v_9:
 {1.0}adp + {1.0}pep = {1.0}atp + {1.0}pyr
  kPK*adp*pep


#parameters
KaldDHAP = 0.170995
KaldF6P = 1.0
KaldGAP = 0.0522437
KbisP = 0.000406744
Kgap = 0.83766
KiADP = 1.14174
Knadp = 0.271013
Knadph = 0.0735253
KpgkADP = 0.0848049
KpgkATP = 9.68404
KpgkP3G = 0.541454
KpgkbisP = 5.59
Kpi = 408.523
KtimDHAP = 0.811626
KtimGAP = 0.245106
KtimP3G = 0.399822
KtimPEP = 0.660519
Vmarev = 24.054
VmfALD = 2.64096
VmfPGK = 17.3221
VmfTIM = 239.835
Vmfor = 20.6526
VmrPGK = 30.04
VmrTIM = 239.414
kGDH = 10.0
kPK = 10.0
kdbpg = 1.05824
kddhap = 0.0225
kdgap = 0.0559
n = 1.56416
pi = 0.0
protALDPase = 0.0036
protGAPDH = 0.0421
protPGK = 0.0034
protTIM = 0.00085
ga = 0.0
default_compartment = 1.0


#initial values
adp = 0.0
atp = 10.0
bpg = 0.0
dhap = 0.03
f6p = 0.12
gap = 0.06
glc = 10.0
nadp = 0.0
nadph = 0.2
p3g = 1.8
pep = 5.0
pyr = 0.4


#assignment rules