begin name
demin1
end name

begin reactions
v1 {1.0}$EXTERNAL = {1.0}AAb
v10in {1.0}LTC4b = {1.0}$EXTERNAL
v10out {1.0}$EXTERNAL = {1.0}LTC4bout
v11 {1.0}LTC4bout = {1.0}LTD4b
v12 {1.0}LTD4b = {1.0}LTE4b
v13 {1.0}LTC4bout = {1.0}$EXTERNAL
v14 {1.0}LTD4b = {1.0}$EXTERNAL
v15 {1.0}LTE4b = {1.0}$EXTERNAL
v16 {1.0}EOb = {1.0}EOib
v17 {1.0}EOib = {1.0}EOab
v18 {1.0}EOab = {1.0}EOb
v19 {1.0}EOb = {1.0}$EXTERNAL
v2 {1.0}AAb = {1.0}$EXTERNAL
v20 {1.0}EOab = {1.0}$EXTERNAL
v21 {1.0}EOib = {1.0}EOiaw
v22 {1.0}EOab = {1.0}EOaaw
v23 {1.0}EOb = {1.0}EOaw
v24 {1.0}$EXTERNAL = {1.0}Hnb
v25 {1.0}Hnb = {1.0}$EXTERNAL
v26 {1.0}$EXTERNAL = {1.0}ILb
v27 {1.0}ILb = {1.0}$EXTERNAL
v28 {1.0}ILb = {1.0}ILbm
v29 {1.0}$EXTERNAL = {1.0}EObm
v3 {1.0}$EXTERNAL = {1.0}HPETEb
v30 {1.0}EObm = {1.0}EOb
v31 {1.0}$EXTERNAL = {1.0}AAaw
v32 {1.0}AAaw = {1.0}$EXTERNAL
v33 {1.0}$EXTERNAL = {1.0}HPETEaw
v34 {1.0}HPETEaw = {1.0}HETEaw
v35 {1.0}HPETEaw = {1.0}HETEaw
v36 {1.0}$EXTERNAL = {1.0}LTA4aw
v37 {1.0}LTA4aw = {1.0}LTC4aw
v38 {1.0}HETEaw = {1.0}$EXTERNAL
v39 {1.0}LTA4aw = {1.0}$EXTERNAL
v4 {1.0}HPETEb = {1.0}HETEb
v40in {1.0}LTC4aw = {1.0}$EXTERNAL
v40out {1.0}$EXTERNAL = {1.0}LTC4awout
v41 {1.0}LTC4awout = {1.0}LTD4aw
v42 {1.0}LTD4aw = {1.0}LTE4aw
v43 {1.0}LTE4aw = {1.0}LTE4b
v44 {1.0}LTD4aw = {1.0}LTD4b
v45 {1.0}LTC4awout = {1.0}LTC4bout
v46 {1.0}EOaw = {1.0}EOiaw
v47 {1.0}EOiaw = {1.0}EOaaw
v48 {1.0}EOaaw = {1.0}EOaw
v49 {1.0}EOaw = {1.0}$EXTERNAL
v5 {1.0}HPETEb = {1.0}HETEb
v50 {1.0}EOaaw = {1.0}$EXTERNAL
v51 {1.0}$EXTERNAL = {1.0}Hnaw
v52 {1.0}Hnaw = {1.0}Hnb
v53 {1.0}$EXTERNAL = {1.0}ILaw
v54 {1.0}ILaw = {1.0}ILb
v55 {1.0}LTE4aw = {1.0}$EXTERNAL
v56 {1.0}LTD4aw = {1.0}$EXTERNAL
v57 {1.0}LTC4awout = {1.0}$EXTERNAL
v58 {1.0}ILaw = {1.0}$EXTERNAL
v59 {1.0}Hnaw = {1.0}$EXTERNAL
v6 {1.0}$EXTERNAL = {1.0}LTA4b
v60 {1.0}ZFintes = {1.0}ZFblood
v61 {1.0}ZFblood = {1.0}ZFairways
v62 {1.0}ZFblood = {1.0}$EXTERNAL
v63 {1.0}MLintes = {1.0}MLblood
v64 {1.0}MLblood = {1.0}$EXTERNAL
v7 {1.0}LTA4b = {1.0}LTC4b
v8 {1.0}HETEb = {1.0}$EXTERNAL
v9 {1.0}LTA4b = {1.0}$EXTERNAL
end reactions

