In silico trial for the assessment of givinostat dose adjustment rules based on the management of key hematological parameters in polycythemia vera patients.

Polycythemia vera (PV) is a chronic myeloproliferative neoplasm characterized by excessive levels of platelets (PLT), white blood cells (WBC), and hematocrit (HCT). Givinostat (ITF2357) is a potent histone-deacetylase inhibitor that showed a good safety/efficacy profile in PV patients during phase I/II studies. A phase III clinical trial had been planned and an adaptive dosing protocol had been proposed where givinostat dose is iteratively adjusted every 28 days (one cycle) based on PLT, WBC, and HCT. As support, a simulation platform to evaluate and refine the proposed givinostat dose adjustment rules was developed. A population pharmacokinetic/pharmacodynamic model predicting the givinostat effects on PLT, WBC, and HCT in PV patients was developed and integrated with a control algorithm implementing the adaptive dosing protocol. Ten in silico trials in ten virtual PV patient populations were simulated 500 times. Considering an eight-treatment cycle horizon, reducing/increasing the givinostat daily dose by 25 mg/day step resulted in a higher percentage of patients with a complete hematological response (CHR), that is, PLT ≤400 × 109 /L, WBC ≤10 × 109 /L, and HCT < 45% without phlebotomies in the last three cycles, and a lower percentage of patients with grade II toxicity events compared with 50 mg/day adjustment steps. After the eighth cycle, 85% of patients were predicted to receive a dose ≥100 mg/day and 40.90% (95% prediction interval = [34, 48.05]) to show a CHR. These results were confirmed at the end of 12th, 18th, and 24th cycles, showing a stability of the response between the eighth and 24th cycles.

Population pharmacokinetic-pharmacodynamic analysis of givinostat.

BACKGROUND: Givinostat (ITF2357), an oral, synthetic histone deacetylase inhibitor, significantly improved all histological muscle biopsy parameters in a Phase II study in boys with Duchenne muscular dystrophy (DMD). RESEARCH DESIGN AND METHODS: A population pharmacokinetic (PK) model, including seven clinical studies, was developed to explore the effect of covariates on givinostat PK. The final model was qualified to simulate pediatric dosing recommendations. A PK/pharmacodynamic (PD) model was developed to simulate the link between givinostat plasma concentration and platelet time course in 10-70-kg children following 6 months of givinostat 20-70 mg twice daily. RESULTS: A two-compartment model, with first-order input with lag and first-order elimination from the central compartment, described givinostat PK, demonstrating increasing apparent clearance with increasing body weight. The PK/PD model well-described platelet count time course. Weight-based dosing (arithmetic mean systemic exposure of 554-641 ng·h/mL) produced an average platelet count decrease from baseline of 45% with maximum decrease within 28 days. After 1 week and 6 months, ~1% and ~14-15% of patients, respectively, had a platelet count <75 × 109/L. CONCLUSIONS: Based on these data, givinostat dosing will be body weight adjusted and include monitoring of platelet counts to support efficacy and safety in a Phase III DMD study.

A Dynamic Energy Budget (DEB) based modeling framework to describe tumor-in-host growth inhibition and cachexia onset during anticancer treatment in in vivo xenograft studies.

Cancer anorexia-cachexia syndrome (CACS) is a very severe complication of cancer for which an adequate therapeutic strategy has not yet been defined. Recently, a notable number of new animal models of human CACS has been made available for translational purposes. Under the assumption that tumor-induced adaptations of host metabolism and tumor-host energetic competition play a major role in CACS (together with possible toxicities induced by the anticancer treatment), we developed a new Dynamic Energy Budget (DEB)-based framework, modeling tumor-in-host growth dynamics and cachexia onset in preclinical animal models during anticancer treatments. The tumor-in-host modeling approach was successfully applied on a multitude of in vivo preclinical studies involving different host species, tumor cell lines, type of anticancer agents and experimental settings among which standard xenograft studies. Obtained results strongly suggested the adoption of the tumor-in-host DEB-based approach in the preclinical oncological setting for a joint assessment of drug efficacy and toxicity and for a better design of the experiments. Further applications of the DEB-based approach to the context of poly-targeted combination therapy, anti-cachectic treatments and preclinical to clinical translation are under investigation with extremely encouraging preliminary results.

