begin name
mcauley1
end name

begin reactions
reaction1 {1.0}$species1 = {1.0}species2
reaction10 {1.0}species7 = {1.0}species2
reaction11 {1.0}$species12 = {1.0}species7
reaction12 {1.0}species7 = {1.0}species13
reaction13 {1.0}species13 = {1.0}species7
reaction14 {1.0}$species16 = {1.0}species10
reaction15 {1.0}species7 = {1.0}species17
reaction16 {1.0}$species19 = {1.0}species18
reaction17 {1.0}species18 = {1.0}species20
reaction18 {1.0}species17 = {1.0}species7
reaction19 {1.0}species17 = {1.0}species21
reaction2 {1.0}$species3 = {1.0}species2
reaction20 {1.0}species21 = {1.0}species7
reaction21 {1.0}species21 = {1.0}species23
reaction22 {1.0}species23 = {1.0}species7
reaction23 {1.0}species23 = {1.0}species7
reaction24 {1.0}species23 = {1.0}species11
reaction25 {1.0}species23 = {1.0}species11
reaction26 {1.0}$species26 = {1.0}species25
reaction27 {1.0}species25 = {1.0}species27
reaction28 {1.0}species11 = {1.0}species28
reaction29 {1.0}species28 = {1.0}species11
reaction3 {1.0}species4 = {1.0}species5
reaction30 {1.0}species11 = {1.0}species29
reaction31 {1.0}species11 + {1.0}species10 = {1.0}species30
reaction32 {1.0}$species32 = {1.0}species11
reaction33 {1.0}species30 = {1.0}species17
reaction34 {1.0}species30 = {1.0}species23
reaction35 {1.0}species30 = {1.0}species7
reaction4 {1.0}species5 = {1.0}species4
reaction5 {1.0}species5 = {1.0}species6
reaction6 {1.0}species7 = {1.0}species4
reaction7 {1.0}species2 = {1.0}species7
reaction8 {1.0}species2 = {1.0}species8
reaction9 {1.0}$species9 = {1.0}species10
end reactions

begin rate equations
reaction1 = multiplier*reaction1k1*species1
reaction10 = function6[reaction10BCRmax, reaction10BCRt, species7[t], reaction10BS]
reaction11 = function7[reaction11HCSmax, species7[t], reaction11HCSt, reaction11HS]
reaction12 = function8[reaction12k9, species14[t], species7[t]]
reaction13 = function9[reaction13k10, species15[t], species13[t]]
reaction14 = function10[reaction14k11, species11[t]]
reaction15 = reaction15k1*species7[t]
reaction16 = function11[reaction16khrs, species19, species7[t]]
reaction17 = reaction17k1*species18[t]
reaction18 = reaction18k1*species17[t]
reaction19 = function12[reaction19k15, species17[t], species22[t]]
reaction2 = function1[reaction2ICSmax, species2[t], reaction2ICt, reaction2IS]
reaction20 = reaction20k1*species21[t]
reaction21 = function13[reaction21k17, species21[t], species24[t]]
reaction22 = function14[reaction22k18, species23[t], species18[t]]
reaction23 = reaction23k1*species23[t]
reaction24 = function15[reaction24k20, species25[t], species23[t]]
reaction25 = reaction25k1*species23[t]
reaction26 = function16[reaction26kprs, species26, species11[t]]
reaction27 = reaction27k1*species25[t]
reaction28 = function17[reaction28k23, species14[t], species11[t]]
reaction29 = function18[reaction29k24, species15[t], species28[t]]
reaction3 = reaction3k1*species4[t]
reaction30 = reaction30k1*species11[t]
reaction31 = function19[reaction31k26, species11[t], species10[t], species31[t]]
reaction32 = function20[reaction32PCSmax, species11[t], reaction32PPCt, reaction32PCSS]
reaction33 = function21[reaction33k27, species30[t], species33[t]]
reaction34 = function22[reaction34k28, species30[t], species33[t]]
reaction35 = function23[reaction35k29, species30[t], species34[t]]
reaction4 = reaction4k1*species5[t]
reaction5 = reaction5k1*species5[t]
reaction6 = function2[reaction6k5, species7[t], species4[t]]
reaction7 = function3[reaction7k6, species2[t], species5[t]]
reaction8 = function4[reaction8k7, species2[t], species5[t]]
reaction9 = function5[reaction9k8, species11[t]]
end rate equations

