Indirect-response model for the time course of leukopenia with anticancer drugs

Supramolecular Biomaterials for Cancer Immunotherapy

Cancer immunotherapy has achieved tremendous successful clinical results and obtained historic victories in tumor treatments. However, great limitations associated with feeble immune responses and serious adverse effects still cannot be neglected due to the complicated multifactorial etiology and pathologic microenvironment in tumors. The rapid development of nanomedical science and material science has facilitated the advanced progress of engineering biomaterials to tackle critical issues. The supramolecular biomaterials with flexible and modular structures have exhibited unparalleled advantages of high cargo-loading efficiency, excellent biocompatibility, and diversiform immunomodulatory activity, thereby providing a powerful weapon for cancer immunotherapy. In past decades, supramolecular biomaterials were extensively explored as versatile delivery platforms for immunotherapeutic agents or designed to interact with the key moleculars in immune system in a precise and controllable manner. In this review, we focused on the crucial role of supramolecular biomaterials in the modulation of pivotal steps during tumor immunotherapy, including antigen delivery and presentation, T lymphocyte activation, tumor-associated macrophage elimination and repolarization, and myeloid-derived suppressor cell depletion. Based on extensive research, we explored the current limitations and development prospects of supramolecular biomaterials in cancer immunotherapy.

Pharmacokinetic-pharmacodynamic modelling of the hypoglycaemic effect of pulsatile administration of human insulin in rats

The relationship between the plasma insulin (INS) concentration–time course and plasma glucose concentration–time course during and after pulsatile INS administration to rats was characterized using a pharmacokinetic–pharmacodynamic (PK–PD) model. A total INS dose of 0.5 IU/kg was intravenously injected in 2 to 20 pulses over a 2-h period. Compared with the single bolus administration, the area under the effect-time curve (AUE) increased depending on the number of pulses, and the AUEs for more than four pulses plateaued at a significantly larger value, which was similar to that after the infusion of a total of 0.5 IU/kg of INS over 2 h. No increase in plasma INS concentration occurred after pulsatile administration. Two indirect response models primarily reflecting the receptor-binding process (IR model) or glucose transporter 4 (GLUT4) translocation (GT model) were applied to describe the PK–PD relationship after single intravenous bolus administration of INS. These models could not explain the observed data after pulsatile administration. However, the IR-GT model, which was a combination of the IR and GT models, successfully explained the effects of pulsatile administration and intravenous infusion. These results indicate that the receptor-binding process and GLUT4 translocation are responsible for the change in AUE after pulsatile administration.

Recent advances in physiologically based pharmacokinetic and pharmacodynamic models for anticancer nanomedicines

Nanoparticles (NPs) have distinct pharmacokinetic (PK) properties and can potentially improve the absorption, distribution, metabolism, and elimination (ADME) of small-molecule drugs loaded therein. Owing to the unwanted toxicities of anticancer agents in healthy organs and tissues, their precise delivery to the tumor is an essential requirement. There have been numerous advancements in the development of nanomedicines for cancer therapy. Physiologically based PK (PBPK) models serve as excellent tools for describing and predicting the ADME properties and the efficacy and toxicity of drugs, in combination with pharmacodynamic (PD) models. The recent preliminary application of these modeling approaches to NPs demonstrated their potential benefits in research and development processes relevant to the ADME and pharmacodynamics of NPs and nanomedicines. Here, we comprehensively review the pharmacokinetics of NPs, the developed PBPK models for anticancer NPs, and the developed PD model for anticancer agents.

Mathematical models for cytarabine-derived myelosuppression in acute myeloid leukaemia

We investigate the personalisation and prediction accuracy of mathematical models for white blood cell (WBC) count dynamics during consolidation treatment using intermediate or high-dose cytarabine (Ara-C) in acute myeloid leukaemia (AML). Ara-C is the clinically most relevant cytotoxic agent for AML treatment. We extend a mathematical model of myelosuppression and a pharmacokinetic model of Ara-C with different hypotheses of Ara-C's pharmacodynamic effects. We cross-validate the 12 model variations using dense WBC count measurements from 23 AML patients. Surprisingly, the prediction accuracy remains satisfactory in each of the models despite different modelling hypotheses. Therefore, we compare average clinical and calculated WBC recovery times for different Ara-C schedules as a successful methodology for model discrimination. As a result, a new hypothesis of a secondary pharmacodynamic effect on the proliferation rate seems plausible. Furthermore, we demonstrate the impact of treatment timing on subsequent nadir values based on personalised predictions as a possibility for influencing/controlling myelosuppression.

Understanding Hematological Toxicities Using Mathematical Modeling

Balancing antitumor efficacy with toxicity is a significant challenge, and drug-induced myelosuppression is a common dose-limiting toxicity of cancer treatments. Mathematical modeling has proven to be a powerful ally in this field, scaling results from animal models to humans, and designing optimized treatment regimens. Here we outline existing mathematical approaches for studying bone marrow toxicity, identify gaps in current understanding, and make future recommendations to advance this vital field of safety research further.

Anticancer and cardio-protective effects of liposomal doxorubicin in the treatment of breast cancer

Breast cancer (BC) is a highly prevalent disease, accounting for the second highest number of cancer-related mortalities worldwide. The anthracycline doxorubicin (DOX), isolated from Streptomyces peucetius var. caesius, is a potent chemotherapeutic drug that is successfully used to treat various forms of liquid and solid tumors and is currently approved to treat BC. DOX exerts its effects by intercalation into DNA and inhibition of topoisomerases I and II, causing damage to DNA and the formation of reactive oxygen species (ROS), resulting in the activation of caspases, which ultimately leads to apoptosis. Unfortunately, DOX also can cause cardiotoxicity, with patients only allowed a cumulative lifetime dose of 550 mg/m². Efforts to decrease cardiotoxicity and to increase the blood circulation time of DOX led to the US Food and Drug Administration (FDA) approval of a PEGylated liposomal formulation (L-DOX), Doxil^® (known internationally as Caelyx^®). Both exhibit better cardiovascular safety profiles; however, they are not currently FDA approved for the treatment of metastatic BC. Here, we provide detailed insights into the mechanism of action of L-DOX and its most common side effects and highlight results of its use in clinical trials for the treatment of BC as single agent and in combination with other commonly used chemotherapeutics.

Administration of temozolomide: Comparison of conventional and metronomic chemotherapy regimens

PURPOSE

We compare the Maximum Tolerated Dose (MTD) and Metronomic Chemotherapy (MC) protocols for temozolomide administration. We develop an innovative methodology for characterizing optimal chemotherapy regimens.

METHODS

