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

A framework for precision "dosing" of mental healthcare services: algorithm development and clinical pilot

BACKGROUND

One in five adults in the US experience mental illness and over half of these adults do not receive treatment. In addition to the access gap, few innovations have been reported for ensuring the right level of mental healthcare service is available at the right time for individual patients.

METHODS

Historical observational clinical data was leveraged from a virtual healthcare system. We conceptualize mental healthcare services themselves as therapeutic interventions and develop a prototype computational framework to estimate their potential longitudinal impacts on depressive symptom severity, which is then used to assess new treatment schedules and delivered to clinicians via a dashboard. We operationally define this process as "session dosing": 497 patients who started treatment with severe symptoms of depression between November 2020 and October 2021 were used for modeling. Subsequently, 22 mental health providers participated in a 5-week clinical quality improvement (QI) pilot, where they utilized the prototype dashboard in treatment planning with 126 patients.

RESULTS

The developed framework was able to resolve patient symptom fluctuations from their treatment schedules: 77% of the modeling dataset fit criteria for using the individual fits for subsequent clinical planning where five anecdotal profile types were identified that presented different clinical opportunities. Based on initial quality thresholds for model fits, 88% of those individuals were identified as adequate for session optimization planning using the developed dashboard, while 12% supported more thorough treatment planning (e.g. different treatment modalities). In the clinical pilot, 90% of clinicians reported using the dashboard a few times or more per member. Although most clinicians (67.5%) either rarely or never used the dashboard to change session types, numerous other discussions were enabled, and opportunities for automating session recommendations were identified.

CONCLUSIONS

It is possible to model and identify the extent to which mental healthcare services can resolve depressive symptom severity fluctuations. Implementation of one such prototype framework in a real-world clinic represents an advancement in mental healthcare treatment planning; however, investigations to assess which clinical endpoints are impacted by this technology, and the best way to incorporate such frameworks into clinical workflows, are needed and are actively being pursued.

A comprehensive regulatory and industry review of modeling and simulation practices in oncology clinical drug development

Exposure-response (E-R) analyses are an integral component in the development of oncology products. Characterizing the relationship between drug exposure metrics and response allows the sponsor to use modeling and simulation to address both internal and external drug development questions (e.g., optimal dose, frequency of administration, dose adjustments for special populations). This white paper is the output of an industry-government collaboration among scientists with broad experience in E-R modeling as part of regulatory submissions. The goal of this white paper is to provide guidance on what the preferred methods for E-R analysis in oncology clinical drug development are and what metrics of exposure should be considered.

Leveraging modeling and simulation to optimize the therapeutic window for epigenetic modifier drugs

Dysregulated epigenetic processes can lead to altered gene expression and give rise to malignant transformation and tumorigenesis. Epigenetic drugs aim to revert the phenotype of cancer cells to normally functioning cells, and are developed and applied to treat both hematological and solid cancers. Despite this promising therapeutic avenue, the successful development of epigenetic modulators has been challenging. We argue that besides identifying the right responder patient population, the selection of an optimized dosing regimen is equally important. For the majority of epigenetic modulators, hematological adverse effects such as thrombocytopenia, anemia or neutropenia are frequently observed and may limit their therapeutic potential. Therefore, one of the key challenges is to identify a dosing regimen that maximizes drug efficacy and minimizes toxicity. This requires a good understanding of the quantitative relationship between the administered dose, the drug exposure and the magnitude and duration of drug response related to safety and efficacy. With case examples, we highlight how modeling and simulation has been successfully applied to address those questions. As an outlook, we suggest the combination of efficacy and safety prediction models that capture the quantitative, mechanistic relationships governing the balance between their safety and efficacy dynamics. A stepwise approach for its implementation is presented. Utilizing in silico explorations, the impact of dosing regimen on the therapeutic window can be explored. This will serve as a basis to select the most promising dosing regimen that maximizes efficacy while minimizing adverse effects and to increase the probability of success for the given epigenetic drug.

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.

Unraveling Disease Pathophysiology with Mathematical Modeling

