#fung1

FIX: HOAc_E

#reactions
R_Acs:
 $pool = {1.0}Acs
  compartment * (alpha2 * pow(AcP / Kg2, n) / (1 + pow(AcP / Kg2, n)) + alpha0)
R_LacI:
 $pool = {1.0}LacI
  compartment * (alpha1 * pow(AcP / Kg1, n) / (1 + pow(AcP / Kg1, n)) + alpha0)
R_Pta:
 $pool = {1.0}Pta
  alpha3 / (1 + pow(LacI / Kg3, n)) + alpha0
R_dAcs:
 {1.0}Acs = $pool
  compartment * kd * Acs
R_dLacI:
 {1.0}LacI = $pool
  compartment * kd * LacI
R_dPta:
 {1.0}Pta = $pool
  compartment * kd * Pta
V_Ace:
 {1.0}OAc = {1.0}HOAc
  compartment * C * (OAc * H - Keq * HOAc)
V_Ack:
 {1.0}AcP = {1.0}OAc
  compartment * (kAck_f * AcP - kAck_r * OAc)
V_Acs:
 {1.0}OAc = {1.0}AcCoA
  compartment * k2 * Acs * OAc / (KM2 + OAc)
V_Pta:
 {1.0}AcCoA = {1.0}AcP
  compartment * k1 * Pta * AcCoA / (KM1 + AcCoA)
V_TCA:
 {1.0}AcCoA = $pool
  compartment * kTCA * AcCoA
V_gly:
 $pool = {1.0}AcCoA
  compartment * S0
V_out:
 {1.0}HOAc = {1.0}HOAc_E
  compartment * k3 * (HOAc - HOAc_E)


#parameters
C = 100.0
H = 1e-07
KM1 = 0.06
KM2 = 0.1
Keq = 0.0005
Kg1 = 10.0
Kg2 = 10.0
Kg3 = 0.001
S0 = 0.5
alpha0 = 0.0
alpha1 = 0.1
alpha2 = 2.0
alpha3 = 2.0
k1 = 80.0
k2 = 0.8
k3 = 0.01
kAck_f = 1.0
kAck_r = 1.0
kTCA = 10.0
kd = 0.06
n = 2.0
HOAc_E = 0.0
compartment = 1.0


#initial values
AcCoA = 0.0
AcP = 0.0
Acs = 0.0
HOAc = 0.0
LacI = 0.0
OAc = 0.0
Pta = 0.0


#assignment rules