# Pharmacokinetic and -dynamic modelling of G-CSF derivatives in humans

- Markus Scholz
^{1, 2}Email author, - Sibylle Schirm
^{1}, - Marcus Wetzler
^{1}, - Christoph Engel
^{1, 2}and - Markus Loeffler
^{1, 2}

**9**:32

https://doi.org/10.1186/1742-4682-9-32

© Scholz et al.; licensee BioMed Central Ltd. 2012

**Received: **17 February 2012

**Accepted: **12 June 2012

**Published: **30 July 2012

## Abstract

### Background

The human granulocyte colony-stimulating factor (G-CSF) is routinely applied to support recovery of granulopoiesis during the course of cytotoxic chemotherapies. However, optimal use of the drug is largely unknown. We showed in the past that a biomathematical compartment model of human granulopoiesis can be used to make clinically relevant predictions regarding new, yet untested chemotherapy regimen. In the present paper, we aim to extend this model by a detailed pharmacokinetic and -dynamic modelling of two commonly used G-CSF derivatives Filgrastim and Pegfilgrastim.

### Results

Model equations are based on our physiological understanding of the drugs which are delayed absorption of G-CSF when applied to the subcutaneous tissue, dose-dependent bioavailability, unspecific first order elimination, specific elimination in dependence on granulocyte counts and reversible protein binding. Pharmacokinetic differences between Filgrastim and Pegfilgrastim were modelled as different parameter sets. Our former cell-kinetic model of granulopoiesis was essentially preserved, except for a few additional assumptions and simplifications. We assumed a delayed action of G-CSF on the bone marrow, a delayed action of chemotherapy and differences between Filgrastim and Pegfilgrastim with respect to stimulation potency of the bone marrow. Additionally, we incorporated a model of combined action of Pegfilgrastim and Filgrastim or endogenous G-CSF which interact via concurrent receptor binding. Unknown pharmacokinetic or cell-kinetic parameters were determined by fitting the predictions of the model to available datasets of G-CSF applications, chemotherapy applications or combinations of it. Data were either extracted from the literature or were received from cooperating clinical study groups. Model predictions fitted well to both, datasets used for parameter estimation and validation scenarios as well. A unique set of parameters was identified which is valid for all scenarios considered. Differences in pharmacokinetic parameter estimates between Filgrastim and Pegfilgrastim were biologically plausible throughout.

### Conclusion

We conclude that we established a comprehensive biomathematical model to explain the dynamics of granulopoiesis under chemotherapy and applications of two different G-CSF derivatives. We aim to apply the model to a large variety of chemotherapy regimen in the future in order to optimize corresponding G-CSF schedules or to individualize G-CSF treatment according to the granulotoxic risk of a patient.

## Keywords

## Introduction and background

### Background

The human granulocyte colony-stimulating factor (G-CSF) is routinely applied in various cancer chemotherapy regimen in order to ameliorate or prevent neutropenia caused by the unspecific toxicity of the drugs used[1–4]. G-CSF proved to be highly potent in stimulating granulopoiesis via several modes of action such as mitotic activation of granulopoietic progenitors and precursors, accelerated maturation of bone marrow cell stages and increased release of mature bone marrow cells[5–10].

In case of conventional (non-myeloablative) chemotherapies, the haematopoietic system usually recovers without further medication. But, G-CSF can significantly speed up this process allowing dose- and time-intensifications of multi-cycle chemotherapies[11, 12]. While platelets and red blood cells can show cumulative toxicity during the course of intensified regimen, appropriate G-CSF prophylaxis often results in a complete recovery of circulating granulocytes within one therapy cycle, i.e. within two or three weeks[4, 13].

Several pharmaceutical derivatives of G-CSF are available now. The first generation of G-CSF pharmaceuticals were recombinant derivatives such as Filgrastim (non-glycosylated) or Lenograstim (glycosylated). Both derivatives are virtually identical to endogenously produced G-CSF[14–16] but Filgrastim is more frequently used in clinical trials. Filgrastim is eliminated by both, renal elimination and specific degradation mediated by G-CSF receptors or neutrophil elastase[17–22]. This results in a short half-life *in vivo* requiring multiple injections during one cycle of chemotherapy. As next generation G-CSF derivative, Pegfilgrastim (pegylated Filgrastim) was developed in order to improve the pharmaceutical properties of Filgrastim. Indeed, Pegfilgrastim shows a remarkably prolonged half-life in vivo mainly due to reduced renal clearance[23–25]. Therefore, only one injection (with fixed dose) is required during one chemotherapy cycle. On the other hand, pegylation of proteins also might reduce the receptor binding affinity, and withit, the efficacy of the drug[26–29]. But for Pegfilgrastim, this effect appears to be less important than the gain in half-life, since it is generally believed that a single injection of Pegfilgrastim is at least as effective as multiple injections of Filgrastim in treating neutropenia[2, 30–33]. There are some ongoing efforts to further improve the pharmacokinetic properties of G-CSF derivatives by additional pegylations (e.g. Maxy-G34,[29]). Although the application of pegylated G-CSF is much more convenient for both patients and clinicians, we believe that Filgrastim will not completely be replaced since it can be applied more individually e.g. in dependence on the neutropenic risk of a patient or in cases when granulopoietic recovery is insufficient while pegylated G-CSF was applied[34]. Additionally, for the purpose of stem cell mobilization, Filgrastim is not inferior compared to Pegfilgrastim but has less severe side-effects[35].

In view of the highly differing pharmacokinetic properties of the available G-CSF derivatives, we constructed pharmacokinetic models of the G-CSF derivatives Filgrastim, Pegfilgrastim and the novel Maxy-G34 in mice and rats[29, 36]. Our aim was to identify basic pharmacokinetic model mechanisms especially with respect to the degradation of G-CSF *in vivo* and to compare the resulting pharmacokinetic parameters between the G-CSF derivatives. We now aim to translate these model insights to humans.

Effectiveness of G-CSF treatment depends on many variable therapy parameters such as applied chemotherapy, individual factors, G-CSF derivative used, and especially, its dosing and timing schedule[4, 34, 37, 38]. Chemotherapy induced neutropenia and G-CSF induced granulocytosis via different modes of action in combination with a strong specific elimination of G-CSF mediated by circulating granulocytes results in complex dynamics of both G-CSF serum concentration and circulating granulocytes as well. In consequence, the optimization of G-CSF treatment is a non-trivial task and cannot be performed solely on the basis of clinical trials. We showed in the past that biomathematical cell-kinetic models of granulopoiesis under chemotherapy and G-CSF support are useful to optimize chemotherapy regimen regarding granulotoxicity[39–41]. Our former model already included a preliminary model of Filgrastim application. In the present paper, we update our model with respect to improved pharmacokinetic and dynamic modelling of Filgrastim and Pegfilgrastim based on our models in mice. Additionally, our cell-kinetic model has been improved by a more elaborated model of chemotherapy action. The resulting model is now able to explain the time courses of granulocytes and G-CSF serum concentrations for virtually all datasets published in the literature. We also discuss how the model can be used to optimize G-CSF scheduling of chemotherapies.

### Structure of the human model of granulopoiesis under chemotherapy with G-CSF support

*S*(pluripotent stem cells),

*CG*(colony forming units of granulocytes and macrophages),

*PGB*(proliferating granulopoietic blasts),

*MGB*(maturing granulopoietic blasts - subdivided into metamyelocytes (

*G*4), banded granulocytes (

*G*5) and segmented granulocytes (

*G*6)) and

*GRA*(circulating granulocytes). At this, the efflux of one compartment equals the influx of the subsequent compartment. The system is highly regulated via growth factor mediated feedback loops. The most important one is G-CSF which regulates the compartments

*CG*,

*PGB*and

*MGB*, but not

*S*[5–10, 42, 43]. Modes of action comprise improvement of proliferation, acceleration of maturation and improvement of the release of mature blood cells from bone marrow to blood. The latter one is also denoted as

*postmitotic amplification*in the following. Production and consumption of G-CSF are regulated by mature cells. We also modelled a subcutaneous compartment in which G-CSF is usually injected. Chemotherapy induces an instantaneous depletion of bone marrow cell stages which is specific for each cell stage and dependent on the applied drugs and drug doses. The cell-kinetic model is essentially the same as presented and discussed in[40], except for a few changes which we will discuss later. Basic model structure is shown in Figure1. A complete set of model equations and parameters is presented in the Additional file1.

