Skip to main content

Modeling the effects of a Staphylococcal Enterotoxin B (SEB) on the apoptosis pathway



The lack of detailed understanding of the mechanism of action of many biowarfare agents poses an immediate challenge to biodefense efforts. Many potential bioweapons have been shown to affect the cellular pathways controlling apoptosis [14]. For example, pathogen-produced exotoxins such as Staphylococcal Enterotoxin B (SEB) and Anthrax Lethal Factor (LF) have been shown to disrupt the Fas-mediated apoptotic pathway [2, 4]. To evaluate how these agents affect these pathways it is first necessary to understand the dynamics of a normally functioning apoptosis network. This can then serve as a baseline against which a pathogen perturbed system can be compared. Such comparisons can expose both the proteins most susceptible to alteration by the agent as well as the most critical reaction rates to better instill control on a biological network.


We explore this through the modeling and simulation of the Fas-mediated apoptotic pathway under normal and SEB influenced conditions. We stimulated human Jurkat cells with an anti-Fas antibody in the presence and absence of SEB and determined the relative levels of seven proteins involved in the core pathway at five time points following exposure. These levels were used to impute relative rate constants and build a quantitative model consisting of a series of ordinary differential equations (ODEs) that simulate the network under both normal and pathogen-influenced conditions. Experimental results show that cells exposed to SEB exhibit an increase in the rate of executioner caspase expression (and subsequently apoptosis) of 1 hour 43 minutes (± 14 minutes), as compared to cells undergoing normal cell death.


Our model accurately reflects these results and reveals intervention points that can be altered to restore SEB-influenced system dynamics back to levels within the range of normal conditions.


Apoptosis is a naturally occurring mechanism by which a cell undergoes programmed death in response to external or internal signals. It has been shown that many potential biowarfare agents including Staphylococcal Enterotoxin B (SEB) improperly activate the pathways controlling this process in a variety of cell types [14]. These agents target specific components of these pathways, inducing a cascade of reactions that modify the rate at which apoptosis occurs. The large scale cell death that ensues brings about a state of unrecoverable shock.

Since the induction of apoptosis underlies the lethality of SEB, this pathway presents a good target for therapeutic intervention. Previous results have shown that SEB treatment affects the expression of proteins involved in the Fas-mediated apoptotic pathway [2, 4]. However, the mechanisms by which these agents interact with this pathway, and the cellular components most advantageous for potential therapies, is unclear. In order to establish a better understanding of these interactions, it is first necessary to develop a model of the pathway as it functions under normal cellular conditions. This model can then serve as a baseline against which experimentally perturbed systems can be compared and key intermediates determined.

The model of receptor-mediated apoptosis is represented in Figure 1 and summarized as follows. The cascade begins with extracellular Fas ligand (FasL) binding to the cell surface receptor Fas/CD95, thus initializing aggregation and activation [5]. The adapter protein FADD binds to the Fas receptor through its death domain (DD) and transfers the apoptotic signal to the inactive forms of caspase-8 (procaspase-8) and caspase-3 (procaspse-3). This causes procaspase-8 to undergo a self proteolytic cleavage reaction resulting in its activation, as well as small traces of proteolytic cleavage of procaspase-3. Once activated, caspase-8 is instrumental in the subsequent activation of multiple other proteins [6]. Of particular interest is the subsequent proteolytic cleavage of procaspase-3 by caspase-8. Activation of caspase-3 rapidly results in DNA degradation and inactivation of other cellular functions [5]. The activation of executioner caspase-3 is also one of the final stages in cell death, thus providing an endpoint marker for verification of apoptosis.

Figure 1
figure 1

Diagram of Fas-mediated apoptosis. Model pathway representation of Fas-mediated apoptosis (modified from Fussenegger et al., 2000).

Figure 2
figure 2

Maximum rate constant range comparison. Resulting minimum aggregate distances between experimental and simulated data for 150 initial value seedings (Monte Carlo simulation). Each curve represents the aggregated error (summated simulated vs. experimental) for the proteins in the model (y-axis) as a function of the trial (x-axis) using increasing maximum allowable relative rate constants. This allows the relative rate constants in the model to be constrained to the smallest range of values.

A level of control over this process is provided by a family of proteins called inhibitors of apoptosis (IAPs). IAPs can intervene at the executioner caspase point to suppress cell death [79]. Leveraging this regulatory reaction can increase our understanding of the system dynamics to instil regulation points and control mechanisms.

Our main objectives within this exposition are two-fold. We developed a parsimonious model representation of Fas-mediated apoptosis under normal cellular conditions in an attempt to (1) reduce the complexity of the biological network to the most critical protein intermediates and (2) elucidate the most influential reaction rates for controlling an SEB-exposed system and reverting the system back to a normal cellular state. Our model [see Additional file 1] provides an accurate representation of the pathway dynamics, as evaluated by consistency with experimental data.


Model of the normally functioning Fas pathway

