begin name
benson1
end name

begin reactions
absorp {1.0}PFMgut = {1.0}PFMp
dist {1.0}PFMp = {1.0}PFMr
elim {1.0}PFMp = {1.0}$EXTERNAL
vAUEb {1.0}Ab = {1.0}$EXTERNAL
vAUEr {1.0}Ar = {1.0}$EXTERNAL
vAbm {1.0}Ab = {1.0}Am
vAdegrb {1.0}Ab = {1.0}$EXTERNAL
vAdegrm {1.0}Am = {1.0}$EXTERNAL
vAdegrr {1.0}Ar = {1.0}$EXTERNAL
vAmp {1.0}Am = {1.0}Ap
vArp {1.0}Ar = {1.0}Ap
vAsynb {1.0}NAPEb = {1.0}Ab
vAsynr {1.0}NAPEr = {1.0}Ar
vFAAHdegrb {1.0}FAAHb = {1.0}$EXTERNAL
vFAAHdegrm {1.0}FAAHm = {1.0}$EXTERNAL
vFAAHdegrr {1.0}FAAHr = {1.0}$EXTERNAL
vFAAHinhb {1.0}FAAHb = {1.0}FAAHinhb
vFAAHinhdegrb {1.0}FAAHinhb = {1.0}$EXTERNAL
vFAAHinhdegrm {1.0}FAAHinhm = {1.0}$EXTERNAL
vFAAHinhdegrr {1.0}FAAHinhr = {1.0}$EXTERNAL
vFAAHinhm {1.0}FAAHm = {1.0}FAAHinhm
vFAAHinhr {1.0}FAAHr = {1.0}FAAHinhr
vFAAHsynb {1.0}$EXTERNAL = {1.0}FAAHb
vFAAHsynm {1.0}$EXTERNAL = {1.0}FAAHm
vFAAHsynr {1.0}$EXTERNAL = {1.0}FAAHr
vLUEb {1.0}Lb = {1.0}$EXTERNAL
vLUEr {1.0}Lr = {1.0}$EXTERNAL
vLbm {1.0}Lb = {1.0}Lm
vLdegrb {1.0}Lb = {1.0}$EXTERNAL
vLdegrm {1.0}Lm = {1.0}$EXTERNAL
vLdegrr {1.0}Lr = {1.0}$EXTERNAL
vLmp {1.0}Lm = {1.0}Lp
vLrp {1.0}Lr = {1.0}Lp
vLsynb {1.0}NLPEb = {1.0}Lb
vLsynr {1.0}NLPEr = {1.0}Lr
vNAPEsynb {1.0}$EXTERNAL = {1.0}NAPEb
vNAPEsynr {1.0}$EXTERNAL = {1.0}NAPEr
vNLPEsynb {1.0}$EXTERNAL = {1.0}NLPEb
vNLPEsynr {1.0}$EXTERNAL = {1.0}NLPEr
vNOPEsynb {1.0}$EXTERNAL = {1.0}NOPEb
vNOPEsynr {1.0}$EXTERNAL = {1.0}NOPEr
vNPPEsynb {1.0}$EXTERNAL = {1.0}NPPEb
vNPPEsynr {1.0}$EXTERNAL = {1.0}NPPEr
vNSPEsynb {1.0}$EXTERNAL = {1.0}NSPEb
vNSPEsynr {1.0}$EXTERNAL = {1.0}NSPEr
vOUEb {1.0}Ob = {1.0}$EXTERNAL
vOUEr {1.0}Or = {1.0}$EXTERNAL
vObm {1.0}Ob = {1.0}Om
vOdegrb {1.0}Ob = {1.0}$EXTERNAL
vOdegrm {1.0}Om = {1.0}$EXTERNAL
vOdegrr {1.0}Or = {1.0}$EXTERNAL
vOrp {1.0}Or = {1.0}Op
vOsynb {1.0}NOPEb = {1.0}Ob
vOsynr {1.0}NOPEr = {1.0}Or
vPUEb {1.0}Pb = {1.0}$EXTERNAL
vPUEr {1.0}Pr = {1.0}$EXTERNAL
vPbm {1.0}Pb = {1.0}Pm
vPdegrb {1.0}Pb = {1.0}$EXTERNAL
vPdegrm {1.0}Pm = {1.0}$EXTERNAL
vPdegrr {1.0}Pr = {1.0}$EXTERNAL
vPmp {1.0}Pm = {1.0}Pp
vPrp {1.0}Pr = {1.0}Pp
vPsynb {1.0}NPPEb = {1.0}Pb
vPsynr {1.0}NPPEr = {1.0}Pr
vSUEb {1.0}Sb = {1.0}$EXTERNAL
vSUEr {1.0}Sr = {1.0}$EXTERNAL
vSbm {1.0}Sb = {1.0}Sm
vSdegrb {1.0}Sb = {1.0}$EXTERNAL
vSdegrm {1.0}Sm = {1.0}$EXTERNAL
vSdegrr {1.0}Sr = {1.0}$EXTERNAL
vSmp {1.0}Sm = {1.0}Sp
vSrp {1.0}Sr = {1.0}Sp
vSsynb {1.0}NSPEb = {1.0}Sb
vSsynr {1.0}NSPEr = {1.0}Sr
vomp {1.0}Om = {1.0}Op
end reactions

