begin name
kouril3
end name

begin reactions
v1 {1.0}atp + {1.0}p3g = {1.0}adp + {1.0}bpg
v2 {1.0}bpg + {1.0}nadph = {1.0}gap + {1.0}nadp
v3 {1.0}gap = {1.0}dhap
v4 {1.0}dhap + {1.0}gap = {1.0}f6p
v5 {1.0}gap = {1.0}$EXTERNAL
v6 {1.0}dhap = {1.0}$EXTERNAL
v7 {1.0}bpg = {1.0}$EXTERNAL
v8 {1.0}glc + {1.0}nadp = {1.0}nadph + {1.0}$ga
v9 {1.0}adp + {1.0}pep = {1.0}atp + {1.0}pyr
end reactions

begin rate equations
v1 = (protPGK*(-((VmrPGK*adp[t]*bpg[t])/(KpgkADP*KpgkbisP)) + (VmfPGK*atp[t]*p3g[t])/(KpgkATP*KpgkP3G)))/((1 + adp[t]/KiADP)*(1 + ((1 + adp[t]/KpgkADP)*bpg[t])/KpgkbisP + ((1 + atp[t]/KpgkATP)*p3g[t])/KpgkP3G))
v2 = (protGAPDH*(bpg[t]/KbisP + (pi*gap[t])/(Kgap*Kpi))^(-1 + n)*(-((pi*Vmarev*gap[t]*nadp[t])/(Kgap*Knadp*Kpi)) + (Vmfor*bpg[t]*nadph[t])/(KbisP*Knadph)))/((1 + (pi/Kpi)^n + (gap[t]/Kgap)^n + (bpg[t]/KbisP + (pi*gap[t])/(Kgap*Kpi))^n)*(1 + nadp[t]/Knadp + nadph[t]/Knadph))
v3 = (protTIM*(-((VmrTIM*dhap[t])/KtimDHAP) + (VmfTIM*gap[t])/KtimGAP))/(1 + dhap[t]/KtimDHAP + gap[t]/KtimGAP + p3g[t]/KtimP3G + pep[t]/KtimPEP)
v4 = (protALDPase*VmfALD*dhap[t]*gap[t])/(KaldDHAP*KaldGAP*(1 + dhap[t]/KaldDHAP + gap[t]/KaldGAP + (dhap[t]*gap[t])/(KaldDHAP*KaldGAP)))
v5 = kdgap*gap[t]
v6 = kddhap*dhap[t]
v7 = kdbpg*bpg[t]
v8 = kGDH*glc[t]*nadp[t]
v9 = kPK*adp[t]*pep[t]
end rate equations

begin 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
defaultcompartment = 1.0
end parameters

begin initial conditions
adp[0] = adpi
atp[0] = atpi
bpg[0] = bpgi
dhap[0] = dhapi
f6p[0] = f6pi
gap[0] = gapi
glc[0] = glci
nadp[0] = nadpi
nadph[0] = nadphi
p3g[0] = p3gi
pep[0] = pepi
pyr[0] = pyri
end initial conditions

begin initial values
adpi = 0.0
atpi = 10.0
bpgi = 0.0
dhapi = 0.03
f6pi = 0.12
gapi = 0.06
glci = 10.0
nadpi = 0.0
nadphi = 0.2
p3gi = 1.8
pepi = 5.0
pyri = 0.4
end initial values

begin assignment rules
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
v_1 urn:miriam:kegg.reaction:R01512
v_2 urn:miriam:kegg.reaction:R01061
v_3 urn:miriam:kegg.reaction:R01015
v_4 urn:miriam:kegg.reaction:R01068
end process annotations

begin species annotations
adp http://identifiers.org/kegg.compound/C00008
adp http://identifiers.org/obo.chebi/CHEBI%3A16761
atp urn:miriam:obo.chebi:CHEBI:15422
atp urn:miriam:kegg.compound:C00002
bpg http://identifiers.org/obo.chebi/CHEBI%3A28907
bpg http://identifiers.org/kegg.compound/C00236
dhap http://identifiers.org/obo.chebi/CHEBI%3A16108
dhap http://identifiers.org/kegg.compound/C00111
f6p http://identifiers.org/kegg.compound/C00085
f6p http://identifiers.org/obo.chebi/CHEBI%3A15946
gap urn:miriam:obo.chebi:CHEBI:17138
gap urn:miriam:kegg.compound:C00661
glc http://identifiers.org/obo.chebi/CHEBI%3A17234
glc http://identifiers.org/kegg.compound/C00031
nadp http://identifiers.org/kegg.compound/C00006
nadp http://identifiers.org/obo.chebi/CHEBI%3A18009
nadph http://identifiers.org/kegg.compound/C00005
nadph http://identifiers.org/obo.chebi/CHEBI%3A16474
p3g http://identifiers.org/kegg.compound/C00197
p3g http://identifiers.org/obo.chebi/CHEBI%3A17794
pep http://identifiers.org/kegg.compound/C00074
pep http://identifiers.org/obo.chebi/CHEBI%3A44897
pyr http://identifiers.org/kegg.compound/C00022
pyr http://identifiers.org/obo.chebi/CHEBI%3A32816
end species annotations

begin units
end units