In order to develop a model of Fas-mediated apoptosis under what we identify as 'normal' conditions, we stimulated human Jurkat cells with an anti-Fas antibody and determined the relative levels of seven proteins involved in the core pathway by western blot at 0, 2, 4, 8, and 12 hours after exposure. It should be noted that the introduction of a signal (the anti-Fas antibody in this case) for induction of apoptosis is artificial as compared to completely normal cellular conditions. However the primary comparison seeks to identify differences between non-pathogen induced and pathogen induced apoptosis. The measured protein levels were used to impute relative rate constants and build a quantitative model consisting of a series of ordinary differential equations (ODEs) that simulate the network under normal conditions. Flow cytometry was also conducted on the stimulated Jurkat cells at the time points specified to verify the presence of cell death and determine a difference in early apoptotic events between a normally functioning Fas pathway and one under the influence of SEB exposure (examined in the next section).

Under these conditions, the simulated results are in agreement with the experimental data. All proteins assayed, with exception to procaspase-8, correlate with data in our model with coefficient values of at least 0.70 (Table 1). From Figure 3a, the Fas receptor is shown to immediately decrease in level after activation with the anti-Fas activating antibody, while FADD increases with a strong inverse relationship (r = -0.72 between FADD and the Fas receptor). The simulated results from the inactive and active caspase-3 proteins reach successive peaks, with a time difference of 1.71 hours (Figure 3c), suggestive of the time required for proteolytic cleavage of the precursor form and subsequent activation. The mean protein levels from the experimental data follow the same pattern, with a time difference of 1.96 ± 1.63 hours between inactive and active caspase-3. This difference between average peaks in inactive and active caspase-3 for experimental and simulated results amounts to approximately 15 minutes. Such a difference is significantly small, when examined in the context of the biological reproducibility for each experimental curve. It is interesting to note the constitutive protein level of the X-linked inhibitor of apoptosis protein (xIAP) (Figure 3d). This protein has been shown previously to directly inhibit caspase-3 [10]. The regulation between xIAP and caspase-3 is examined in detail in the Reaction rates space section. The strong similarities between the experimental and simulated data indicate that this model describes this Fas-mediated mechanism.

Table 1 Correlation between experimental and simulated data The p-values are calculated for a one-sided alternative hypothesis, assuming a positive association between the experimental and simulated curves for each protein (H0: β= 0 vs. H1: β > 0).
Figure 3
figure 3

Simulation results under normal cellularconditions. Simulated protein profiles for the seven proteins assayed. (a) Fas receptor and FADD protein. (b) Inactive (procaspase-8) and active (capase-8) caspase-8 proteins. Note that the two proteins are plotted on the same graph with different scales. (c) Inactive (procaspase-3) and active (capase-3) caspase-3 proteins. Note that the two proteins are plotted on the same graph with different scales. (d) xIAP protein.

The simulated procaspase-8 protein is least correlated with experimental results due to the abnormally large peak in level at ~ 4.5 hours, followed by a steep shift in slope direction after this point. Procaspase-8 is expressed at large protein levels, as compared to procaspase-3. Such a result is likely due to transcriptional activation of the procaspase-8 gene. Since this particular caspase is also involved in the internally-mediated (mitochondrial) mechanism of apoptosis, larger concentrations might be necessary for activation of competing protein cascades, not captured within this model representation. The simulated procaspase-8 curve resembles the experimental results up to the apex of the curve. After this point, the simulated curve continues to increase with positive slope, whereas the experimental curve decreased rapidly. The modeled representation of procaspase-8 was calculated to assume a protein expression pattern that both proceeds and resembles caspase-8, similar to the relationship exhibited between procaspase-3 and caspase-3. So this unexpected discrepancy between simulated and experimental results for procaspase-8 is likely the result of a transcriptional effect, not captured within the model representation.

Model of the Fas pathway under the influence of SEB exposure

In order to develop a model of this pathway under SEB-influenced conditions, the previously described experiments were performed with the addition of recombinant SEB. Compared to normal conditions, the experimental data shows significant differences (using a threshold of p < 0.05) in protein level patterns for three of the seven critical proteins – procaspase-8, caspase-3, and xIAP. In addition, flow cytometry results for the SEB-exposed cells indicate a positive shift in density for annexin-stained cells, as compared to cells stimulated with only the anti-Fas antibody at ~ 8 hours (data not shown). This distribution difference in the former illustrates how exposure to SEB (and anti-Fas antibody) induces apoptosis at a greater rate than cells stimulated with the anti-Fas antibody. Table 2 contains the T2 statistics, F-values, degrees of freedom (df), and p-values for the protein comparisons.

Table 2 Summary of Hotelling's T2 two-sample test. For a size α = 0.05 test, there is a significant (denoted with a '*' symbol) difference between normal and SEB-exposed conditions for procaspase-8, caspase-3, and xIAP.

When visualizing the curve patterns for these three proteins, the differences between conditions are apparent (Figure 4). The inactive procaspase-8 protein reaches maximum levels at 4.41 ± 0.10 hours under normal cellular conditions (Figure 4a), whereas in SEB-exposed cells, this protein reaches maximum levels at 2.20 ± 0.22 hours (Figure 4b). This is a time lead difference of 2 hours 13 minutes (± 14 minutes).

