#akman1



#reactions
E1Fdeg:
 {1.0}E1F = $pool
  gam1 * E1F
E1Fpdeg:
 {1.0}E1Fp = $pool
  gam1p * E1Fp
E1Fptrl:
 $pool = {1.0}E1Fp
  a3p * MF
E1Ftrl:
 $pool = {1.0}E1F
  a3 * MF
E1Wdeg:
 {1.0}E1W = $pool
  gam2 * E1W
E1Wtrl:
 $pool = {1.0}E1W
  (a6 + a7 * (PF + PFp)) * MW
E2Fdeg:
 {1.0}E2F = $pool
  gam1 * E2F
E2Fpdeg:
 {1.0}E2Fp = $pool
  gam1p * E2Fp
E2Fptrl:
 {1.0}E1Fp = {1.0}E2Fp
  f1p * E1Fp
E2Ftrl:
 {1.0}E1F = {1.0}E2F
  f1 * E1F
E2Wdeg:
 {1.0}E2W = $pool
  gam2 * E2W
E2Wtrl:
 {1.0}E1W = {1.0}E2W
  f2 * E1W
MFdeg:
 {1.0}MF = $pool
  d1 * MF / (MF + b5)
MFtrn:
 $pool = {1.0}MF
  a1 * pow(PWL, n) / ((1 + (PF + PFp) / b1) * (pow(PWL, n) + pow(b2, n))) + a2 * pow(PW, m) / ((1 + (PF + PFp) / b3) * (pow(PW, m) + pow(b4, m)))
MWdeg:
 {1.0}MW = $pool
  d3 * MW / (MW + b8)
MWtrn:
 $pool = {1.0}MW
  a4 + a5 * pow(PWL, k) / (pow(PWL, k) + pow(b7, k))
PFdeg:
 {1.0}PF = $pool
  d2 * PF
PFpdeg:
 {1.0}PFp = $pool
  d2p * PFp
PFptrl:
 {1.0}E2Fp = {1.0}PFp
  f1p * E2Fp
PFtrl:
 {1.0}E2F = {1.0}PF
  f1 * E2F
PWLdeg:
 {1.0}PWL = $pool
  d5 * PWL / (PWL + b10)
PWdeg:
 {1.0}PW = $pool
  d4 * PW / (PW + b9)
PWtrl:
 {1.0}E2W = {1.0}PW
  f2 * E2W
PWtrs:
 {1.0}PW = {1.0}PWL
  r1 * amp * PW * (1 + tanh(24 * (time - 24 * floor(time / 24) - dawn))) * (1 - tanh(24 * (time - 24 * floor(time / 24) - dusk))) / 4 - r2 * PWL


#parameters
a1 = 24.9736
a2 = 3.59797
a3 = 0.2834
a3p = 0.34578
a4 = 0.46227
a5 = 0.02917
a6 = 0.20695
a7 = 3.02856
amp = 0.0
b1 = 0.00209
b10 = 93.03664
b2 = 2.13476
b3 = 0.08039
b4 = 0.45859
b5 = 0.13056
b6 = 0.0
b7 = 2.96739
b8 = 0.11167
b9 = 81.10381
d1 = 1.43549
d2 = 0.21251
d2p = 0.18191
d3 = 0.50309
d4 = 3.36641
d5 = 0.41085
dawn = 12.0
dusk = 24.0
f1 = 0.09292
f1p = 0.09588
f2 = 0.14979
gam1 = 0.34603
gam1p = 0.40119
gam2 = 0.00028
k = 2.18234
m = 1.34979
n = 1.02419
r1 = 2.71574
r2 = 35.40005
cytosol = 1.0
nucleus = 1.0


#initial values
E1F = 0.43076
E1Fp = 0.45583
E1W = 5.84748
E2F = 0.10647
E2Fp = 0.09872
E2W = 5.70265
MF = 0.6935
MW = 1.2689
PF = 0.06565
PFp = 0.07719
PW = 26.4393
PWL = 0.0


#assignment rules
!F WC1_tot = PW + PWL + E1W + E2W
!F Frq_tot = sFrq_tot + lFrq_tot
!F lFrq_tot = PFp + E2Fp + E1Fp
!F sFrq_tot = PF + E2F + E1F