Population modeling of filgrastim PK-PD in healthy adults following intravenous and subcutaneous administrations

Semimechanistic pharmacokinetic-pharmacodynamic model of tripegfilgrastim for pediatric patients after chemotherapy

Tripegfilgrastim is a long-acting granulocyte colony-stimulating factor (G-CSF) that has been used to prevent chemotherapy-induced neutropenia in adults. This study aimed to establish a pharmacokinetic (PK)-pharmacodynamic (PD) model to explore the impact of chemotherapy and tripegfilgrastim on absolute neutrophil counts (ANCs) and to further propose a fixed-dose regimen in pediatric patients. Because neutrophils affect the clearance of tripegfilgrastim, the semimechanistic PK-PD model was developed simultaneously by using data from 40 healthy adults and 27 pediatric patients with solid tumors. Tripegfilgrastim PK and ANC dynamics were described with a pharmacodynamics-mediated drug disposition model assuming quasi-equilibrium with five transit compartments mimicking neutrophil granulopoiesis. The effect of chemotherapy on neutrophils was included by stimulating the elimination of the G-CSF receptor at the mitotic cells. Healthy adult and pediatric patients showed significantly different value for dissociation constant of the tripegfilgrastim-G-CSF receptor complex (Kd ) and apparent volume of distribution (Vd /F). Patients treated with chemotherapy had a higher Vd /F and 62% lower Kd than healthy adults. As the age increased, the absorption rate of tripegfilgrastim was decreased. Body weight affected the G-CSF receptor-mediated internalization of tripegfilgrastim, and the baseline ANC value impacted the production rate of G-CSF receptors. Simulations from the developed model suggested that 1.5, 2.5, 4, and 6 mg single subcutaneous tripegfilgrastim doses for the respective weight groups of 10-20, 21-30, 31-44, and more than 45 kg significantly reduced the duration of Grade 4 neutropenia similar to tripegfilgrastim weight-based treatment with 100 μg/kg.

A Pharmacometric Model to Predict Chemotherapy-Induced Myelosuppression and Associated Risk Factors in Non-Small Cell Lung Cancer

Chemotherapy often induces severe neutropenia due to the myelosuppressive effect. While predictive pharmacokinetic (PK)/pharmacodynamic (PD) models of absolute neutrophil count (ANC) after anticancer drug administrations have been developed, their deployments to routine clinics have been limited due to the unavailability of PK data and sparseness of PD (or ANC) data. Here, we sought to develop a model describing temporal changes of ANC in non-small cell lung cancer patients receiving (i) combined chemotherapy of paclitaxel and cisplatin and (ii) granulocyte colony stimulating factor (G-CSF) treatment when needed, under such limited circumstances. Maturation of myelocytes into blood neutrophils was described by transit compartments with negative feedback. The K-PD model was employed for drug effects with drug concentration unavailable and the constant model for G-CSF effects. The fitted model exhibited reasonable goodness of fit and parameter estimates. Covariate analyses revealed that ANC decreased in those without diabetes mellitus and female patients. Using the final model obtained, an R Shiny web-based application was developed, which can visualize predicted ANC profiles and associated risk of severe neutropenia for a new patient. Our model and application can be used as a supportive tool to identify patients at the risk of grade 4 neutropenia early and suggest dose reduction.

A prospective study of filgrastim pharmacokinetics in morbidly obese patients compared with non-obese controls

INTRODUCTION

Filgrastim is a human granulocyte colony-stimulating factor (G-CSF). There are limited data on dosing filgrastim in obesity. The objective of this study was to compare filgrastim pharmacokinetic parameters for morbidly obese and non-obese patients after a single subcutaneous dose of filgrastim dosed per actual body weight.

METHODS

This prospective, matched-pair study (NCT01719432) included patients ≥18 years of age, receiving filgrastim at 5 μg/kg with a weight >190% of their ideal body weight (IBW) for "morbidly obese" patients or within 80%-124% of IBW for matched-control patients. The control group was prospectively matched for age (within 10 years), degree of neutropenia, and gender. Filgrastim doses were not rounded to vial size, to allow more accurate assessment of exposure. Blood samples were collected at 0 (prior to dose), 2, 4, 6, 8, 12, and 24 h after the first subcutaneous administration of filgrastim.

