#saavedra

FIX: Glu PPhi Phi etoh glycogen n3POHpyr

#reactions
v_1:
 {1.0}Glu + {1.0}ATP = {1.0}ADP + {1.0}G6P
  HKrevparVmf/(HKrevparalfa*HKrevparKa*HKrevparKb)*(Glu*ATP-G6P*ADP/HKrevparKeq)/(1+Glu/HKrevparKa+ATP/HKrevparKb+Glu*ATP/(HKrevparalfa*HKrevparKa*HKrevparKb)+AMP/HKrevparKi+Glu*AMP/(HKrevparbeta*HKrevparKa*HKrevparKi)+AMP*ATP/(HKrevpargama*HKrevparKi*HKrevparKb)+G6P/HKrevparKp+ADP/HKrevparKq+G6P*ADP/(HKrevparalfa*HKrevparKp*HKrevparKq))
v_10:
 {1.0}NAD + {1.0}n3PG = {1.0}NADH + {1.0}n3POHpyr
  n3PGDHparVm*NAD*n3PG/(n3PGDHparKa*n3PGDHparKb+n3PGDHparKb*NAD+NAD*n3PG)
v_11:
 {1.0}n3PG = {1.0}n2PG
  (PGAMparVf*n3PG/PGAMparKms-PGAMparVr*n2PG/PGAMparKmp)/(1+n3PG/PGAMparKms+n2PG/PGAMparKmp)
v_12:
 {1.0}n2PG = {1.0}PEP
  (ENOparVf*n2PG/ENOparKms-ENOparVr*PEP/ENOparKmp)/(1+n2PG/ENOparKms+PEP/ENOparKmp)
v_13:
 {1.0}PEP + {1.0}AMP + {1.0}PPhi = {1.0}pyr + {1.0}ATP + {1.0}Phi
  (PPDKparVmf*(PEP*AMP*PPhi/(PPDKparalfa*PPDKparKa*PPDKparKb*PPDKparKc))-PPDKparVmr*(pyr*ATP*Phi/(PPDKparbeta*PPDKparKp*PPDKparKq*PPDKparKr)))/(1+PEP/PPDKparKa+PEP*AMP/(PPDKparalfa*PPDKparKa*PPDKparKb)+PEP*AMP*PPhi/(PPDKparalfa*PPDKparbeta*PPDKparKa*PPDKparKb*PPDKparKc)+pyr*ATP*Phi/(PPDKparalfa*PPDKparbeta*PPDKparKp*PPDKparKq*PPDKparKr)+ATP*Phi/(PPDKparalfa*PPDKparKq*PPDKparKr)+Phi/PPDKparKr)
v_14:
 {1.0}NADH + {1.0}pyr = {1.0}NAD + {1.0}acald
  (PFORALDHparVmf*(NADH*pyr/(PFORALDHparalfa*PFORALDHparKa*PFORALDHparKb))-PFORALDHparVmr*(acald*NAD/(PFORALDHparbeta*PFORALDHparKp*PFORALDHparKq)))/(1+NADH/PFORALDHparKa+pyr/PFORALDHparKb+NADH*pyr/(PFORALDHparalfa*PFORALDHparKa*PFORALDHparKb)+acald*NAD/(PFORALDHparbeta*PFORALDHparKp*PFORALDHparKq)+acald/PFORALDHparKp+NAD/PFORALDHparKq)
v_15:
 {1.0}NADH + {1.0}acald = {1.0}NAD + {1.0}etoh
  ADHparVm*NADH*acald/(ADHparKa*ADHparKb+ADHparKb*NADH+NADH*acald)
v_16:
 {2.0}ADP = {1.0}AMP + {1.0}ATP
  (AKpark1*ADP**2 - AKpark2*ATP*AMP)
v_17:
 {1.0}ATP = {1.0}ADP
  ATPasesparK*ATP
v_18:
 {1.0}ATP = {1.0}AMP + {1.0}PPhi
  PPhisynthesispark*ATP
