Pharmacokinetic and dynamic modelling of GCSF derivatives in humans.  
Jump to Full Text  
MedLine Citation:

PMID: 22846180 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

BACKGROUND: The human granulocyte colonystimulating factor (GCSF) 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 GCSF derivatives Filgrastim and Pegfilgrastim. RESULTS: Model equations are based on our physiological understanding of the drugs which are delayed absorption of GCSF when applied to the subcutaneous tissue, dosedependent 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 cellkinetic model of granulopoiesis was essentially preserved, except for a few additional assumptions and simplifications. We assumed a delayed action of GCSF 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 GCSF which interact via concurrent receptor binding. Unknown pharmacokinetic or cellkinetic parameters were determined by fitting the predictions of the model to available datasets of GCSF 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 GCSF derivatives. We aim to apply the model to a large variety of chemotherapy regimen in the future in order to optimize corresponding GCSF schedules or to individualize GCSF treatment according to the granulotoxic risk of a patient. 
Authors:

Markus Scholz; Sibylle Schirm; Marcus Wetzler; Christoph Engel; Markus Loeffler 
Publication Detail:

Type: Journal Article; Research Support, NonU.S. Gov't Date: 20120730 
Journal Detail:

Title: Theoretical biology & medical modelling Volume: 9 ISSN: 17424682 ISO Abbreviation: Theor Biol Med Model Publication Date: 2012 
Date Detail:

Created Date: 20121128 Completed Date: 20130621 Revised Date: 20130712 
Medline Journal Info:

Nlm Unique ID: 101224383 Medline TA: Theor Biol Med Model Country: England 
Other Details:

Languages: eng Pagination: 32 Citation Subset: IM 
Affiliation:

Institute for Medical Informatics, Statistics and Epidemiology, University of Leipzig, Haertelstrasse 1618, 04107 Leipzig, Germany. markus.scholz@imise.unileipzig.de 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:

Bone Marrow
/
drug effects Granulocyte ColonyStimulating Factor / pharmacokinetics*, pharmacology Granulocytes / drug effects Hematopoiesis / drug effects Humans Models, Biological* Recombinant Proteins / pharmacokinetics, pharmacology 
Chemical  
Reg. No./Substance:

0/Recombinant Proteins; 0/pegfilgrastim; 121181531/Filgrastim; 143011727/Granulocyte ColonyStimulating Factor 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): Theor Biol Med Model Journal ID (isoabbrev): Theor Biol Med Model ISSN: 17424682 Publisher: BioMed Central 
Article Information Download PDF Copyright ©2012 Scholz et al.; licensee BioMed Central Ltd. openaccess: Received Day: 17 Month: 2 Year: 2012 Accepted Day: 12 Month: 6 Year: 2012 collection publication date: Year: 2012 Electronic publication date: Day: 30 Month: 7 Year: 2012 Volume: 9First Page: 32 Last Page: 32 PubMed Id: 22846180 ID: 3507764 Publisher Id: 17424682932 DOI: 10.1186/17424682932 
Pharmacokinetic and dynamic modelling of GCSF derivatives in humans  
Markus Scholz12  Email: markus.scholz@imise.unileipzig.de 
Sibylle Schirm1  Email: sschirm@imise.unileipzig.de 
Marcus Wetzler1  Email: marcus.wetzler@imise.unileipzig.de 
Christoph Engel12  Email: christoph.engel@imise.unileipzig.de 
Markus Loeffler12  Email: markus.loeffler@imise.unileipzig.de 
1Institute for Medical Informatics, Statistics and Epidemiology, University of Leipzig, Haertelstrasse 1618, 04107 Leipzig, Germany 

