begin name
locke1
end name

begin reactions
four12 {1.0}cLc = {1.0}$EXTERNAL
four15 {1.0}cTc = {1.0}$EXTERNAL
four18 {1.0}cXc = {1.0}$EXTERNAL
four21 {1.0}cYc = {1.0}$EXTERNAL
one11 {1.0}$EXTERNAL = {1.0}cLm
one12 {1.0}$EXTERNAL = {1.0}cLc
one13 {1.0}cLc = {1.0}cLn
one14 {1.0}$EXTERNAL = {1.0}cTm
one15 {1.0}$EXTERNAL = {1.0}cTc
one16 {1.0}cTc = {1.0}cTn
one17 {1.0}$EXTERNAL = {1.0}cXm
one18 {1.0}$EXTERNAL = {1.0}cXc
one19 {1.0}cXc = {1.0}cXn
one20 {1.0}$EXTERNAL = {1.0}cYm
one21 {1.0}$EXTERNAL = {1.0}cYc
one22 {1.0}cYc = {1.0}cYn
one23 {1.0}$EXTERNAL = {1.0}cPn
three11 {1.0}cLm = {1.0}$EXTERNAL
three12 {1.0}cLn = {1.0}cLc
three13 {1.0}cLn = {1.0}$EXTERNAL
three15 {1.0}cTn = {1.0}cTc
three16 {1.0}cTn = {1.0}$EXTERNAL
three18 {1.0}cXn = {1.0}cXc
three19 {1.0}cXn = {1.0}$EXTERNAL
three21 {1.0}cYn = {1.0}cYc
three22 {1.0}cYn = {1.0}$EXTERNAL
three23 {1.0}cPn = {1.0}$EXTERNAL
two11 {1.0}$EXTERNAL = {1.0}cLm
two14 {1.0}cTm = {1.0}$EXTERNAL
two17 {1.0}cXm = {1.0}$EXTERNAL
two20 {1.0}cYm = {1.0}$EXTERNAL
two23 {1.0}cPn = {1.0}$EXTERNAL
end reactions

begin rate equations
four12 = (compartment*m2*cLc[t])/(k2 + cLc[t])
four15 = (compartment*cTc[t]*(m6 + m5*(1 - LD[t, dayLength])))/(k5 + cTc[t])
four18 = (compartment*m10*cXc[t])/(k8 + cXc[t])
four21 = (compartment*m13*cYc[t])/(k11 + cYc[t])
one11 = compartment*q1*cPn[t]*LD[t, dayLength]
one12 = compartment*p1*cLm[t]
one13 = compartment*r1*cLc[t]
one14 = (compartment*g3^c*n2*cYn[t]^b)/((g3^c + cLn[t]^c)*(g2^b + cYn[t]^b))
one15 = compartment*p2*cTm[t]
one16 = compartment*r3*cTc[t]
one17 = (compartment*n3*cTn[t]^d)/(g4^d + cTn[t]^d)
one18 = compartment*p3*cXm[t]
one19 = compartment*r5*cXc[t]
one20 = (compartment*g6^f*(q2*cPn[t]*LD[t, dayLength] + (g5^e*(n5 + n4*LD[t, dayLength]))/(g5^e + cTn[t]^e)))/(g6^f + cLn[t]^f)
one21 = compartment*p4*cYm[t]
one22 = compartment*r7*cYc[t]
one23 = compartment*p5*(1 - LD[t, dayLength])
three11 = (compartment*m1*cLm[t])/(k1 + cLm[t])
three12 = compartment*r2*cLn[t]
three13 = (compartment*m3*cLn[t])/(k3 + cLn[t])
three15 = compartment*r4*cTn[t]
three16 = (compartment*cTn[t]*(m8 + m7*(1 - LD[t, dayLength])))/(k6 + cTn[t])
three18 = compartment*r6*cXn[t]
three19 = (compartment*m11*cXn[t])/(k9 + cXn[t])
three21 = compartment*r8*cYn[t]
three22 = (compartment*m14*cYn[t])/(k12 + cYn[t])
three23 = compartment*q3*cPn[t]*LD[t, dayLength]
two11 = (compartment*n1*cXn[t]^a)/(g1^a + cXn[t]^a)
two14 = (compartment*m4*cTm[t])/(k4 + cTm[t])
two17 = (compartment*m9*cXm[t])/(k7 + cXm[t])
two20 = (compartment*m12*cYm[t])/(k10 + cYm[t])
two23 = (compartment*m15*cPn[t])/(k13 + cPn[t])
end rate equations

begin parameters
CP = 24.0
Lmax = 1.0
Lmin = 0.0
a = 3.3064
b = 1.0258
c = 1.0258
d = 1.4422
dayLength = 12.0
e = 3.6064
f = 1.0237
g1 = 0.876738488
g2 = 0.036805783
g3 = 0.26593318
g4 = 0.538811228
g5 = 1.17803247
g6 = 0.064455137
k1 = 1.817
k10 = 1.7303
k11 = 1.8258
k12 = 1.8066
k13 = 1.2
k2 = 1.5644
k3 = 1.2765
k4 = 2.5734
k5 = 2.7454
k6 = 0.4033
k7 = 6.5585
k8 = 0.6632
k9 = 17.1111
m1 = 1.5283
m10 = 0.2179
m11 = 3.3442
m12 = 4.297
m13 = 0.1347
m14 = 0.6114
m15 = 1.2
m2 = 20.44
m3 = 3.6888
m4 = 3.8231
m5 = 0.0013
m6 = 3.1741
m7 = 0.0492
m8 = 4.0424
m9 = 10.1132
n1 = 5.1694
n2 = 3.0087
n3 = 0.2431
n4 = 0.0857
n5 = 0.1649
p1 = 0.8295
p2 = 4.324
p3 = 2.147
p4 = 0.2485
p5 = 0.5
q1 = 2.4514
q2 = 2.40178
q3 = 1.0
r1 = 16.8363
r2 = 0.1687
r3 = 0.3166
r4 = 2.1509
r5 = 1.0352
r6 = 3.3017
r7 = 2.2123
r8 = 0.2002
two11Fch0 = 8.0
four12Fch0 = 8.0
compartment = 1.0
end parameters

begin initial conditions
cLc[0] = cLci
cLm[0] = cLmi
cLn[0] = cLni
cPn[0] = cPni
cTc[0] = cTci
cTm[0] = cTmi
cTn[0] = cTni
cXc[0] = cXci
cXm[0] = cXmi
cXn[0] = cXni
cYc[0] = cYci
cYm[0] = cYmi
cYn[0] = cYni
end initial conditions

begin initial values
cLci = 0.015
cLmi = 0.539
cLni = 0.0855
cPni = 0.835
cTci = 8.7
cTmi = 0.456
cTni = 0.4
cXci = 1.34
cXmi = 0.08
cXni = 0.4
cYci = 0.004
cYmi = 0.037
cYni = 0.013
end initial values

begin assignment rules
end assignment rules

begin function definitions
LD[tod_,length_] := Ceiling[Sin[0.001 + (Pi*tod)/length]/2]
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
cLc http://identifiers.org/uniprot/O81713
cLn http://identifiers.org/uniprot/O81713
cTc http://identifiers.org/uniprot/Q9LKL2
cTn http://identifiers.org/uniprot/Q9LKL2
end species annotations

begin units
end units