#tyson8

FIX: S

#reactions
r0:
 $pool = {1.0}R
  env * k0 * Ep
r1:
 $pool = {1.0}R
  env * k1 * S
r2:
 {1.0}R = $pool
  env * k2 * R
r2_prime:
 {1.0}R = $pool
  env * k2_prime * R * X
r3:
 $pool = $pool
  env * k3 * R * E / (Km3 + E)
r4:
 $pool = $pool
  env * k4 * Ep / (Km4 + Ep)
r5:
 $pool = {1.0}X
  env * k5 * R
r6:
 {1.0}X = $pool
  env * k6 * X


#parameters
Et = 1.0
J3 = 0.3
J4 = 0.3
k0 = 4.0
k1 = 1.0
k2 = 1.0
k2_prime = 1.0
k3 = 1.0
k4 = 1.0
k5 = 0.1
k6 = 0.075
S = 0.0
env = 1.0


#initial values
R = 0.0
X = 0.0


#assignment rules
!F Km3 = J3 * Et
!F Ep = goldbeter_koshland(k3 * R, k4, J3, J4) * Et
!F E = Et - Ep
!F Km4 = J4 * Et