#conradie2



#reactions
v1:
 {1.0}ERG = $pool
  k16 * ERG
v10:
 {1.0}CA = {1.0}P27
  K30 * CA * CDc20
v11:
 {1.0}CE = {1.0}CYCE + {1.0}P27
  K25R * CE
v12:
 {1.0}CA = {1.0}CYCA + {1.0}P27
  K25R * CA
v13:
 {1.0}CE = {1.0}P27
  V8 * CE
v14:
 {1.0}CYCE = $pool
  V8 * CYCE
v15:
 {1.0}P27 = $pool
  V6 * P27
v16:
 {1.0}CE = {1.0}CYCE
  V6 * CE
v17:
 {1.0}CD = {1.0}CYCD
  V6 * CD
v18:
 {1.0}CA = {1.0}CYCA
  V6 * CA
v19:
 {1.0}CYCB = $pool
  V2 * CYCB
v2:
 {1.0}DRG = $pool
  k18 * DRG
v20:
 $pool = {1.0}CDh1
  (K3a + K3 * CDc20) * (1 - CDh1) / (1 + J3 - CDh1)
v21:
 {1.0}CDh1 = $pool
  V4 * CDh1 / (J4 + CDh1)
v22:
 {1.0}PPX = $pool
  K34 * PPX
v23:
 $pool = {1.0}IEP
  K31 * CYCB * (1 - IEP) / (1 + J31 - IEP)
v24:
 {1.0}IEP = $pool
  K32 * IEP * PPX / (J32 + IEP)
v25:
 {1.0}CDc20T = $pool
  K12 * CDc20T
v26:
 $pool = {1.0}CDc20
  K13 * (-CDc20 + CDc20T) * IEP / (J13 - CDc20 + CDc20T)
v27:
 {1.0}CDc20 = $pool
  K14 * CDc20 / (J14 + CDc20)
v28:
 {1.0}CDc20 = $pool
  K12 * CDc20
v29:
 {1.0}var5 = {1.0}var1 + {1.0}var2
  K20 * (LA * CYCA + LB * CYCB + LD * (CD + CYCD) + LE * CYCE) * var5
v3:
 {1.0}CD = {1.0}P27
  K10 * CD
v30:
 {1.0}var6 = {1.0}var1 + {1.0}var3
  K20 * (LA * CYCA + LB * CYCB + LD * (CD + CYCD) + LE * CYCE) * var6
v31:
 $pool = {1.0}GM
  K27 * MASS * r31switch
v32:
 {1.0}GM = $pool
  K28 * GM
v33:
 $pool = {1.0}MASS
  eps * MU * GM
v34:
 $pool = {1.0}ERG
  eps * k15 / (1 + pow(DRG, 2) / pow(J15, 2))
v35:
 $pool = {1.0}CDc20T
  eps * (K11a + K11 * CYCB)
v36:
 $pool = {1.0}CYCA
  eps * K29 * MASS * var2
v37:
 $pool = {1.0}PPX
  eps * K33
v38:
 $pool = {1.0}CYCE
  eps * (K7a + K7 * var2)
v39:
 $pool = {1.0}CYCD
  eps * K9 * DRG
v4:
 {1.0}CYCD = $pool
  K10 * CYCD
v40:
 $pool = {1.0}P27
  eps * K5
v41:
 $pool = {1.0}DRG
  eps * (k17 * pow(DRG, 2) / (pow(J17, 2) * (1 + pow(DRG, 2) / pow(J17, 2))) + k17a * ERG)
v42:
 $pool = {1.0}CYCB
  eps * (K1a + K1 * pow(CYCB, 2) / (pow(J1, 2) * (1 + pow(CYCB, 2) / pow(J1, 2))))
v43:
 {1.0}var4 = {1.0}var1
  K20 * (LA * CYCA + LB * CYCB + LD * (CD + CYCD) + LE * CYCE) * var4
v44:
 {1.0}var1 = {1.0}var4
  (K19 * PP1A + K19a * (PP1T - PP1A)) * var1
v45:
 {1.0}var5 = {1.0}var2 + {1.0}var4
  K26R * var5
v46:
 {1.0}var2 = {1.0}var3
  (K23a + K23 * (CYCA + CYCB)) * var2
v47:
 {1.0}var3 = {1.0}var2
  K22 * var3
v48:
 {1.0}var2 + {1.0}var4 = {1.0}var5
  K26 * var2 * var4