begin rate equations
absorp = kabsPFM*MD
dist = koutPFM*PFMp[t] - kinPFM*PFMr[t]
elim = klinearPFM*PFMp[t] + (VmPFM*PFMp[t])/(VssPFM*(KmPFM + PFMp[t]/VssPFM))
vAUEb = BRAIN*bFAAHBrain*kclA*Ab[t]
vAUEr = cNAAAROB*kclA*Ar[t]
vAbm = (ktrmpA*MEC*(Ab[t] - Am[t]))/(KmpmA + Ab[t] + Am[t])
vAdegrb = (aFAAHA*BRAIN*kcatFAAH*Ab[t]*FAAHb[t])/(FAAHDb*KmFAAHA)
vAdegrm = (aFAAHA*kcatFAAH*MEC*Am[t]*FAAHm[t])/(FAAHDm*KmFAAHA)
vAdegrr = (aFAAHA*kcatFAAH*ROB*Ar[t]*FAAHr[t])/(FAAHDr*KmFAAHA)
vAmp = (ktrmpA*MEC*(Am[t] - KtrpmA*Ap[t]))/(KmpmA + Am[t] + Ap[t])
vArp = (ktrrp*PLASMA*(-(KtrprA*Ap[t]) + Ar[t]))/(KmpmA + Ap[t] + Ar[t])
vAsynb = (BRAIN*kNAPE*PLDb*NAPEb[t])/(denb*KmNAPE)
vAsynr = (kNAPE*PLDr*ROB*NAPEr[t])/(denr*KmNAPE)
vFAAHdegrb = BRAIN*kdegFAAH*FAAHb[t]
vFAAHdegrm = kdegFAAH*MEC*FAAHm[t]
vFAAHdegrr = kdegFAAH*ROB*FAAHr[t]
vFAAHinhb = BRAIN*kinh*PFb*FAAHb[t]
vFAAHinhdegrb = BRAIN*kdegFAAH*FAAHinhb[t]
vFAAHinhdegrm = kdegFAAH*MEC*FAAHinhm[t]
vFAAHinhdegrr = kdegFAAH*ROB*FAAHinhr[t]
vFAAHinhm = kinh*MEC*PFm*FAAHm[t]
vFAAHinhr = kinh*PFr*ROB*FAAHr[t]
vFAAHsynb = BRAIN*bFAAHBrain*FAAHt*kdegFAAH
vFAAHsynm = bFAAHMEC*FAAHt*kdegFAAH*MEC
vFAAHsynr = cFAAHROB*FAAHt*kdegFAAH
vLUEb = BRAIN*bFAAHBrain*kclL*Lb[t]
vLUEr = cNAAAROB*kclL*Lr[t]
vLbm = ktrmpL*MEC*(Lb[t] - Lm[t])
vLdegrb = (aFAAHL*BRAIN*kcatFAAH*FAAHb[t]*Lb[t])/(FAAHDb*KmFAAHL)
vLdegrm = (aFAAHL*kcatFAAH*MEC*FAAHm[t]*Lm[t])/(FAAHDm*KmFAAHL)
vLdegrr = (aFAAHL*kcatFAAH*ROB*FAAHr[t]*Lr[t])/(FAAHDr*KmFAAHL)
vLmp = ktrmpL*MEC*(Lm[t] - KtrpmL*Lp[t])
vLrp = ktrrp*PLASMA*(-(KtrprL*Lp[t]) + Lr[t])
vLsynb = (BRAIN*kNLPE*PLDb*NLPEb[t])/(denb*KmNLPE)
vLsynr = (kNLPE*PLDr*ROB*NLPEr[t])/(denr*KmNLPE)
vNAPEsynb = aNATA*BRAIN*bNATBrain*pA*VmaxNAT
vNAPEsynr = aNATA*cNATROB*pA*VmaxNAT
vNLPEsynb = aNATL*BRAIN*bNATBrain*pL*VmaxNAT
vNLPEsynr = aNATL*cNATROB*pL*VmaxNAT
vNOPEsynb = aNATO*BRAIN*bNATBrain*pO*VmaxNAT
vNOPEsynr = aNATO*cNATROB*pO*VmaxNAT
vNPPEsynb = aNATP*BRAIN*bNATBrain*pP*VmaxNAT
vNPPEsynr = aNATP*cNATROB*pP*VmaxNAT
vNSPEsynb = aNATS*BRAIN*bNATBrain*pS*VmaxNAT
vNSPEsynr = aNATS*cNATROB*pS*VmaxNAT
vOUEb = BRAIN*bFAAHBrain*kclO*Ob[t]
vOUEr = cNAAAROB*kclO*Or[t]
vObm = ktrmpO*MEC*(Ob[t] - Om[t])
vOdegrb = (aFAAHO*BRAIN*kcatFAAH*FAAHb[t]*Ob[t])/(FAAHDb*KmFAAHO)
vOdegrm = (aFAAHO*kcatFAAH*MEC*FAAHm[t]*Om[t])/(FAAHDm*KmFAAHO)
vOdegrr = (aFAAHO*kcatFAAH*ROB*FAAHr[t]*Or[t])/(FAAHDr*KmFAAHO)
vOrp = ktrrp*PLASMA*(-(KtrprO*Op[t]) + Or[t])
vOsynb = (BRAIN*kNOPE*PLDb*NOPEb[t])/(denb*KmNOPE)
vOsynr = (kNOPE*PLDr*ROB*NOPEr[t])/(denr*KmNOPE)
vPUEb = BRAIN*bFAAHBrain*kclP*Pb[t]
vPUEr = cNAAAROB*kclP*Pr[t]
vPbm = ktrmpP*MEC*(Pb[t] - Pm[t])
vPdegrb = (aFAAHP*BRAIN*kcatFAAH*FAAHb[t]*Pb[t])/(FAAHDb*KmFAAHP)
vPdegrm = (aFAAHP*kcatFAAH*MEC*FAAHm[t]*Pm[t])/(FAAHDm*KmFAAHP)
vPdegrr = (aFAAHP*kcatFAAH*ROB*FAAHr[t]*Pr[t])/(FAAHDr*KmFAAHP)
vPmp = ktrmpP*MEC*(Pm[t] - KtrpmP*Pp[t])
vPrp = ktrrp*PLASMA*(-(KtrprP*Pp[t]) + Pr[t])
vPsynb = (BRAIN*kNPPE*PLDb*NPPEb[t])/(denb*KmNPPE)
vPsynr = (kNPPE*PLDr*ROB*NPPEr[t])/(denr*KmNPPE)
vSUEb = BRAIN*bFAAHBrain*kclS*Sb[t]
vSUEr = cNAAAROB*kclS*Sr[t]
vSbm = ktrmpS*MEC*(Sb[t] - Sm[t])
vSdegrb = (aFAAHS*BRAIN*kcatFAAH*FAAHb[t]*Sb[t])/(FAAHDb*KmFAAHS)
vSdegrm = (aFAAHS*kcatFAAH*MEC*FAAHm[t]*Sm[t])/(FAAHDm*KmFAAHS)
vSdegrr = (aFAAHS*kcatFAAH*ROB*FAAHr[t]*Sr[t])/(FAAHDr*KmFAAHS)
vSmp = ktrmpS*MEC*(Sm[t] - KtrpmS*Sp[t])
vSrp = ktrrp*PLASMA*(-(KtrprS*Sp[t]) + Sr[t])
vSsynb = (BRAIN*kNSPE*PLDb*NSPEb[t])/(denb*KmNSPE)
vSsynr = (kNSPE*PLDr*ROB*NSPEr[t])/(denr*KmNSPE)
vomp = ktrmpO*MEC*(Om[t] - KtrpmO*Op[t])
end rate equations