Figure 4
figure 4

Experimental data comparison between normal andSEB-exposed conditions. (a) Procaspase-8 protein under normal cellular conditions. (b) Procaspase-8 protein under SEB-exposed conditions. (c) Caspase-3 protein under normal cellular conditions. (d) Caspase-3 protein under SEB-exposed conditions. (e) xIAP protein under normal cellular conditions. (f) xIAP protein under SEB-exposed conditions. Points in each plot are from a random sample from experimental values. Red dashed lines are placed at the time point where the peak occurs for each curve while blue lines represent 95% point-wise and simultaneous Bonferonni confidence bounds.

The active caspase-3 protein shows a similar pattern to procaspase-8. Under normal cellular conditions, this protein reaches a maximum level at 4.65 ± 0.18 hours (Figure 4c) (approximately 14 minutes after procaspase-8), whereas under conditions of SEB exposure, caspase-3 reaches a maximum level at 2.94 ± 0.16 hours (Figure 4d). This is a time lead difference of 1 hour 43 minutes (± 14 minutes). The result that under SEB exposure, both caspase-3 and procaspase-8 show accelerated rates by approximately 2 hours implies that these two proteins are coupled to some mechanism induced by SEB. This association, however, is likely related to another pathway and not to this Fas-mediated model, since the intermediates in the cascade (procaspase-3 and active caspase-8) do not exhibit a similar accelerated pattern.

The xIAP protein seems to demonstrate the largest apparent difference of the three. Under normal cellular conditions, this protein is at a constitutive level (considering confidence bounds-Figure 4e), whereas under SEB-exposed conditions, this protein is down-regulated, or essentially degrading at a fast rate (Figure 4f). The maximum level for xIAP under SEB-exposed conditions is at t = 0. Previous work has demonstrated xIAP to be a direct inhibitor of caspase-3. Our experimental results support this finding in that under normal conditions, xIAP is at a constitutive level. However, under conditions of SEB exposure, xIAP is inhibited, or simply negatively affected, thereby reducing the repressive effect on caspase-3 and consequently increasing the rate of caspase-3. This result implies two notions: (1) SEB influences the degradation of xIAP, thereby both releasing a regulatory effect and increasing the rate of caspase-3, and subsequently, apoptosis, and (2) under SEB-exposed conditions, increasing the rate of xIAP can be a potential solution to restoring the level of caspase-3 and the normal rate of apoptosis. Such an implication is evaluated next for possible solutions within the rate constant space.

Reaction rate space

To understand the conditions under which caspase-3 is both expressed normally and at an increased rate (under the influence of SEB exposure), the relative kinetic rate space is sampled to define regions of normal and perturbed cellular conditions. As was explained previously, caspase-3 is one of the final proteins leading to DNA damage and cell death in the Fas-mediated mechanism, so understanding the parameters under which this protein is synthesized provides control over the outcome of apoptosis. Such boundaries are necessary to identify those areas within the plane associated with the system behaviour. This perspective on the allowable rate space provides a regulatory mechanism to revert a perturbed cellular system back to a stable normal state.

Two possible scenarios were examined within this analysis. First, the two rate constants from the primary regulatory reaction between xIAP and caspase-3 were simulated for solutions under normal cellular conditions (Figure 1-ultimate reaction). This criterion encompasses normal caspase-3 synthesis, where the expression peak occurs at 4.65 hours and xIAP is either constitutively, or simply positively expressed (Figures 4c,e). The second scenario simulates the rate constants where the peak in expression of caspase-3 is increased from 4.65 hours and the expression of xIAP is a decreasing function (Figures 4d,f), resembling either degradation or decreasing expression. Possible example solutions to these two scenarios are provided in Figure 5. For each plot, the bold line represents the same average protein profile as those observed in Figure 4. The dashed lines represent examples of potential solutions that are simulated. Note that in the dashed line solutions, the protein levels that are simulated can vary, however, the peak rates and linear trends cannot.

Figure 5
figure 5

Caspase-3 and xIAP under normal and simulated SEBconditions. (a) Caspase-3 potential solutions where a peak in expression occurs at 4.65 ± 0.18 hours. (b) xIAP potential solutions where the protein is either constitutively expressed or simply positively expressed after t = 0. (c) Simulated potential solutions for caspase-3, where the rate of synthesis is increased from 4.65 hours. Note that the level of expression can fluctuate with an increased rate. (d) Simulated potential solutions for xIAP, where the protein has negative expression levels after t = 0. The multiple dashed lines for each plot represent the various degrees of potential protein levels that can occur, while maintaining the expected pattern.

