#novak1



#reactions
R1:
 $pool = {1.0}cyclin
  k1AA
R10:
 {1.0}p_dimer = $pool
  k2 * p_dimer
R11:
 {1.0}p_dimer_p = {1.0}p_dimer
  kinh * p_dimer_p
R12:
 {1.0}dimer_p = {1.0}p_dimer_p
  kwee * dimer_p
R13:
 {1.0}p_dimer_p = {1.0}dimer_p
  k25 * p_dimer_p
R14:
 {1.0}p_dimer_p = $pool
  k2 * p_dimer_p
R15:
 {1.0}dimer_p = $pool
  k2 * dimer_p
R17:
 $pool = {1.0}cdc25_p
  ka * dimer_p * (total_cdc25 - cdc25_p) / (K_a + total_cdc25 - cdc25_p)
R18:
 {1.0}cdc25_p = $pool
  kbPPase * cdc25_p / (K_b + cdc25_p)
R19:
 $pool = {1.0}wee1_p
  ke * dimer_p * (total_wee1 - wee1_p) / (K_e + total_wee1 - wee1_p)
R2:
 {1.0}cyclin = $pool
  k2 * cyclin
R20:
 {1.0}wee1_p = $pool
  kfPPase * wee1_p / (K_f + wee1_p)
R21:
 $pool = {1.0}IE_p
  kg * dimer_p * (total_IE - IE_p) / (K_g + total_IE - IE_p)
R22:
 {1.0}IE_p = $pool
  khPPAse * IE_p / (K_h + IE_p)
R23:
 $pool = {1.0}UbE_star
  kc * IE_p * (total_UbE - UbE_star) / (K_c + total_UbE - UbE_star)
R24:
 {1.0}UbE_star = $pool
  kd_anti_IE * UbE_star / (K_d + UbE_star)
R3:
 {1.0}cyclin = {1.0}dimer
  k3 * cyclin * cdc2
R4:
 {1.0}dimer_p = {1.0}dimer
  kinh * dimer_p
R5:
 {1.0}dimer = {1.0}p_dimer
  kwee * dimer
R6:
 {1.0}dimer = {1.0}dimer_p
  kcak * dimer
R7:
 {1.0}dimer = $pool
  k2 * dimer
R8:
 {1.0}p_dimer = {1.0}dimer
  k25 * p_dimer
R9:
 {1.0}p_dimer = {1.0}p_dimer_p
  kcak * p_dimer


#parameters
K_a = 0.1
K_b = 0.1
K_c = 0.01
K_d = 0.01
K_e = 0.3
K_f = 0.3
K_g = 0.01
K_h = 0.01
V25_double_prime = 2.0
V25_prime = 0.1
V2_double_prime = 1.0
V2_prime = 0.015
Vwee_double_prime = 1.0
Vwee_prime = 0.1
k1AA = 1.0
k3 = 0.01
ka = 0.01
kbPPase = 0.125
kc = 0.1
kcak = 0.25
kd_anti_IE = 0.095
ke = 0.0133
kfPPase = 0.1
kg = 0.0065
khPPAse = 0.087
kinh = 0.025
total_IE = 1.0
total_UbE = 1.0
total_cdc2 = 100.0
total_cdc25 = 1.0
total_wee1 = 1.0
cytoplasm = 1.0


#initial values
IE_p = 0.0
UbE_star = 0.0
cdc25_p = 0.0
cyclin = 100.0
dimer = 0.0
dimer_p = 0.0
p_dimer = 0.0
p_dimer_p = 0.0
wee1_p = 0.0


#assignment rules
!F wee1 = total_wee1 - wee1_p
!F cdc25 = total_cdc25 - cdc25_p
!F cdc2 = total_cdc2 - (dimer + p_dimer + p_dimer_p + dimer_p)
!F UbE = total_UbE - UbE_star
!F IE = total_IE - IE_p
!F k2 = V2_prime * (total_UbE - UbE_star) + V2_double_prime * UbE_star
!F kwee = Vwee_prime * wee1_p + Vwee_double_prime * (total_wee1 - wee1_p)
!F Y15P = p_dimer + p_dimer_p
!F total_cyclin = cyclin + dimer + dimer_p + p_dimer + p_dimer_p
!F k25 = V25_prime * (total_cdc25 - cdc25_p) + V25_double_prime * cdc25_p