Skip to main content
  • Research article
  • Open access
  • Published:

Modelling chemotherapy effects on granulopoiesis



Although the growth-factor G-CSF is widely used to prevent granulotoxic side effects of cytotoxic chemotherapies, its optimal use is still unknown since treatment outcome depends on many parameters such as dosing and timing of chemotherapies, pharmaceutical derivative of G-CSF used and individual risk factors. We showed in the past that a pharmacokinetic and –dynamic model of G-CSF and human granulopoiesis can be used to predict the performance of yet untested G-CSF schedules. However, only a single chemotherapy was considered so far.

In the present paper, we propose a comprehensive model of chemotherapy toxicity and combine it with our cell kinetic model of granulopoiesis. Major assumptions are: proportionality of cell numbers and cell loss, delayed action of chemotherapy, drug, drug-dose and cell stage specific toxicities, no interaction of drugs and higher toxicity of drugs at the first time of application. Correspondingly, chemotherapies can be characterized by a set of toxicity parameters which can be estimated by fitting the predictions of our model to clinical time series data of patients under therapy. Data were either extracted from the literature or were received from cooperating clinical study groups.


Model assumptions proved to be feasible in explaining granulotoxicity of 10 different chemotherapeutic drugs or drug-combinations applied in 33 different schedules with and without G-CSF. Risk groups of granulotoxicity were traced back to differences in toxicity parameters.


We established a comprehensive model of combined G-CSF and chemotherapy action in humans which allows us to predict and compare the outcome of alternative G-CSF schedules. We aim to apply the model in different clinical contexts to optimize and individualize G-CSF treatment.


The effectivity of antineoplastic chemotherapy of some cancer types, such as lymphomas or breast cancer, depends on dose density of applied cytostatic drugs [1]-[6]. Dose density is defined as the amount of drug given per body surface per time unit (mg/m2/week) [7],[8]. It has been shown that a decrease in dose density such as treatment delays or dose reductions, can have negative impact on remission rates, recurrence rates and overall survival rates [9]-[15].

Physicians are frequently forced to reduce dose density due to serious chemotherapy-associated side effects, of which neutropenia, i.e. a reduction of white blood cells, is the most common one. Because neutrophils are an essential part of the nonspecific immune system, neutropenic patients are prone to bacterial and fungal infections, frequently resulting in an increased need of antibiotics, prolonged hospitalization and a higher risk of therapy discontinuation [16]-[22].

To ameliorate neutropenia, the recombinant haematopoietic growth factor G-CSF (granulocyte colony stimulating factor) is routinely applied. It is a major requirement to make dose-dense therapies feasible. G-CSF increases the mitotic activity, accelerates the maturation of different immature granuloid precursor cells in the bone marrow and increases the release of mature granulocytes into blood [23],[24]. Nowadays, a variety of pharmaceutical derivatives of G-CSF are available differing in pharmacokinetic and –dynamic properties. In consequence, combined chemotherapy and G-CSF treatment result in complex dynamics of granulocytes due to the interaction of G-CSF pharmacokinetics, G-CSF induced granulocytosis via different mechanisms and chemotherapy induced cell destruction.

We recently established a biomathematical model of G-CSF applications in humans explaining numerous scenarios of G-CSF applications of the most frequently used derivatives Filgrastim and Pegfilgrastim into healthy volunteers and first simple chemotherapies [25]. However, the large variety of chemotherapies supported by G-CSF in current clinical practice is not yet covered.

In the present article we introduce a substantially refined model of chemotherapy action applicable for different diseases and risk groups. The major objective of our model is to allow predictions of the dynamics of granulocytes after combined chemotherapy and G-CSF applications especially regarding effects of alternative, yet untested G-CSF treatment schedules on neutropenia.


We recently developed a comprehensive model of pharmacokinetics and pharmacodynamics of Filgrastim and Pegfilgrastim [25]. We briefly sketch this model in the following:

PK/PD model structure and basic properties

Figure 1 depicts the structure of the PK/PD-model. Granulopoiesis is divided into five distinct cell compartments, representing haematopoietic stem cells (S), early progenitors (CG), proliferating precursors (PGB), maturing precursors (MGB) and mature granulocytes in circulation (GRA). Dynamics of each compartment are described by ordinary differential equations (ODE). The system is regulated by several feed-back loops, mostly mediated by the cytokine G-CSF. Endogenous G-CSF production is regulated by the cell demand of the system.

Figure 1
figure 1

Structure of the granulopoiesis model under chemotherapy and G-CSF treatment. Boxes represent major cell- or cytokine compartments of the model: S = haematopoietic stem cells, CG = granulopoietic progenitors, PGB = granulopoietic precursors, MGB = maturing granulopoietic precursors in bone marrow, GRA = mature granulocytes in blood, LY = lymphocytes. We modelled two G-CSF derivatives (fil = Filgrastim, peg = Pegfilgrastim). Arrows represent cell/cytokine fluxes and interactions. CX represents the strength of chemotherapy (see below).

All granulopoietic cells originate from the stem cell compartment. Cell division and differentiation rates in S are regulated in such a way that self-maintenance gains priority in case of low stem cell numbers. Cells committed to granulopoietic lineages enter the subsequent compartment CG, which represents the most immature cell stage committed to granulopoiesis. The next compartment, PGB, represents the mitotic granuloid precursors (myeloblasts, promyelocytes and myelocytes). Compartment MGB represents all non-mitotic precursors (metamyelocytes, banded cells and mature neutrophils) in the bone marrow. The final compartment, GRA, comprises all mature neutrophils in the peripheral blood. Reductions of lymphocytes LY due to chemotherapy are also modelled.

Without chemotherapy, changes of compartment sizes are determined by balance equations of cell influx, cell production and cell efflux or degradation:

d C t d t = C in t A t C t T t

where d C t d t represents the changing rate of compartment size, C in (t) represents the cell influx rate from the preceding compartment, A(t) the amplification of cell numbers and T(t) is the average time of a cell residing in the compartment (transition time).

Filgrastim and Pegfilgrastim are assumed to have different pharmacodynamic properties in the model, whereas Filgrastim and endogenous G-CSF are assumed to be undistinguishable. Amplification rate and transition time in PGB increase with G-CSF serum concentration. In contrast, the transition time in MGB and the apoptosis rate are reduced with increasing G-CSF concentrations. Details of Filgrastim and Pegfilgrastim pharmacokinetics, pharmacodynamics, and corresponding regulatory processes are described elsewhere [25].

On the basis of this baseline model, we aim at developing a more comprehensive model of cytotoxic chemotherapy action on granulopoiesis. Corresponding assumptions and equations are explained in the following. Some detailed information can be found in [26].

Modeling chemotherapy

Next we present and discuss our assumptions reagarding chemotherapy modelling in detail.

Assumption 1 (injection)

Injection of chemotherapy is modeled by pulse functions according to the applied dosing and timing schedules. In general, each drug is modeled by a separate pulse function, where the length of the pulse corresponds to the injection time and the amplitude is normalized in such a way that the area under the curve after a single injection equals one (see equations below).

Assumption 2 (delayed action)

We assumed a delayed maximum of chemotherapy damage after injection. This is modeled in a phenomenological rather than mechanistic way by a set of concatenated first order transitions resulting in a delayed maximum after injection (see [27]).

d Ψ drug i t d t = Ψ drug _ out i 1 t k Delay drug Ψ t o x i t ,i=1,,4,with
Ψ drug 0 t = CHEMO drug t ,
CHEMO drug t = i = 1 N cycle ( Hv ( t t i )Hv t t i t inf )/ t inf
Ψ drug _ out i t = k Delay drug Ψ drug i t ,i=1,,4,

where Hv is the Heavyside function: Hv t = 1 : t > 0 0 : t 0 , t i are the time points of chemotherapy applications and tinf is the infusion time. Thus, function CHEMO represents the chemotherapy schedule. In summary, Ψ drug _ out 4 t represents the strength of the (delayed) toxic effect.

Assumption 3 (toxicity)

Drug, drug-dose and cell-stage specific toxicity functions are derived from Ψ drug _ out 4 t by multiplications with specific toxicity values: K X drug t

Ψ drug X t = K X drug t Ψ drug _ out 4 t

The quantities K X drug are called toxicity coefficients in the following.

Assumption 4 (first cycle effect)

The term “first cycle effect” refers to increased toxicity of chemotherapeutic drugs when applied for the first time. Accordingly, we assumed a ‘first cycle effect’ by multiplying the toxicity of the first chemotherapy cycle by a factor f c  ≥ 1. Hence,

K X drug t = f fc drug K X drug K X drug if t < t 2 drug else

Assumption 5 (no interactions between drugs)

Most chemotherapy regimens consist of multiple drugs administered simultaneously. If toxicity functions of single drugs are available, the resulting total toxicity is obtained by adding these functions. Thus, in general, no interactions between single drugs are assumed:

Ψ total X t = drug Ψ drug X

where “drug” summarizes all drugs applied in combination. The overall process of defining toxicity functions is sketched in Figure 2.

Figure 2
figure 2

Construction of toxicity functions exemplified by the high-CHOEP regimen and its toxic effect on stem cells. The Heavyside functions for a single application of the combination of Cyclophosphamide, Doxorubicin and Vincristin (at time point 0) and for single Etoposide (applications at time points 0, 1 and 2) are determined (first column), delayed (second column) and multiplied with the corresponding stem cell toxicity factors (third column). Finally, the functions were added resulting in the overall toxicity function.

