begin name
elowitz1
end name

begin reactions
Reaction1 {1.0}X = {1.0}$EXTERNAL
Reaction10 {1.0}$EXTERNAL = {1.0}X
Reaction11 {1.0}$EXTERNAL = {1.0}Y
Reaction12 {1.0}$EXTERNAL = {1.0}Z
Reaction2 {1.0}Y = {1.0}$EXTERNAL
Reaction3 {1.0}Z = {1.0}$EXTERNAL
Reaction4 {1.0}$EXTERNAL = {1.0}PX
Reaction5 {1.0}$EXTERNAL = {1.0}PY
Reaction6 {1.0}$EXTERNAL = {1.0}PZ
Reaction7 {1.0}PX = {1.0}$EXTERNAL
Reaction8 {1.0}PY = {1.0}$EXTERNAL
Reaction9 {1.0}PZ = {1.0}$EXTERNAL
end reactions

begin rate equations
Reaction1 = kdmRNA*X[t]
Reaction10 = a0tr + (atr*KM^n)/(KM^n + PZ[t]^n)
Reaction11 = a0tr + (atr*KM^n)/(KM^n + PX[t]^n)
Reaction12 = a0tr + (atr*KM^n)/(KM^n + PY[t]^n)
Reaction2 = kdmRNA*Y[t]
Reaction3 = kdmRNA*Z[t]
Reaction4 = ktl*X[t]
Reaction5 = ktl*Y[t]
Reaction6 = ktl*Z[t]
Reaction7 = kdprot*PX[t]
Reaction8 = kdprot*PY[t]
Reaction9 = kdprot*PZ[t]
end rate equations

begin parameters
KM = 40.0
eff = 20.0
n = 2.0
ps0 = 0.0005
psa = 0.5
taumRNA = 2.0
tauprot = 10.0
cell = 1.0
end parameters

begin initial conditions
PX[0] = PXi
PY[0] = PYi
PZ[0] = PZi
X[0] = Xi
Y[0] = Yi
Z[0] = Zi
end initial conditions

begin initial values
PXi = 0.0
PYi = 0.0
PZi = 0.0
Xi = 0.0
Yi = 20.0
Zi = 0.0
end initial values

begin assignment rules
alpha0 := (a0tr*eff*tauprot)/(KM*Log[2])
alpha := (atr*eff*tauprot)/(KM*Log[2])
kdmRNA := Log[2]/taumRNA
kdprot := Log[2]/tauprot
a0tr := 60*ps0
atr := 60*(-ps0 + psa)
ktl := eff/tave
beta := taumRNA/tauprot
tave := taumRNA/Log[2]
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
PX http://identifiers.org/uniprot/P03023
PY http://identifiers.org/uniprot/P04483
PZ http://identifiers.org/uniprot/P03034
end species annotations

begin units
end units