begin name
weimann1
end name

begin reactions
ActiveBMAL1degradation {1.0}y7 = {1.0}$EXTERNAL
BMAL1activation {1.0}y6 = {1.0}y7
BMAL1deactivation {1.0}y7 = {1.0}y6
BMAL1nuclearexport {1.0}y6 = {1.0}y5
BMAL1nuclearimport {1.0}y5 = {1.0}y6
BMAL1translation {1.0}$EXTERNAL = {1.0}y5
Bmal1mRNAdegradation {1.0}y4 = {1.0}$EXTERNAL
Bmal1transcription {1.0}$EXTERNAL = {1.0}y4
cytoplasmicBMAL1degradation {1.0}y5 = {1.0}$EXTERNAL
cytoplasmicper2crycomplexdegradation {1.0}y2 = {1.0}$EXTERNAL
nuclearBMAL1degradation {1.0}y6 = {1.0}$EXTERNAL
nuclearper2crycomplexdegradation {1.0}y3 = {1.0}$EXTERNAL
per2crycomplexformation {1.0}$EXTERNAL = {1.0}y2
per2crymRNAdegradation {1.0}y1 = {1.0}$EXTERNAL
per2crynuclearexport {1.0}y3 = {1.0}y2
per2crynuclearimport {1.0}y2 = {1.0}y3
per2crytranscription {1.0}$EXTERNAL = {1.0}y1
end reactions

begin rate equations
ActiveBMAL1degradation = k7d*Nucleus*y7[t]
BMAL1activation = k6a*Nucleus*y6[t]
BMAL1deactivation = k7a*Nucleus*y7[t]
BMAL1nuclearexport = k6t*Nucleus*y6[t]
BMAL1nuclearimport = Cytoplasm*k5t*y5[t]
BMAL1translation = Cytoplasm*k5b*y4[t]
Bmal1mRNAdegradation = Cytoplasm*k4d*y4[t]
Bmal1transcription = Cytoplasm*transBmal1
cytoplasmicBMAL1degradation = Cytoplasm*k5d*y5[t]
cytoplasmicper2crycomplexdegradation = Cytoplasm*k2d*y2[t]
nuclearBMAL1degradation = k6d*Nucleus*y6[t]
nuclearper2crycomplexdegradation = k3d*Nucleus*y3[t]
per2crycomplexformation = Cytoplasm*k2b*y1[t]^q
per2crymRNAdegradation = Cytoplasm*k1d*y1[t]
per2crynuclearexport = k3t*Nucleus*y3[t]
per2crynuclearimport = Cytoplasm*k2t*y2[t]
per2crytranscription = Cytoplasm*transper2cry
end rate equations

begin parameters
c = 0.01
hillcoeff = 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
end parameters

begin initial conditions
y1[0] = y1i
y2[0] = y2i
y3[0] = y3i
y4[0] = y4i
y5[0] = y5i
y6[0] = y6i
y7[0] = y7i
end initial conditions

begin initial values
y1i = 0.2
y2i = 0.0
y3i = 1.1
y4i = 0.8
y5i = 1.0
y6i = 1.0
y7i = 1.05
end initial values

begin assignment rules
y5y6y7 := y5[t] + y6[t] + y7[t]
transBmal1 := (v4b*y3[t]^r)/(k4b^r + y3[t]^r)
transper2cry := (v1b*(c + y7[t]))/(c + k1b*(1 + (y3[t]/k1i)^hillcoeff) + y7[t])
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units