2LIFE  Leipzig Research Center for Civilization Diseases, University of Leipzig, PhilippRosenthalStrasse 27, 04103 Leipzig, Germany 
The human granulocyte colonystimulating factor (GCSF) 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}]. GCSF 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 (nonmyeloablative) chemotherapies, the haematopoietic system usually recovers without further medication. But, GCSF can significantly speed up this process allowing dose and timeintensifications of multicycle chemotherapies [^{11},^{12}]. While platelets and red blood cells can show cumulative toxicity during the course of intensified regimen, appropriate GCSF 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 GCSF are available now. The first generation of GCSF pharmaceuticals were recombinant derivatives such as Filgrastim (nonglycosylated) or Lenograstim (glycosylated). Both derivatives are virtually identical to endogenously produced GCSF [^{14}^{}^{16}] but Filgrastim is more frequently used in clinical trials. Filgrastim is eliminated by both, renal elimination and specific degradation mediated by GCSF receptors or neutrophil elastase [^{17}^{}^{22}]. This results in a short halflife in vivo requiring multiple injections during one cycle of chemotherapy. As next generation GCSF derivative, Pegfilgrastim (pegylated Filgrastim) was developed in order to improve the pharmaceutical properties of Filgrastim. Indeed, Pegfilgrastim shows a remarkably prolonged halflife 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 halflife, 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 GCSF derivatives by additional pegylations (e.g. MaxyG34, [^{29}]). Although the application of pegylated GCSF 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 GCSF was applied [^{34}]. Additionally, for the purpose of stem cell mobilization, Filgrastim is not inferior compared to Pegfilgrastim but has less severe sideeffects [^{35}].
In view of the highly differing pharmacokinetic properties of the available GCSF derivatives, we constructed pharmacokinetic models of the GCSF derivatives Filgrastim, Pegfilgrastim and the novel MaxyG34 in mice and rats [^{29},^{36}]. Our aim was to identify basic pharmacokinetic model mechanisms especially with respect to the degradation of GCSF in vivo and to compare the resulting pharmacokinetic parameters between the GCSF derivatives. We now aim to translate these model insights to humans.
Effectiveness of GCSF treatment depends on many variable therapy parameters such as applied chemotherapy, individual factors, GCSF derivative used, and especially, its dosing and timing schedule [^{4},^{34},^{37},^{38}]. Chemotherapy induced neutropenia and GCSF induced granulocytosis via different modes of action in combination with a strong specific elimination of GCSF mediated by circulating granulocytes results in complex dynamics of both GCSF serum concentration and circulating granulocytes as well. In consequence, the optimization of GCSF treatment is a nontrivial task and cannot be performed solely on the basis of clinical trials. We showed in the past that biomathematical cellkinetic models of granulopoiesis under chemotherapy and GCSF 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 cellkinetic 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 GCSF serum concentrations for virtually all datasets published in the literature. We also discuss how the model can be used to optimize GCSF scheduling of chemotherapies.
Our cell kinetic model of granulopoiesis is an ordinary differential equations system modelling the timedependent content of and the fluxes between the following cell compartments: S (pluripotent stem cells), CG (colony forming units of granulocytes and macrophages), PGB (proliferating granulopoietic blasts), MGB (maturing granulopoietic blasts  subdivided into metamyelocytes (G4), banded granulocytes (G5) and segmented granulocytes (G6)) 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 GCSF 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 GCSF are regulated by mature cells. We also modelled a subcutaneous compartment in which GCSF 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 cellkinetic 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 Figure 1. A complete set of model equations and parameters is presented in the Additional file 1.
For unperturbed granulopoiesis, the model is autonomous and has a single fixedpoint (steadystate). This fixedpoint appears to be stable for the parameter set which we propose later, i.e. transient perturbations result in damped oscillations until the fixedpoint is reestablished. Permanent perturbations such as constant GCSF stimulation or chemotherapy damage result in new fixedpoints 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 steadystate at the beginning of any treatment. That is, initial values equals normal values.
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}].
Most of our cell kinetic parameters such as amplification or transition times are regulated between a minimum and a maximum value by a so called Zfunction, which is a function of another quantity such as circulating GCSF. We use the following class of Zfunctions:
(1)
Y=ZX=Ymax−Ymax−Ymine−lnYmax−YminYmax−YnorXbYforYmin<Ynor<YmaxorYmin>Ynor>Ymax,YnorforYmin=Ynor=Ymax 
where X is the regulator of the quantity Y and Y^{min}, Y^{nor}, Y^{max}, b_{Y }are the parameters of the Zfunction for the regulation of Y . The parameter b_{Y }defines the steepness of the function and is called sensitivity parameter in the following.
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 GCSF or chemotherapy actions.
Our pharmacokinetic model of GCSF is essentially based on models developed for mice and rats [^{29},^{36}]. The underlying model assumptions are as follows:
1. The pharmacokinetic model contains three compartments in which GCSF is present: The subcutaneous compartment CG−CSFsc in which GCSF is usually injected, the central compartment CG−CSFcent in which GCSF is haematologically active and a peripheral compartment CG−CSFper representing reversible bindings of GCSF (e.g. protein binding [^{44}]).
2. Subcutaneously injected GCSF results in a delayed influx of GCSF 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 GCSF 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 GCSF 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 GCSF derivatives is dosedependent [^{29}], we assume that a part of the applied GCSF is removed from the subcutaneous compartment without entering the central compartment. This is modelled by a MichaelisMenten kinetic within the first subcompartment of the subcutaneous tissue.
6. GCSF 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 MichaelisMenten 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 GCSF 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 GCSF receptors are within the bone marrow which dynamics are somewhat different from the dynamics of mature granulocytes.
7. Differences in GCSF derivatives are model by different model parameters rather than differences in model structure. In view of the high similarity of Filgrastim and endogenous GCSF [^{14}^{}^{16}], we assume that the pharmacokinetic and dynamic parameters of Filgrastim and endogenous GCSF 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 GCSF pharmaceuticals in order to model observed differences in GCSF serum dynamics after subcutaneous applications [^{23}^{}^{25}]. Pharmacodynamik differences are modelled by different parameterizations of the GCSF mediated regulatory mechanisms (Zfunctions). The latter is motivated by experimental results suggesting a reduced receptor binding affinity of pegylated GCSF [^{26}^{}^{29}].
Based on these assumptions, we formulate the pharmacokinetic model equations in the next section. A schematic structure of the model can be found in Figure 2.
According to [^{40}], the relative GCSF production PG−CSFendo is a Zfunction of the relative content of segmented granulocytes in bone marrow and granulocytes in circulation.
(2)
PG−CSFendo=ZωG6CG6+ωGRACGRAωG6CG6nor+ωGRACGRAnor 
where ω_{G6} and ω_{GRA} are weighting parameters. It holds that PG−CSFendo_nor=1 whereas all other parameters of the Zfunctions are free parameters.
Exogenous GCSF applications are modelled by an injection function. Let
Hv(t)=0:t≤01:t>0 
be the Heavisidefunction, then the injection function reads as follows:
(3)
PG−CSFexo=∑i=1LdG−CSFtiHvt−ti−Hvt−ti−tinftinf 
where t_{i }≥ 0 (i = 1,…,L) are the time points of GCSF 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 }= 5s).
The injection function is specific for each GCSF 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 PG−CSFexo_sc_fil, PG−CSFexo_iv_fil for subcutaneous and intravenous Filgrastim injections and PG−CSFexo_sc_peg and PG−CSFexo_iv_peg for subcutaneous and intravenous Pegfilgrastim injections respectively.
The subcutaneous compartment is divided into two subcompartments sc_1 and sc_2 where the efflux of the first compartment is the influx to the second compartment. GCSF is applied to the first subcompartment (second term of (4)). In the first subcompartment there is a dosedependent loss of GCSF modelled by a MichaelisMenten kinetic (third term of (4)). For Filgrastim injections it holds that
(4)
ddtCG−CSFsc_1=PG−CSFexo_sc_fil−kscFCG−CSFsc_1−vmaxFCG−CSFsc_1kmF+CG−CSFsc_1 
(5)
ddtCG−CSFsc_2=kscFCG−CSFsc_1−CG−CSFsc_2 
with the initial values CG−CSFsc_1(0)=CG−CSFsc_2(0)=0. For Pegfilgrastim injections the first term of the righthand side of (4) is substituted by PG−CSFexo_peg_sc. Likewise, the Filgrastim parameters kscF, vmaxF and kmF are substituted by corresponding Pegfilgrastim parameters.
For Filgrastim injections it holds that
(6)
ddtCG−CSFcent=PG−CSFrefPG−CSFendo+PG−CSFexo_fil_iv+kscFCG−CSFsc_2−kuFCG−CSFcent−kcpFCG−CSFcent+kpcFCG−CSFper−vmaxGRA_FCG−CSFcentkmGRA_F+CG−CSFcentCGRArel. 
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 (MichaelisMenten kinetic) which is proportional to the relative granulocyte concentration CGRArel. The corresponding equation for Pegfilgrastim is the same except for the endogenous production, which is zero, the intravenous injection function which is substituted by PG−CSFexo_peg_iv, the parameters and the initial value which is again zero.
The relative endogenous production function of GCSF is multiplied by a constant PG−CSFref in order to adjust the normal value of CG−CSFcent to a reference amount of GCSF which is the product of the distribution volume VDF and the reference GCSF serum concentration CG−CSFcent_ref which will be determined later on the basis of measurements available from the literature. Assume that CG−CSFcent(0)=CG−CSFcent_nor=VDFCG−CSFcent_ref, the parameter PG−CSFref can be calculated exploiting the steadystate conditions of equation (6), i.e. PG−CSFendo(0)=1, PG−CSFexo_fil_iv(0)=0 and −kcpFCG−CSFcent(0)+kpcFCG−CSFper(0)=0:
(7)
PG−CSFref=VDFCG−CSFcent_refkuF+vmaxGRA_FkmGRA_F+VDFCG−CSFcent_ref 
For both GCSF derivatives, we have
(8)
ddtCG−CSFper=kcpCG−CSFcent−kpcCG−CSFper 
(9)
CG−CSFper(0)=CG−CSFper_nor=VDkcpkpcCG−CSFcent_ref 
where the parameters k_{cp}, k_{pc} and V_{D} are specific for Filgrastim and Pegfilgrastim respectively.
The pharmacodynamic model describes the dynamics of the bone marrow cell stages, circulating granulocytes, GCSF, corresponding regulations and the action of chemotherapy as well. In Table 1 we present the major model compartments and its regulatory features. A complete set of equations can be found in the Additional file 1.
We used the same pharmacodynamic model of GCSF as presented in [^{40}] except for a few simplifications and additional assumptions which we will discuss now:
The cytokine GMCSF (granulocyte macrophage colonystimulating factor) was no longer considered in order to simplify the model. Since the endogenous productions of GMCSF and GCSF in our former model were both related to the demand of bone marrow cell stages, the time courses of endogenous GCSF and GMCSF after external perturbations are similar, making the cytokines undistinguishable from the modelling point of view (compare [^{40}]). Therefore, only GCSF was considered in the present model. It replaces GMCSF regarding the regulation of the CG compartment. This simplification is also biologically plausible since GCSF 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 GCSF via a Zfunction. 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 GCSF receptor binding affinities [^{26}^{}^{29}], we assume different regulatory Zfunctions for both derivatives. But we assumed the same pharmacodynamic parameters for Filgrastim and endogenously produced GCSF.
While in our former model version the GCSF concentration instantaneously affected the value of the Zfunctions, we now introduce a time delay regarding GCSF action. This was motivated by transduction network analyses which revealed a delayed response of the transcriptome to GCSF stimulations [^{57}]. The delay is modelled by a cascade of four subcompartments (see section “Basic model mechanisms”). The efflux of the last subcompartment is the delayed GCSF concentration CG−CSFcent_rel_del which is the new argument of our Zfunctions.
(10)
ddtCG−CSF1=CG−CSFcent_rel−DG−CSF·CG−CSF1 
(11)
ddtCG−CSFi=DG−CSFCG−CSFi−1−CG−CSFii=2,…,4 
with CG−CSFcent_rel_del=DG−CSF·CG−CSF4. The delay parameter D_{G−CSF }is specific for Pegfilgrastim and Filgrastim but constant for all Zfunctions. Hence, (1) reads
(12)
Y=ZYCG−CSFcent_rel_del 
for all quantities Y regulated by GCSF. At this, the normal value of the Zfunction Y^{nor }refers to CG−CSFnor for endogenous GCSF or Filgrastim and to an absolute amount of 1μg for Pegfilgrastim.
In case of Pegfilgrastim injections, Pegfilgrastim and endogenous GCSF compete with respect to receptor binding. To model this process, the Zfunctions of Pegfilgrastim and Filgrastim were added using a weighting factor ω_{P }which is again a Zfunction of the quotient of the two doses with minimum 0 and maximum 1.
(13)
ωP=ZωPCG−CSFcent_rel_del_peg/CG−CSFcent_rel_del_filwithωPmin=0andωPmax=1 
(14)
ZY=ωP·ZYCG−CSFcent_rel_del_peg+1−ωP·ZYCG−CSFcent_rel_del_fil 
where Y is an arbitrary regulated quantity such as a transition time or an amplification. For all these quantities we assumed the same Zfunction 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.
Since our model is a model of granulopoiesis, we can only make predictions for absolute neutrophil counts. However, in clinical practice often only leukocytes are available. In our former model version, we assumed proportionality of leukocytes and absolute neutrophil count which is only roughly correct [^{38},^{40}]. To be more precise, we now calculate the leukocyte count as the sum of lymphocytes and granulocytes. To avoid a full model of lymphopoiesis, we modelled the reduced lymphocyte count under chemotherapy by an exponential function of the corresponding toxicity function.
(15)
CWBC(t)≈cLYexp−ΨLY(t)+cGRACGRA(t)CGRAnor 
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.)
Since the precise structure of the toxicity function depends on the schedule of the chemotherapy, we demonstrate the construction of the toxicity functions using six cycles of CHOP therapy with a cycle duration of 14 days (6xCHOP14) as example. During CHOP14 therapy, the cytotoxic drugs cyclophosphamide, doxorubicine and vincristine are applied concomitantly at the first day of each cycle. Since contributions of these drugs to granulotoxicity cannot be separated, we assume a unique toxicity function for this drug combination. Hence, the chemotherapy injection function Ψ6×CHOP−14inj reads as follows:
Ψ6×CHOP−14inj=∑i=16Hvt−14d·i−1−Hvt−14d·i−1−t6×CHOP−14inft6×CHOP−14inf 
which is in analogy to (3). Here, no dose parameter is required since intensity of damage is defined by the toxicity parameters later. We set t6×CHOP−14inf=1d. The delayed action of chemotherapy can now be modelled in analogy to (10), (11):
ddtΨ6×CHOP−14X_1=Ψ6×CHOP−14inj−D6×CHOP−14X·Ψ6×CHOP−14X_1ddtΨ6×CHOP−14X_i=D6×CHOP−14X·Ψ6×CHOP−14X_(i−1)−Ψ6×CHOP−14X_i 
where 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:
Ψ6×CHOP−14LY=kLY·D6×CHOP−14LY·Ψ6×CHOP−14LY_4Ψ6×CHOP−14GRA_Y=kY·D6×CHOP−14GRA·Ψ6×CHOP−14GRA_4 
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.
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.
To address these challenges, we established the following stepwise fitting procedure keeping parameters identified at the previous step constant.
1. Pharmacokinetic model parameters were determined on the basis of available cytokine dynamics after GCSF 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 (GCSF 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 file 1. Not all data sets were used to fit parameters. A few data sets were kept in reserve in order to validate the model.
Data sets were collected from the literature by an extensive search. For our modelling purposes, closemeshed time series of GCSF 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 GCSF in healthy volunteers and conventional chemotherapies of different diseases with or without GCSF prophylaxis.
Additionally, we can rely on own clinical trials data for which one of us (Markus Loeffler) is the responsible biostatistician or for which we have cooperation agreements. Leukocyte raw data under chemotherapy are available from published studies of the German High Grade NonHodgkinLymphoma 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 Table 2.
As mentioned above, unknown parameters of the model were determined by fitting the predictions of the model to available clinical data minimizing the L^{1} distance between logarithmized model and logarithmized median of data. More precisely, we have
(16)
∫t0t1log(fmodel(t,k))−logfdata(t)dt→mink, 
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 logf_{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 nondeterministic 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 selfadapting 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 GCSF 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}].
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 file 1.
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.
Given the model equations of section “Pharmacokinetic model equations”, we determined the pharmacokinetic parameters for Filgrastim and Pegfilgrastim separately by fitting the model to time series of GCSF serum concentrations after single or multiple application of one of the two GCSF derivatives. Parameter estimates are shown in Table 3.
Combined with our pharmacodynamic model (see section “The new pharmacodynamic model”), these parameter estimates resulted in a good fit of all model scenarios. Examples are presented in Figure 4. A complete list of all scenarios is presented in the Additional file 1.
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 GCSF 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 GCSF injections. Examples are shown in Figure 4. All other scenarios of GCSF application can be found in the Additional file 1. 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 file 1.
Pharmacodynamic differences between Filgrastim and Pegfilgrastim can be traced back to differences of the regulation functions. In general, compared to Pegfilgrastim, the regulation functions of Filgrastim express a higher sensitivity regarding changes of the GCSF concentration, that is, greater slopes and higher values under maximum stimulation (see Figure 5 for an example).
In contrast to our former model of granulopoiesis, we assumed a delayed effect of GCSF (assumption 4 in section “New pharmacodynamic model assumptions”). Consequences of this assumption are studied in Figure 6 on the basis of model simulations of single Filgrastim injections.
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}]).
Chemotherapy was modelled by a transient depletion of bone marrow cell stages. In contrast to our former model of granulopoiesis, we assumed a delay of the bone marrow depletion. The effect of this delay is studied in Figure 7.
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).
Data of the CHOP regimen were utilized to fit both, the pharmacokinetic and dynamic model and the set of specific toxicity parameters as well. This set of toxicity parameters was valid for all GCSF schedules applied as supportive therapy for CHOP. Results of these scenarios are shown in Figure 8.
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 file 1.
Since we assumed a simplified model of lymphocyte toxicity under chemotherapy, we can estimate the ratio of granulocytes and leukocytes under therapy, offering a possibility to validate the model. Results of the CHOP14 regimen with Filgrastim application at day 4 to 13 are displayed in Figure 9. This ratio was estimated to be clearly not constant varying between 68% and 98%.
A few datasets were kept in reserve in order to validate our model. The phase 1 data of Varki 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 GCSF serum levels as well. Both scenarios fit well with our model prediction (see Figure 10). No additional parameter fittings were performed to model these scenarios.
In the present paper, we developed an ordinary differential equations model of human granulopoiesis under chemotherapy and GCSF 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 GCSF derivatives (Filgrastim and Pegfilgrastim). At this, the pharmacokinetic model was adopted from similar models developed for GCSF 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 GCSF into healthy volunteers or patients under CHOP chemotherapy. We were able to describe the differences between the GCSF 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 GCSF applications. Published models comprise for example pure pharmacokinetic models [^{25},^{46}], pharmacokinetic and dynamic models of GCSF application on the cellular level [^{49}], in healthy volunteers [^{50},^{51},^{70}], for the treatment of cyclic neutropenia [^{71}], for highdose 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 GCSF 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. GCSF 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 GCSF 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 firstorder differential equations rather than a fixed time delay. We showed in the past that this kind of modelling is equivalent to a Gammadistributed 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 GCSF has a dosedependent 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 GCSF molecules, we modelled a first order transition between the blood compartment and a peripheral compartment [^{46}].
Our model assumptions regarding endogenous production of GCSF are speculative. We assumed that the production is regulated between a minimum and a maximum value in dependence on bone marrow cellularity [^{47}]. In steadystate, the production is constant in order to sustain a fixed serum concentration determined by averaged data from the literature.
Degradation of GCSF was modelled by two independent processes, an unspecific renal clearance modelled by a first order transition and a specific degradation via neutrophil elastase or GCSF 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 MichaelisMenten 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 GCSF receptors are also present in bone marrow progenitors and precursors [^{42}], it appears to be less appropriate for the receptormediated 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 GCSF by both, mature neutrophils and precursors which did not significantly improve the quality of our model predictions.
Variable numbers of GCSF receptors per cell were observed in dependence on the GCSF 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 GCSF 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 GCSF derivative, namely MaxyG34, which is a novel GCSF 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 GCSF, it was necessary to attach a pharmacodynamic model of GCSF applications. Since there is some evidence that the pegylations of the Pegfilgrastim molecule interact with its binding affinity to the GCSF receptor [^{26}^{}^{29}], we decided to assume different regulation functions for Filgrastim and Pegfilgrastim. Hence, Zfunctions 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 Zfunctions for Filgrastim and endogenous GCSF. Due to the fact that Filgrastim/endogenous GCSF and Pegfilgrastim were assumed to be different regarding Zfunctions, it was necessary to merge the superimposing effects of Pegfilgrastim and endogenous GCSF 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 GCSF in the system. If the Pegfilgrastim concentration is high or low, then the system is mainly influenced by the Zfunction of Pegfilgrastim or Filgrastim respectively. Although this assumption is plausible, the complete regulation mechanism via the combined Zfunctions 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 GCSF 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 semiexplicite 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 GCSF 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 GCSF 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 GCSF 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 cellkinetic and toxicity parameters is available. Rough ranges for transition times and amplification rates in steadystate or under stimulation by GCSF were obtained from the literature. But especially values under minimal stimulation, the sensitivity parameters of the Zfunctions 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 GCSF serum levels after application in combination with granulocyte or leukocyte counts after chemotherapy and different GCSF 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 GCSF 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 GCSF 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 cellkinetic, pharmacodynamic and toxicity parameters also resulted in a good fit of the time courses of granulocytes and leukocytes after application of chemotherapy, GCSF 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 daywise scale but might be unable to explain shortterm or transient effects of GCSF applications e.g. on the scale of hours. Modelling such shortterm effects would require a better database, since almost all available time courses of granulocytes and leukocytes were measured at most at a daywise scale. However, in order to make predictions regarding the efficiency of different GCSF schedules, we are also more interested in the longterm dynamics of granulocytes and leukocytes in the course of the therapy rather than shortterm 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 MaxyG34 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 GCSF 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 GCSF 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 GCSF 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 GCSF serum concentrations, bone marrow cell stages and mature cell counts in circulation under different GCSF schedules, allowing to optimize GCSF treatment. We will exploit the clinically relevant applications of our model in the near future.
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 GCSF serum concentrations, granulocytes and leukocytes of patients treated with GCSF 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 GCSF, bone marrow cell stages and circulating granulocytes or leukocytes of yet untested GCSF schedules. The model is currently applied in the planning phase of clinical trials in order to optimize GCSF treatment.
The authors declare that they have no competing interests.
Model development: S.S., M.W., C.E., M.S. Parameter estimation and model simulations: S.S. Paper writing: M.S. All authors contributed to the discussion and the paper writing. All authors read and approved the final manuscript.
Supplement Material. Complete list of model equations, complete list of model parameters, additional model and data comparisons, sensitivity analysis [^{8},^{40},^{64},^{78}^{}^{80}].
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.
References
Crawford J,Pegfilgrastim administered once per cycle reduces incidence of chemotherapyinduced neutropeniaDrugsYear: 200262Suppl 1899812479597  
Dale D,Current management of chemotherapyinduced neutropenia: the role of colonystimulating factorsSemin OncolYear: 2003303914508714  
Siena S,Secondino S,Giannetta L,Carminati O,Pedrazzoli P,Optimising management of neutropenia and anaemia in cancer chemotherapyadvances in cytokine therapyCrit Rev Oncol HematolYear: 200348S39s4714563520  
Wunderlich A,Kloess M,Reiser M,Rudolph C,Truemper L,Bittner S,Schmalenberg H,Schmits R,Pfreundschuh M,Loeffler M,German HighGrade NonHodgkin’s Lymphoma Study Group (DSHNHL): Practicability and acute haematological toxicity of 2 and 3weekly CHOP and CHOEP chemotherapy for aggressive nonHodgkin’s lymphoma: results from the NHLB trial of the German HighGrade NonHodgkin’s Lymphoma Study Group (DSHNHL)Ann OncolYear: 20031488189312796026  
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 colonystimulating factor: effects on normal and leukemic myeloid cellsScienceYear: 198623261652420009  
Begley CG,Nicola NA,Metcalf D,Proliferation of normal human promyelocytes and myelocytes after a single pulse stimulation by purified GMCSF or GCSFBloodYear: 1988716406452449922  
Lord BI,Bronchud MH,Owens S,Chang J,Howell A,Souza L,Dexter TM,The kinetics of human granulopoiesis following treatment with granulocyte colonystimulating factor in vivoProc Natl Acad Sci USAYear: 198986949995032480603  
Mackey MC,Aprikyan AA,Dale DC,The rate of apoptosis in post mitotic neutrophil precursors of normal and neutropenic humansCell ProlifYear: 200336273412558658  
Kim HK,De La Luz Sierra M,Williams CK,Gulino AV,Tosato G,GCSF downregulation of CXCR4 expression identified as a mechanism for mobilization of myeloid cellsBloodYear: 200610881282016537807  
Christopher MJ,Link DC,Regulation of neutrophil homeostasisCurr Opin HematolYear: 2007143817133093  
Diehl V,Franklin J,Pfreundschuh M,Lathan B,Paulus U,Hasenclever D,Tesch H,Herrmann R,Dorken B,MullerHermelink HK,Duhmke E,Loeffler M,Standard and increaseddose BEACOPP chemotherapy compared with COPPABVD for advanced Hodgkin’s diseaseN Engl J MedYear: 20033482386239512802024  
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 HighGrade NonHodgkin’s Lymphoma Study GroupTwoweekly or 3weekly chop chemotherapy with or without etoposide for the treatment of elderly patients with aggressive lymphomas: results of the NHLB2 trial of the DSHNHLBloodYear: 2004104363464115016643  
Engel C,Loeffler M,Schmitz S,Tesch H,Diehl V,Acute hematologic toxicity and practicability of doseintensified BEACOPP chemotherapy for advanced stage Hodgkin’s disease. German Hodgkin’s Lymphoma Study Group (GHSG)Ann OncolYear: 200011911051411061603  
Frampton JE,Lee CR,Faulds D,Filgrastim. A review of its pharmacological properties and therapeutic efficacy in neutropeniaDrugsYear: 1994485731607530630  
Frampton JE,Yarker YE,Goa KL,Lenograstim. A review of its pharmacological properties and therapeutic efficacy in neutropenia and related clinical settingsDrugsYear: 1995495767937541335  
Houston AC,Stevens LA,Cour V,Pharmacokinetics of glycosylated recombinant human granulocyte colonystimulating factor (lenograstim) in healthy male volunteersBr J Clin PharmacolYear: 19994727928410215752  
Tanaka H,Tokiwa T,Influence of renal and hepatic failure on the pharmacokinetics of recombinant human granulocyte colonystimulating factor (KRN8601) in the ratCancer ResYear: 199050661566191698539  
Khwaja A,Carver J,Jones HM,Paterson D,Linch DC,Expression and dynamic modulation of the human granulocyte colonystimulating factor receptor in immature and differentiated myeloid cellsBr J HaematolYear: 1993852542597506564  
Ericson SG,Gao H,Gericke GH,Lewis LD,The role of polymorphonuclear neutrophils (PMNs) in clearance of granulocyte colonystimulating factor (GCSF) in vivo and in vitroExp HematolYear: 199725131313259406990  
El Ouriaghli F,Fujiwara H,Melenhorst JJ,Sconocchia G,Hensel N,Barrett AJ,Neutrophil elastase enzymatically antagonizes the in vitro action of GCSF: implications for the regulation of granulopoiesisBloodYear: 20031011752175812393522  
Hunter MG,Druhan LJ,Massullo PR,Avalos BR,Proteolytic cleavage of granulocyte colonystimulating factor and its receptor by neutrophil elastase induces growth inhibition and decreased cell surface expression of the granulocyte colonystimulating factor receptorAm J HematolYear: 20037414915514587040  
KottoKome AC,Fox SE,Lu W,Yang BB,Christensen RD,Calhoun DA,Evidence that the granulocyte colonystimulating factor (GCSF) receptor plays a role in the pharmacokinetics of GCSF and PegGCSF using a GCSFR KO modelPharmacol ResYear: 200450555815082029  
Tanaka H,SatakeIshikawa R,Ishikawa M,Matsuki S,Asano K,Pharmacokinetics of recombinant human granulocyte colonystimulating factor conjugated to polyethylene glycol in ratsCancer ResYear: 199151371037141712246  
Crawford J,Clinical uses of pegylated pharmaceuticals in oncologyCancer Treat RevYear: 200228Suppl A71112173409  
Yang BB,Lum PK,Hayashi MM,Roskos LK,Polyethylene glycol modification of filgrastim results in decreased renal clearance of the protein in ratsJ Pharm SciYear: 2004931367137315067712  
Harris JM,Chess RB,Effect of pegylation on pharmaceuticalsNat Rev Drug DiscovYear: 2003221422112612647  
Lowenhaupt K,Wang PJ,Horan T,Lauffenburger DA,Sarkar CA,Parsing the effects of binding, signaling, and trafficking on the mitogenic potencies of granulocyte colonystimulating factor analoguesBiotechnol ProgYear: 20031995596412790662  
Veronese FM,Mero A,The impact of PEGylation on biological therapiesBioDrugsYear: 20082231532918778113  
Scholz M,Engel C,Apt D,Sankar SL,Goldstein E,Loeffler M,Pharmacokinetic and Pharmacodynamic modelling of the novel human GCSF derivative MaxyG34 and Pegfilgrastim in the ratCell ProlifYear: 20094268233719780759  
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 oncepercycle pegfilgrastim and daily injection filgrastim in chemotherapyinduced neutropenia: a multicenter dosefinding study in women with breast cancerAnn OncolYear: 20021390390912123336  
Grigg A,SolalCeligny 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 nonHodgkin’s lymphomaLeuk LymphomaYear: 2003441503150814565651  
Vose JM,Crump M,Lazarus H,Emmanouilides C,Schenkein D,Moore J,Frankel S,Flinn I,Lovelace W,Hackett J,Liang BC,Randomized, Multicenter, OpenLabel Study of Pegfilgrastim Compared With Daily Filgrastim After Chemotherapy for LymphomaJ Clin OncolYear: 200321351451912560443  
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 metaanalysis of randomized controlled trialsCurr Med Res OpinYear: 2007232283229517697451  
Ziepert M,Schmits R,Trumper L,Pfreundschuh M,Loeffler M,Prognostic factors for hematotoxicity of chemotherapy in aggressive nonHodgkin’s lymphomaAnn OncolYear: 20081975276218048382  
Kroschinsky F,Holig K,Ehninger G,The role of pegfilgrastim in mobilization of hematopoietic stem cellsTransfus Apher SciYear: 20083832374418490197  
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 cyclophosphamideinduced granulocytopaeniaCell ProlifYear: 20094268132219689472  
Hartmann F,Zeynalova S,Nickenig C,Reiser M,Lengfelder E,Duerk H,de Witt M,Schubert J,Loeffler M,Pfreundschuh M,Pegfilgrastim (PegF) on day 4 of (R)CHOP14 chemotherapy compared to day 2 in elderly patients with diffuse large Bcell lymphoma (DLBCL): Results of a randomized trial of the German highgrade nonHodgkin’s lymphoma study group (DSHNHL)J Clin Oncol, ASCO Ann Meeting Proc Part IYear: 20072518S19511  
Scholz M,Ackermann M,Emmrich F,Loeffler M,Kamprad M,Effectiveness of cytopenia prophylaxis for different filgrastim and pegfilgrastim schedules in a chemotherapy mouse modelBiologics: Targets & TherapyYear: 200932737  
Engel C,Scholz M,Loeffler M,A computational model of human granulopoiesis to simulate the hematotoxic effects of multicycle polychemotherapyBloodYear: 20041042323233115226178  
Scholz M,Engel C,Loeffler M,Modelling human granulopoiesis under polychemotherapy with GCSF supportJ Math BiolYear: 20055039743915614553  
Scholz M,Engel C,Loeffler M,Modelbased design of chemotherapeutic regimens that account for heterogeneity in leucopoeniaBr J HaematolYear: 200613272373516487172  
Tsuji K,Ebihara Y,Expression of GCSF receptor on myeloid progenitorsLeuk LymphomaYear: 20014261351711911419  
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 functionJ Clin InvestYear: 1998101112456679616217  
Kuwabara T,Uchimura T,Takai K,Kobayashi H,Kobayashi S,Sugiyama Y,Saturable uptake of a recombinant human granulocyte colonystimulating factor derivative, nartograstim, by the bone marrow and spleen of rats in vivoJ Pharmacol Exp TherYear: 1995273111411227540687  
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 modelDrug Metab DisposYear: 2007352211221717875672  
Kuwabara T,Kobayashi S,Sugiyama Y,Pharmacokinetics and pharmacodynamics of a recombinant human granulocyte colony stimulating factorDrug Metab RevYear: 1996286256588959393  
Lenhoff S,Rosberg B,Olofsson T,Granulocyte interactions with GMCSF and GCSF secretion by endothelial cells and monocytesEur Cytokine NetworkYear: 1999104525532  
Layton JE,Hockman H,Sheridan WP,Morstyn G,Evidence for a novel in vivo control mechanism of granulopoiesis: mature cellrelated control of a regulatory growth factorBloodYear: 198974130313072475185  
Sarkar CA,Lauffenburger DA,Celllevel pharmacokinetic model of granulocyte colonystimulating factor: implications for ligand lifetime and potency in vivoMol PharmacolYear: 20036314715812488547  
Roskos LK,Lum P,Lockbaum P,Schwab G,Yang BB,Pharmacokinetic / pharmacodynamic modeling of pegfilgrastim in healthy subjectsJ Clin PharmacolYear: 20064674775716809800  
Wang B,Ludden TM,Cheung EN,Schwab GG,Roskos LK,Population pharmacokineticpharmacodynamic modeling of filgrastim (rmetHuGCSF) in healthy volunteersJ Pharmacokinet PharmacodynYear: 20012832134211677930  
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 colonystimulating factorBloodYear: 1999932506251410194429  
Levesque JP,Takamatsu Y,Nilsson SK,Haylock DN,Simmons PJ,Vascular cell adhesion molecule1 (CD106) is cleaved by neutrophil proteases in the bone marrow following hematopoietic progenitor cell mobilization by granulocyte colonystimulating factorBloodYear: 2001981289129711520773  
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 colonystimulating factor after autologous blood stem cell transplantationExp HematolYear: 199523149715028542937  
Steinman RA,Tweardy DJ,Granulocyte colonystimulating factor receptor mRNA upregulation is an immediate early marker of myeloid differentiation and exhibits dysfunctional regulation in leukemic cellsBloodYear: 1994831191278274731  
Tkatch LS,Rubin KA,Ziegler SF,Tweardy DJ,Modulation of human GCSF receptor mRNA and protein in normal and leukemic myeloid cells by GCSF and retinoic acidJ Leukoc BiolYear: 1995579649717540644  
Srinivasa SP,Doshi PD,Extracellular signalregulated kinase and p38 mitogenactivated protein kinase pathways cooperate in mediating cytokineinduced proliferation of a leukemic cell lineLeukemiaYear: 20021622445311840291  
Lohrmann HP,Schreml W,Cytotoxic Drugs and the Granulopoietic SystemYear: 1982Berlin: Springer Verlag  
Scholz M,Gross A,Loeffler M,A biomathematical model of human thrombopoiesis under chemotherapyJ Theor BiolYear: 2010264228730020083124  
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 gcsf mutant (Ro 258315) in healthy volunteers: comparison with single and multiple daily doses of filgrastimAm J HematolYear: 200166424525111279634  
Borleffs JC,Bosschaert M,Vrehen HM,Schneider MM,van Strijp J,Small MK,Borkett KM,Effect of escalating doses of recombinant human granulocyte colonystimulating factor (filgrastim) on circulating neutrophils in healthy subjectsClin TherYear: 19982047227369737832  
Chatta GS,Price TH,Allen RC,Dale DC,Effects of in vivo recombinant methionyl human granulocyte colonystimulating factor on the neutrophil response and peripheral blood colonyforming cells in healthy young and elderly adult volunteersBloodYear: 1994849292329297524759  
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 (AVI014) in healthy subjects: a firstin human, single dose, controlled studyBMC Clin PharmacolYear: 200992  
Johnston E,Crawford J,Blackwell S,Bjurstrom T,Lockbaum P,Roskos L,Yang BB,Gardner S,MillerMessana MA,Shoemaker D,Garst J,Schwab G,Randomized, doseescalation study of sd/01 compared with daily filgrastim in patients receiving chemotherapyJ Clin OncolYear: 200018132522252810893282  
George S,Yunus F,Case D,Yang BB,Hackett J,Shogan JE,Meza LA,Neumann TA,Liang BC,Fixeddose pegfilgrastim is safe and allows neutrophil recovery in patients with nonhodgkin’s lymphomaLeuk LymphomaYear: 200344101691169614692520  
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 HighGrade NonHodgkin’s Lymphoma Study GroupSix versus eight cycles of biweekly CHOP14 with or without rituximab in elderly patients with aggressive CD20+ Bcell lymphomas: a randomised controlled trial (RICOVER60)Lancet OncolYear: 20089210511618226581  
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 HighGrade NonHodgkin Lymphoma Study GroupRandomized comparison of pegfilgrastim day 4 versus day 2 for the prevention of chemotherapyinduced leukocytopeniaAnn OncolYear: 20112281872187721292644  
Rechenberg I,Evolutionsstrategie ’94Year: 1994Stuttgart: frommannholzboog  
Schwefel HP,Evolution strategies: A family of nonlinear optimization techniques based on imitating some principles of organic evolutionAnn Oper ResYear: 1984165167  
Vainstein V,Ginosar Y,Shoham M,Ranmar DO,Ianovski A,Agur Z,The complex effect of granulocyte colonystimulating factor on human granulopoiesis analyzed by a new physiologicallybased mathematical modelJ Theor BiolYear: 2005234331132715784267  
Foley C,Bernard S,Mackey MC,Costeffective GCSF therapy strategies for cyclical neutropenia: mathematical modelling based hypothesesJ Theor BiolYear: 2006238475476316115650  
Ostby I,Rusten LS,Kvalheim G,Grottum P,A mathematical model for reconstitution of granulopoiesis after high dose chemotherapy with autologous stem cell transplantationJ Math BiolYear: 200347210113612883857  
Ostby I,Kvalheim G,Rusten LS,Grottum P,Mathematical modeling of granulocyte reconstitution after highdose chemotherapy with stem cell support: effect of posttransplant GCSF treatmentJ Theor BiolYear: 20042311698315363930  
Shochat E,RomKedar V,Segel LA,GCSF control of neutrophils dynamics in the bloodBull Math BiolYear: 2007697229933817554586  
Foley C,Mackey MC,Mathematical model for GCSF administration after chemotherapyJ Theor BiolYear: 20092571274419007795  
Waladkhani AR,Pegfilgrastim: a recent advance in the prophylaxis of chemotherapyinduced neutropeniaEur J Cancer CareYear: 200413371379  
Wiczling P,Lowe P,Pigeolet E,Luedicke F,Balser S,Krzyzanski W,Population pharmacokinetic modelling of filgrastim in healthy adults following intravenous and subcutaneous administrationsClin PharmacokinetYear: 200948128172619902989  
Wichmann HE,Loeffler M,Mathematical Modeling of Cell Proliferation: Stem Cell Regulation in HemopoiesisYear: 1985Boca Raton: sCRC Press  
Schmitz S,Franke H,Loeffler M,Wichmann HE,Diehl V,Model analysis of the contrasting effects of GMCSF and GCSF treatment on peripheral blood neutrophils observed in three patients with childhoodonset cyclic neutropenia BritJ HaematolYear: 199695616625  
Dale DC,Fauci AS,Wolff SM,Alternateday prednisoneTnE J MedYear: 19742912211541158 
Figures
[Figure ID: F1] 
Figure 1
(Basic structure of the cellkinetic model of granulopoiesis). Major model compartments describing granulopoietic cell stages are S (pluripotent stem cells), CG (colony forming units of granulocytes and macrophages), PGB (proliferating granulopoietic blasts), MGB (maturing granulopoietic blasts  subdivided into metamyelocytes (G4),banded granulocytes (G5) and segmented granulocytes (G6)) and GRA (circulating granulocytes). The system is regulated by feedback loops. A major loop is mediated by GCSF which is produced endogenously but can also be applied subcutaneously. Chemotherapy (CX) induces acute cell loss. The model is essentially the same as in [^{40}]. 
[Figure ID: F2] 
Figure 2
(Model structure of the pharmacokinetic model of GCSF). The major compartments, cytokine fluxes and regulations are presented (MM = Michaelis Menten kinetic). The subcutaneous compartment is divided into two subcompartments with first order transition. 
[Figure ID: F3] 
Figure 3
(Estimated bioavailability of subcutaneously injected Filgrastim or Pegfilgrastim based on systematic model simulations): Bioavailability was estimated by calculating GCSF amounts absorbed by the central compartment relative and the total amount of subcutaneously injected GCSF (xaxis). Due to the modelled loss in the subcutaneous tissue, the bioavailability is dosedependent. Circles indicate estimates for pharmaceutically available doses of 300μg and 480μg of Filgrastim or 6000 μg of Pegfilgrastim respectively [^{12},^{67}]. 
[Figure ID: F4] 
Figure 4
(Comparison of model and data for GCSF applications): Comparison of model and selected datasets of single Filgrastim injections (scenario A), multiple Filgrastim injections (scenario B) and single Pegfilgrastim injections of different doses (scenarios C and D). For each scenario, we present the time courses of ANC and GCSF, respectively. A complete list of scenarios can be found in the Additional file 1. 
[Figure ID: F5] 
Figure 5
(Regulation functions of Filgrastim and Pegfilgrastim): Comparison of Filgrastim and Pegfilgrastim with respect to the regulation function of the amplification in PGB. The circle marks the value under steadystate conditions. 
[Figure ID: F6] 
Figure 6
(Effect of GCSF delay): Effect of the delay parameter of GCSF action on cellcounts of specific cell compartments. 
[Figure ID: F7] 
Figure 7
(Effect of chemotherapy delay): Effect of the delay parameter of the chemotherapy studied for the CHOP regimen. 
[Figure ID: F8] 
Figure 8
(Comparison of model and data for chemotherapy scenarios): Comparison of model and data for the CHOP21 regimen and time intensified CHOP14 regimen supported by various Filgrastim and Pegfilgrastim schedules. The solid line is the model prediction. Dots are patient medians at corresponding time points and the grey lines mark the interquartile range of the data. All scenarios are based on the same model parameters. 
[Figure ID: F9] 
Figure 9
(Ratio of granulocytes and leukocytes): Based on model simulations, the ratio of granulocytes and leukocytes under CHOP14 chemotherapy is predicted. 
[Figure ID: F10] 
Figure 10
(Validation of model): Validation of the model on the basis of two datasets not used for model fitting. Solid line is the model prediction. Dots and dotted lines are the data and the interpolated data respectively. 
Tables
Major compartments of the pharmacokinetic model and corresponding regulations
Compartment  Regulations  Regulator 