Assumption 6 (cell loss)

Cytotoxic drugs cause a depletion of bone marrow cells. The loss rate is proportional to the number of cells in each compartment (first order kinetics). The overall toxicity function Ψ total X t defined above serves as proportionality factor, i.e.:

d C X t d t = C X in t A t C X t T t Ψ total X t C X t

where C X ist the content of compartment X.

Assumption 7 (risk groups)

Risk groups of patients with differing toxic response can be described by different sets of corresponding toxicity parameters. Motivated by the observation that G-CSF response does not differ between elderly and younger patients [28], we assumed the same cell kinetic parameters among risk groups.

Assumption 8 (lymphopoiesis)

Depletion of lymphocytes is modelled by the following simple equation:

C W B C t = C A N C t + C L Y e Ψ L Y t

C WBC , C ANC and C LY are the concentrations of white blood cells, neutrophils and lymphocytes in peripheral blood. The factor e Ψ L Y t quantifies the overall cytoreductive effect of chemotherapy on lymphocytes. Ψ LY (t) is analogously defined as toxicities regarding granulopoiesis. Note that lymphocyte dynamics are not explicitly modelled except for this chemotherapy effect. Thus, C LY  = 3000/μl is constant, i.e. the normal concentration of lymphocytes.

Assumption 9 (Prednisone)

Prednisone is a chemotherapeutic drug without cytotoxic effect. It is often applied for the treatment of malignant lymphoma (e.g., CHOP, CHOEP, highCHOEP and BEACOPP, see Table 1, “Summary of modelled chemotherapies”). It is well known that prednisone temporarily increases granulocyte counts caused by temporarily prolonged half-life [29]-[31].

Table 1 Summary of modelled chemotherapies: we present names of therapies and corresponding dosage of drugs

Assumption 10 (Cell kinetic parameters are un-affected by chemotherapy)

Besides the above mentioned (toxic) effects of chemotherapy, it is assumed that all cell-kinetic parameters of our granulopoiesis model remain unchanged. Specifically, we do not assume irreversible deterioration of granulopoiesis or reduced G-CSF response of cell stages during or after chemotherapy.

Clinical data

Since the toxic effects of chemotherapies on bone marrow cannot directly be observed, at least one clinical dataset of patients under therapy is required to estimate it. Different schedules of the chemotherapy and different schedules of supportive G-CSF treatment provide additional information which can be used to improve toxicity estimates or to validate model predictions. Data of G-CSF application into healthy volunteers were already used to validate our pharmacokinetic and –dynamic model of G-CSF [25].

Data of 10 different chemotherapies used to treat Hodgkin’s lymphoma (HD), non-Hodgkin’s disease (NHL), breast cancer (BC) and non-small cell lung cancer (NSCLC) are available (Table 1), either from literature or from cooperating clinical study groups (German Hodgkin’s Lymphoma Study Group (Professor Engert), German High-Grade Non-Hodgkin’s Lymphoma Study Group (Professor Pfreundschuh), German Breast Group (Professor von Minckwitz)). Considering different cytotoxic drug and G-CSF schedules, data of 33 different chemotherapies are available (Table 2).

Table 2 Clinical data sets used for modelling: we present disease, chemotherapy protocol and corresponding G-CSF schedules

Data sets comprise time series data of G-CSF serum concentrations, ANC or WBC of patients under therapy. For modelling issues, we used patient’s medians throughout.

Parameter estimation

Pharmakokinetic and pharmakodynamic parameters for Filgrastim and Pegfilgrastim are described elsewhere and remained unchanged in the present work [25]. In our model, the toxicity of a chemotherapy regimen is characterised by a set of cell-stage and drug specific toxicity parameters (for S, CG, PGB, MGB, LY) and a drug specific delay parameter. Parameter estimation was realised using an algorithm based on evolutionary strategies. Evolutionary strategies are stochastic algorithms used for numerical optimization [48]. The cost function to calculate model fitness was defined as:

t 0 t 1 | log f model t , k log f data t |dtmin,

where t 0 is the time of the first data point, t 1 is the last data point, f model (t, k) is the solution of the model equation system for the granulocyte compartment at the time of t (t 0 ≤ t ≤ t 1) based on the parameter set k = {k 1, …, k n } and f data (t) is the linearly interpolated data curve. Agreement of logarithms was pursued since cell counts are usually log-normally distributed.

We split the data of the NHL trial (CHOP-like chemotherapies of high-grade non-Hodgkin’s disease) into young and elderly patients to account for risk specific toxicities (chemotherapy assumption 7, see above). The toxicity parameters were estimated in a stepwise manner starting with simple chemotherapies which require only a few parameters estimates. More complex chemotherapies were modelled by estimating toxicity parameter sets for yet unconsidered drugs or drug combinations. Toxicity parameters estimated in earlier fitting steps were kept constant throughout the fitting process. The parameters for the first cycle effect and the two delay parameters are kept constant for different dose levels and for young and elderly patients as well. If drugs are always applied in combinations, it is impossible to separate the toxic effects of its components. In these cases, a single set of toxicity parameters was estimated for the combination.

In more detail:

  1. 1.

    We estimated CHOP parameters separately for elderly patients (scenarios 13, 15, 21, 30, 31 of Table 2, chemotherapy dosings can be found in Table 1), and young patients (scenarios 14, 20, 33). We assumed higher toxicity for elderly patients [49].

  2. 2.

    Using CHOEP data sets and the parameters found in step 1, we determined parameter settings for Etoposide 100 mg/m2 for young (scenarios 16, 23) and elderly patients (scenarios 17, 24).

  3. 3.

    Using BEACOPP basis data sets 18 and 26 and the parameters for Etoposide 100 for young patients estimated in step 2 we determined parameter settings for the combination of Cyclophosphamide 650 mg/m2, Doxorubicin 25 mg/m2 and Vincristine 2 mg, with the constraint that the parameter values of this combination must be smaller than those for CHOP young (because of lower or equal dosage of Cyclophosphamide, Doxorubicine and Vincristine compared to CHOP). Parameters for Bleomycin 10 mg/m2 and Procarbacine 100 mg/m2 were also determined.

  4. 4.

    With the data set 27 (BEACOPP escalated) the parameter settings for the combination of Cyclophosphamide 1250 mg/m2, Doxorubicin 35 mg/m2 and Vincristine 2 mg, and for Etoposide 200 mg/m2 were estimated with the constraint that the parameter values must be larger than those estimated for BEACOPP basis.

  5. 5.

    Taking the high-CHOEP data sets 3 and 25, we estimated parameters for the combination Cyclophosphamide 1400 mg/m2, Doxorubicin 32.5 mg/m2 and Vincristine 2 mg, and for Etoposide 175 mg/m2 with the constraint that parameter values must be larger than for CHOEP young.

  6. 6.

    Independently of the previous settings, parameters are determined for Doxorubicin 60 mg/m2 and Docetaxel 75 mg/m2 using data sets 1, 4–8.

  7. 7.

    Using data sets 2, 9–11, parameters are determined for the combination of Carboplatin and Paclitaxel 225 mg/m2.

  8. 8.

    Using simultaneously the data sets of E-T-C (data set 28) and EC-T (data set 19), the parameter settings for Epirubicin (dose 90 mg/m2 or 150 mg/m2), Paclitaxel (dose 175 mg/m2 or 225 mg/m2) and Cyclophosphamide (dose 600 mg/m2 or 2500 mg/m2) were determined with the constraint that lower doses have lower values of toxicity parameters.

  9. 9.

    With the ESHAP data set 29, parameter settings for Etoposide 40 mg/m2, Cytarabine 2000 mg/m2 and Cisplatin 25 mg/m2 were determined.

Three scenarios were not used for parameter estimation and served as model validation: WBC data from non Hodgkin lymphoma patients treated with CHOP-14 and Filgrastim on day 6–12 (data set 22), ANC data of patients with relapsed or persistent HD or NHL, treated with ESHAP and Pegfilgrastim 100 μg/kg on day 6 (data set 32) and WBC and G-CSF serum level data from non Hodgkin lymphoma patients treatet with CHOP-14 and Pegfilgrastim 6000 μg on day 2 (data set 12).

Quantification of myelotoxicity

In order to compare toxicity of different chemotherapy and G-CSF scenarios, it is necessary to quantify the degree of reduction of granulocytes during the course of the therapy. There is evidence that the risk for infectious complications in neutropenic patients depends on the depth as well as on the duration of neutropenia [50],[51]. Therefore, we defined the area between a certain threshold and the model curve below the threshold (AOC) as an appropriate summary measure for neutropenia/leukopenia or severity of reduction of other bone marrow cell stages. We used 2.000/μl and 4.000/μl as thresholds for total cell counts of neutrophils and leukocytes respectively. For normalized cell counts we always use the steady-state value 1 as threshold. The AOC was either used to compare overall toxicity between schedules or served as a target measure for optimizing G-CSF schedules.

Technical implementation

The model equations were programmed and solved on a standard personal computer using the numeric computation software Matlab (R2007b) and the integrated Simulink toolbox v7.0 (The MathWorks, Natick, MA). Model simulations were performed by numerical integration of the ODE system. For our model, evaluation of functions is expensive. Therefore, we used the variable step solver from Adams and Bashford (ode113).