v_19:
 {1.0}NAD = {1.0}NADH
  (DHasespark1*NAD-DHasespark2*NADH)
v_2:
 {1.0}glycogen + {1.0}Phi = {1.0}G6P
  glycogendegradationpark*glycogen*Phi
v_3:
 {1.0}ATP + {1.0}G6P = {1.0}ADP + {1.0}PPhi + {1.0}glycogen
  glycogensynthesispark*G6P*ATP
v_4:
 {1.0}G6P = {1.0}F6P
  (HPhiparVf*G6P/HPhiparKms-HPhiparVr*F6P/HPhiparKmp)/(1+G6P/HPhiparKms+F6P/HPhiparKmp)
v_5:
 {1.0}PPhi + {1.0}F6P = {1.0}FBP + {1.0}Phi
  (PPhiPFKparVmf*(F6P*PPhi/(PPhiPFKparalfa*PPhiPFKparKa*PPhiPFKparKb))-PPhiPFKparVmr*(FBP*Phi/(PPhiPFKparbeta*PPhiPFKparKp*PPhiPFKparKq)))/(1+F6P/PPhiPFKparKa+PPhi/PPhiPFKparKb+F6P*PPhi/(PPhiPFKparalfa*PPhiPFKparKa*PPhiPFKparKb)+FBP*Phi/(PPhiPFKparbeta*PPhiPFKparKp*PPhiPFKparKq)+FBP/PPhiPFKparKp+Phi/PPhiPFKparKq)
v_6:
 {1.0}FBP = {1.0}G3P + {1.0}DHAP
  (ALDOparVmf*(FBP/ALDOparKs)-ALDOparVmr*(DHAP*G3P/(ALDOparKp+ALDOparKq)))/(1+FBP/ALDOparKs+DHAP/ALDOparKp+G3P/ALDOparKq+DHAP*G3P/(ALDOparKp*ALDOparKq))
v_7:
 {1.0}DHAP = {1.0}G3P
  (TPhiparVf*DHAP/TPhiparKms-TPhiparVr*G3P/TPhiparKmp)/(1+DHAP/TPhiparKms+G3P/TPhiparKmp)
v_8:
 {1.0}G3P + {1.0}NAD = {1.0}NADH + {1.0}n13BPG
  (GAPDHparVmf*(NAD*G3P/(GAPDHparalfa*GAPDHparKa*GAPDHparKb))-GAPDHparVmr*(n13BPG*NADH/(GAPDHparbeta*GAPDHparKp*GAPDHparKq)))/(1+NAD/GAPDHparKa+G3P/GAPDHparKb+NAD*G3P/(GAPDHparalfa*GAPDHparKa*GAPDHparKb)+n13BPG*NADH/(GAPDHparbeta*GAPDHparKp*GAPDHparKq)+n13BPG/GAPDHparKp+NADH/GAPDHparKq)
v_9:
 {1.0}ADP + {1.0}n13BPG = {1.0}n3PG + {1.0}ATP
  (PGKparVmf*(n13BPG*ADP/(PGKparalfa*PGKparKa*PGKparKb))-PGKparVmr*(n3PG*ATP/(PGKparbeta*PGKparKp*PGKparKq)))/(1+n13BPG/PGKparKa+ADP/PGKparKb+n13BPG*ADP/(PGKparalfa*PGKparKa*PGKparKb)+n3PG*ATP/(PGKparbeta*PGKparKp*PGKparKq)+n3PG/PGKparKp+ATP/PGKparKq)


#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
default_compartment = 1.0


#initial values
ADP = 3.3
AMP = 1.6
ATP = 5.0
DHAP = 1.15
F6P = 1.1
FBP = 0.43
G3P = 0.36
G6P = 6.2
NAD = 1.5
NADH = 0.05
PEP = 1.0
acald = 0.001
n13BPG = 1.0
n2PG = 1.0
n3PG = 1.0
pyr = 0.92


#assignment rules