(* Generated by JWS Online *) (* This is an experimental feature of JWS Online. Please report any mistakes.*) (* Note that the following notable SBML entities or features are not supported in notebook outputyet: *) (* Events *) (* Constraints *) (* Units and UnitDefinitions *) (* AlgebraicRules *) (* conversionFactors *) variables = { APC\[LetterSpace]PS\[LetterSpace]VIIIa\[LetterSpace]l[t], APC\[LetterSpace]PS\[LetterSpace]Va\[LetterSpace]l[t], APC\[LetterSpace]PS\[LetterSpace]l[t], APC\[LetterSpace]f[t], APC\[LetterSpace]l[t], AT\[LetterSpace]f[t], II\[LetterSpace]f[t], II\[LetterSpace]l[t], IIa\[LetterSpace]AT\[LetterSpace]f[t], IIa\[LetterSpace]TM\[LetterSpace]PC\[LetterSpace]l[t], IIa\[LetterSpace]TM\[LetterSpace]l[t], IIa\[LetterSpace]f[t], IX\[LetterSpace]f[t], IX\[LetterSpace]l[t], IXa\[LetterSpace]AT\[LetterSpace]f[t], IXa\[LetterSpace]VIIIa\[LetterSpace]X\[LetterSpace]l[t], IXa\[LetterSpace]VIIIa\[LetterSpace]l[t], IXa\[LetterSpace]f[t], IXa\[LetterSpace]l[t], LIPID[t], PC\[LetterSpace]f[t], PC\[LetterSpace]l[t], PS\[LetterSpace]f[t], PS\[LetterSpace]l[t], TFPI\[LetterSpace]Xa\[LetterSpace]TF\[LetterSpace]VIIa\[LetterSpace]l[t], TFPI\[LetterSpace]Xa\[LetterSpace]l[t], TFPI\[LetterSpace]f[t], TF\[LetterSpace]VII\[LetterSpace]Xa\[LetterSpace]l[t], TF\[LetterSpace]VII\[LetterSpace]l[t], TF\[LetterSpace]VIIa\[LetterSpace]IX\[LetterSpace]l[t], TF\[LetterSpace]VIIa\[LetterSpace]IXa\[LetterSpace]l[t], TF\[LetterSpace]VIIa\[LetterSpace]X\[LetterSpace]l[t], TF\[LetterSpace]VIIa\[LetterSpace]Xa\[LetterSpace]l[t], TF\[LetterSpace]VIIa\[LetterSpace]l[t], TF\[LetterSpace]l[t], TM\[LetterSpace]l[t], VIII\[LetterSpace]IIa\[LetterSpace]l[t], VIII\[LetterSpace]Xa\[LetterSpace]l[t], VIII\[LetterSpace]f[t], VIII\[LetterSpace]l[t], VIII\[LetterSpace]mIIa\[LetterSpace]l[t], VIIIa\[LetterSpace]f[t], VIIIa\[LetterSpace]l[t], VIIIai\[LetterSpace]f[t], VIIIai\[LetterSpace]l[t], VII\[LetterSpace]Xa\[LetterSpace]l[t], VII\[LetterSpace]f[t], VII\[LetterSpace]l[t], VIIa\[LetterSpace]f[t], VIIa\[LetterSpace]l[t], V\[LetterSpace]IIa\[LetterSpace]l[t], V\[LetterSpace]Xa\[LetterSpace]l[t], V\[LetterSpace]f[t], V\[LetterSpace]l[t], V\[LetterSpace]mIIa\[LetterSpace]l[t], Va\[LetterSpace]f[t], Va\[LetterSpace]l[t], Vai\[LetterSpace]f[t], Vai\[LetterSpace]l[t], XI\[LetterSpace]IIa\[LetterSpace]l[t], XI\[LetterSpace]f[t], XIa\[LetterSpace]IX\[LetterSpace]l[t], XIa\[LetterSpace]l[t], X\[LetterSpace]f[t], X\[LetterSpace]l[t], Xa\[LetterSpace]AT\[LetterSpace]f[t], Xa\[LetterSpace]Va\[LetterSpace]II\[LetterSpace]l[t], Xa\[LetterSpace]Va\[LetterSpace]l[t], Xa\[LetterSpace]Va\[LetterSpace]mIIa\[LetterSpace]l[t], Xa\[LetterSpace]f[t], Xa\[LetterSpace]l[t], mIIa\[LetterSpace]AT\[LetterSpace]l[t], mIIa\[LetterSpace]f[t], mIIa\[LetterSpace]l[t] }; initialValues = { APC\[LetterSpace]PS\[LetterSpace]VIIIa\[LetterSpace]l[0] == 0.0, APC\[LetterSpace]PS\[LetterSpace]Va\[LetterSpace]l[0] == 0.0, APC\[LetterSpace]PS\[LetterSpace]l[0] == 0.0, APC\[LetterSpace]f[0] == 0.0, APC\[LetterSpace]l[0] == 0.0, AT\[LetterSpace]f[0] == 3400.0, II\[LetterSpace]f[0] == 1400.0, II\[LetterSpace]l[0] == 0.0, IIa\[LetterSpace]AT\[LetterSpace]f[0] == 0.0, IIa\[LetterSpace]TM\[LetterSpace]PC\[LetterSpace]l[0] == 0.0, IIa\[LetterSpace]TM\[LetterSpace]l[0] == 0.0, IIa\[LetterSpace]f[0] == 0.0, IX\[LetterSpace]f[0] == 90.0, IX\[LetterSpace]l[0] == 0.0, IXa\[LetterSpace]AT\[LetterSpace]f[0] == 0.0, IXa\[LetterSpace]VIIIa\[LetterSpace]X\[LetterSpace]l[0] == 0.0, IXa\[LetterSpace]VIIIa\[LetterSpace]l[0] == 0.0, IXa\[LetterSpace]f[0] == 0.0, IXa\[LetterSpace]l[0] == 0.0, LIPID[0] == 849079.0, PC\[LetterSpace]f[0] == 60.0, PC\[LetterSpace]l[0] == 0.0, PS\[LetterSpace]f[0] == 300.0, PS\[LetterSpace]l[0] == 0.0, TFPI\[LetterSpace]Xa\[LetterSpace]TF\[LetterSpace]VIIa\[LetterSpace]l[0] == 0.0, TFPI\[LetterSpace]Xa\[LetterSpace]l[0] == 0.0, TFPI\[LetterSpace]f[0] == 2.5, TF\[LetterSpace]VII\[LetterSpace]Xa\[LetterSpace]l[0] == 0.0, TF\[LetterSpace]VII\[LetterSpace]l[0] == 0.0, TF\[LetterSpace]VIIa\[LetterSpace]IX\[LetterSpace]l[0] == 0.0, TF\[LetterSpace]VIIa\[LetterSpace]IXa\[LetterSpace]l[0] == 0.0, TF\[LetterSpace]VIIa\[LetterSpace]X\[LetterSpace]l[0] == 0.0, TF\[LetterSpace]VIIa\[LetterSpace]Xa\[LetterSpace]l[0] == 0.0, TF\[LetterSpace]VIIa\[LetterSpace]l[0] == 0.0, TF\[LetterSpace]l[0] == 0.005, TM\[LetterSpace]l[0] == 1.0, VIII\[LetterSpace]IIa\[LetterSpace]l[0] == 0.0, VIII\[LetterSpace]Xa\[LetterSpace]l[0] == 0.0, VIII\[LetterSpace]f[0] == 0.7, VIII\[LetterSpace]l[0] == 0.0, VIII\[LetterSpace]mIIa\[LetterSpace]l[0] == 0.0, VIIIa\[LetterSpace]f[0] == 0.0, VIIIa\[LetterSpace]l[0] == 0.0, VIIIai\[LetterSpace]f[0] == 0.0, VIIIai\[LetterSpace]l[0] == 0.0, VII\[LetterSpace]Xa\[LetterSpace]l[0] == 0.0, VII\[LetterSpace]f[0] == 10.0, VII\[LetterSpace]l[0] == 0.0, VIIa\[LetterSpace]f[0] == 0.1, VIIa\[LetterSpace]l[0] == 0.0, V\[LetterSpace]IIa\[LetterSpace]l[0] == 0.0, V\[LetterSpace]Xa\[LetterSpace]l[0] == 0.0, V\[LetterSpace]f[0] == 20.0, V\[LetterSpace]l[0] == 0.0, V\[LetterSpace]mIIa\[LetterSpace]l[0] == 0.0, Va\[LetterSpace]f[0] == 0.0, Va\[LetterSpace]l[0] == 0.0, Vai\[LetterSpace]f[0] == 0.0, Vai\[LetterSpace]l[0] == 0.0, XI\[LetterSpace]IIa\[LetterSpace]l[0] == 0.0, XI\[LetterSpace]f[0] == 30.0, XIa\[LetterSpace]IX\[LetterSpace]l[0] == 0.0, XIa\[LetterSpace]l[0] == 0.0, X\[LetterSpace]f[0] == 170.0, X\[LetterSpace]l[0] == 0.0, Xa\[LetterSpace]AT\[LetterSpace]f[0] == 0.0, Xa\[LetterSpace]Va\[LetterSpace]II\[LetterSpace]l[0] == 0.0, Xa\[LetterSpace]Va\[LetterSpace]l[0] == 0.0, Xa\[LetterSpace]Va\[LetterSpace]mIIa\[LetterSpace]l[0] == 0.0, Xa\[LetterSpace]f[0] == 0.0, Xa\[LetterSpace]l[0] == 0.0, mIIa\[LetterSpace]AT\[LetterSpace]l[0] == 0.0, mIIa\[LetterSpace]f[0] == 0.0, mIIa\[LetterSpace]l[0] == 0.0 }; rates = { LB1, LB10, LB11, LB12, LB13, LB14, LB15, LB16, LB17, LB2, LB3, LB4, LB5, LB6, LB7, LB8, LB9, R1, R10, R10b, R11, R11b, R12, R12b, R13, R14, R15, R15b, R16, R16b, R17, R17b, R18, R18b, R19, R19b, R2, R20, R21, R22, R23, R24, R25, R25b, R26, R26b, R27, R28, R28b, R29, R3, R30, R31, R31b, R3b, R4, R4b, R4c, R5, R5b, R6, R7, R8, R8b, R9, R9b }; rateEquations = { LB1 -> compartment*(-(koffII*II\[LetterSpace]l[t]) + (konII*II\[LetterSpace]f[t]*LIPID[t])/nva), LB10 -> compartment*(-(koffIXa*IXa\[LetterSpace]l[t]) + (konIXa*IXa\[LetterSpace]f[t]*LIPID[t])/nva), LB11 -> compartment*((konX*LIPID[t]*X\[LetterSpace]f[t])/nva - koffX*X\[LetterSpace]l[t]), LB12 -> (konXa*LIPID[t]*Xa\[LetterSpace]f[t])/nva - koffXa*Xa\[LetterSpace]l[t], LB13 -> compartment*(-(koffAPC*APC\[LetterSpace]l[t]) + (konAPC*APC\[LetterSpace]f[t]*LIPID[t])/nva), LB14 -> compartment*((konPS*LIPID[t]*PS\[LetterSpace]f[t])/nva - koffPS*PS\[LetterSpace]l[t]), LB15 -> compartment*((konVIIIai*LIPID[t]*VIIIai\[LetterSpace]f[t])/nva - koffVIIIai*VIIIai\[LetterSpace]l[t]), LB16 -> compartment*((konVai*LIPID[t]*Vai\[LetterSpace]f[t])/nva - koffVai*Vai\[LetterSpace]l[t]), LB17 -> compartment*((konPC*LIPID[t]*PC\[LetterSpace]f[t])/nva - koffPC*PC\[LetterSpace]l[t]), LB2 -> compartment*((konmIIa*LIPID[t]*mIIa\[LetterSpace]f[t])/nva - koffmIIa*mIIa\[LetterSpace]l[t]), LB3 -> compartment*((konV*LIPID[t]*V\[LetterSpace]f[t])/nva - koffV*V\[LetterSpace]l[t]), LB4 -> compartment*((konVa*LIPID[t]*Va\[LetterSpace]f[t])/nva - koffVa*Va\[LetterSpace]l[t]), LB5 -> compartment*((konVII*LIPID[t]*VII\[LetterSpace]f[t])/nva - koffVII*VII\[LetterSpace]l[t]), LB6 -> compartment*((konVIIa*LIPID[t]*VIIa\[LetterSpace]f[t])/nva - koffVIIa*VIIa\[LetterSpace]l[t]), LB7 -> compartment*((konVIII*LIPID[t]*VIII\[LetterSpace]f[t])/nva - koffVIII*VIII\[LetterSpace]l[t]), LB8 -> compartment*((konVIIIa*LIPID[t]*VIIIa\[LetterSpace]f[t])/nva - koffVIIIa*VIIIa\[LetterSpace]l[t]), LB9 -> compartment*(-(koffIX*IX\[LetterSpace]l[t]) + (konIX*IX\[LetterSpace]f[t]*LIPID[t])/nva), R1 -> compartment*(-(k2*TF\[LetterSpace]VIIa\[LetterSpace]l[t]) + k1*TF\[LetterSpace]l[t]*VIIa\[LetterSpace]l[t]), R10 -> -(k25*VIII\[LetterSpace]Xa\[LetterSpace]l[t]) + k24*VIII\[LetterSpace]l[t]*Xa\[LetterSpace]l[t], R10b -> compartment*k26*VIII\[LetterSpace]Xa\[LetterSpace]l[t], R11 -> compartment*(-(k28*V\[LetterSpace]IIa\[LetterSpace]l[t]) + k27*IIa\[LetterSpace]f[t]*V\[LetterSpace]l[t]), R11b -> compartment*k29*V\[LetterSpace]IIa\[LetterSpace]l[t], R12 -> compartment*(-(k31*VIII\[LetterSpace]IIa\[LetterSpace]l[t]) + k30*IIa\[LetterSpace]f[t]*VIII\[LetterSpace]l[t]), R12b -> compartment*k32*VIII\[LetterSpace]IIa\[LetterSpace]l[t], R13 -> compartment*(-(k34*Xa\[LetterSpace]Va\[LetterSpace]II\[LetterSpace]l[t]) + k33*II\[LetterSpace]l[t]*Xa\[LetterSpace]Va\[LetterSpace]l[t]), R14 -> compartment*(k35*mIIa\[LetterSpace]l[t]*Xa\[LetterSpace]Va\[LetterSpace]l[t] - k36*Xa\[LetterSpace]Va\[LetterSpace]mIIa\[LetterSpace]l[t]), R15 -> compartment*k37*Xa\[LetterSpace]Va\[LetterSpace]II\[LetterSpace]l[t], R15b -> compartment*k38*Xa\[LetterSpace]Va\[LetterSpace]mIIa\[LetterSpace]l[t], R16 -> compartment*(-(k40*VII\[LetterSpace]Xa\[LetterSpace]l[t]) + k39*VII\[LetterSpace]l[t]*Xa\[LetterSpace]l[t]), R16b -> compartment*k41*VII\[LetterSpace]Xa\[LetterSpace]l[t], R17 -> compartment*(k42*IIa\[LetterSpace]f[t]*XI\[LetterSpace]f[t] - k43*XI\[LetterSpace]IIa\[LetterSpace]l[t]), R17b -> compartment*k44*XI\[LetterSpace]IIa\[LetterSpace]l[t], R18 -> compartment*(-(k46*APC\[LetterSpace]PS\[LetterSpace]VIIIa\[LetterSpace]l[t]) + k45*APC\[LetterSpace]PS\[LetterSpace]l[t]*VIIIa\[LetterSpace]l[t]), R18b -> compartment*k47*APC\[LetterSpace]PS\[LetterSpace]VIIIa\[LetterSpace]l[t], R19 -> compartment*(-(k49*APC\[LetterSpace]PS\[LetterSpace]Va\[LetterSpace]l[t]) + k48*APC\[LetterSpace]PS\[LetterSpace]l[t]*Va\[LetterSpace]l[t]), R19b -> compartment*k50*APC\[LetterSpace]PS\[LetterSpace]Va\[LetterSpace]l[t], R2 -> compartment*(-(k4*TF\[LetterSpace]VII\[LetterSpace]l[t]) + k3*TF\[LetterSpace]l[t]*VII\[LetterSpace]l[t]), R20 -> compartment*(-(k52*TFPI\[LetterSpace]Xa\[LetterSpace]l[t]) + k51*TFPI\[LetterSpace]f[t]*Xa\[LetterSpace]f[t]), R21 -> compartment*(-(k54*TFPI\[LetterSpace]Xa\[LetterSpace]TF\[LetterSpace]VIIa\[LetterSpace]l[t]) + k53*TFPI\[LetterSpace]Xa\[LetterSpace]l[t]*TF\[LetterSpace]VIIa\[LetterSpace]l[t]), R22 -> compartment*k55*AT\[LetterSpace]f[t]*IXa\[LetterSpace]f[t], R23 -> compartment*k56*AT\[LetterSpace]f[t]*Xa\[LetterSpace]f[t], R24 -> compartment*k57*AT\[LetterSpace]f[t]*IIa\[LetterSpace]f[t], R25 -> compartment*(k58*mIIa\[LetterSpace]l[t]*V\[LetterSpace]l[t] - k59*V\[LetterSpace]mIIa\[LetterSpace]l[t]), R25b -> compartment*k60*V\[LetterSpace]mIIa\[LetterSpace]l[t], R26 -> compartment*(k61*mIIa\[LetterSpace]l[t]*VIII\[LetterSpace]l[t] - k62*VIII\[LetterSpace]mIIa\[LetterSpace]l[t]), R26b -> compartment*k63*VIII\[LetterSpace]mIIa\[LetterSpace]l[t], R27 -> compartment*(-(k65*IIa\[LetterSpace]TM\[LetterSpace]l[t]) + k64*IIa\[LetterSpace]f[t]*TM\[LetterSpace]l[t]), R28 -> compartment*(-(k67*IIa\[LetterSpace]TM\[LetterSpace]PC\[LetterSpace]l[t]) + k66*IIa\[LetterSpace]TM\[LetterSpace]l[t]*PC\[LetterSpace]l[t]), R28b -> compartment*k68*IIa\[LetterSpace]TM\[LetterSpace]PC\[LetterSpace]l[t], R29 -> compartment*k69*AT\[LetterSpace]f[t]*mIIa\[LetterSpace]f[t], R3 -> compartment*(-(k6*TF\[LetterSpace]VIIa\[LetterSpace]IX\[LetterSpace]l[t]) + k5*IX\[LetterSpace]l[t]*TF\[LetterSpace]VIIa\[LetterSpace]l[t]), R30 -> compartment*(-(k71*APC\[LetterSpace]PS\[LetterSpace]l[t]) + k70*APC\[LetterSpace]l[t]*PS\[LetterSpace]l[t]), R31 -> compartment*(-(k73*XIa\[LetterSpace]IX\[LetterSpace]l[t]) + k72*IX\[LetterSpace]l[t]*XIa\[LetterSpace]l[t]), R31b -> compartment*k74*XIa\[LetterSpace]IX\[LetterSpace]l[t], R3b -> compartment*k7*TF\[LetterSpace]VIIa\[LetterSpace]IX\[LetterSpace]l[t], R4 -> compartment*(-(k9*TF\[LetterSpace]VIIa\[LetterSpace]X\[LetterSpace]l[t]) + k8*TF\[LetterSpace]VIIa\[LetterSpace]l[t]*X\[LetterSpace]l[t]), R4b -> compartment*k10*TF\[LetterSpace]VIIa\[LetterSpace]X\[LetterSpace]l[t], R4c -> compartment*k75*TF\[LetterSpace]VIIa\[LetterSpace]Xa\[LetterSpace]l[t], R5 -> compartment*(-(k12*TF\[LetterSpace]VII\[LetterSpace]Xa\[LetterSpace]l[t]) + k11*TF\[LetterSpace]VII\[LetterSpace]l[t]*Xa\[LetterSpace]l[t]), R5b -> compartment*k13*TF\[LetterSpace]VII\[LetterSpace]Xa\[LetterSpace]l[t], R6 -> compartment*(-(k15*IXa\[LetterSpace]VIIIa\[LetterSpace]l[t]) + k14*IXa\[LetterSpace]l[t]*VIIIa\[LetterSpace]l[t]), R7 -> compartment*(k16*Va\[LetterSpace]l[t]*Xa\[LetterSpace]l[t] - k17*Xa\[LetterSpace]Va\[LetterSpace]l[t]), R8 -> compartment*(-(k19*IXa\[LetterSpace]VIIIa\[LetterSpace]X\[LetterSpace]l[t]) + k18*IXa\[LetterSpace]VIIIa\[LetterSpace]l[t]*X\[LetterSpace]l[t]), R8b -> compartment*k20*IXa\[LetterSpace]VIIIa\[LetterSpace]X\[LetterSpace]l[t], R9 -> compartment*(-(k22*V\[LetterSpace]Xa\[LetterSpace]l[t]) + k21*V\[LetterSpace]l[t]*Xa\[LetterSpace]l[t]), R9b -> compartment*k23*V\[LetterSpace]Xa\[LetterSpace]l[t] }; parameters = { k1 -> 0.5, k10 -> 1.5, k11 -> 0.05, k12 -> 44.8, k13 -> 15.2, k14 -> 0.1, k15 -> 0.2, k16 -> 1.0, k17 -> 1.0, k18 -> 0.1, k19 -> 10.7, k2 -> 0.005, k20 -> 8.3, k21 -> 0.1, k22 -> 1.0, k23 -> 0.043, k24 -> 0.1, k25 -> 2.1, k26 -> 0.023, k27 -> 0.1, k28 -> 6.94, k29 -> 0.23, k3 -> 0.005, k30 -> 0.1, k31 -> 13.8, k32 -> 0.9, k33 -> 0.1, k34 -> 100.0, k35 -> 0.1, k36 -> 66.0, k37 -> 13.0, k38 -> 15.0, k39 -> 0.05, k4 -> 0.005, k40 -> 44.8, k41 -> 15.2, k42 -> 0.1, k43 -> 10.0, k44 -> 1.43, k45 -> 0.1, k46 -> 1.6, k47 -> 0.4, k48 -> 0.1, k49 -> 1.6, k5 -> 0.01, k50 -> 0.4, k51 -> 0.016, k52 -> 0.00033, k53 -> 0.01, k54 -> 0.0011, k55 -> 4.9*^-07, k56 -> 2.3*^-06, k57 -> 6.83*^-05, k58 -> 0.1, k59 -> 6.94, k6 -> 2.09, k60 -> 1.035, k61 -> 0.1, k62 -> 13.8, k63 -> 0.9, k64 -> 1.0, k65 -> 0.5, k66 -> 0.1, k67 -> 6.4, k68 -> 3.6, k69 -> 6.83*^-06, k7 -> 0.34, k70 -> 0.1, k71 -> 0.5, k72 -> 0.01, k73 -> 1.417, k74 -> 0.183, k75 -> 1.0, k8 -> 0.1, k9 -> 32.5, koffAPC -> 3.5, koffII -> 1.0, koffIX -> 0.115, koffIXa -> 0.115, koffPC -> 11.5, koffPS -> 0.2, koffV -> 0.145, koffVII -> 0.66, koffVIII -> 0.1, koffVIIIa -> 0.335, koffVIIIai -> 0.335, koffVIIa -> 0.227, koffVa -> 0.17, koffVai -> 0.17, koffX -> 1.9, koffXa -> 3.3, koffmIIa -> 0.475, konAPC -> 0.05, konII -> 0.0043, konIX -> 0.05, konIXa -> 0.05, konPC -> 0.05, konPS -> 0.05, konV -> 0.05, konVII -> 0.05, konVIII -> 0.05, konVIIIa -> 0.05, konVIIIai -> 0.05, konVIIa -> 0.05, konVa -> 0.057, konVai -> 0.057, konX -> 0.01, konXa -> 0.029, konmIIa -> 0.05, nva -> 100.0, compartment -> 1.0 }; assignments = { }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { APC\[LetterSpace]PS\[LetterSpace]VIIIa\[LetterSpace]l'[t] == 1.0*R18 -1.0*R18b, APC\[LetterSpace]PS\[LetterSpace]Va\[LetterSpace]l'[t] == 1.0*R19 -1.0*R19b, APC\[LetterSpace]PS\[LetterSpace]l'[t] == 1.0*R18b +1.0*R19b +1.0*R30 -1.0*R18 -1.0*R19, APC\[LetterSpace]f'[t] == -1.0*LB13, APC\[LetterSpace]l'[t] == 1.0*LB13 +1.0*R28b -1.0*R30, AT\[LetterSpace]f'[t] == -1.0*R22 -1.0*R23 -1.0*R24 -1.0*R29, II\[LetterSpace]f'[t] == -1.0*LB1, II\[LetterSpace]l'[t] == 1.0*LB1 -1.0*R13, IIa\[LetterSpace]AT\[LetterSpace]f'[t] == 1.0*R24 , IIa\[LetterSpace]TM\[LetterSpace]PC\[LetterSpace]l'[t] == 1.0*R28 -1.0*R28b, IIa\[LetterSpace]TM\[LetterSpace]l'[t] == 1.0*R27 +1.0*R28b -1.0*R28, IIa\[LetterSpace]f'[t] == 1.0*R11b +1.0*R12b +1.0*R15b +1.0*R17b -1.0*R11 -1.0*R12 -1.0*R17 -1.0*R24 -1.0*R27, IX\[LetterSpace]f'[t] == -1.0*LB9, IX\[LetterSpace]l'[t] == 1.0*LB9 -1.0*R3 -1.0*R31, IXa\[LetterSpace]AT\[LetterSpace]f'[t] == 1.0*R22 , IXa\[LetterSpace]VIIIa\[LetterSpace]X\[LetterSpace]l'[t] == 1.0*R8 -1.0*R8b, IXa\[LetterSpace]VIIIa\[LetterSpace]l'[t] == 1.0*R6 +1.0*R8b -1.0*R8, IXa\[LetterSpace]f'[t] == -1.0*LB10 -1.0*R22, IXa\[LetterSpace]l'[t] == 1.0*LB10 +1.0*R3b +1.0*R31b -1.0*R6, LIPID'[t] == -100.0*LB1 -100.0*LB2 -100.0*LB3 -100.0*LB4 -100.0*LB5 -100.0*LB6 -100.0*LB7 -100.0*LB8 -100.0*LB9 -100.0*LB10 -100.0*LB11 -100.0*LB12 -100.0*LB13 -100.0*LB14 -100.0*LB15 -100.0*LB16 -100.0*LB17, PC\[LetterSpace]f'[t] == -1.0*LB17, PC\[LetterSpace]l'[t] == 1.0*LB17 -1.0*R28, PS\[LetterSpace]f'[t] == -1.0*LB14, PS\[LetterSpace]l'[t] == 1.0*LB14 -1.0*R30, TFPI\[LetterSpace]Xa\[LetterSpace]TF\[LetterSpace]VIIa\[LetterSpace]l'[t] == 1.0*R21 , TFPI\[LetterSpace]Xa\[LetterSpace]l'[t] == 1.0*R20 -1.0*R21, TFPI\[LetterSpace]f'[t] == -1.0*R20, TF\[LetterSpace]VII\[LetterSpace]Xa\[LetterSpace]l'[t] == 1.0*R5 -1.0*R5b, TF\[LetterSpace]VII\[LetterSpace]l'[t] == 1.0*R2 -1.0*R5, TF\[LetterSpace]VIIa\[LetterSpace]IX\[LetterSpace]l'[t] == 1.0*R3 -1.0*R3b, TF\[LetterSpace]VIIa\[LetterSpace]IXa\[LetterSpace]l'[t] == 0.0 , TF\[LetterSpace]VIIa\[LetterSpace]X\[LetterSpace]l'[t] == 1.0*R4 -1.0*R4b, TF\[LetterSpace]VIIa\[LetterSpace]Xa\[LetterSpace]l'[t] == 1.0*R4b -1.0*R4c, TF\[LetterSpace]VIIa\[LetterSpace]l'[t] == 1.0*R1 +1.0*R3b +1.0*R4c +1.0*R5b -1.0*R3 -1.0*R4 -1.0*R21, TF\[LetterSpace]l'[t] == -1.0*R1 -1.0*R2, TM\[LetterSpace]l'[t] == -1.0*R27, VIII\[LetterSpace]IIa\[LetterSpace]l'[t] == 1.0*R12 -1.0*R12b, VIII\[LetterSpace]Xa\[LetterSpace]l'[t] == 1.0*R10 -1.0*R10b, VIII\[LetterSpace]f'[t] == -1.0*LB7, VIII\[LetterSpace]l'[t] == 1.0*LB7 -1.0*R10 -1.0*R12 -1.0*R26, VIII\[LetterSpace]mIIa\[LetterSpace]l'[t] == 1.0*R26 -1.0*R26b, VIIIa\[LetterSpace]f'[t] == -1.0*LB8, VIIIa\[LetterSpace]l'[t] == 1.0*LB8 +1.0*R10b +1.0*R12b +1.0*R26b -1.0*R6 -1.0*R18, VIIIai\[LetterSpace]f'[t] == -1.0*LB15, VIIIai\[LetterSpace]l'[t] == 1.0*LB15 +1.0*R18b , VII\[LetterSpace]Xa\[LetterSpace]l'[t] == 1.0*R16 -1.0*R16b, VII\[LetterSpace]f'[t] == -1.0*LB5, VII\[LetterSpace]l'[t] == 1.0*LB5 -1.0*R2 -1.0*R16, VIIa\[LetterSpace]f'[t] == -1.0*LB6, VIIa\[LetterSpace]l'[t] == 1.0*LB6 +1.0*R16b -1.0*R1, V\[LetterSpace]IIa\[LetterSpace]l'[t] == 1.0*R11 -1.0*R11b, V\[LetterSpace]Xa\[LetterSpace]l'[t] == 1.0*R9 -1.0*R9b, V\[LetterSpace]f'[t] == -1.0*LB3, V\[LetterSpace]l'[t] == 1.0*LB3 -1.0*R9 -1.0*R11 -1.0*R25, V\[LetterSpace]mIIa\[LetterSpace]l'[t] == 1.0*R25 -1.0*R25b, Va\[LetterSpace]f'[t] == -1.0*LB4, Va\[LetterSpace]l'[t] == 1.0*LB4 +1.0*R9b +1.0*R11b +1.0*R25b -1.0*R7 -1.0*R19, Vai\[LetterSpace]f'[t] == -1.0*LB16, Vai\[LetterSpace]l'[t] == 1.0*LB16 +1.0*R19b , XI\[LetterSpace]IIa\[LetterSpace]l'[t] == 1.0*R17 -1.0*R17b, XI\[LetterSpace]f'[t] == -1.0*R17, XIa\[LetterSpace]IX\[LetterSpace]l'[t] == 1.0*R31 -1.0*R31b, XIa\[LetterSpace]l'[t] == 1.0*R17b +1.0*R31b -1.0*R31, X\[LetterSpace]f'[t] == -1.0*LB11, X\[LetterSpace]l'[t] == 1.0*LB11 -1.0*R4 -1.0*R8, Xa\[LetterSpace]AT\[LetterSpace]f'[t] == 1.0*R23 , Xa\[LetterSpace]Va\[LetterSpace]II\[LetterSpace]l'[t] == 1.0*R13 -1.0*R15, Xa\[LetterSpace]Va\[LetterSpace]l'[t] == 1.0*R7 +1.0*R15b -1.0*R13 -1.0*R14, Xa\[LetterSpace]Va\[LetterSpace]mIIa\[LetterSpace]l'[t] == 1.0*R14 +1.0*R15 -1.0*R15b, Xa\[LetterSpace]f'[t] == -1.0*LB12 -1.0*R20 -1.0*R23, Xa\[LetterSpace]l'[t] == 1.0*LB12 +1.0*R4c +1.0*R5b +1.0*R8b +1.0*R9b +1.0*R10b +1.0*R16b -1.0*R5 -1.0*R7 -1.0*R9 -1.0*R10 -1.0*R16, mIIa\[LetterSpace]AT\[LetterSpace]l'[t] == 1.0*R29 , mIIa\[LetterSpace]f'[t] == -1.0*LB2 -1.0*R29, mIIa\[LetterSpace]l'[t] == 1.0*LB2 +1.0*R25b +1.0*R26b -1.0*R14 -1.0*R25 -1.0*R26 }; timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}]; (* Steady-state solution initialized with result of time evolution *) findRootEquations = odes /.D[_[t],t]->0; findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2]; steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100] fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations (* Plot the time evolution of the variables *) plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]