begin name
novak
end name

begin reactions
v1 {1.0}ERG = {1.0}$p
v10 {1.0}CA = {1.0}P27
v11 {1.0}CE = {1.0}P27 + {1.0}CYCE
v12 {1.0}CA = {1.0}P27 + {1.0}CYCA
v13 {1.0}CE = {1.0}P27
v14 {1.0}CYCE = {1.0}$p
v15 {1.0}P27 = {1.0}$p
v16 {1.0}CE = {1.0}CYCE
v17 {1.0}CD = {1.0}CYCD
v18 {1.0}CA = {1.0}CYCA
v19 {1.0}CYCB = {1.0}$p
v2 {1.0}DRG = {1.0}$p
v20 {1.0}$s = {1.0}CDh1
v21 {1.0}CDh1 = {1.0}$p
v22 {1.0}PPX = {1.0}$p
v23 {1.0}$s = {1.0}IEP
v24 {1.0}IEP = {1.0}$p
v25 {1.0}CDc20T = {1.0}$p
v26 {1.0}$s = {1.0}CDc20
v27 {1.0}CDc20 = {1.0}$p
v28 {1.0}CDc20 = {1.0}$p
v29 {1.0}$s = {1.0}E2F
v3 {1.0}CD = {1.0}P27
v30 {1.0}E2F = {1.0}$p
v31 {1.0}$s = {1.0}GM
v32 {1.0}GM = {1.0}$p
v33 {1.0}$s = {1.0}MASS
v34 {1.0}$s = {1.0}ERG
v35 {1.0}$s = {1.0}CDc20T
v36 {1.0}$s = {1.0}CYCA
v37 {1.0}$s = {1.0}PPX
v38 {1.0}$s = {1.0}CYCE
v39 {1.0}$s = {1.0}CYCD
v4 {1.0}CYCD = {1.0}$p
v40 {1.0}$s = {1.0}P27
v41 {1.0}$s = {1.0}DRG
v42 {1.0}$s = {1.0}CYCB
v5 {1.0}P27 + {1.0}CYCE = {1.0}CE
v6 {1.0}P27 + {1.0}CYCA = {1.0}CA
v7 {1.0}P27 + {1.0}CYCD = {1.0}CD
v8 {1.0}CD = {1.0}P27 + {1.0}CYCD
v9 {1.0}CYCA = {1.0}$p
end reactions

begin rate equations
v1 = k16*ERG[t]
v10 = K30*CA[t]*CDc20[t]
v11 = K25R*CE[t]
v12 = K25R*CA[t]
v13 = v8*CE[t]
v14 = v8*CYCE[t]
v15 = v6*P27[t]
v16 = v6*CE[t]
v17 = v6*CD[t]
v18 = v6*CA[t]
v19 = v2*CYCB[t]
v2 = k18*DRG[t]
v20 = ((K3a + K3*CDc20[t])*(1 - CDh1[t]))/(1 + J3 - CDh1[t])
v21 = (v4*CDh1[t])/(J4 + CDh1[t])
v22 = K34*PPX[t]
v23 = (K31*CYCB[t]*(1 - IEP[t]))/(1 + J31 - IEP[t])
v24 = (K32*IEP[t]*PPX[t])/(J32 + IEP[t])
v25 = K12*CDc20T[t]
v26 = (K13*(-CDc20[t] + CDc20T[t])*IEP[t])/(J13 - CDc20[t] + CDc20T[t])
v27 = (K14*CDc20[t])/(J14 + CDc20[t])
v28 = K12*CDc20[t]
v29 = K22*(E2FT - E2F[t])
v3 = K10*CD[t]
v30 = (K23a + K23*(CYCA[t] + CYCB[t]))*E2F[t]
v31 = K27*MASS[t]*Piecewise[{{0, RBH/RBT > 0.8}}, 1]
v32 = K28*GM[t]
v33 = eps*MU*GM[t]
v34 = (eps*k15)/(1 + DRG[t]^2/J15^2)
v35 = eps*(K11a + K11*CYCB[t])
v36 = E2FA*eps*K29*MASS[t]
v37 = eps*K33
v38 = eps*(E2FA*K7 + K7a)
v39 = eps*K9*DRG[t]
v4 = K10*CYCD[t]
v40 = eps*K5
v41 = eps*((k17*DRG[t]^2)/(J17^2*(1 + DRG[t]^2/J17^2)) + k17a*ERG[t])
v42 = eps*(K1a + (K1*CYCB[t]^2)/(J1^2*(1 + CYCB[t]^2/J1^2)))
v5 = K25*CYCE[t]*P27[t]
v6 = K25*CYCA[t]*P27[t]
v7 = k24*CYCD[t]*P27[t]
v8 = k24r*CD[t]
v9 = K30*CDc20[t]*CYCA[t]
end rate equations

