#goldbeter4



#reactions
R0:
 $pool = {1.0}Z
  cytosol * v0
R1:
 $pool = {1.0}Z
  cytosol * v1 * beta
R2:
 {1.0}Z = {1.0}Y
  cytosol * (Vm2 * pow(Z, n) / (pow(K2, n) + pow(Z, n)))
R3:
 {1.0}Y = {1.0}Z
  store * (Vm3 * pow(Y, m) * pow(Z, p) / ((pow(Kr, m) + pow(Y, m)) * (pow(Ka, p) + pow(Z, p))))
R_eff:
 {1.0}Z = $pool
  cytosol * k * Z
Rf:
 {1.0}Y = {1.0}Z
  store * kf * Y


#parameters
K2 = 1.0
Ka = 0.9
Kr = 2.0
Vm2 = 65.0
Vm3 = 500.0
beta = 0.301
k = 10.0
kf = 1.0
m = 2.0
n = 2.0
p = 4.0
v0 = 1.0
v1 = 7.3
cytosol = 1.0
store = 1.0


#initial values
Y = 1.6
Z = 0.15


#assignment rules