begin name
rehm1
end name

begin reactions
Apoformation {1.0}PC9 + {1.0}APAF1 = {1.0}C9
CytCrelease {1.0}CytCmit = {1.0}CytCcell
SMACrelease {1.0}SMACmito = {1.0}SMAC
reaction1 {1.0}$EXTERNAL = {1.0}PC3
reaction10 {1.0}C3 + {1.0}BIR12 = {1.0}BIR12C3
reaction11 {1.0}C3 + {1.0}XIAP = {1.0}BIR12 + {1.0}BIR3R + {1.0}C3
reaction12 {1.0}C3 + {1.0}XIAPC9 = {1.0}BIR12 + {1.0}BIR3RC9 + {1.0}C3
reaction13 {1.0}C3 + {1.0}XIAPC3 = {1.0}C3 + {1.0}BIR12C3 + {1.0}BIR3R
reaction14 {1.0}C3 + {1.0}XIAPp2frag = {1.0}C3 + {1.0}BIR12 + {1.0}BIR3Rp2frag
reaction15 {1.0}C3 + {1.0}XIAPp2fragC3 = {1.0}C3 + {1.0}BIR12C3 + {1.0}BIR3Rp2frag
reaction16 {1.0}C3 + {1.0}XIAPC9C3 = {1.0}C3 + {1.0}BIR12C3 + {1.0}BIR3RC9
reaction17 {1.0}C3 + {1.0}XIAP2SMAC = {1.0}C3 + {1.0}BIR12SMAC + {1.0}BIR3RSMAC
reaction18 {1.0}C3 + {1.0}XIAPC9C3 = {1.0}C3 + {1.0}C9P + {1.0}XIAPp2fragC3
reaction19 {1.0}C3 + {1.0}XIAPC9 = {1.0}C3 + {1.0}C9P + {1.0}XIAPp2frag
reaction2 {1.0}$EXTERNAL = {1.0}XIAP
reaction20 {1.0}C3 + {1.0}BIR3RC9 = {1.0}C3 + {1.0}BIR3Rp2frag + {1.0}C9P
reaction21 {1.0}C9 + {1.0}XIAP = {1.0}XIAPC9
reaction22 {1.0}C9 + {1.0}XIAPC3 = {1.0}XIAPC9C3
reaction23 {1.0}C9 + {1.0}BIR3R = {1.0}BIR3RC9
reaction24 {1.0}BIR3Rp2frag = {1.0}BIR3R
reaction25 {1.0}XIAPp2frag = {1.0}XIAP
reaction26 {1.0}XIAP + {2.0}SMAC = {1.0}XIAP2SMAC
reaction27 {1.0}XIAPC9 + {2.0}SMAC = {1.0}XIAP2SMAC + {1.0}C9
reaction28 {1.0}XIAPC3 + {2.0}SMAC = {1.0}XIAP2SMAC + {1.0}C3
reaction29 {1.0}XIAPC9C3 + {2.0}SMAC = {1.0}XIAP2SMAC + {1.0}C3 + {1.0}C9
reaction3 {1.0}C9 + {1.0}PC3 = {1.0}C3 + {1.0}C9
reaction30 {1.0}BIR12 + {1.0}SMAC = {1.0}BIR12SMAC
reaction31 {1.0}BIR3R + {1.0}SMAC = {1.0}BIR3RSMAC
reaction32 {1.0}BIR12C3 + {1.0}SMAC = {1.0}BIR12SMAC + {1.0}C3
reaction33 {1.0}BIR3RC9 + {1.0}SMAC = {1.0}BIR3RSMAC + {1.0}C9
reaction34 {1.0}XIAPp2frag + {2.0}SMAC = {1.0}XIAPp2frag2SMAC
reaction35 {1.0}C9P = {1.0}$EXTERNAL
reaction36 {1.0}C9 = {1.0}$EXTERNAL
reaction37 {1.0}C3 = {1.0}$EXTERNAL
reaction38 {1.0}XIAPC3 = {1.0}$EXTERNAL
reaction39 {1.0}XIAPC9C3 = {1.0}$EXTERNAL
reaction4 {1.0}C9 + {1.0}C3 = {1.0}C9P + {1.0}C3
reaction40 {1.0}XIAPC9 = {1.0}$EXTERNAL
reaction41 {1.0}XIAPp2frag = {1.0}$EXTERNAL
reaction42 {1.0}XIAPp2fragC3 = {1.0}$EXTERNAL
reaction43 {1.0}XIAPp2frag2SMAC = {1.0}$EXTERNAL
reaction44 {1.0}XIAP2SMAC = {1.0}$EXTERNAL
reaction45 {1.0}BIR12 = {1.0}$EXTERNAL
reaction46 {1.0}BIR3R = {1.0}$EXTERNAL
reaction47 {1.0}BIR12SMAC = {1.0}$EXTERNAL
reaction48 {1.0}BIR3RSMAC = {1.0}$EXTERNAL
reaction49 {1.0}BIR12C3 = {1.0}$EXTERNAL
reaction5 {1.0}C9P + {1.0}PC3 = {1.0}C9P + {1.0}C3
reaction50 {1.0}BIR3RC9 = {1.0}$EXTERNAL
reaction51 {1.0}BIR3Rp2frag = {1.0}$EXTERNAL
reaction52 {1.0}SMAC = {1.0}$EXTERNAL
reaction53 {1.0}C3 + {1.0}Substrate = {1.0}C3 + {1.0}ClvgPrds
reaction6 {1.0}PC3 + {1.0}C3 = {2.0}C3
reaction7 {1.0}C3 + {1.0}XIAP = {1.0}XIAPC3
reaction8 {1.0}C3 + {1.0}XIAPC9 = {1.0}XIAPC9C3
reaction9 {1.0}C3 + {1.0}XIAPp2frag = {1.0}XIAPp2fragC3
end reactions

