begin name
bruggeman2
end name

begin reactions
v1 {1.0}$S = {1.0}x1
v2 {1.0}x1 = {1.0}x2
v3 {1.0}x2 = {1.0}x3
v4 {1.0}x3 = {1.0}$P
v5 {1.0}$mRsource = {1.0}mR
v6 {1.0}mR = {1.0}$mRsink
v7 {1.0}$esource = {1.0}e
v8 {1.0}e = {1.0}$esink
end reactions

begin rate equations
v1 = (k1*S*e[t]*(1 - x1[t]/(Keq1*S))*(((1 + S/K1s + x1[t]/K1x1)/(1 + (cs*S)/K1s + (cp*x1[t])/K1x1))^n + (cs*L*(1 + S/K1s + x1[t]/K1x1)*(1 + x3[t]/K1x3)^n)/(1 + (cs*S)/K1s + (cp*x1[t])/K1x1)))/(K1s*(1 + S/K1s + x1[t]/K1x1)*(((1 + S/K1s + x1[t]/K1x1)/(1 + (cs*S)/K1s + (cp*x1[t])/K1x1))^n + L*(1 + x3[t]/K1x3)^n))
v2 = (k2*e[t]*x1[t]*(1 - x2[t]/(Keq2*x1[t])))/(K2x1*(1 + x1[t]/K2x1 + x2[t]/K2x2))
v3 = (k3*e[t]*x2[t]*(1 - x3[t]/(Keq3*x2[t])))/(K3x2*(1 + x2[t]/K3x2 + x3[t]/K3x3))
v4 = (V4*(1 - P/(Keq4*x3[t]))*x3[t])/(K4x3*(1 + P/K4P + x3[t]/K4x3))
v5 = (1 + (x3[t]/Ktrx3)^ntr)^(-1)
v6 = ktrscdeg*mR[t]
v7 = ktrnlsyn*mR[t]
v8 = ktrnldeg*e[t]
end rate equations

begin parameters
K1s = 0.05
K1x1 = 1000.0
K1x3 = 1.0
K2x1 = 1.0
K2x2 = 1.0
K3x2 = 1.0
K3x3 = 10.0
K4P = 1.0
K4x3 = 0.4
Keq1 = 100.0
Keq2 = 10.0
Keq3 = 1.0
Keq4 = 10.0
Ktrx3 = 0.5
L = 0.8
V4 = 100.0
cp = 0.0001
cs = 0.0001
k1 = 100.0
k2 = 50.0
k3 = 75.0
ktrnldeg = 0.5
ktrnlsyn = 1.0
ktrscdeg = 0.5
n = 3.0
ntr = 2.0
P = 1.0
S = 1.0
esink = 1.0
esource = 1.0
mRsink = 1.0
mRsource = 1.0
defaultcompartment = 1.0
end parameters

begin initial conditions
e[0] = ei
mR[0] = mRi
x1[0] = x1i
x2[0] = x2i
x3[0] = x3i
end initial conditions

begin initial values
ei = 1.0
mRi = 1.0
x1i = 1.0
x2i = 1.0
x3i = 1.0
end initial values

begin assignment rules
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