#vernoux2



#reactions
reac_DIA:
 {1.0}A + {1.0}I = {1.0}D_IA
  k_IA * I * A - kprime_IA * D_IA
reac_DIAxA:
 {1.0}D_IA = {1.0}A
  gamma_I * d_I * K_aux * aux / (K_aux * aux + 1) * D_IA
reac_DII:
 {1.0}I + {1.0}I = {1.0}D_II
  k_II * I * I - kprime_II * D_II
reac_DIIxI:
 {1.0}D_II = {1.0}I
  gamma_I * d_I * K_aux * aux / (K_aux * aux + 1) * D_II
reac_degrA:
 {1.0}A = $pool
  d_A * A
reac_degrDIA:
 {1.0}D_IA = $pool
  d_IA * D_IA
reac_degrDII:
 {1.0}D_II = $pool
  d_II * D_II
reac_degrI:
 {1.0}I = $pool
  gamma_I * d_I * K_aux * aux / (K_aux * aux + 1) * I
reac_degrR:
 {1.0}R = $pool
  d_r * R
reac_prodA:
 $pool = {1.0}A
  pi_A
reac_prodI:
 $pool = {1.0}I
  pi_I * R
reac_prodR:
 $pool = {1.0}R
  (1 + f_c / B_d * A * (1 + w_A * f_A * A / B_d)) / (1 + A / B_d * (1 + w_A * A / B_d) + w_I * A * I / (K_IA * B_d) + w_D * D_IA / B_d + k_Am)


#parameters
B_d = 100.0
K_IA = 10.0
K_II = 10.0
K_aux = 1.0
d_A = 0.003
d_I = 0.05
d_IA = 0.003
d_II = 0.003
d_r = 0.007
f_A = 10.0
f_c = 10.0
gamma_I = 10.0
k_Am = 10.0
k_IA = 1.0
k_II = 1.0
pi_A = 1.0
pi_I = 1.0
w_A = 10.0
w_D = 10.0
w_I = 10.0
compartment_0000001 = 1.0


#initial values
A = 10.0
D_IA = 10.0
D_II = 10.0
I = 10.0
R = 1.0


#assignment rules
!F aux_basal = 1.0 / (K_aux * (gamma_I - 1.0))
!F kprime_II = K_II * k_II
!F kprime_IA = K_IA * k_IA
!F aux = 1.11 * sin(2.0 * 3.1416 / 800.0 * (time - 200.0)) + 1.11