Modeling has enabled fundamental advances in our understanding of the mechanisms of health and disease for centuries, since at least the time of William Harvey almost 500 years ago. Recent technological advances in molecular methods, computation, and imaging generate optimism that mathematical modeling will enable the biomedical research community to accelerate its efforts in unraveling the molecular, cellular, tissue-, and organ-level processes that maintain health, predispose to disease, and determine response to treatment. In this review, we discuss some of the roles of mathematical modeling in the study of human physiology and pathophysiology and some challenges and opportunities in general and in two specific areas: in vivo modeling of pulmonary function and in vitro modeling of blood cell populations.

Model-based approach for methoxy polyethylene glycol-epoetin beta drug development in paediatric patients with anaemia of chronic kidney disease

AIMS

Methoxy polyethylene glycol-epoetin beta (continuous erythropoietin receptor activator, C.E.R.A.) is used for the treatment of anaemia in adults with chronic kidney disease (CKD). Patients treated with shorter-acting erythropoiesis-stimulating agents up to three times weekly can be switched to once-monthly C.E.R.A.. Doses can be adjusted on a monthly basis based on haemoglobin (Hb) levels during the preceding period. A model-based approach was applied to optimise C.E.R.A. development, more specifically the confirmatory trial of the paediatric plan.

METHODS

Pharmacokinetic and pharmacodynamic data from a phase II paediatric study and phase II and III adult studies were analysed together using modelling and simulation to determine the pharmacokinetic/pharmacodynamic characteristics of C.E.R.A. in a broad population. Model-based simulations of C.E.R.A. treatment outcomes in paediatric patients were performed, notably when administered subcutaneously and compared to clinical and real-world data.

RESULTS

Age and body weight explained differences in pharmacokinetics, while the pharmacodynamic characteristics of C.E.R.A. were similar between adult and paediatric populations. Simulated Hb levels (mean and 95% prediction interval 10.9 [10.6, 11.2] g dL^-1 ) and C.E.R.A. doses (median and 95% prediction interval 105 [72, 159] μg) 20 weeks after switching to subcutaneous C.E.R.A. were confirmed by observed real-world data from International Pediatric Dialysis Network registries (mean Hb was 10.8 g dL^-1 and median C.E.R.A. dose was 100 μg).

CONCLUSIONS

These analyses have facilitated optimisation of the C.E.R.A. development programme in paediatric patients with anaemia of CKD to provide this patient population with faster access to the drug while avoiding unnecessary clinical trial exposure and related monitoring burden in children.

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.

Toxicodynetics: A new discipline in clinical toxicology

OBJECTIVES

Regarding the different disciplines that encompass the pharmacology and the toxicology, none is specifically dedicated to the description and analysis of the time-course of relevant toxic effects both in experimental and clinical studies. The lack of a discipline devoted to this major field in toxicology results in misconception and even in errors by clinicians.

MATERIAL AND METHODS

Review of the basic different disciplines that encompass pharmacology toxicology and comparing with the description of the time-course of effects in conditions in which toxicological analysis was not performed or with limited analytical evidence.

RESULTS

Review of the literature clearly shows how misleading is the current extrapolation of toxicokinetic data to the description of the time-course of toxic effects.

CONCLUSION

A new discipline entitled toxicodynetics should be developed aiming at a more systematic description of the time-course of effects in acute human and experimental poisonings. Toxicodynetics might help emergency physicians in risk assessment when facing a poisoning and contribute to a better assessment of quality control of data collected by poison control centres. Toxicodynetics would also allow a quantitative approach to the clinical effects resulting from drug-drug interaction.

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.

Lifespan based pharmacokinetic-pharmacodynamic model of tumor growth inhibition by anticancer therapeutics