RESULTS

A total of 30 patients were enrolled in this prospective pharmacokinetic study, with 15 patients assigned to each arm. Non-compartmental analysis showed that the systemic clearance (Cl) was 0.111 ± 0.041 ml/min in the morbidly obese group versus 0.124 ± 0.045 ml/min in the non-obese group (p = 0.44). Additionally, the mean area under the curve (AUC_0-24h ) was 49.3 ± 13.9 ng/ml × min in the morbidly obese group versus 46.3 ± 16.8 ng/mL x min in the non-obese group (p = 0.6). No differences were seen in maximum concentrations (C_max ) between the two groups (morbidly obese: 48.1 ± 14.7 ng/ml vs. non-obese: 49.2 ± 20.7 ng/ml (p = 0.87)). The morbidly obese group had a numerically higher, but not statistically significant, increase in time to maximum concentration (T_max ) compared to the non-obese group (544 ± 145 min vs 436 ± 156 min (p = 0.06), respectively).

CONCLUSION

Calculating subcutaneous filgrastim doses using actual body weight appears to produce similar systemic exposure in morbidly obese and non-obese patients with severe neutropenia.

A Quantitative Systems Pharmacology Framework for Optimal Doxorubicin Granulocyte Colony-Stimulating Factor Regimens in Triple-Negative Breast Cancer

INTRODUCTION

To mitigate the risk of neutropenia during chemotherapy treatment of triple-negative breast cancer, prophylactic and supportive therapy with granulocyte colony-stimulating factor (G-CSF) is administered concomitant to chemotherapy. The proper timing of combined chemotherapy and G-CSF is crucial for treatment outcomes.

METHODS

Leveraging our established mathematical model of neutrophil production by G-CSF, we developed quantitative systems pharmacology (QSP) framework to investigate how modulating chemotherapy dose frequency and intensity can maximize antitumour effects. To establish schedules that best control tumour size while minimizing neutropenia, we combined Gompertzian tumour growth with pharmacokinetic/pharmacodynamic models of doxorubicin and G-CSF, and our QSP model of neutrophil production.

RESULTS

We optimized a range of chemotherapeutic cycle lengths and dose sizes to establish regimens that simultaneously reduced tumour burden while minimizing neutropenia. Our results suggest that cytotoxic chemotherapy with doxorubicin 45 mg/m2 every 14 days provides effective control of tumour growth while mitigating neutropenic risks.

CONCLUSION

This work suggests future avenues for optimal regimens of chemotherapy with prophylactic G-CSF support. Importantly, the algorithmic approach that we developed can aid in balancing the anticancer and the neutropenic effects of both drugs, and therefore contributes to rational considerations in clinical decision-making in triple-negative breast cancer.

A whole-body circulatory neutrophil model with application to predicting clinical neutropenia from in vitro studies

A circulatory model of granulopoiesis and its regulation is presented that includes neutrophil trafficking in the lungs, liver, spleen, bone marrow, lymph nodes, and blood. In each organ, neutrophils undergo transendothelial migration from vascular to interstitial space, clearance due to apoptosis, and recycling via the lymphatic flow. The model includes cell cycling of progenitor cells in the bone marrow, granulocyte colony-stimulating factor (G-CSF) kinetics and its neutrophil regulatory action, as well as neutrophil margination in the blood. From previously reported studies, 111 In-labeled neutrophil kinetic data in the blood and sampled organs were used to estimate the organ trafficking parameters in the model. The model was further developed and evaluated using absolute neutrophil count (ANC), band cell, and segmented neutrophil time course data from healthy volunteers following four dose levels of pegfilgrastim (r2  = 0.77-0.99), along with ANC time course responses following filgrastim (r2  = 0.96). The baseline values of various cell types in bone marrow and blood, as well as G-CSF concentration in the blood, predicted by the model are consistent with available literature reports. After incorporating the mechanism of action of both paclitaxel and carboplatin, as determined from an in vitro bone marrow studies, the model reliably predicted the observed ANC time course following paclitaxel plus carboplatin observed in a phase I trial of 46 patients (r2  = 0.70). The circulatory neutrophil model may provide a mechanistic framework for predicting multi-organ neutrophil homeostasis and dynamics in response to therapeutic agents that target neutrophil dynamics and trafficking in different organs.

