ia1_bind_forward

ia1_bind_forward

IPTG + ia1_ActiveTF > ia1_InactiveTF

ia1_bind_reverse

ia1_bind_reverse

ia1_InactiveTF > IPTG + ia1_ActiveTF

pp1_v1

pp1_v1

∅ > proSLS1

pp1_v2

pp1_v2

proSLS1 > ∅

pp2_v1

pp2_v1

∅ > proAUR1

pp2_v2

pp2_v2

proAUR1 > ∅

pp3_v1

pp3_v1

∅ > proSLS4

pp3_v2

pp3_v2

proSLS4 > ∅

pp4_v1

pp4_v1

∅ > proLACI

pp4_v2

pp4_v2

proLACI > ∅

pp5_v1

pp5_v1

∅ > proLAMDAR

pp5_v2

pp5_v2

proLAMDAR > ∅

pp6_v1

pp6_v1

∅ > proTETR

pp6_v2

pp6_v2

proTETR > ∅

IPTG*

ia1_ActiveTF*

ia1_InactiveTF*

p1*

p2*

p3*

pAUR1*

pSLS1*

pSLS4*

proAUR1*

proLACI*

proLAMDAR*

proSLS1*

proSLS4*

proTETR*

ter1*

ter2*

ter3*

ter4*

ter5*

ter6*

Assignment rules

rs6 = 1.0 / (1.0 + pow(proTETR / tr4_Kd, tr4_h))

rs4 = 1.0 / (1.0 + pow(proLAMDAR / tr5_Kd, tr5_h))

AUR1 = pAUR1_strength * 1.0 / (1.0 + pow(proSLS1 / tr1_Kd, tr1_h))

LACI = p2_strength * 1.0 / (1.0 + pow(proLAMDAR / tr5_Kd, tr5_h))

SLS4 = pSLS4_strength * 1.0 / (1.0 + pow(proAUR1 / tr2_Kd, tr2_h))

rs5 = 1.0 / (1.0 + pow(proLACI / tr3_Kd, tr3_h))

rs1 = 1.0 / (1.0 + pow(proSLS1 / tr1_Kd, tr1_h))

rs2 = 1.0 / (1.0 + pow(proAUR1 / tr2_Kd, tr2_h))

ope1 = 1.0 / (1.0 + pow(ia1_ActiveTF / ia1_repression_Kd, ia1_repression_h))

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))

as1 = (1.0 + pow(proSLS4 / ta1_Kd, ta1_h) - 1.0) / (1.0 + pow(proSLS4 / ta1_Kd, ta1_h))

TETR = p3_strength * 1.0 / (1.0 + pow(proLACI / tr3_Kd, tr3_h))

LAMDAR = p1_strength * 1.0 / (1.0 + pow(proTETR / tr4_Kd, tr4_h))

rs3 = 1.0 / (1.0 + pow(proLACI / tr6_Kd, tr6_h))

Function definitions

Note that constraints are not enforced in simulations. It remains the responsibility of the user to verify that simulation results satisfy these constraints.