#mandlik1



#reactions
ia1_bind_forward:
 {1.0}IPTG + {1.0}ia1_ActiveTF = {1.0}ia1_InactiveTF
  IPTG * ia1_ActiveTF
ia1_bind_reverse:
 {1.0}ia1_InactiveTF = {1.0}IPTG + {1.0}ia1_ActiveTF
  ia1_Kd * ia1_InactiveTF
pp1_v1:
 $pool = {1.0}proSLS1
  1
pp1_v2:
 {1.0}proSLS1 = $pool
  proSLS1_degradation_rate * proSLS1
pp2_v1:
 $pool = {1.0}proAUR1
  pp2_translation_rate * pAUR1_strength * AUR1
pp2_v2:
 {1.0}proAUR1 = $pool
  proAUR1_degradation_rate * proAUR1
pp3_v1:
 $pool = {1.0}proSLS4
  pp3_translation_rate * pSLS4_strength * SLS4
pp3_v2:
 {1.0}proSLS4 = $pool
  proSLS4_degradation_rate * proSLS4
pp4_v1:
 $pool = {1.0}proLACI
  pp4_translation_rate * p2_strength * LACI
pp4_v2:
 {1.0}proLACI = $pool
  proLACI_degradation_rate * proLACI
pp5_v1:
 $pool = {1.0}proLAMDAR
  pp5_translation_rate * p1_strength * LAMDAR
pp5_v2:
 {1.0}proLAMDAR = $pool
  proLAMDAR_degradation_rate * proLAMDAR
pp6_v1:
 $pool = {1.0}proTETR
  pp6_translation_rate * p3_strength * TETR
pp6_v2:
 {1.0}proTETR = $pool
  proTETR_degradation_rate * proTETR


#parameters
ia1_Kd = 0.001
ia1_repression_Kd = 1.0
ia1_repression_h = 2.0
p1_strength = 2.5
p2_strength = 2.4
p3_strength = 2.5
pAUR1_strength = 2.014
pSLS1_strength = 2.5
pSLS4_strength = 2.57
pp2_translation_rate = 1.0
pp3_translation_rate = 1.0
pp4_translation_rate = 1.0
pp5_translation_rate = 1.0
pp6_translation_rate = 1.0
proAUR1_degradation_rate = 0.1
proLACI_degradation_rate = 0.07
proLAMDAR_degradation_rate = 0.0415
proSLS1_degradation_rate = 0.07
proSLS4_degradation_rate = 0.04
proTETR_degradation_rate = 0.08
ta1_Kd = 1.0
ta1_h = 2.0
tr1_Kd = 1.0
tr1_h = 2.0
tr2_Kd = 1.0
tr2_h = 2.0
tr3_Kd = 1.0
tr3_h = 2.0
tr4_Kd = 1.0
tr4_h = 2.0
tr5_Kd = 1.0
tr5_h = 2.0
tr6_Kd = 1.0
tr6_h = 2.0
DefaultCompartment = 1.0


#initial values
IPTG = 10000.0
ia1_ActiveTF = 100.0
ia1_InactiveTF = 100.0
p1 = 10.0
p2 = 10.0
p3 = 10.0
pAUR1 = 10.0
pSLS1 = 10.0
pSLS4 = 10.0
proAUR1 = 0.0
proLACI = 0.0
proLAMDAR = 0.0
proSLS1 = 0.0
proSLS4 = 0.0
proTETR = 0.0
ter1 = 10.0
ter2 = 10.0
ter3 = 10.0
ter4 = 10.0
ter5 = 10.0
ter6 = 10.0


#assignment rules
!F rs6 = 1.0 / (1.0 + pow(proTETR / tr4_Kd, tr4_h))
!F rs4 = 1.0 / (1.0 + pow(proLAMDAR / tr5_Kd, tr5_h))
!F AUR1 = pAUR1_strength * 1.0 / (1.0 + pow(proSLS1 / tr1_Kd, tr1_h))
!F LACI = p2_strength * 1.0 / (1.0 + pow(proLAMDAR / tr5_Kd, tr5_h))
!F SLS4 = pSLS4_strength * 1.0 / (1.0 + pow(proAUR1 / tr2_Kd, tr2_h))
!F rs5 = 1.0 / (1.0 + pow(proLACI / tr3_Kd, tr3_h))
!F rs1 = 1.0 / (1.0 + pow(proSLS1 / tr1_Kd, tr1_h))
!F rs2 = 1.0 / (1.0 + pow(proAUR1 / tr2_Kd, tr2_h))
!F ope1 = 1.0 / (1.0 + pow(ia1_ActiveTF / ia1_repression_Kd, ia1_repression_h))
!F SLS1 = pSLS1_strength * (1.0 + pow(proSLS4 / ta1_Kd, ta1_h) - 1.0) / (1.0 + pow(proSLS4 / ta1_Kd, ta1_h)) * 1.0 / (1.0 + pow(ia1_ActiveTF / ia1_repression_Kd, ia1_repression_h)) * 1.0 / (1.0 + pow(proLACI / tr6_Kd, tr6_h))
!F as1 = (1.0 + pow(proSLS4 / ta1_Kd, ta1_h) - 1.0) / (1.0 + pow(proSLS4 / ta1_Kd, ta1_h))
!F TETR = p3_strength * 1.0 / (1.0 + pow(proLACI / tr3_Kd, tr3_h))
!F LAMDAR = p1_strength * 1.0 / (1.0 + pow(proTETR / tr4_Kd, tr4_h))
!F rs3 = 1.0 / (1.0 + pow(proLACI / tr6_Kd, tr6_h))