Continuous Intravenous Administration of Granulocyte-Colony-Stimulating Factors-A Breakthrough in the Treatment of Cancer Patients with Febrile Neutropenia

Background: Febrile neutropenia (FN) remains one of the most challenging problems in medical oncology and is a very severe side effect of chemotherapy. Its late consequences, when it is recurrent or of a severe grade, are dose reduction and therapy delays. Current guidelines allow the administration of granulocyte-colony-stimulating factors (G-CSF) for profound FN (except for the case when a pegylated form of G-CSF is administrated with prophylactic intention) in addition to antibiotics and supportive care. Methods: This is a prospective study that included 96 patients with confirmed malignancy, treated with chemotherapy, who developed FN during their oncological therapy, and were hospitalized. They received standard treatment plus a dose of G-CSF of 16 µg/Kg/day IV continuous infusion. Results: The gender distribution was almost symmetrical: Male patients made up 48.96% and 51.04% were female patients, with no significance on recovery from FN (p = 1.00). The patients who received prophylactic G-CSF made up 20.21%, but this was not a predictive or prognostic factor for the recovery time from aplasia (p = 0.34). The median chemotherapy line where patients with FN were included was two and the number of previous chemotherapy cycles before FN was three. The median serological number of neutrophils (PMN) was 450/mm³ and leucocytes (WBC) 1875/mm³ at the time of FN. Ten patients possess PMN less than 100/mm³. The median time to recovery was 25.5 h for 96 included patients, with one failure in which the patient possessed grade 5 FN. Predictive factors for shorter recovery time were lower levels of C reactive protein (p < 0.001) and procalcitonin (p = 0.002) upon hospital admission and higher WBC (p = 0.006) and PMN (p < 0.001) at the time of the provoking cycle of chemotherapy for FN. The best chance for a shorter duration of FN was a short history of chemotherapy regarding the number of cycles) (p < 0.0001). Conclusions: Continuous IV administration of G-CSF could be an alternative salvage treatment for patients with profound febrile neutropenia, with a very fast recovery time for neutrophiles.

A mathematical model of viral oncology as an immuno-oncology instigator

We develop and analyse a mathematical model of tumour-immune interaction that explicitly incorporates heterogeneity in tumour cell cycle duration by using a distributed delay differential equation. We derive a necessary and sufficient condition for local stability of the cancer-free equilibrium in which the amount of tumour-immune interaction completely characterizes disease progression. Consistent with the immunoediting hypothesis, we show that decreasing tumour-immune interaction leads to tumour expansion. Finally, by simulating the mathematical model, we show that the strength of tumour-immune interaction determines the long-term success or failure of viral therapy.

Predicting Chemotherapy-Induced Neutropenia and Granulocyte Colony-Stimulating Factor Response Using Model-Based In Vitro to Clinical Translation

The ability to predict the incidence of chemotherapy-induced neutropenia in early drug development can inform risk monitoring and mitigation strategies, as well as decisions on advancing compounds to clinical trials. In this report, a physiological model of granulopoiesis that incorporates the drug's mechanism of action on cell cycle proliferation of bone marrow progenitor cells was extended to include the action of the cytotoxic agents paclitaxel, carboplatin, doxorubicin, and gemcitabine. In vitro bone marrow studies were conducted with each compound, and results were used to determine the model's drug effect parameters. Population simulations were performed to predict the absolute neutrophil count (ANC) and incidence of neutropenia for each compound, which were compared to results reported in the literature. In addition, using the single agent in vitro study results, the model was able to predict ANC time course in response to paclitaxel plus carboplatin in combination, which compared favorably to the results reported in a phase 1 clinical trial of 46 patients (r2 = 0.70). Model simulations were used to compare the relative risk (RR) of neutropenia in patients with high baseline ANCs for five chemotherapeutic regimens: doxorubicin (RR = 0.59), paclitaxel plus carboplatin combination (RR = 0.079), carboplatin (RR = 0.047), paclitaxel (RR = 0.031), and gemcitabine (RR = 0.013). Finally, the model was applied to quantify the reduced incidence of neutropenia with coadministration of pegfilgrastim or filgrastim, for both paclitaxel and the combination of paclitaxel plus carboplatin. The model provides a framework for predicting clinical neutropenia using in vitro bone marrow studies of anticancer agents that may guide drug development decisions.