We use a PK/PD model based on Faivre et al. (2013) for the pharmacokinetics of temozolomide, as well as the pharmacodynamics of its efficacy. For toxicity, which is measured by the nadir of the normalized absolute neutrophil count, we formalize the myelosuppression effect of temozolomide with the physiological model of Panetta et al. (2003b). We introduce a multi-criteria tool for comparing protocols along their efficacy and toxicity dimensions.

RESULTS

We show that the toxicity of the MC regimen proposed by Faivre et al. (2013) can greatly be reduced without affecting its efficacy, while the standard MTD protocol efficacy cannot be improved without impairing its toxicity. We also show that for any acceptable toxicity level, the optimal protocol remains closely related to standard MTD.

CONCLUSIONS

Overall, our new method enables a rich comparison between protocols along multiple dimensions. We can rank protocols for temozolomide administration. It is a first step toward building optimal individual protocols.

Vaccination during myeloid cell depletion by cancer chemotherapy fosters robust T cell responses

Therapeutic vaccination with human papillomavirus type 16 synthetic long peptides (HPV16-SLPs) results in T cell-mediated regression of HPV16-induced premalignant lesions but fails to install clinically effective immunity in patients with HPV16-positive cervical cancer. We explored whether HPV16-SLP vaccination can be combined with standard carboplatin and paclitaxel chemotherapy to improve immunity and which time point would be optimal for vaccination. This was studied in the HPV16 E6/E7-positive TC-1 mouse tumor model and in patients with advanced cervical cancer. In mice and patients, the presence of a progressing tumor was associated with abnormal frequencies of circulating myeloid cells. Treatment of TC-1-bearing mice with chemotherapy and therapeutic vaccination resulted in superior survival and was directly related to a chemotherapy-mediated altered composition of the myeloid cell population in the blood and tumor. Chemotherapy had no effect on tumor-specific T cell responses. In advanced cervical cancer patients, carboplatin-paclitaxel also normalized the abnormal numbers of circulating myeloid cells, and this was associated with increased T cell reactivity to recall antigens. The effect was most pronounced starting 2 weeks after the second cycle of chemotherapy, providing an optimal immunological window for vaccination. This was validated with a single dose of HPV16-SLP vaccine given in this time window. The resulting proliferative HPV16-specific T cell responses were unusually strong and were retained after all cycles of chemotherapy. In conclusion, carboplatin-paclitaxel therapy fosters vigorous vaccine-induced T cell responses when vaccination is given after chemotherapy and has reset the tumor-induced abnormal myeloid cell composition to normal values.

Therapeutic drug monitoring of cancer chemotherapy

Therapeutic drug monitoring is not routinely used for chemotherapy agents. There are Several reasons, but one major drawback is the lack of established therapeutic Concentration ranges. Combination chemotherapy makes the establishment of Therapeutic ranges for individual drugs difficult, the concentration-effect relationship for a single drug may not be the same as when that drug is used in a drug combination. Pharmacokinetic optimization protocols for many classes of cytotoxic compounds exist in specialized centers, and some of these protocols are now part of large multicentre trials. Nonetheless, TDM clearly has the potential to improve the clinical use of chemotherapy gents, most of which have very narrow therapeutic indices and highly variable pharmacokinetics. A substantial body of literature accumulating during the past 15 years demonstrates relationships between systemic exposure to various chemotherapy agents and their toxic or therapeutic effects. This article reviews TDM concepts in addition to tools based on pharmacokinetic modeling of chemotherapy agents. The administered dose of chemotherapy agents is sometimes adjusted individually using either a priori or a posteriori methods. These models can only be applied by using the same dose and schedule as the original study. Bayesian estimation offers more flexibility in blood sampling times and, owing to its precision and to the amount of information provided is the method of choice for ensuring that a given patient benefits from the desired systemic exposure. Moreover, the role and application of Pharmacogenetics as a tool for individualizing chemotherapy is discussed highlighting the agents and mechanisms that have been well studied and defined and their relevance to clinical practice. Finally, this paper address issues critical to the optimal use of TDM in a clinical setting, and the role of clinical pharmacist in this regard. In addition, it discusses future developments in this field that can contribute to improving cancer chemotherapy In terms of patient outcome and survival. J Oncol Pharm Practice (2007) 13: 207—221.

Pharmacodynamic modeling of adverse effects of anti-cancer drug treatment

Purpose

Adverse effects related to anti-cancer drug treatment influence patient’s quality of life, have an impact on the realized dosing regimen, and can hamper response to treatment. Quantitative models that relate drug exposure to the dynamics of adverse effects have been developed and proven to be very instrumental to optimize dosing schedules. The aims of this review were (i) to provide a perspective of how adverse effects of anti-cancer drugs are modeled and (ii) to report several model structures of adverse effect models that describe relationships between drug concentrations and toxicities.

Methods

Various quantitative pharmacodynamic models that model adverse effects of anti-cancer drug treatment were reviewed.

Results

Quantitative models describing relationships between drug exposure and myelosuppression, cardiotoxicity, and graded adverse effects like fatigue, hand-foot syndrome (HFS), rash, and diarrhea have been presented for different anti-cancer agents, including their clinical applicability.

Conclusions

Mathematical modeling of adverse effects proved to be a helpful tool to improve clinical management and support decision-making (especially in establishment of the optimal dosing regimen) in drug development. The reported models can be used as templates for modeling a variety of anti-cancer-induced adverse effects to further optimize therapy.

The effect of "Jelly" CdTe QD uptake on RAW264.7 monocytes: immune responses and cell fate study

Encapsulation of Quantum Dots (QDs) has become an essential factor which regulates particles cytotoxicity, as well as physical and chemical stability. Negatively charged cellular membranes have a great affinity to nanoparticles with surface molecules carrying positive charge, hence creating perfect conditions for fast and aggressive intracellular penetration. The preference for non-charged outer shells is topical in QD design and various applications. In the current paper we develop gelatination as a prominent coating approach to create neutrally passivated QDs with improved biocompatibility. We have revealed the trends in particle's uptake, accumulation, intracellular localisation and retaining time as well as RAW264.7 monocyte cell fate and immune responses. Also the difference in particle endocytosis kinetics and dynamics has been shown to depend on the QD core size. The intracellular QD content along with cell responses at the population level was quantified by flow cytometry.

Pharmacodynamic models of age-structured cell populations

The purpose of this work is to review basic pharmacodynamic (PD) models describing drug effects on cell populations and expand them to age-structured models using the theory of physiologically structured populations. The plasma drug concentrations are interpreted as the environment affecting the cell production and mortality rates. An explicit solution to model equations provides the age density distribution that serves to establish a relationship between the cell lifespan distribution and the hazard of cell removal. Given the lifespan distributions, the age distributions for most commonly applied PD models of cell responses including basic cell turnover, transit compartments, and basic lifespan models have been derived both for the baseline conditions and drug treatment. The steady-state age distribution for basic indirect response models is exponential, and it is uniform for the basic lifespan model. As an example of more complex cell population, the age distribution of human red blood cells has been simulated based on a recent model of red blood cell survival. The age distribution for cells in the transit compartment model is the sum of the gamma functions. Means and variances of age distributions for all discussed models were calculated. A brief discussion of numerical challenges and possible future model developments is presented.

