#dupont2



#reactions
Protein_Phosphorylation:
 $pool = {1.0}W_star
  cytosol * (vp / Wt) * (vk / vp * ((1 - W_star) / (K1 + 1 - W_star)) - W_star / (K2 + W_star))
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(Kp, 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(K_A, p) + pow(Z, p))))
R_eff:
 {1.0}Z = $pool
  cytosol * k * Z
Rf:
 {1.0}Y = {1.0}Z
  store * kf * Y


#parameters
K1 = 0.01
K2 = 0.01
K_A = 0.9
Ka = 2.5
Kp = 1.0
Kr = 2.0
Vm2 = 65.0
Vm3 = 500.0
k = 10.0
kf = 1.0
m = 2.0
n = 2.0
p = 4.0
q = 1.0
v0 = 1.0
v1_beta = 2.7
vMK = 20.0
vp = 2.5
cytosol = 1.0
store = 1.0


#initial values
W_star = 0.0
Wt = 1.0
Y = 1.6
Z = 0.15


#assignment rules
!F vk = vMK * pow(Z, q) / (pow(Ka, q) + pow(Z, q))