begin name
saavedra
end name

begin reactions
v1 {1.0}$Glu + {1.0}ATP = {1.0}ADP + {1.0}G6P
v10 {1.0}NAD + {1.0}n3PG = {1.0}NADH + {1.0}$n3POHpyr
v11 {1.0}n3PG = {1.0}n2PG
v12 {1.0}n2PG = {1.0}PEP
v13 {1.0}PEP + {1.0}AMP + {1.0}$PPhi = {1.0}pyr + {1.0}ATP + {1.0}$Phi
v14 {1.0}NADH + {1.0}pyr = {1.0}NAD + {1.0}acald
v15 {1.0}NADH + {1.0}acald = {1.0}NAD + {1.0}$etoh
v16 {2.0}ADP = {1.0}AMP + {1.0}ATP
v17 {1.0}ATP = {1.0}ADP
v18 {1.0}ATP = {1.0}AMP + {1.0}$PPhi
v19 {1.0}NAD = {1.0}NADH
v2 {1.0}$glycogen + {1.0}$Phi = {1.0}G6P
v3 {1.0}ATP + {1.0}G6P = {1.0}ADP + {1.0}$PPhi + {1.0}$glycogen
v4 {1.0}G6P = {1.0}F6P
v5 {1.0}$PPhi + {1.0}F6P = {1.0}FBP + {1.0}$Phi
v6 {1.0}FBP = {1.0}G3P + {1.0}DHAP
v7 {1.0}DHAP = {1.0}G3P
v8 {1.0}G3P + {1.0}NAD = {1.0}NADH + {1.0}n13BPG
v9 {1.0}ADP + {1.0}n13BPG = {1.0}n3PG + {1.0}ATP
end reactions