Impact of granulocyte colony-stimulating factor on FOLFIRINOX-induced neutropenia prevention: A population pharmacokinetic/pharmacodynamic approach

AIMS

Granulocyte colony-stimulating factor (G-CSF) is frequently prescribed to prevent chemotherapy-induced neutropenia, but the administration schedule remains empirical in case of bimonthly chemotherapy such as FOLFIRINOX regimen. This pharmacokinetic/pharmacodynamic (PK/PD) study was performed to determine the effect of different G-CSF regimens on the incidence and duration of neutropenia following FOLFIRINOX administration in order to propose an optimal G-CSF dosing schedule.

METHODS

A population PK/PD model was developed to describe individual neutrophil time course from absolute neutrophil counts (ANC) obtained in 40 advanced cancer patients receiving FOLFIRINOX regimen. The structural model considered ANC dynamics, neutropenic effect of cytotoxics and the stimulating effect of G-CSF on neutrophils. Final model estimates were used to simulate different G-CSF dosing schedules for 1000 virtual subjects. The incidence and duration of neutropenia were then calculated for different G-CSF dosing schedules.

RESULTS

The final model successfully described the myelosuppressive effect induced by the 3 cytotoxics for all patients. Simulations showed that pegfilgrastim administration reduced the risk of severe neutropenia by 22.9% for subjects with low ANC at the start of chemotherapy. Median duration in this group was also shortened by 3.1 days when compared to absence of G-CSF. Delayed G-CSF administration was responsible for higher incidence and longer duration of neutropenia compared to absence of administration.

CONCLUSION

The PK/PD model well described our population's ANC data. Simulations showed that pegylated-G-CSF administration 24 hours after the end of chemotherapy seems to be the optimal schedule to reduce FOLFIRINOX-induced neutropenia. We also underline the potential negative effect of G-CSF maladministration.

Model-Based Optimal AML Consolidation Treatment

OBJECTIVE

Neutropenia is an adverse event commonly arising during intensive chemotherapy of acute myeloid leukemia (AML). It is often associated with infectious complications. Mathematical modeling, simulation, and optimization of the treatment process would be a valuable tool to support clinical decision making, potentially resulting in less severe side effects and deeper remissions. However, until now, there has been no validated mathematical model available to simulate the effect of chemotherapy treatment on white blood cell (WBC) counts and leukemic cells simultaneously.

METHODS

We developed a population pharmacokinetic/pharmacodynamic (PK/PD) model combining a myelosuppression model considering endogenous granulocyte-colony stimulating factor (G-CSF), a PK model for cytarabine (Ara-C), a subcutaneous absorption model for exogenous G-CSF, and a two-compartment model for leukemic blasts. This model was fitted to data of 44 AML patients during consolidation therapy with a novel Ara-C plus G-CSF schedule from a phase II controlled clinical trial. Additionally, we were able to optimize treatment schedules with respect to disease progression, WBC nadirs, and the amount of Ara-C and G-CSF.

RESULTS

The developed PK/PD model provided good prediction accuracies and an interpretation of the interaction between WBCs, G-CSF, and blasts. For 14 patients (those with available bone marrow blast counts), we achieved a median 4.2-fold higher WBC count at nadir, which is the most critical time during consolidation therapy. The simulation results showed that relative bone marrow blast counts remained below the clinically important threshold of 5%, with a median of 60% reduction in Ara-C.

CONCLUSION

These in silico findings demonstrate the benefits of optimized treatment schedules for AML patients.

SIGNIFICANCE

Until 2017, no new drug had been approved for the treatment of AML, fostering the optimal use of currently available drugs.

A physiological model of granulopoiesis to predict clinical drug induced neutropenia from in vitro bone marrow studies: with application to a cell cycle inhibitor