Semimechanistic cell-cycle type-based pharmacokinetic/pharmacodynamic model of chemotherapy-induced neutropenic effects of diflomotecan under different dosing schedules

The current work integrates cell-cycle dynamics occurring in the bone marrow compartment as a key element in the structure of a semimechanistic pharmacokinetic/pharmacodynamic model for neutropenic effects, aiming to describe, with the same set of system- and drug-related parameters, longitudinal data of neutropenia gathered after the administration of the anticancer drug diflomotecan (9,10-difluoro-homocamptothecin) under different dosing schedules to patients (n = 111) with advanced solid tumors. To achieve such an objective, the general framework of the neutropenia models was expanded, including one additional physiologic process resembling cell cycle dynamics. The main assumptions of the proposed model are as follows: within the stem cell compartment, proliferative and quiescent cells coexist, and only cells in the proliferative condition are sensitive to drug effects and capable of following the maturation chain. Cell cycle dynamics were characterized by two new parameters, FProl (the fraction of proliferative [Prol] cells that enters into the maturation chain) and kcycle (first-order rate constant governing cell cycle dynamics within the stem cell compartment). Both model parameters were identifiable as indicated by the results from a bootstrap analysis, and their estimates were supported by date from the literature. The estimates of FProl and kcycle were 0.58 and 1.94 day(-1), respectively. The new model could properly describe the neutropenic effects of diflomotecan after very different dosing scenarios, and can be used to explore the potential impact of dosing schedule dependencies on neutropenia prediction.

Mechanism-based model characterizing bidirectional interaction between PEGylated liposomal CKD-602 (S-CKD602) and monocytes in cancer patients

S-CKD602 is a PEGylated liposomal formulation of CKD-602, a potent topoisomerase I inhibitor. The objective of this study was to characterize the bidirectional pharmacokinetic-pharmacodynamic (PK-PD) interaction between S-CKD602 and monocytes. Plasma concentrations of encapsulated CKD-602 and monocytes counts from 45 patients with solid tumors were collected following intravenous administration of S-CKD602 in the phase I study. The PK-PD models were developed and fit simultaneously to the PK-PD data, using NONMEM(®). The monocytopenia after administration of S-CKD602 was described by direct toxicity to monocytes in a mechanism-based model, and by direct toxicity to progenitor cells in bone marrow in a myelosuppression-based model. The nonlinear PK disposition of S-CKD602 was described by linear degradation and irreversible binding to monocytes in the mechanism-based model, and Michaelis-Menten kinetics in the myelosuppression-based model. The mechanism-based PK-PD model characterized the nonlinear PK disposition, and the bidirectional PK-PD interaction between S-CKD602 and monocytes.

Comparison of different semi-mechanistic models for chemotherapy-related neutropenia: application to BI 2536 a Plk-1 inhibitor

PURPOSE

The aim of this investigation was to compare the performance of a commonly used semi-mechanistic model for drug-related neutropenia with other semi-mechanistic models published in the literature.

METHODS

After their implementation in NONMEM VI, five semi-mechanistic models were assessed using the pharmacokinetic and absolute neutrophil count data obtained from 95 patients with non-small cell lung cancer receiving either 200 mg on day 1 or 50 or 60 mg on days 1, 2 and 3 of a 21-day treatment course with the new Plk-1 inhibitor BI 2536. The model performance was compared by means of predictive (visual and numerical) checks, precision in the parameter estimates and objective function-based measures. Details of model parameterization, model stability and run times are also provided.

RESULTS

The time course of the drug plasma concentrations was described by a three compartment model with a first-order elimination rate. With respect to neutropenia, all models were successfully implemented in NONMEM and provided reasonable fits for the median (although not all models described all percentiles of the data well), and in general precise parameter estimates.

CONCLUSION

In the current evaluation performed in a single drug, none of the models showed superior performance compared to the most commonly used model first described by Friberg et al. (J Clin Oncol 20:4713-4721, 2002).

Semi-mechanistic population pharmacokinetic/pharmacodynamic model for neutropenia following therapy with the Plk-1 inhibitor BI 2536 and its application in clinical development

PURPOSE

(1) To describe the neutropenic response of BI 2536 a polo-like kinase 1 inhibitor in patients with cancer using a semi-mechanistic model. (2) To explore by simulations (a) the neutropenic effects for the maximum tolerated dose (MTD) and the dose at which dose-limiting toxicity occurred, (b) the possibility to reduce the cycle duration without increasing neutropenia substantially, and (c) the impact of the initial absolute neutrophil count (ANC) on the degree of neutropenia for different doses.

EXPERIMENTAL DESIGN

BI 2536 was administered as intravenous infusion over 60 min in the dose range from 25 to 250 mg. Three different administration schedules were explored: (a) day 1, (b) days 1, 2, and 3 or (c) days 1 and 8 within a 3 week treatment cycle. BI 2536 plasma concentrations and ANC obtained during the first treatment cycle from 104 patients were analysed using the population approach with NONMEM VI.

RESULTS

Neutropenia was described by a semi-mechanistic model resembling proliferation at the stem cell compartment, maturation, degradation, and homeostatic regulation. BI 2536 acts decreasing proliferation rate. Simulations showed that (1) all MTD doses showed an acceptable risk of neutropenia, (2) when BI 2536 is given as 200 mg single administration, cycle duration can be reduced from 3 to 2 weeks, and (3) baseline ANC might be considered to individualise the dose of BI 2536.

CONCLUSIONS

A semi-mechanistic population model was applied to describe the neutropenic effects of BI 2536. The model was used for simulations to support further clinical development.

Predictive ability of a semi-mechanistic model for neutropenia in the development of novel anti-cancer agents: two case studies

In cancer chemotherapy neutropenia is a common dose-limiting toxicity. An ability to predict the neutropenic effects of cytotoxic agents based on proposed trial designs and models conditioned on previous studies would be valuable. The aim of this study was to evaluate the ability of a semi-mechanistic pharmacokinetic/pharmacodynamic (PK/PD) model for myelosuppression to predict the neutropenia observed in Phase I clinical studies, based on parameter estimates obtained from prior trials. Pharmacokinetic and neutropenia data from 5 clinical trials for diflomotecan and from 4 clinical trials for indisulam were used. Data were analyzed and simulations were performed using the population approach with NONMEM VI. Parameter sets were estimated under the following scenarios: (a) data from each trial independently, (b) pooled data from all clinical trials and (c) pooled data from trials performed before the tested trial. Model performance in each of the scenarios was evaluated by means of predictive (visual and numerical) checks. The semi-mechanistic PK/PD model for neutropenia showed adequate predictive ability for both anti-cancer agents. For diflomotecan, similar predictions were obtained for the three scenarios. For indisulam predictions were better when based on data from the specific study, however when the model parameters were conditioned on data from trials performed prior to a specific study, similar predictions of the drug related-neutropenia profiles and descriptors were obtained as when all data were used. This work provides further indication that modeling and simulation tools can be applied in the early stages of drug development to optimize future trials.