begin parameters
AG2b = 0.0
Dose = 10.0
ED50 = 0.53
EmaxPFM = 0.773
FAAHt = 78.0
Gut = 1.65
H = 0.48
Heart = 0.31
KdAG2 = 3424.0
KdCB1A = 239.2
KiA = 230.0
KiL = 1000.0
KiO = 240.0
KiP = 6700.0
KiS = 840.0
Kidney = 0.28
KmFAAHA = 8200.0
KmFAAHL = 10800.0
KmFAAHO = 52200.0
KmFAAHP = 543000.0
KmFAAHS = 10000.0
KmNAPE = 2800.0
KmNLPE = 1000.0
KmNOPE = 2900.0
KmNPPE = 3300.0
KmNSPE = 3400.0
KmPFM = 26.1
KmpmA = 1.0
KpbPF = 1.3
KpmPF = 1.3
KprPF = 1.5
KtrpmA = 1.89
KtrpmL = 2.77
KtrpmO = 9.07
KtrpmP = 2.65
KtrpmS = 30.01
KtrprA = 0.62
KtrprL = 0.89
KtrprO = 2.8
KtrprP = 0.85
KtrprS = 9.19
LEU = 0.025
LIVER = 1.69
Leucocytes = 0.025
Lungs = 1.172
MA = 347.5
ML = 323.5
MO = 325.5
MP = 299.5
MPF = 455.4
MS = 321.5
Muscles = 35.0
PLDb = 10000000.0
PLDr = 10000000.0
Pancreas = 0.077
RPF = 0.6
Spleen = 0.192
Testis = 0.036
Thymus = 0.029
VmPFM = 1511.0
VmaxNAT = 300.0
VssPFM = 58.328
aFAAHA = 1.0
aFAAHL = 1.15
aFAAHO = 5.7
aFAAHP = 37.8
aFAAHS = 1.0
aNATA = 1.0
aNATL = 8.6
aNATO = 13.0
aNATP = 0.42
aNATS = 1.0
bFAAHBrain = 0.197
bFAAHGut = 0.034
bFAAHKidney = 0.069
bFAAHLeucocytes = 0.0
bFAAHLiver = 1.0
bFAAHLungs = 0.032
bFAAHMEC = 0.137
bFAAHSpleen = 0.03
bFAAHTestis = 0.126
bNAAABrain = 0.6
bNAAAGut = 0.2
bNAAAHeart = 0.2
bNAAAKidney = 0.6
bNAAALiver = 1.0
bNAAALungs = 14.0
bNAAASpleen = 8.0
bNAAATestis = 0.6
bNAAAThymus = 4.0
bNATBrain = 1.667
bNATHeart = 1.0
bNATKidney = 0.667
bNATLeucocytes = 0.0
bNATLungs = 0.033
bNATMuscles = 0.333
bNATPancreas = 0.333
bNATTestis = 0.667
e1 = 0.0
kNAPE = 202.0
kNLPE = 100.0
kNOPE = 230.0
kNPPE = 270.0
kNSPE = 280.0
kdegFAAH = 0.0051
kinh = 1.1
kabsPFM = 2.2
kcatFAAH = 18000.0
kclA = 1.74
kclL = 1.25
kclO = 2.5
kclP = 2.61
kclS = 1.2
kinPFM = 0.117
klinearPFM = 0.0803
koutPFM = 0.18
ktrmpA = 150.0
ktrmpL = 0.0
ktrmpO = 10.0
ktrmpP = 10.0
ktrmpS = 10.0
ktrrp = 100.0
pA = 0.051
pL = 0.016
pO = 0.098
pP = 0.615
pS = 0.191
t = 0.0
BRAIN = 1.45000004768372
Default = 1.0
MEC = "1.49999996210681e-05"
PLASMA = 2.6489999294281
ROB = 65.3000030517578
end parameters