Accurate prediction of tumor growth is critical in modeling the effects of anti-tumor agents. Popular models of tumor growth inhibition (TGI) generally offer empirical description of tumor growth. We propose a lifespan-based tumor growth inhibition (LS TGI) model that describes tumor growth in a xenograft mouse model, on the basis of cellular lifespan T. At the end of the lifespan, cells divide, and to account for tumor burden on growth, we introduce a cell division efficiency function that is negatively affected by tumor size. The LS TGI model capability to describe dynamic growth characteristics is similar to many empirical TGI models. Our model describes anti-cancer drug effect as a dose-dependent shift of proliferating tumor cells into a non-proliferating population that die after an altered lifespan T_A. Sensitivity analysis indicated that all model parameters are identifiable. The model was validated through case studies of xenograft mouse tumor growth. Data from paclitaxel mediated tumor inhibition was well described by the LS TGI model, and model parameters were estimated with high precision. A study involving a protein casein kinase 2 inhibitor, AZ968, contained tumor growth data that only exhibited linear growth kinetics. The LS TGI model accurately described the linear growth data and estimated the potency of AZ968 that was very similar to the estimate from an established TGI model. In the case study of AZD1208, a pan-Pim inhibitor, the doubling time was not estimable from the control data. By fixing the parameter to the reported in vitro value of the tumor cell doubling time, the model was still able to fit the data well and estimated the remaining parameters with high precision. We have developed a mechanistic model that describes tumor growth based on cell division and has the flexibility to describe tumor data with diverse growth kinetics.

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.

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

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.

Erythropoietic response to endogenous erythropoietin in premature very low birth weight infants

Despite the common occurrence of anemia in very low birth weight (VLBW) infants, the erythropoiesis and Hb production rates and their relationship to plasma erythropoietin (EPO) concentrations remain unknown in these subjects. To determine these quantities, all blood removed by phlebotomy and administered by red blood cell (RBC) transfusion over the first 30 days of life was recorded in 14 ventilated VLBW infants born at 24 to 28 weeks of gestation. Discarded blood from frequent clinically ordered laboratory blood samples was used to construct plasma EPO, Hb, and RBC concentration-time profiles for each infant. A pharmacodynamic Hb mass balance model that accounted for the dynamic hematological conditions experienced by these infants was simultaneously fitted to the plasma EPO, Hb, and RBC concentrations from each individual subject, while accounting for subject growth. Based on the model estimates, an average of 4.69 g of Hb was produced over the first 30 days of life, compared with 5.97 g removed by phlebotomies and 12.3 g administered by transfusions. These high transfusion amounts were consistent with a relatively short RBC life span and rapidly expanding blood volume with infant growth. The estimated mean body weight-scaled Hb production rate dropped nearly 3-fold after birth to 0.144 g/day x (kg)(3/4). Although only estimated in a subset of the subjects, the mean plasma EPO EC(50) of 28.5 mU/ml and maximal Hb production rate (E(max)) indicated that a severalfold increase in Hb production rate could be achieved with only a modest increase in plasma EPO concentrations.

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.

Basic pharmacodynamic models for agents that alter the lifespan distribution of natural cells

A new class of basic indirect pharmacodynamic models for agents that alter the loss of natural cells based on a lifespan concept are presented. The lifespan indirect response (LIDR) models assume that cells (R) are produced at a constant rate (k(in)), survive during a certain duration T(R), and finally are lost. The rate of cell loss is equal to the production rate but is delayed by T(R). A therapeutic agent can increase or decrease the baseline cell lifespan to a new cell lifespan, T(D), by temporally changing the proportion of cells belonging to the two modes of the lifespan distribution. Therefore, the change of lifespan at time t is described according to the Hill function, H(C(t)), with capacity (E(max)) and sensitivity (EC(50)), and the pharmacokinetic function C(t). A one-compartment cell model was examined through simulations to describe the role of pharmacokinetics, pharmacodynamics and cell properties for the cases where the drug increases (T(D) > T(R)) or decreases (T(D) < T(R)) the cell lifespan. The area under the effect curve (AUCE) and explicit solutions of LIDR models for large doses were derived. The applicability of the model was further illustrated using the effects of recombinant human erythropoietin (rHuEPO) on reticulocytes. The cases of both stimulation of the proliferation of bone marrow progenitor cells and the increase of reticulocyte lifespans were used to describe mean data from healthy subjects who received single subcutaneous doses of rHuEPO ranging from 20 to 160 kIU. rHuEPO is about 4.5-fold less potent in increasing reticulocyte survival than in stimulating the precursor production. A maximum increase of 4.1 days in the mean reticulocyte lifespan was estimated and the effect duration on the lifespan distribution was dose dependent. LIDR models share similar properties with basic indirect response models describing drug stimulation or inhibition of the response loss rate with the exception of the presence of a lag time and a dose independent peak time. The current concept can be applied to describe the pharmacodynamic effects of agents affecting survival of hematopoietic cell populations yielding realistic physiological parameters.

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.