Within the model, the primary regulatory reaction between xIAP and caspase-3 is dictated by two relative rate constants, allowing the plane to be readily visualized within two dimensions (Figure 6). This extraction of the solution space for the two rates under both conditions is represented with green circles for normal states and blue crosses for SEB-simulated states, where each point represents a rate value for the forward (k15) and reverse (k16) reaction. It is interesting to note the peak pattern for the possible rate values necessary to achieve SEB-exposed conditions (blue crosses). Such points are indicative of caspase-3 expression at an increased rate and a decreasing expression function for xIAP. As the k16 rate values increase, the range of k15 values begins to decrease, for allowable solutions to the SEB-exposed condition. Though the allowable rates for the model were restricted to the range of 1.0 × 10-9 to 18 (justified in the Methods section), the y-axis (representative of k16) was simulated for values up to 70, to identify the highest point for the SEB-simulated solution space. The rate space under normal cellular conditions is restricted to two small areas within the plane: very small values of k15 and k16, and large values of k16 for k15 limited to the range of 5.5 to 9. Simulations using the rate values within these two small regions provide results associated with normal cellular conditions for caspase-3 and xIAP, as determined by agreement with experimental data. In attempt to revert a perturbed system from SEB-exposed to normal cellular conditions, this representation can guide the necessary directionality and magnitude of rate adjustments. For example, at relative rate values of k15 = 8 and k16 = 30, the system exhibits behaviour indicative of SEB exposure, however, an increase in the reverse rate constant (k16) to values greater than 40 can restore the expression of both caspase-3 and xIAP to behaviours simulating a normal system. This result is counterintuitive when first examining this reaction, mainly because an increase in the reverse reaction rate, k16 would imply an increase in the rate of caspase-3 synthesis, rather than a reduction. However, it is important to remember that k16 is the sole synthesis parameter for xIAP and the level of xIAP is but a single component of one degradation term for caspase-3. So k16 values that occur greater than a certain threshold (values greater than 40 in this example) dramatically increase xIAP to high levels, which simultaneously increases the degradation term of caspase-3, thus slowing down the rate. To simulate this result within a cellular system, this is analogous to increasing the cellular concentration of xIAP to a point where apoptosis can be slowed down, indicative of normal cellular conditions.

Figure 6
figure 6

Rate space for caspase-3 and xIAP synthesis forsingle regulatory reaction. Plot of simulated solutions to normal conditions and SEB-exposed conditions within varied values of relative rate constants k16 vs. k15. Normally distributed random noise (jitter) has been added to the green points to reduce the overlap and illustrate the frequency of simulations.


We have constructed a parsimonious model representation of Fas-mediated apoptosis under both normal and SEB-influenced cellular conditions. Our simulated results agree with experimental data as evaluated by correlation. In addition, experimental results demonstrate an increased rate of apoptosis and significant differences for three proteins (procaspase-8, caspase-3, and xIAP) between normal and SEB-exposed conditions.

In addition to the direct effect of SEB on the Fas pathway, there are several other factors that could lead to this acceleration of apoptosis. First, several factors could be causing an increase in FasL concentration. Specifically, both SEB treatment of various cell types (including Jurkats) and activation of the Fas receptor [13, 14] have been shown to cause the de novo production and secretion of soluble FasL (sFasL). Also, SEB treatment has been shown to induce the production of several members of the (tumor necrosis factor) TNF superfamily such as TNFα and tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) [11]. Increased TNFα levels have been shown to lead to proteolytic cleavage of existing FasL from the cell surface [12]. This increase in sFasL would create a positive feedback loop thus increasing the rate of apoptosis.

Second, the increase in TRAIL expression has been shown to initiate Jurkat cell apoptosis through the DR5 receptor [15]. In addition to promoting apoptosis, activation of this parallel pathway eventually leads to an increase in caspase-8 expression and may partially explain our results (discussed further below) [15]. The addition of other pathways into our minimal model will be the subject of further work.

Finally, the induction of fratricide apoptosis in which activated T cells begin inducing apoptosis in neighboring cell as a result of SEB exposure may play a role in this rate increase as well [16].

It is interesting to note that the increased rate for procaspase-8 is not complimented by an increase in proteolytic cleavage and subsequent activation of caspase-8. This implies that SEB is inducing the increased rate on the gene expression point, since the inactive caspase demonstrates a significant difference in the presence of SEB exposure, as opposed to the active caspase.

This accelerated rate is likely the result of SEB affecting several signal transduction pathways that regulate caspase-8 expression. Studies have shown that SEB stimulation of Jurkat cells causes them to begin secreting various cytokines such as interferon gamma (IFNγ) within 30 minutes following exposure [1719] The release of IFN-γ can then increase the rate of transcription of caspase-8 in an interferon regulatory factor 1 (IRF1) dependent manner. [18, 20, 21].

The accelerated rate of caspase-3 and degradation of xIAP under SEB influence imply a repressive regulation of xIAP under normal conditions, such that SEB exposure releases the repressive effect and increases the rate of caspase-3 and subsequently, apoptosis. SEB has been shown to induce the activation of caspase-3 [12], however, the combined effect on the interaction (indirect or direct) between SEB and xIAP/caspase-3 is not well understood. Under the assumption of a mitochondrial mechanism (internally mediated apoptosis), SEB could inhibit xIAP similar to how synthetic Smac/DIABLO has been demonstrated to, thereby leading to an increased activity of caspase-3 and caspase-9 in etoposide-induced cells [22]. However, since autocrine release of fas was observed, and previous work has shown the SEB treatment induces FasL release [12], the mechanism is more likely not primarily mediated through the mitochondrial pathway, and instead, fas-mediated.