A Tumor-in-Host DEB-Based Approach for Modeling Cachexia and Bevacizumab Resistance.

Adequate energy intake and homeostasis are fundamental for the appropriate growth and maintenance of an organism; the presence of a tumor can break this equilibrium. Tumor energy requests can lead to extreme weight loss in animals and cachexia in cancer patients. Angiogenesis inhibitors, acting on tumor vascularization, counteract this tumor-host energy imbalance, with significant results in preclinical models and more limited results in the clinic. Current pharmacokinetic-pharmacodynamic models mainly focus on the antiangiogenic effects on tumor growth but do not provide information about host conditions. A model that can predict energetic conditions that provide significant tumor growth inhibition with acceptable host body weight reduction is therefore needed. We developed a new tumor-in-host dynamic energy budget (DEB)-based model to account for the cytostatic activity of antiangiogenic treatments. Drug effect was implemented as an inhibition of the energy fraction subtracted from the host by the tumor. The model was tested on seven xenograft experiments involving bevacizumab and three different tumor cell lines. The model successfully predicted tumor and host body growth data, providing a quantitative measurement of drug potency and tumor-related cachexia. The inclusion of a hypoxia-triggered resistance mechanism enabled investigation of the decreased efficacy frequently observed with prolonged bevacizumab treatments. In conclusion, the tumor-in-host DEB-based approach has been extended to account for the effect of bevacizumab. The resistance model predicts the response to different administration protocols and, for the first time, the impact of tumor-related cachexia in different cell lines. Finally, the physiologic base of the model strongly suggests its use in translational human research. SIGNIFICANCE: A mathematical model describes tumor growth in animal models, taking into consideration the energy balance involving both the growth of tumor and the physiologic functions of the host.

Histological effects of givinostat in boys with Duchenne muscular dystrophy.

Duchenne Muscular Dystrophy (DMD) is caused by mutations in the dystrophin gene leading to dystrophin deficiency, muscle fiber degeneration and progressive fibrotic replacement of muscles. Givinostat, a histone deacetylase (HDAC) inhibitor, significantly reduced fibrosis and promoted compensatory muscle regeneration in mdx mice. This study was conducted to evaluate whether the beneficial histological effects of Givinostat could be extended to DMD boys. Twenty ambulant DMD boys aged 7 to <11 years on stable corticosteroid treatment were enrolled in the study and treated for ≥12 months with Givinostat. A muscle biopsy was collected at the beginning and at the end of treatment to evaluate the amount of muscle and fibrotic tissue. Histological effects were the primary objectives of the study. Treatment with Givinostat significantly increased the fraction of muscle tissue in the biopsies and reduced the amount of fibrotic tissue. It also substantially reduced tissue necrosis and fatty replacement. Overall the drug was safe and tolerated. Improvement in functional tests was not observed in this study, but the sample size of the study was not sufficient to draw definitive conclusions. This study showed that treatment with Givinostat for more than 1 year significantly counteracted histological disease progression in ambulant DMD boys aged 7 to 10 years.

Modeling of human tumor xenografts and dose rationale in oncology.

Xenograft models are commonly used in oncology drug development. Although there are discussions about their ability to generate meaningful data for the translation from animal to humans, it appears that better data quality and better design of the preclinical experiments, together with appropriate data analysis approaches could make these data more informative for clinical development. An approach based on mathematical modeling is necessary to derive experiment-independent parameters which can be linked with clinically relevant endpoints. Moreover, the inclusion of biomarkers as predictors of efficacy is a key step towards a more general mechanism-based strategy.

Predictive pharmacokinetic-pharmacodynamic modeling of tumor growth after administration of an anti-angiogenic agent, bevacizumab, as single-agent and combination therapy in tumor xenografts.