Neutropenia is one of the most common dose-limiting toxocities associated with anticancer drug therapy. The ability to predict the probability and severity of neutropenia based on in vitro studies of drugs in early drug development will aid in advancing safe and efficacious compounds to human testing. Toward this end, a physiological model of granulopoiesis and its regulation is presented that includes the bone marrow progenitor cell cycle, allowing for a mechanistic representation of the action of relevant anticancer drugs based on in vitro studies. Model development used data from previously reported tracer kinetic studies of granulocyte disposition in healthy humans to characterize the dynamics of neutrophil margination in the presence of endogenous granulocyte-colony stimulating factor (G-CSF). In addition, previously published data from healthy volunteers following pegfilgrastim and filgrastim were used to quantify the regulatory effects of support G-CSF therapies on granulopoiesis. The model was evaluated for the cell cycle inhibitor palbociclib, using an in vitro system of human bone marrow mononuclear cells to quantify the action of palbociclib on proliferating progenitor cells, including its inhibitory effect on G1 to S phase transition. The in vitro results were incorporated into the physiological model of granulopoiesis and used to predict the time course of absolute neutrophil count (ANC) and the incidence of neutropenia observed in three previously reported clinical trials of palbociclib. The model was able to predict grade 3 and 4 neutropenia due to palbociclib treatment with 86% accuracy based on in vitro data.

Mathematical analysis and drug exposure evaluation of pharmacokinetic models with endogenous production and simultaneous first-order and Michaelis-Menten elimination: the case of single dose

Drugs with an additional endogenous source often exhibit simultaneous first-order and Michaelis-Menten elimination and are becoming quite common in pharmacokinetic modeling. In this paper, we investigate the case of single dose intravenous bolus administration for the one-compartment model. Relying on a formerly introduced transcendent function, we were able to analytically express the concentration time course of this model and provide the pharmacokinetic interpretation of its components. Using the concept of the corrected concentration, the mathematical expressions for the partial and total areas under the concentration time curve (AUC) were also given. The impact on the corrected concentration and AUC is discussed as well as the relative contribution of the exogenous part in presence of endogenous production. The present findings theoretically elucidate several pharmacokinetic issues for the considered drug compounds and provide guidance for the rational estimation of their pharmacokinetic parameters.

Pharmacokinetic-pharmacodynamic modelling of neutrophil response to G-CSF in healthy subjects and patients with chemotherapy-induced neutropenia

AIM

The objective of the present study was to use pharmacokinetic-pharmacodynamic modelling to characterize the effects of chemotherapy on the granulopoietic system and to predict the absolute neutrophil counts (ANCs) for patients with chemotherapy-induced neutropenia treated with filgrastim and pegfilgrastim.

METHODS

Data were extracted from 10 phase I-III studies conducted in 110 healthy adults, and 618 adult and 52 paediatric patients on chemotherapy following administration of filgrastim or pegfilgrastim. The structural model accounted for ANC dynamics and the effects of filgrastim and pegfilgrastim, chemotherapy and corticosteroids. The impact of neutrophils on drug disposition was based on a drug receptor-binding model that assumed quasi-equilibrium and stimulation of the production and maturation of neutrophils upon treatment. The chemotherapy and corticosteroid effects were represented by kinetic-pharmacodynamic-type models, where chemotherapy stimulated elimination of neutrophil precursors at the mitotic stage, and corticosteroids stimulated neutrophil production.

RESULTS

The systemic half-lives of filgrastim (2.6 h) and pegfilgrastim (10.1 h) were as expected. The effective half-life of chemotherapy was 9.6 h, with a 2-day killing effect. The rate of receptor elimination from mitotic compartments exhibited extreme interindividual variability (% coefficient of variation >200), suggesting marked differences in sensitivity to chemotherapy effects on ANCs. The stimulatory effects of pegfilgrastim were significantly greater than those of filgrastim. Model qualification confirmed the predictive capability of this model.

CONCLUSION

This qualified model simulates the time course of ANC in the absence or presence of chemotherapy and predicts nadir, time to nadir and time of recovery from different grades of neutropenia upon treatment with filgrastim and pegfilgrastim.

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.

AFIR: A Dimensionless Potency Metric for Characterizing the Activity of Monoclonal Antibodies