Simulation results also demonstrate a means to instil control on an SEB-exposed system. Utilizing the primary inhibitory reaction between caspase-3 and xIAP, the relative rate space was computed to illustrate values for converting a perturbed system back to normal cellular conditions. Such results indicate that at specific increased concentrations of xIAP, the rate of apoptosis can be reduced from an SEB-exposed rate to a normal cellular rate.


In future work, we plan to expand the model to account for the protein intermediates that are observed to have an effect on procaspase-8, in addition to the indirect connection between increased rates of procapase-8 and caspase-3 under the influence of SEB.


Flow cytometry

After the cells became confluent in free RPMI1640 with 1.0 × 106/ml cells, they were split equally into three flasks. One flask was incubated with anti-Fas (0.01 ug/ml) and the other was incubated with the combination of anti-Fas (0.01 ug/ml) and SEB (1 ug/ml). The control group was untreated. Cells were then aspirated out from three flasks after 2, 4, and 8 hour of incubations. They were washed with cold (2–8°C) 1X PBS, and centrifuged at 500 × g for 5 min at room temperature before labeled with Annexin V-FITC, PI, or a combination of Annexin V-FITC and PI (amounts of reagents added are according to the manufacture protocol, TACS Annexin V-FITC Apoptosis Detection Kit, R&D systems), respectively. Cells were analyzed in a Becton Dickinson FACSCaliber flow cytometer using Cell Quest software. At least 10,000 viable cells per condition were analyzed.

Western blot analysis

Cell lines and cell culture media were purchased from ATCC (Manassas, VA). Cells were maintained in RPMI medium supplemented with 10% FBS and antibiotics. The Fas pathway was activated by the addition of 0.01 μg/ml of an anti-Fas activating antibody (clone DX2, BD Pharmingen, San Diego, CA) for normal cellular conditions and 1 μg/ml of SEB (Sigma Chemical, St. Louis, MO), in addition to the anti- Fas antibody, for perturbed conditions within the time range indicated (0–12 hours).

All blots were performed using standard procedures. In summary: following treatment, cells were lysed using M-PER lysis buffer (Pierce, Rockford IL) and the protein concentration was determined by BCA assay. A mass of 50 μg of total protein was loaded on a 10% polyacrylamide gel and separated by electrophoresis. Proteins were transferred to a nitrocellulose membrane and probed with specific antibodies (below) for 1 hour at room temperature. Blots were developed using ImmunoStar ECL reagent (BioRad, Hercules, CA) and imaged using a VersaDoc Imaging System (BioRad, Hercules, CA). Relative protein concentrations were determined using QuantityOne software (BioRad, Hercules, CA). Initial concentration values for the seven proteins assayed were simultaneously determined on the same nitrocellulose membrane.

The antibodies (Table 3) were obtained from Santa Cruz Biotechnology (Santa Cruz, CA) or R&D Systems (Minneapolis, MN) and used at 1 μg/ml.

Table 3 Antibody information

Relative protein level curves

Mean protein levels were calculated from a minimum of three independent replicates for each protein, ranging to nine replicates for some proteins. The five points for each protein curve (five time states: 0, 2, 4, 8, and 12 hours) were increased to 50 points with linear interpolation. Curve kurtosis was then reduced with a local polynomial regression fitting (loess smoothing). The 95% point-wise and simultaneous Bonferroni confidence intervals were calculated from the variance for each curve. All error terms and propagation of error given for time point intervals are calculated with the standard error of the mean.

Ordinary differential equations (ODEs)

Initial value ODEs were determined from the model reactions (see Appendix) for each protein/protein complex. The equations were solved using a classical Runge-Kutta 4th order integration function available in the 'odesolve' package of R [23]. The protein with the minimum initial level (active caspase-3) was used as the basal value for the protein complexes at (t = 0). This was determined to support the assumption that proteolytic activation of caspases (caspase-3 and caspase8) and protein complex formation occurs at minimum levels (thus at low concentrations) at initial time points.


The optimization problem seeks to minimize the aggregate Euclidean distance between the experimental and simulated data given as:

( 1.1 ) d a = j = 1 n ( i = 1 t ( x j i y j i ) 2 ) 1 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqabeqacaaabaWaaeWaaeaacqaIXaqmcqGGUaGlcqaIXaqmaiaawIcacaGLPaaaaeaacqWGKbazdaWgaaWcbaGaemyyaegabeaakiabg2da9maaqahabaGaeiikaGYaaabCaeaacqGGOaakcqWG4baEdaWgaaWcbaGaemOAaOMaemyAaKgabeaakiabgkHiTiabdMha5naaBaaaleaacqWGQbGAcqWGPbqAaeqaaOGaeiykaKYaaWbaaSqabeaacqaIYaGmaaGccqGGPaqkdaahaaWcbeqaamaalmaameaacqaIXaqmaeaacqaIYaGmaaaaaaWcbaGaemyAaKMaeyypa0JaeGymaedabaGaemiDaqhaniabggHiLdaaleaacqWGQbGAcqGH9aqpcqaIXaqmaeaacqWGUbGBa0GaeyyeIuoaaaaaaa@5354@