begin initial conditions
Ab[0] = Abi
Am[0] = Ami
Ap[0] = Api
Ar[0] = Ari
FAAHb[0] = FAAHbi
FAAHm[0] = FAAHmi
FAAHr[0] = FAAHri
FAAHinhb[0] = FAAHinhbi
FAAHinhm[0] = FAAHinhmi
FAAHinhr[0] = FAAHinhri
Lb[0] = Lbi
Lm[0] = Lmi
Lp[0] = Lpi
Lr[0] = Lri
NAPEb[0] = NAPEbi
NAPEr[0] = NAPEri
NLPEb[0] = NLPEbi
NLPEr[0] = NLPEri
NOPEb[0] = NOPEbi
NOPEr[0] = NOPEri
NPPEb[0] = NPPEbi
NPPEr[0] = NPPEri
NSPEb[0] = NSPEbi
NSPEr[0] = NSPEri
Ob[0] = Obi
Om[0] = Omi
Op[0] = Opi
Or[0] = Ori
PFMgut[0] = PFMguti
PFMp[0] = PFMpi
PFMr[0] = PFMri
Pb[0] = Pbi
Pm[0] = Pmi
Pp[0] = Ppi
Pr[0] = Pri
Sb[0] = Sbi
Sm[0] = Smi
Sp[0] = Spi
Sr[0] = Sri
end initial conditions