begin parameters
multiplier = 1.0
species1 = 304.0
species12 = 0.0
species16 = 0.0
species19 = 600.0
species26 = 575.16
species3 = 0.0
species32 = 0.0
species9 = 0.0
reaction10BCRmax = 2000.0
reaction1k1 = 1.0
reaction2ICSmax = 100.0
reaction2ICt = 3120.0
reaction2IS = 5.0
reaction3k1 = 6.0
reaction4k1 = 4.29
reaction5k1 = 0.856
reaction6k5 = 2.66
reaction7k6 = 0.0005286
reaction8k7 = 0.0005286
reaction9k8 = 0.0005
reaction10BCRt = 55326.0
reaction10BS = 5.0
reaction11HCSmax = 500.0
reaction11HCSt = 93925.0
reaction11HS = 5.0
reaction12k9 = 1.0
reaction13k10 = 5.998
reaction14k11 = 0.005
reaction15k1 = 0.016
reaction16khrs = 100.0
reaction17k1 = 0.01
reaction18k1 = 0.0496
reaction19k15 = 0.43
reaction20k1 = 0.054
reaction21k17 = 0.38
reaction22k18 = 0.068
reaction23k1 = 0.005
reaction24k20 = 0.00675
reaction25k1 = "5e-06"
reaction26kprs = 100.0
reaction27k1 = 0.01
reaction28k23 = 0.017386
reaction29k24 = 0.1068
reaction30k1 = 0.0005
reaction31k26 = "1.5e-05"
reaction32PCSmax = 500.0
reaction32PPCt = 80342.0
reaction32PCSS = 5.0
reaction33k27 = 0.01
reaction34k28 = 0.001
reaction35k29 = 0.05
Excreted = 1.0
HepaticTissue = 1.0
Intake = 1.0
Intestine = 1.0
PeripheralTissue = 1.0
Plasma = 1.0
end parameters

begin initial conditions
species10[0] = species10i
species11[0] = species11i
species13[0] = species13i
species14[0] = species14i
species15[0] = species15i
species17[0] = species17i
species18[0] = species18i
species2[0] = species2i
species20[0] = species20i
species21[0] = species21i
species22[0] = species22i
species23[0] = species23i
species24[0] = species24i
species25[0] = species25i
species27[0] = species27i
species28[0] = species28i
species29[0] = species29i
species30[0] = species30i
species31[0] = species31i
species33[0] = species33i
species34[0] = species34i
species4[0] = species4i
species5[0] = species5i
species6[0] = species6i
species7[0] = species7i
species8[0] = species8i
end initial conditions

begin initial values
species10i = 100.0
species11i = 57516.0
species13i = 10000.0
species14i = 100.0
species15i = 100.0
species17i = 20.0
species18i = 100.0
species2i = 3150.0
species20i = 0.0
species21i = 20.0
species22i = 100.0
species23i = 100.0
species24i = 100.0
species25i = 100.0
species27i = 0.0
species28i = 9363.0
species29i = 0.0
species30i = 45.0
species31i = 100.0
species33i = 100.0
species34i = 100.0
species4i = 400.0
species5i = 467.0
species6i = 0.0
species7i = 60000.0
species8i = 0.0
end initial values

begin assignment rules
end assignment rules

begin function definitions
function4[k7_,IC_,IBS_] := IBS*IC*k7
function5[k8_,PFC_] := k8*PFC
function6[BCRmax_,BCRt_,HFC_,BS_] := BCRmax/(1 + (BCRt/HFC)^BS)
function10[k11_,PFC_] := k11*PFC
function8[k9_,ACAT_,HFC_] := ACAT*HFC*k9
function16[kprs_,PLDLRsS_,PFC_] := (kprs*PLDLRsS)/PFC
function15[k20_,PLDLRs_,LDLC_] := k20*LDLC*PLDLRs
function14[k18_,LDLC_,HLDLRs_] := HLDLRs*k18*LDLC
function13[k17_,IDLC_,HSL_] := HSL*IDLC*k17
function12[k15_,VLDLC_,LPL_] := k15*LPL*VLDLC
function11[khrs_,HLDLRsS_,HFC_] := (HLDLRsS*khrs)/HFC
function1[ICSmax_,IC_,ICt_,IS_] := ICSmax/(1 + (IC/ICt)^IS)
function2[k5_,HFC_,HBS_] := (HFC*k5)/HBS
function3[k6_,IC_,IBS_] := IBS*IC*k6
function7[HCSmax_,HFC_,HCSt_,HS_] := HCSmax/(1 + (HFC/HCSt)^HS)
function9[k10_,CEH_,HCE_] := CEH*HCE*k10
function23[k29_,HDLC_,SRB1_] := HDLC*k29*SRB1
function22[k28_,HDLC_,CETP_] := CETP*HDLC*k28
function21[k27_,HDLC_,CETP_] := CETP*HDLC*k27
function20[PCSmax_,PFC_,PPCt_,PCSS_] := PCSmax/(1 + (PFC/PPCt)^PCSS)
function19[k26_,PFC_,NHDL_,LCAT_] := k26*LCAT*NHDL*PFC
function18[k24_,CEH_,PCE_] := CEH*k24*PCE
function17[k23_,ACAT_,PFC_] := ACAT*k23*PFC
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units