For monoclonal antibody (mAb) drugs, soluble targets may accumulate several thousand fold after binding to the drug. Time course data of mAb and total target is often collected and, although free target is more closely related to clinical effect, it is difficult to measure. Therefore, mathematical models of this data are used to predict target engagement. In this article, a “potency factor” is introduced as an approximation for the model‐predicted target inhibition. This potency factor is defined to be the time‐Averaged Free target concentration to Initial target concentration Ratio (AFIR), and it depends on three key quantities: the average drug concentration at steady state; the binding affinity; and the degree of target accumulation. AFIR provides the intuition for how changes in dosing regimen and binding affinity affect target capture and AFIR can be used to predict the druggability of new targets and the expected benefits of more potent, second‐generation mAbs.

A cell-level model of pharmacodynamics-mediated drug disposition

We aimed to develop a cell-level pharmacodynamics-mediated drug disposition (PDMDD) model to analyze in vivo systems where the PD response to a drug has an appreciable effect on the pharmacokinetics (PK). An existing cellular level model of PD stimulation was combined with the standard target-mediated drug disposition (TMDD) model and the resulting model structure was parametrically identifiable from typical in vivo PK and PD data. The PD model of the cell population was controlled by the production rate k in and elimination rate k out which could be stimulated or inhibited by the number of bound receptors on a single cell. Simulations were performed to assess the impact of single and repeated dosing on the total drug clearance. The clinical utility of the cell-level PDMDD model was demonstrated by fitting published data on the stimulatory effects of filgrastim on absolute neutrophil counts in healthy subjects. We postulated repeated dosing as a means of detecting and quantifying PDMDD as a single dose might not be sufficient to elicit the cellular response capable of altering the receptor pool to visibly affect drug disposition. In the absence of any PD effect, the model reduces down to the standard TMDD model. The applications of this model can be readily extended to include chemotherapy-induced cytopenias affecting clearance of endogenous hematopoietic growth factors, different monoclonal antibodies and immunogenicity effects on PK.

A Mathematical Model of Granulopoiesis Incorporating the Negative Feedback Dynamics and Kinetics of G-CSF/Neutrophil Binding and Internalization

We develop a physiological model of granulopoiesis which includes explicit modelling of the kinetics of the cytokine granulocyte colony-stimulating factor (G-CSF) incorporating both the freely circulating concentration and the concentration of the cytokine bound to mature neutrophils. G-CSF concentrations are used to directly regulate neutrophil production, with the rate of differentiation of stem cells to neutrophil precursors, the effective proliferation rate in mitosis, the maturation time, and the release rate from the mature marrow reservoir into circulation all dependent on the level of G-CSF in the system. The dependence of the maturation time on the cytokine concentration introduces a state-dependent delay into our differential equation model, and we show how this is derived from an age-structured partial differential equation model of the mitosis and maturation and also detail the derivation of the rest of our model. The model and its estimated parameters are shown to successfully predict the neutrophil and G-CSF responses to a variety of treatment scenarios, including the combined administration of chemotherapy and exogenous G-CSF. This concomitant treatment was reproduced without any additional fitting to characterize drug-drug interactions.

Neutrophil dynamics during concurrent chemotherapy and G-CSF administration: Mathematical modelling guides dose optimisation to minimise neutropenia

The choice of chemotherapy regimens is often constrained by the patient's tolerance to the side effects of chemotherapeutic agents. This dose-limiting issue is a major concern in dose regimen design, which is typically focused on maximising drug benefits. Chemotherapy-induced neutropenia is one of the most prevalent toxic effects patients experience and frequently threatens the efficient use of chemotherapy. In response, granulocyte colony-stimulating factor (G-CSF) is co-administered during chemotherapy to stimulate neutrophil production, increase neutrophil counts, and hopefully avoid neutropenia. Its clinical use is, however, largely dictated by trial and error processes. Based on up-to-date knowledge and rational considerations, we develop a physiologically realistic model to mathematically characterise the neutrophil production in the bone marrow which we then integrate with pharmacokinetic and pharmacodynamic (PKPD) models of a chemotherapeutic agent and an exogenous form of G-CSF (recombinant human G-CSF, or rhG-CSF). In this work, model parameters represent the average values for a general patient and are extracted from the literature or estimated from available data. The dose effect predicted by the model is confirmed through previously published data. Using our model, we were able to determine clinically relevant dosing regimens that advantageously reduce the number of rhG-CSF administrations compared to original studies while significantly improving the neutropenia status. More particularly, we determine that it could be beneficial to delay the first administration of rhG-CSF to day seven post-chemotherapy and reduce the number of administrations from ten to three or four for a patient undergoing 14-day periodic chemotherapy.