PURPOSE: Pharmacokinetic-pharmacodynamic (PK-PD) models able to predict the action of anticancer compounds in tumor xenografts have an important impact on drug development. In case of anti-angiogenic compounds, many of the available models show difficulties in their applications, as they are based on a cell kill hypothesis, while these drugs act on the tumor vascularization, without a direct tumor cell kill effect. For this reason, a PK-PD model able to describe the tumor growth modulation following treatment with a cytostatic therapy, as opposed to a cytotoxic treatment, is proposed here. METHODS: Untreated tumor growth was described using an exponential growth phase followed by a linear one. The effect of anti-angiogenic compounds was implemented using an inhibitory effect on the growth function. The model was tested on a number of experiments in tumor-bearing mice given the anti-angiogenic drug bevacizumab either alone or in combination with another investigational compound. Nonlinear regression techniques were used for estimating the model parameters. RESULTS: The model successfully captured the tumor growth data following different bevacizumab dosing regimens, allowing to estimate experiment-independent parameters. A combination model was also developed under a 'no-interaction' assumption to assess the effect of the combination of bevacizumab with a target-oriented agent. The observation of a significant difference between model-predicted and observed tumor growth curves was suggestive of the presence of a pharmacological interaction that was further accommodated into the model. CONCLUSIONS: This approach can be used for optimizing the design of preclinical experiments. With all the inherent limitations, the estimated experiment-independent model parameters can be used to provide useful indications for the single-agent and combination regimens to be explored in the subsequent development phases.

A4S: a user-friendly graphical tool for pharmacokinetic and pharmacodynamic (PK/PD) simulation.

Effective communication of PK/PD principles and results in a biomedical research environment remains a significant challenge which can result in lack of buy-in and engagement from scientists outside the modeling and simulation communities. In our view, one of the barriers in this area is a lack of user-friendly tools which allow "non experts" to use PK/PD models without the need to develop technical skills and expertise in advanced mathematical principles and specialist software. The costs of commercial software may also prevent large-scale distribution. One attempt to address this issue internally in our research organizations has resulted in the development of the A4S ("Accelera for Sandwich") software, which is a simple-to-use, menu-drive Matlab-based PK/PD simulator targeted at biomedical researchers with little PK/PD experience.

NMS-P937, an orally available, specific small-molecule polo-like kinase 1 inhibitor with antitumor activity in solid and hematologic malignancies.

Polo-like kinase 1 (PLK1) is a serine/threonine protein kinase considered to be the master player of cell-cycle regulation during mitosis. It is indeed involved in centrosome maturation, bipolar spindle formation, chromosome separation, and cytokinesis. PLK1 is overexpressed in a variety of human tumors and its overexpression often correlates with poor prognosis. Although five different PLKs are described in humans, depletion or inhibition of kinase activity of PLK1 is sufficient to induce cell-cycle arrest and apoptosis in cancer cell lines and in xenograft tumor models. NMS-P937 is a novel, orally available PLK1-specific inhibitor. The compound shows high potency in proliferation assays having low nanomolar activity on a large number of cell lines, both from solid and hematologic tumors. NMS-P937 potently causes a mitotic cell-cycle arrest followed by apoptosis in cancer cell lines and inhibits xenograft tumor growth with clear PLK1-related mechanism of action at well-tolerated doses in mice after oral administration. In addition, NMS-P937 shows potential for combination in clinical settings with approved cytotoxic drugs, causing tumor regression in HT29 human colon adenocarcinoma xenografts upon combination with irinotecan and prolonged survival of animals in a disseminated model of acute myelogenous leukemia in combination with cytarabine. NMS-P937, with its favorable pharmacologic parameters, good oral bioavailability in rodent and nonrodent species, and proven antitumor activity in different preclinical models using a variety of dosing regimens, potentially provides a high degree of flexibility in dosing schedules and warrants investigation in clinical settings.

Identification of 4,5-dihydro-1H-pyrazolo[4,3-h]quinazoline derivatives as a new class of orally and selective Polo-like kinase 1 inhibitors.

Polo-like kinase 1 (Plk1) is a fundamental regulator of mitotic progression whose overexpression is often associated with oncogenesis and therefore is recognized as an attractive therapeutic target in the treatment of proliferative diseases. Here we discuss the structure-activity relationship of the 4,5-dihydro-1H-pyrazolo[4,3-h]quinazoline class of compounds that emerged from a high throughput screening (HTS) campaign as potent inhibitors of Plk1 kinase. Furthermore, we describe the discovery of 49, 8-{[2-methoxy-5-(4-methylpiperazin-1-yl)phenyl]amino}-1-methyl-4,5-dihydro-1H-pyrazolo[4,3-h]quinazoline-3-carboxamide, as a highly potent and specific ATP mimetic inhibitor of Plk1 (IC(50) = 0.007 microM) as well as its crystal structure in complex with the methylated Plk1(36-345) construct. Compound 49 was active in cell proliferation against different tumor cell lines with IC(50) values in the submicromolar range and active in vivo in the HCT116 xenograft model where it showed 82% tumor growth inhibition after repeated oral administration.

