#novak

FIX: p s

#reactions
v_1:
 {1.0}ERG = {1.0}p
  k16*ERG
v_10:
 {1.0}CA = {1.0}P27
  K30*CDc20*CA
v_11:
 {1.0}CE = {1.0}P27 + {1.0}CYCE
  K25R*CE
v_12:
 {1.0}CA = {1.0}P27 + {1.0}CYCA
  K25R*CA
v_13:
 {1.0}CE = {1.0}P27
  v_8*CE
v_14:
 {1.0}CYCE = {1.0}p
  v_8*CYCE
v_15:
 {1.0}P27 = {1.0}p
  v_6*P27
v_16:
 {1.0}CE = {1.0}CYCE
  v_6*CE
v_17:
 {1.0}CD = {1.0}CYCD
  v_6*CD
v_18:
 {1.0}CA = {1.0}CYCA
  v_6*CA
v_19:
 {1.0}CYCB = {1.0}p
  v_2*CYCB
v_2:
 {1.0}DRG = {1.0}p
  k18*DRG
v_20:
 {1.0}s = {1.0}CDh1
  (K3a + K3*CDc20)*(1 - CDh1)/(J3 + 1 - CDh1)
v_21:
 {1.0}CDh1 = {1.0}p
  v_4*CDh1/(J4 + CDh1)
v_22:
 {1.0}PPX = {1.0}p
  K34*PPX
v_23:
 {1.0}s = {1.0}IEP
  K31*CYCB*(1 - IEP)/(J31 + 1 - IEP)
v_24:
 {1.0}IEP = {1.0}p
  K32*PPX*IEP/(J32 + IEP)
v_25:
 {1.0}CDc20T = {1.0}p
  K12*CDc20T
v_26:
 {1.0}s = {1.0}CDc20
  K13*IEP*(CDc20T - CDc20)/(J13 + CDc20T - CDc20)
v_27:
 {1.0}CDc20 = {1.0}p
  K14*CDc20/(J14 + CDc20)
v_28:
 {1.0}CDc20 = {1.0}p
  K12*CDc20
v_29:
 {1.0}s = {1.0}E2F
  K22*(E2FT - E2F)
v_3:
 {1.0}CD = {1.0}P27
  K10*CD
v_30:
 {1.0}E2F = {1.0}p
  (K23a + K23*(CYCA + CYCB))*E2F
v_31:
 {1.0}s = {1.0}GM
  K27*MASS*piecewise(0, RBH/RBT > 0.8, 1)
v_32:
 {1.0}GM = {1.0}p
  K28*GM
v_33:
 {1.0}s = {1.0}MASS
  eps*MU*GM
v_34:
 {1.0}s = {1.0}ERG
  eps*k15/(1 + (DRG/J15)**2)
v_35:
 {1.0}s = {1.0}CDc20T
  eps*(K11a + K11*CYCB)
v_36:
 {1.0}s = {1.0}CYCA
  eps*K29*E2FA*MASS
v_37:
 {1.0}s = {1.0}PPX
  eps*K33
v_38:
 {1.0}s = {1.0}CYCE
  eps*(K7a + K7*E2FA)
v_39:
 {1.0}s = {1.0}CYCD
  eps*K9*DRG
v_4:
 {1.0}CYCD = {1.0}p
  K10*CYCD
v_40:
 {1.0}s = {1.0}P27
  eps*K5
v_41:
 {1.0}s = {1.0}DRG
  eps*(k17a*ERG + k17*(DRG/J17)**2/(1 + (DRG/J17)**2))
v_42:
 {1.0}s = {1.0}CYCB
  eps*(K1a + K1*(CYCB/J1)**2/(1 + (CYCB/J1)**2))
v_5:
 {1.0}P27 + {1.0}CYCE = {1.0}CE
  K25*P27*CYCE
v_6:
 {1.0}P27 + {1.0}CYCA = {1.0}CA
  K25*P27*CYCA
v_7:
 {1.0}P27 + {1.0}CYCD = {1.0}CD
  k24*CYCD*P27
v_8:
 {1.0}CD = {1.0}P27 + {1.0}CYCD
  k24r*CD
v_9:
 {1.0}CYCA = {1.0}p
  K30*CDc20*CYCA


#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
default_compartment = 1.0


#initial values
CA = 0.0434486
CD = 0.0488772
CDc20 = 1e-07
CDc20T = 0.006466
CDh1 = 0.98702
CE = 0.1380924
CYCA = 0.3411
CYCB = 0.00514
CYCD = 0.401389
CYCE = 1.09959
DRG = 0.9005328
E2F = 4.25991
ERG = 0.0121808
GM = 1.032033
IEP = 1.98e-05
MASS = 1.28486
P27 = 0.011877
PPX = 1.0


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