#tyson12

FIX: S

#reactions
r0:
 $pool = {1.0}R
  env * k0
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 * E
r3:
 $pool = $pool
  env * k3 * Ep / (Km3 + Ep)
r4:
 $pool = $pool
  env * k4 * R * E / (Km4 + E)


#parameters
Et = 1.0
J3 = 0.05
J4 = 0.05
k0 = 0.0
k1 = 0.05
k2 = 0.1
k2_prime = 0.5
k3 = 0.2
k4 = 1.0
S = 0.0
env = 1.0


#initial values
R = 0.0


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