begin rate equations
v1 = Default*((PL*PLA2Ca*VmaxPLA2)/(KmPLA2APC + PL) - (VCoA*AAb[t])/(KmCoAAA + AAb[t]))
v10in = Default*Kd50*LTC4b[t]
v10out = 1000000*Default*Kd50*VLTCCB*LTC4b[t]
v11 = fupLT*kggt*VdLTC*LTC4bout[t]
v12 = fupLT*kdp*VdLTD*LTD4b[t]
v13 = fupLT*kltcltdel*VdLTC*LTC4bout[t]
v14 = fupLT*kltcltdel*VdLTD*LTD4b[t]
v15 = fupLT*(kacet + klteel)*VdLTE*LTE4b[t]
v16 = (ca*ka*OLb^hact*VB*EOb[t])/(EC50act^hact + OLb^hact)
v17 = ca*kEOm*VB*EOib[t]
v18 = ca*kia*VB*EOab[t]
v19 = kEOd*VB*EOb[t]
v2 = (Default*FLO3tb*klo*AAb[t])/KAA
v20 = kEOad*VB*EOab[t]
v21 = (kEOtbaw*Recoccupmigr^hmigr*VB*EOib[t])/(EC50migr^hmigr + Recoccupmigr^hmigr)
v22 = (kEOtbaw*Recoccupmigr^hmigr*VB*EOab[t])/(EC50migr^hmigr + Recoccupmigr^hmigr)
v23 = (kEOtbaw*Recoccupmigr^hmigr*VB*EOb[t])/(EC50migr^hmigr + Recoccupmigr^hmigr)
v24 = kHnp*VB*(EOab[t] + EOb[t] + EOib[t])
v25 = fupHn*kHnd*VdHn*Hnb[t]
v26 = kILp*VB*EOab[t]
v27 = kILd*VdIL5*ILb[t]
v28 = kILtbbm*(ILb[t] - ILbm[t]) - JBMlymfl*ILbm[t]
v29 = VBM*(k1min + (k1*ILbm[t]^hmatur)/(Km1^hmatur + ILbm[t]^hmatur))
v3 = Default*(1 + AAb[t]/KiAA)*(FLO5HPb*k3 - FLO3tb*k3*HPETEb[t])
v30 = ca*kEOtbmb*VBM*EObm[t]
v31 = Default*((PL*PLA2Ca*VmaxPLA2)/(KmPLA2APC + PL) - (VCoA*AAaw[t])/(KmCoAAA + AAaw[t]))
v32 = (Default*FLO3taw*klo*AAaw[t])/KAA
v33 = Default*(1 + AAaw[t]/KiAA)*(FLO5HPaw*k3 - FLO3taw*k3*HPETEaw[t])
v34 = Default*r1*(-((FLO3aw*HETEaw[t])/Keox) + FLO2aw*HPETEaw[t])
v35 = (Baw*Default*GPx)/Aaw
v36 = Default*FLO5HPaw*kltasyn
v37 = (Default*nomLTCsaw)/denLTCsaw
v38 = (Bhedhaw*Default*HEDH5)/Ahedhaw
v39 = Default*Kd12*LTA4aw[t]
v4 = Default*r1*(-((FLO3b*HETEb[t])/Keox) + FLO2b*HPETEb[t])
v40in = Default*Kd50*LTC4aw[t]
v40out = 1000000*Default*Kd50*VLTCCAW*LTC4aw[t]
v41 = kggt*VdAWLTC*LTC4awout[t]
v42 = kdp*VdAWLTD*LTD4aw[t]
v43 = QAWblf*RLTEB*((RLTEAW*LTE4aw[t])/KpLTEAW - LTE4b[t])
v44 = QAWblf*RLTDB*((RLTDAW*LTD4aw[t])/KpLTDAW - LTD4b[t])
v45 = QAWblf*RLTCB*((RLTCAW*LTC4awout[t])/KpLTCAW - LTC4bout[t])
v46 = (ca*ka*OLaw^hact*VAW*EOaw[t])/(EC50act^hact + OLaw^hact)
v47 = ca*kEOm*VAW*EOiaw[t]
v48 = ca*kia*VAW*EOaaw[t]
v49 = kEOd*VAW*EOaw[t]
v5 = (Bb*Default*GPx)/Ab
v50 = kEOad*VAW*EOaaw[t]
v51 = kHnp*VAW*(EOaw[t] + EOaaw[t] + EOiaw[t])
v52 = QAWblf*RHnB*((RHnAW*Hnaw[t])/KpHnAW - Hnb[t])
v53 = kILp*VAW*EOaaw[t]
v54 = JAWlymfl*ILaw[t] + kILtawb*(ILaw[t] - ILb[t])
v55 = (kacet + klteel)*VdAWLTE*LTE4aw[t]
v56 = kltcltdel*VdAWLTD*LTD4aw[t]
v57 = kltcltdel*VdAWLTC*LTC4awout[t]
v58 = kILd*VAW*ILaw[t]
v59 = kHnd*VdAWHn*Hnaw[t]
v6 = Default*FLO5HPb*kltasyn
v60 = Default*kabszf*((1000*DOSEzf*(1 - a + a*ftzf)*Fzf*oral)/MZF + ZFintes[t])
v61 = QAWblf*RZFB*(-((RZFAW*ZFairways[t])/KpZFAW) + ZFblood[t])
v62 = kelimzf*VdZF*ZFblood[t]
v63 = Default*kabsml*((1000000000*DOSEml*(1 - a + a*ftml)*Fml*oral)/MML + MLintes[t])
v64 = kelimml*VdML*MLblood[t]
v7 = (Default*nomLTCsb)/denLTCsb
v8 = (Bhedhb*Default*HEDH5)/Ahedhb
v9 = Default*Kd12*LTA4b[t]
end rate equations

