#schaber

FIX: p s

#reactions
v_1:
 {1.0}Ste5 + {1.0}Ste11 = {1.0}Ste5Ste11
  k1 * Ste5 * Ste11
v_10:
 {1.0}Ste5Ste11GbgKss1 = {1.0}Ste5Ste11GbgKss1P
  k10*Ste5Ste11GbgKss1
v_11:
 {1.0}Ste5Ste11GbgKss1P = {1.0}Ste5Ste11GbgP + {1.0}Kss1PP
  k11*Ste5Ste11GbgKss1P
v_12:
 {1.0}Ste5Ste11GbgP + {1.0}Kss1 = {1.0}Ste5Ste11GbgKss1P
  k12*Ste5Ste11GbgP*Kss1
v_13:
 {1.0}Ste11 = {1.0}Ste11P
  k13*Ste11*beta
v_14:
 {1.0}Ste11P = {1.0}Ste11
  k14*Ste11P
v_15:
 {1.0}Kss1 = {1.0}Kss1PP
  k15*Kss1*Ste11P + k30*Kss1*Ste11Ubi
v_16:
 {1.0}Kss1PP = {1.0}Kss1
  k16*Kss1PP + k28*Kss1PP*Fus3PP
v_17:
 {1.0}Ste12Kss1 = {2.0}Kss1 + {1.0}Ste12
  k17*Ste12Kss1
v_18:
 {2.0}Kss1 + {1.0}Ste12 = {1.0}Ste12Kss1
  k18*Kss1*Ste12
v_19:
 {1.0}Ste12 = {1.0}Ste12P
  k19*Ste12*Fus3PP + k29*Ste12*Kss1PP
v_2:
 {1.0}Ste5Ste11 + {1.0}Gbg = {1.0}Ste5Ste11Gbg
  k2*Ste5Ste11*Gbg*alpha
v_20:
 {1.0}s = {1.0}PREP
  k20*Ste12P
v_21:
 {1.0}Ste12TeSte5Kss1 = {1.0}Kss1 + {1.0}Ste12TeSte5
  k21*Ste12TeSte5Kss1
v_22:
 {1.0}Kss1 + {1.0}Ste12TeSte5 = {1.0}Ste12TeSte5Kss1
  k22*Kss1*Ste12TeSte5
v_23:
 {1.0}Ste12TeSte5 = {1.0}Ste12TeSte5P
  k23*Ste12TeSte5*Kss1PP
v_24:
 {1.0}Ste12TeSte5 = {1.0}p
  k24*Ste12TeSte5*Fus3PP
v_25:
 {1.0}s = {1.0}FREP
  k25*Ste12TeSte5P
v_26:
 {1.0}Fus3PP = {1.0}Fus3
  k26*Fus3PP
v_27:
 {1.0}Ste5Ste11 = {1.0}Ste5 + {1.0}Ste11
  k27*Ste5Ste11
v_28:
 {1.0}Ste12P = {1.0}Ste12
  k31*Ste12P
v_29:
 {1.0}PREP = {1.0}p
  k32*PREP
v_3:
 {1.0}Ste5Ste11Gbg + {1.0}Fus3 = {1.0}Ste5Ste11GbgFus3
  k3*Ste5Ste11Gbg*Fus3
v_30:
 {1.0}Ste12TeSte5P = {1.0}Ste12TeSte5
  k33*Ste12TeSte5P
v_31:
 {1.0}FREP = {1.0}p
  k34*FREP
v_4:
 {1.0}Ste5Ste11GbgFus3 = {1.0}Ste5Ste11GbgFus3P
  k4*Ste5Ste11GbgFus3
v_5:
 {1.0}Ste5Ste11GbgFus3P = {1.0}Fus3PP + {1.0}Ste5Ste11GbgP
  k5*Ste5Ste11GbgFus3P
v_6:
 {1.0}Fus3 + {1.0}Ste5Ste11GbgP = {1.0}Ste5Ste11GbgFus3P
  k6*Fus3*Ste5Ste11GbgP
v_7:
 {1.0}Ste5Ste11GbgP = {1.0}Gbg + {1.0}Ste11Ubi + {1.0}Ste5
  k7*Ste5Ste11GbgP
v_8:
 {1.0}Ste11Ubi = {1.0}p
  k8*Ste11Ubi
v_9:
 {1.0}Ste5Ste11Gbg + {1.0}Kss1 = {1.0}Ste5Ste11GbgKss1
  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
default_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*(piecewise(1,(time - betat)*(betae - time)>=0,0)*betaA*(1 - exp(-(time - betat)/betas)) + piecewise(1,(time - betae)>=0,0)*betaA*exp(-(time - betae)/betas))
!F alpha = alphastim*(piecewise(1,(time - alphat)*(alphae - time)>=0,0)*alphaA*(1 - exp(-(time - alphat)/alphas)) + piecewise(1,time - alphae>=0,0)*alphaA*exp(-(time - alphae)/alphas))