begin name
tyson7
end name

begin reactions
CKIdegradation {1.0}CKIt = {1.0}$EXTERNAL
CKItphosphorilationviaSK {1.0}CKIt = {1.0}$EXTERNAL
CKItsynthesis {1.0}$EXTERNAL = {1.0}CKIt
Cdc20activation {1.0}$EXTERNAL = {1.0}Cdc20a
Cdc20adegradation {1.0}Cdc20a = {1.0}$EXTERNAL
Cdc20ainhibition {1.0}Cdc20a = {1.0}$EXTERNAL
Cdc20tdeg {1.0}Cdc20t = {1.0}$EXTERNAL
Cdc20tsynthesis {1.0}$EXTERNAL = {1.0}Cdc20t
Cdh1degradation {1.0}Cdh1 = {1.0}$EXTERNAL
Cdh1synthesis {1.0}$EXTERNAL = {1.0}Cdh1
CycBdegradation {1.0}CycBt = {1.0}$EXTERNAL
CycBdegradationviaCdh1 {1.0}CycBt = {1.0}$EXTERNAL
CycBtsynthesis {1.0}$EXTERNAL = {1.0}CycBt
CycBtdegradationviaCdc20a {1.0}CycBt = {1.0}$EXTERNAL
IEPdegradation {1.0}IEP = {1.0}$EXTERNAL
IEPsynthesis {1.0}$EXTERNAL = {1.0}IEP
SKdegradation {1.0}SK = {1.0}$EXTERNAL
SKsynthesis {1.0}$EXTERNAL = {1.0}SK
eq7 {1.0}CKIt = {1.0}$EXTERNAL
growth {1.0}$EXTERNAL = {1.0}m
end reactions

begin rate equations
CKIdegradation = k12p*CKIt[t]
CKItphosphorilationviaSK = k12pp*CKIt[t]*SK[t]
CKItsynthesis = k11
Cdc20activation = (k7*(-Cdc20a[t] + Cdc20t[t])*IEP[t])/(J7 - Cdc20a[t] + Cdc20t[t])
Cdc20adegradation = k6*Cdc20a[t]
Cdc20ainhibition = (k8*Mad*Cdc20a[t])/(J8 + Cdc20a[t])
Cdc20tdeg = k6*Cdc20t[t]
Cdc20tsynthesis = k5p + (k5pp*((CycB*m[t])/J5)^n)/(1 + ((CycB*m[t])/J5)^n)
Cdh1degradation = (CycB*k4*Cdh1[t]*m[t] + k4p*Cdh1[t]*SK[t])/(J4 + Cdh1[t])
Cdh1synthesis = ((k3p + k3pp*Cdc20a[t])*(1 - Cdh1[t]))/(1 + J3 - Cdh1[t])
CycBdegradation = k2p*CycBt[t]
CycBdegradationviaCdh1 = k2pp*Cdh1[t]*CycBt[t]
CycBtsynthesis = k1
CycBtdegradationviaCdc20a = k2ppp*Cdc20a[t]*CycBt[t]
IEPdegradation = k10*IEP[t]
IEPsynthesis = CycB*k9*(1 - IEP[t])*m[t]
SKdegradation = k14*SK[t]
SKsynthesis = k13*TF
eq7 = CycB*k12ppp*CKIt[t]*m[t]
growth = mu*m[t]*(1 - m[t]/mmax)
end rate equations

begin 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
end parameters

begin initial conditions
CKIt[0] = CKIti
Cdc20a[0] = Cdc20ai
Cdc20t[0] = Cdc20ti
Cdh1[0] = Cdh1i
CycBt[0] = CycBti
IEP[0] = IEPi
SK[0] = SKi
m[0] = mi
end initial conditions

begin initial values
CKIti = 0.001
Cdc20ai = 0.001
Cdc20ti = 0.001
Cdh1i = 0.001
CycBti = 0.001
IEPi = 0.001
SKi = 0.001
mi = 0.5
end initial values

begin assignment rules
CycB := CycBt[t] - (2*CKIt[t]*CycBt[t])/(Keq^(-1) + CKIt[t] + CycBt[t] + Sqrt[-4*CKIt[t]*CycBt[t] + (Keq^(-1) + CKIt[t] + CycBt[t])^2])
Mad := 1
TF := GK[k15p*m[t] + k15pp*SK[t], k16p + CycB*k16pp*m[t], J15, J16]
Trimer := (2*CKIt[t]*CycBt[t])/(Keq^(-1) + CKIt[t] + CycBt[t] + Sqrt[-4*CKIt[t]*CycBt[t] + (Keq^(-1) + CKIt[t] + CycBt[t])^2])
end assignment rules

begin function definitions
GK[A1_,A2_,A3_,A4_] := (2*A1*A4)/(-A1 + A2 + A2*A3 + A1*A4 + Sqrt[-4*A1*(-A1 + A2)*A4 + (-A1 + A2 + A2*A3 + A1*A4)^2])
end function definitions

begin events
unnamed = (CycB < 0.1, {m[t] -> m[t]/2}, 0)
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units