#weimann1



#reactions
Active_BMAL1_degradation:
 {1.0}y7 = $pool
  Nucleus * k7d * y7
BMAL1_activation:
 {1.0}y6 = {1.0}y7
  Nucleus * k6a * y6
BMAL1_deactivation:
 {1.0}y7 = {1.0}y6
  Nucleus * k7a * y7
BMAL1_nuclear_export:
 {1.0}y6 = {1.0}y5
  Nucleus * k6t * y6
BMAL1_nuclear_import:
 {1.0}y5 = {1.0}y6
  Cytoplasm * k5t * y5
BMAL1_translation:
 $pool = {1.0}y5
  Cytoplasm * k5b * y4
Bmal1_mRNA_degradation:
 {1.0}y4 = $pool
  Cytoplasm * k4d * y4
Bmal1_transcription:
 $pool = {1.0}y4
  Cytoplasm * trans_Bmal1
cytoplasmic_BMAL1_degradation:
 {1.0}y5 = $pool
  Cytoplasm * k5d * y5
cytoplasmic_per2_cry_complex_degradation:
 {1.0}y2 = $pool
  Cytoplasm * k2d * y2
nuclear_BMAL1_degradation:
 {1.0}y6 = $pool
  Nucleus * k6d * y6
nuclear_per2_cry_complex_degradation:
 {1.0}y3 = $pool
  Nucleus * k3d * y3
per2_cry_complex_formation:
 $pool = {1.0}y2
  Cytoplasm * k2b * pow(y1, q)
per2_cry_mRNA_degradation:
 {1.0}y1 = $pool
  Cytoplasm * k1d * y1
per2_cry_nuclear_export:
 {1.0}y3 = {1.0}y2
  Nucleus * k3t * y3
per2_cry_nuclear_import:
 {1.0}y2 = {1.0}y3
  Cytoplasm * k2t * y2
per2_cry_transcription:
 $pool = {1.0}y1
  Cytoplasm * trans_per2_cry


#parameters
c = 0.01
hill_coeff = 8.0
k1b = 1.0
k1d = 0.12
k1i = 0.56
k2b = 0.3
k2d = 0.05
k2t = 0.24
k3d = 0.12
k3t = 0.02
k4b = 2.16
k4d = 0.75
k5b = 0.24
k5d = 0.06
k5t = 0.45
k6a = 0.09
k6d = 0.12
k6t = 0.06
k7a = 0.003
k7d = 0.09
q = 2.0
r = 3.0
v1b = 9.0
v4b = 3.6
Cytoplasm = 1.0
Nucleus = 1.0


#initial values
y1 = 0.2
y2 = 0.0
y3 = 1.1
y4 = 0.8
y5 = 1.0
y6 = 1.0
y7 = 1.05


#assignment rules
!F y5_y6_y7 = y5 + y6 + y7
!F trans_Bmal1 = v4b * pow(y3, r) / (pow(k4b, r) + pow(y3, r))
!F trans_per2_cry = v1b * (y7 + c) / (k1b * (1.0 + pow(y3 / k1i, hill_coeff)) + y7 + c)