begin rate equations
Apoformation = (apolim*cell*Log[2])/thApop
CytCrelease = (cell*CytCmit[t]*Log[2])/thCytC
SMACrelease = (cell*Log[2]*SMACmito[t])/thSMAC
reaction1 = cell*(k1 - k1r*PC3[t])
reaction10 = cell*(-(k10r*BIR12C3[t]) + k10*BIR12[t]*C3[t])
reaction11 = cell*k11*C3[t]*XIAP[t]
reaction12 = cell*k12*C3[t]*XIAPC9[t]
reaction13 = cell*k13*C3[t]*XIAPC3[t]
reaction14 = cell*k14*C3[t]*XIAPp2frag[t]
reaction15 = cell*k15*C3[t]*XIAPp2fragC3[t]
reaction16 = cell*k16*C3[t]*XIAPC9C3[t]
reaction17 = cell*k17*C3[t]*XIAP2SMAC[t]
reaction18 = cell*k18*C3[t]*XIAPC9C3[t]
reaction19 = cell*k19*C3[t]*XIAPC9[t]
reaction2 = cell*(k2 - k2r*XIAP[t])
reaction20 = cell*k20*BIR3RC9[t]*C3[t]
reaction21 = cell*(k21*C9[t]*XIAP[t] - k21r*XIAPC9[t])
reaction22 = cell*(k22*C9[t]*XIAPC3[t] - k22r*XIAPC9C3[t])
reaction23 = cell*(-(k23r*BIR3RC9[t]) + k23*BIR3R[t]*C9[t])
reaction24 = cell*k24*BIR3Rp2frag[t]
reaction25 = cell*k25*XIAPp2frag[t]
reaction26 = cell*(k26*SMAC[t]^2*XIAP[t] - k26r*XIAP2SMAC[t])
reaction27 = cell*(-(k27r*C9[t]*XIAP2SMAC[t]) + k27*SMAC[t]^2*XIAPC9[t])
reaction28 = cell*(-(k28r*C3[t]*XIAP2SMAC[t]) + k28*SMAC[t]^2*XIAPC3[t])
reaction29 = cell*(-(k29r*C3[t]*C9[t]*XIAP2SMAC[t]) + k29*SMAC[t]^2*XIAPC9C3[t])
reaction3 = cell*k3*C9[t]*PC3[t]
reaction30 = cell*(-(k30r*BIR12SMAC[t]) + k30*BIR12[t]*SMAC[t])
reaction31 = cell*(-(k31r*BIR3RSMAC[t]) + k31*BIR3R[t]*SMAC[t])
reaction32 = cell*(-(k32r*BIR12SMAC[t]*C3[t]) + k32*BIR12C3[t]*SMAC[t])
reaction33 = cell*(-(k33r*BIR3RSMAC[t]*C9[t]) + k33*BIR3RC9[t]*SMAC[t])
reaction34 = cell*(k34*SMAC[t]^2*XIAPp2frag[t] - k34r*XIAPp2frag2SMAC[t])
reaction35 = cell*k35*C9P[t]
reaction36 = cell*k36*C9[t]
reaction37 = cell*k37*C3[t]
reaction38 = cell*k38*XIAPC3[t]
reaction39 = cell*k39*XIAPC9C3[t]
reaction4 = cell*k4*C3[t]*C9[t]
reaction40 = cell*k40*XIAPC9[t]
reaction41 = cell*k41*XIAPp2frag[t]
reaction42 = cell*k42*XIAPp2fragC3[t]
reaction43 = cell*k43*XIAPp2frag2SMAC[t]
reaction44 = cell*k44*XIAP2SMAC[t]
reaction45 = cell*k45*BIR12[t]
reaction46 = cell*k46*BIR3R[t]
reaction47 = cell*k47*BIR12SMAC[t]
reaction48 = cell*k48*BIR3RSMAC[t]
reaction49 = cell*k49*BIR12C3[t]
reaction5 = cell*k5*C9P[t]*PC3[t]
reaction50 = cell*k50*BIR3RC9[t]
reaction51 = cell*k51*BIR3Rp2frag[t]
reaction52 = cell*k52*SMAC[t]
reaction53 = cell*k53*C3[t]*Substrate[t]
reaction6 = cell*k6*C3[t]*PC3[t]
reaction7 = cell*(k7*C3[t]*XIAP[t] - k7r*XIAPC3[t])
reaction8 = cell*(k8*C3[t]*XIAPC9[t] - k8r*XIAPC9C3[t])
reaction9 = cell*(k9*C3[t]*XIAPp2frag[t] - k9r*XIAPp2fragC3[t])
end rate equations