begin rate equations
v1 = (HKrevparVmf*(Glu*ATP[t] - (ADP[t]*G6P[t])/HKrevparKeq))/(HKrevparalfa*HKrevparKa*HKrevparKb*(1 + Glu/HKrevparKa + ADP[t]/HKrevparKq + AMP[t]/HKrevparKi + (Glu*AMP[t])/(HKrevparbeta*HKrevparKa*HKrevparKi) + ATP[t]/HKrevparKb + (Glu*ATP[t])/(HKrevparalfa*HKrevparKa*HKrevparKb) + (AMP[t]*ATP[t])/(HKrevpargama*HKrevparKb*HKrevparKi) + G6P[t]/HKrevparKp + (ADP[t]*G6P[t])/(HKrevparalfa*HKrevparKp*HKrevparKq)))
v10 = (n3PGDHparVm*n3PG[t]*NAD[t])/(n3PGDHparKa*n3PGDHparKb + n3PGDHparKb*NAD[t] + n3PG[t]*NAD[t])
v11 = (-((PGAMparVr*n2PG[t])/PGAMparKmp) + (PGAMparVf*n3PG[t])/PGAMparKms)/(1 + n2PG[t]/PGAMparKmp + n3PG[t]/PGAMparKms)
v12 = ((ENOparVf*n2PG[t])/ENOparKms - (ENOparVr*PEP[t])/ENOparKmp)/(1 + n2PG[t]/ENOparKms + PEP[t]/ENOparKmp)
v13 = ((PPDKparVmf*PPhi*AMP[t]*PEP[t])/(PPDKparalfa*PPDKparKa*PPDKparKb*PPDKparKc) - (Phi*PPDKparVmr*ATP[t]*pyr[t])/(PPDKparbeta*PPDKparKp*PPDKparKq*PPDKparKr))/(1 + Phi/PPDKparKr + (Phi*ATP[t])/(PPDKparalfa*PPDKparKq*PPDKparKr) + PEP[t]/PPDKparKa + (AMP[t]*PEP[t])/(PPDKparalfa*PPDKparKa*PPDKparKb) + (PPhi*AMP[t]*PEP[t])/(PPDKparalfa*PPDKparbeta*PPDKparKa*PPDKparKb*PPDKparKc) + (Phi*ATP[t]*pyr[t])/(PPDKparalfa*PPDKparbeta*PPDKparKp*PPDKparKq*PPDKparKr))
v14 = (-((PFORALDHparVmr*acald[t]*NAD[t])/(PFORALDHparbeta*PFORALDHparKp*PFORALDHparKq)) + (PFORALDHparVmf*NADH[t]*pyr[t])/(PFORALDHparalfa*PFORALDHparKa*PFORALDHparKb))/(1 + acald[t]/PFORALDHparKp + NAD[t]/PFORALDHparKq + (acald[t]*NAD[t])/(PFORALDHparbeta*PFORALDHparKp*PFORALDHparKq) + NADH[t]/PFORALDHparKa + pyr[t]/PFORALDHparKb + (NADH[t]*pyr[t])/(PFORALDHparalfa*PFORALDHparKa*PFORALDHparKb))
v15 = (ADHparVm*acald[t]*NADH[t])/(ADHparKa*ADHparKb + ADHparKb*NADH[t] + acald[t]*NADH[t])
v16 = AKpark1*ADP[t]^2 - AKpark2*AMP[t]*ATP[t]
v17 = ATPasesparK*ATP[t]
v18 = PPhisynthesispark*ATP[t]
v19 = DHasespark1*NAD[t] - DHasespark2*NADH[t]
v2 = glycogen*glycogendegradationpark*Phi
v3 = glycogensynthesispark*ATP[t]*G6P[t]
v4 = (-((HPhiparVr*F6P[t])/HPhiparKmp) + (HPhiparVf*G6P[t])/HPhiparKms)/(1 + F6P[t]/HPhiparKmp + G6P[t]/HPhiparKms)
v5 = ((PPhi*PPhiPFKparVmf*F6P[t])/(PPhiPFKparalfa*PPhiPFKparKa*PPhiPFKparKb) - (Phi*PPhiPFKparVmr*FBP[t])/(PPhiPFKparbeta*PPhiPFKparKp*PPhiPFKparKq))/(1 + PPhi/PPhiPFKparKb + Phi/PPhiPFKparKq + F6P[t]/PPhiPFKparKa + (PPhi*F6P[t])/(PPhiPFKparalfa*PPhiPFKparKa*PPhiPFKparKb) + FBP[t]/PPhiPFKparKp + (Phi*FBP[t])/(PPhiPFKparbeta*PPhiPFKparKp*PPhiPFKparKq))
v6 = ((ALDOparVmf*FBP[t])/ALDOparKs - (ALDOparVmr*DHAP[t]*G3P[t])/(ALDOparKp + ALDOparKq))/(1 + DHAP[t]/ALDOparKp + FBP[t]/ALDOparKs + G3P[t]/ALDOparKq + (DHAP[t]*G3P[t])/(ALDOparKp*ALDOparKq))
v7 = ((TPhiparVf*DHAP[t])/TPhiparKms - (TPhiparVr*G3P[t])/TPhiparKmp)/(1 + DHAP[t]/TPhiparKms + G3P[t]/TPhiparKmp)
v8 = ((GAPDHparVmf*G3P[t]*NAD[t])/(GAPDHparalfa*GAPDHparKa*GAPDHparKb) - (GAPDHparVmr*n13BPG[t]*NADH[t])/(GAPDHparbeta*GAPDHparKp*GAPDHparKq))/(1 + G3P[t]/GAPDHparKb + n13BPG[t]/GAPDHparKp + NAD[t]/GAPDHparKa + (G3P[t]*NAD[t])/(GAPDHparalfa*GAPDHparKa*GAPDHparKb) + NADH[t]/GAPDHparKq + (n13BPG[t]*NADH[t])/(GAPDHparbeta*GAPDHparKp*GAPDHparKq))
v9 = ((PGKparVmf*ADP[t]*n13BPG[t])/(PGKparalfa*PGKparKa*PGKparKb) - (PGKparVmr*ATP[t]*n3PG[t])/(PGKparbeta*PGKparKp*PGKparKq))/(1 + ADP[t]/PGKparKb + ATP[t]/PGKparKq + n13BPG[t]/PGKparKa + (ADP[t]*n13BPG[t])/(PGKparalfa*PGKparKa*PGKparKb) + n3PG[t]/PGKparKp + (ATP[t]*n3PG[t])/(PGKparbeta*PGKparKp*PGKparKq))
end rate equations