Applying our model, we simulated 14 different chemotherapy scenarios (TA, CP, CHOP-14, CHOP-21, CHOEP-14, CHOEP-21, high CHOEP-14, high CHOEP-21, BEACOPP-14, BEACOPP-21, BEACOPP escalated, EC-T, E-T-C, and ESHAP) including 14 different drugs. Taking into account individual risk groups, we estimated a total of 12 different parameter sets. Considering different schedules of Pegfilgrastim and Filgrastim, 33 scenarios were modelled. First, we study the qualitative behaviour of our resulting chemotherapy model.

Qualitative behaviour of the chemotherapy model

In Figure 3, we studied the behaviour of our chemotherapy model on the basis of simplified chemotherapy actions. The estimated parameter set for CHOP in elderly patients was considered for this purpose. At first, the effect of an isolated stem cell kill imposed by the CHOP chemotherapy is simulated (Figure 3A). As a result, the stem cells are diminished quickly, while the later cell stages decrease with some delay. After mild oscillations, the cell counts approach normal levels.

Figure 3
figure 3

Toxicity of CHOP chemotherapy. Time courses of normalised cell counts of different cell stages after a single application of CHOP chemotherapy. A: CHOP effect only on stem cells. B: CHOP effect only on PGB. C: CHOP effect on all cell stages.

Next, we simulated an isolated CHOP chemotherapy effect on the compartment PGB alone (Figure 3B). PGB decrease immediately, and, after certain delay, other cell stages are reduced too. After oscillation, the cell counts return to normal levels similar to the isolated stem cell kill.

In Figure 3C we show the results of CHOP toxicity affecting all cell stages. This toxic effect is equivalent to later simulations of clinical scenarios involving CHOP. The figure shows that due to the combined toxicity on all cell stages, the compartment CG is most seriously affected. Recovery of the system takes much more time than in the above mentioned scenarios.

Simulations of simple chemotherapies

According to step 1 of our estimation procedure, we fitted parameters of simple chemotherapy scenarios first. Simple chemotherapies refer to those comprising either a small number of different cytotoxic drugs or drug combinations applied at the same time. This applies for data sets 1–12, 13–17, 20, 21, 23–25, 30, 31 and 33, where only one or two toxicity parameters per cell stage are required to describe the therapy. As example, a comparison of model and data for the CHOP and CHOEP young scenarios with and without G-CSF treatment can be found in Figure 4. Other scenarios can be found in the appendix, Additional file 1: Figures A3-A8. Estimated parameter sets (Table 3) resulted in a good agreement of model and data.

Figure 4
figure 4

Simulation results for CHOP-21, CHOP-14, CHOEP-21, CHOEP-14, younger patients. We show simulated cell counts for CHOP and CHOEP chemotherapy. Dots represent patient medians, grey lines represent interquartile range of patient data, black squares correspond to chemotherapy administrations, “+” correspond to days with G-CSF-injections. Clinical data originate from our collaborating clinical trials group [39], see Table 2.

Table 3 Toxicity parameters: Each drug or drug-combination is characterized by a set of eight parameters

Chemotherapy model: more complex chemotherapy simulations

In the next step, more complex chemotherapies containing a higher number of cytotoxic drugs or more complex schedules are modeled (scenarios 18, 26–29). If applicable, toxicity estimates of drugs or drug combinations established in the previous section were kept constant. Comparisons of model and data for selected scenarios can be found in Figure 5. All other scenarios can be found in the appendix (Additional file 1: Figures A9-A10).

Figure 5
figure 5

Simulation results of selected complex chemotherapies. We present results for the two breast cancer therapies, EC-T and ETC with Filgrastim on days 3–10 (first row). Note that in these schedules, chemotherapeutic drugs differ between cycles: For EC-T the drugs epirubicine and cyclophosphamide where applied in combination in the first four cycles. The single drug paclitaxel was applied for the last four cycles. For ETC, the single drug epirubicine was applied in three cycles followed by three cycles of paclitaxel and three cycles of cyclophosphamide. We also present two therapies of advanced Hodgkin’s lymphoma, BEACOPP-21 and BEACOPP escalated with Filgrastim on days 8–15, in which multiple drugs are administered at different time points per cycle (second row). Dosages can be found in Table 1. Dots represent patient medians, grey lines represent interquartile range of patient data, squares represent chemotherapy administrations, “+” denote G-CSF-injections.

Quantification of chemotherapy toxicity

As can be seen from Figures 4 and 5 and those of other scenarios presented in the appendix, our model assumptions regarding chemotherapy action and corresponding toxicity parameters resulted in a reasonable fit of almost all scenarios considered.

Our toxicity parameters can be interpreted as the strength of chemotherapy damage on the respective cell stage. We now discuss and interprete these parameters in more detail. An overview for different drugs and drug doses can be found in Table 3.

Sensitivity analysis (Additional file 1: Figures A1, A2) revealed that among bone marrow toxicities, those estimated for stem cells showed the highest precision in most scenarios. The LY toxicity and the delay parameters are sensitive too.

Figure 6 shows the relation between stem cell toxicity and resulting WBC toxicity. The correlation (Spearman) is 0.88, i.e. the stem cell toxicity is a good predictor of the overall toxicity. Additional file 1: Figure A11 of the supplement material shows the correlation of MGB AOC and WBC AOC. In Additional file 1: Figure A12 the correlation of ANC AOC and WBC AOC is depicted. Both are highly correlated too (MGB AOC vs. WBC AOC r = 0.94, ANC AOC vs. WBC AOC r = 0.92).

Figure 6
figure 6

Correlation of stem cell toxicity and peripheral toxicity. We determined cumulative toxicities of stem cells and mature blood cells for each drug or drug combination considered. Toxicities are expressed in terms of AOC of normalized cell counts applying the steady-state value 1 as threshold. AOC is calculated over 28 days. Only a single injection of chemotherapy was simulated for this purpose. The unit of AOC is “d”. A good correlation between stem cell and peripheral toxicity can be observed (Spearman’s r = 0.88). Toxicity relations between schedules are plausible. Note that schedules may be attributed to different groups of patients.

For Cyclophosphamide (with or without Doxorubicin and Vincristin) we estimated a high stem cell toxicity in agreement with the literature (e.g. [52]). For Etoposide, we obtained a rather low stem cell toxicity, even for increased dose levels, and higher toxicity to later cell stages. This also complies with the literature (e.g. [2]). Due to the rather small haematotoxic influence of Vincristine, we abstained from determining a separate parameter set for Vincristine [52].

Two anthracyclines were considered, Doxorubicine and Epirubicine. The first one was always applied in combination with other drugs, namely with Cyclophsphamide for therapies of lymphoma diseases and with Docetaxel for the TA regimen as adjuvant breast cancer therapy. Therefore, no separate parameter set of Doxorubicine alone could be derived. In contrast, Epirubicine was applied as single drug in the ETC therapy of breast cancer patients. This allows us to derive a separate set of toxicity parameters for two dose levels of Epirubicine (90 and 150 mg/m2) showing a considerable, dose-dependend stem cell toxicity in agreement with the literature [52].

According to our assumption 5, cytotoxic drugs are assumed to contribute to overall toxicity independently of each other. This does not apply for the combination of Carboplatin and Paclitaxel for which it is known that the combination is less toxic than the single drugs [53]. Therefore, a new set of toxicity parameters was determined for this drug combination, which indeed resulted in lower estimates than for Paclitaxel alone (see Table 3).

We assumed that risk groups of haematotoxicity can be traced back to differences in toxicity parameters (assumption 7). This assumption worked fine if comparing the toxicity outcomes of young and elderly patients treated with CHOP or CHOEP chemotherapies. For both risk groups the agreement of model and data is fine while corresponding toxicity parameters are higher in elderly patients.


Data sets not used for parameter fitting served as validation scenarios of our model. This requires that the corresponding chemotherapy parameters were determined on the basis of other data sets. Scenarios 12, 22, 32 fulfill these requirements. Figure 7 shows the agreement of model and data for validation scenarios.

Figure 7
figure 7

Validation scenarios. We compare model results with clinical data from validation scenarios: A: CHOP with Pegfilgrastim 6000 μg on day 2 (data: [37]), B: ESHAP with Pegfilgrastim on day 5 (data: [45]), C: CHOP with 480 μg Filgrastim on cycle-days 6–12, for elderly patients, data: [43]. Dots represent patient medians, squares represent the chemotherapy administrations, + are time points with G-CSF-injections.

Model predictions

A key feature of our model is that it allows simulations of alternative G-CSF schedules and its effects on overall leukotoxicity. We demonstrate this on the basis of the CHOP regimen for elderly patients: Using toxicity parameters estimated for CHOP based on the G-CSF schedules presented above, we modified the starting time and the duration of Filgrastim treatment. Comparisons of schedules can be performed by calculating the AOC of the simulation results. Two examples of simulated G-CSF schedules are presented in Figure 8. The regimen day 2–8 results in clearly inferior AOC than the current standard (G-CSF at day 3–12). In contrast, for the schedule day 5–13 we predict a better AOC than the current standard even though the number of injections is reduced.

Figure 8
figure 8

Analysis of CHOP with Filgrastim (age > 60 years). Simulated cell counts for CHOP-14, with Filgrastim 480 μg on days 5–13 (black line), days 2–8 (grey line) and days 3–12. Filgrastim at day 2–8 results in particularly low leukocyte counts. Better results are obtained by the Filgrastim application on days 3–12 or days 5–13.