Closed form solutions and dominant elimination pathways of simultaneous first-order and Michaelis-Menten kinetics

The current study aims to provide the closed form solutions of one-compartment open models exhibiting simultaneous linear and nonlinear Michaelis-Menten elimination kinetics for single- and multiple-dose intravenous bolus administrations. It can be shown that the elimination half-time ([Formula: see text]) has a dose-dependent property and is upper-bounded by [Formula: see text] of the first-order elimination model. We further analytically distinguish the dominant role of different elimination pathways in terms of model parameters. Moreover, for the case of multiple-dose intravenous bolus administration, the existence and local stability of the periodic solution at steady state are established. The closed form solutions of the models are obtained through a newly introduced function motivated by the Lambert W function.

Pharmacokinetics and pharmacodynamics of FSK0808 and Gran after single intravenous drip administration or single subcutaneous administration: comparative study in healthy Japanese adult male subjects

FSK0808 is a recombinant human granulocyte colony-stimulating factor developed by Fuji Pharma Co., Ltd and Mochida Pharmaceutical Co., Ltd. as a biosimilar product of Gran®. We verified the pharmacokinetic/pharmacodynamic equivalence of FSK0808 and commercially available Gran® by a randomized crossover study of single intravenous dose (200 µg/m(2)) and single subcutaneous dose (400 µg/m(2)) in healthy Japanese adult male subjects. According to the bioequivalence guidelines, the area under the blood concentration - time curve by 48 hours after administration (AUC0-48) in a single intravenous drip (IVD) study, and AUC0-48 and maximum blood concentration (Cmax) in a single subcutaneous (SC) dose study were used as primary endpoints, and the pharmacodynamic parameters including absolute neutrophil count (ANC) or number of CD34 positive cells (CD34(+) cells) as secondary endpoints. The safety was evaluated based on the characteristics and incidence of adverse reactions. As a result, the 90% confidence interval (CI) of the difference in mean value for AUC0-48 among drugs ranged from log(0.8) to log(1.25), in the IVD study, and those for Cmax and AUC0-48 were within the range of log(0.8)-log(1.25) in the SC study. Those for secondary endpoints were all within the range of log(0.8)-log(1.25). Thus, the pharmacokinetics/pharmacodynamics of both drugs were considered equivalent for all routes of administration, and the profiles of adverse reactions were also very similar.

Basic concepts in population modeling, simulation, and model-based drug development: part 3-introduction to pharmacodynamic modeling methods

Population pharmacodynamic (PD) models describe the time course of drug effects, relating exposure to response, and providing a more robust understanding of drug action than single assessments. PD models can test alternative dose regimens through simulation, allowing for informed assessment of potential dose regimens and study designs. This is the third paper in a three-part series, providing an introduction into methods for developing and evaluating population PD models. Example files are available in the Supplementary Data.

T-lymphocyte perturbation following large-scale apheresis and hematopoietic stem cell transplantation in HIV-infected individuals

Analysis and mathematical modeling of T-lymphocyte perturbation following administration of granulocyte colony stimulating factor (G-CSF) and two large-scale aphereses are reported. 74 HIV-1 positive antiretroviral-treated individuals were infused with gene- or sham-transduced CD34+ hematopoietic stem cells (HSC) in a Phase II clinical trial. T cell numbers were examined in four phases: 1) during steady state; 2) increases in peripheral blood (PB) following G-CSF administration; 3) depletion post-aphereses and 4) reconstitution post HSC infusion. The present analysis provides the first direct estimate of CD4+ T cell distribution and trafficking in HIV-infected individuals on stable HAART, indicating that CD4+ T lymphocytes in PB represent 5.5% of the pool of CD4+ T lymphocytes that traffic to PB.

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.