Pharmacodynamic model of interleukin-21 effects on red blood cells in cynomolgus monkeys

Interleukin-21 (IL-21) is a novel cytokine that is currently under clinical investigations as a potential anti-cancer agent. Like many other anti-cancer agents, including other interleukins, IL-21 is seen to produce a broad range of biological effects that may be related to both efficacy and safety of treatment. The present analysis investigates the observed pharmacodynamics effects on red blood cells following various treatment schedules of human IL-21 administrated to cynomolgus monkeys. These effects are described by a novel non-linear mixed-effects model that enabled separation of drug effects and sampling effects, the latter believed to be due partly to blood loss and partly to stress induced haemolysis in connection with blood sampling. Two different studies with a total of 9 different treatment groups of cynomolgus monkeys were used for model development. In conclusion, the model describes the IL-21 induced drop in red blood cells to be (1) caused by removal rather than suppression of production, consistent with increased reticulocyte concentration, and (2) considerably delayed compared to dosing, i.e. not related to the drop in red blood cells observed immediately post dose. It is believed that the structural model presented here can be used for other types of drug induced loss of red blood cells, whereas the mechanism for sampling related blood loss is relevant for investigations of anaemia in all pharmacological studies with smaller animals.

Population cell life span models for effects of drugs following indirect mechanisms of action

Pharmacokinetic/pharmacodynamic (PK/PD) models for hematological drug effects exist that assume that cells are produced by a zero- or first-order process, survive for a specific duration (cell lifespan), and then are lost. Due to the fact that delay differential equations (DDE) are needed for cell lifespan models, their software implementation is not straightforward. Our objective is to demonstrate methods to implement three different cell lifespan models for dealing with hematological drug effects and to evaluate the performance of NONMEM to estimate the model parameters. For the basic lifespan indirect response (LIDR) model, cells are produced by a zero-order process and removed due to senescence. The modified LIDR model adds a precursor pool. The LIDR model of cytotoxicity assumes a three-pool indirect model to account for the cell proliferation with capacity-limited cytotoxicity followed by maturation, and removal from the circulation. A numerical method (method of steps) implementing DDE in NONMEM was introduced. Simulation followed by estimation was used to evaluate NONMEM performance and the impact of the minimization algorithm (first-order method vs. first-order conditional estimation method) and the model for residual variability on the estimates of the population parameters. The FOCE method combined with log-transformation of data was found to be superior. This report provides methodology that will assist in application of population methods for assessing hematological responses to various types of drugs.

Interactions of prednisolone and other immunosuppressants used in dual treatment of systemic lupus erythematosus in lymphocyte proliferation assays

Systemic lupus erythematosus is an autoimmune disease primarily affecting women. Currently, systemic lupus erythematosus therapy is suboptimal due to adverse effects of immunosuppressants, particularly corticosteroids. This study determines the single effects of prednisolone, dehydroepiandrosterone, bromocriptine, tamoxifen, mycophenolic acid, 2-chloro-2'-deoxyadenosine, azathioprine, and chloroquine on lectin-stimulated proliferation of human T lymphocytes, as well as determining whether there are interactions in the joint effects of prednisolone and these agents. The T lymphocytes from the whole blood of 10 middle-aged women were stimulated by phytohemagglutinin and cultured with varying drug concentrations. The Hill function was used to evaluate single-drug response data. Isobolograms were constructed to qualitatively analyze interactions. Parametric analysis based on competitive and noncompetitive interaction models was further applied to quantify the joint interactions and predict steroid-sparing potential. The surface interaction parameter (psi) estimated from parametric analysis was in concordance with isobolographic inspection for all interactions studied. All interactions favored the noncompetitive model. Results suggest that dehydroepiandrosterone is additive in its effect with prednisolone, whereas tamoxifen interacts synergistically, both providing steroid-sparing effects. Novel immuno-suppressants such as mycophenolic acid may still provide added pharmacologic benefit during therapy despite a slight antagonistic interaction with prednisolone. These studies help rationalize actual or potential use of other drugs with prednisolone in the treatment of systemic lupus erythematosus.