Semi-mechanistic model for neutropenia after high dose of chemotherapy in breast cancer patients

PURPOSE

To describe the absolute neutrophil counts (ANC) profile in breast cancer patients receiving high-dose of chemotherapy and peripheral blood stem-cells (PBSC) transplantation.

METHODS

Data from 41 subjects receiving cyclophosphamide, thiotepa and carboplatin were used to develop the ANC model consisting of a drug-sensitive progenitor cell compartment, linked to the peripheral blood compartment, through three transition compartments. PBSC were incorporated into the first transit compartment following a zero-order process, k(in), and the rebound effect was explained by a feedback mechanism. A 'kinetics of drug action' model was used to quantify the HDC effect on the progenitor cells according to a linear function, with a slope (alpha).

RESULTS

The typical of the ANC at baseline (Circ(0)), mean transit time (MTT), feedback parameter (gamma), k(in) and alpha were estimated to be 5,610 x 10(6)/L, 3.25 days, 0.145, 0.954 cell/kg/day and 2.50 h/U, respectively. rHuG-CSF shortens the MTT by 92% and increases the mitotic activity by 120%. Bootstrap analysis, visual predictive check and numerical predictive checks evidenced accurate prediction of the ANC nadir, time to ANC nadir and time to grade 4 neutropenia recovery.

CONCLUSION

The time course of neutropenia following high-dose of chemotherapy and PBSC transplantation was accurately predicted. Higher amount of CD34+ cells in the PBSC transplantation and earlier administration rHuG-CSF were associated with faster haematological recovery.

Multiple-pool cell lifespan models for neutropenia to assess the population pharmacodynamics of unbound paclitaxel from two formulations in cancer patients

PURPOSE

Our objective was to build a mechanism-based pharmacodynamic model for the time course of neutropenia in cancer patients following paclitaxel treatment with a tocopherol-based Cremophor-free formulation (Tocosol Paclitaxel) and Cremophor EL-formulated paclitaxel (Taxol).

METHODS

A randomized two-way crossover trial was performed with 35 adult patients who received 175 mg/m(2) paclitaxel as either 15 min (Tocosol Paclitaxel) or 3 h (Taxol) intravenous infusions. Paclitaxel concentrations were measured by LC-MS/MS. NONMEM VI was used for population pharmacodynamics.

RESULTS

The cytotoxic effect on neutrophils was described by four mechanism-based models predicated on known properties of paclitaxel that used unbound concentrations in the central, deep peripheral or an intracellular compartment as forcing functions. Tocosol Paclitaxel was estimated to release 9.8% of the dose directly into the deep peripheral compartment (DPC). All models provided reasonable fitting of neutropenic effects. The model with the best predictive performance assumed that this dose fraction was released into 22.5% of the DPC which included the site of toxicity. The second-order cytotoxic rate constant was 0.00211 mL/ng per hour (variability: 52% CV). The relative exposure at the site of toxicity was 2.21 +/- 0.41 times (average +/- SD) larger for Tocosol Paclitaxel compared to Taxol. Lifespan was 11.0 days for progenitor cells, 1.95 days for maturating cells, and 4.38 days for neutrophils. Total drug exposure in blood explained half of the variance in nadir to baseline neutrophil count ratio.

CONCLUSIONS

The relative exposure of unbound paclitaxel at the site of toxicity was twice as large for Tocosol Paclitaxel compared to Taxol. The proposed mechanism-based models explained the extent and time course of neutropenia jointly for both formulations.

Using pharmacokinetic and pharmacodynamic modeling and simulation to evaluate importance of schedule in topotecan therapy for pediatric neuroblastoma

PURPOSE

The study aims to use mathematical modeling and simulation to assess the relative contribution of topotecan systemic exposure and scheduling in the activity and myelosuppression of topotecan in pediatric patients with neuroblastoma.

EXPERIMENTAL DESIGN

Pharmacokinetic and pharmacodynamic data were obtained from a phase II study for pediatric patients with high-risk neuroblastoma. The topotecan dosage was individualized to attain a topotecan lactone area under the plasma concentration-time curve between 80 and 120 ng/mL h and given over a protracted schedule (i.e., 10 days). Four mathematical models describing topotecan pharmacokinetics, tumor growth, and neutrophil and platelet dynamics were developed. The models were combined to simulate and compare different topotecan treatment strategies with respect to systemic exposure and schedule.

RESULTS

The median change in tumor volume was significantly different between schedules (5% increase for D x 5 versus 60% decrease for D x 5 x 2; P < 0.0001) when administering the same total systemic exposure. Whereas protracted schedules showed increased neutropenia (median of 7 versus 12 days below an absolute neutrophil count of 500/microL; P < 0.0001) and thrombocytopenia (median of 3 versus 10 days below a platelet count of 20,000/microL; P < 0.00001), simulations showed that delays in topotecan therapy would not be required. Simulations showed that an increase in topotecan exposure on the D x 5 schedule by 2.4-fold resulted in a modest decrease in tumor volume (i.e., median percentage change tumor volume of 24% versus 3%).

CONCLUSIONS

The present mathematical model gave an innovative approach to determining relevant topotecan schedules for possible evaluation in the clinic, which could lead to improved tumor response with minimized toxicities.

Pharmacodynamic model for chemotherapy-induced anemia in rats