We calculated the WBC AOC for 6 cycles of the CHOP-14 regimen, administered to elderly patients, in dependence on different Filgrastim doses and injection numbers starting on day seven (Figure 9). Best results are predicted, if Filgrastim injections are applied from day seven up to the end of the therapy cycle (see also [54]). Increasing G-CSF dose results only in marginal improvements.

Figure 9
figure 9

Modified Filgrastim schedules for CHOP-14 in elderly patients. Predicted WBC AOC (applying a threshold of 4000/μl, calculated over 84 days, unit of AOC is 1000/μl*d) under CHOP-14 with Filgrastim: G-CSF injections start at day 7 in each cycle. We modified the number of Filgrastim injections and its doses. The color scale on the right corresponds to AOC values.


Conventional cytotoxic chemotherapy plays a major role in cancer therapy. Development of intensified regimen improved the outcome of several diseases [39]-[41],[55],[56] but is limited by toxic side effects. A major, frequently dose-limiting side effect is granulotoxicity which is routinely treated with the growth factor G-CSF. A variety of pharmaceutical derivatives are available, which differ greatly in pharmacokinetic and -dynamic properties. Furthermore, outcome of growth factor treatment depends on many factors such as chemotherapy drugs used, drug doses, growth-factor derivatives and individual risk factors [57],[58]. Due to this variety of variable therapy parameters, identification of optimal growth-factor schedules cannot be performed solely on the basis of clinical trials.

We recently developed a model of the pharmacokinetic and –dynamic action of the G-CSF derivatives Filgrastim and Pegfilgrastim under conventional poly-chemotherapy [25]. We also showed that the model successfully predicts the outcome of alternative G-CSF schedules [54]. However, so far only a single simple chemotherapy schedule was considered for which the data base was most comprehensive, namely the CHOP regimen used to treat high-grade non-Hodgkin’s lymphoma diseases. The major purpose of the present work is to extend the applicability of our model considering a broad range of conventional chemotherapy schedules. This requires the construction of a comprehensive model of chemotherapy action on the granulopoietic system. Making a number of biologically plausible assumptions and translating them into differential equations allowed us to predict granulocyte and leukocyte dynamics of virtually all chemotherapy scenarios with published time series data of granulocytes and leukocytes (33 scenarios comprising 10 different chemotherapies). Modelling of chemotherapies essentially requires estimation of dose, drug and cell-stage specific toxicity parameters. In consequence, our model can easily be applied to novel chemotherapy scenarios for which time series data are available allowing estimation of these parameters. We showed how the model can be used to systematically explore the outcomes of alternative G-CSF schedules for a chemotherapy for which toxicity parameters are available.

Ongoing efforts to model haematopoiesis under chemotherapy and growth-factor applications are considerable [59]-[84]. Most newer models consider G-CSF as the major stimulant of granulopoiesis, and account for corresponding intracellular mechanisms, as well as for receptor binding kinetics and endocytic ligand depletion [62],[64],[85],[86]. Shochat et al. [63], and Foley et al. [59], proposed models considering both, stimulating effects of G-CSF as well as the cell depleting effects of chemotherapy. However, published models usually consider selected chemotherapy regimens. So far it has not been shown that these model concepts are valid for a broad range of chemotherapies and schedules [59],[63]. First attempts to predict the performance of alternative G-CSF schedules on the basis of these models were performed [87].

In order to construct a comprehensive model of chemotherapy action on granulopoiesis, we made the following assumptions and translated them into differential equations:

Delayed toxicity

It is assumed that the cell depleting effect of chemotherapy starts immediately after drug application. The maximum is reached after some time delay. This assumption is motivated by available time series data of murine bone marrow cellularity after a variety of cytotoxic drug applications often showing a maximum response to chemotherapy treatment at later time points even if the underlying drugs are quickly metabolized in vivo [52]. This phenomenon can be explained for example by delayed apoptosis of cells after damage, e.g. at time when cells entering their next cell cycle.

Cell type specific toxicity

It is assumed that chemotherapy acts cell type specific which is supported by numerous experimental data [52]. This implies that toxicity parameters are assumed to be dose, drug and cell-stage specific. Most of our data sets comprise leukocyte counts instead of neutrophil profiles. To account for this fact, we accompanied our cell-kinetic model of granulopoiesis by a simple model of lymphotoxicity. This is motivated by differing dynamics of granulocytes and lymphocytes observed in chemotherapy-treated mice [88]. Further clinical evidence is provided by a trial with breast cancer patients undergoing polychemotherapy. In this study, suppressed B- and T-cell populations were still present at times when the absolute neutrophil count had returned to normal or even higher than normal values [89]. These findings are in good agreement with the results of our parameter estimation (prolonged toxicity for lymphocytes). However, to precisely quantify lymphocyte toxicity profiles, more detailed differential blood counts of patients undergoing chemotherapy would be required than currently available.

First cycle effect

There is some evidence that the first application of chemotherapeutic drugs results in higher toxicity [81],[90]-[92]. We modelled this effect in a phenomenologic way by multiplying toxicity parameters with a factor ≥ 1 at time of first application of chemotherapy.

Toxicity of drug combinations

To estimate the overall toxicity of a drug combination, we achieved satisfying results by adding the toxicity parameters of single substances or groups of substances. However, there is evidence that some drug combinations interact in a paradoxical way, resulting in an overall toxicity that is smaller than the toxicity of either one of the single agents. For example, thrombocytopenia tends to be significantly less pronounced in patients treated with carboplatin when combined with paclitaxel (which by itself causes thrombocytopenia, too) [53],[93],[94]. The reason for the platelet sparing effect of this combination is still unknown. Therefore, it is possible that the toxicity of a drug combination cannot be simply derived by adding the toxicity parameters of their components determined in previous studies. In these cases, new parameter fittings are required.

Action of chemotherapy on granulopoiesis

We modeled chemotherapy to act cytotoxic, rather than cytostatic. As result, cells in our model were removed from the compartments directly, whereas cell kinetic properties (amplification and transit time) are not affected. However, it is well known that many drugs act cytostatically, for example by disrupting cellular metabolism or DNA synthesis. Hence, onset and severity of chemotherapy associated myelotoxicity depend on the cell cycle. However, in our modeling framework, it is neither possible nor necessary to distinguish between cytostatic and cytotoxic effects since both result in reduced cell numbers within a relatively small time frame. Closely meshed time series data of bone marrow cell stages would be required for a more detailed modelling of this issue which however cannot be established for humans.

Modelling risk groups

Numerous clinical risk factors regarding toxic response of patients are known such as age, sex and general health status [49]. We hypothesize that this heterogeneity can be traced back to different sets of toxicity parameters rather than cell kinetic parameters of granulopoiesis [28],[95]. This assumption allowed us for example to stratify patients into risk groups described by risk-specific toxicity parameters [58]. So far, differences between younger and elderly patients could be successfully explained by higher toxicity parameters as can be seen on our toxicity parameters for etoposide and the combination of cyclophosphamide, doxorubicine and vincristine applied in CHOP-like chemotherapy regimens. For both risk groups, corresponding parameter estimates resulted in good explanation of clinical data.

Overall, our model assumptions proved to be feasible for modelling almost all published time series data after a large variety of chemotherapies and schedules.

A few scenarios, for which complete sets of toxicity parameters are available, were used to successfully validate the model. We qualitatively compared derived toxicity parameters between schedules and received clinically plausible results: We estimated for example that Bleomycin (10 mg/m2), Procarbacine (100 mg/m2) and low dose Etoposide (100 mg/m2) have low granulotoxicity according to clinical experiences. In contrast, BEACOPP escalated and high-CHOEP are among the most toxic therapies in agreement with high percentages of grade 3 and 4 leucopoenia observed in these patients [43],[49].

We observed that stem cell toxicities of drugs or drug combinations correlate well with resulting leukocyte toxicity (r = 0.88, Figure 6) indicating that the parameter of stem cell toxicity is the most sensitive of our parameters characterizing chemotherapy toxicity. In contrast, the parameters for CG, PGB and MGB are less well characterized, in general. For these three cell stages, we have to acknowledge that higher toxicity at later stages can somewhat be compensated with lower toxicity at earlier cell stages and vice versa.

Finally, we demonstrated how the model could be used to make clinically relevant predictions regarding the outcome of alternative growth-factor schedules after chemotherapy. This requires that the toxicity parameters of the considered therapy are available. Then, the model can be used to simulate and compare alternative growth-factor schedules. Since hard clinical endpoints such as febrile neutropenia, use of antibiotics or length of stay in hospital cannot be addressed by our modelling, it was necessary to use surrogate markers in order to compare efficacy of G-CSF prophylaxis between schedules. We used the area between model curve and the line of 2.000/μl neutrophils or 4.000/μl leukocytes for this purpose. Other measures of relative toxicity are discussed elsewhere [95].

Finally, we have to acknowledge that the present model only allows median predictions while critical time-courses are clinically more relevant and therapy-limiting. Although this aspect is not yet covered, there is a clear perspective towards modelling individual data either by fitting parameter sets for patient risk groups or by assuming distributions of model parameters. Accordingly, we plan to extend our model and apply it in order to support improvement and individualisation of G-CSF therapies. Since we have the clear intention to apply our model in clinical contexts, we also plan public release of our model software in the near future.


