begin name
leroux
end name

begin reactions
v1 {1.0}e + {1.0}$A = {1.0}eA
v10 {1.0}eAQ + {1.0}$B = {1.0}eABQ
v11 {1.0}eQ + {1.0}$B = {1.0}eBQ
v12 {1.0}eBQ + {1.0}$A = {1.0}eABQ
v13 {1.0}eABQ = {1.0}P + {1.0}XQ
v14 {1.0}XQ = {1.0}eR
v15 {1.0}eR = {1.0}e + {1.0}$R
v16 {1.0}XQ + {1.0}$B = {1.0}XQB
v17 {1.0}XQ + {1.0}$R = {1.0}XQR
v19 {1.0}e + {1.0}$C = {1.0}eC
v2 {1.0}eA + {1.0}$B = {1.0}eAB
v20 {1.0}eC + {1.0}$A = {1.0}eAC
v21 {1.0}eAC + {1.0}$B = {1.0}eABC
v22 {1.0}eC + {1.0}$B = {1.0}eBC
v23 {1.0}eBC + {1.0}$A = {1.0}eABC
v24 {1.0}eA + {1.0}$C = {1.0}eAC
v25 {1.0}eB + {1.0}$C = {1.0}eBC
v26 {1.0}eAB + {1.0}$C = {1.0}eABC
v27 {1.0}eABC = {1.0}P + {1.0}XC
v28 {1.0}XC = {1.0}eeQ
v29 {1.0}eeQ = {1.0}e + {1.0}Q
v3 {1.0}e + {1.0}$B = {1.0}eB
v30 {1.0}XC + {1.0}$B = {1.0}XCB
v31 {1.0}XC + {1.0}$R = {1.0}XCR
v32 {1.0}P = {1.0}$S
v4 {1.0}eB + {1.0}$A = {1.0}eAB
v5 {1.0}eAB + {1.0}Q = {1.0}eABQ
v6 {1.0}eA + {1.0}Q = {1.0}eAQ
v7 {1.0}eB + {1.0}Q = {1.0}eBQ
v8 {1.0}e + {1.0}Q = {1.0}eQ
v9 {1.0}eQ + {1.0}$A = {1.0}eAQ
end reactions

begin rate equations
v1 = (A*k1*e[t] - k1a*eA[t])/Cell
v10 = (-(k10a*eABQ[t]) + B*k10*eAQ[t])/Cell
v11 = (-(k11a*eBQ[t]) + B*k11*eQ[t])/Cell
v12 = (-(k12a*eABQ[t]) + A*k12*eBQ[t])/Cell
v13 = (k13*eABQ[t])/Cell
v14 = (-(k14a*eR[t]) + k14*XQ[t])/Cell
v15 = (-(k15a*R*e[t]) + k15*eR[t])/Cell
v16 = (B*k16*XQ[t] - k16a*XQB[t])/Cell
v17 = (k17*R*XQ[t] - k17a*XQR[t])/Cell
v19 = (C*k19*e[t] - k19a*eC[t])/Cell
v2 = (B*k2*eA[t] - k2a*eAB[t])/Cell
v20 = (-(k20a*eAC[t]) + A*k20*eC[t])/Cell
v21 = (-(k21a*eABC[t]) + B*k21*eAC[t])/Cell
v22 = (-(k22a*eBC[t]) + B*k22*eC[t])/Cell
v23 = (-(k23a*eABC[t]) + A*k23*eBC[t])/Cell
v24 = (C*k24*eA[t] - k24a*eAC[t])/Cell
v25 = (C*k25*eB[t] - k25a*eBC[t])/Cell
v26 = (C*k26*eAB[t] - k26a*eABC[t])/Cell
v27 = (k27*eABC[t])/Cell
v28 = (-(k28a*eeQ[t]) + k28*XC[t])/Cell
v29 = (k29*eeQ[t] - k29a*e[t]*Q[t])/Cell
v3 = (B*k3*e[t] - k3a*eB[t])/Cell
v30 = (B*k30*XC[t] - k30a*XCB[t])/Cell
v31 = (k31*R*XC[t] - k31a*XCR[t])/Cell
v32 = (k32*P[t])/Cell
v4 = (-(k4a*eAB[t]) + A*k4*eB[t])/Cell
v5 = (-(k5a*eABQ[t]) + k5*eAB[t]*Q[t])/Cell
v6 = (-(k6a*eAQ[t]) + k6*eA[t]*Q[t])/Cell
v7 = (-(k7a*eBQ[t]) + k7*eB[t]*Q[t])/Cell
v8 = (-(k8a*eQ[t]) + k8*e[t]*Q[t])/Cell
v9 = (-(k9a*eAQ[t]) + A*k9*eQ[t])/Cell
end rate equations