begin parameters
Ca = 1.0
CaFEVex = 10000.0
CaoFEV = 10000.0
DOSEml = 0.0
DOSEzf = 0.0
EC50MLFEV = 500.0
EC50act = 0.75
EC50migr = 0.115
EtLTCs = 1.0
FLOa = 1.5
Fml = 0.660688
Fzf = 0.082
GPx = 1.6
GSHaw = 1000.0
GSHb = 5000.0
GSpoolaw = 10000.0
GSpoolb = 10000.0
HEDH5 = 0.5
JAWlymfl = 0.00115
JBMlymfl = 0.00049
KAA = 10.74959
KCa2 = 14.36738
KCa3 = 7116.527
KCaFEV = 150000.0
KGSH = 744.9176
KLTA = 1.75951
KLTAGSH = 1696.6
KLTC = 0.1951215
KPLA2Ca = 0.1
Kd12 = 0.007
Kd50 = 0.43
KdZ = 20.0
KdHnFEV = 6300000.0
KdILmigr = 50.0
KdLT = 1000.0
KdLTEmigr = 50000.0
KdLTR1FEV = 1000.0
KdLTR2FEV = 10000.0
KdLT2 = 10000.0
KdgpxGSSG = 0.07154222
KdgpxHETE5 = 6.043446
KdhedhHETE5 = 0.331696
KdhedhNADP = 2.895899
KdhedhNADPH = 2.685581
KdhedhoxoETE5 = 1.667019
Keox = 99.99979
Kered = "5.761955e-07"
KiAA = 551.8748
KiHETE = 0.5408177
KiMLEOa = 500.0
Km1 = 2.0
KmCoAAA = 0.005
KmPLA2APC = 20.0
KmgpxGSH = 600.0
KmgpxHPETE5 = 5.974381
KpHnAW = 3950.0
KpLTCAW = 0.22
KpLTDAW = 0.22
KpLTEAW = 0.22
KpZFAW = 0.204
KpZFIW = 0.631
LOHaw = 0.0
LOHb = 0.0
LOOHaw = 20.0
LOOHb = 1.0
MAXFEV = 4.94
MLairways = 0.0
MLex = 0.0
MHn = 111.2
MIL5 = 45000.0
MLTC = 625.8
MLTD = 496.7
MLTE = 439.6
MML = 586.18
MZF = 236.0
NADPaw = 2900.0
NADPb = 2000.0
NPpoolaw = 3000.0
NPpoolb = 3000.0
NApmole = 602000000000.0
OLbex = 0.0
PL = 110.0
QAWblf = 5.23
R1portionEOa = 0.39
R1portionFEV = 0.016
R2portionFEV = 0.037
RFEV = 3.007
RHnAW = 5130.0
RHnB = 141.0
RLTCAW = 1.4
RLTCB = 0.538
RLTDAW = 1.4
RLTDB = 0.538
RLTEAW = 1.4
RLTEB = 0.538
RZFAW = 2.96
RZFB = 0.533
RinrelaxFEV = 2.073
T = 1440.0
TSN0 = 15.5
VCoA = 350.0
VmaxPLA2 = 450.0
a = 1.0
al = 1.0
ca = 10.0
diamEO = 0.00012
fupHn = 0.77
fupLT = 0.16
fupML = 0.004
fupZF = 0.069
hact = 3.0
hmatur = 1.0
hmigr = 3.0
k1 = "1e-06"
k1min = "1.6e-07"
k3 = 34.0
k3 = 263640.0
kEOad = 0.00015
kEOd = 0.0003
kEOm = 10.0
kEOtbaw = 0.04
kEOtbmb = 0.02
kHnd = 0.033
kHnp = 18000000000.0
kILd = 0.0046
kILp = 16.0
kILtawb = 0.05
kILtbbm = 0.001
kLTCsback = 1003.709
kLTCsfow = 1068016.0
kabsml = 0.012
kabszf = 0.018
kacet = 0.002703885
kdp = 0.067
kelimml = 0.00225
kelimzf = 0.004
kfeveff = 3000000.0
kggt = 0.1
kgpxcat = 0.4884995
khedh1 = 88.33768
khedh2 = 1724.404
khedh3 = 31.49749
khedh4 = 8.077531
klo = 4642.68
kltasyn = 54420.0
kltcltdel = 0.1
klteel = 0.04
kox = 0.0002658
kox2 = 67.2
kred = 0.0002538
kred2 = "4.428e-05"
ka = 500.0
ki = 25000.0
kia = 0.001
nFEV = 1.0
naEOLTCsyn = 0.0
npi = 3.14
oral = 1.0
oxoETEaw = 0.0
oxoETEb = 0.0
portionmigr = 0.3
zfinh = 0.0
Default = 1.0
VAW = 0.209999993443489
VB = 2.80999994277954
VBM = 0.824999988079071
VdAWHn = 5640.0
VdAWLTC = 1.53999996185303
VdAWLTD = 1.53999996185303
VdAWLTE = 1.53999996185303
VdAWZF = 3.25999999046326
VdHn = 78100.0
VdIL5 = 10.1999998092651
VdLTC = 9.47999954223633
VdLTD = 9.47999954223633
VdLTE = 9.47999954223633
VdML = 10.0
VdZF = 9.44999980926514
end parameters