We successfully developed a bio-mathematical model of granulopiesis under chemotherapy and applications of the G-CSF derivatives Filgrastim and Pegfilgrastim. Our model is able to simulate neutrophil and leukocyte profiles in the peripheral blood under various chemotherapies, with and without Filgrastim or Pegfilgrastim. The model consistently explains available clinical data, and can be used to predict the performance of alternative G-CSF schedules.

Additional file


  1. Bonadonna G, Valagussa P, Moliterni A, Zambetti M, Brambilla C: Adjuvant cyclophosphamide, methotrexate, and fluorouracil in node-positive breast cancer: the results of 20 years of follow-up. N Engl J Med. 1995, 332 (14): 901-906. 10.1056/NEJM199504063321401.

    Article  CAS  PubMed  Google Scholar 

  2. Budman DRBD, Cirrincione CT, Henderson IC, Wood WC, Weiss RB, Ferree CR, Muss HB, Green MR, Norton L, Frei E: Dose and dose intensity as determinants of outcome in the adjuvant treatment of breast cancer. The Cancer and Leukemia Group B. J Natl Cancer Inst. 1998, 90: 1205-1211. 10.1093/jnci/90.16.1205.

    Article  CAS  PubMed  Google Scholar 

  3. Kwak LW, Halpern J, Olshen RA, Horning SJ: Prognostic significance of actual dose intensity in diffuse large-cell lymphoma: results of a tree-structured survival analysis. J Clin Oncol. 1990, 8 (6): 963-977.

    CAS  PubMed  Google Scholar 

  4. Epelbaum R, Faraggi D, Ben-Arie Y, Ben-Shahar M, Haim N, Ron Y, Robinson E, Cohen Y: Survival of diffuse large cell lymphoma. Cancer. 1990, 66 (6): 1124-1129. 10.1002/1097-0142(19900915)66:6<1124::AID-CNCR2820660608>3.0.CO;2-T.

    Article  CAS  PubMed  Google Scholar 

  5. Lepage E, Gisselbrecht C, Haioun C, Sebban C, Tilly H, Bosly A, Morel P, Herbrecht R, Reyes F, Coiffier B: Prognostic significance of received relative dose intensity in non-Hodgkin’s lymphoma patients: application to LNH-87 protocol. The GELA. (Groupe d’Etude des Lymphomes de l’Adulte). Ann Oncol. 1993, 4 (8): 651-656.

    CAS  PubMed  Google Scholar 

  6. Dale DC: Colony-stimulating factors for the management of neutropenia in cancer patients. Drugs. 2002, 62 (Suppl 1): 1-15. 10.2165/00003495-200262001-00001.

    Article  CAS  PubMed  Google Scholar 

  7. Hryniuk W: Dosage parameters in chemotherapy of breast cancer. Breast Dis. 2001, 14: 21-30.

    CAS  PubMed  Google Scholar 

  8. Hryniuk WM, Goodyear M: The calculation of received dose intensity. J Clin Oncol. 1990, 8 (12): 1935-1937.

    CAS  PubMed  Google Scholar 

  9. Smith TJ, Khatcheressian J, Lyman GH, Ozer H, Armitage JO, Balducci L, Bennett CL, Cantor SB, Crawford J, Cross SJ, Demetri G, Desch CE, Pizzo PA, Schiffer CA, Schwartzberg L, Somerfield MR, Somlo G, Wade JC, Wade JL, Winn RJ, Wozniak AJ, Wolff AC: 2006 update of recommendations for the use of white blood cell growth factors: An evidence-based clinical practice guideline. J Clin Oncol. 2006, 24 (19): 3187-3205. 10.1200/JCO.2006.06.4451.

    Article  CAS  PubMed  Google Scholar 

  10. Gisselbrecht C, Haioun C, Lepage E, Bastion Y, Tilly H, Bosly A, Dupriez B, Marit G, Herbrecht R, Deconinck E, Marolleau JP, Yver A, Dabouz-Harrouche F, Coiffier B, Reyes F: Placebo-controlled phase III study of lenograstim (glycosylated recombinant human granulocyte colony-stimulating factor) in aggressive non-Hodgkin’s lymphoma: factors influencing chemotherapy administration. Groupe d’Etude des Lymphomes de l’Adulte. Leuk Lymphoma. 1997, 25 (3–4): 289-300.

    CAS  PubMed  Google Scholar 

  11. Sonneveld P, de Ridder M, van der Lelie H, Nieuwenhuis K, Schouten H, Mulder A, van Reijswoud I, Hop W, Lowenberg B: Comparison of doxorubicin and mitoxantrone in the treatment of elderly patients with advanced diffuse non-Hodgkin’s lymphoma using CHOP versus CNOP chemotherapy. J Clin Oncol. 1995, 13: 2530-2539.

    CAS  PubMed  Google Scholar 

  12. Tirelli U: CHOP is the standard regimen in patients > or = 70 years of age with intermediate-grade and high-grade non-Hodgkin’s lymphoma: results of a randomized study of the European Organization for Research and Treatment of Cancer Lymphoma Cooperative Study Group. J Clin Oncol. 1998, 16: 27-34.

    CAS  PubMed  Google Scholar 

  13. Osby E, Hagberg H, Kvaløy S, Teerenhovi L, Anderson H, Cavallin-Stahl E, Holte H, Myhre J, Pertovaara H, Björkholm M: CHOP is superior to CNOP in elderly patients with aggressive lymphoma while outcome is unaffected by filgrastim treatment: results of a Nordic Lymphoma Group randomized trial. Blood. 2003, 101: 3840-3848. 10.1182/blood-2002-10-3238.

    Article  CAS  PubMed  Google Scholar 

  14. Meyer RM, Browman GP, Samosh ML, Benger AM, Bryant-Lukosius D, Wilson WE, Frank GL, Leber BF, Sternbach MS, Foster GA: Randomized phase II comparison of standard CHOP with weekly CHOP in elderly patients with non-Hodgkin’s lymphoma. J Clin Oncol. 1995, 13 (9): 2386-2393.

    CAS  PubMed  Google Scholar 

  15. Marangolo M, Bengala C, Conte PF, Danova M, Pronzato P, Rosti G, Sagrada P: Dose and outcome: the hurdle of neutropenia (Review). Oncol Rep. 2006, 16 (2): 233-248.

    CAS  PubMed  Google Scholar 

  16. Briones MA, Josephson CD, Hillyer CD: Granulocyte transfusion: revisited. Curr Hematol Rep. 2003, 2 (6): 522-527.

    PubMed  Google Scholar 

  17. Dale D: Current management of chemotherapy-induced neutropenia: the role of colony-stimulating factors. Semin Oncol. 2003, 30 (4 Suppl 13): 3-9. 10.1016/S0093-7754(03)00326-9.

    Article  CAS  PubMed  Google Scholar 

  18. Mitchell PL, Morland B, Stevens MC, Dick G, Easlea D, Meyer LC, Pinkerton CR: Granulocyte colony-stimulating factor in established febrile neutropenia: a randomized study of pediatric patients. J Clin Oncol. 1997, 15 (3): 1163-1170.

    CAS  PubMed  Google Scholar 

  19. Lyman GH, Kuderer NM: Filgrastim in patients with neutropenia: potential effects on quality of life. Drugs. 2002, 62 (Suppl 1): 65-78. 10.2165/00003495-200262001-00005.

    Article  CAS  PubMed  Google Scholar 

  20. Garcia-Carbonero R, Mayordomo JI, Tornamira MV, Lopez-Brea M, Rueda A, Guillem V, Arcediano A, Yubero A, Ribera F, Gomez C, Trés A, Pérez-Gracia JL, Lumbreras C, Hornedo J, Cortés-Funes H, Paz-Ares L: Granulocyte colony-stimulating factor in the treatment of high-risk febrile neutropenia: a multicenter randomized trial. J Natl Cancer Inst. 2001, 93 (1): 31-38. 10.1093/jnci/93.1.31.

    Article  CAS  PubMed  Google Scholar 

  21. Maher DW, Lieschke GJ, Green M, Bishop J, Stuart-Harris R, Wolf M, Sheridan WP, Kefford RF, Cebon J, Olver I, McKendrick J, Toner G, Bradstock K, Lieschke M, Cruickshank S, Tomita DK, Hoffman EW, Fox RM, Morstyn G: Filgrastim in patients with chemotherapy-induced febrile neutropenia. A double-blind, placebo-controlled trial. Ann Intern Med. 1994, 121 (7): 492-501. 10.7326/0003-4819-121-7-199410010-00004.

    Article  CAS  PubMed  Google Scholar 

  22. Siena S, Secondino S, Giannetta L, Carminati O, Pedrazzoli P: Optimising management of neutropenia and anaemia in cancer chemotherapy-advances in cytokine therapy. Crit Rev Oncol Hematol. 2003, 48 (Suppl): S39-47. 10.1016/j.critrevonc.2003.05.002.

    Article  PubMed  Google Scholar 

  23. Schmitz S, Franke H, Brusis J, Wichmann HE: Quantification of the Cell Kinetic Effects of G-Csf Using a Model of Human Granulopoiesis. Exp Hematol. 1993, 21 (6): 755-760.

    CAS  PubMed  Google Scholar 

  24. Lord BI, Bronchud MH, Owens S, Chang J, Howell A, Souza L, Dexter TM: The kinetics of human granulopoiesis following treatment with granulocyte colony-stimulating factor in vivo. Proc Natl Acad Sci U S A. 1989, 86 (23): 9499-9503. 10.1073/pnas.86.23.9499.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Scholz M, Schirm S, Wetzler M, Engel C, Loeffler M: Pharmacokinetic and -dynamic modelling of G-CSF derivatives in humans. Theor Biol Med Model. 2012, 9: 32-10.1186/1742-4682-9-32.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Wetzler M: Mathematisches Kompartimentmodell der menschlichen Granulopoese zur Beschreibung der pharmakokinetischen und pharmakodynamischen Eigenschaften der G-CSF-Derivate Filgrastim und Pegfilgrastim unter Chemotherapie. 2012, Leipzig, Leipzig

    Google Scholar 

  27. Schirm S, Engel C, Loeffler M, Scholz M: A biomathematical model of human erythropoiesis under erythropoietin and chemotherapy administration. PLoS One. 2013, 8 (6): e65630-10.1371/journal.pone.0065630.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Chatta GS, Price TH, Allen RC, Dale DC: Effects of in vivo recombinant methionyl human granulocyte colony-stimulating factor on the neutrophil response and peripheral blood colony-forming cells in healthy young and elderly adult volunteers. Blood. 1994, 84 (9): 2923-2929.

    CAS  PubMed  Google Scholar 

  29. Bishop CR, Athens JW, Boggs DR, Warner HR, Cartwright GE, Wintrobe MM: Leukokinetic studies. 13. A non-steady-state kinetic evaluation of the mechanism of cortisone-induced granulocytosis. J Clin Invest. 1968, 47 (2): 249-260. 10.1172/JCI105721.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Dale DC, Fauci AS, Wolff SM: Alternate-day prednisone. Leukocyte kinetics and susceptibility to infections. N Engl J Med. 1974, 291 (22): 1154-1158. 10.1056/NEJM197411282912203.

    Article  CAS  PubMed  Google Scholar 

  31. Karle H, Hansen NE: Neutrophil kinetics in corticosteroid induced neutrophilia measured by plasma lysozyme. Scand J Haematol. 1975, 14 (3): 190-195. 10.1111/j.1600-0609.1975.tb00316.x.

    Article  CAS  PubMed  Google Scholar 

  32. Holmes FA, Jones SE, O’Shaughnessy J, Vukelja S, George T, Savin M, Richards D, Glaspy J, Meza L, Cohen G, Dhami M, Budman DR, Hackett J, Brassard M, Yang BB, Liang BC: Comparable efficacy and safety profiles of once-per-cycle pegfilgrastim and daily injection filgrastim in chemotherapy-induced neutropenia: a multicenter dose-finding study in women with breast cancer. Ann Oncol. 2002, 13 (6): 903-909. 10.1093/annonc/mdf130.

    Article  CAS  PubMed  Google Scholar 

  33. Johnston E, Crawford J, Blackwell S, Bjurstrom T, Lockbaum P, Roskos L, Yang BB, Gardner S, Miller-Messana MA, Shoemaker D, Garst J, Schwab G: Randomized, dose-escalation study of SD/01 compared with daily filgrastim in patients receiving chemotherapy. J Clin Oncol. 2000, 18 (13): 2522-2528.

    CAS  PubMed  Google Scholar 

  34. Trumper L, Zwick C, Ziepert M, Hohloch K, Schmits R, Mohren M, Liersch R, Bentz M, Graeven U, Wruck U, Hoffmann M, Metzner B, Hasenclever D, Loeffler M, Pfreundschuh M: Dose-escalated CHOEP for the treatment of young patients with aggressive non-Hodgkin’s lymphoma: I. A randomized dose escalation and feasibility study with bi- and tri-weekly regimens. Ann Oncol. 2008, 19 (3): 538-544. 10.1093/annonc/mdm497.

    Article  CAS  PubMed  Google Scholar 

  35. Zamboni WC: Pharmacokinetics of pegfilgrastim. Pharmacotherapy. 2003, 23 (8 Pt 2): 9S-14S. 10.1592/phco.23.9.9S.32888.

    Article  CAS  PubMed  Google Scholar 

  36. Yowell SL, Blackwell S: Novel effects with polyethylene glycol modified pharmaceuticals. Cancer Treat Rev. 2002, 28: 3-6. 10.1016/S0305-7372(02)80002-0.

    Article  CAS  PubMed  Google Scholar 

  37. George S, Yunus F, Case D, Yang BB, Hackett J, Shogan JE, Meza LA, Neumann TA, Liang BC: Fixed-dose pegfilgrastim is safe and allows neutrophil recovery in patients with non-Hodgkin’s lymphoma. Leuk Lymphoma. 2003, 44 (10): 1691-1696. 10.1080/1042819031000063462.

    Article  CAS  PubMed  Google Scholar 

  38. Mey UJ, Maier A, Schmidt-Wolf IG, Ziske C, Forstbauer H, Banat GA, Reber M, Strehl JW, Gorschlueter M: Pegfilgrastim as hematopoietic support for dose-dense chemoimmunotherapy with R-CHOP-14 as first-line therapy in elderly patients with diffuse large B cell lymphoma. Support Care Cancer. 2007, 15 (7): 877-884. 10.1007/s00520-006-0201-z.

    Article  PubMed  Google Scholar 

  39. Pfreundschuh M, Trumper L, Kloess M, Schmits R, Feller AC, Rudolph C, Reiser M, Hossfeld DK, Metzner B, Hasenclever D, Schmitz N, Glass B, Rübe C, Loeffler M: Two-weekly or 3-weekly CHOP chemotherapy with or without etoposide for the treatment of young patients with good-prognosis (normal LDH) aggressive lymphomas: results of the NHL-B1 trial of the DSHNHL. Blood. 2004, 104 (3): 626-633. 10.1182/blood-2003-06-2094.

    Article  CAS  PubMed  Google Scholar 

  40. Pfreundschuh M, Trumper L, Kloess M, Schmits R, Feller AC, Rube C, Rudolph C, Reiser M, Hossfeld DK, Eimermacher H, Hasenclever D, Schmitz N, Loeffler M: Two-weekly or 3-weekly CHOP chemotherapy with or without etoposide for the treatment of elderly patients with aggressive lymphomas: results of the NHL-B2 trial of the DSHNHL. Blood. 2004, 104 (3): 634-641. 10.1182/blood-2003-06-2095.

    Article  CAS  PubMed  Google Scholar 

  41. Diehl V, Franklin J, Pfreundschuh M, Lathan B, Paulus U, Hasenclever D, Tesch H, Herrmann R, Dörken B, Müller-Hermelink HK, Dühmke E, Loeffler M: Standard and increased-dose BEACOPP chemotherapy compared with COPP-ABVD for advanced Hodgkin’s disease. New Engl J Med 2003, 348(24):2386–2395.,

  42. Möbus VJC, Lück HJ: Intense dose-dense sequential chemotherapy with epirubicin, paclitaxel and cyclophosphamide compared with conventionally scheduled chemotherapy in high-risk primary breast cancer (4+ LN): mature results of an AGO-phase-III study. J Clin Oncol. 2010, 28: 2874-2880. 10.1200/JCO.2009.24.7643.

    Article  Google Scholar 

  43. Pfreundschuh M, Schubert J, Ziepert M, Schmits R, Mohren M, Lengfelder E, Reiser M, Nickenig C, Clemens M, Peter N, Bokemeyer C, Eimermacher H, Ho A, Hoffmann M, Mertelsmann R, Trümper L, Balleisen L, Liersch R, Metzner B, Hartmann F, Glass B, Poeschel V, Schmitz N, Ruebe C, Feller AC, Loeffler M: Six versus eight cycles of bi-weekly CHOP-14 with or without rituximab in elderly patients with aggressive CD20+ B-cell lymphomas: a randomised controlled trial (RICOVER-60). Lancet Oncol. 2008, 9 (2): 105-116. 10.1016/S1470-2045(08)70002-0.

    Article  CAS  PubMed  Google Scholar 

  44. Sieber M, Bredenfeld H, Josting A, Reineke T, Rueffer U, Koch T, Naumann R, Boissevain F, Koch P, Worst P, Soekler M, Eich H, Müller-Hermelink HK, Franklin J, Paulus U, Wolf J, Engert A, Diehl V: 14-day variant of the bleomycin, etoposide, doxorubicin, cyclophosphamide, vincristine, procarbazine, and prednisone regimen in advanced-stage Hodgkin’s lymphoma: Results of a pilot study of the German Hodgkin’s Lymphoma Study Group. J Clin Oncol. 2003, 21 (9): 1734-1739. 10.1200/JCO.2003.06.028.

    Article  CAS  PubMed  Google Scholar 

  45. Vose JM, Crump M, Lazarus H, Emmanouilides C, Schenkein D, Moore J, Frankel S, Flinn I, Lovelace W, Hackett J, Liang BC: Randomized, multicenter, open-label study of pegfilgrastim compared with daily filgrastim after chemotherapy for lymphoma. J Clin Oncol. 2003, 21 (3): 514-519. 10.1200/JCO.2003.03.040.

    Article  CAS  PubMed  Google Scholar 

  46. Zwick C, Hartmann F, Zeynalova S, Poschel V, Nickenig C, Reiser M, Lengfelder E, Peter N, Schlimok G, Schubert J, Schmitz N, Loeffler M, Pfreundschuh M: Randomized comparison of pegfilgrastim day 4 versus day 2 for the prevention of chemotherapy-induced leukocytopenia. Ann Oncol. 2011, 22 (8): 1872-1877. 10.1093/annonc/mdq674.

    Article  CAS  PubMed  Google Scholar 

  47. Brusamolino E, Rusconi C, Montalbetti L, Gargantini L, Uziel L, Pinotti G, Fava S, Rigacci L, Pagnucco G, Pascutto C, Morra E, Lazzarino M: Dose-dense R-CHOP-14 supported by pegfilgrastim in patients with diffuse large B-cell lymphoma: a phase II study of feasibility and toxicity. Haematologica. 2006, 91 (4): 496-502.

    CAS  PubMed  Google Scholar 

  48. Rechenberg I: Evolutionsstrategie 94. 1994, Frommann-Holzboog, Stuttgart

    Google Scholar 

  49. Wunderlich A, Kloess M, Reiser M, Rudolph C, Truemper L, Bittner S, Schmalenberg H, Schmits R, Pfreundschuh M, Loeffler M: Practicability and acute haematological toxicity of 2- and 3-weekly CHOP and CHOEP chemotherapy for aggressive non-Hodgkin’s lymphoma: results from the NHL-B trial of the German High-Grade Non-Hodgkin’s Lymphoma Study Group (DSHNHL). Ann Oncol. 2003, 14 (6): 881-893. 10.1093/annonc/mdg249.

    Article  CAS  PubMed  Google Scholar 

  50. Bodey GP: Infection in cancer patients, A continuing association. Am J Med. 1986, 81 (1A): 11-26. 10.1016/0002-9343(86)90510-3.

    Article  CAS  PubMed  Google Scholar 

  51. Bodey GP, Buckley M, Sathe YS, Freireic E: Quantitative Relationships between Circulating Leukocytes and Infection in Patients with Acute Leukemia. Ann Intern Med. 1966, 328 (2): 328-340. 10.7326/0003-4819-64-2-328.

    Article  Google Scholar 

  52. Lohrmann H-P, Schreml W: Cytotoxic Drugs and the Granulopoietic System. 1982, Springer, Berlin

    Book  Google Scholar 

  53. Guminski A, Harnett P, de Fazio A: Carboplatin and paclitaxel interact antagonistically in a megakaryoblast cell line–a potential mechanism for paclitaxel-mediated sparing of carboplatininduced thrombocytopenia. Cancer Chemother Pharmacol. 2001, 48 (3): 229-334. 10.1007/s002800100279.

    Article  CAS  PubMed  Google Scholar 

  54. Zeynalova S, Ziepert M, Scholz M, Schirm S, Zwick C, Pfreundschuh M, Loeffler M, German High-Grade Non-Hodgkin Lymphoma Study G: Comparison and modelling of pegylated or unpegylated G-CSF schedules in CHOP-14 regimen of elderly patients with aggressive B-cell lymphoma. Ann Hematol. 2013, 92 (12): 1641-1652. 10.1007/s00277-013-1842-x.

    Article  CAS  PubMed  Google Scholar 

  55. Untch M, Fasching PA, Konecny GE, von Koch F, Conrad U, Fett W, Kurzeder C, Luck HJ, Stickeler E, Urbaczyk H, Liedtke B, Salat C, Harbeck N, Müller V, Schmidt M, Hasmüller S, Lenhard M, Schuster T, Nekljudova V, Lebeau A, Loibl S, von Minckwitz G: PREPARE trial: a randomized phase III trial comparing preoperative, dose-dense, dose-intensified chemotherapy with epirubicin, paclitaxel and CMF versus a standard-dosed epirubicin/cyclophosphamide followed by paclitaxel +/− darbepoetin alfa in primary breast cancer--results at the time of surger. Ann Oncol. 2011, 22 (9): 1988-1998. 10.1093/annonc/mdq709.

    Article  CAS  PubMed  Google Scholar 

  56. Untch M, von Minckwitz G, Konecny GE, Conrad U, Fett W, Kurzeder C, Luck HJ, Stickeler E, Urbaczyk H, Liedtke B, Beckmann MW, Salat C, Harbeck N, Müller V, Schmidt M, Hasmüller S, Lenhard M, Nekljudova V, Lebeau A, Loibl S, Fasching PA: PREPARE trial: a randomized phase III trial comparing preoperative, dose-dense, dose-intensified chemotherapy with epirubicin, paclitaxel, and CMF versus a standard-dosed epirubicin-cyclophosphamide followed by paclitaxel with or without darbepoetin alfa in primary breast cancer--outcome on prognosi. Ann Oncol. 2011, 22 (9): 1999-2006. 10.1093/annonc/mdq713.

    Article  CAS  PubMed  Google Scholar 

  57. Ziepert M, Hasenclever D, Kuhnt E, Glass B, Schmitz N, Pfreundschuh M, Loeffler M: Standard International prognostic index remains a valid predictor of outcome for patients with aggressive CD20+ B-cell lymphoma in the rituximab era. J Clin Oncol. 2010, 28 (14): 2373-2380. 10.1200/JCO.2009.26.2493.

    Article  CAS  PubMed  Google Scholar 

  58. Ziepert M, Schmits R, Trumper L, Pfreundschuh M, Loeffler M, German High-Grade Non-Hodgkin’s Lymphoma Study G: Prognostic factors for hematotoxicity of chemotherapy in aggressive non-Hodgkin’s lymphoma. Ann Oncol. 2008, 19 (4): 752-762. 10.1093/annonc/mdm541.

    Article  CAS  PubMed  Google Scholar 

  59. Foley C, Mackey MC: Mathematical model for G-CSF administration after chemotherapy. J Theor Biol. 2009, 257 (1): 27-44. 10.1016/j.jtbi.2008.09.043.

    Article  CAS  PubMed  Google Scholar 

  60. Ostby I, Kvalheim G, Rusten LS, Grottum P: Mathematical modeling of granulocyte reconstitution after high-dose chemotherapy with stem cell support: effect of post-transplant G-CSF treatment. J Theor Biol. 2004, 231 (1): 69-83. 10.1016/j.jtbi.2004.05.010.

    Article  CAS  PubMed  Google Scholar 

  61. Ostby I, Rusten LS, Kvalheim G, Grottum P: A mathematical model for reconstitution of granulopoiesis after high dose chemotherapy with autologous stem cell transplantation. J Math Biol. 2003, 47 (2): 101-136. 10.1007/s00285-003-0198-6.

    Article  PubMed  Google Scholar 

  62. Roskos LK, Lum P, Lockbaum P, Schwab G, Yang BB: Pharmacokinetic/pharmacodynamic modeling of pegfilgrastim in healthy subjects. J Clin Pharmacol. 2006, 46 (7): 747-757. 10.1177/0091270006288731.

    Article  CAS  PubMed  Google Scholar 

  63. Shochat E, Rom-Kedar V, Segel LA: G-CSF control of neutrophils dynamics in the blood. Bull Math Biol. 2007, 69 (7): 2299-2338. 10.1007/s11538-007-9221-1.

    Article  CAS  PubMed  Google Scholar 

  64. Vainstein V, Ginosar Y, Shoham M, Ranmar DO, Ianovski A, Agur Z: The complex effect of granulocyte colony-stimulating factor on human granulopoiesis analyzed by a new physiologically-based mathematical model. J Theor Biol. 2005, 234 (3): 311-327. 10.1016/j.jtbi.2004.11.026.

    Article  CAS  PubMed  Google Scholar 

  65. Wang B, Ludden TM, Cheung EN, Schwab GG, Roskos LK: Population pharmacokinetic-pharmacodynamic modeling of filgrastim (r-metHuG-CSF) in healthy volunteers. J Pharmacokinet Pharmacodyn. 2001, 28 (4): 321-342. 10.1023/A:1011534529622.

    Article  CAS  PubMed  Google Scholar 

  66. Wichmann HE, Loeffler M, Herkenrath P, Gerhardts MD, Wesselborg C, Wulff H: Mathematical models in hematology. Klin Wochenschr. 1983, 61 (19): 935-940. 10.1007/BF01550265.

    Article  CAS  PubMed  Google Scholar 

  67. Wichmann HE: Mathematical Modeling of Cell Proliferation: Stem Cell Regulation in Hemopoiesis. Vol 1: Model Description, Irradiation, Erythropoietic Stimulation. In., vol. 1. 1985, Boca Raton, Florida, Crc Pr I Llc

    Google Scholar 

  68. Meyer G: Charakterisierung der zellkinetischen Wirkungen bei exogener Applikation von Erythropoetin auf die Erythropoese des Menschen mit Hilfe eines mathematischen Kompartimentmodells. 2004, Leipzig University, Leipzig

    Google Scholar 

  69. Banks HT, Cole CE, Schlosser PM, Tran HT: Modeling and optimal regulation of erythropoiesis subject to benzene intoxication. Math Biosci Eng. 2004, 1 (1): 15-48. 10.3934/mbe.2004.1.15.

    Article  CAS  PubMed  Google Scholar 

  70. Belair J, Mackey MC, Mahaffy JM: Age-structured and two-delay models for erythropoiesis. Math Biosci. 1995, 128 (1–2): 317-346. 10.1016/0025-5564(94)00078-E.

    Article  CAS  PubMed  Google Scholar 

  71. Brooks G, Provencher G, Lei J, Mackey MC: Neutrophil dynamics after chemotherapy and G-CSF: the role of pharmacokinetics in shaping the response. J Theor Biol. 2012, 315: 97-109. 10.1016/j.jtbi.2012.08.028.

    Article  CAS  PubMed  Google Scholar 

  72. Colijn C, Foley C, Mackey MC: G-CSF treatment of canine cyclical neutropenia: a comprehensive mathematical model. Exp Hematol. 2007, 35 (6): 898-907. 10.1016/j.exphem.2007.02.015.

    Article  CAS  PubMed  Google Scholar 

  73. Colijn C, Mackey MC: A mathematical model of hematopoiesis–I. Periodic chronic myelogenous leukemia. J Theor Biol. 2005, 237 (2): 117-132. 10.1016/j.jtbi.2005.03.033.

    Article  PubMed  Google Scholar 

  74. Colijn C, Mackey MC: A mathematical model of hematopoiesis: II. Cyclical neutropenia. J Theor Biol. 2005, 237 (2): 133-146. 10.1016/j.jtbi.2005.03.034.

    Article  PubMed  Google Scholar 

  75. Crauste F, Demin I, Gandrillon O, Volpert V: Mathematical study of feedback control roles and relevance in stress erythropoiesis. J Theor Biol. 2010, 263 (3): 303-316. 10.1016/j.jtbi.2009.12.026.

    Article  PubMed  Google Scholar 

  76. Engel C, Scholz M, Loeffler M: A computational model of human granulopoiesis to simulate the hematotoxic effects of multicycle polychemotherapy. Blood. 2004, 104 (8): 2323-2331. 10.1182/blood-2004-01-0306.

    Article  CAS  PubMed  Google Scholar 

  77. Glauche I, Cross M, Loeffler M, Roeder I: Lineage specification of hematopoietic stem cells: mathematical modeling and biological implications. Stem Cells. 2007, 25 (7): 1791-1799. 10.1634/stemcells.2007-0025.

    Article  CAS  PubMed  Google Scholar 

  78. Krinner A, Roeder I, Loeffler M, Scholz M: Merging concepts - coupling an agent-based model of hematopoietic stem cells with an ODE model of granulopoiesis. BMC Syst Biol. 2013, 7: 117-10.1186/1752-0509-7-117.

    Article  PubMed Central  PubMed  Google Scholar 

  79. Krzyzanski W, Jusko WJ, Wacholtz MC, Minton N, Cheung WK: Pharmacokinetic and pharmacodynamic modeling of recombinant human erythropoietin after multiple subcutaneous doses in healthy subjects. Eur J Pharm Sci. 2005, 26 (3–4): 295-306. 10.1016/j.ejps.2005.06.010.

    Article  CAS  PubMed  Google Scholar 

  80. Obeyesekere MN, Berry RW, Spicer PP, Korbling M: A mathematical model of haemopoiesis as exemplified by CD34 cell mobilization into the peripheral blood. Cell Prolif. 2004, 37 (4): 279-294. 10.1111/j.1365-2184.2004.00312.x.

    Article  CAS  PubMed  Google Scholar 

  81. Scholz M, Engel C, Loeffler M: Modelling human granulopoiesis under poly-chemotherapy with G-CSF support. J Math Biol. 2005, 50 (4): 397-439. 10.1007/s00285-004-0295-1.

    Article  CAS  PubMed  Google Scholar 

  82. Scholz M, Gross A, Loeffler M: A biomathematical model of human thrombopoiesis under chemotherapy. J Theor Biol. 2010, 264 (2): 287-300. 10.1016/j.jtbi.2009.12.032.

    Article  CAS  PubMed  Google Scholar 

  83. Zhuge C, Lei J, Mackey MC: Neutrophil dynamics in response to chemotherapy and G-CSF. J Theor Biol. 2012, 293: 111-120. 10.1016/j.jtbi.2011.10.017.

    Article  CAS  PubMed  Google Scholar 

  84. Foley C, Bernard S, Mackey MC: Cost-effective G-CSF therapy strategies for cyclical neutropenia: mathematical modelling based hypotheses. J Theor Biol. 2006, 238 (4): 754-763. 10.1016/j.jtbi.2005.06.021.

    Article  CAS  PubMed  Google Scholar 

  85. Sarkar CA, Lauffenburger DA: Cell-level pharmacokinetic model of granulocyte colony-stimulating factor: implications for ligand lifetime and potency in vivo. Mol Pharmacol. 2003, 63 (1): 147-158. 10.1124/mol.63.1.147.

    Article  CAS  PubMed  Google Scholar 

  86. Wiczling P, Lowe P, Pigeolet E, Ludicke F, Balser S, Krzyzanski W: Population pharmacokinetic modelling of filgrastim in healthy adults following intravenous and subcutaneous administrations. Clin Pharmacokinet. 2009, 48 (12): 817-826. 10.2165/11318090-000000000-00000.

    Article  CAS  PubMed  Google Scholar 

  87. Shochat E, Rom-Kedar V: Novel strategies for granulocyte colony-stimulating factor treatment of severe prolonged neutropenia suggested by mathematical modeling. Clin Cancer Res. 2008, 14 (20): 6354-6363. 10.1158/1078-0432.CCR-08-0807.

    Article  CAS  PubMed  Google Scholar 

  88. Scholz M, Ackermann M, Emmrich F, Loeffler M, Kamprad M: Effectiveness of cytopenia prophylaxis for different filgrastim and pegfilgrastim schedules in a chemotherapy mouse model. Biologics. 2009, 2009 (3): 27-37.

    Google Scholar 

  89. Talmadge JE, Jackson JD, Borgeson CD, Perry GA: Differential Recovery of Polymorphonuclear Neutrophils, B-Cell and T-Cell Subpopulations in the Thymus, Bone-Marrow, Spleen and Blood of Mice Following Split-Dose Polychemotherapy. Cancer Immunol Immun. 1994, 39 (1): 59-67. 10.1007/BF01517182.

    Article  CAS  Google Scholar 

  90. Cullen MH, Billingham LJ, Gaunt CH, Steven NM: Rational selection of patients for antibacterial prophylaxis after chemotherapy. J Clin Oncol. 2007, 25 (30): 4821-4828. 10.1200/JCO.2006.08.7395.

    Article  PubMed  Google Scholar 

  91. Ozer H: The timing of chemotherapy-induced neutropenia and its clinical and economic impact. Oncology. 2006, 20: 11-15.

    PubMed  Google Scholar 

  92. Rabinowitz AP, Weiner NJ, Tronic BS, Fridman M, Liberman RF, Delgado DJ: Severe neutropenia in CHOP occurs most frequently in cycle 1: A predictive model. Leuk Lymphoma. 2006, 47 (5): 853-858. 10.1080/10428190500404316.

    Article  PubMed  Google Scholar 

  93. Fujiwara K: The platelet-sparing effect of paclitaxel is not related to changes in the pharmacokinetics of carboplatin. Cancer Chemother Pharmacol. 2001, 47: 22-26. 10.1007/s002800000212.

    Article  CAS  PubMed  Google Scholar 

  94. Loi S, Rischin D, Michael M, Yuen K, Stokes K, Ellis AG, Millward MJ, Webster LK: A randomized cross-over trial to determine the effect of Cremophor EL on the pharmacodynamics and pharmacokinetics of carboplatin chemotherapy. Cancer Chemoth Pharm. 2004, 54 (5): 407-414. 10.1007/s00280-004-0792-3.

    Article  CAS  Google Scholar 

  95. Scholz M, Engel C, Loeffler M: German High Grade Non-Hodgkin Lymphoma Study G. Br J Haematol. 2006, 132 (6): 723-735. 10.1111/j.1365-2141.2005.05957.x.

    Article  CAS  PubMed  Google Scholar 