Population pharmacokinetic and pharmacodynamic model-based comparability assessment of a recombinant human Epoetin Alfa and the Biosimilar HX575

The aim of this study was to develop an integrated pharmacokinetic and pharmacodynamic (PK/PD) model and assess the comparability between epoetin alfa HEXAL/Binocrit (HX575) and a comparator epoetin alfa by a model-based approach. PK/PD data-including serum drug concentrations, reticulocyte counts, red blood cells, and hemoglobin levels-were obtained from 2 clinical studies. In sum, 149 healthy men received multiple intravenous or subcutaneous doses of HX575 (100 IU/kg) and the comparator 3 times a week for 4 weeks. A population model based on pharmacodynamics-mediated drug disposition and cell maturation processes was used to characterize the PK/PD data for the 2 drugs. Simulations showed that due to target amount changes, total clearance may increase up to 2.4-fold as compared with the baseline. Further simulations suggested that once-weekly and thrice-weekly subcutaneous dosing regimens would result in similar efficacy. The findings from the model-based analysis were consistent with previous results using the standard noncompartmental approach demonstrating PK/PD comparability between HX575 and comparator. However, due to complexity of the PK/PD model, control of random effects was not straightforward. Whereas population PK/PD model-based analyses are suited for studying complex biological systems, such models have their limitations (statistical), and their comparability results should be interpreted carefully.

Theoretical considerations of target-mediated drug disposition models: simplifications and approximations

PURPOSE

To clarify relationships among various types of target-mediated disposition (TMD) models including the Michaelis-Menten, quasi-steady-state (Qss), and rapid binding models and propose measures for the closeness of some models as approximations to the general TMD model (Mager and Jusko, J Pharmacokinet Pharmacodyn 28(6):507-532, 2001).

METHODS

Based on the classic singular perturbation theory by selecting appropriate scales of time, we derive requirements with which the Michaelis-Menten and Qss models are suitable approximations. Under the Qss assumption we show that other simplifications of the general TMD model can be similarly obtained as the Michaelis-Menten and Qss models. We compare these models by simulations using known application examples.

RESULTS

The Michaelis-Menten and Qss models are direct simplifications of the general TMD model and, moreover, suitable approximations if certain specific requirements on the parameters are met.

CONCLUSIONS

As a first attempt to quantify the closeness of some simplifications to the general TMD model, our work should provide a more rigorous basis for the theoretical and practical research of TMD models, which are important for investigating the pharmacokinetic-pharmacodynamic relationships of many biological compounds.

Interpretation of transit compartments pharmacodynamic models as lifespan based indirect response models

Transit compartments (TC) models are used to describe pharmacodynamic responses that involve drug action on cells undergoing differentiation and maturation. Such pharmacodynamic systems can also be described by lifespan based indirect response (LIDR) models. The purpose of this report is to investigate conditions under which the transit compartments models can be considered a special case of LIDR models. An integral representation of a solution to TC model has been used to determine the lifespan distribution for cell population described by this model. The distribution served as a basis for definition of new LIDRE (lifespan based indirect response with an effect on the lifespan distribution) models. Time courses of responses described by both types of models were simulated for a monoexponential pharmacokinetic function. The limit response was calculated as the number of transit compartments approached infinity. The difference between the limit response and TC responses were evaluated by computer simulations using MATLAB 7.7. TC models are a special case of LIDR models with the lifespan distribution described by the gamma function. If drug affects only the production of cells, then the cell lifespan distribution is time invariant. In this case an increase in the number of compartments results in a basic LIDR model with a point lifespan distribution. When the drug inhibits or stimulates cell aging, the cell lifespan distribution becomes time dependent revealing a new mechanism for drug effect on the gamma probability density function. The TC model with a large number of transit compartments converges to an LIDRE model. The limit LIDR models are approximated by the TC models when the number of compartments is at least 5. A moderate improvement in the approximation is observed if this number exceeds 20. The lifespan distribution for a cell population described by a TC model is described by the gamma probability density function. A drug affects this distribution only if it stimulates or inhibits the rate of cell maturation. If the number of transit compartments increases, then the TC model converges to a new type of LIDR model.