#bai1



#reactions
E2Fdegradationviacellcycleprogression_1:
 {1.0}theta_1 = $pool
  cell * dtheta_1 * (X_1 / (qtheta_1 + X_1)) * theta_1
_1:
 $pool = {1.0}D_1
  cell * aD_1 * (GF_1 / (pow(k1_1, -1) + GF_1))
_2:
 {1.0}RS_1 = {1.0}theta_1
  cell * (pD_1 * RS_1 * D_1 / (qD_1 + RS_1 + D_1))
cellcycleprogression_1:
 $pool = {1.0}X_1
  cell * (aX_1 * E_1 + f_1 * theta_1 + g_1 * pow(X_1, 2) * E_1)
cellcycleprogressionslow_1:
 {1.0}X_1 = $pool
  cell * dX_1 * X_1
cyclinD_1:
 {1.0}D_1 = $pool
  cell * dD_1 * E_1 * D_1
cyclinEdegradation_1:
 {1.0}E_1 = $pool
  cell * dE_1 * X_1 * E_1
cyclinEsynthesis_1:
 $pool = {1.0}E_1
  cell * aE_1 * (GF_1 / (pow(k2_1, -1) + GF_1) + aF_1 * theta_1)
null2_1:
 {1.0}RS_1 = {1.0}theta_1
  cell * (pE_1 * RS_1 * E_1 / (qE_1 + RS_1 + E_1))
null3_1:
 $pool = {1.0}theta_1
  cell * atheta_1 * (GF_1 / (pow(k3_1, -1) + GF_1) + fC_1_1 * theta_1)
pRBdeplation_1:
 {1.0}R_1 + {1.0}theta_1 = {1.0}RS_1
  cell * Mass_Action_2_1(pS_1, R_1, theta_1)
pRBsynthesis_1:
 $pool = {1.0}R_1
  cell * (pX_1 * (RT_1 - RS_1 - R_1) * X_1 / (qX_1 + RT_1 - RS_1 - R_1 + X_1))


#parameters
GF_1 = 6.3
RT_1 = 2.5
aD_1 = 0.4
aE_1 = 0.16
aF_1 = 0.9
aX_1 = 0.08
atheta_1 = 0.05
dD_1 = 0.4
dE_1 = 0.2
dX_1 = 1.04
dtheta_1 = 0.12
fC_1_1 = 0.003
f_1 = 0.35
g_1 = 0.528
k1_1 = 0.05
k2_1 = 1000.0
k3_1 = 1.5
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
qtheta_1 = 0.3
cell = 1.0


#initial values
D_1 = 0.0
E_1 = 0.0
RS_1 = 0.0
R_1 = 0.0
X_1 = 0.0
theta_1 = 0.0


#assignment rules
!F Rb_phos = RT_1 - RS_1 - R_1