#goldbeter3



#reactions
reaction1:
 $pool = {1.0}C
  cell * vi
reaction2:
 {1.0}C = $pool
  C * cell * kd
reaction3:
 {1.0}C = $pool
  C * cell * vd * X * pow(C + Kd, -1)
reaction4:
 {1.0}MI = {1.0}M
  cell * MI * V1 * pow(K1 + MI, -1)
reaction5:
 {1.0}M = {1.0}MI
  cell * M * V2 * pow(K2 + M, -1)
reaction6:
 {1.0}XI = {1.0}X
  cell * V3 * XI * pow(K3 + XI, -1)
reaction7:
 {1.0}X = {1.0}XI
  cell * V4 * X * pow(K4 + X, -1)


#parameters
Kc = 0.5
VM1 = 3.0
VM3 = 1.0
reaction3_vd = 0.25
reaction1_vi = 0.025
reaction2_kd = 0.01
reaction3_Kd = 0.02
reaction4_K1 = 0.005
reaction5_V2 = 1.5
reaction5_K2 = 0.005
reaction6_K3 = 0.005
reaction7_K4 = 0.005
reaction7_V4 = 0.5
cell = 1.0


#initial values
C = 0.01
M = 0.01
MI = 0.99
X = 0.01
XI = 0.99


#assignment rules
!F V3 = M * VM3
!F V1 = C * VM1 * pow(C + Kc, -1.0)