begin initial values
Abi = 0.7493309
Ami = 0.97761
Api = 0.8740574
Ari = 0.5419204
FAAHbi = 15.366
FAAHmi = 10.686
FAAHri = 2.165868
FAAHinhbi = 0.0
FAAHinhmi = 0.0
FAAHinhri = 0.0
Lbi = 2.319571
Lmi = 0.0
Lpi = 1.916254
Lri = 1.705466
NAPEbi = 3.879041e-05
NAPEri = 4.241633e-06
NLPEbi = 7.550331e-05
NLPEri = 8.256095e-06
NOPEbi = 0.0008814287
NOPEri = 9.638198e-05
NPPEbi = 0.0001732296
NPPEri = 1.894222e-05
NSPEbi = 0.0001272629
NSPEri = 1.391587e-05
Obi = 20.77858
Omi = 16.3219
Opi = 5.085073
Ori = 14.23822
PFMguti = 0.0
PFMpi = 0.0
PFMri = 0.0
Pbi = 6.541209
Pmi = 5.809415
Ppi = 4.849307
Pri = 4.121915
Sbi = 3.427807
Smi = 2.968774
Spi = 0.273772
Sri = 2.515968
end initial values

begin assignment rules
tid := 0.041666666666666664*t
CB1b := (AG2b/KdAG2 + Ab[t]/KdCB1A)/(1 + AG2b/KdAG2 + Ab[t]/KdCB1A)
SGp := (MS*Sp[t])/1000
LGp := (ML*Lp[t])/1000
PGp := (MP*Pp[t])/1000
PFGp := PFMp[t]/(1000*VssPFM)
Fr := (1500.*FAAHr[t])/(FAAHinhr[t] + FAAHr[t])
FAAHDb := 1 + Ab[t]/KmFAAHA + Lb[t]/KmFAAHL + Ob[t]/KmFAAHO + Pb[t]/KmFAAHP + Sb[t]/KmFAAHS
PFm := KpmPF*PFp
PFr := KprPF*PFp
denb := 1 + slag1b + slag2b
MD := 1000000*Dose*FPFM + PFMgut[t]
PFb := KpbPF*PFp
FAAHDr := 1 + Ar[t]/KmFAAHA + Lr[t]/KmFAAHL + Or[t]/KmFAAHO + Pr[t]/KmFAAHP + Sr[t]/KmFAAHS
denr := 1 + slag1r + slag2r
cFAAHROB := bFAAHGut*Gut + bFAAHKidney*Kidney + bFAAHLeucocytes*Leucocytes + bFAAHLiver*LIVER + bFAAHLungs*Lungs + bFAAHSpleen*Spleen + bFAAHTestis*Testis
AGp := (MA*Ap[t])/1000
OGp := (MO*Op[t])/1000
FAAHDm := 1 + Am[t]/KmFAAHA + Lm[t]/KmFAAHL + Om[t]/KmFAAHO + Pm[t]/KmFAAHP + Sm[t]/KmFAAHS
cNAAAROB := bNAAAGut*Gut + bNAAAHeart*Heart + bNAAAKidney*Kidney + bNAAALiver*LIVER + bNAAALungs*Lungs + bNAAASpleen*Spleen + bNAAAThymus*Testis^2 + bNAAAThymus*Thymus
slag2r := Ar[t]/KiA + Lr[t]/KiL + Or[t]/KiO + Pr[t]/KiP + Sr[t]/KiS
slag1r := NAPEr[t]/KmNAPE + NLPEr[t]/KmNLPE + NOPEr[t]/KmNOPE + NPPEr[t]/KmNPPE + NSPEr[t]/KmNSPE
cNATROB := bNATHeart*Heart + bNATKidney*Kidney + bNATLeucocytes*Leucocytes + bNATLungs*Lungs + bNATMuscles*Muscles + bNATPancreas*Pancreas + bNATTestis*Testis
slag2b := Ab[t]/KiA + Lb[t]/KiL + Ob[t]/KiO + Pb[t]/KiP + Sb[t]/KiS
slag1b := NAPEb[t]/KmNAPE + NLPEb[t]/KmNLPE + NOPEb[t]/KmNOPE + NPPEb[t]/KmNPPE + NSPEb[t]/KmNSPE
PFp := PFMp[t]/(MPF*VssPFM)
FPFM := (Dose*EmaxPFM)/(Dose + ED50)
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