begin initial conditions
AAaw[0] = AAawi
AAb[0] = AAbi
EOaaw[0] = EOaawi
EOab[0] = EOabi
EOaw[0] = EOawi
EOb[0] = EObi
EObm[0] = EObmi
EOiaw[0] = EOiawi
EOib[0] = EOibi
HETEaw[0] = HETEawi
HETEb[0] = HETEbi
HPETEaw[0] = HPETEawi
HPETEb[0] = HPETEbi
Hnaw[0] = Hnawi
Hnb[0] = Hnbi
ILaw[0] = ILawi
ILb[0] = ILbi
ILbm[0] = ILbmi
LTA4aw[0] = LTA4awi
LTA4b[0] = LTA4bi
LTC4aw[0] = LTC4awi
LTC4awout[0] = LTC4awouti
LTC4b[0] = LTC4bi
LTC4bout[0] = LTC4bouti
LTD4aw[0] = LTD4awi
LTD4b[0] = LTD4bi
LTE4aw[0] = LTE4awi
LTE4b[0] = LTE4bi
MLblood[0] = MLbloodi
MLintes[0] = MLintesi
ZFairways[0] = ZFairwaysi
ZFblood[0] = ZFbloodi
ZFintes[0] = ZFintesi
end initial conditions

begin initial values
AAawi = 0.1304887
AAbi = 0.2890944
EOaawi = 0.001358713
EOabi = 0.0005383185
EOawi = 1.176558e-06
EObi = 4.68114e-07
EObmi = 1.637391e-06
EOiawi = 1.358698e-07
EOibi = 5.479176e-08
HETEawi = 1.960926
HETEbi = 5.037703
HPETEawi = 1.780604
HPETEbi = 0.2965837
Hnawi = 14162.15
Hnbi = 14995.88
ILawi = 0.6627439
ILbi = 0.5994857
ILbmi = 0.4023394
LTA4awi = 41.69257
LTA4bi = 1.081167
LTC4awi = 6.806687
LTC4awouti = 359.9744
LTC4bi = 0.8869873
LTC4bouti = 2168.571
LTD4awi = 205.7602
LTD4bi = 1308.488
LTE4awi = 322.6366
LTE4bi = 2053.137
MLbloodi = 0.0
MLintesi = 0.0
ZFairwaysi = 0.0
ZFbloodi = 0.0
ZFintesi = 0.0
end initial values

