#bruggeman2

FIX: P S esink esource mRsink mRsource

#reactions
v_1:
 {1.0}S = {1.0}x1
  (e*k1*S*(1 - x1/(Keq1*S))*(((1 + S/K1s + x1/K1x1)/(1 + (cs*S)/K1s + (cp*x1)/K1x1))**n + (cs*L*(1 + S/K1s + x1/K1x1)*(1 + x3/K1x3)**n)/(1 + (cs*S)/K1s + (cp*x1)/K1x1)))/(K1s*(1 + S/K1s + x1/K1x1)*(((1 + S/K1s + x1/K1x1)/(1 + (cs*S)/K1s + (cp*x1)/K1x1))**n + L*(1 + x3/K1x3)**n))
v_2:
 {1.0}x1 = {1.0}x2
  (e*k2*x1*(1 - x2/(Keq2*x1)))/(K2x1*(1 + x1/K2x1 + x2/K2x2))
v_3:
 {1.0}x2 = {1.0}x3
  (e*k3*x2*(1 - x3/(Keq3*x2)))/(K3x2*(1 + x2/K3x2 + x3/K3x3))
v_4:
 {1.0}x3 = {1.0}P
  (V4*(1 - P/(Keq4*x3))*x3)/(K4x3*(1 + P/K4P + x3/K4x3))
v_5:
 {1.0}mRsource = {1.0}mR
  (1 + (x3/Ktrx3)**ntr)**(-1)
v_6:
 {1.0}mR = {1.0}mRsink
  ktrscdeg*mR
v_7:
 {1.0}esource = {1.0}e
  ktrnlsyn*mR
v_8:
 {1.0}e = {1.0}esink
  e*ktrnldeg


#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
default_compartment = 1.0


#initial values
e = 1.0
mR = 1.0
x1 = 1.0
x2 = 1.0
x3 = 1.0


#assignment rules