Modelling the genesis and treatment of cancer: the potential role of physiologically based pharmacodynamics.

Physiologically based modelling of pharmacodynamics/toxicodynamics requires an a priori knowledge on the underlying mechanisms causing toxicity or causing the disease. In the context of cancer, the objective of the expert meeting was to discuss the molecular understanding of the disease, modelling approaches used so far to describe the process, preclinical models of cancer treatment and to evaluate modelling approaches developed based on improved knowledge. Molecular events in cancerogenesis can be detected using 'omics' technology, a tool applied in experimental carcinogenesis, but also for diagnostics and prognosis. The molecular understanding forms the basis for new drugs, for example targeting protein kinases specifically expressed in cancer. At present, empirical preclinical models of tumour growth are in great use as the development of physiological models is cost and resource intensive. Although a major challenge in PKPD modelling in oncology patients is the complexity of the system, based in part on preclinical models, successful models have been constructed describing the mechanism of action and providing a tool to establish levels of biomarker associated with efficacy and assisting in defining biologically effective dose range selection for first dose in man. To follow the concentration in the tumour compartment enables to link kinetics and dynamics. In order to obtain a reliable model of tumour growth dynamics and drug effects, specific aspects of the modelling of the concentration-effect relationship in cancer treatment that need to be accounted for include: the physiological/circadian rhythms of the cell cycle; the treatment with combinations and the need to optimally choose appropriate combinations of the multiple agents to study; and the schedule dependence of the response in the clinical situation.

A phase I dose-escalation study of danusertib (PHA-739358) administered as a 24-hour infusion with and without granulocyte colony-stimulating factor in a 14-day cycle in patients with advanced solid tumors.

PURPOSE: This study was conducted to assess the safety, tolerability, pharmacokinetics, and pharmacodynamics of the i.v. pan-aurora kinase inhibitor PHA-739358, danusertib, in patients with advanced solid tumors. EXPERIMENTAL DESIGN: In part 1, patients received escalating doses of danusertib (24-hour infusion every 14 days) without filgrastim (granulocyte colony-stimulating factor, G-CSF). Febrile neutropenia was the dose-limiting toxicity without G-CSF. Further dose escalation was done in part 2 with G-CSF. Blood samples were collected for danusertib pharmacokinetics and pharmacodynamics. Skin biopsies were collected to assess histone H3 phosphorylation (pH3). RESULTS: Fifty-six patients were treated, 40 in part 1 and 16 in part 2. Febrile neutropenia was the dose-limiting toxicity in part 1 without G-CSF. Most other adverse events were grade 1 to 2, occurring at doses >or=360 mg/m(2) with similar incidence in parts 1 and 2. The maximum tolerated dose without G-CSF is 500 mg/m(2). The recommended phase 2 dose in part 2 with G-CSF is 750 mg/m(2). Danusertib showed dose-proportional pharmacokinetics in parts 1 and 2 with a median half-life of 18 to 26 hours. pH3 modulation in skin biopsies was observed at >or=500 mg/m(2). One patient with refractory small cell lung cancer (1,000 mg/m(2) with G-CSF) had an objective response lasting 23 weeks. One patient with refractory ovarian cancer had 27% tumor regression and 30% CA125 decline. CONCLUSIONS: Danusertib was well tolerated with target inhibition in skin at >or=500 mg/m(2). Preliminary evidence of antitumor activity, including a partial response and several occurrences of prolonged stable disease, was seen across a variety of advanced refractory cancers. Phase II studies are ongoing.

A model-based approach to the in vitro evaluation of anticancer activity.