Anticancer agents often cause bone marrow toxicity resulting in progressive anemia which may influence the therapeutic effects of erythropoietic-stimulating agents. The objective of this study was to develop a pharmacodynamic (PD) model to describe chemotherapy-induced anemia in rats. Anemia was induced in male Wistar rats with a single intravenous (i.v.) injection of 60 mg/kg carboplatin. Hematological responses including reticulocytes, red blood cells (RBC), hemoglobin, and endogenous rat erythropoietin (EPO) were measured for up to 4 weeks. A catenary, lifespan-based, indirect response model served as a basic PD model to represent erythroid cellular populations in the bone marrow and blood involved in erythropoiesis. The model assumed that actively proliferating progenitor cells in the bone marrow are sensitive to anti-cancer agents and subject to an irreversible removal process. The removal rate of the target cells is proportional to drug activity concentrations and the cell numbers. An additional RBC loss from the circulation resulting from thrombocytopenia was described by a first-order process. The turnover process of rat EPO and EPO-mediated feedback inhibition mechanism regulated by hemoglobin changes were incorporated. Reticulocyte counts decreased rapidly and reached a nadir by day 3 after administration of carboplatin and returned to the baseline by day 13. This was followed by a gradual increase and the rebound peak occurred at about day 15. The hemoglobin nadir was approximately 9 g/dl observed at about 11-13 days compared to its normal value of 13 g/dl and hemoglobin returned to the baseline by day 30. The increase in endogenous rat EPO mirrored inversely hemoglobin changes and the maximum increase was observed soon after the hemoglobin nadir. The carboplatin-treated rats exhibited progressive anemia. The proposed model adequately described the time course of hematological changes after carboplatin in rats and can be a useful tool to explore potential strategies for the management of anemia caused by chemotherapy.

Contribution of modelling chemotherapy-induced hematological toxicity for clinical practice

Anticancer chemotherapies are responsible for numerous adverse events. Among these, hematological toxicity is one of the main causes for ending treatment. These toxicities decrease production of red blood cells (anemia), production of white blood cells (neutropenia or granulocytopenia), and production of platelets (thrombocytopenia), which may be life-threatening to the patient. Preventing such discontinuation would be valuable for treating patients more effectively. In order to achieve this goal, numerous mathematical and physiological or semiphysiological models have been developed. The complexity of models has increased over the years, from empiric E(max) models to mechanistic models including physiological mechanisms such as feedback control. This review discusses several approaches of modelling hematological toxicities illustrated with some examples: pharmacodynamic models for the hematological toxicity of 5-fluorouracil, epirubicin, melphalan, paclitaxel, topotecan, and indisulam.

Population Pharmacokinetic-Pharmacodynamic Model for Neutropenia with Patient Subgroup Identification: Comparison across Anticancer Drugs

PURPOSE

Cancer chemotherapy, although based on body surface area, often causes unpredictable myelosuppression, especially severe neutropenia. The aim of this study was to evaluate qualitatively and quantitatively the influence of patient-specific characteristics on the neutrophil concentration-time course, to identify patient subgroups, and to compare covariates on system-related pharmacodynamic variable between drugs.

EXPERIMENTAL DESIGN

Drug and neutrophil concentration, demographic, and clinical chemistry data of several trials with docetaxel (637 patients), paclitaxel (45 patients), etoposide (71 patients), or topotecan (191 patients) were included in the covariate analysis of a physiology-based pharmacokinetic-pharmacodynamic neutropenia model. Comparisons of covariate relations across drugs were made.

RESULTS

A population model incorporating four to five relevant patient factors for each drug to explain variability in the degree and duration of neutropenia has been developed. Sex, previous anticancer therapy, performance status, height, binding partners, or liver enzymes influenced system-related variables and alpha1-acid glycoprotein, albumin, bilirubin, concomitant cytotoxic agents, or administration route changed drug-specific variables. Overall, female and pretreated patients had a lower baseline neutrophil concentration. Across-drug comparison revealed that several covariates (e.g., age) had minor (clinically irrelevant) influences but consistently shifted the pharmacodynamic variable in the same direction.

CONCLUSIONS

These mechanistic models, including patient characteristics that influence drug-specific parameters, form the rationale basis for more tailored dosing of individual patients or subgroups to minimize the risk of infection and thus might contribute to a more successful therapy. In addition, nonsignificant or clinically irrelevant relations on system-related parameters suggest that these covariates could be negligible in clinical trails and daily use.

Semi-physiological model describing the hematological toxicity of the anti-cancer agent indisulam

Indisulam (N-(3-chloro-7-indolyl)-1,4-benzenedisulfonamide, GOAL, E7070) is a novel anti-cancer drug currently in phase II clinical development for the treatment of solid tumors. Phase I dose-escalation studies were conducted comparing four treatment schedules. Neutropenia and thrombocytopenia were dose limiting in all schedules. The aim of this study was to describe the extent and the time course of the hematological toxicity and its possible schedule dependency using a semi-physiological model. Data from 142 patients were analyzed using NONMEM. The semi-physiological model comprised a progenitor blood cell compartment, linked to the central circulation compartment, through 3 transition compartments representing the maturation chain in the bone marrow. Plasma concentrations of the drug were assumed to reduce the proliferation rate in the progenitor compartment according to a linear function. A feedback mechanism was included in the model representing the rebound effect of endogenous growth factors. The model was validated using a posterior predictive check. The model adequately described the extent and time course of neutropenia and thrombocytopenia. The mean transition time (MTT, i.e. maturation time in bone marrow) of neutrophils was increased by 47% in patients who received indisulam as a weekly dose administered for four out of every six weeks. For platelets, MTT was increased by 33% in patients who received this schedule and also in patients who received a continuous 120-h infusion. The validation procedure indicated that the model adequately predicts the nadir value of neutrophils and platelets and the time to reach this nadir. A semi-physiological model was successfully applied to describe the time course and extent of the neutropenia and thrombocytopenia after indisulam administration for four treatment schedules.

Comparison of pharmacokinetics and pharmacodynamics of docetaxel and Cisplatin in elderly and non-elderly patients: why is toxicity increased in elderly patients?

PURPOSE

Following phase I studies of docetaxel and cisplatin in patients with non-small-cell lung cancer, the recommended doses of docetaxel were different for elderly (> or = 75 years) and non-elderly (< 75 years) patients. To elucidate the mechanism of the difference, the pharmacokinetics of docetaxel and cisplatin were investigated in two phase II studies separately conducted in elderly and non-elderly patients.

PATIENTS AND METHODS

Twenty-seven elderly and 25 non-elderly patients were treated with three weekly administrations of docetaxel and cisplatin every 4 weeks. Doses of docetaxel were 20 and 35 mg/m(2) for elderly and non-elderly patients, respectively. All patients received 25 mg/m(2) of cisplatin. The pharmacokinetics and pharmacodynamics of docetaxel and cisplatin were compared in elderly and non-elderly patients.

RESULTS

There were no differences in pharmacokinetics of docetaxel or cisplatin between elderly versus non-elderly patients with regard to clearance and volume of distribution. In the pharmacodynamic analysis, neutropenia was positively correlated with the area under the concentration-time curve for docetaxel but not for cisplatin. In evaluating the relationship between neutropenia and the area under the concentration-time curve of docetaxel, elderly patients experienced greater neutropenia than those predicted by a pharmacodynamic model developed in non-elderly patients; the residual for prediction of the percent change in neutrophil count was -11.2% (95% CI, -21.8 to -0.5%).

CONCLUSION