S

proliferative fraction

bone marrow content


selfrenewal probability

bone marrow content

CG

proliferative fraction

bone marrow content


amplification

GCSF


transition time

GCSF

PGB

amplification

GCSF


transition time

GCSF

MGB

postmitotic amplification

GCSF


transition time

GCSF

GRA

turnover



GCSF

endogenous production

late bone marrow cell stages


specific degradation

GRA

external applications   
Most of these compartments are mediated by GCSF. A complete set of model equations can be found in the Additional file 1.
Data sets utilized for the establishment and validation of the granulopoiesis model
Type of data  Disease  GCSF schedules  Chemotherapy  References 

Phase I studies with Filgrastim

none

single application 3, 5, 10 μg/kg

none

[^{31}]

none

single application 5, 10 μg/kg

none

[^{60}]



none

10 applications 75, 150, 300, 600 μg

none

[^{61}]


none

14 applications 30, 300 μg

none

[^{62}]


none

single application 4, 8 μg/kg

none

[^{63}]

Phase I studies with Pegfilgrastim

none

30, 60, 100, 300 μg/kg

none

[^{50}]

Phase II studies with Pegfilgrastim

LuCa

30, 100, 300 μg/kg

none

[^{64}]


NHL

6000 μg, day 2

CHOP

[^{65}]