begin parameters
APAF1tot = APAF1
C3tot = PC3
C9tot = PC9
CytCtot = CytCmit
SMACtot = SMACmito
Sbreakdown = 0.01
XIAPini = XIAP
k1 = 0.000468
k10 = 156.0
k10r = 0.144
k11 = 12.0
k12 = 12.0
k13 = 12.0
k14 = 12.0
k15 = 12.0
k16 = 12.0
k17 = 12.0
k18 = 12.0
k19 = 12.0
k1r = 0.0039
k2 = 0.0007308
k20 = 12.0
k21 = 156.0
k21r = 0.144
k22 = 156.0
k22r = 0.144
k23 = 156.0
k23r = 0.144
k24 = 0.0
k25 = 0.0
k26 = 420.0
k26r = 0.133
k27 = 420.0
k27r = 156.0
k28 = 420.0
k28r = 156.0
k29 = 420.0
k29r = 156.0
k2r = 0.0116
k3 = 6.0
k30 = 4.45
k30r = 31.9
k31 = 0.33
k31r = 14.2
k32 = 4.45
k32r = 156.0
k33 = 0.33
k33r = 156.0
k34 = 420.0
k34r = 156.0
k35 = 0.0058
k36 = 0.0058
k37 = 0.0058
k38 = 0.0347
k39 = 0.0347
k4 = 12.0
k40 = 0.0347
k41 = 0.0058
k42 = 0.0347
k43 = 0.0347
k44 = 0.0347
k45 = 0.0058
k46 = 0.0347
k47 = 0.0058
k48 = 0.0347
k49 = 0.0058
k5 = 48.0
k50 = 0.0058
k51 = 0.0347
k52 = 0.0058
k53 = 12.0
k6 = 2.4
k7 = 156.0
k7r = 0.144
k8 = 156.0
k8r = 0.144
k9 = 156.0
k9r = 0.144
thApop = 2.3
thCytC = 1.5
thSMAC = 7.0
cell = 1.0
mito = 1.0
end parameters

