#leloup1



#reactions
Mp_degradation:
 {1.0}Mp = $pool
  Cell * k_d * Mp + Cell * V_mP * Mp / (K_mP + Mp)
Mp_production:
 $pool = {1.0}Mp
  Cell * v_sP * pow(K_IP, n) / (pow(K_IP, n) + pow(Cn, n))
Mt_degradation:
 {1.0}Mt = $pool
  Cell * k_d * Mt + Cell * V_mT * Mt / (K_mT + Mt)
Mt_production:
 $pool = {1.0}Mt
  Cell * V_sT * pow(K_IT, n) / (pow(K_IT, n) + pow(Cn, n))
P0_degradation:
 {1.0}P0 = $pool
  Cell * k_d * P0
P0_production:
 $pool = {1.0}P0
  Cell * k_sP * Mp
P0_to_P1:
 {1.0}P0 = {1.0}P1
  Cell * V_1P * P0 / (K1_P + P0)
P1_degradation:
 {1.0}P1 = $pool
  Cell * k_d * P1
P1_to_P0:
 {1.0}P1 = {1.0}P0
  Cell * V_2P * P1 / (K_2P + P1)
P1_to_P2:
 {1.0}P1 = {1.0}P2
  Cell * V_3P * P1 / (K_3P + P1)
P2_degradation:
 {1.0}P2 = $pool
  Cell * k_d * P2 + Cell * V_dP * P2 / (K_dP + P2)
P2_to_P1:
 {1.0}P2 = {1.0}P1
  Cell * V_4P * P2 / (K_4P + P2)
PT_complex_degradation:
 {1.0}CC = $pool
  Cell * k_dC * CC
PT_complex_formation:
 {1.0}P2 + {1.0}T2 = {1.0}CC
  Cell * k3 * P2 * T2 - Cell * k4 * CC
PT_complex_nucleation:
 {1.0}CC = {1.0}Cn
  Cell * k1 * CC - compartment_0000002 * k2 * Cn
PTnucl_complex_degradation:
 {1.0}Cn = $pool
  compartment_0000002 * k_dN * Cn
T0_degradation:
 {1.0}T0 = $pool
  Cell * k_d * T0
T0_production:
 $pool = {1.0}T0
  Cell * k_sT * Mt
T0_to_T1:
 {1.0}T0 = {1.0}T1
  Cell * V_1T * T0 / (K_1T + T0)
T1_degradation:
 {1.0}T1 = $pool
  Cell * k_d * T1
T1_to_T0:
 {1.0}T1 = {1.0}T0
  Cell * V_2T * T1 / (K_2T + T1)
T1_to_T2:
 {1.0}T1 = {1.0}T2
  Cell * V_3T * T1 / (K_3T + T1)
T2_degradation:
 {1.0}T2 = $pool
  Cell * k_d * T2 + Cell * V_dT * T2 / (K_dT + T2)
T2_to_T1:
 {1.0}T2 = {1.0}T1
  Cell * V_4T * T2 / (K_4T + T2)


#parameters
V_dT = 2.0
V_mT = 0.7
P1_to_P2_K_3P = 2.0
P0_to_P1_K1_P = 2.0
P0_to_P1_V_1P = 8.0
T0_to_T1_K_1T = 2.0
T0_to_T1_V_1T = 8.0
P1_to_P0_K_2P = 2.0
P1_to_P0_V_2P = 1.0
T1_to_T0_K_2T = 2.0
T1_to_T0_V_2T = 1.0
P1_to_P2_V_3P = 8.0
T1_to_T2_K_3T = 2.0
T1_to_T2_V_3T = 8.0
P2_to_P1_K_4P = 2.0
P2_to_P1_V_4P = 1.0
T2_to_T1_K_4T = 2.0
T2_to_T1_V_4T = 1.0
P0_degradation_k_d = 0.01
T0_degradation_k_d = 0.01
P1_degradation_k_d = 0.01
T1_degradation_k_d = 0.01
P2_degradation_k_d = 0.01
P2_degradation_V_dP = 2.0
P2_degradation_K_dP = 0.2
T2_degradation_k_d = 0.01
T2_degradation_K_dT = 0.2
PT_complex_formation_k3 = 1.2
PT_complex_formation_k4 = 0.6
PT_complex_nucleation_k1 = 0.6
PT_complex_nucleation_k2 = 0.2
PT_complex_degradation_k_dC = 0.01
PTnucl_complex_degradation_k_dN = 0.01
Mp_production_v_sP = 1.0
Mp_production_K_IP = 1.0
Mp_production_n = 4.0
Mt_production_V_sT = 1.0
Mt_production_K_IT = 1.0
Mt_production_n = 4.0
P0_production_k_sP = 0.9
T0_production_k_sT = 0.9
Mp_degradation_k_d = 0.01
Mp_degradation_V_mP = 0.7
Mp_degradation_K_mP = 0.2
Mt_degradation_k_d = 0.01
Mt_degradation_K_mT = 0.2
Cell = 1.0
compartment_0000002 = 1.0


#initial values
CC = 0.0
Cn = 0.0
Mp = 0.0
Mt = 0.0
P0 = 0.0
P1 = 0.0
P2 = 0.0
T0 = 0.0
T1 = 0.0
T2 = 0.0


#assignment rules
!F Tt = CC + Cn + T0 + T1 + T2
!F Pt = CC + Cn + P0 + P1 + P2