PURPOSE: The use of in vitro screening tests for characterizing the activity of anticancer agents is a standard practice in oncology research and development. In these studies, human A2780 ovarian carcinoma cells cultured in plates are exposed to different concentrations of the compounds for different periods of time. Their anticancer activity is then quantified in terms of EC(50) comparing the number of metabolically active cells present in the treated and the control arms at specified time points. The major concern of this methodology is the observed dependency of the EC(50) on the experimental design in terms of duration of exposure. This dependency could affect the efficacy ranking of the compounds, causing possible biases especially in the screening phase, when compound selection is the primary purpose of the in vitro analysis. To overcome this problem, the applicability of a modeling approach to these in vitro studies was evaluated. METHODS: The model, consisting of a system of ordinary differential equations, represents the growth of tumor cells using a few identifiable and biologically relevant parameters related to cell proliferation dynamics and drug action. In particular, the potency of the compounds can be measured by a unique and drug-specific parameter that is essentially independent of drug concentration and exposure time. Parameter values were estimated using weighted nonlinear least squares. RESULTS: The model was able to adequately describe the growth of tumor cells at different experimental conditions. The approach was validated both on commercial drugs and discovery candidate compounds. In addition, from this model the relationship between EC(50) and the exposure time was derived in an analytic form. CONCLUSIONS: The proposed approach provides a new tool for predicting and/or simulating cell responses to different treatments with useful indications for optimizing in vitro experimental designs. The estimated potency parameter values obtained from different compounds can be used for an immediate ranking of anticancer activity.

A minimal model of tumor growth inhibition.

The preclinical development of antitumor drugs greatly benefits from the availability of models capable of predicting tumor growth as a function of the drug administration schedule. For being of practical use, such models should be simple enough to be identifiable from standard experiments conducted on animals. In the present paper, a stochastic model is derived from a set of minimal assumptions formulated at cellular level. Tumor cells are divided in two groups: proliferating and nonproliferating. The probability that a proliferating cell generates a new cell is a function of the tumor weight. The probability that a proliferating cell becomes nonproliferating is a function of the plasma drug concentration. The time-to-death of a nonproliferating cell is a random variable whose distribution reflects the nondeterministic delay between drug action and cell death. The evolution of the expected value of tumor weight obeys two differential equations (an ordinary and a partial differential one), whereas the variance is negligible. Therefore, the tumor growth dynamics can be well approximated by the deterministic evolution of its expected value. The tumor growth inhibition model, which is a lumped parameter model that in the last few years has been successfully applied to several antitumor drugs, is shown to be a special case of the minimal model presented here.

PHA-739358, a potent inhibitor of Aurora kinases with a selective target inhibition profile relevant to cancer.

PHA-739358 is a small-molecule 3-aminopyrazole derivative with strong activity against Aurora kinases and cross-reactivities with some receptor tyrosine kinases relevant for cancer. PHA-739358 inhibits all Aurora kinase family members and shows a dominant Aurora B kinase inhibition-related cellular phenotype and mechanism of action in cells in vitro and in vivo. p53 status-dependent endoreduplication is observed upon treatment of cells with PHA-739358, and phosphorylation of histone H3 in Ser(10) is inhibited. The compound has significant antitumor activity in different xenografts and spontaneous and transgenic animal tumor models and shows a favorable pharmacokinetic and safety profile. In vivo target modulation is observed as assessed by the inhibition of the phosphorylation of histone H3, which has been validated preclinically as a candidate biomarker for the clinical phase. Pharmacokinetics/pharmacodynamics modeling was used to define drug potency and to support the prediction of active clinical doses and schedules. We conclude that PHA-739358, which is currently tested in clinical trials, has great therapeutic potential in anticancer therapy in a wide range of cancers.

1,4,5,6-tetrahydropyrrolo[3,4-c]pyrazoles: identification of a potent Aurora kinase inhibitor with a favorable antitumor kinase inhibition profile.

The optimization of a series of 5-phenylacetyl 1,4,5,6-tetrahydropyrrolo[3,4-c]pyrazole derivatives toward the inhibition of Aurora kinases led to the identification of compound 9d. This is a potent inhibitor of Aurora kinases that also shows low nanomolar potency against additional anticancer kinase targets. Based on its high antiproliferative activity on different cancer cell lines, favorable chemico-physical and pharmacokinetic properties, and high efficacy in in vivo tumor models, compound 9d was ultimately selected for further development.

A predictive model for exemestane pharmacokinetics/pharmacodynamics incorporating the effect of food and formulation.

AIMS: Exemestane (Aromasin) is an irreversible aromatase inactivator used for the treatment of postmenopausal women with advanced breast cancer. The objective of this study was to evaluate the effect of formulation comparing a sugar-coated tablet (SCT) with a suspension and food on the pharmacokinetics (PK) and pharmacodynamics (PD) with respect to plasma estrone sulphate (E1S) concentrations of exemestane, using a PK/PD approach. METHODS: This was an open, three-period, randomized, crossover study. Twelve healthy postmenopausal women received single oral doses of 25 mg exemestane as a SCT after fasting or food and as a suspension after fasting. Exemestane and E1S concentrations were determined before and up to 14 days after drug administration. Population analysis was performed in two steps: (i) a compartmental PK model was selected incorporating the effect of food and formulation; (ii) conditional on the PK model, a PD model was developed employing indirect response models. Model selection was performed using standard statistical tests. Validation and assessment of the predictive capability of the selected model was performed using real test data sets obtained from the literature. RESULTS: A three-compartment model with first-order elimination rate best described exemestane disposition (k12 0.454, k21 0.158, k13 0.174, k31 0.016 and k 0.738 h(-1)). Absorption was described by a mono-exponential function [ka 2.3 (SCT after fasting), 1.1 (SCT after food) and 7.6 h(-1) (suspension); lag time 0.2 h]. The PD model assumed that E1S plasma concentrations are determined by a zero-order synthesis rate (6.5 pg ml(-1) h(-1)) and a first-order elimination constant (0.032 h(-1)). Exemestane inhibited E1S synthesis with a C50 value of 22.1 pg ml(-1). The mean population estimates were used to simulate the administration of different doses of the drug (0.5, 1, 2.5, 5 and 25 mg day(-1)). The model predictions were in agreement with historical data. CONCLUSIONS: Exemestane absorption is influenced by the formulation of the drug and by food, but its disposition is independent of both. PK differences do no translate into clinically important differences in the PD. The PK/PD model developed was able to predict successfully the response to different doses and administration schedules with respect to oestrogen suppression.

The effects of degree of hepatic or renal impairment on the pharmacokinetics of exemestane in postmenopausal women.

PURPOSE: Two studies were conducted to compare the pharmacokinetics and tolerability of exemestane in postmenopausal subjects with various degrees of impairment of hepatic or renal function with those in healthy postmenopausal subjects. METHODS: All subjects were postmenopausal females. In study 1, nine subjects had normal hepatic function (Child-Pugh grade A), and nine had moderately (grade B) and eight severely (grade C) impaired hepatic function. In study 2, six subjects had normal renal function, and six moderately (creatinine clearance, CrCL, 30-60 ml/min per 1.73 m(2)) and seven severely (CrCL <29 ml/min per 1.73 m(2)) impaired renal function. Each subject took a single oral dose of 25 mg exemestane. Samples of plasma (to 168 h after dosing) and urine (to 72 h in study 1, or 96 h in study 2) were taken for pharmacokinetic analysis. Safety and tolerability were assessed by monitoring vital signs, laboratory safety tests, ECG and adverse events. RESULTS: Exposure to exemestane was increased two- to threefold in patients with hepatic impairment. Thus, the geometric mean AUC(0- infinity ) values were 41.71 (90% CI 32.2 to 54.0), 99.02 (76.5 to 128.2) and 118.59 ng.h/ml (90.2 to 156.0) in healthy subjects, and in patients with moderate and severe hepatic impairment, respectively ( P<0.01). C(max) also increased twofold. Compared with healthy subjects, patients with hepatic impairment had lower apparent oral clearance and apparent volume of distribution of exemestane. Renal impairment was also associated with two- to threefold increases in AUC(0- infinity ): 34.64 (90% CI 23.9 to 50.2), 94.10 (64.9 to 136.4) and 65.52 ng.h/ml (46.5 to 92.4) in healthy subjects, and in patients with moderate and severe hepatic impairment, respectively ( P<0.05). C(max) did not change significantly. Apparent oral clearance was directly correlated with CrCL ( r(2)=0.43). Exemestane was tolerated well, with no safety concerns. CONCLUSIONS: Oral clearance of exemestane was reduced in the presence of significant hepatic or renal disease. However, in view of the relatively large safety margin and the mild nature of the side effects of exemestane, the therapeutic implications of these changes in pharmacokinetics are considered minor and of no clinical significance.