begin parameters
E2FA = 0.0
E2FT = 5.0
FB = 2.0
FE = 25.0
GA = 0.3
GB = 1.0
GE = 0.0
HA = 0.5
HB = 1.0
HE = 0.5
J1 = 0.1
J13 = 0.005
J14 = 0.005
J15 = 0.1
J17 = 0.3
J3 = 0.01
J31 = 0.01
J32 = 0.01
J4 = 0.01
J8 = 0.1
K1 = 0.06
K10 = 0.5
K11 = 0.15
K11a = 0.0
K12 = 0.15
K13 = 0.5
K14 = 0.25
K19 = 2.0
K19a = 0.0
K1a = 0.01
K2 = 2.0
K20 = 1.0
K21 = 1.0
K22 = 0.1
K23 = 0.1
K23a = 0.0005
K25 = 100.0
K25R = 1.0
K26 = 1000.0
K26R = 20.0
K27 = 0.02
K28 = 0.02
K29 = 0.005
K2a = 0.005
K2aa = 0.1
K3 = 14.0
K30 = 2.0
K31 = 0.07
K32 = 0.18
K33 = 0.005
K34 = 0.005
K3a = 0.75
K4 = 4.0
K5 = 2.0
K6 = 10.0
K6a = 1.0
K7 = 0.06
K7a = 0.0
K8 = 0.2
K8a = 0.01
K9 = 0.25
LA = 3.0
LB = 5.0
LD = 3.3
LE = 5.0
MU = 0.0061
PP1T = 1.0
RBT = 10.0
YB = 0.05
YE = 1.0
eps = 1.0
k15 = 0.025
k16 = 0.025
k17 = 1.0
k17a = 0.035
k18 = 1.0
k24 = 100.0
k24r = 1.0
p = 1.0
s = 1.0
defaultcompartment = 1.0
end parameters

begin initial conditions
CA[0] = CAi
CD[0] = CDi
CDc20[0] = CDc20i
CDc20T[0] = CDc20Ti
CDh1[0] = CDh1i
CE[0] = CEi
CYCA[0] = CYCAi
CYCB[0] = CYCBi
CYCD[0] = CYCDi
CYCE[0] = CYCEi
DRG[0] = DRGi
E2F[0] = E2Fi
ERG[0] = ERGi
GM[0] = GMi
IEP[0] = IEPi
MASS[0] = MASSi
P27[0] = P27i
PPX[0] = PPXi
end initial conditions

begin initial values
CAi = 0.0434486
CDi = 0.0488772
CDc20i = 1e-07
CDc20Ti = 0.006466
CDh1i = 0.98702
CEi = 0.1380924
CYCAi = 0.3411
CYCBi = 0.00514
CYCDi = 0.401389
CYCEi = 1.09959
DRGi = 0.9005328
E2Fi = 4.25991
ERGi = 0.0121808
GMi = 1.032033
IEPi = 1.98e-05
MASSi = 1.28486
P27i = 0.011877
PPXi = 1.0
end initial values

begin assignment rules
L := (K26R + K20*(LA*CYCA[t] + LB*CYCB[t] + LD*CYCD[t] + LE*CYCE[t]))/K26
E2RBC := (2*E2FT*RBH)/(E2FT + L + RBH + (-4*E2FT*RBH + (E2FT + L + RBH)^2)^0.5)
RBH := RBT/(1 + (K20*(LA*CYCA[t] + LB*CYCB[t] + LD*CYCD[t] + LE*CYCE[t]))/(K19*PP1A + K19a*(-PP1A + PP1T)))
PP1A := PP1T/(1 + K21*(FB*CYCB[t] + FE*(CYCA[t] + CYCE[t])))
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units