#obeyesekere1



#reactions
_1:
 $pool = {1.0}D_1
  aD_1 * (k_1 * GF_1 / (1 + k_1 * GF_1))
_7:
 {1.0}RS_1 = $pool
  pE_1 * RS_1 * E_1 / (qE_1 + RS_1 + E_1)
cycECDK2degradation_1:
 {1.0}E_1 = $pool
  dE_1 * X_1 * E_1
cyclebreak_1:
 {1.0}X_1 = $pool
  dX_1 * X_1
cycleprogression_1:
 $pool = {1.0}X_1
  aX_1 * E_1 + f_1 * E2F_1 + g_1 * pow(X_1, 2) * E_1
cyclinCDK4degradation_1:
 {1.0}D_1 = $pool
  dD_1 * E_1 * D_1
cyclin_1:
 $pool = {1.0}E_1
  aE_1 * (1 + af_1 * E2F_1)
pRBE2Fcomplexassociation_1:
 {1.0}R_1 = {1.0}RS_1
  pS_1 * E2F_1 * R_1
pRBE2FcomplexdeassociationviacycDCDK4_1:
 {1.0}RS_1 = $pool
  pD_1 * RS_1 * D_1 / (qD_1 + RS_1 + D_1)
pRBpdephosphorylation_1:
 $pool = {1.0}R_1
  pX_1 * RP_1 * X_1 / (qX_1 + RP_1 + X_1)


#parameters
GF_1 = 6.3
RT_1 = 2.5
aD_1 = 0.4
aE_1 = 0.16
aX_1 = 0.08
af_1 = 0.9
dD_1 = 0.4
dE_1 = 0.2
dX_1 = 1.04
f_1 = 0.2
g_1 = 0.528
k_1 = 0.05
pD_1 = 0.48
pE_1 = 0.096
pS_1 = 0.6
pX_1 = 0.48
qD_1 = 0.6
qE_1 = 0.6
qX_1 = 0.8
theta_1 = 1.5
cell_1 = 1.0


#initial values
D_1 = 0.1
E_1 = 0.6
RS_1 = 1.0
R_1 = 0.5
X_1 = 0.7


#assignment rules
!F RP_1 = RT_1 - RS_1 - R_1
!F E2F_1 = theta_1 - RS_1
!F unpho_RB = R_1 + RS_1