begin parameters
ADHparKa = 0.05
ADHparKb = 0.15
ADHparVm = 0.4
AKpark1 = 2.26
AKpark2 = 1.0
ALDOparKp = 0.21
ALDOparKq = 0.26
ALDOparKs = 0.028
ALDOparVmf = 0.16
ALDOparVmr = 0.284
ATPasesparK = 0.002
DHasespark1 = 0.1
DHasespark2 = 1.5
ENOparKmp = 0.102
ENOparKms = 0.06
ENOparVf = 0.508
ENOparVr = 0.103
GAPDHparKa = 0.083
GAPDHparKb = 0.043
GAPDHparKp = 0.016
GAPDHparKq = 0.06
GAPDHparVmf = 0.405
GAPDHparVmr = 3.68
GAPDHparalfa = 1.0
GAPDHparbeta = 1.0
HKrevparKa = 0.025
HKrevparKb = 0.121
HKrevparKeq = 656.0
HKrevparKi = 0.036
HKrevparKp = 1.0
HKrevparKq = 0.235
HKrevparVmf = 0.095
HKrevparalfa = 1.0
HKrevparbeta = 1.0
HKrevpargama = 1.0
HPhiparKmp = 0.46
HPhiparKms = 0.61
HPhiparVf = 0.233
HPhiparVr = 0.206
PFORALDHparKa = 0.18
PFORALDHparKb = 0.14
PFORALDHparKp = 0.1
PFORALDHparKq = 0.2
PFORALDHparVmf = 0.3
PFORALDHparVmr = 0.01
PFORALDHparalfa = 1.0
PFORALDHparbeta = 1.0
PGAMparKmp = 0.106
PGAMparKms = 0.83
PGAMparVf = 0.116
PGAMparVr = 0.104
PGKparKa = 0.125
PGKparKb = 0.04
PGKparKp = 0.505
PGKparKq = 0.061
PGKparVmf = 3.182
PGKparVmr = 1.742
PGKparalfa = 1.0
PGKparbeta = 1.0
PPDKparKa = 0.03
PPDKparKb = 0.002
PPDKparKc = 0.091
PPDKparKp = 0.3
PPDKparKq = 0.28
PPDKparKr = 1.0
PPDKparVmf = 0.304
PPDKparVmr = 0.019
PPDKparalfa = 0.01
PPDKparbeta = 2.0
PPhiPFKparKa = 0.695
PPhiPFKparKb = 0.38
PPhiPFKparKp = 0.109
PPhiPFKparKq = 2.3
PPhiPFKparVmf = 0.213
PPhiPFKparVmr = 0.346
PPhiPFKparalfa = 0.5
PPhiPFKparbeta = 2.0
PPhisynthesispark = 0.0065
TPhiparKmp = 0.32
TPhiparKms = 0.445
TPhiparVf = 4.366
TPhiparVr = 6.098
glycogendegradationpark = 0.0001
glycogensynthesispark = 0.0015
n3PGDHparKa = 0.087
n3PGDHparKb = 0.212
n3PGDHparVm = 0.01
Glu = 5.0
PPhi = 0.45
Phi = 5.0
etoh = 10.0
glycogen = 1.0
n3POHpyr = 1.0
defaultcompartment = 1.0
end parameters

begin initial conditions
ADP[0] = ADPi
AMP[0] = AMPi
ATP[0] = ATPi
DHAP[0] = DHAPi
F6P[0] = F6Pi
FBP[0] = FBPi
G3P[0] = G3Pi
G6P[0] = G6Pi
NAD[0] = NADi
NADH[0] = NADHi
PEP[0] = PEPi
acald[0] = acaldi
n13BPG[0] = n13BPGi
n2PG[0] = n2PGi
n3PG[0] = n3PGi
pyr[0] = pyri
end initial conditions

begin initial values
ADPi = 3.3
AMPi = 1.6
ATPi = 5.0
DHAPi = 1.15
F6Pi = 1.1
FBPi = 0.43
G3Pi = 0.36
G6Pi = 6.2
NADi = 1.5
NADHi = 0.05
PEPi = 1.0
acaldi = 0.001
n13BPGi = 1.0
n2PGi = 1.0
n3PGi = 1.0
pyri = 0.92
end initial values

begin assignment rules
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units