begin parameters
Cell = 1.0
eXTeRNAL = 0.0
k1 = 13.4953706168796
k10 = 261.19
k10a = 7289.2804823
k11 = 6210.12345627707
k11a = 315.68381975
k12 = 29.9958100432
k12a = 1330.84779006
k13 = 2.27895805686
k14 = 543.476605024
k14a = "1e-10"
k15 = 34.9267638563
k15a = 4941.84508678
k16 = 488.901949443
k16a = 2199.40597503
k17 = 48.5266386368
k17a = 82.5043501462
k19 = 0.01739668189544
k19a = 868.323556017581
k1a = 134.2939838
k2 = 6317.116733
k20 = "2.005267865347e-06"
k20a = "1e-10"
k21 = 94.6031417823562
k21a = 56632.7173679014
k22 = 897.73729266064
k22a = 26.8002418488192
k23 = "1e-10"
k23a = "1e-10"
k24 = 54684.4190465446
k24a = 13678.3788270169
k25 = "8.61285612366646e-14"
k25a = "1e-10"
k26 = 0.0474957203508858
k26a = 1058.42481633921
k27 = 15.1770330320787
k28 = 455.465653318972
k28a = 290.422551587577
k29 = 30247.7736950703
k29a = 200.14320986275
k2a = 42.4469248983
k3 = 223.401281347
k30 = 234.236580254189
k30a = 1218.73611055286
k31 = "1e-10"
k31a = 420.551496677673
k32 = 14253.4745502918
k3a = 286.7064474
k4 = 762.2665647
k4a = 39.71496125
k5 = 22.91729291
k5a = 1465.04965
k6 = "6.4969834151108e-09"
k6a = "1e-10"
k7 = 7067.22002146676
k7a = 530.538761235
k8 = 125.296650629147
k8a = 237.469413308
k9 = 369.688932945
k9a = 29.8763965658
A = 2300.0
B = 300.0
C = 8000.0
R = 0.0
S = 0.0
defaultcompartment = 1.0
end parameters

begin initial conditions
P[0] = Pi
Q[0] = Qi
XC[0] = XCi
XCB[0] = XCBi
XCR[0] = XCRi
XQ[0] = XQi
XQB[0] = XQBi
XQR[0] = XQRi
e[0] = ei
eA[0] = eAi
eAB[0] = eABi
eABC[0] = eABCi
eABQ[0] = eABQi
eAC[0] = eACi
eAQ[0] = eAQi
eB[0] = eBi
eBC[0] = eBCi
eBQ[0] = eBQi
eC[0] = eCi
eQ[0] = eQi
eR[0] = eRi
eeQ[0] = eeQi
end initial conditions

begin initial values
Pi = 0.0
Qi = 0.0
XCi = 0.0
XCBi = 0.0
XCRi = 0.0
XQi = 0.0
XQBi = 0.0
XQRi = 0.0
ei = 1.0
eAi = 0.0
eABi = 0.0
eABCi = 0.0
eABQi = 0.0
eACi = 0.0
eAQi = 0.0
eBi = 0.0
eBCi = 0.0
eBQi = 0.0
eCi = 0.0
eQi = 0.0
eRi = 0.0
eeQi = 0.0
end initial values

begin assignment rules
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