Phase III studies with CX and w/wo Filgrastim

NHL

no GCSF

CHOP21*

[^{12}]

NHL

480 μg, day 4–13

CHOP14*

[^{12}]


NHL

480 μg, day 6–12

CHOP14*

[^{66}]


Phase III studies with CX + Peg  NHL  6000 μg, day 2, 4  CHOP14*  [^{67}] 
Studies with access to raw data are indicated with an asterisk. CX = chemotherapy, LuCa = lung cancer, NHL = highgrade nonHodgkin’s lymphoma.
Pharmacokinetic parameters of Filgrastim and Pegfilgrastim
Parameter  Meaning  Filgrastim  Pegfilgrastim 

k_{sc}

subcutaneous absorption [h^{−1}]

0.161

0.107

k_{m}

MichaelisMenten constant of

34.7

5.5


subcutaneous elimination [μg]



v_{max}

Maximum of subcutaneous

67.3

16.5


elimination [h^{−1}]



k_{u}

unspecific elimination [h^{−1}]

0.441

0.087

kmGRA

MichaelisMenten constant of

22.4

30.8


specific elimination [μg]



vmaxGRA

Maximum of specific elimination [h^{−1}]

4.77

5.16

k_{cp}

transition central to peripheral [h^{−1}]

0.000

0.075

k_{pc}

transition peripheral to central [h^{−1}]



0.548

V_{D}

distribution volume [l]

1.156

4.091

CG−CSFcent_ref  reference GCSF serum concentration μgl  0.02 
Compared to Pegfilgrastim, we estimated that Filgrastim is more easily absorbed from the subcutaneous compartment, has a lower bioavailability (see Figure 3) and a higher specific and unspecific elimination. Reversible binding is neglectable for Filgrastim but not for Pegfilgrastim. The distribution volume is higher for Pegfilgrastim than for Filgrastim.
Article Categories:
Keywords: Chemotherapy, Filgrastim, Granulopoiesis, Haematotoxicity, Leucopenia, Pegfilgrastim. 
Previous Document: CYP1A1 MspI polymorphism and acute myeloid leukemia risk: metaanalyses based on 5018 subjects.
Next Document: Superimposed Preeclampsia in Women with Chronic Kidney Disease.