#schaber1

FIX: p s

#reactions
v1:
 {1.0}Ste5 + {1.0}Ste11 = {1.0}Ste5Ste11
  compartment * k1 * Ste5 * Ste11
v10:
 {1.0}Ste5Ste11GbgKss1 = {1.0}Ste5Ste11GbgKss1P
  compartment * k10 * Ste5Ste11GbgKss1
v11:
 {1.0}Ste5Ste11GbgKss1P = {1.0}Ste5Ste11GbgP + {1.0}Kss1PP
  compartment * k11 * Ste5Ste11GbgKss1P
v12:
 {1.0}Ste5Ste11GbgP + {1.0}Kss1 = {1.0}Ste5Ste11GbgKss1P
  compartment * k12 * Ste5Ste11GbgP * Kss1
v13:
 {1.0}Ste11 = {1.0}Ste11P
  compartment * k13 * Ste11 * beta
v14:
 {1.0}Ste11P = {1.0}Ste11
  compartment * k14 * Ste11P
v15:
 {1.0}Kss1 = {1.0}Kss1PP
  compartment * (k15 * Kss1 * Ste11P + k30 * Kss1 * Ste11Ubi)
v16:
 {1.0}Kss1PP = {1.0}Kss1
  compartment * (k16 * Kss1PP + k28 * Kss1PP * Fus3PP)
v17:
 {1.0}Ste12Kss1 = {2.0}Kss1 + {1.0}Ste12
  compartment * k17 * Ste12Kss1
v18:
 {2.0}Kss1 + {1.0}Ste12 = {1.0}Ste12Kss1
  compartment * k18 * Kss1 * Ste12
v19:
 {1.0}Ste12 = {1.0}Ste12P
  compartment * (k19 * Ste12 * Fus3PP + k29 * Ste12 * Kss1PP)
v2:
 {1.0}Ste5Ste11 + {1.0}Gbg = {1.0}Ste5Ste11Gbg
  compartment * k2 * Ste5Ste11 * Gbg * alpha
v20:
 {1.0}s = {1.0}PREP
  compartment * k20 * Ste12P
v21:
 {1.0}Ste12TeSte5Kss1 = {1.0}Kss1 + {1.0}Ste12TeSte5
  compartment * k21 * Ste12TeSte5Kss1
v22:
 {1.0}Kss1 + {1.0}Ste12TeSte5 = {1.0}Ste12TeSte5Kss1
  compartment * k22 * Kss1 * Ste12TeSte5
v23:
 {1.0}Ste12TeSte5 = {1.0}Ste12TeSte5P
  compartment * k23 * Ste12TeSte5 * Kss1PP
v24:
 {1.0}Ste12TeSte5 = {1.0}p
  compartment * k24 * Ste12TeSte5 * Fus3PP
v25:
 {1.0}s = {1.0}FREP
  compartment * k25 * Ste12TeSte5P
v26:
 {1.0}Fus3PP = {1.0}Fus3
  compartment * k26 * Fus3PP
v27:
 {1.0}Ste5Ste11 = {1.0}Ste5 + {1.0}Ste11
  compartment * k27 * Ste5Ste11
v28:
 {1.0}Ste12P = {1.0}Ste12
  compartment * k31 * Ste12P
v29:
 {1.0}PREP = {1.0}p
  compartment * k32 * PREP
v3:
 {1.0}Ste5Ste11Gbg + {1.0}Fus3 = {1.0}Ste5Ste11GbgFus3
  compartment * k3 * Ste5Ste11Gbg * Fus3
v30:
 {1.0}Ste12TeSte5P = {1.0}Ste12TeSte5
  compartment * k33 * Ste12TeSte5P
v31:
 {1.0}FREP = {1.0}p
  compartment * k34 * FREP
v4:
 {1.0}Ste5Ste11GbgFus3 = {1.0}Ste5Ste11GbgFus3P
  compartment * k4 * Ste5Ste11GbgFus3
v5:
 {1.0}Ste5Ste11GbgFus3P = {1.0}Fus3PP + {1.0}Ste5Ste11GbgP
  compartment * k5 * Ste5Ste11GbgFus3P
v6:
 {1.0}Fus3 + {1.0}Ste5Ste11GbgP = {1.0}Ste5Ste11GbgFus3P
  compartment * k6 * Fus3 * Ste5Ste11GbgP
v7:
 {1.0}Ste5 + {1.0}Ste5Ste11GbgP = {1.0}Gbg + {1.0}Ste11Ubi
  compartment * k7 * Ste5Ste11GbgP
v8:
 {1.0}Ste11Ubi = {1.0}p
  compartment * k8 * Ste11Ubi
v9:
 {1.0}Ste5Ste11Gbg + {1.0}Kss1 = {1.0}Ste5Ste11GbgKss1
  compartment * k9 * Ste5Ste11Gbg * Kss1


#parameters
alphaA = 1.0
alphae = 10.0
alphas = 2.0
alphastim = 1.0
alphat = 0.0
betaA = 1.0
betae = 360.0
betas = 20.0
betastim = 1.0
betat = 0.0
k1 = 0.01
k10 = 1.0
k11 = 1.0
k12 = 1.0
k13 = 1.0
k14 = 0.1
k15 = 0.1
k16 = 0.1
k17 = 1.0
k18 = 10.0
k19 = 1.0
k2 = 0.01
k20 = 1.0
k21 = 1.0
k22 = 1.0
k23 = 1.0
k24 = 0.01
k25 = 1.0
k26 = 0.1
k27 = 1.0
k28 = 0.01
k29 = 0.01
k3 = 1.0
k30 = 0.1
k31 = 1.0
k32 = 1.0
k33 = 1.0
k34 = 1.0
k4 = 1.0
k5 = 1.0
k6 = 1.0
k7 = 10.0
k8 = 0.1
k9 = 1.0
p = 0.0
s = 0.0
compartment = 1.0


#initial values
FREP = 0.0
Fus3 = 217.0
Fus3PP = 0.0
Gbg = 53.0
Kss1 = 54.4
Kss1PP = 0.0
PREP = 0.0
Ste11 = 13.3
Ste11P = 0.0
Ste11Ubi = 0.0
Ste12 = 0.07
Ste12Kss1 = 35.9
Ste12P = 0.0
Ste12TeSte5 = 0.25
Ste12TeSte5Kss1 = 13.7
Ste12TeSte5P = 0.0
Ste5 = 42.3
Ste5Ste11 = 5.6
Ste5Ste11Gbg = 0.0
Ste5Ste11GbgFus3 = 0.0
Ste5Ste11GbgFus3P = 0.0
Ste5Ste11GbgKss1 = 0.0
Ste5Ste11GbgKss1P = 0.0
Ste5Ste11GbgP = 0.0


#assignment rules
!F beta = betastim * betaA * piecewise(1.0 - exp(-(time - betat) / betas), and(geq(time, betat), leq(time, betae)), piecewise(exp(-(time - betae) / betas), gt(time, betae), 0.0))
!F alpha = alphastim * piecewise(alphaA * (1.0 - exp(-(time - alphat) / alphas)), and(geq(time, alphat), leq(time, alphae)), piecewise(alphaA * exp(-(time - alphat) / alphas), geq(time, alphae), 0.0))