The pharmacokinetics of docetaxel and unchanged cisplatin were not different between elderly and non-elderly patients. The elderly patients were more sensitive to docetaxel exposure than the non-elderly patients, resulting in the different recommended doses for the phase II studies.

Pharmacokinetic/Pharmacodynamic Analysis of Neutrophil Proliferation Induced by rhG-CSF in Patients Receiving Antineoplastic Drugs

In the present study, we analyzed the effect of recombinant human granulocyte colony-stimulating factor (rhG-CSF) on neutrophil counts in cancer patients undergoing chemotherapy using a previously developed pharmacokinetic/pharmacodynamic model.(7)) The time profiles of neutrophil counts in blood after repeated administration of rhG-CSF to lung cancer patients undergoing chemotherapy could be analyzed by this model by considering the inhibition of neutrophil production by antineoplastic drugs. Although deviation was observed between the predicted and observed neutrophil counts in ovarian cancer patients, it may be possible to use this model for determining a rational dosage regimen of rhG-CSF for patients undergoing chemotherapy.

Effect of valspodar on the pharmacokinetics of unbound paclitaxel

The aim of this multicenter study was to determine whether valspodar (Amdray; code designation, SDZ PSC 833), a potent P-glycoprotein (P-gp) inhibitor, affects the pharmacokinetics of unbound paclitaxel (Cu). Data were obtained from 31 patients with advanced breast cancer. Thirteen patients were treated with paclitaxel alone (3-h infusion at 175 mg/m2) and another 18 received paclitaxel (3-h infusion at 70 mg/m2) in combination with a 21-day cycle of oral valspodar (5 mg/kg given four times a day) starting 1 day before administration of paclitaxel. Serial blood samples were taken in the first course and Cu in plasma determined using equilibrium dialysis with a [G-3H]paclitaxel tracer. The apparent clearance of Cu was not significantly different between the two groups, with mean +/- standard deviation (+/- SD) values of 230 +/- 56.0 and 202 +/- 49.9 L/h/m2 in the absence and presence of valspodar, respectively (P = 0.17). The volume of Cu distribution was slightly larger in the presence of valspodar (1160 +/- 474 vs. 1620 +/- 552 L/m2; P = 0.025), which contributed to a minor difference in the terminal disposition half-life (6.12 +/- 3.42 vs. 8.50 +/- 2.06 h; P = 0.028). These data indicate that (i) valspodar lacks the significant interaction with paclitaxel observed previously with other P-gp modulators, (ii) the majority of the increased toxicity of the combination does not appear to be attributable to increased levels of Cu, and (iii) provide further evidence of the conjecture that the plasma concentration of paclitaxel may not be an appropriate measure to monitor the impact of P-gp inhibition.

A mechanistic mathematical model of temozolomide myelosuppression in children with high-grade gliomas

Temozolomide (TMZ) is currently being evaluated for the treatment of high-grade gliomas in children. Myelosuppression (the suppression of bone marrow activity) is the dose-limiting toxicity for TMZ in adults and children. Empirical methods (i.e. relations between the percent change in absolute neutrophil count (ANC) and the area under the plasma concentration curve (AUC) of TMZ or its active metabolite MTIC) showed poor results when attempting to describe myelosuppression from serial data derived during TMZ therapy in a Phase II study of children with high-grade glioma. Therefore, to improve our understanding of the myelosuppressive effects of TMZ and MTIC in children we developed a mechanistic mathematical model. The model describes the progression of neutrophils from their production in the bone marrow to their release in the plasma. Included in the model are the feedback effects of granulocyte colony stimulating factor (G-CSF), which stimulates neutrophil production when there is a decrease in circulating neutrophils. The model is fit to serial ANC measurements obtained after TMZ dosing and it is able to explain, among other things, the lag in ANC reduction following a dose of TMZ, the ANC nadir, and the 'rebound effect' observed where the ANC recovers to levels greater than that observed pre-TMZ dose. This model will be useful for the prospective design of clinical trials of TMZ in children with cancer.

Comparison of neutropenia in a randomized, crossover trial of 3-, 6-, and 24-h infusions of paclitaxel

OBJECTIVE

Paclitaxel is most commonly infused over 3 hs rather than the original schedule of 24 h as the briefer infusion duration results in greater convenience, similar efficacy, significantly less myelosuppression, and less cost. While differences in toxicity between 3- and 24-h infusions are well described, there is little information about the effect of modest prolongation of infusion duration, which is often employed in patients who develop hypersensitivity reactions. To assess whether prolonging a 3-h infusion significantly increases the degree of neutropenia, we reviewed our data from a randomized, crossover trial of 3-h versus 6-h versus 24-h regimens of paclitaxel.

METHODS

Results from 12 patients who were randomized to receive one 3-h, one 6-h and one 24-h infusion of paclitaxel in varied sequences during their first three cycles of treatment were analysed. The blood counts were monitored closely throughout each cycle of treatment.

RESULTS

Crossover trial methodology was used to assess the differences in the degree of neutropenia caused by the three different infusion regimens. The 24-h infusion regimen resulted in significantly worse neutropenia than the 3- or 6-h infusion regimens. There was no statistically significant difference between the 3- and 6-h infusion regimens with respect to all endpoints. The estimated mean difference in the duration of grade 3 or 4 neutropenia between the 3- and 6-h infusion regimens (6 h - 3-h) was 1.1 day (95% CI: -0.9, 3.2), and for grade 4 neutropenia, the estimated mean difference in the duration was 0.8 day (95% CI: -0.4, 2.0).

CONCLUSIONS

Increasing the duration of paclitaxel infusions from 3 to 6 h does not result in a statistically significant increase in the degree of neutropenia. Any additional neutropenia is likely to be of brief duration.

Pharmacokinetic-pharmacodynamic modeling of methotrexate-induced toxicity in mice

The prediction of chemotherapeutic efficacy is complicated by "protocol dependencies" in dose-effect and dose-toxicity relationships. It has been proposed that pharmacokinetic-pharmacodynamic mathematical models may allow characterization of chemotherapeutic protocol dependencies, and may facilitate the prediction of chemotherapeutic efficacy; however, few demonstrations exist in the literature. The present study examines the pharmacokinetics and toxicodynamics of methotrexate (MTX), a commonly used anticancer agent, after intraperitoneal (i.p.) administration to mice. MTX was administered via bolus or infusion (24, 72, and 168 h), at doses of 2.5-1000 mg/kg. MTX plasma and peritoneal pharmacokinetics were characterized through standard noncompartmental and compartmental techniques. Body weight loss was used as a measure of MTX-induced toxicity. We found that MTX pharmacokinetics were independent of dose (over a range of 3-600 mg/kg) and independent of dosing mode (i.e., i.p. bolus vs. i.p. infusion). However, MTX-induced toxicity was shown to be highly dependent on the dosing protocol used. For example, the maximally tolerated dose (i.e., the dose related to a mean body weight loss of 10%) was 200-fold greater after bolus administration relative to that observed for 72-h infusion (760 mg/kg vs. 3.8 mg/kg). This profound protocol dependence in the relationship between MTX-induced toxicity and MTX exposure was characterized through the use of a time-dissociated pharmacokinetic-pharmacodynamic model (median prediction error: 3.9%).

Mechanistic models for myelosuppression

As myelosuppression is the dose-limiting toxicity for most chemotherapeutic drugs, modelers attempt to find relationships between drug and toxicity to optimize treatment. Mechanistic models, i.e. models based on physiology and pharmacology, are preferable over empirical models, as prior information can be utilized and as they generally are more reliable for extrapolations. To account for different dosing-regimens and possible schedule-dependent effects, the whole concentration-time profile should be used as input into the pharmacokinetic-pharmacodynamic model. It is also of importance to model the whole time course of myelosuppression to be able to predict both the degree and duration of toxicity as well as consecutive courses of therapy. A handful of (semi)-mechanistic pharmacokinetic-pharmacodynamic models with the above properties have been developed and are reviewed. Ideally, a model of myelosuppression should separate drug-specific parameters from system related parameters to be applicable across drugs and useful under different clinical settings. Introduction of mechanistic models of myelosuppression in the design and evaluation of clinical trials can guide in the decision of optimal sampling times, contribute to knowledge of optimal doses and treatment regimens at an earlier time point and identify sub-groups of patients at a high risk of myelosuppression.

Diversity of mechanism-based pharmacodynamic models

Pharmacodynamics is the study of the time course of pharmacological effects of drugs. The field of pharmacodynamic modeling has made many advances, due in part to the relatively recent development of basic and extended mechanism-based models. The purpose of this article is to describe the classic as well as contemporary approaches, with an emphasis on pertinent equations and salient model features. In addition, current methods of integrating various system complexities into these models are discussed. Future pharmacodynamic models will most likely reflect an assembly of the basic components outlined in this review.

Population pharmacokinetics/pharmacodynamics relationships of an anticancer drug

This paper proposes a method for studying the toxicity of an anticancer drug with a delayed effect. The goal is to predict a dosage regimen with controlled toxicity. To this end, a semi-physiological model is used. A limit of toxicity is demonstrated, which is intrinsic to the model. It reduces the effect of high drug concentrations. This limit explains the mixed behaviour of the drug: time-dependence and concentration-dependence, according to the dose actually administered. A population analysis is performed to estimate the parameters of the model, and to predict a safe dosage regimen.

Model of chemotherapy-induced myelosuppression with parameter consistency across drugs

PURPOSE

To develop a semimechanistic pharmacokinetic-pharmacodynamic model describing chemotherapy-induced myelosuppression through drug-specific parameters and system-related parameters, which are common to all drugs.

PATIENTS AND METHODS

Patient leukocyte and neutrophil data after administration of docetaxel, paclitaxel, and etoposide were used to develop the model, which was also applied to myelosuppression data from 2'-deoxy-2'-methylidenecytidine (DMDC), irinotecan (CPT-11), and vinflunine administrations. The model consisted of a proliferating compartment that was sensitive to drugs, three transit compartments that represented maturation, and a compartment of circulating blood cells. Three system-related parameters were estimated: baseline, mean transit time, and a feedback parameter. Drug concentration-time profiles affected the proliferation of sensitive cells by either an inhibitory linear model or an inhibitory E(max) model. To evaluate the model, system-related parameters were fixed to the same values for all drugs, which were based on the results from the estimations, and only drug-specific parameters were estimated. All modeling was performed using NONMEM software.

RESULTS

For all investigated drugs, the model successfully described myelosuppression. Consecutive courses and different schedules of administration were also well characterized. Similar system-related parameter estimates were obtained for the different drugs and also for leukocytes compared with neutrophils. In addition, when system-related parameters were fixed, the model well characterized chemotherapy-induced myelosuppression for the different drugs.

CONCLUSION

This model predicted myelosuppression after administration of one of several different chemotherapeutic drugs. In addition, with fixed system-related parameters to proposed values, and only drug-related parameters estimated, myelosuppression can be predicted. We propose that this model can be a useful tool in the development of anticancer drugs and therapies.

Multiple-pool cell lifespan model of hematologic effects of anticancer agents

The leukopenic effects of anticancer agents are described using a semi-physiologic multiple-pool cell lifespan model. The time course of myelosuppression in relation to the drug concentration vs. time profile was characterized using a three pool indirect model. The proliferation and maturation stages of myeloid cells in the bone marrow and cell removal from the circulation were quantitated with a cell life-span concept. Drug effects were assumed to take place in the bone marrow based on irreversible linear or capacity-limited cytotoxicity. Mathematical derivations and computer simulations (Adapt II) were used to examine the properties of the model. Data from the literature were also analyzed. Cell response profiles after therapy typically exhibit a lag period, reduction to a nadir, and return to baseline. The predicted values of the time periods of granulopoiesis were 10-14 days for proliferation, and 1-6 days for maturation of progenitor cells in the bone marrow. The proposed irreversible mechanism of cell killing by anticancer drugs explains previously observed relationships between leukocyte nadir counts and exposure to the drug and/or duration of drug concentrations above some threshold level. The model was applied to literature data for paclitaxel and etoposide effects on leukocyte counts. The predicted value of KC50 for paclitaxel ranged from 0.004 to 0.2 microgram/mL and for etoposide 2 micrograms/mL. The present model accounts for drug-induced leukopenia using a physiologic cell production and loss model and irreversible cytotoxicity in a precursor pool.

Irinotecan pharmacokinetics-pharmacodynamics: the clinical relevance of prolonged exposure to SN-38

We have shown previously that the terminal disposition half-life of SN-38, the active metabolite of irinotecan, is much longer than earlier thought. Currently, it is not known whether this prolonged exposure has any relevance toward SN-38-induced toxicity. Here, we found that SN-38 concentrations present in human plasma for up to 3 weeks after a single irinotecan infusion induce significant cytotoxicity in vitro. Using pharmacokinetic data from 26 patients, with sampling up to 500 h, relationships were evaluated between systemic exposure (AUC) to SN-38 and the per cent decrease in absolute neutrophil count (ANC) at nadir, or by taking the entire time course of ANC into account (AOC). The time course of SN-38 concentrations (AUC(500 h)) was significantly related to this AOC (P<0.001). Based on these findings, a new limited-sampling model was developed for SN-38 AUC(500 h) using only two timed samples: AUC(500 h)=(6.588 x C(2.5 h))+(146.4 x C(49.5 h))+15.53, where C(2.5 h) and C(49.5 h) are plasma concentrations at 2.5 and 49.5 h after start of infusion, respectively. The use of this limited-sampling model may open up historic databases to retrospectively obtain information about SN-38-induced toxicity in patients treated with irinotecan.

In vitro toxicity of ET-743 and aplidine, two marine-derived antineoplastics, on human bone marrow haematopoietic progenitors. comparison with the clinical results

Ecteinascidine-743 (ET-743) and aplidine are two marine-derived antineoplastics currently in phase II development. With the aim of evaluating whether in vitro haematopoietic studies can predict the toxicity of these two drugs in patients, human bone marrow (BM) samples were incubated with these drugs under conditions which mimicked the administration exposures used in the clinics. As it was observed in different cancer cell lines, ET-743 was more toxic on an equimolar basis in human hematopoietic progenitors (inhibitory concentration reducing the viability to 50% after 24 h exposures; IC50(24h): 10-50 nM) compared with doxorubicin (IC50(24h) values: 280-460 nM), used as a control anticancer drug. In contrast to the high haematotoxic effects observed for ET-743, similar IC values were obtained for aplidine (IC50(24h): 150-530 nM) compared with doxorubicin. For both ET-743 and aplidine, the megakaryocytic progenitor was the most sensitive, compared with the other haematopoietic progenitors (IC50 values were 3- to 5-fold lower in the CFU-Megs compared with the CFU-GMs). The observation that the Cmax observed in patients treated with the aplidine maximum tolerated dose (MTD) (7.1 nM) was 21-75 fold lower than the IC50(24h) value observed for the different haematopoietic progenitors is highly consistent with the lack of haematotoxicity observed in patients treated with this drug. In the case of ET-743, differences between the Cmax value corresponding to the MTD (2.6 nM) and the in vitro IC50 values corresponding to the different progenitors were much lower (4-19-fold), also consistent with the haematotoxicity that was observed in patients treated at recommended doses (RDs) and MTDs. Although CFU-Megs were more sensitive than CFU-GM progenitors to ET-743 in vitro, clinical data showed that neutropenic events were more frequent than thrombocytopenic episodes. Aiming to further improve the predictive value of in vitro IC values corresponding to the different haematopoietic progenitors, additional refinement parameters derived from pharmacokinetic and animal studies are proposed.

Pharmacodynamic modeling of chemotherapeutic effects: application of a transit compartment model to characterize methotrexate effects in vitro

The time course of chemotherapeutic effect is often delayed relative to the time course of chemotherapeutic exposure. In many cases, this delay is difficult to characterize mathematically through the use of standard pharmacodynamic models. In the present work, we investigated the relationship between methotrexate (MTX) exposure and the time course of MTX effects on tumor cell growth in culture. Two cancer cell lines, Ehrlich ascites cells and sarcoma 180 cells, were exposed for 24 hours to MTX concentrations that varied more than 700-fold (0.19-140 micro g/mL). Viable cells were counted on days 1, 3, 5, 7, 9, 11, 13, 15, 17, 20, 22, and 24 for Ehrlich ascites cells and on days 1, 2, 3, 5, 7, 9, 11, 13, 14, 15, 17, 19, and 21 for sarcoma 180 cells, through the use of a tetrazolium assay. Although MTX was removed 24 hours after application, cell numbers reached nadir values more than 100 hours after MTX exposure. Data from each cell line were fitted to 3 pharmacodynamic models of chemotherapeutic cell killing: a cell cycle phase-specific model, a phase-nonspecific model, and a transit compartment model (based on the general model recently reported by Mager and Jusko, Clin Pharmacol Ther. 70:210-216, 2001). The transit compartment model captured the data much more accurately than the standard pharmacodynamic models, with correlation coefficients ranging from 0.86 to 0.999. This report shows the successful application of a transit compartment model for characterization of the complex time course of chemotherapeutic effects; such models may be very useful in the development of optimization strategies for cancer chemotherapy.

Mechanism-based pharmacokinetic model for paclitaxel

PURPOSE

To create a model based on known mechanisms of paclitaxel distribution that could describe the pharmacokinetics (PK) of total and unbound plasma concentrations, as well as blood concentrations. In addition, to investigate the relationship between exposure, based on unbound and total concentrations, and neutropenia.

PATIENTS AND METHODS

Paclitaxel and Cremophor EL (CrEL) concentrations were obtained from 23 female and three male patients (50 courses in total) with different cancer types that received paclitaxel (Taxol; Bristol-Myers Squibb Co, Princeton, NJ) (135 to 225 mg/m(2)) as 3- or 24-hour intravenous infusions. Seven of the patients received combination therapy with doxorubicin or cisplatin. The population PK model was built to fit three types of data simultaneously: unbound, total plasma, and blood concentrations. The area under the curve, threshold, and general models were used to relate neutrophil survival fraction from 19 patients (29 courses in total) to exposure based on unbound and total plasma concentration, respectively.

RESULTS

The PK model included a linear three-compartment model for unbound concentration, binding directly proportional to CrEL, linear and nonlinear binding to plasma proteins, and linear and nonlinear binding to blood cells. The threshold model best described the PK/pharmacodynamic (PD) relationship for total concentration. No distinction could be made between the models for unbound drug.

CONCLUSION

Earlier PK models for paclitaxel have been empirical. This study shows that a mechanistic model can be used to describe the nonlinear PK of paclitaxel. There is an indication that the PK/PD relationship is not the same for unbound and total plasma concentrations.

Paclitaxel pharmacodynamics: application of a mechanism-based neutropenia model

Antineoplastic agents exert adverse effects that impact both dose and scheduling of drug administration. Our objective was to develop a quantitative relationship between paclitaxel (taxol) exposure and pharmacodynamic endpoints, such as neutropenia or body weight loss. Paclitaxel in liposomes or Cremophor EL was administered to rats at doses of 20 or 40 mg/kg. Body weight and absolute neutrophil count were determined daily. The decrease in body weight was greater for paclitaxel in Cremophor EL than for liposomal paclitaxel, but hematological toxicity was similar. The hematological data was fit using a pharmacodynamic model to investigate the temporal delay between drug exposure and neutropenia. From the model, the lifespan of neutrophils (T(N)), of surviving precursor cells in bone marrow (T(P)), and a killing rate constant (K) were determined. The values of T(N), T(P), and K for liposomal paclitaxel were 95 h, 82 h, and 0.735 (microM h)(-1), respectively, and for paclitaxel in Cremophor EL, 86 h, 78 h, and 0.475 (microM h)(-1), respectively. Simulations of various doses indicated a dependency of the neutropenia time course on paclitaxel exposure. The entire time course of changes in neutrophil count is more informative than a single measurement if myelosuppression is prolonged and at a level associated with increased incidence of clinical adverse effects.