Download references


This research was supported by the BMBF-grant on Medical Systems Biology “HaematoSys” (BMBF-FKZ 0315452). We further thank the German Hodgkin’s Lymphoma Study Group (Volker Diehl, chairman), the German High-Grade Non-Hodgkin’s Lymphoma Study Group (Michael Pfreundschuh, chairman) and the German Breast Group (Gunther von Minckwitz, chairman) for kind permission to use clinical trial data.

SS, MS were funded by the Leipzig Interdisciplinary Research Cluster of Genetic Factors, Clinical Phenotypes and Environment (LIFE Center, Universität Leipzig). LIFE is funded by means of the European Union, by the European Regional Development Fund (ERFD), the European Social Fund and by means of the Free State of Saxony within the framework of the excellence initiative.

We acknowledge support from the German Research Foundation (DFG) and Leipzig University within the program of Open Access Publishing.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Markus Scholz.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Developed the model: SS, CE, MS. Conceived and designed model simulations: SS, MS. Performed model simulations: SS. Wrote the article: SS, MS. Contributed to paper writing: CE. Contributed to discussion: ML. Final approval of manuscript: all authors.

Electronic supplementary material


Additional file 1: Modelling chemotherapy effects on granulopoiesis: Supplement Material. The file GraPaper-2-Appendix.pdf contains sensitivity analysis and further simulation results. (PDF 267 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Schirm, S., Engel, C., Loeffler, M. et al. Modelling chemotherapy effects on granulopoiesis. BMC Syst Biol 8, 138 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: