begin name
mandlik1
end name

begin reactions
ia1bindforward {1.0}IPTG + {1.0}ia1ActiveTF = {1.0}ia1InactiveTF
ia1bindreverse {1.0}ia1InactiveTF = {1.0}IPTG + {1.0}ia1ActiveTF
pp1v1 {1.0}$EXTERNAL = {1.0}proSLS1
pp1v2 {1.0}proSLS1 = {1.0}$EXTERNAL
pp2v1 {1.0}$EXTERNAL = {1.0}proAUR1
pp2v2 {1.0}proAUR1 = {1.0}$EXTERNAL
pp3v1 {1.0}$EXTERNAL = {1.0}proSLS4
pp3v2 {1.0}proSLS4 = {1.0}$EXTERNAL
pp4v1 {1.0}$EXTERNAL = {1.0}proLACI
pp4v2 {1.0}proLACI = {1.0}$EXTERNAL
pp5v1 {1.0}$EXTERNAL = {1.0}proLAMDAR
pp5v2 {1.0}proLAMDAR = {1.0}$EXTERNAL
pp6v1 {1.0}$EXTERNAL = {1.0}proTETR
pp6v2 {1.0}proTETR = {1.0}$EXTERNAL
end reactions

begin rate equations
ia1bindforward = ia1ActiveTF[t]*IPTG[t]
ia1bindreverse = ia1Kd*ia1InactiveTF[t]
pp1v1 = 1
pp1v2 = proSLS1degradationrate*proSLS1[t]
pp2v1 = AUR1*pAUR1strength*pp2translationrate
pp2v2 = proAUR1degradationrate*proAUR1[t]
pp3v1 = pp3translationrate*pSLS4strength*SLS4
pp3v2 = proSLS4degradationrate*proSLS4[t]
pp4v1 = LACI*p2strength*pp4translationrate
pp4v2 = proLACIdegradationrate*proLACI[t]
pp5v1 = LAMDAR*p1strength*pp5translationrate
pp5v2 = proLAMDARdegradationrate*proLAMDAR[t]
pp6v1 = p3strength*pp6translationrate*TETR
pp6v2 = proTETRdegradationrate*proTETR[t]
end rate equations

begin parameters
ia1Kd = 0.001
ia1repressionKd = 1.0
ia1repressionh = 2.0
p1strength = 2.5
p2strength = 2.4
p3strength = 2.5
pAUR1strength = 2.014
pSLS1strength = 2.5
pSLS4strength = 2.57
pp2translationrate = 1.0
pp3translationrate = 1.0
pp4translationrate = 1.0
pp5translationrate = 1.0
pp6translationrate = 1.0
proAUR1degradationrate = 0.1
proLACIdegradationrate = 0.07
proLAMDARdegradationrate = 0.0415
proSLS1degradationrate = 0.07
proSLS4degradationrate = 0.04
proTETRdegradationrate = 0.08
ta1Kd = 1.0
ta1h = 2.0
tr1Kd = 1.0
tr1h = 2.0
tr2Kd = 1.0
tr2h = 2.0
tr3Kd = 1.0
tr3h = 2.0
tr4Kd = 1.0
tr4h = 2.0
tr5Kd = 1.0
tr5h = 2.0
tr6Kd = 1.0
tr6h = 2.0
DefaultCompartment = 1.0
end parameters

begin initial conditions
IPTG[0] = IPTGi
ia1ActiveTF[0] = ia1ActiveTFi
ia1InactiveTF[0] = ia1InactiveTFi
p1[0] = p1i
p2[0] = p2i
p3[0] = p3i
pAUR1[0] = pAUR1i
pSLS1[0] = pSLS1i
pSLS4[0] = pSLS4i
proAUR1[0] = proAUR1i
proLACI[0] = proLACIi
proLAMDAR[0] = proLAMDARi
proSLS1[0] = proSLS1i
proSLS4[0] = proSLS4i
proTETR[0] = proTETRi
ter1[0] = ter1i
ter2[0] = ter2i
ter3[0] = ter3i
ter4[0] = ter4i
ter5[0] = ter5i
ter6[0] = ter6i
end initial conditions

begin initial values
IPTGi = 10000.0
ia1ActiveTFi = 100.0
ia1InactiveTFi = 100.0
p1i = 10.0
p2i = 10.0
p3i = 10.0
pAUR1i = 10.0
pSLS1i = 10.0
pSLS4i = 10.0
proAUR1i = 0.0
proLACIi = 0.0
proLAMDARi = 0.0
proSLS1i = 0.0
proSLS4i = 0.0
proTETRi = 0.0
ter1i = 10.0
ter2i = 10.0
ter3i = 10.0
ter4i = 10.0
ter5i = 10.0
ter6i = 10.0
end initial values

begin assignment rules
rs6 := (1 + (proTETR[t]/tr4Kd)^tr4h)^(-1)
rs4 := (1 + (proLAMDAR[t]/tr5Kd)^tr5h)^(-1)
AUR1 := pAUR1strength/(1 + (proSLS1[t]/tr1Kd)^tr1h)
LACI := p2strength/(1 + (proLAMDAR[t]/tr5Kd)^tr5h)
SLS4 := pSLS4strength/(1 + (proAUR1[t]/tr2Kd)^tr2h)
rs5 := (1 + (proLACI[t]/tr3Kd)^tr3h)^(-1)
rs1 := (1 + (proSLS1[t]/tr1Kd)^tr1h)^(-1)
rs2 := (1 + (proAUR1[t]/tr2Kd)^tr2h)^(-1)
ope1 := (1 + (ia1ActiveTF[t]/ia1repressionKd)^ia1repressionh)^(-1)
SLS1 := (pSLS1strength*(proSLS4[t]/ta1Kd)^ta1h)/((1 + (ia1ActiveTF[t]/ia1repressionKd)^ia1repressionh)*(1 + (proLACI[t]/tr6Kd)^tr6h)*(1 + (proSLS4[t]/ta1Kd)^ta1h))
as1 := (proSLS4[t]/ta1Kd)^ta1h/(1 + (proSLS4[t]/ta1Kd)^ta1h)
TETR := p3strength/(1 + (proLACI[t]/tr3Kd)^tr3h)
LAMDAR := p1strength/(1 + (proTETR[t]/tr4Kd)^tr4h)
rs3 := (1 + (proLACI[t]/tr6Kd)^tr6h)^(-1)
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units