// Created by libAntimony v2.8.0
function activitycalc(x, slopea, yintercepta)
yintercepta + x*slopea;
end
model *BIOMD0000000436()
// Compartments and Species:
compartment c1;
species LPS in c1, DG in c1;
substanceOnly species $FA in c1, HETE in c1, PGH2 in c1, PGE2 in c1, PGF2a in c1;
substanceOnly species PGD2 in c1, PGJ2 in c1, dPGJ2 in c1, AA in c1, GPCho in c1;
substanceOnly species dPGD2 in c1;
// Assignment Rules:
timevalue := time + t0;
LPSactivity := activitycalc(timevalue, LPSslope, LPSint);
DGactivity := activitycalc(timevalue, DGslope, DGint);
GPChoratio := activitycalc(timevalue, GPChoslope, GPChoint)/GPChointzero;
// Reactions:
r2: $FA => AA; k2*FA;
r3: DG => AA; k3*DGactivity*DGperDNA*DNA_;
r7: GPCho => AA; k7*GPChoratio*GPCho;
r5: GPCho => AA; (k5*DGactivity*GPChoratio*GPCho)/DNA_;
r11: AA => PGH2; k11*AA;
r10: AA => PGH2; (k10*DGactivity*AA)/DNA_;
r8: AA => HETE; k8*AA;
r4: AA => ; k4*AA;
r9: HETE => ; k9*HETE;
r13: PGH2 => PGE2; k13*PGH2;
r14: PGE2 => ; k14*PGE2;
r15: PGH2 => PGF2a; k15*PGH2;
r16: PGF2a => ; k16*PGF2a;
r17: PGH2 => PGD2; k17*PGH2;
r18: PGD2 => PGJ2; k18*PGD2;
r19: PGD2 => dPGD2; k19*PGD2;
r20: dPGD2 => ; k20*dPGD2;
r21: PGJ2 => dPGJ2; k21*PGJ2;
r22: dPGJ2 => ; k22*dPGJ2;
r1: $FA => AA; k1*onepmol*LPSactivity;
r6: GPCho => AA; k6*GPCho*GPChoratio*LPSactivity;
r12: AA => PGH2; k12*AA*LPSactivity;
// Events:
_E0: at timevalue > t0: GPChoslope = GPChoslopezero, GPChoint = GPChointzero, DGslope = DGslopezero, DGint = DGintzero, LPSslope = LPSslopezero, LPSint = LPSintzero;
_E1: at timevalue >= t1: GPChoslope = GPChoslope1, GPChoint = GPChoint1, DGslope = DGslope1, DGint = DGint1, LPSslope = LPSslope1, LPSint = LPSint1;
_E2: at timevalue >= t2: GPChoslope = GPChoslope2, GPChoint = GPChoint2, DGslope = DGslope2, DGint = DGint2;
_E3: at timevalue >= t3: GPChoslope = GPChoslope3, GPChoint = GPChoint3, DGslope = DGslope3, DGint = DGint3, LPSslope = LPSslope3, LPSint = LPSint3;
_E4: at timevalue >= t4: GPChoslope = GPChoslope4, GPChoint = GPChoint4, DGslope = DGslope4, DGint = DGint4;
_E5: at timevalue >= t5: GPChoslope = GPChoslope5, GPChoint = GPChoint5, DGslope = DGslope5, DGint = DGint5;
_E6: at timevalue >= t6: GPChoslope = GPChoslope6, GPChoint = GPChoint6, DGslope = DGslope6, DGint = DGint6;
// Species initializations:
FA = 1/c1;
FA has pmol_per_litre;
HETE = 0;
HETE has pmol_per_litre;
PGH2 = 0;
PGH2 has pmol_per_litre;
PGE2 = 0;
PGE2 has pmol_per_litre;
PGF2a = 0;
PGF2a has pmol_per_litre;
PGD2 = 0;
PGJ2 = 0;
PGJ2 has pmol_per_litre;
dPGJ2 = 0;
dPGJ2 has pmol_per_litre;
AA = 25/c1;
AA has pmol_per_litre;
LPS = 0.5/c1;
LPS has dimensionless_per_litre;
DG = 0;
DG has dimensionless_per_litre;
GPCho = 1/c1;
GPCho has pmol_per_litre;
dPGD2 = 0;
dPGD2 has pmol_per_litre;
// Compartment initializations:
c1 = 1;
c1 has litre;
timevalue has time_unit;
// Variable initializations:
t0 = 0;
t0 has time_unit;
LPSint = 0;
LPSint has dimensionless;
LPSintzero = 0;
LPSintzero has dimensionless;
LPSslope = 0;
LPSslope has perh;
LPSslopezero = 2;
LPSslopezero has perh;
DGint = 1;
DGint has dimensionless;
DGintzero = 1;
DGintzero has dimensionless;
DGslope = 1;
DGslope has perh;
DGslopezero = 0.3014;
DGslopezero has perh;
GPChoint = 315.3137;
GPChoint has pmol;
GPChointzero = 315.3137;
GPChointzero has pmol;
GPChoslope = -158.43;
GPChoslope has pmolperh;
GPChoslopezero = -158.43;
GPChoslopezero has pmolperh;
t1 = 0.5;
t1 has time_unit;
LPSint1 = 1.3335;
LPSint1 has dimensionless;
LPSslope1 = -0.6667;
LPSslope1 has perh;
DGint1 = 1.2292;
DGint1 has dimensionless;
DGslope1 = -0.1572;
DGslope1 has perh;
GPChoint1 = 179.96;
GPChoint1 has pmol;
GPChoslope1 = 112.28;
GPChoslope1 has pmolperh;
t2 = 1;
t2 has time_unit;
DGint2 = 1.0282;
DGint2 has dimensionless;
DGslope2 = 0.0439;
DGslope2 has perh;
GPChoint2 = 267.93;
GPChoint2 has pmol;
GPChoslope2 = 24.31;
GPChoslope2 has pmolperh;
t3 = 2;
t3 has time_unit;
LPSint3 = 0;
LPSint3 has dimensionless;
LPSslope3 = 0;
LPSslope3 has perh;
DGint3 = 1.2987;
DGint3 has dimensionless;
DGslope3 = -0.0913;
DGslope3 has perh;
GPChoint3 = 396.41;
GPChoint3 has pmol;
GPChoslope3 = -39.93;
GPChoslope3 has pmolperh;
t4 = 4;
t4 has time_unit;
DGint4 = 0.543;
DGint4 has dimensionless;
DGslope4 = 0.0976;
DGslope4 has perh;
GPChoint4 = 214.77;
GPChoint4 has pmol;
GPChoslope4 = 5.48;
GPChoslope4 has pmolperh;
t5 = 8;
t5 has time_unit;
DGint5 = 0.7549;
DGint5 has dimensionless;
DGslope5 = 0.0711;
DGslope5 has perh;
GPChoint5 = 41.93;
GPChoint5 has pmol;
GPChoslope5 = 27.08;
GPChoslope5 has pmolperh;
t6 = 12;
t6 has time_unit;
DGint6 = 0.5927;
DGint6 has dimensionless;
DGslope6 = 0.0846;
DGslope6 has perh;
GPChoint6 = 425.43;
GPChoint6 has pmol;
GPChoslope6 = -4.874;
GPChoslope6 has pmolperh;
k1 = 355.637;
k1 has perh;
k2 = 1e-015;
k2 has perh;
k3 = 1e-015;
k3 has perh;
k4 = 1e-015;
k4 has perh;
k5 = 1e-015;
k5 has ugDNAperh;
k6 = 0.33;
k6 has perh;
k7 = 1e-015;
k7 has perh;
k8 = 0.007;
k8 has perh;
k9 = 0.187;
k9 has perh;
k10 = 0.024;
k10 has ugDNAperh;
k11 = 0.111;
k11 has perh;
k12 = 0.098;
k12 has perh;
k13 = 0.204;
k13 has perh;
k14 = 1e-015;
k14 has perh;
k15 = 0.061;
k15 has perh;
k16 = 1e-015;
k16 has perh;
k17 = 3.116;
k17 has perh;
k18 = 0.054;
k18 has perh;
k19 = 0.029;
k19 has perh;
k20 = 0.014;
k20 has perh;
k21 = 0.034;
k21 has perh;
k22 = 0.116;
k22 has perh;
DNA_ = 1;
DNA_ has ug;
DGperDNA = 1;
DGperDNA has pmolperugDNA;
onepmol = 1;
onepmol has pmol;
zero = 0;
zero has dimensionless;
LPSactivity has dimensionless;
DGactivity has dimensionless;
GPChoratio has dimensionless;
// Other declarations:
var timevalue, LPSint, LPSslope, DGint, DGslope, GPChoint, GPChoslope, DGperDNA;
var LPSactivity, DGactivity, GPChoratio;
const c1, t0, LPSintzero, LPSslopezero, DGintzero, DGslopezero, GPChointzero;
const GPChoslopezero, t1, LPSint1, LPSslope1, DGint1, DGslope1, GPChoint1;
const GPChoslope1, t2, DGint2, DGslope2, GPChoint2, GPChoslope2, t3, LPSint3;
const LPSslope3, DGint3, DGslope3, GPChoint3, GPChoslope3, t4, DGint4, DGslope4;
const GPChoint4, GPChoslope4, t5, DGint5, DGslope5, GPChoint5, GPChoslope5;
const t6, DGint6, DGslope6, GPChoint6, GPChoslope6, k1, k2, k3, k4, k5;
const k6, k7, k8, k9, k10, k11, k12, k13, k14, k15, k16, k17, k18, k19;
const k20, k21, k22, DNA_, onepmol, zero;
// Unit definitions:
unit time_unit = 3600 second;
unit perh = 1 / 0.0002778 second;
unit pmolperugDNA = 1e-12 mole / 1e-6 gram;
unit ug = 1e-6 gram;
unit pmol = 1e-12 mole;
unit pmolperh = 1e-12 mole / 0.0002778 second;
unit ugDNAperh = 1e-6 gram / 3600 second;
unit pmol_per_litre = 1e-12 mole / litre;
unit dimensionless_per_litre = dimensionless / litre;
// Display Names:
AA is "Arachidonic acid";
end