#tyson7



#reactions
CKIdegradation:
 {1.0}CKIt = $pool
  k12p * CKIt
CKItphosphorilationviaSK:
 {1.0}CKIt = $pool
  k12pp * SK * CKIt
CKItsynthesis:
 $pool = {1.0}CKIt
  k11
Cdc20activation:
 $pool = {1.0}Cdc20a
  k7 * IEP * (Cdc20t - Cdc20a) / (J7 + Cdc20t - Cdc20a)
Cdc20adegradation:
 {1.0}Cdc20a = $pool
  k6 * Cdc20a
Cdc20ainhibition:
 {1.0}Cdc20a = $pool
  k8 * Mad * Cdc20a / (J8 + Cdc20a)
Cdc20t_deg:
 {1.0}Cdc20t = $pool
  k6 * Cdc20t
Cdc20tsynthesis:
 $pool = {1.0}Cdc20t
  k5p + k5pp * pow(CycB * m / J5, n) / (1 + pow(CycB * m / J5, n))
Cdh1degradation:
 {1.0}Cdh1 = $pool
  (k4p * SK * Cdh1 + k4 * m * CycB * Cdh1) / (J4 + Cdh1)
Cdh1synthesis:
 $pool = {1.0}Cdh1
  (k3p + k3pp * Cdc20a) * (1 - Cdh1) / (J3 + 1 - Cdh1)
CycBdegradation:
 {1.0}CycBt = $pool
  k2p * CycBt
CycBdegradationviaCdh1:
 {1.0}CycBt = $pool
  k2pp * Cdh1 * CycBt
CycBt_synthesis:
 $pool = {1.0}CycBt
  k1
CycBtdegradationviaCdc20a:
 {1.0}CycBt = $pool
  k2ppp * Cdc20a * CycBt
IEPdegradation:
 {1.0}IEP = $pool
  k10 * IEP
IEPsynthesis:
 $pool = {1.0}IEP
  k9 * m * CycB * (1 - IEP)
SKdegradation:
 {1.0}SK = $pool
  k14 * SK
SKsynthesis:
 $pool = {1.0}SK
  k13 * TF
eq_7:
 {1.0}CKIt = $pool
  k12ppp * m * CycB * CKIt
growth:
 $pool = {1.0}m
  mu * m * (1 - m / mmax)


#parameters
J15 = 0.01
J16 = 0.01
J3 = 0.04
J4 = 0.04
J5 = 0.3
J7 = 0.001
J8 = 0.001
Keq = 1000.0
k1 = 0.04
k10 = 0.02
k11 = 1.0
k12p = 0.2
k12pp = 50.0
k12ppp = 100.0
k13 = 1.0
k14 = 1.0
k15p = 1.5
k15pp = 0.05
k16p = 1.0
k16pp = 3.0
k2p = 0.04
k2pp = 1.0
k2ppp = 1.0
k3p = 1.0
k3pp = 10.0
k4 = 35.0
k4p = 2.0
k5p = 0.005
k5pp = 0.2
k6 = 0.1
k7 = 1.0
k8 = 0.5
k9 = 0.1
mmax = 10.0
mu = 0.005
n = 4.0
cell = 1.0


#initial values
CKIt = 0.001
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 - 2.0 * CycBt * CKIt / (CycBt + CKIt + 1.0 / Keq + pow(pow(CycBt + CKIt + 1.0 / Keq, 2.0) - 4.0 * CycBt * CKIt, 1.0 / 2.0))
!F Mad = 1.0
!F TF = GK(k15p * m + k15pp * SK, k16p + k16pp * m * CycB, J15, J16)
!F Trimer = 2.0 * CycBt * CKIt / (CycBt + CKIt + 1.0 / Keq + pow(pow(CycBt + CKIt + 1.0 / Keq, 2.0) - 4.0 * CycBt * CKIt, 1.0 / 2.0))