v49:
 {1.0}var6 = {1.0}var3 + {1.0}var4
  K26R * var6
v5:
 {1.0}CYCE + {1.0}P27 = {1.0}CE
  K25 * CYCE * P27
v50:
 {1.0}var3 + {1.0}var4 = {1.0}var6
  K26 * var3 * var4
v51:
 {1.0}var6 = {1.0}var5
  K22 * var6
v52:
 {1.0}var5 = {1.0}var6
  (K23a + K23 * (CYCA + CYCB)) * var5
v6:
 {1.0}CYCA + {1.0}P27 = {1.0}CA
  K25 * CYCA * P27
v7:
 {1.0}CYCD + {1.0}P27 = {1.0}CD
  k24 * CYCD * P27
v8:
 {1.0}CD = {1.0}CYCD + {1.0}P27
  k24r * CD
v9:
 {1.0}CYCA = $pool
  K30 * CDc20 * CYCA


#parameters
E2FT = 5.0
FB = 2.0
FE = 25.0
Flag = 1.0
GA = 0.3
GB = 1.0
GE = 0.0
HA = 0.5
HB = 1.0
HE = 0.5
J8 = 0.1
K10 = 5.0
K12 = 1.5
K2 = 20.0
K20 = 10.0
K21 = 1.0
K22 = 1.0
K23 = 1.0
K23a = 0.005
K25 = 1000.0
K25R = 10.0
K26 = 10000.0
K26R = 200.0
K2a = 0.05
K2aa = 1.0
K30 = 20.0
K4 = 40.0
K6 = 100.0
K6a = 10.0
K8 = 2.0
K8a = 0.1
LA = 3.0
LB = 5.0
LD = 3.3
LE = 5.0
PP1T = 1.0
RBT = 10.0
YB = 0.05
YE = 1.0
eps = 1.0
r31switch = 1.0
v1_k16 = 0.25
v2_k18 = 10.0
v7_k24 = 1000.0
v8_k24r = 10.0
v20_J3 = 0.01
v20_K3 = 140.0
v20_K3a = 7.5
v21_J4 = 0.01
v22_K34 = 0.05
v23_J31 = 0.01
v23_K31 = 0.7
v24_J32 = 0.01
v24_K32 = 1.8
v26_J13 = 0.005
v26_K13 = 5.0
v27_J14 = 0.005
v27_K14 = 2.5
v31_K27 = 0.2
v32_K28 = 0.2
v33_MU = 0.061
v34_J15 = 0.1
v34_k15 = 0.25
v35_K11 = 1.5
v35_K11a = 0.0
v36_K29 = 0.05
v37_K33 = 0.05
v38_K7 = 0.6
v38_K7a = 0.0
v39_K9 = 2.5
v40_K5 = 20.0
v41_J17 = 0.3
v41_k17 = 10.0
v41_k17a = 0.35
v42_J1 = 0.1
v42_K1 = 0.6
v42_K1a = 0.1
v44_K19 = 20.0
v44_K19a = 0.0
cell = 1.0


#initial values
CA = 0.0356927
CD = 0.010976
CDc20 = 0.00220177
CDc20T = 2.36733
CDh1 = 0.000653278
CE = 0.000542587
CYCA = 1.4094
CYCB = 2.72898
CYCD = 0.43929
CYCE = 0.0229112
DRG = 0.900533
ERG = 0.0121809
GM = 1.35565
IEP = 0.154655
MASS = 1.68776
P27 = 0.00922806
PPX = 1.0
var1 = 9.97574
var2 = 0.989986
var3 = 3.98594
var4 = 0.000190871
var5 = 0.00478911
var6 = 0.0192822


#assignment rules
!F P27T = CA + CD + CE + P27
!F CYCAT = CA + CYCA
!F CYCDT = CD + CYCD
!F CYCET = CE + CYCE
!F V8 = K8a + K8 * (YB * CYCB + YE * (CYCA + CYCE)) / (CYCET + J8)
!F V6 = K6a + K6 * (HA * CYCA + HB * CYCB + HE * CYCE)
!F V4 = K4 * (GA * CYCA + GB * CYCB + GE * CYCE)
!F V2 = K2aa * CDc20 + K2a * (1.0 - CDh1) + K2 * CDh1
!F PP1A = PP1T / (1.0 + K21 * (FB * CYCB + FE * (CYCA + CYCE)))