where x and y are simulated and experimental protein concentrations, respectively, n is the total number of proteins, and t is the total number of times states. The 16 relative kinetic rate constants (parameters) were computed using a two-phase minimization approach over multiple trials. Due to the dependence between initial value seeding and the minimized function (local optima), a Monte Carlo simulation [24] was used to iteratively seed the initial parameter values to begin each optimization. The distance function was first minimized with a Newton-type non-linear procedure available in the 'stats' package of R [23]. The parameter values determined at convergence of the distance function were then seeded into a subsequent quasi-Newton minimization method with box constraints [25]. This method constrains the parameter values to a specified range in the minimization of the distance function. This final aggregate distance was stored in a list for each simulated trial.

For each iteration in the Monte Carlo simulation, the parameter values were sampled from a random uniform distribution with a lower bound near zero and an upper bound of λ. The λ variable was increased from an initial value of nine to a maximum value of λ = 27. This procedure was implemented for determination of the minimum allowable parameter value range with the minimum distance function. Such a method constrains the rate constants to a small range of values. The results are represented in Figure 2. Each line corresponds to a different maximum parameter value λ over a series of 150 minimization trials. When the 16 parameter values are confined to the range of near-zero (1.0 × 10-9) to λ = 9, the minimum aggregate distance obtained over 150 trials is 45.48. However, as the maximum allowable parameter value is increased to λ = 18, the minimum aggregate distance is reduced to 44.46. Additionally, all maximum parameter values λ ≥ 18 converge to the same minimum aggregate distance. This result indicates that the distance function is not reduced any further with a maximum parameter value greater than λ = 18. As such, the range of rate constants was limited to values of λ = 1.0 × 10-9 to λ = 18.

Hotelling's two-sample T2 test

For assignment of a statistical summary to the difference between each protein under normal and SEB-exposed conditions, a T2 test was calculated, using all replicates for each protein. This multivariate analog to the Student's t-test is implemented to account for the correlation structure between the t time state variables. The squared distance is given by:

(2.1)     T2 = n ( X ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWGybawgaqeaaaa@2DFD@ - μ0)T (S)-1 ( X ¯ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWGybawgaqeaaaa@2DFD@ - μ0), where

( 2.2 ) X ¯ ( t × 1 ) = 1 n j = 1 n X i MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqabeqacaaabaWaaeWaaeaacqaIYaGmcqGGUaGlcqaIYaGmaiaawIcacaGLPaaaaeaacuWGybawgaqeamaaBaaaleaacqGGOaakcqWG0baDcqGHxdaTcqaIXaqmcqGGPaqkaeqaaaaakiabg2da9maalaaabaGaeGymaedabaGaemOBa4gaamaaqahabaGaemiwaG1aaSbaaSqaaiabdMgaPbqabaaabaGaemOAaOMaeyypa0JaeGymaedabaGaemOBa4ganiabggHiLdaaaa@45D5@

, and

( 2.3 ) S ( t × t ) = 1 n 1 j = 1 n ( X i X ¯ ) ( X i X ¯ ) T MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqabeqacaaabaWaaeWaaeaacqaIYaGmcqGGUaGlcqaIZaWmaiaawIcacaGLPaaaaeaacqWGtbWudaWgaaWcbaGaeiikaGIaemiDaqNaey41aqRaemiDaqNaeiykaKcabeaaaaGccqGH9aqpdaWcaaqaaiabigdaXaqaaiabd6gaUjabgkHiTiabigdaXaaadaaeWbqaaiabcIcaOiabdIfaynaaBaaaleaacqWGPbqAaeqaaOGaeyOeI0IafmiwaGLbaebacqGGPaqkcqGGOaakcqWGybawdaWgaaWcbaGaemyAaKgabeaakiabgkHiTiqbdIfayzaaraGaeiykaKYaaWbaaSqabeaacqqGubavaaaabaGaemOAaOMaeyypa0JaeGymaedabaGaemOBa4ganiabggHiLdaaaa@5423@


( 2.4 ) μ 0 ( t × 1 ) = [ μ 10 μ 20 . . μ t 0 ] MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqaaeqacaaabaWaaeWaaeaacqaIYaGmcqGGUaGlcqaI0aanaiaawIcacaGLPaaaaeaaiiGacqWF8oqBdaWgaaWcbaGaeGimaaZaaSbaaWqaaiabcIcaOiabdsha0jabgEna0kabigdaXiabcMcaPaqabaaaleqaaOGaeyypa0ZaamWaaeaafaqaaeqbbaaaaeaacqWF8oqBdaWgaaWcbaGaeGymaeJaeGimaadabeaaaOqaaiab=X7aTnaaBaaaleaacqaIYaGmcqaIWaamaeqaaaGcbaGaeiOla4cabaGaeiOla4cabaGae8hVd02aaSbaaSqaaiabdsha0jabicdaWaqabaaaaaGccaGLBbGaayzxaaaaaaaa@4AF4@

