#tyson5

FIX: ext

#reactions
v_1:
 {1.0}ext = {1.0}CycBT
  k1 + k0*SK
v_10:
 {1.0}ext = {1.0}IEP
  k9*m*CycB*(1 - IEP)
v_11:
 {1.0}IEP = {1.0}ext
  k10*IEP
v_12:
 {1.0}ext = {1.0}m
  0.005*m*(1 - m/10)
v_13:
 {1.0}ext = {1.0}SK
  1*TF + k0*SK
v_14:
 {1.0}SK = {1.0}ext
  1*SK
v_2:
 {1.0}CycBT = {1.0}ext
  (k21 + k22*Cdh1 + 1*Cdc20A)*CycBT
v_3:
 {1.0}ext = {1.0}Cdh1
  (k31 + k32*Cdc20A)*(1 - Cdh1)/(J3 + 1 - Cdh1)
v_4:
 {1.0}Cdh1 = {1.0}ext
  (2*SK + k4*m*CycB)*Cdh1/(J4 + Cdh1)
v_5:
 $pool = {1.0}Cdc20T
  k51 + k52*(CycB*m/J5)**n/(1 + (CycB*m/J5)**n)
v_6:
 {1.0}Cdc20T = {1.0}ext
  k6*(Cdc20T - Cdc20A)
v_7:
 {1.0}ext = {1.0}Cdc20A
  k7*IEP*(Cdc20T - Cdc20A)/(J7 + Cdc20T - Cdc20A)
v_8:
 {1.0}Cdc20A = {1.0}ext
  k8*MAD*Cdc20A/(J8 + Cdc20A)
v_9:
 {1.0}Cdc20A + {1.0}Cdc20T = {1.0}ext
  k6*Cdc20A


#parameters
J3 = 0.04
J4 = 0.04
J5 = 0.3
J7 = 0.001
J8 = 0.001
MAD = 1.0
k0 = 0.0
k1 = 0.04
k10 = 0.02
k21 = 0.04
k22 = 1.0
k31 = 1.0
k32 = 10.0
k4 = 35.0
k51 = 0.005
k52 = 0.2
k6 = 0.1
k7 = 1.0
k8 = 0.5
k9 = 0.1
mu = 0.01
n = 4.0
ext = 0.0
default_compartment = 1.0


#initial values
Cdc20A = 0.001
Cdc20T = 0.001
Cdh1 = 0.001
CycBT = 0.001
IEP = 0.001
SK = 0.001
m = 0.5


#assignment rules
!F CycB = CycBT
!F TF = GK((1.5*m + 0.05*SK), (1 + 3*m*CycB), 0.01, 0.01)