begin initial conditions
APAF1[0] = APAF1i
BIR12[0] = BIR12i
BIR12C3[0] = BIR12C3i
BIR12SMAC[0] = BIR12SMACi
BIR3R[0] = BIR3Ri
BIR3RC9[0] = BIR3RC9i
BIR3RSMAC[0] = BIR3RSMACi
BIR3Rp2frag[0] = BIR3Rp2fragi
C3[0] = C3i
C9[0] = C9i
C9P[0] = C9Pi
ClvgPrds[0] = ClvgPrdsi
CytCcell[0] = CytCcelli
CytCmit[0] = CytCmiti
PC3[0] = PC3i
PC9[0] = PC9i
SMAC[0] = SMACi
SMACmito[0] = SMACmitoi
Substrate[0] = Substratei
XIAP[0] = XIAPi
XIAP2SMAC[0] = XIAP2SMACi
XIAPC3[0] = XIAPC3i
XIAPC9[0] = XIAPC9i
XIAPC9C3[0] = XIAPC9C3i
XIAPp2frag[0] = XIAPp2fragi
XIAPp2frag2SMAC[0] = XIAPp2frag2SMACi
XIAPp2fragC3[0] = XIAPp2fragC3i
end initial conditions

begin initial values
APAF1i = 3.372
BIR12i = 0.0
BIR12C3i = 0.0
BIR12SMACi = 0.0
BIR3Ri = 0.0
BIR3RC9i = 0.0
BIR3RSMACi = 0.0
BIR3Rp2fragi = 0.0
C3i = 0.0
C9i = 0.0
C9Pi = 0.0
ClvgPrdsi = 0.0
CytCcelli = 0.0
CytCmiti = 10.0
PC3i = 0.12
PC9i = 0.03
SMACi = 0.0
SMACmitoi = 0.126
Substratei = 1.0
XIAPi = 0.063
XIAP2SMACi = 0.0
XIAPC3i = 0.0
XIAPC9i = 0.0
XIAPC9C3i = 0.0
XIAPp2fragi = 0.0
XIAPp2frag2SMACi = 0.0
XIAPp2fragC3i = 0.0
end initial values

begin assignment rules
BIR3Rp2fragnorm := BIR3Rp2frag[t]/XIAPini
BIR3RC9norm := BIR3RC9[t]/XIAPini
BIR12C3norm := BIR12C3[t]/XIAPini
BIR3Rnorm := BIR3R[t]/XIAPini
BIR12norm := BIR12[t]/XIAPini
XIAPp2fragnorm := (XIAPp2frag[t] + XIAPp2fragC3[t])/XIAPini
XIAPC9norm := (XIAPC9[t] + XIAPC9C3[t])/XIAPini
XIAPC3norm := (XIAPC3[t] + XIAPC9C3[t])/XIAPini
XIAPnorm := XIAP[t]/XIAPini
Apoform := 1 - PC9[t]/C9tot
SMACrel := 1 - SMACmito[t]/SMACtot
PC3norm := PC3[t]/C3tot
C3norm := C3[t]/C3tot
C9Pnorm := C9P[t]/C9tot
C9norm := C9[t]/C9tot
XIAP2SMACnorm := XIAP2SMAC[t]/XIAPini
SMACnorm := SMAC[t]/XIAPini
apolim := Piecewise[{{PC9[t], PC9[t] < APAF1[t]}}, APAF1[t]]
end assignment rules

begin function definitions
end function definitions

begin events
Production_Breakdown = (Substrate[t] < Sbreakdown, {k1 -> 0, k1r -> 0, k2 -> 0, k2r -> 0, k38 -> 0.0058, k39 -> 0.0058, k40 -> 0.0058, k43 -> 0.0058, k44 -> 0.0058, k46 -> 0.0058, k48 -> 0.0058, k50 -> 0.0058, k51 -> 0.0058}, 0)
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units