begin assignment rules
ftml := Ceiling[t/T]
ZFbloodconc := fupZF*ZFblood[t]
ZFairwaysconc := ZFairways[t]
MLbloodconc := fupML*MLblood[t]
MLairwaysconc := fupML*MLblood[t]
wEO := (diamEO^3*npi)/6
VaCB := NApmole*VB*wEO*EOab[t]
VaCAW := NApmole*VAW*wEO*EOaaw[t]
GSSGb := (-GSHb + GSpoolb)/2
PLA2D := 1 + Ca/KPLA2Ca
KiAAAA := KiAA
KiHPETEAA := KiAA
OOHb := LOOHb + HPETEb[t]
OHb := LOHb + HETEb[t]
r1 := kox + (Ca*kox2)/KCa2
r2 := kred + (Ca*kred2)/KCa2
RELFLO5b := ((klo*AAb[t])/KAA + k3*(1 + AAb[t]/KiAA)*HPETEb[t])/(kltasyn + k3*(1 + AAb[t]/KiHPETEAA))
FLO3b := FLOa/DFLOab
FLO2b := FLO3b*REDOXb
FLO3tb := FLO3b*(1 + Ca/KCa3)
FLO2tb := FLO2b*(1 + Ca/KCa2)
FLO5HPb := FLO3tb*RELFLO5b
Cb := (GSHb^2*HPETEb[t])/(KmgpxGSH^2*KmgpxHPETE5)
deltaLTCsb := 1 + GSHb/KGSH + LTA4b[t]/KLTA + (GSHb*LTA4b[t])/(KLTA*KLTAGSH) + LTC4b[t]/KLTC
Chedhb := ((khedh2*NADPb)/KdhedhNADP + (khedh3*oxoETEb)/KdhedhoxoETE5)*(1 + NADPHb/KdhedhNADPH + HETEb[t]/KdhedhHETE5)
GSSGaw := (-GSHaw + GSpoolaw)/2
NADPHaw := -NADPaw + NPpoolaw
OOHaw := LOOHaw + HPETEaw[t]
OHaw := LOHaw + HETEaw[t]
RELFLO5aw := ((klo*AAaw[t])/KAA + k3*(1 + AAaw[t]/KiAA)*HPETEaw[t])/(kltasyn + k3*(1 + AAaw[t]/KiHPETEAA))
FLO3aw := FLOa/DFLOaaw
FLO2aw := FLO3aw*REDOXaw
FLO3taw := FLO3aw*(1 + Ca/KCa3)
FLO2taw := FLO2aw*(1 + Ca/KCa2)
FLO5HPaw := FLO3taw*RELFLO5aw
Caw := (GSHaw^2*HPETEaw[t])/(KmgpxGSH^2*KmgpxHPETE5)
Baw := (GSHaw^2*kgpxcat*HPETEaw[t])/(KmgpxGSH^2*KmgpxHPETE5)
denLTCsaw := deltaLTCsaw*KGSH*KLTAGSH*KLTC
Chedhaw := ((khedh2*NADPaw)/KdhedhNADP + (khedh3*oxoETEaw)/KdhedhoxoETE5)*(1 + NADPHaw/KdhedhNADPH + HETEaw[t]/KdhedhHETE5)
Bhedhaw := -((khedh3*khedh4*NADPHaw*oxoETEaw)/(KdhedhNADPH*KdhedhoxoETE5)) + (khedh1*khedh2*NADPaw*HETEaw[t])/(KdhedhHETE5*KdhedhNADP)
Recoccupmigr := ((1 - portionmigr)*ILb[t])/(KdILmigr + ILb[t]) + (fupLT*portionmigr*LTE4b[t])/(KdLTEmigr + fupLT*LTE4b[t])
rinFEV := -((RFEV - RinrelaxFEV)*Sqrt[RFEV/routFEV]) + routFEV
FEV1 := (MAXFEV*rinFEV^2)/RinrelaxFEV^2
timehour := t/60
timeday := 0.0006944444444444445*t
NEObm := 6.02*^8*EObm[t]
NEOb := 6.02*^8*EOb[t]
NEOib := 6.02*^8*EOib[t]
EObtot := EOab[t] + EOb[t] + EOib[t]
NEObtot := 6.02*^8*(EOab[t] + EOb[t] + EOib[t])
NEOaw := 6.02*^8*EOaw[t]
NEOiaw := 6.02*^8*EOiaw[t]
NEOaaw := 6.02*^8*EOaaw[t]
EOawtot := EOaw[t] + EOaaw[t] + EOiaw[t]
NEOawtot := 6.02*^8*(EOaw[t] + EOaaw[t] + EOiaw[t])
EObtotperss := 187617.26078799248*(-0.000533 + EObtot)
NEOawtotperc := 0.00012213964225298783*NEOawtot
NEObtotperc := 0.000023441047289734392*NEObtot
LTC4bpM := 1000000*LTC4b[t]
NEOa := (NEOaaw*VAW + NEOab*VB)/(VAW + VB)
NEO := (NEOawtot*VAW + NEObtot*VB)/(VAW + VB)
NEOawperc := 0.00012213964225298783*NEOawtot
Hnawperc := 0.00706114955514758*Hnaw[t]
NEOperc := 0.0002787394288071625*NEO
OLASM := ((1 - R1portionFEV - R2portionFEV)*Hnaw[t])/(KdHnFEV*(1 + Hnaw[t]/KdHnFEV)) + (R1portionFEV*LTD4aw[t])/(KdLTR1FEV*(1 + MLairwaysconc/EC50MLFEV + LTD4aw[t]/KdLTR1FEV)) + (R2portionFEV*(LTC4awout[t] + LTD4aw[t]))/(KdLTR2FEV*(1 + (LTC4awout[t] + LTD4aw[t])/KdLTR2FEV))
CaFEVLTR2 := (kfeveff*(((1 - R1portionFEV - R2portionFEV)*Hnaw[t])/(KdHnFEV*(1 + Hnaw[t]/KdHnFEV)) + (R2portionFEV*(LTC4awout[t] + LTD4aw[t]))/(KdLTR2FEV*(1 + (LTC4awout[t] + LTD4aw[t])/KdLTR2FEV))))/(-CaoFEV + CaFEV)
CaFEVLTR1 := (kfeveff*R1portionFEV*LTD4aw[t])/((-CaoFEV + CaFEV)*KdLTR1FEV*(1 + MLairwaysconc/EC50MLFEV + LTD4aw[t]/KdLTR1FEV))
LTD4awpers := 0.4860031104199067*LTD4aw[t]
MLuM := MLblood[t]/1000000
LTD4bfree := fupLT*LTD4b[t]
LTsawpg := ((MLTC*LTC4awout[t])/fupLT + (MLTD*LTD4aw[t])/fupLT + (MLTE*LTE4aw[t])/fupLT)/1000
NEOab := 6.02*^8*EOab[t]
FEV1percent := 28.34467120181406*(-3.528 + FEV1)
TSNHn := (CaFEV^nFEV*(RFEV - (RFEV*(RFEV - RinrelaxFEV)^2)^(1/3))*TSN0)/(CaFEV^nFEV + KCaFEV^nFEV)
TSN := (CaFEVex^nFEV*(RFEV - (RFEV*(RFEV - RinrelaxFEV)^2)^(1/3))*TSN0)/(CaFEVex^nFEV + KCaFEV^nFEV)
routFEV := RFEV - (CaFEV^nFEV*(RFEV - (RFEV*(RFEV - RinrelaxFEV)^2)^(1/3)))/(CaFEV^nFEV + KCaFEV^nFEV)
CaFEV := CaoFEV + kfeveff*(((1 - R1portionFEV - R2portionFEV)*Hnaw[t])/(KdHnFEV*(1 + Hnaw[t]/KdHnFEV)) + (R1portionFEV*LTD4aw[t])/(KdLTR1FEV*(1 + MLairwaysconc/EC50MLFEV + LTD4aw[t]/KdLTR1FEV)) + (R2portionFEV*(LTC4awout[t] + LTD4aw[t]))/(KdLTR2FEV*(1 + (LTC4awout[t] + LTD4aw[t])/KdLTR2FEV)))
OLaw := (R1portionEOa*LTD4aw[t])/(KdLT*(1 + MLairwaysconc/KiMLEOa + LTD4aw[t]/KdLT)) + ((1 - R1portionEOa)*(LTC4awout[t] + LTD4aw[t]))/(KdLT2*(1 + (LTC4awout[t] + LTD4aw[t])/KdLT2))
OLb := (fupLT*R1portionEOa*LTD4b[t])/(KdLT*(1 + MLbloodconc/KiMLEOa + (fupLT*LTD4b[t])/KdLT)) + (fupLT*(1 - R1portionEOa)*(LTC4bout[t] + LTD4b[t]))/(KdLT2*(1 + (fupLT*(LTC4bout[t] + LTD4b[t]))/KdLT2))
Ahedhaw := Chedhaw + (1 + NADPaw/KdhedhNADP + oxoETEaw/KdhedhoxoETE5)*((khedh4*NADPHaw)/KdhedhNADPH + (khedh1*HETEaw[t])/KdhedhHETE5)
nomLTCsaw := EtLTCs*(GSHaw*KLTC*kLTCsfow*LTA4aw[t] - KGSH*KLTAGSH*kLTCsback*LTC4aw[t])
deltaLTCsaw := 1 + GSHaw/KGSH + LTA4aw[t]/KLTA + (GSHaw*LTA4aw[t])/(KLTA*KLTAGSH) + LTC4aw[t]/KLTC
Aaw := Caw + (GSHaw^2*(1 + GSSGaw/KdgpxGSSG))/KmgpxGSH^2 + ((1 + HETEaw[t]/KdgpxHETE5)*HPETEaw[t])/KmgpxHPETE5
DFLOaaw := (1 + Ca/KCa3)*(1 + ((1 + Ca/KCa2)*REDOXaw)/(1 + Ca/KCa3) + (al*ZFairwaysconc)/KdZ + (AAaw[t]*(1 + AAaw[t]/KiAAAA))/KAA + RELFLO5aw*(1 + AAaw[t]/KiHPETEAA) + HETEaw[t]/KiHETE + (AAaw[t]*(1 + HETEaw[t]/KAAHETE))/KiAA)
REDOXaw := ((OHaw*r1)/Keox + r2/Kered + (al*ki*(1 + Ca/KCa3)*ZFairwaysconc)/KdZ)/(OOHaw*r1 + r2)
Bhedhb := -((khedh3*khedh4*NADPHb*oxoETEb)/(KdhedhNADPH*KdhedhoxoETE5)) + (khedh1*khedh2*NADPb*HETEb[t])/(KdhedhHETE5*KdhedhNADP)
Ahedhb := Chedhb + (1 + NADPb/KdhedhNADP + oxoETEb/KdhedhoxoETE5)*((khedh4*NADPHb)/KdhedhNADPH + (khedh1*HETEb[t])/KdhedhHETE5)
denLTCsb := deltaLTCsb*KGSH*KLTAGSH*KLTC
nomLTCsb := EtLTCs*(GSHb*KLTC*kLTCsfow*LTA4b[t] - KGSH*KLTAGSH*kLTCsback*LTC4b[t])
Bb := (GSHb^2*kgpxcat*HPETEb[t])/(KmgpxGSH^2*KmgpxHPETE5)
Ab := Cb + (GSHb^2*(1 + GSSGb/KdgpxGSSG))/KmgpxGSH^2 + ((1 + HETEb[t]/KdgpxHETE5)*HPETEb[t])/KmgpxHPETE5
DFLOab := (1 + Ca/KCa3)*(1 + ((1 + Ca/KCa2)*REDOXb)/(1 + Ca/KCa3) + (al*ZFbloodconc)/KdZ + (AAb[t]*(1 + AAb[t]/KiAAAA))/KAA + RELFLO5b*(1 + AAb[t]/KiHPETEAA) + HETEb[t]/KiHETE + (AAb[t]*(1 + HETEb[t]/KAAHETE))/KiAA)
REDOXb := ((OHb*r1)/Keox + r2/Kered + (al*ki*(1 + Ca/KCa3)*ZFbloodconc)/KdZ)/(OOHb*r1 + r2)
KAAHETE := KiHETE
PLA2Ca := Ca/(KPLA2Ca*PLA2D)
NADPHb := -NADPb + NPpoolb
VLTCCAW := NApmole*VAW*wEO*(EOaaw[t] + naEOLTCsyn*(EOaw[t] + EOiaw[t]))
VCAW := NApmole*VAW*wEO*(EOaw[t] + EOaaw[t] + EOiaw[t])
VLTCCB := NApmole*VB*wEO*(EOab[t] + naEOLTCsyn*(EOb[t] + EOib[t]))
VCB := NApmole*VB*wEO*(EOab[t] + EOb[t] + EOib[t])
ftzf := Ceiling[(-960. + t)/T] + Ceiling[(-660. + t)/T] + Ceiling[(-360. + t)/T] + Ceiling[t/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