such that at large values of T2, the hypothesis H0: μ = μ0 is rejected. The T2 statistic is distributed as ( n 1 ) t ( n t ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaamaabmaabaGaemOBa4MaeyOeI0IaeGymaedacaGLOaGaayzkaaGaemiDaqhabaWaaeWaaeaacqWGUbGBcqGHsislcqWG0baDaiaawIcacaGLPaaaaaaaaa@3844@ Ft,n-t, where Ft,n-t, is a random variable that follows an F – distribution with t and n - t degrees of freedom [26]. For each comparison, Bartlett's test was computed to verify the homogeneity of the covariance matrices


  1. Agrawal A, Pulendran B: Anthrax lethal toxin: a weapon of multisystem destruction. Cell Mol Life Sci. 2004, 61 (22): 2859-65. 10.1007/s00018-004-4251-4.

    Article  CAS  PubMed  Google Scholar 

  2. Bonfoco E, Stuart PM, Brunner T, Lin T, Griffith TS, Gao Y, Nakajima H, Henkart PA, Ferguson TA, Green DR: Inducible nonlymphoid expression of Fas ligand is responsible for superantigen-induced peripheral deletion of T cells. Immunity. 1998, 9 (5): 711-720. 10.1016/S1074-7613(00)80668-8.

    Article  CAS  PubMed  Google Scholar 

  3. Popov SG, Villasmil R, Bernardi J, Grene E, Cardwell J, Wu A, Alibek D, Bailey C, Alibek K: Lethal toxin of Bacillus anthracis causes apoptosis of macrophages. Biochem Biophys Res Commun. 2002, 293 (1): 349-55. 10.1016/S0006-291X(02)00227-9.

    Article  CAS  PubMed  Google Scholar 

  4. Sohn MH, Kim JW, Kim WK, Jang GC, Kim KE: Staphylococcal enterotoxin B upregulates fas-mediated apoptosis of peripheral blood mononuclear cells in childhood atopic dermatitis. Scand J Immunol. 2003, 57 (1): 62-67. 10.1046/j.1365-3083.2003.01183.x.

    Article  CAS  PubMed  Google Scholar 

  5. Jacobson MD, McCarthy N: Apoptosis the molecular biology of programmed cell death. 2002, New York: Oxford University Press

    Google Scholar 

  6. Lin Y, Devin A, Rodriguez Y, Liu ZG: Cleavage of the death domain kinase RIP by Caspase-8 prompts TNF-induced apoptosis. Genes Dev. 1999, 13 (19): 2514-2526. 10.1101/gad.13.19.2514.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Scott FL, Denault J, Riedl R, Shin H, Renatus M, Salvesen GS: XIAP inhibits caspase-3 and -7 using two binding sites: evolutionarily conserved mechanism of IAPs. EMBO. 2005, 24 (3): 645-655. 10.1038/sj.emboj.7600544.

    Article  CAS  Google Scholar 

  8. Silke J, Ekert PG, Day CL, Hawkins CJ, Baca M, Chew J, Pakusch M, Verhagen AM, Vaux DL: Direct inhibition of caspase3 is dispensable for the anti-apoptotic activity of XIAP. EMBO. 2001, 20 (12): 3114-3123. 10.1093/emboj/20.12.3114.

    Article  CAS  Google Scholar 

  9. Fussenegger M, Bailey JE, Varner J: A mathematical model of caspase function in apoptosis. Nature Biotechnology. 2000, 18: 768-774. 10.1038/81208.

    Article  CAS  PubMed  Google Scholar 

  10. Verhagen A, Coulson E, Vaux D: Inhibitor of apoptosis proteins and their relatives: IAPs and other BIRPs. Genome Biol. 2001, 2 (7): 1-reviews3009. 10.1186/gb-2001-2-7-reviews3009.

    Article  Google Scholar 

  11. Droin NM, Pinkoski MJ, Dejardin E, Green DR: Egr family members regulate nonlymphoid expression of Fas ligand, TRAIL, and tumor necrosis factor during immune responses. Mol Cell Biol. 2003, 7638-47. 10.1128/MCB.23.21.7638-7647.2003. 21

  12. Dhein J, Walczak H, Baumler C, Debatin KM, Krammer PH: Autocrine T-cell suicide mediated by APO-1/(Fas/CD95). Nature. 1995, 373 (6513): 385-6. 10.1038/373438a0.

    Article  Google Scholar 

  13. Ettinger R, Panka DJ, Wang JK, Stanger BZ, Ju ST, Marshak-Rothstein A: Fas ligand-mediated cytotoxicity is directly responsible for apoptosis of normal CD4+ T cells responding to a bacterial superantigen. J Immunol. 1995, 154: 4302-4308.

    CAS  PubMed  Google Scholar 

  14. Sohn MH, Kim JW, Kim WK, Jang GC, Kim KE: Staphylococcalenterotoxin B upregulates fas-mediated apoptosis of peripheral blood mononuclear cells in childhood atopic dermatitis. Scand J Immunol. 2003, 57 (1): 62-7. 10.1046/j.1365-3083.2003.01183.x.

    Article  CAS  PubMed  Google Scholar 

  15. Jang YJ, Park KS, Chung HY, Kim HI: Analysis of the phenotypes of Jurkat clones with different TRAIL-sensitivities. Cancer Lett. 2003, 194 (1): 107-17. 10.1016/S0304-3835(02)00680-8.

    Article  CAS  PubMed  Google Scholar 

  16. Gorak-Stolinska P, Kemeny DM, Noble A: Activation-induced cell death in human T cells is a suicidal process regulated by cell density but superantigen induces T cell fratricide. Cell Immunol. 2002, 19 (2): 98-107. 10.1016/S0008-8749(02)00598-1.

    Article  Google Scholar 

  17. Maeurer M, Zitvogel L, Elder E, Storkus WJ, Lotze MT: Human intestinal V delta 1+ T cells obtained from patients with colon cancer respond exclusively to SEB but not to SEA. Nat Immun. 1995, 14 (4): 188-197.

    CAS  PubMed  Google Scholar 

  18. Shinohara F, Sato K, Suzuki M, Takoda H, Echigo S, Rikiishi H: Role of caspase in CD80 expression of superantigen-stimulated monocytes. Int J Mol Med. 2004, 14 (2): 241-6.

    CAS  PubMed  Google Scholar 

  19. Miethke T, Wahl C, Heeg K, Echtenacher B, Krammer PH, Wagner H: T cell-mediated lethal shock triggered in mice by the superantigen staphylococcal enterotoxin B: critical role of tumor necrosis factor. J Exp Med. 1992, 175 (1): 91-8. 10.1084/jem.175.1.91.

    Article  CAS  PubMed  Google Scholar 

  20. Refaeli Y, Van Parijs L, Alexander SI, Abbas AK: Interferon gamma is required for activation-induced death of T lymphocytes. J Exp Med. 2002, 196 (7): 999-1005. 10.1084/jem.20020666.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Ruiz-Ruz C, Ruiz de Almodovar C, Rodriguez A, Ortiz-Ferron G, Redondo JM, Lopez-Revas A: The up-regulation of human caspase 8 by interferon-gamma in breast tumor cells requires the induction and action of the transcription factor interferon regulatory factor-1. J Biol Chem. 2004, 279 (19): 19712-19720. 10.1074/jbc.M313023200.

    Article  Google Scholar 

  22. Mizukawa K, Kawamura A, Sasayama T, Tanaka K, Kamei M, Sasaki M, Kohmura E: Synthetic Smac peptide enhances the effect ofetoposide-induced apoptosis in human glioblastoma celllines. J Neurooncol. 2006, Mar 31,

  23. R Development Core Team: R: A Language and Environment for Statistical Computing. 2004, R Foundation for Statistical Computing, Vienna, Austria

    Google Scholar 

  24. Metropolis , Nicholas , Ulam : The Monte Carlo method. Journal of the American Statistical Association. 1949, 44 (247): 335-341. 10.2307/2280232.

    Article  CAS  PubMed  Google Scholar 

  25. Byrd RH, Lu P, Nocedal J, Zhu C: A limited memory algorithm for bound constrained optimization. SIAM J Scientific Computing. 1995, 16: 1190-1208. 10.1137/0916069.

    Article  Google Scholar 

  26. Johnson R, Wichern D: Applied Multivariate Statistical Analysis. 2002, New Jersey: Prentice Hall

    Google Scholar 

Download references


The authors would like to thank the anonymous reviewers of this manuscript for their valuable criticism, as well as the members of the BioSpice community for inspiring the research focus of this work.

Author information

Authors and Affiliations


Corresponding author

Correspondence to John Dileo.

Additional information

Authors' contributions

The modeling and simulation was conducted primarily by BWH and WEC. The experimental results were obtained by HBS and JD. BWH, JD, WEC, and HBS drafted the manuscript. The critical assessment of results and scientific discussions were contributed by OJP, JCF, RH, and MJ. All authors read and approved the final manuscript.

Brandon W Higgs, John Dileo, Wenling E Chang, Olivia J Peters and Jordan C Feidler contributed equally to this work.

Electronic supplementary material


Additional File 1: Model equations. Series of 13 initial value ODEs used for the model of Fas-mediated apoptosis. (PDF 12 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Higgs, B.W., Dileo, J., Chang, W.E. et al. Modeling the effects of a Staphylococcal Enterotoxin B (SEB) on the apoptosis pathway. BMC Microbiol 6, 48 (2006).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: