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

Liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) methods for the therapeutic drug monitoring of cytotoxic anticancer drugs: An update

In the era of precision medicine, there is increasing evidence that conventional cytotoxic agents may be suitable candidates for therapeutic drug monitoring (TDM)- guided drug dosage adjustments and patient's tailored personalization of non-selective chemotherapies. To that end, many liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) assays have been developed for the quantification of conventional cytotoxic anticancer chemotherapies, that have been comprehensively and critically reviewed. The use of stable isotopically labelled internal standards (IS) of cytotoxic drugs was strikingly uncommon, accounting for only 48 % of the methods found, although their use could possible to suitably circumvent patients' samples matrix effects variability. Furthermore, this approach would increase the reliability of cytotoxic drug quantification in highly multi-mediated cancer patients with complex fluctuating pathophysiological and clinical conditions. LC-MS/MS assays can accommodate multiplexed analyses of cytotoxic drugs with optimal selectivity and specificity as well as short analytical times and, when using stable-isotopically labelled IS for quantification, provide concentrations measurements with a high degree of certainty. However, there are still organisational, pharmacological, and medical constraints to tackle before TDM of cytotoxic drugs can be more largely adopted in the clinics for contributing to our ever-lasting quest to improve cancer treatment outcomes.

Drug Exposure to Establish Pharmacokinetic-Response Relationships in Oncology

In the oncology field, understanding the relationship between the dose administered and the exerted effect is particularly important because of the narrow therapeutic index associated with anti-cancer drugs and the high interpatient variability. Therefore, in this review, we provide a critical perspective of the different methods of characterising treatment exposure in the oncology setting. The increasing number of modelling applications in oncology reflects the applicability and the impact of pharmacometrics on all phases of the drug development process and patient management as well. Pharmacometric modelling is a worthy component within the current paradigm of model-based drug development, but pharmacometric modelling techniques are also accessible for the clinician in the optimisation of current oncology therapies. Consequently, the application of population models in a hospital setting by generating close collaborations between physicians and pharmacometricians is highly recommended, providing a systematic means of developing and assessing model-based metrics as 'drivers' for various responses to treatments, which can then be evaluated as predictors for treatment success. Characterising the key determinants of variability in exposure is of particular importance for anticancer agents, as efficacy and toxicity are associated with exposure. We present the different strategies to describe and predict drug exposure that can be applied depending on the data available, with the objective of obtaining the most useful information in the patients' favour throughout the full drug cycle. Therefore, the objective of the present article is to review the different approaches used to characterise a patient's exposure to oncology drugs, which will result in a better understanding of the time course of the response and the magnitude of interpatient variability.

A physiologically based pharmacokinetic - pharmacodynamic modelling approach to predict incidence of neutropenia as a result of drug-drug interactions of paclitaxel in cancer patients

Paclitaxel is the backbone of standard chemotherapeutic regimens used in a number of malignancies and is frequently given with concomitant medications. Newly developed oncolytic agents, including tyrosine kinase inhibitors are often shown to be CYP3A4 and P-gp inhibitors. The aim of this study was to develop a PBPK model for intravenously administered paclitaxel in order to predict the incidence of neutropenia and to estimate the DDI potential as a victim drug. The dose-dependent effects on paclitaxel plasma protein binding, volume of distribution and drug clearance were considered for dose levels of 80 mg/m², 135 mg/m² and 175 mg/m². A pharmacodynamics model that incorporate the impact of paclitaxel on the neutrophil was developed. The relative metabolic clearance via CYP3A4 and CYP2C8, the renal clearance as well as P-gp mediated biliary clearance were incorporated in the model in order to assess the neutropenia in the presence of DDI. The developed PBPK-PD model was able to recover the drop in neutrophils observed after the administration of 175mg/m2 of paclitaxel over a 3-h duration. The mean nadir observed was 1.9 × 10⁹ neutrophils/L and was attained after 10 days of treatment, and a fraction of 47% of the population was predicted to have at some point a neutropenia including 12% with severe neutropenia. In the case of concomitant administration of ketoconazole, 39% of the population was predicted to suffer from severe neutropenia. In summary, PBPK-PD modeling allows a priori prediction of DDIs and safety events involving complex combination therapies which are often utilized in an oncology setting.

Effects of paclitaxel and resveratrol on blood characteristics in rabbits

We investigated the adverse effects of paclitaxel (PAC) on blood characteristics including the percentage of alpha naphthyl acetate esterase (ANAE)-positive peripheral blood lymphocytes, white cell count, number of neutrophil nuclear lobe and, erythrocyte diameter, and the protective effects of resveratrol (RES) on these characteristics. Male New Zealand rabbits were assigned to four groups. The control group (C) was given 40 ml normal saline, the PAC group was given 5 mg/kg PAC in 40 ml normal saline, the RES group was given 4 mg/kg RES in 40 ml normal saline, and the PAC + RES group was given 5 mg/kg PAC + 4 mg/kg RES in 40 ml normal saline. All injections were performed intravenously once/week for 8 weeks. PAC caused an increased total percentage of lymphocytes and monocytes, which was associated with neutropenia. The PAC group showed a right shift in the lobulation curve of neutrophil nuclei together with a decreased percentage of ANAE-positive lymphocytes. RES reduced the percentage of ANAE-positive lymphocytes, slightly increased the percentage of neutrophil leukocytes and reduced the total lymphocyte count. Administration of RES combined with PAC prevented, while PAC induced, neutropenia and lymphocytosis, and increased the percentage of ANAE-positive lymphocytes.

Towards Quantitative Systems Pharmacology Models of Chemotherapy-Induced Neutropenia

Neutropenia is a serious toxic complication of chemotherapeutic treatment. For years, mathematical models have been developed to better predict hematological outcomes during chemotherapy in both the traditional pharmaceutical sciences and mathematical biology disciplines. An increasing number of quantitative systems pharmacology (QSP) models that combine systems approaches, physiology, and pharmacokinetics/pharmacodynamics have been successfully developed. Here, I detail the shift towards QSP efforts, emphasizing the importance of incorporating systems-level physiological considerations in pharmacometrics.

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.

Natural Products as a Vital Source for the Discovery of Cancer Chemotherapeutic and Chemopreventive Agents

Throughout history, natural products have played a dominant role in the treatment of human ailments. For example, the legendary discovery of penicillin transformed global existence. Presently, natural products comprise a large portion of current-day pharmaceutical agents, most notably in the area of cancer therapy. Examples include Taxol, vinblastine, and camptothecin. These structurally unique agents function by novel mechanisms of action; isolation from natural sources is the only plausible method that could have led to their discovery. In addition to terrestrial plants as sources for starting materials, the marine environment (e.g., ecteinascidin 743, halichondrin B, and dolastatins), microbes (e.g., bleomycin, doxorubicin, and staurosporin), and slime molds (e.g., epothilone B) have yielded remarkable cancer chemotherapeutic agents. Irrespective of these advances, cancer remains a leading cause of death worldwide. Undoubtedly, the prevention of human cancer is highly preferable to treatment. Cancer chemoprevention, the use of vaccines or pharmaceutical agents to inhibit, retard, or reverse the process of carcinogenesis, is another important approach for easing this formidable public health burden. Similar to cancer chemotherapeutic agents, natural products play an important role in this field. There are many examples, including dietary phytochemicals such as sulforaphane and phenethyl isothiocyanate (cruciferous vegetables) and resveratrol (grapes and grape products). Overall, natural product research is a powerful approach for discovering biologically active compounds with unique structures and mechanisms of action. Given the unfathomable diversity of nature, it is reasonable to suggest that chemical leads can be generated that are capable of interacting with most or possibly all therapeutic targets.

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.

Clinical population pharmacokinetics and toxicodynamics of linezolid

Thrombocytopenia is a common side effect of linezolid, an oxazolidinone antibiotic often used to treat multidrug-resistant Gram-positive bacterial infections. Various risk factors have been suggested, including linezolid dose and duration of therapy, baseline platelet counts, and renal dysfunction; still, the mechanisms behind this potentially treatment-limiting toxicity are largely unknown. A clinical study was conducted to investigate the relationship between linezolid pharmacokinetics and toxicodynamics and inform strategies to prevent and manage linezolid-associated toxicity. Forty-one patients received 42 separate treatment courses of linezolid (600 mg every 12 h). A new mechanism-based, population pharmacokinetic/toxicodynamic model was developed to describe the time course of plasma linezolid concentrations and platelets. A linezolid concentration of 8.06 mg/liter (101% between-patient variability) inhibited the synthesis of platelet precursor cells by 50%. Simulations predicted treatment durations of 5 and 7 days to carry a substantially lower risk than 10- to 28-day therapy for platelet nadirs of <100 ×10(9)/liter. The risk for toxicity did not differ noticeably between 14 and 28 days of therapy and was significantly higher for patients with lower baseline platelet counts. Due to the increased risk of toxicity after longer durations of linezolid therapy and large between-patient variability, close monitoring of patients for development of toxicity is important. Dose individualization based on plasma linezolid concentration profiles and platelet counts should be considered to minimize linezolid-associated thrombocytopenia. Overall, oxazolidinone therapy over 5 to 7 days even at relatively high doses was predicted to be as safe as 10-day therapy of 600 mg linezolid every 12 h.

Prediction of nanoparticle prodrug metabolism by pharmacokinetic modeling of biliary excretion

Pharmacokinetic modeling and simulation is a powerful tool for the prediction of drug concentrations in the absence of analytical techniques that allow for direct quantification. The present study applied this modeling approach to determine active drug release from a nanoparticle prodrug formulation. A comparative pharmacokinetic study of a nanoscale micellar docetaxel (DTX) prodrug, Procet 8, and commercial DTX formulation, Taxotere, was conducted in bile duct cannulated rats. The nanoscale (~40nm) size of the Procet 8 formulation resulted in confinement within the plasma space and high prodrug plasma concentrations. Ex vivo prodrug hydrolysis during plasma sample preparation resulted in unacceptable error that precluded direct measurement of DTX concentrations. Pharmacokinetic modeling of Taxotere and Procet 8 plasma concentrations, and their associated biliary metabolites, allowed for prediction of the DTX concentration profile and DTX bioavailability, and thereby evaluation of Procet 8 metabolism. Procet 8 plasma decay and in vitro plasma hydrolytic rates were identical, suggesting that systemic clearance of the prodrug was primarily metabolic. The Procet 8 and Taxotere plasma profiles, and associated docetaxel hydroxy-tert-butyl carbamate (HDTX) metabolite biliary excretion, were best fit by a two compartment model, with both linear and non-linear DTX clearance, and first order Procet 8 hydrolysis. The model estimated HDTX clearance rate agreed with in vitro literature values, supporting the predictability of the proposed model. Model simulation at the 10mg DTX equivalent/kg dose level predicted DTX formation rate-limited kinetics and a peak plasma DTX concentration of 39ng/mL at 4h for Procet 8, in comparison to 2826ng/mL for Taxotere. As a result of nonlinear DTX clearance, the DTX AUCinf for the Procet 8 formulation was predicted to be 2.6 times lower than Taxotere (775 vs. 2017h×ng/mL, respectively), resulting in an absolute bioavailability estimate of 38%. As DTX clearance in man is considered linear, this low bioavailability is likely species-dependent. These data support the use of pharmacokinetic modeling and simulation in cases of complex formulations, where analytical methods for direct measurement of free (released) drug concentrations are unavailable. Uses of such models may include interpretation of preclinical toxicology studies, selection of first in man dosing regimens, and PK/PD model development.

Meta-analysis of nanoparticulate paclitaxel delivery system pharmacokinetics and model prediction of associated neutropenia

PURPOSE

Nanoparticulate paclitaxel carriers have entered clinical evaluation as alternatives to the Cremophor-based standard Taxol(®) (Cre-pac). Their pharmacokinetics (PK) is complex, and factors influencing their pharmacodynamics (PD) are poorly understood. We aimed to develop a unified quantitative model for 4 paclitaxel carriers that captures systems-level PK, predicts micro-scale PK processes, and permits correlations between carrier properties and observed PD.

METHODS

Data consisting of 54 PK profiles and 574 observations were extracted from 20 clinical studies investigating Cre-pac, albumin-(A-pac), liposome-(L-pac), and tocopherol-(T-pac) nanocarriers. A population-PK approach was used for data analysis. All datasets were simultaneously fitted to produce a unified model. Model-based simulations explored relationships between predicted PK and myelosuppression for each formulation.

RESULTS

The final model employed nonlinear drug-binding mechanisms to describe Cre-pac and a delayed-release model for A-pac, L-pac, and T-pac. Estimated drug-release rate constants (h(-1)): Cre-pac (5.19), L-pac (1.26), A-pac (0.72), T-pac (0.74). Simulations of equivalent dosing schemes ranked neutropenia severity (highest to lowest): T-pac~Cre-pac>L-pac~A-pac and predicted remarkably well the clinically-observed relationships between neutropenia and free drug exposure relative to a threshold concentration.

CONCLUSIONS

Paclitaxel disposition was well-described for all formulations. The derived model predicts toxicodynamics of diverse paclitaxel carriers.

Lifespan based indirect response models

In the field of hematology, several mechanism-based pharmacokinetic-pharmacodynamic models have been developed to understand the dynamics of several blood cell populations under different clinical conditions while accounting for the essential underlying principles of pharmacology, physiology and pathology. In general, a population of blood cells is basically controlled by two processes: the cell production and cell loss. The assumption that each cell exits the population when its lifespan expires implies that the cell loss rate is equal to the cell production rate delayed by the lifespan and justifies the use of delayed differential equations for compartmental modeling. This review is focused on lifespan models based on delayed differential equations and presents the structure and properties of the basic lifespan indirect response (LIDR) models for drugs affecting cell production or cell lifespan distribution. The LIDR models for drugs affecting the precursor cell production or decreasing the precursor cell population are also presented and their properties are discussed. The interpretation of transit compartment models as LIDR models is reviewed as the basis for introducing a new LIDR for drugs affecting the cell lifespan distribution. Finally, the applications and limitations of the LIDR models are discussed.

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).

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.

Maitake beta-glucan promotes recovery of leukocytes and myeloid cell function in peripheral blood from paclitaxel hematotoxicity

Bone marrow myelotoxicity is a major limitation of chemotherapy. While granulocyte colony stimulating factor (G-CSF) treatment is effective, alternative approaches to support hematopoietic recovery are sought. We previously found that a beta-glucan extract from maitake mushroom Grifola frondosa (MBG) enhanced colony forming unit-granulocyte monocyte (CFU-GM) activity of mouse bone marrow and human hematopoietic progenitor cells (HPC), stimulated G-CSF production and spared HPC from doxorubicin toxicity in vitro. This investigation assessed the effects of MBG on leukocyte recovery and granulocyte/monocyte function in vivo after dose intensive paclitaxel (Ptx) in a normal mouse. After a cumulative dose of Ptx (90-120 mg/kg) given to B6D2F1mice, daily oral MBG (4 or 6 mg/kg), intravenous G-CSF (80 microg/kg) or Ptx alone were compared for effects on the dynamics of leukocyte recovery in blood, CFU-GM activity in bone marrow and spleen, and granulocyte/monocyte production of reactive oxygen species (ROS). Leukocyte counts declined less in Ptx + MBG mice compared to Ptx-alone (p = 0.024) or Ptx + G-CSF treatment (p = 0.031). Lymphocyte levels were higher after Ptx + MBG but not Ptx + G-CSF treatment compared to Ptx alone (p < 0.01). MBG increased CFU-GM activity in bone marrow and spleen (p < 0.001, p = 0.002) 2 days after Ptx. After two additional days (Ptx post-day 4), MBG restored granulocyte/monocyte ROS response to normal levels compared to Ptx-alone and increased ROS response compared to Ptx-alone or Ptx + G-CSF (p < 0.01, both). The studies indicate that oral MBG promoted maturation of HPC to become functionally active myeloid cells and enhanced peripheral blood leukocyte recovery after chemotoxic bone marrow injury.

Differential pharmacodynamic effects of paclitaxel formulations in an intracranial rat brain tumor model

Nano- and microparticulate carriers can exert a beneficial impact on the pharmacodynamics of anticancer agents. To investigate the relationships between carrier and antitumor pharmacodynamics, paclitaxel incorporated in liposomes (L-pac) was compared with the clinical standard formulated in Cremophor-EL/ethanol (Cre-pac) in a rat model of advanced primary brain cancer. Three maximum-tolerated-dose regimens given by intravenous administration were investigated: 50 mg/kg on day 8 (d8) after implantation of 9L gliosarcoma tumors; 40 mg/kg on d8 and d15; 20 mg/kg on d8, d11, and d15. Body weight change and neutropenia were assessed as pharmacodynamic markers of toxicity. The pharmacodynamic markers of antitumor efficacy were increase in lifespan (ILS) and tumor volume progression, measured noninvasively by magnetic resonance imaging. At equivalent doses, neutropenia was similar for both formulations, but weight loss was more severe for Cre-pac. No regimen of Cre-pac extended survival, whereas L-pac at 40 mg/kg x2 doses was well tolerated and mediated 26% ILS (p < 0.0002) compared with controls. L-pac at a lower cumulative dose (20 mg/kg x3) was even more effective (40% ILS; p < 0.0001). In striking contrast, the identical regimen of Cre-pac was lethal. Development of a novel semimechanistic pharmacodynamic model permitted quantitative hypothesis testing with the tumor volume progression data, and suggested the existence of a transient treatment effect that was consistent with sensitization or "priming" of tumors by more frequent L-pac dosing schedules. Therefore, improved antitumor responses of carrier-based paclitaxel formulations can arise both from dose escalation, because of reduced toxicity, and from novel carrier-mediated alterations of antitumor pharmacodynamic effects.