For unperturbed granulopoiesis, the model is autonomous and has a single fixed-point (steady-state). This fixed-point appears to be stable for the parameter set which we propose later, i.e. transient perturbations result in damped oscillations until the fixed-point is re-established. Permanent perturbations such as constant G-CSF stimulation or chemotherapy damage result in new fixed-points of increased or depressed granulopoiesis, respectively.

However, stable oscillations of the system can occur for alternative parameter sets especially if parameters of stem cell regulation are changed.

In general, we assume that the system is in steady-state at the beginning of any treatment. That is, initial values equals normal values.

### Basic model mechanisms

In this section, we briefly describe two important regulation mechanisms of our granulopoiesis model which are needed in the following. A more detailed discussion of all regulation mechanisms of the model is given in[40].

#### Z-function

where *X* is the regulator of the quantity *Y* and *Y*^{
min
}, *Y*^{
nor
}, *Y*^{
max
}, *b*_{
Y
}are the parameters of the Z-function for the regulation of *Y* . The parameter *b*_{
Y
}defines the steepness of the function and is called *sensitivity parameter* in the following.

#### Modelling of delays

The maturation of cells and the transition between maturing compartments is neither a random first order transition without dependence of cell age (which is equivalent to an exponential distribution of transit times) nor a transition with fixed time delay. In[40] we showed that a cascade of first order transitions results in a Gamma distribution of the transition times. More precisely, if we divide a compartment with transit time *T* in *N* subcompartments connected by first order transitions with time *T*/*N*, the resulting distribution is a Gamma distribution with expectation *T* and variance *T*^{2}/*N*. Hence, the number of subcompartments corresponds to a variance estimate regarding the transition time.

The method of dividing compartments into subcompartments in order to introduce an element of time delay will later be adopted for the modelling of delayed G-CSF or chemotherapy actions.

## Refinement of the model

### Pharmacokinetic model assumptions

- 1.
The pharmacokinetic model contains three compartments in which G-CSF is present: The subcutaneous compartment ${C}_{G-\mathrm{CSF}}^{\mathrm{sc}}$ in which G-CSF is usually injected, the central compartment ${C}_{G-\mathrm{CSF}}^{\mathrm{cent}}$ in which G-CSF is haematologically active and a peripheral compartment ${C}_{G-\mathrm{CSF}}^{\mathrm{per}}$ representing reversible bindings of G-CSF (e.g. protein binding [44]).

- 2.
Subcutaneously injected G-CSF results in a delayed influx of G-CSF into the central compartment caused e.g. by lymphatic absorption [45]. The delay is modelled by the division of the subcutaneous compartment into two subcompartments.

- 3.
Transitions between central and peripheral compartment are reversible and are modelled by a first order kinetic in both directions [46].

- 4.
Endogenous G-CSF is produced by endothelial cells [47]. The production is regulated by the demand of mature granulocytes. We implemented a phenomenological rather than mechanistic modelling of this principle. At this, the production of endogenous G-CSF is modelled as a function of the content of the final bone marrow compartment and circulating granulocytes. This is in complete analogy to former versions of our model [40].

- 5.
Since the bioavailability of G-CSF derivatives is dose-dependent [29], we assume that a part of the applied G-CSF is removed from the subcutaneous compartment without entering the central compartment. This is modelled by a Michaelis-Menten kinetic within the first subcompartment of the subcutaneous tissue.

- 6.
G-CSF is irreversibly removed from the central compartment by two independent processes: An unspecific renal elimination which is modelled by a first order kinetic [24, 46] and a specific degradation mediated by the number of circulating granulocytes. The latter one is modelled by a Michaelis-Menten kinetic which is assumed to be proportional to the number of circulating granulocytes [36, 48–51]. Two key mechanisms are discussed for the specific degradation which are cleavage by neutrophil elastase [20, 21, 52, 53] and G-CSF receptor binding and internalization [19, 22, 54–56]. Since neutrophil elastase is mainly produced by granulocytes [20], proportionality of degradation with granulocyte count can be assumed for the first mechanism. For the second mechanism, the proportionality assumption is less clear since the majority of G-CSF receptors are within the bone marrow which dynamics are somewhat different from the dynamics of mature granulocytes.

- 7.
Differences in G-CSF derivatives are model by different model parameters rather than differences in model structure. In view of the high similarity of Filgrastim and endogenous G-CSF [14–16], we assume that the pharmacokinetic and -dynamic parameters of Filgrastim and endogenous G-CSF are the same. In contrast, we assume differences between Filgrastim and Pegfilgrastim for some of the parameters. More precisely, pharmacokinetic differences were assumed with respect to absorption, distribution and degradation of the G-CSF pharmaceuticals in order to model observed differences in G-CSF serum dynamics after subcutaneous applications [23–25]. Pharmacodynamik differences are modelled by different parameterizations of the G-CSF mediated regulatory mechanisms (Z-functions). The latter is motivated by experimental results suggesting a reduced receptor binding affinity of pegylated G-CSF [26–29].

### Pharmacokinetic model equations

#### Endogenous production

where *ω*_{G 6} and *ω*_{
GRA
} are weighting parameters. It holds that${P}_{G-\mathit{\text{CSF}}}^{\mathit{\text{endo}}\text{\_}\mathit{\text{nor}}}=1$ whereas all other parameters of the Z-functions are free parameters.

#### Exogenous G-CSF application

where *t*_{
i
}≥ 0 (*i* = 1,…,*L*) are the time points of G-CSF injections and *d*_{G−CSF}(*t*_{
i
}) are the corresponding doses (in *μg*). The parameter *t*_{
inf
}is the duration of the injection which we assumed to be constant (*t*_{
inf
}= 5*s*).

The injection function is specific for each G-CSF derivative and site of injection (subcutaneous, intravenous). That is, for concurrent Filgrastim and Pegfilgrastim applications injections at concurrent sites one needs a maximum of four injection functions${P}_{G-\mathit{\text{CSF}}}^{\mathit{\text{exo}}\_\mathit{\text{sc}}\_\mathit{\text{fil}}}$,${P}_{G-\mathit{\text{CSF}}}^{\mathit{\text{exo}}\_\mathit{\text{iv}}\_\mathit{\text{fil}}}$ for subcutaneous and intravenous Filgrastim injections and${P}_{G-\mathit{\text{CSF}}}^{\mathit{\text{exo}}\_\mathit{\text{sc}}\_\mathit{\text{peg}}}$ and${P}_{G-\mathit{\text{CSF}}}^{\mathit{\text{exo}}\_\mathit{\text{iv}}\_\mathit{\text{peg}}}$ for subcutaneous and intravenous Pegfilgrastim injections respectively.

#### Subcutaneous compartment

*sc*_1 and

*sc*_2 where the efflux of the first compartment is the influx to the second compartment. G-CSF is applied to the first subcompartment (second term of (4)). In the first subcompartment there is a dose-dependent loss of G-CSF modelled by a Michaelis-Menten kinetic (third term of (4)). For Filgrastim injections it holds that

with the initial values${C}_{G-\mathit{\text{CSF}}}^{\mathit{\text{sc}}\_1}\left(0\right)\phantom{\rule{2.77626pt}{0ex}}\phantom{\rule{2.22144pt}{0ex}}=\phantom{\rule{2.77626pt}{0ex}}\phantom{\rule{2.22144pt}{0ex}}{C}_{G-\mathit{\text{CSF}}}^{\mathit{\text{sc}}\_2}\left(0\right)\phantom{\rule{2.77626pt}{0ex}}\phantom{\rule{2.22144pt}{0ex}}=\phantom{\rule{2.77626pt}{0ex}}\phantom{\rule{2.22144pt}{0ex}}0$. For Pegfilgrastim injections the first term of the right-hand side of (4) is substituted by${P}_{G-\mathit{\text{CSF}}}^{\mathit{\text{exo}}\_\mathit{\text{peg}}\_\mathit{\text{sc}}}$. Likewise, the Filgrastim parameters${k}_{\mathit{\text{sc}}}^{F}$,${v}_{\mathit{\text{max}}}^{F}$ and${k}_{m}^{F}$ are substituted by corresponding Pegfilgrastim parameters.

#### Central compartment

The balance equation (6) contains terms in the following order: The endogenous production, a potential intravenous injection, the influx from the subcutaneous compartment, the unspecific elimination, the efflux to the peripheral compartment, the influx from the peripheral compartment and the specific elimination (Michaelis-Menten kinetic) which is proportional to the relative granulocyte concentration${C}_{\mathit{\text{GRA}}}^{\mathit{\text{rel}}}$. The corresponding equation for Pegfilgrastim is the same except for the endogenous production, which is zero, the intravenous injection function which is substituted by${P}_{G-\mathit{\text{CSF}}}^{\mathit{\text{exo}}\_\mathit{\text{peg}}\_\mathit{\text{iv}}}$, the parameters and the initial value which is again zero.

#### Peripheral compartment

where the parameters *k*_{
cp
}, *k*_{
pc
} and *V*_{
D
} are specific for Filgrastim and Pegfilgrastim respectively.

### New pharmacodynamic model assumptions

**Major compartments of the pharmacokinetic model and corresponding regulations**

Compartment | Regulations | Regulator |
---|---|---|

| proliferative fraction | bone marrow content |

self-renewal probability | bone marrow content | |

| proliferative fraction | bone marrow content |

amplification | G-CSF | |

transition time | G-CSF | |

| amplification | G-CSF |

transition time | G-CSF | |

| post-mitotic amplification | G-CSF |

transition time | G-CSF | |

| turn-over | - |

| endogenous production | late bone marrow cell stages |

specific degradation |
| |

external applications | - |

We used the same pharmacodynamic model of G-CSF as presented in[40] except for a few simplifications and additional assumptions which we will discuss now:

The cytokine GM-CSF (granulocyte macrophage colony-stimulating factor) was no longer considered in order to simplify the model. Since the endogenous productions of GM-CSF and G-CSF in our former model were both related to the demand of bone marrow cell stages, the time courses of endogenous G-CSF and GM-CSF after external perturbations are similar, making the cytokines undistinguishable from the modelling point of view (compare[40]). Therefore, only G-CSF was considered in the present model. It replaces GM-CSF regarding the regulation of the *CG* compartment. This simplification is also biologically plausible since G-CSF receptors are also expressed at myeloid progenitors[42].

The transition times *T*_{
CG
}and *T*_{
PGB
}were constant in our former version of the model but are now dependent on G-CSF via a Z-function. This model assumption was made in order to account for the increased number of mitoses in these compartments. But, due to the lack of data, it must be considered as speculative.

Since Pegfilgrastim and Filgrastim are supposed to have different G-CSF receptor binding affinities[26–29], we assume different regulatory Z-functions for both derivatives. But we assumed the same pharmacodynamic parameters for Filgrastim and endogenously produced G-CSF.

*D*

_{G−CSF}is specific for Pegfilgrastim and Filgrastim but constant for all Z-functions. Hence, (1) reads

for all quantities *Y* regulated by G-CSF. At this, the normal value of the Z-function *Y*^{
nor
}refers to${C}_{G-\mathit{\text{CSF}}}^{\mathit{\text{nor}}}$ for endogenous G-CSF or Filgrastim and to an absolute amount of 1*μ* g for Pegfilgrastim.

*ω*

_{ P }which is again a Z-function of the quotient of the two doses with minimum 0 and maximum 1.

where *Y* is an arbitrary regulated quantity such as a transition time or an amplification. For all these quantities we assumed the same Z-function of the weighting parameter *ω*_{
P
}.

In the previous version of our model, chemotherapy was modelled by an instantaneous depletion of bone marrow cells lasting exactly one day. Since independence of cytotoxic action of the single chemotherapy components was assumed, the effect of chemotherapy could be modelled by a step function with a step width of one day. Although the metabolism of cytotoxic drugs is usually fast, the nadir of bone marrow cell stages is typically reached a few days after the application[58]. To account for this observation, we delayed the toxic effect of chemotherapy applications by a cascade of four subcompartments in complete analogy to (10), (11) resulting in a delayed toxicity function *Ψ*^{
GRA
}. The delay parameter is specific for the cytotoxic drugs used. Hence, two parameters are required to define the toxic effect of a chemotherapeutic drug to a single cell line: a delay parameter and a scaling factor of the toxicity function. While the delay is only specific for the applied chemotherapeutic drug but constant for all cell stages, the scaling factor is specific for both.

where *c*_{
LY
}= 3000 cells per *μl* and *c*_{
GRA
}= 4000 cells per *μl* are the normal concentrations of lymphocytes and granulocytes respectively. *Ψ*_{
LY
}is the toxicity function for lymphocytes which is analogously defined as the toxicity functions of granulopoiesis (see 6.)

### Construction of toxicity functions

*i*= 2,…,4,

*X*∈{

*LY*,

*GRA*}, i.e. the delay parameter

*D*is specific for each chemotherapy and different for the toxicity functions of granulocytes and lymphocytes. Now, the toxicity functions are defined as:

where *k* are the toxicity parameters and *Y* *ε*{*S*,*CG*,*PGB*,*MGB*}, i.e. the toxicity function is specific for the different cell stages of granulopoiesis. If drugs are applied at different schedules, corresponding toxicity functions are added. The sketched principle can easily be generalized to derive toxicity functions of arbitrary chemotherapy schedules.

## Model calibration, parameter estimation and validation

### Estimation of parameters

A main goal of our study is to construct pharmacokinetic models of Filgrastim and Pegfilgrastim. Based on the pharmacokinetic model equations presented above, a set of unknown pharmacokinetic parameters needs to be determined.

Since detailed bone marrow data of human granulopoiesis are not available, most of the bone marrow parameters are only known up to a certain range or are completely unknown. Furthermore, model parameters regarding sensitivity of regulatory mechanisms (sensitivity parameters) have no direct biological measurable equivalent. Finally, we want to apply the model to chemotherapy settings requiring a quantification of corresponding toxicity and delay parameters. Many of the parameters of the present model version were adopted from an earlier version of the model especially if they are not very sensitive regarding model behaviour (compare[40]). But the inclusion of new regulatory mechanisms (see section “New pharmacodynamic model assumptions”) made some adaptations of model parameters necessary.

- 1.
Pharmacokinetic model parameters were determined on the basis of available cytokine dynamics after G-CSF application. In order to model the specific elimination mechanism, we imprinted the corresponding data of granulocyte dynamics at this stage of modelling. We obtained a unique parameter set which is valid for all dosing and timing schedules of all scenarios considered.

- 2.
Pharmacodynamic parameters were determined by fitting the predictions of the model to available granulocyte and leukocyte dynamics of different scenarios (G-CSF application and simple chemotherapies for which the number of chemotherapy parameters to be fitted is relatively low). The resulting parameter set is valid for all scenarios with and without chemotherapy applications. Since the stem cell compartment is the basis of all of our models of haematopoietic lineages, we decided to keep corresponding parameters constant and as presented in [40, 59].

- 3.
Afterwards, toxicity parameters of more complex chemotherapies can be estimated. More details of chemotherapy modelling, parameter estimation and exploration of patients with different risk of haematotoxicity can be found in a separate publication of our group (Wetzler et al., to appear).

A complete parameter list for our model is provided in the Additional file1. Not all data sets were used to fit parameters. A few data sets were kept in reserve in order to validate the model.

### Available data sets

Data sets were collected from the literature by an extensive search. For our modelling purposes, close-meshed time series of G-CSF and absolute neutrophil counts (ANC) or leukocytes are especially valuable. Corresponding data were extracted from the publications as precise as possible using automated tools. Data sets for which no means or medians of the patients could be retrieved were neglected. Data sets comprise single or multiple applications of G-CSF in healthy volunteers and conventional chemotherapies of different diseases with or without G-CSF prophylaxis.

*German High Grade Non-Hodgkin-Lymphoma Study Group*. These studies were conducted in accordance with the Declaration of Helsinki. Corresponding protocols were approved by the ethics review committee of each participating center. Written informed consent was obtained from the patient for publication of this report and any accompanying images. An overview of data used for model fitting and validation is given in Table2.

**Data sets utilized for the establishment and validation of the granulopoiesis model**

Type of data | Disease | G-CSF schedules | Chemotherapy | References |
---|---|---|---|---|

Phase I studies with Filgrastim | none | single application 3, 5, 10 | none | [31] |

none | single application 5, 10 | none | [60] | |

none | 10 applications 75, 150, 300, 600 | none | [61] | |

none | 14 applications 30, 300 | none | [62] | |

none | single application 4, 8 | none | [63] | |

Phase I studies with Pegfilgrastim | none | 30, 60, 100, 300 | none | [50] |

Phase II studies with Pegfilgrastim | LuCa | 30, 100, 300 | none | [64] |

NHL | 6000 | CHOP | [65] | |

Phase III studies with CX and w/wo Filgrastim | NHL | no G-CSF | CHOP-21* | [12] |

NHL | 480 | CHOP-14* | [12] | |

NHL | 480 | CHOP-14* | [66] | |

Phase III studies with CX + Peg | NHL | 6000 | CHOP-14* | [67] |

### Fitting procedure

*L*

^{1}distance between logarithmized model and logarithmized median of data. More precisely, we have

where *f*_{
model
}(*t*,**k**) is the solution of the model equation system for the granulocyte compartment at the time *t* based on the parameter set **k** = *k*_{1},…*k*_{
n
}. For each scenario, *t*_{0} and *t*_{1} describe the first and the last time point for which data are available. To obtain the curve log*f*_{
data
}(*t*), the logarithms of the patients medians were linearly interpolated. Logarithms of data were used to provide an optimal fit of the nadir phase of cell counts. In the following, the left hand side of equation (16) is referred to as the *fitness function*.

As in our previous papers, evolutionary strategies were used for the numerical solution of the optimisation problem. This method is especially suitable for our problem since it requires a minimum of computationally expensive calculations of the fitness function, it can deal with a large number of free parameters with only a linear growth in effort and it is the only chance to obtain a global optimum as good as possible.

Evolutionary strategies are non-deterministic optimization methods which are based on the principles of evolution (mutation by chance, reproduction, realization of phenotypes and survival of the fittest). For mathematical optimisation, parameter settings were taken instead of livings, that is parental parameter settings are changed by chance (mutation), combined to form new parameter settings (reproduction) and were used to solve the model equation system (realization). The parameter settings for which one obtains a good agreement between the model prediction and the data were taken to create the next generation of parameter settings (survival of the fittest). The fitness function is a measure for this agreement. We used a (1+3) evolutionary strategy with self-adapting mutation step size most of the time. That is, one possibly immortal parent creates three children in each step. See also[68, 69] for further details of evolutionary strategies.

Fitting of chemotherapy schedules requires additional parameters with respect to the chemotherapy toxicity. The delay parameter of the toxicity is specific for each drug but is the same for all bone marrow cell stages. Four toxicity parameters are required to model the cell stage specific toxicities. Another parameter represents increased toxicity for the first chemotherapy application (first cycle effect). Finally, two parameters represent the toxicity to the lymphopoetic system (eq. (15)). In general, modelling a chemotherapy required several sets of these parameters in order to model all drugs or drug combination with different schedules. However, for the purpose of model calibration, we only considered simple chemotherapy schedules in order to reduce the number of additional unknown parameters to be fitted. The CHOP regimens are based on the application of three cytotoxic drugs (cyclophosphamide, doxorubicine and vincristine) at the same time. Hence, only one set of chemotherapy parameters was assumed to model the effects of this drug combination. Additionally, for these regimens different G-CSF schedules are available, which is especially useful for the calibration of our pharmacokinetic and -dynamic models of Filgrastim and Pegfilgrastim.

Not all of our data sets were utilized for parameter fitting. A subset of data sets was kept in reserve in order to validate the model. Data sets used for model validation comprise the data of[63, 65].

### Sensitivity analysis

Since our model contains several parameters which are speculative or unknown only up to a certain range, we peformed an extensive sensitivity analysis of all parameters. For this purpose, parameters were increased or decreased by 2.5% and the corresponding change of the fitness function was determined. At this, only affected model scenarios were considered. The changes of the fitness function were plotted as bar diagrams for each parameter in order to facilitate comparisons. Figures are shown in the Additional file1.

### Simulation and numerical methods

Our model has been programmed with MATLAB 7.5.0.342 (R2007b) and SIMULINK toolbox (The MathWorks Inc., Natick, MA, USA). Simulations of the model were performed by numerical integration of the equation system using the variable step solver from Adams and Bashford implemented in the SIMULINK toolbox.

## Results

### Pharmacokinetic model of Filgrastim and Pegfilgrastim

**Pharmacokinetic parameters of Filgrastim and Pegfilgrastim**

Parameter | Meaning | Filgrastim | Pegfilgrastim |
---|---|---|---|

| subcutaneous absorption [ | 0.161 | 0.107 |

| Michaelis-Menten constant of | 34.7 | 5.5 |

subcutaneous elimination [ | |||

| Maximum of subcutaneous | 67.3 | 16.5 |

elimination [ | |||

| unspecific elimination [ | 0.441 | 0.087 |

${k}_{m}^{\mathit{\text{GRA}}}$ | Michaelis-Menten constant of | 22.4 | 30.8 |

specific elimination [ | |||

${v}_{\mathit{\text{max}}}^{\mathit{\text{GRA}}}$ | Maximum of specific elimination [ | 4.77 | 5.16 |

| transition central to peripheral [ | 0.000 | 0.075 |

| transition peripheral to central [ | - | 0.548 |

| distribution volume [ | 1.156 | 4.091 |

${C}_{G-\mathit{\text{CSF}}}^{\mathit{\text{cent}}\_\mathit{\text{ref}}}$ | reference G-CSF serum concentration$\left[\frac{\mu g}{l}\right]$ | 0.02 |

### The new pharmacodynamic model

Unknown pharmacodynamic parameters were determined by fitting the predictions of our model to available time courses of ANC or leukocytes after single or multiple injections of G-CSF or chemotherapy. All scenarios presented in section “Available data sets” were used for this simultaneous fitting process except for those reserved for model validation (compare section “Fitting procedure”). Fitted parameters resulted in a good fit of all model scenarios except for time points shortly after G-CSF injections. Examples are shown in Figure4. All other scenarios of G-CSF application can be found in the Additional file1. Chemotherapy scenarios are presented in the next section.

The most sensitive pharmacodynamic parameters influencing the behaviour of the model are parameters regarding the CG compartment (regulation of the proliferative fraction, transition time and amplification), the amplification in the PGB compartment and the postmitotic amplification. Sensitivity parameters of the regulation functions are generally less sensitive except for the sensitivity parameter of the regulation of the postmitotic amplification. A complete list of the results of the sensitivity analysis can be found in the Additional file1.

The estimated delay is moderate in size resulting in a small shift of the time course of cell stages. This shift is negligible for *MGB* but more pronounced for the granulocyte compartment which can be explained by the postmitotic amplification mechanism (see[40]).

### The new chemotherapy model

The delay resulted in a later occurrence and a reduced depth of the nadir of leukocytes. The delay was assumed to be different for different chemotherapeutic drugs or drug combinations (see Wetzler *et al.* for further details).

Comparison of model and data revealed a good agreement. For almost all time points, the model curve is within the interquartile range of the data. Chemotherapy and delay parameters for the CHOP regimen can be found in the Additional file1.

### Validation of the model

*et al.*[63] were not used for model fitting as well as the data of CHOP chemotherapy under Pegfilgrastim treatment of George

*et al.*[65]. Compared to the CHOP data used for model fitting, the data of George

*et al.*[65] comprise G-CSF serum levels as well. Both scenarios fit well with our model prediction (see Figure10). No additional parameter fittings were performed to model these scenarios.

## Discussion

In the present paper, we developed an ordinary differential equations model of human granulopoiesis under chemotherapy and G-CSF support. The model was built on the basis of a former model of granulopoiesis of our group which now has been improved primarily by the incorporation of a detailed pharmacokinetic and -dynamic model of two G-CSF derivatives (Filgrastim and Pegfilgrastim). At this, the pharmacokinetic model was adopted from similar models developed for G-CSF applications in mice and rats developed by our group. Unknown model parameters were obtained by fitting the predictions of the model to available datasets. The combined pharmacokinetic and -dynamic model correctly predicts the time course of a variety of datasets comprising single or multiple injections of G-CSF into healthy volunteers or patients under CHOP chemotherapy. We were able to describe the differences between the G-CSF derivatives by a set of different pharmacokinetic- and -dynamic parameters. The model was validated on the basis of datasets not used for model fitting.

The presented model is by far not the first attempt to model granulopoiesis or G-CSF applications. Published models comprise for example pure pharmacokinetic models[25, 46], pharmacokinetic and -dynamic models of G-CSF application on the cellular level[49], in healthy volunteers[50, 51, 70], for the treatment of cyclic neutropenia[71], for high-dose chemotherapy with stem cell transplantation[72, 73] or for conventional chemotherapy patients[74, 75]. We developed a model of human granulopoiesis under chemotherapy in the past including a preliminary model of Filgrastim application[40]. To the best of our knowledge, there is so far no granulopoiesis model of humans under conventional chemotherapy comprising a detailed pharmakokinetic- and dynamic model of the two G-CSF derivatives Filgrastim and Pegfilgrastim. This combination allows us to derive clinically meaningful applications of the model.

As mentioned, the presented model was based on a former model of our group. This model was based on biologically plausible assumptions regarding the production of mature granulocytes via a cascade of bone marrow cell stages, the action of chemotherapy and the action of growth factor mediated feedback loops. Equations describe the fluxes between cell compartments. G-CSF was modelled as the major regulatory element of both, the transition time and the amplification within the compartments. Chemotherapy was modelled by an instantaneous and transient cell loss of all bone marrow cell stages. Since these basic model assumptions were intensively discussed in our former paper[40] we will focus on our new model assumptions and parameters in the following.

The major improvement of our model is the incorporation of a detailed pharmacokinetic model of the two G-CSF derivatives Filgrastim and Pegfilgrastim. Both are widely used in clinical practice in order to ameliorate leukopoenia during cancer chemotherapy. The pharmacokinetic model was constructed in complete analogy to the pharmacokinetic models which we developed for mice and rats recently[29, 36]. That is, we made the same physiological assumptions and used the same model equations but different parameters. Furthermore, we used the same model equations for both drugs assuming that pharmacokinetic differences between the drugs can be traced back to different parameters rather than different mechanisms of action.

The drugs were typically injected into the subcutaneous tissue resulting in a delayed absorption by the circulating blood compartment probably via lymphatic absorption[45]. The delay was modelled by a set of concatenated first-order differential equations rather than a fixed time delay. We showed in the past that this kind of modelling is equivalent to a Gamma-distributed transition time, which is biologically plausible. At this, the variance is determined by the number of subcompartments[40]. However, this variance appears to be of lesser importance for the model behaviour. In analogy to[36], we observed that a number of subcompartments between two and ten would also work well. To reduce the computational in the present model, we used the smaller number.

Data collected in mouse and rat experiments suggested that subcutaneously injected G-CSF has a dose-dependent bioavailability[29, 36]. Therefore, we introduced a loss term into the equations of the subcutaneous compartment. Since there is some evidence of reversible protein binding of G-CSF molecules, we modelled a first order transition between the blood compartment and a peripheral compartment[46].

Our model assumptions regarding endogenous production of G-CSF are speculative. We assumed that the production is regulated between a minimum and a maximum value in dependence on bone marrow cellularity[47]. In steady-state, the production is constant in order to sustain a fixed serum concentration determined by averaged data from the literature.

Degradation of G-CSF was modelled by two independent processes, an unspecific renal clearance modelled by a first order transition and a specific degradation via neutrophil elastase or G-CSF receptors. All three degradation mechanisms are biologically well understood but their relative importance is unknown (see discussed literature in section “Pharmacokinetic model assumptions”). The specific degradation was modelled by a Michaelis-Menten kinetic which was assumed to be proportional to the number of granulocytes. For the degradation mechanism due to neutrophil elastase, this assumption seems to be appropriate[21]. On the other hand, since G-CSF receptors are also present in bone marrow progenitors and precursors[42], it appears to be less appropriate for the receptor-mediated clearance mechanism. Nevertheless, we assumed proportionality as the most parsimonious model resulting in a good agreement of predictions and data. The assumption also worked well for the pharmacokinetic modelling in mice[36] but not for rats[29]. Moreover, we experimented with alternative model assumptions assuming consumption of G-CSF by both, mature neutrophils and precursors which did not significantly improve the quality of our model predictions.

Variable numbers of G-CSF receptors per cell were observed in dependence on the G-CSF level[56]. Modelling this obervation would require additional assumptions and parameters. We decided to skip this in the current version of the model in view of the relatively good quality of model predictions in the clinical scenarios considered so far.

The model equations worked well to explain the dynamics of G-CSF serum concentrations after single or multiple injections of Filgrastim or Pegfilgrastim in healthy or diseased people. The same equations worked also well for a third G-CSF derivative, namely Maxy-G34, which is a novel G-CSF derivative currently under development by Maxygen Inc. However, due to a confidentiality agreement with Maxygen Inc., the results are not shown in the present paper.

In order to make predictions regarding the response of granulopoiesis after the application of G-CSF, it was necessary to attach a pharmacodynamic model of G-CSF applications. Since there is some evidence that the pegylations of the Pegfilgrastim molecule interact with its binding affinity to the G-CSF receptor[26–29], we decided to assume different regulation functions for Filgrastim and Pegfilgrastim. Hence, Z-functions of the transition times and amplifications in *CG*, *PGB* and *MGB* are assumed to be different for Filgrastim and Pegfilgrastim. On the other hand, we assumed the same Z-functions for Filgrastim and endogenous G-CSF. Due to the fact that Filgrastim/endogenous G-CSF and Pegfilgrastim were assumed to be different regarding Z-functions, it was necessary to merge the superimposing effects of Pegfilgrastim and endogenous G-CSF or concurrent Pegfilgrastim and Filgrastim applications as well. This was solved by a weighting factor which is regulated between zero and one in dependence on the ratio of Pegfilgrastim and Filgrastim or endogenous G-CSF in the system. If the Pegfilgrastim concentration is high or low, then the system is mainly influenced by the Z-function of Pegfilgrastim or Filgrastim respectively. Although this assumption is plausible, the complete regulation mechanism via the combined Z-functions must be considered as speculative, since especially the shape of the regulation functions can hardly be observed or measured.

Another speculative mechanism introduced into our model update was the delayed effect of G-CSF action. By model fitting, we estimated that the corresponding delay time is about 6h which appears to be in the right order of magnitude compared to the dynamics presented in[57]. However, the overall impact of the delay on model dynamics is limited. At least for the scenarios considered in the present paper, it is not critical for the quality of the agreement of model and data.

Furthermore, some adjustments were performed regarding chemotherapy modelling. Instead of an instantaneous cytotoxic effect of chemotherapy assumed in our former model version, we now assume a delayed effect to account for available data of the dynamics of bone marrow cell stages after chemotherapy applications in mice[58]. The delay parameter was assumed to be constant for all cell stages but specific for the applied drugs or drug combinations and was modelled by a cascade of first order transitions. This modelling is rather a phenomenological than a mechanistic approach since the delay is caused by many factors such as toxification of the applied drugs at different time scales, transient cell cycle arrests of cells and delayed apoptosis of cells due to irreversible damage[58].

Another improvement of our model is due to a semi-explicite modelling of lymphocyte toxicity. This was necessary in order to apply the model to a sufficiently large dataset of time courses of both, granulocytes and leukocytes as well. In our former model, we assumed proportionality of granulocytes and leukocytes during G-CSF application or chemotherapy which is only roughly correct[38] and further unpublished data of our chemotherapy studies). To avoid a complete cell kinetic model of lymphopoiesis requiring a large set of new and unknown model parameters, lymphocyte counts were modelled by a separate simple characteristic. We assumed no effect of G-CSF on lymphopoiesis but a toxic effect of chemotherapy modelled by an exponential depletion of lymphocytes according to a (delayed) chemotherapy toxicity function. This toxicity is again specific for the drugs or drug combinations applied. The resulting model was able to explain the time courses of leukocytes and granulocytes under G-CSF or chemotherapy adequately within the framework of one model.

The model is based on a relatively large set of parameters. Due to missing bone marrow data during chemotherapy and growth factor application, only a limited knowledge regarding the required cell-kinetic and toxicity parameters is available. Rough ranges for transition times and amplification rates in steady-state or under stimulation by G-CSF were obtained from the literature. But especially values under minimal stimulation, the sensitivity parameters of the Z-functions and the toxicity parameters are not available from literature data. Hence, many model parameters were determined indirectly by fitting the predictions of the model to available datasets. For this purpose, we collected a set of suitable data from the literature and clinical trials for which we have access to raw data. Densely measured time courses of G-CSF serum levels after application in combination with granulocyte or leukocyte counts after chemotherapy and different G-CSF schedules are especially useful. Data of patients were pooled to construct a model that fits to the median of patients. A unique parameter set was identified which is valid for all scenarios considered. No adjustments were performed in order to fit single scenarios. Not all datasets were used for model fitting enabling an opportunity for model validation.

Despite of the utilization of several datasets comprising different G-CSF dosing and timing schedules with and without chemotherapy, there remained a large uncertainty regarding parameter estimates, and consequently, the current parameter settings must be considered as preliminary. This is especially true for parameters with a low impact on our fitness function in the scenarios considered as demonstrated by our sensitivity analysis. Additionally, the toxicity parameters show some degree of dependence in the sense that a higher toxicity at one cell stage can to some degree be compensated by a lower toxicity at a subsequent cell stage and vice versa. Consequently, further datasets and validation scenarios are required to improve the confidence regarding our parameter settings.

The estimates of our pharmacokinetic parameters resulted in a good fit of all time series data of G-CSF serum concentration for both Filgrastim and Pegfilgrastim applications as well. The estimated values fit well to our biological understanding of the drugs. Due to the pegylation of the drug, it was expected that the unspecific renal clearance is significantly reduced for Pegfilgrastim, which is in agreement with our parameter estimates. We also estimated a reduced specific degradation of Pegfilgrastim which could be explained by a reduced receptor binding affinity or hydrophilic properties of pegylated molecules[26, 27, 76]. We also made the same observation for our pharmacokinetic models constructed in mice and rats[29, 36]. Protein binding was estimated to be almost negligible for Filgrastim but important for Pegfilgrastim in agreement with our observations in mice[36]. Finally, the estimated distribution volumes are in rough agreement with findings of other authors[50, 77].

The estimates of our cell-kinetic, pharmacodynamic and toxicity parameters also resulted in a good fit of the time courses of granulocytes and leukocytes after application of chemotherapy, G-CSF or combinations of it. Possible exceptions are cell counts measured shortly after the first application of Filgrastim, which seem to be underestimated in some scenarios. We conclude that the model works well on a day-wise scale but might be unable to explain short-term or transient effects of G-CSF applications e.g. on the scale of hours. Modelling such short-term effects would require a better database, since almost all available time courses of granulocytes and leukocytes were measured at most at a day-wise scale. However, in order to make predictions regarding the efficiency of different G-CSF schedules, we are also more interested in the long-term dynamics of granulocytes and leukocytes in the course of the therapy rather than short-term effects after single injections.

Estimates of the pharmacodynamic parameters of Filgrastim and Pegfilgrastim suggest that Filgrastim has a higher potency to stimulate the bone marrow. This is in agreement with our biological understanding that pegylation reduces the receptor binding affinity. An analogous observation was made for Pegfilgrastim and the novel drug Maxy-G34 which has even more pegylation sites than Pegfilgrastim[29].

Although modelling of chemotherapy was not the primary goal of the present paper, it was necessary to model at least a few conventional chemotherapy regimen to study the pharmacokinetic properties of the G-CSF derivatives under granulopenic conditions. Chemotherapy was modelled as a transient delayed toxic effect on all bone marrow cell stages. Corresponding toxicity parameters are specific for the bone marrow cell stages and the drugs and drug concentrations used. For the development of our model, we used the data of the most simple chemotherapy regimen CHOP. With our toxicity parameter estimates one obtains a good fit of all CHOP regimen with different G-CSF schedules.

Other conventional chemotherapies were modelled by assuming different toxicity parameters but the same cell kinetic model. Corresponding model simulations will be demonstrated in a separate paper of our group comprising about 20 different chemotherapy scenarios. Generally, the model can be applied to arbitrary conventional chemotherapy regimens for which data of leukocyte or granulocyte time courses are available for at least one G-CSF scheduling. Based on these data, sets of toxicity parameters of corresponding chemotherapies can be estimated. Using these parameters, it is possible to make clinically relevant predictions regarding the time course of G-CSF serum concentrations, bone marrow cell stages and mature cell counts in circulation under different G-CSF schedules, allowing to optimize G-CSF treatment. We will exploit the clinically relevant applications of our model in the near future.

## Conclusions

We established a human pharmacokinetic and -dynamic model of Filgrastim and Pegfilgrastim applications under cytotoxic chemotherapy. The model is able to explain a large number of clinical time series data of G-CSF serum concentrations, granulocytes and leukocytes of patients treated with G-CSF and with or without chemotherapy. A unique parameter set valid for all scenarios was established by fitting the predictions of the model to clinical data. The model was validated on a set of scenarios not used for parameter fitting. Differences between Filgrastim and Pegfilgrastim could be traced back to biologically plausible differences in parameter estimates. Effects of chemotherapy can be quantified by a set of toxicity parameters. Given these toxicity parameters, the model can be used to simulate the dynamics of G-CSF, bone marrow cell stages and circulating granulocytes or leukocytes of yet untested G-CSF schedules. The model is currently applied in the planning phase of clinical trials in order to optimize G-CSF treatment.

## Declarations

### Acknowledgements

S.S. and M.W. were funded by a grant of the Federal Ministry of Education and Research of the Federal Republic of Germany (“Haematosys”, BMBF / PTJ0315452A). M.S. was funded by LIFE – Leipzig Research Center for Civilization Diseases, University of Leipzig. LIFE is funded by means of the European Union, by the European Regional Development Fund (ERDF) and by means of the Free State of Saxony within the framework of the excellence initiative.

## Authors’ Affiliations

## References

- Crawford J: Pegfilgrastim administered once per cycle reduces incidence of chemotherapy-induced neutropenia. Drugs. 2002, 62 (Suppl 1): 89-98.View ArticlePubMedGoogle Scholar
- Dale D: Current management of chemotherapy-induced neutropenia: the role of colony-stimulating factors. Semin Oncol. 2003, 30: 3-9.View ArticlePubMedGoogle Scholar
- 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: S39-s47.View ArticlePubMedGoogle Scholar
- Wunderlich A, Kloess M, Reiser M, Rudolph C, Truemper L, Bittner S, Schmalenberg H, Schmits R, Pfreundschuh M, Loeffler M: German High-Grade Non-Hodgkin’s Lymphoma Study Group (DSHNHL): 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: 881-893.View ArticlePubMedGoogle Scholar
- Souza LM, Boone TC, Gabrilove J, Lai PH, Zsebo KM, Murdock DC, Chazin VR, Bruszewski J, Lu H, Chen KK, Barendt J, Platzer E, Moore MAS, Mertelsmann R, Welte K: Recombinant human granulocyte colony-stimulating factor: effects on normal and leukemic myeloid cells. Science. 1986, 232: 61-65.View ArticlePubMedGoogle Scholar
- Begley CG, Nicola NA, Metcalf D: Proliferation of normal human promyelocytes and myelocytes after a single pulse stimulation by purified GM-CSF or G-CSF. Blood. 1988, 71: 640-645.PubMedGoogle Scholar
- 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 USA. 1989, 86: 9499-9503.PubMed CentralView ArticlePubMedGoogle Scholar
- Mackey MC, Aprikyan AA, Dale DC: The rate of apoptosis in post mitotic neutrophil precursors of normal and neutropenic humans. Cell Prolif. 2003, 36: 27-34.View ArticlePubMedGoogle Scholar
- Kim HK, De La Luz Sierra M, Williams CK, Gulino AV, Tosato G: G-CSF down-regulation of CXCR4 expression identified as a mechanism for mobilization of myeloid cells. Blood. 2006, 108: 812-820.PubMed CentralView ArticlePubMedGoogle Scholar
- Christopher MJ, Link DC: Regulation of neutrophil homeostasis. Curr Opin Hematol. 2007, 14: 3-8.View ArticlePubMedGoogle Scholar
- Diehl V, Franklin J, Pfreundschuh M, Lathan B, Paulus U, Hasenclever D, Tesch H, Herrmann R, Dorken B, Muller-Hermelink HK, Duhmke E, Loeffler M: Standard and increased-dose BEACOPP chemotherapy compared with COPP-ABVD for advanced Hodgkin’s disease. N Engl J Med. 2003, 348: 2386-2395.View ArticlePubMedGoogle Scholar
- Pfreundschuh M, Truemper L, Kloess M, Schmits R, Feller AC, Ruebe C, Rudolph C, Reiser M, Hossfeld DK, Eimermacher H, Hasenclever D, Schmitz N, Loeffler M, German High-Grade Non-Hodgkin’s Lymphoma Study Group: 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.View ArticlePubMedGoogle Scholar
- Engel C, Loeffler M, Schmitz S, Tesch H, Diehl V: Acute hematologic toxicity and practicability of dose-intensified BEACOPP chemotherapy for advanced stage Hodgkin’s disease. German Hodgkin’s Lymphoma Study Group (GHSG). Ann Oncol. 2000, 11 (9): 1105-14.View ArticlePubMedGoogle Scholar
- Frampton JE, Lee CR, Faulds D: Filgrastim. A review of its pharmacological properties and therapeutic efficacy in neutropenia. Drugs. 1994, 48 (5): 731-60.View ArticlePubMedGoogle Scholar
- Frampton JE, Yarker YE, Goa KL: Lenograstim. A review of its pharmacological properties and therapeutic efficacy in neutropenia and related clinical settings. Drugs. 1995, 49 (5): 767-93.View ArticlePubMedGoogle Scholar
- Houston AC, Stevens LA, Cour V: Pharmacokinetics of glycosylated recombinant human granulocyte colony-stimulating factor (lenograstim) in healthy male volunteers. Br J Clin Pharmacol. 1999, 47: 279-284.PubMed CentralView ArticlePubMedGoogle Scholar
- Tanaka H, Tokiwa T: Influence of renal and hepatic failure on the pharmacokinetics of recombinant human granulocyte colonystimulating factor (KRN8601) in the rat. Cancer Res. 1990, 50: 6615-6619.PubMedGoogle Scholar
- Khwaja A, Carver J, Jones HM, Paterson D, Linch DC: Expression and dynamic modulation of the human granulocyte colony-stimulating factor receptor in immature and differentiated myeloid cells. Br J Haematol. 1993, 85: 254-259.View ArticlePubMedGoogle Scholar
- Ericson SG, Gao H, Gericke GH, Lewis LD: The role of polymorphonuclear neutrophils (PMNs) in clearance of granulocyte colony-stimulating factor (G-CSF) in vivo and in vitro. Exp Hematol. 1997, 25: 1313-1325.PubMedGoogle Scholar
- El Ouriaghli F, Fujiwara H, Melenhorst JJ, Sconocchia G, Hensel N, Barrett AJ: Neutrophil elastase enzymatically antagonizes the in vitro action of G-CSF: implications for the regulation of granulopoiesis. Blood. 2003, 101: 1752-1758.View ArticlePubMedGoogle Scholar
- Hunter MG, Druhan LJ, Massullo PR, Avalos BR: Proteolytic cleavage of granulocyte colony-stimulating factor and its receptor by neutrophil elastase induces growth inhibition and decreased cell surface expression of the granulocyte colony-stimulating factor receptor. Am J Hematol. 2003, 74: 149-155.View ArticlePubMedGoogle Scholar
- Kotto-Kome AC, Fox SE, Lu W, Yang BB, Christensen RD, Calhoun DA: Evidence that the granulocyte colony-stimulating factor (G-CSF) receptor plays a role in the pharmacokinetics of G-CSF and PegG-CSF using a G-CSF-R KO model. Pharmacol Res. 2004, 50: 55-58.View ArticlePubMedGoogle Scholar
- Tanaka H, Satake-Ishikawa R, Ishikawa M, Matsuki S, Asano K: Pharmacokinetics of recombinant human granulocyte colonystimulating factor conjugated to polyethylene glycol in rats. Cancer Res. 1991, 51: 3710-3714.PubMedGoogle Scholar
- Crawford J: Clinical uses of pegylated pharmaceuticals in oncology. Cancer Treat Rev. 2002, 28 (Suppl A): 7-11.View ArticlePubMedGoogle Scholar
- Yang BB, Lum PK, Hayashi MM, Roskos LK: Polyethylene glycol modification of filgrastim results in decreased renal clearance of the protein in rats. J Pharm Sci. 2004, 93: 1367-1373.View ArticlePubMedGoogle Scholar
- Harris JM, Chess RB: Effect of pegylation on pharmaceuticals. Nat Rev Drug Discov. 2003, 2: 214-221.View ArticlePubMedGoogle Scholar
- Lowenhaupt K, Wang PJ, Horan T, Lauffenburger DA, Sarkar CA: Parsing the effects of binding, signaling, and trafficking on the mitogenic potencies of granulocyte colony-stimulating factor analogues. Biotechnol Prog. 2003, 19: 955-964.View ArticlePubMedGoogle Scholar
- Veronese FM, Mero A: The impact of PEGylation on biological therapies. BioDrugs. 2008, 22: 315-329.View ArticlePubMedGoogle Scholar
- Scholz M, Engel C, Apt D, Sankar SL, Goldstein E, Loeffler M: Pharmacokinetic and Pharmacodynamic modelling of the novel human G-CSF derivative Maxy-G34 and Pegfilgrastim in the rat. Cell Prolif. 2009, 42 (6): 823-37.View ArticlePubMedGoogle Scholar
- 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: 903-909.View ArticlePubMedGoogle Scholar
- Grigg A, Solal-Celigny P, Hoskin P, Taylor K, McMillan A, Forstpointner R, Bacon P, Renwick J, Hiddemann W: Openlabel, randomized study of pegfilgrastim vs. daily filgrastim as an adjunct to chemotherapy in elderly patients with non-Hodgkin’s lymphoma. Leuk Lymphoma. 2003, 44: 1503-1508.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Pinto L, Liu Z, Doan Q, Bernal M, Dubois R, Lyman G: Comparison of pegfilgrastim with filgrastim on febrile neutropenia, grade IV neutropenia and bone pain: a meta-analysis of randomized controlled trials. Curr Med Res Opin. 2007, 23: 2283-2295.View ArticlePubMedGoogle Scholar
- Ziepert M, Schmits R, Trumper L, Pfreundschuh M, Loeffler M: Prognostic factors for hematotoxicity of chemotherapy in aggressive non-Hodgkin’s lymphoma. Ann Oncol. 2008, 19: 752-762.View ArticlePubMedGoogle Scholar
- Kroschinsky F, Holig K, Ehninger G: The role of pegfilgrastim in mobilization of hematopoietic stem cells. Transfus Apher Sci. 2008, 38 (3): 237-44.View ArticlePubMedGoogle Scholar
- Scholz M, Ackermann M, Engel C, Emmrich F, Loeffler M, Kamprad M: A pharmacokinetic model of filgrastim and pegfilgrastim application in normal mice and those with cyclophosphamide-induced granulocytopaenia. Cell Prolif. 2009, 42 (6): 813-22.View ArticlePubMedGoogle Scholar
- Hartmann F, Zeynalova S, Nickenig C, Reiser M, Lengfelder E, Duerk H, de Witt M, Schubert J, Loeffler M, Pfreundschuh M: Peg-filgrastim (Peg-F) on day 4 of (R-)CHOP-14 chemotherapy compared to day 2 in elderly patients with diffuse large B-cell lymphoma (DLBCL): Results of a randomized trial of the German high-grade non-Hodgkin’s lymphoma study group (DSHNHL). J Clin Oncol, ASCO Ann Meeting Proc Part I. 2007, 25 (18S): 19511-Google Scholar
- 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: Targets & Therapy. 2009, 3: 27-37.Google Scholar
- Engel C, Scholz M, Loeffler M: A computational model of human granulopoiesis to simulate the hematotoxic effects of multicycle polychemotherapy. Blood. 2004, 104: 2323-2331.View ArticlePubMedGoogle Scholar
- Scholz M, Engel C, Loeffler M: Modelling human granulopoiesis under poly-chemotherapy with G-CSF support. J Math Biol. 2005, 50: 397-439.View ArticlePubMedGoogle Scholar
- Scholz M, Engel C, Loeffler M: Model-based design of chemotherapeutic regimens that account for heterogeneity in leucopoenia. Br J Haematol. 2006, 132: 723-735.View ArticlePubMedGoogle Scholar
- Tsuji K, Ebihara Y: Expression of G-CSF receptor on myeloid progenitors. Leuk Lymphoma. 2001, 42 (6): 1351-7.View ArticlePubMedGoogle Scholar
- Prosper F, Stroncek D, McCarthy JB, Verfaillie CM: Mobilization and homing of peripheral blood progenitors is related to reversible downregulation of alpha4 beta1 integrin expression and function. J Clin Invest. 1998, 101 (11): 2456-67.PubMed CentralView ArticlePubMedGoogle Scholar
- Kuwabara T, Uchimura T, Takai K, Kobayashi H, Kobayashi S, Sugiyama Y: Saturable uptake of a recombinant human granulocyte colony-stimulating factor derivative, nartograstim, by the bone marrow and spleen of rats in vivo. J Pharmacol Exp Ther. 1995, 273: 1114-1122.PubMedGoogle Scholar
- Kota J, Machavaram KK, McLennan DN, Edwards GA, Porter CJ, Charman SA: Lymphatic absorption of subcutaneously administered proteins: influence of different injection sites on the absorption of darbepoetin alfa using a sheep model. Drug Metab Dispos. 2007, 35: 2211-2217.View ArticlePubMedGoogle Scholar
- Kuwabara T, Kobayashi S, Sugiyama Y: Pharmacokinetics and pharmacodynamics of a recombinant human granulocyte colony stimulating factor. Drug Metab Rev. 1996, 28: 625-658.View ArticlePubMedGoogle Scholar
- Lenhoff S, Rosberg B, Olofsson T: Granulocyte interactions with GM-CSF and G-CSF secretion by endothelial cells and monocytes. Eur Cytokine Network. 1999, 10 (4): 525-532.Google Scholar
- Layton JE, Hockman H, Sheridan WP, Morstyn G: Evidence for a novel in vivo control mechanism of granulopoiesis: mature cell-related control of a regulatory growth factor. Blood. 1989, 74: 1303-1307.PubMedGoogle Scholar
- 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: 147-158.View ArticlePubMedGoogle Scholar
- Roskos LK, Lum P, Lockbaum P, Schwab G, Yang BB: Pharmacokinetic / pharmacodynamic modeling of pegfilgrastim in healthy subjects. J Clin Pharmacol. 2006, 46: 747-757.View ArticlePubMedGoogle Scholar
- 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: 321-342.View ArticlePubMedGoogle Scholar
- Falanga A, Marchetti M, Evangelista V, Manarini S, Oldani E, Giovanelli S, Galbusera M, Cerletti C, Barbui T: Neutrophil activation and hemostatic changes in healthy donors receiving granulocyte colony-stimulating factor. Blood. 1999, 93: 2506-2514.PubMedGoogle Scholar
- Levesque JP, Takamatsu Y, Nilsson SK, Haylock DN, Simmons PJ: Vascular cell adhesion molecule-1 (CD106) is cleaved by neutrophil proteases in the bone marrow following hematopoietic progenitor cell mobilization by granulocyte colony-stimulating factor. Blood. 2001, 98: 1289-1297.View ArticlePubMedGoogle Scholar
- Shimazaki C, Uchiyama H, Fujita N, Araki S, Sudo Y, Yamagata N, Ashihara E, Goto H, Inaba T, Haruyama H: Serum levels of endogenous and exogenous granulocyte colony-stimulating factor after autologous blood stem cell transplantation. Exp Hematol. 1995, 23: 1497-1502.PubMedGoogle Scholar
- Steinman RA, Tweardy DJ: Granulocyte colony-stimulating factor receptor mRNA upregulation is an immediate early marker of myeloid differentiation and exhibits dysfunctional regulation in leukemic cells. Blood. 1994, 83: 119-127.PubMedGoogle Scholar
- Tkatch LS, Rubin KA, Ziegler SF, Tweardy DJ: Modulation of human G-CSF receptor mRNA and protein in normal and leukemic myeloid cells by G-CSF and retinoic acid. J Leukoc Biol. 1995, 57: 964-971.PubMedGoogle Scholar
- Srinivasa SP, Doshi PD: Extracellular signal-regulated kinase and p38 mitogen-activated protein kinase pathways cooperate in mediating cytokine-induced proliferation of a leukemic cell line. Leukemia. 2002, 16 (2): 244-53.View ArticlePubMedGoogle Scholar
- Lohrmann HP, Schreml W: Cytotoxic Drugs and the Granulopoietic System. 1982, Berlin: Springer VerlagView ArticleGoogle Scholar
- Scholz M, Gross A, Loeffler M: A biomathematical model of human thrombopoiesis under chemotherapy. J Theor Biol. 2010, 264 (2): 287-300.View ArticlePubMedGoogle Scholar
- van der Auwera P, Platzer E, Xu ZX, Schulz R, Feugeas O, Capdeville R, Edwards DJ: Pharmacodynamics and pharmacokinetics of single doses of subcutaneous pegylated human g-csf mutant (Ro 25-8315) in healthy volunteers: comparison with single and multiple daily doses of filgrastim. Am J Hematol. 2001, 66 (4): 245-251.View ArticlePubMedGoogle Scholar
- Borleffs JC, Bosschaert M, Vrehen HM, Schneider MM, van Strijp J, Small MK, Borkett KM: Effect of escalating doses of recombinant human granulocyte colony-stimulating factor (filgrastim) on circulating neutrophils in healthy subjects. Clin Ther. 1998, 20 (4): 722-736.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Varki R, Pequignot E, Leavitt MC, Ferber A, Kraft WK: A glycosylated recombinant human granulocyte colony stimulating factor produced in a novel protein production system (AVI-014) in healthy subjects: a first-in human, single dose, controlled study. BMC Clin Pharmacol. 2009, 9 (2):Google Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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, Truemper L, Balleisen L, Liersch R, Metzner B, Hartmann F, Glass B, Poeschel V, Schmitz N, Ruebe C, Feller AC, Loeffler M, German High-Grade Non-Hodgkin’s Lymphoma Study Group: 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.View ArticlePubMedGoogle Scholar
- Zwick C, Hartmann F, Zeynalova S, Poeschel V, Nickenig C, Reiser M, Lengfelder E, Peter N, Schlimok G, Schubert J, Schmitz N, Loeffler M, Pfreundschuh M, German High-Grade Non-Hodgkin Lymphoma Study Group: Randomized comparison of pegfilgrastim day 4 versus day 2 for the prevention of chemotherapy-induced leukocytopenia. Ann Oncol. 2011, 22 (8): 1872-1877.View ArticlePubMedGoogle Scholar
- Rechenberg I: Evolutionsstrategie ’94. 1994, Stuttgart: frommann-holzboogGoogle Scholar
- Schwefel HP: Evolution strategies: A family of nonlinear optimization techniques based on imitating some principles of organic evolution. Ann Oper Res. 1984, 1: 65-167.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Shochat E, Rom-Kedar V, Segel LA: G-CSF control of neutrophils dynamics in the blood. Bull Math Biol. 2007, 69 (7): 2299-338.View ArticlePubMedGoogle Scholar
- Foley C, Mackey MC: Mathematical model for G-CSF administration after chemotherapy. J Theor Biol. 2009, 257 (1): 27-44.View ArticlePubMedGoogle Scholar
- Waladkhani AR: Pegfilgrastim: a recent advance in the prophylaxis of chemotherapy-induced neutropenia. Eur J Cancer Care. 2004, 13: 371-379.View ArticleGoogle Scholar
- Wiczling P, Lowe P, Pigeolet E, Luedicke F, Balser S, Krzyzanski W: Population pharmacokinetic modelling of filgrastim in healthy adults following intravenous and subcutaneous administrations. Clin Pharmacokinet. 2009, 48 (12): 817-26.View ArticlePubMedGoogle Scholar
- Wichmann H-E, Loeffler M: Mathematical Modeling of Cell Proliferation: Stem Cell Regulation in Hemopoiesis. 1985, Boca Raton: sCRC PressGoogle Scholar
- Schmitz S, Franke H, Loeffler M, Wichmann HE, Diehl V: Model analysis of the contrasting effects of GM-CSF and G-CSF treatment on peripheral blood neutrophils observed in three patients with childhood-onset cyclic neutropenia Brit. J Haematol. 1996, 95: 616-625.View ArticleGoogle Scholar
- Dale DC, Fauci AS, Wolff SM: Alternate-day prednisone. TnE J Med. 1974, 291 (22): 1154-1158.Google Scholar

## Copyright

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