PKPD Modeling of Predictors for Adverse Effects and Overall Survival in Sunitinib-Treated Patients With GIST

A Pharmacodynamic Study of Aminoglycosides against Pathogenic E. coli through Monte Carlo Simulation

This research focuses on combating the increasing problem of antimicrobial resistance, especially in Escherichia coli (E. coli), by assessing the efficacy of aminoglycosides. The study specifically addresses the challenge of developing new therapeutic approaches by integrating experimental data with mathematical modeling to better understand the action of aminoglycosides. It involves testing various antibiotics like streptomycin (SMN), kanamycin (KMN), gentamicin (GMN), tobramycin (TMN), and amikacin (AKN) against the O157:H7 strain of E. coli. The study employs a pharmacodynamics (PD) model to analyze how different antibiotic concentrations affect bacterial growth, utilizing minimum inhibitory concentration (MIC) to gauge the effective bactericidal levels of the antibiotics. The study's approach involved transforming bacterial growth rates, as obtained from time-kill curve data, into logarithmic values. A model was then developed to correlate these log-transformed values with their respective responses. To generate additional data points, each value was systematically increased by an increment of 0.1. To simulate real-world variability and randomness in the data, a Gaussian scatter model, characterized by parameters like κ and EC50, was employed. The mathematical modeling was pivotal in uncovering the bactericidal properties of these antibiotics, indicating different PD MIC (zMIC) values for each (SMN: 1.22; KMN: 0.89; GMN: 0.21; TMN: 0.32; AKN: 0.13), which aligned with MIC values obtained through microdilution methods. This innovative blend of experimental and mathematical approaches in the study marks a significant advancement in formulating strategies to combat the growing threat of antimicrobial-resistant E. coli, offering a novel pathway to understand and tackle antimicrobial resistance more effectively.

Assessment of ibrutinib scheduling on leukocyte, lymph node size and blood pressure dynamics in chronic lymphocytic leukemia through pharmacokinetic-pharmacodynamic models

Ibrutinib is a Bruton tyrosine kinase (Btk) inhibitor for treating chronic lymphocytic leukemia (CLL). It has also been associated with hypertension. The optimal dosing schedule for mitigating this adverse effect is currently under discussion. A quantification of relationships between systemic ibrutinib exposure and efficacy (i.e., leukocyte count and sum of the product of perpendicular diameters [SPD] of lymph nodes) and hypertension toxicity (i.e., blood pressure), and their association with overall survival is needed. Here, we present a semi-mechanistic pharmacokinetic-pharmacodynamic modeling framework to characterize such relationships and facilitate dose optimization. Data from a phase Ib/II study were used, including ibrutinib plasma concentrations to derive daily 0-24-h area under the concentration-time curve, leukocyte count, SPD, survival, and blood pressure measurements. A nonlinear mixed effects modeling approach was applied, considering ibrutinib's pharmacological action and CLL cell dynamics. The final framework included (i) an integrated model for SPD and leukocytes consisting of four CLL cell subpopulations with ibrutinib inhibiting phosphorylated Btk production, (ii) a turnover model in which ibrutinib stimulates an increase in blood pressure, and (iii) a competing risk model for dropout and death. Simulations predicted that the approved dosing schedule had a slightly higher efficacy (24-month, progression-free survival [PFS] 98%) than de-escalation schedules (24-month, average PFS ≈ 97%); the latter had, on average, ≈20% lower proportions of patients with hypertension. The developed modeling framework offers an improved understanding of the relationships among ibrutinib exposure, efficacy and toxicity biomarkers. This framework can serve as a platform to assess dosing schedules in a biologically plausible manner.

Optimization of blood pressure measurement practices for pharmacodynamic analyses of tyrosine-kinase inhibitors

Blood pressure measurements form a critical component of adverse event monitoring for tyrosine kinase inhibitors, but might also serve as a biomarker for dose titrations. This study explored the impact of various sources of within-individual variation on blood pressure readings to improve measurement practices and evaluated the utility for individual- and population-level dose selection. A pharmacokinetic-pharmacodynamic modeling framework was created to describe circadian blood pressure changes, inter- and intra-day variability, changes from dipper to non-dipper profiles, and the relationship between drug exposure and blood pressure changes over time. The framework was used to quantitatively evaluate the influence of physiological and pharmacological aspects on blood pressure measurements, as well as to compare measurement techniques, including office-based, home-based, and ambulatory 24-h blood pressure readings. Circadian changes, as well as random intra-day and inter-day variability, were found to be the largest sources of within-individual variation in blood pressure. Office-based and ambulatory 24-h measurements gave rise to potential bias (>5 mmHg), which was mitigated by model-based estimations. Our findings suggest that 5-8 consecutive, home-based, measurements taken at a consistent time around noon, or alternatively within a limited time frame (e.g., 8.00 a.m. to 12.00 p.m. or 12.00 p.m. to 5.00 p.m.), will give rise to the most consistent blood pressure estimates. Blood pressure measurements likely do not represent a sufficiently accurate method for individual-level dose selection, but may be valuable for population-level dose identification. A user-friendly tool has been made available to allow for interactive blood pressure simulations and estimations for the investigated scenarios.

Multistate model for pharmacometric analyses of overall survival in HER2-negative breast cancer patients treated with docetaxel

The aim of this study was to develop a multistate model for overall survival (OS) analysis, based on parametric hazard functions and combined with an investigation of predictors derived from a longitudinal tumor size model on the transition hazards. Different states - stable disease, tumor response, progression, second-line treatment, and death following docetaxel treatment initiation (stable state) in patients with HER2-negative breast cancer (n = 183) were used in model building. Past changes in tumor size prospectively predicts the probability of state changes. The hazard of death after progression was lower for subjects who had longer treatment response (i.e., longer time-to-progression). Young age increased the probability of receiving second-line treatment. The developed multistate model adequately described the transitions between different states and jointly the overall event and survival data. The multistate model allows for simultaneous estimation of transition rates along with their tumor model derived metrics. The metrics were evaluated in a prospective manner so not to cause immortal time bias. Investigation of predictors and characterization of the time to develop response, the duration of response, the progression-free survival, and the OS can be performed in a single multistate modeling exercise. This modeling approach can be applied to other cancer types and therapies to provide a better understanding of efficacy of drug and characterizing different states, thereby facilitating early clinical interventions to improve anticancer therapy.

Relevance of primary lesion location, tumour heterogeneity and genetic mutation demonstrated through tumour growth inhibition and overall survival modelling in metastatic colorectal cancer

AIMS

The aims of this work were to build a semi-mechanistic tumour growth inhibition (TGI) model for metastatic colorectal cancer (mCRC) patients receiving either cetuximab + chemotherapy or chemotherapy alone and to identify early predictors of overall survival (OS).

METHODS

A total of 1716 patients from 4 mCRC clinical studies were included in the analysis. The TGI model was built with 8973 tumour size measurements where the probability of drop-out was also included and modelled as a time-to-event variable using parametric survival models, as it was the case in the OS analysis. The effects of patient- and tumour-related covariates on model parameters were explored.

RESULTS

Chemotherapy and cetuximab effects were included in an additive form in the TGI model. Development of resistance was found to be faster for chemotherapy (drug effect halved at wk 8) compared to cetuximab (drug effect halved at wk 12). KRAS wild-type status and presenting a right-sided primary lesion were related to a 3.5-fold increase in cetuximab drug effect and a 4.7× larger cetuximab resistance, respectively. The early appearance of a new lesion (HR = 4.14), a large tumour size at baseline (HR = 1.62) and tumour heterogeneity (HR = 1.36) were the main predictors of OS.

CONCLUSIONS

Semi-mechanistic TGI and OS models have been developed in a large population of mCRC patients receiving chemotherapy in combination or not with cetuximab. Tumour-related predictors, including a machine learning derived-index of tumour heterogeneity, were linked to changes in drug effect, resistance to treatment or OS, contributing to the understanding of the variability in clinical response.

Pharmacodynamic analysis of hypertension caused by lenvatinib using real-world postmarketing surveillance data

Lenvatinib is a tyrosine kinase inhibitor of the vascular endothelial growth factor receptor used against nonoperative thyroid cancer; however, hypertension is a major dose-limiting side effect. In this study, hypertension caused by lenvatinib was described through a novel population pharmacodynamic model using postmarketing surveillance data obtained in Japan. The model consists of two maximum effect model components based on the (1) concentration of lenvatinib in plasma and (2) cumulative area under the curve of lenvatinib. In addition, antihypertensive drug of either an angiotensin-converting enzyme inhibitor/angiotensin receptor blocker or calcium channel blocker accounted for by lowering effect on diastolic blood pressure. Based on virtual simulations, the combination of antihypertensive drug and dose adjustment of lenvatinib showed a reduction in the probability of grade greater than or equal to 3 hypertension. The present model provides useful guidance in managing hypertension during treatment with lenvatinib in the real-world setting.

Dose individualisation in oncology using chemotherapy-induced neutropenia: Example of docetaxel in non-small cell lung cancer patients

AIMS

Chemotherapy-induced neutropenia has been associated with an increase in overall survival in non-small cell lung cancer patients. Therefore, neutrophil counts could be an interesting biomarker for drug efficacy as well as linked directly to toxicity. For drugs where neutropenia is dose limiting, neutrophil counts might be used for monitoring drug effect and for dosing optimisation.

METHODS

The relationship between drug effect on the first cycle neutrophil counts and patient survival was explored in a Phase III clinical trial where patients with non-small cell lung cancer were treated with docetaxel. Once the association has been established, dosing optimisation was performed for patients with severe toxicities (neutropenia) without compromising drug efficacy (overall survival).

RESULTS

After taking into account baseline prognostic factors, such as Eastern Cooperative Oncology Group performance status, smoking status, liver metastasis, tumour burden, neutrophil counts and albumin levels, a model-predicted drug effect on the first cycle neutrophil counts was strongly associated with patient survival (P = .005). Utilising this relationship in a dose optimisation algorithm, 194 out of 366 patients would have benefited from a dose reduction after the first cycle of docetaxel. The simulated 1-year survival probabilities associated with the original dose and the individualised dose were not different.

CONCLUSION

The strong relationship between drug effect on the first cycle neutrophil counts and patient survival suggests that this variable could be used to individualise dosing, possibly without needing pharmacokinetic samples. The algorithm highlights that doses could be reduced in case of severe haematological toxicities without compromising drug efficacy.

PK/PD modeling based on NO-ET homeostasis for improving management of sunitinib-induced hypertension in rats

Sunitinib is an oral small molecule multitargeted tyrosine kinase inhibitor, which is currently used to treat severe cancers. Clinical research has shown that patients treated with sunitinib develop hypertension. As soon as sunitinib-induced hypertension appears, it is usual to administer anti-hypertension agent. But this treatment may cause acute blood pressure fluctuation which may lead to additional cardiovascular risk. The aim of this study is to establish a mathematical model for managing sunitinib-induced hypertension and blood pressure fluctuation. A mechanism-based PK/PD model was developed based on animal experiments. Then this model was used to perform simulations, thus to propose an anti-hypertension indication, according to which the anti-hypertension treatment might yield relative low-level AUC and fluctuation of blood pressure. The simulation results suggest that the anti-hypertension agent may yield low-level AUC and fluctuation of blood pressure when relative ET-1 level ranges from −15% to 5% and relative NO level is more than 10% compared to control group. Finally, animal experiments were conducted to verify the simulation results. Macitentan (30 mg/kg) was administered based on the above anti-hypertension indication. Compared with the untreated group, the optimized treatment significantly reduced the AUC of blood pressure; meanwhile the fluctuation of blood pressure in optimized treatment group was 70% less than that in immediate treatment group. This work provides a novel model with potential translational value for managing sunitinib-induced hypertension.

Precision Dosing Priority Criteria: Drug, Disease, and Patient Population Variables

The administered dose of a drug modulates whether patients will experience optimal effectiveness, toxicity including death, or no effect at all. Dosing is particularly important for diseases and/or drugs where the drug can decrease severe morbidity or prolong life. Likewise, dosing is important where the drug can cause death or severe morbidity. Since we believe there are many examples where more precise dosing could benefit patients, it is worthwhile to consider how to prioritize drug-disease targets. One key consideration is the quality of information available from which more precise dosing recommendations can be constructed. When a new more precise dosing scheme is created and differs significantly from the approved label, it is important to consider the level of proof necessary to either change the label and/or change clinical practice. The cost and effort needed to provide this proof should also be considered in prioritizing drug-disease precision dosing targets. Although precision dosing is being promoted and has great promise, it is underutilized in many drugs and disease states. Therefore, we believe it is important to consider how more precise dosing is going to be delivered to high priority patients in a timely manner. If better dosing schemes do not change clinical practice resulting in better patient outcomes, then what is the use? This review paper discusses variables to consider when prioritizing precision dosing candidates while highlighting key examples of precision dosing that have been successfully used to improve patient care.

Model-Based Biomarker Selection for Dose Individualization of Tyrosine-Kinase Inhibitors

Tyrosine-kinase inhibitors (TKIs) demonstrate high inter-individual variability with respect to safety and efficacy and would therefore benefit from dose or schedule adjustments. This study investigated the efficacy, safety, and economical aspects of alternative dosing options for sunitinib in gastro-intestinal stromal tumors (GIST) and axitinib in metastatic renal cell carcinoma (mRCC). Dose individualization based on drug concentration, adverse effects, and sVEGFR-3 was explored using a modeling framework connecting pharmacokinetic and pharmacodynamic models, as well as overall survival. Model-based simulations were performed to investigate four different scenarios: (I) the predicted value of high-dose pulsatile schedules to improve clinical outcomes as compared to regular daily dosing, (II) the potential of biomarkers for dose individualizations, such as drug concentrations, toxicity measurements, and the biomarker sVEGFR-3, (III) the cost-effectiveness of biomarker-guided dose-individualizations, and (IV) model-based dosing approaches versus standard sample-based methods to guide dose adjustments in clinical practice. Simulations from the axitinib and sunitinib frameworks suggest that weekly or once every two weeks high-dosing result in lower overall survival in patients with mRCC and GIST, compared to continuous daily dosing. Moreover, sVEGFR-3 appears a safe and cost-effective biomarker to guide dose adjustments and improve overall survival (€36 784.- per QALY). Model-based estimations were for biomarkers in general found to correctly predict dose adjustments similar to or more accurately than single clinical measurements and might therefore guide dose adjustments. A simulation framework represents a rapid and resource saving method to explore various propositions for dose and schedule adjustments of TKIs, while accounting for complicating factors such as circulating biomarker dynamics and inter-or intra-individual variability.

Tumor Size and Overall Survival in Patients With Platinum-Resistant Ovarian Cancer Treated With Chemotherapy and Bevacizumab

Introduction:

Ovarian cancer is now recognized as a constellation of distinct subtypes of neoplasia involving the ovary and related structures. As a consequence of this heterogeneity, the analysis of covariates influencing the overall survival is crucial in this disease segment. In this work, an overall survival model incorporating tumor kinetics metrics in patients with platinum-resistant ovarian cancer was developed from the randomized, open label, phase 3 AURELIA trial.

Methods:

Tumor size data from 361 patients randomly allocated to the bevacizumab + chemotherapy or chemotherapy study arm were collected at baseline and every 8 to 9 weeks until disease progression. Patients continued to be followed for survival after treatment discontinuation. A landmarked Cox proportional hazard survival model was developed to characterize the overall survival distribution.

Results:

Two sets of factors were found to be influential on survival time: those describing the type and severity of disease (Eastern Cooperative Oncology Group [ECOG], Féderation Internationale de Gynécologie et d’Obstétrique [FIGO] stages, presence of ascites) and those summarizing the key features of the tumor kinetic model (tumor shrinkage at week 8 and tumor size at treatment onset). The treatment group was not required in the final model as the drug effect was accounted for in the tumor kinetics model.

Conclusions:

This work has identified both ascites and tumor kinetics metrics as being the 2 most influential factors to explain variability in overall survival in patients with platinum-resistant ovarian cancer.

Modeling Exposure-Driven Adverse Event Time Courses in Oncology Exemplified by Afatinib

Models were developed to characterize the relationship between afatinib exposure and diarrhea and rash/acne adverse event (AE) trajectories, and their predictive ability was assessed. Based on pooled data from seven phase II/III clinical studies including 998 patients, mixed-effects models for ordered categorical data were applied to describe daily AE severity. Clinical trial simulation aided by trial execution models was used for internal and external model evaluation. The final exposure-safety model consisted of longitudinal logistic regression models with first-order Markov elements for both AEs. Drug exposure was included as daily area under the concentration-time curve (AUC), and drug effects on the AEs were correlated. Clinical trial simulation allowed adequate prediction of maximum AE grades and AE severity time courses but overestimated the proportion of AE-dependent dose reductions and discontinuations. Both diarrhea and rash/acne were correlated with afatinib exposure. The developed modeling framework allows a prospective comparison of dosing strategies and study designs with respect to safety.

The impact of tacrolimus exposure on extrarenal adverse effects in adult renal transplant recipients

AIMS

Tacrolimus has been associated with notable extrarenal adverse effects (AEs), which are unpredictable and impact patient morbidity. The association between model-predicted tacrolimus exposure metrics and standardized extrarenal AEs in stable renal transplant recipients was investigated and a limited sampling strategy (LSS) was developed to predict steady-state tacrolimus area under the curve over a 12-h dosing period (AUC_ss,0-12h ).

METHODS

All recipients receiving tacrolimus and mycophenolic acid ≥6 months completed a 12-h cross-sectional observational pharmacokinetic-pharmacodynamic study. Patients were evaluated for the presence of individual and composite gastrointestinal, neurological, and aesthetic AEs during the study visit. The associations between AEs and tacrolimus exposure metrics generated from a published population pharmacokinetic model were investigated using a logistic regression analysis in NONMEM 7.3. An LSS was determined using a Bayesian estimation method with the same patients.

RESULTS

Dose-normalized tacrolimus AUC_ss,0-12h and apparent clearance were independently associated with diarrhoea, dyspepsia, insomnia and neurological AE ratio. Dose-normalized tacrolimus maximum concentration was significantly correlated with skin changes and acne. No AE associations were found with trough concentrations. Using limited sampling at 0, 2h; 0, 1, 4h; and 0, 1, 2, 4h provided a precise and unbiased prediction of tacrolimus AUC (root mean squared prediction error < 10%), which was not well characterized using trough concentrations only (root mean squared prediction error >15%).

CONCLUSIONS

Several AEs (i.e. diarrhoea, dyspepsia, insomnia and neurological AE ratio) were associated with tacrolimus dose normalized AUC_ss,0-12h and clearance. Skin changes and acne were associated with dose-normalized maximum concentrations. To facilitate clinical implementation, a LSS was developed to predict AUC_ss,0-12h values using sparse patient data to efficiently assess projected immunosuppressive exposure and potentially minimize AE manifestations.

A PK/PD Analysis of Circulating Biomarkers and Their Relationship to Tumor Response in Atezolizumab-Treated non-small Cell Lung Cancer Patients

To assess circulating biomarkers as predictors of antitumor response to atezolizumab (anti-programmed death-ligand 1 (PD-L1), Tecentriq) serum pharmacokinetic (PK) and 95 plasma biomarkers were analyzed in 88 patients with relapsed/refractory non-small cell lung cancer (NSCLC) receiving atezolizumab i.v. q3w (10-20 mg/kg) in the PCD4989g phase I clinical trial. Following exploratory analyses, two plasma biomarkers were chosen for further study and correlation with change in tumor size (the sum of the longest diameter) was assessed in a pharmacokinetic/pharmacodynamic (PK/PD) tumor modeling framework. When longitudinal kinetics of biomarkers and tumor size were modeled, tumor shrinkage was found to significantly correlate with area under the curve (AUC), baseline factors (metastatic sites, liver metastases, and smoking status), and relative change in interleukin (IL)-18 level from baseline at day 21 (RCFB_{IL -18,d21} ). Although AUC was a major predictor of tumor shrinkage, the effect was estimated to dissipate with an average half-life of 80 days, whereas RCFB_{IL -18,d21} seemed relevant to the duration of the response.

Mathematical modeling identifies optimum lapatinib dosing schedules for the treatment of glioblastoma patients

Human primary glioblastomas (GBM) often harbor mutations within the epidermal growth factor receptor (EGFR). Treatment of EGFR-mutant GBM cell lines with the EGFR/HER2 tyrosine kinase inhibitor lapatinib can effectively induce cell death in these models. However, EGFR inhibitors have shown little efficacy in the clinic, partly because of inappropriate dosing. Here, we developed a computational approach to model the in vitro cellular dynamics of the EGFR-mutant cell line SF268 in response to different lapatinib concentrations and dosing schedules. We then used this approach to identify an effective treatment strategy within the clinical toxicity limits of lapatinib, and developed a partial differential equation modeling approach to study the in vivo GBM treatment response by taking into account the heterogeneous and diffusive nature of the disease. Despite the inability of lapatinib to induce tumor regressions with a continuous daily schedule, our modeling approach consistently predicts that continuous dosing remains the best clinically feasible strategy for slowing down tumor growth and lowering overall tumor burden, compared to pulsatile schedules currently known to be tolerated, even when considering drug resistance, reduced lapatinib tumor concentrations due to the blood brain barrier, and the phenotypic switch from proliferative to migratory cell phenotypes that occurs in hypoxic microenvironments. Our mathematical modeling and statistical analysis platform provides a rational method for comparing treatment schedules in search for optimal dosing strategies for glioblastoma and other cancer types.

Exposure-Response Analysis of Necitumumab Efficacy in Squamous Non-Small Cell Lung Cancer Patients

We sought to describe the exposure-response relationship of necitumumab efficacy in squamous non-small cell lung cancer patients and evaluate intrinsic and extrinsic patient descriptors that may guide dosing. SQUIRE was a phase III study comparing necitumumab in combination with gemcitabine and cisplatin vs. gemcitabine and cisplatin alone in 1,014 patients. An integrated model for tumor size dynamics and overall survival was developed, where reduction in tumor size results in a decrease in survival hazard. The change in tumor size was characterized using linear growth and first-order shrinkage. Overall survival was described using a combination of a Weibull function and Gompertz function for the hazard, with dynamic tumor size being a predictor for the hazard. Although body weight resulted in higher clearance and lower exposure, simulations showed that an 800 mg flat dose provided optimal response regardless of body weight.

A Minimal Continuous-Time Markov Pharmacometric Model

In this work, an alternative model to discrete-time Markov model (DTMM) or standard continuous-time Markov model (CTMM) for analyzing ordered categorical data with Markov properties is presented: the minimal CTMM (mCTMM). Through a CTMM reparameterization and under the assumption that the transition rate between two consecutive states is independent on the state, the Markov property is expressed through a single parameter, the mean equilibration time, and the steady-state probabilities are described by a proportional odds (PO) model. The mCTMM performance was evaluated and compared to the PO model (ignoring Markov features) and to published Markov models using three real data examples: the four-state fatigue and hand-foot syndrome data in cancer patients initially described by DTMM and the 11-state Likert pain score data in diabetic patients previously analyzed with a count model including Markovian transition probability inflation. The mCTMM better described the data than the PO model, and adequately predicted the average number of transitions per patient and the maximum achieved scores in all examples. As expected, mCTMM could not describe the data as well as more flexible DTMM but required fewer estimated parameters. The mCTMM better fitted Likert data than the count model. The mCTMM enables to explore the effect of potential predictive factors such as drug exposure and covariates, on ordered categorical data, while accounting for Markov features, in cases where DTMM and/or standard CTMM is not applicable or conveniently implemented, e.g., non-uniform time intervals between observations or large number of categories.

A Pharmacometric Framework for Axitinib Exposure, Efficacy, and Safety in Metastatic Renal Cell Carcinoma Patients

The relationships between exposure, biomarkers (vascular endothelial growth factor (VEGF), soluble VEGF receptors (sVEGFR)-1, -2, -3, and soluble stem cell factor receptor (sKIT)), tumor sum of longest diameters (SLD), diastolic blood pressure (dBP), and overall survival (OS) were investigated in a modeling framework. The dataset included 64 metastatic renal cell carcinoma patients (mRCC) treated with oral axitinib. Biomarker timecourses were described by indirect response (IDR) models where axitinib inhibits sVEGFR-1, -2, and -3 production, and VEGF degradation. No effect was identified on sKIT. A tumor model using sVEGFR-3 dynamics as driver predicted SLD data well. An IDR model, with axitinib exposure stimulating the response, characterized dBP increase. In a time-to-event model the SLD timecourse predicted OS better than exposure, biomarker- or dBP-related metrics. This type of framework can be used to relate pharmacokinetics, efficacy, and safety to long-term clinical outcome in mRCC patients treated with VEGFR inhibitors. (ClinicalTrial.gov identifier NCT00569946.).

A Review of Modeling Approaches to Predict Drug Response in Clinical Oncology

Model-based approaches have emerged as important tools for quantitatively understanding temporal relationships between drug dose, concentration, and effect over the course of treatment, and have now become central to optimal drug development and tailored drug treatment. In oncology, the therapeutic index of a chemotherapeutic drug is typically narrow and a full dose-response relationship is not available, often because of treatment failure. Noting the benefits of model-based approaches and the low therapeutic index of oncology drugs, in recent years, modeling approaches have been increasingly used to streamline oncologic drug development through early identification and quantification of dose-response relationships. With this background, this report reviews publications that used model-based approaches to evaluate drug treatment outcome variables in oncology therapeutics, ranging from tumor size dynamics to tumor/biomarker time courses and survival response.

Bridging Sunitinib Exposure to Time-to-Tumor Progression in Hepatocellular Carcinoma Patients With Mathematical Modeling of an Angiogenic Biomarker

Hepatocellular carcinoma (HCC) is third in cancer-related causes of death worldwide and its treatment is a significant unmet medical need. Sunitinib is a selective tyrosine kinase inhibitor of the angiogenic biomarker: soluble vascular endothelial growth factor receptor-2 (sVEGFR2 ). Sunitinib failed its primary overall survival endpoint in patients with advanced HCC in a phase III trial compared to sorafenib. In the present study, pharmacokinetic-pharmacodynamic modeling was used to link drug-exposure to tumor-growth-inhibition (TGI) and time-to-tumor progression (TTP) through sVEGFR2 dynamics. The results suggest that 1) active drug concentration (i.e., sunitinib and its metabolite) inhibits the release of sVEGFR2 and that such inhibition is associated with TGI, and 2) daily sVEGFR2 exposure is likely a reliable predictor for the TTP in HCC patients. Moreover, the model quantitatively links the dynamics of an angiogenesis biomarker to TTP and accurately predicts observed literature-reported results of placebo treatment.

PK-PD modeling of individual lesion FDG-PET response to predict overall survival in patients with sunitinib-treated gastrointestinal stromal tumor

Pharmacometric models were developed to characterize the relationships between lesion-level tumor metabolic activity, as assessed by the maximum standardized uptake value (SUVmax) obtained on [(18)F]-fluorodeoxyglucose (FDG) positron emission tomography (PET), tumor size, and overall survival (OS) in 66 patients with gastrointestinal stromal tumor (GIST) treated with intermittent sunitinib. An indirect response model in which sunitinib stimulates tumor loss best described the typically rapid decrease in SUVmax during on-treatment periods and the recovery during off-treatment periods. Substantial interindividual and interlesion variability were identified in SUVmax baseline and drug sensitivity. A parametric time-to-event model identified the relative change in SUVmax at one week for the lesion with the most pronounced response as a better predictor of OS than tumor size. Based on the proposed modeling framework, early changes in FDG-PET response may serve as predictor for long-term outcome in sunitinib-treated GIST.

Pharmacodynamic modeling of adverse effects of anti-cancer drug treatment

Purpose

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

Methods

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

Results

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

Conclusions

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

Establishing the Quantitative Relationship Between Lanreotide Autogel®, Chromogranin A, and Progression-Free Survival in Patients with Nonfunctioning Gastroenteropancreatic Neuroendocrine Tumors

The objective of this work was to establish the quantitative relationship between Lanreotide Autogel® (LAN) on serum chromogranin A (CgA) and progression-free survival (PFS) in patients with nonfunctioning gastroenteropancreatic neuroendocrine tumors (GEP-NETs) through an integrated pharmacokinetic/pharmacodynamic (PK/PD) model. In CLARINET, a phase III, randomized, double-blind, placebo-controlled study, 204 patients received deep subcutaneous injections of LAN 120 mg (n = 101) or placebo (n = 103) every 4 weeks for 96 weeks. Data for 810 LAN and 1298 CgA serum samples (n = 632 placebo and n = 666 LAN) were used to develop a parametric time-to-event model to relate CgA levels and PFS (76 patients experienced disease progression: n = 49 placebo and n = 27 LAN). LAN serum profiles were described by a one-compartment disposition model. Absorption was characterized by two parallel pathways following first- and zero-order kinetics. As PFS data were considered informative dropouts, CgA and PFS responses were modeled jointly. The LAN-induced decrease in CgA levels was described by an inhibitory E MAX model. Patient age and target lesions at baseline were associated with an increment in baseline CgA. Weibull model distribution showed that decreases in CgA from baseline reduced the hazard of disease progression significantly (P < 0.001). Covariates of tumor location in the pancreas and tumor hepatic tumor load were associated with worse prognosis (P < 0.001). We established a semimechanistic PK/PD model to better understand the effect of LAN on a surrogate endpoint (serum CgA) and ultimately the clinical endpoint (PFS) in treatment-naive patients with nonfunctioning GEP-NETs.

Population pharmacokinetic-pharmacodynamic modelling in oncology: a tool for predicting clinical response

In oncology trials, overall survival (OS) is considered the most reliable and preferred endpoint to evaluate the benefit of drug treatment. Other relevant variables are also collected from patients for a given drug and its indication, and it is important to characterize the dynamic effects and links between these variables in order to improve the speed and efficiency of clinical oncology drug development. However, the drug-induced effects and causal relationships are often difficult to interpret because of temporal differences. To address this, population pharmacokinetic–pharmacodynamic (PKPD) modelling and parametric time-to-event (TTE) models are becoming more frequently applied. Population PKPD and TTE models allow for exploration towards describing the data, understanding the disease and drug action over time, investigating relevance of biomarkers, quantifying patient variability and in designing successful trials. In addition, development of models characterizing both desired and adverse effects in a modelling framework support exploration of risk-benefit of different dosing schedules. In this review, we have summarized population PKPD modelling analyses describing tumour, tumour marker and biomarker responses, as well as adverse effects, from anticancer drug treatment data. Various model-based metrics used to drive PD response and predict OS for oncology drugs and their indications are also discussed.

Bringing Model-Based Prediction to Oncology Clinical Practice: A Review of Pharmacometrics Principles and Applications

Despite much investment and progress, oncology is still an area with significant unmet medical needs, with new therapies and more effective use of current therapies needed. The emergent field of pharmacometrics combines principles from pharmacology (pharmacokinetics [PK] and pharmacodynamics [PD]), statistics, and computational modeling to support drug development and optimize the use of already marketed drugs. Although it has gained a role within drug development, its use in clinical practice remains scarce. The aim of the present study was to review the principal pharmacometric concepts and provide some examples of its use in oncology. Integrated population PK/PD/disease progression models as part of the pharmacometrics platform provide a powerful tool to predict outcomes so that the right dose can be given to the right patient to maximize drug efficacy and reduce drug toxicity. Population models often can be developed with routinely collected medical record data; therefore, we encourage the application of such models in the clinical setting by generating close collaborations between physicians and pharmacometricians.

IMPLICATIONS FOR PRACTICE

The present review details how the emerging field of pharmacometrics can integrate medical record data with predictive pharmacological and statistical models of drug response to optimize and individualize therapies. In order to make this routine practice in the clinic, greater awareness of the potential benefits of the field is required among clinicians, together with closer collaboration between pharmacometricians and clinicians to ensure the requisite data are collected in a suitable format for pharmacometrics analysis.

Population pharmacokinetic-pharmacodynamic modelling of 24-h diastolic ambulatory blood pressure changes mediated by axitinib in patients with metastatic renal cell carcinoma

BACKGROUND

Increased blood pressure (BP) is commonly observed in patients treated with vascular endothelial growth factor pathway inhibitors, including axitinib. Ambulatory BP monitoring (ABPM) and pharmacokinetic data were collected in a randomised, double-blind phase II study of axitinib with or without dose titration in previously untreated patients with metastatic renal cell carcinoma.

OBJECTIVE

Aims of these analyses were to (1) develop a population pharmacokinetic-pharmacodynamic model for describing the relationship between axitinib exposure and changes in diastolic BP (dBP) and (2) simulate changes in dBP with different axitinib dosing regimens.

METHODS

We employed a three-stage modelling approach, which included development of (1) a baseline 24-h ABPM model, (2) a pharmacokinetic model from serial and sparse pharmacokinetic data, and (3) an indirect-response, maximum-effect (Emax) model to evaluate the exposure-driven effect of axitinib on dBP. Simulations (N = 1,000) were performed using the final pharmacokinetic-pharmacodynamic model to evaluate dBP changes on days 4 and 15 of treatment with different axitinib doses.

RESULTS

Baseline ABPM data from 62 patients were best described by 24-h mean dBP and two cosine terms. The final indirect-response Emax model showed good agreement between observed 24-h ABPM data and population and individual predictions. The maximum increase in dBP was 20.8 %, and the axitinib concentration at which 50 % of the maximal increase in dBP was reached was 12.4 ng/mL.

CONCLUSION

Our model adequately describes the relationship between axitinib exposure and dBP increases. Results from these analyses may potentially be applied to infer dBP changes in patients administered axitinib at nonstandard doses.

A Modeling and Simulation Framework for Adverse Events in Erlotinib-Treated Non-Small-Cell Lung Cancer Patients

Treatment with erlotinib, an epidermal growth factor receptor tyrosine kinase inhibitor used for treating non-small-cell lung cancer (NSCLC) and other cancers, is frequently associated with adverse events (AE). We present a modeling and simulation framework for the most common erlotinib-induced AE, rash, and diarrhea, providing insights into erlotinib toxicity. We used the framework to investigate the safety of high-dose erlotinib pulses proposed to limit acquired resistance while treating NSCLC. Continuous-time Markov models were developed using rash and diarrhea AE data from 39 NSCLC patients treated with erlotinib (150 mg/day). Exposure and different covariates were investigated as predictors of variability. Rash was also tested as a survival predictor. Models developed were used in a simulation analysis to compare the toxicities of different regimens, including the previously mentioned pulsed strategy. Probabilities of experiencing rash or diarrhea were found to be highest early during treatment. Rash, but not diarrhea, was positively correlated with erlotinib exposure. In contrast with some common understandings, radiotherapy decreased transitioning to higher rash grades by 81% (p < 0.01), and experiencing rash was not correlated with positive survival outcomes. Model simulations predicted that the proposed pulsed regimen (1600 mg/week + 50 mg/day remaining week days) results in a maximum of 20% of the patients suffering from severe rash throughout the treatment course in comparison to 12% when treated with standard dosing (150 mg/day). In conclusion, the framework demonstrated that radiotherapy attenuates erlotinib-induced rash, providing an opportunity to use radiotherapy and erlotinib together, and demonstrated the tolerability of high-dose pulses intended to address acquired resistance to erlotinib.

Mathematical models of the steps involved in the systemic delivery of a chemotherapeutic to a solid tumor: From circulation to survival

The efficacy of an intravenously administered chemotherapeutic for treatment of a solid tumor is dependent on a sequence of steps, including circulation, extravasation by the enhanced permeability and retention effect, transport in the tumor microenvironment, the mechanism of cellular uptake and trafficking, and the mechanism of drug action. These steps are coupled since the time dependent concentration in circulation determines the concentration and distribution in the tumor microenvironment, and hence the amount taken up by individual cells within the tumor. Models have been developed for each of the steps in the delivery process although their predictive power remains limited. Advances in our understanding of the steps in the delivery process will result in refined models with improvements in predictive power and ultimately allow the development of integrated models that link systemic administration of a drug to the probability of survival. Integrated models that predict outcomes based on patient specific data could be used to select the optimum therapeutic regimens. Here we present an overview of current models for the steps in the delivery process and highlight knowledge gaps that are key to developing integrated models.

Antiangiogenic agents in gynecological cancer: State of art and perspectives of clinical research

Vascular endothelial growth factor [VEGF] pathway, which plays a key role in angiogenesis, may be blocked by either extracellular interference with VEGF itself (bevacizumab [BEV] or aflibercept), or intracytoplasmic inhibition of VEGF receptor (pazopanib, nintedanib, cediranid, sunitinib and sorafenib). An alternative approach is represented by trebananib, a fusion protein that prevents the interaction of angiopoietin [Ang]-1 and Ang-2 with Tie2 receptor on vascular endothelium. The combination of antiangiogenic agents, especially BEV, and chemotherapy is a rational therapeutic option for primary or recurrent ovarian carcinoma. However, it will be difficult to accept that it represents the new standard treatment, until biological characterization of ovarian carcinoma has not identified subsets of tumors with different responsiveness to BEV. Anti-angiogenesis is an interesting target also for recurrent cervical or endometrial cancer, but nowadays the use of anti-angiogenic agents in these malignancies should be reserved to patients enrolled in clinical trials.

Early Prediction of Disease Progression in Small Cell Lung Cancer: Toward Model-Based Personalized Medicine in Oncology

Predictive biomarkers can play a key role in individualized disease monitoring. Unfortunately, the use of biomarkers in clinical settings has thus far been limited. We have previously shown that mechanism-based pharmacokinetic/pharmacodynamic modeling enables integration of nonvalidated biomarker data to provide predictive model-based biomarkers for response classification. The biomarker model we developed incorporates an underlying latent variable (disease) representing (unobserved) tumor size dynamics, which is assumed to drive biomarker production and to be influenced by exposure to treatment. Here, we show that by integrating CT scan data, the population model can be expanded to include patient outcome. Moreover, we show that in conjunction with routine medical monitoring data, the population model can support accurate individual predictions of outcome. Our combined model predicts that a change in disease of 29.2% (relative standard error 20%) between two consecutive CT scans (i.e., 6-8 weeks) gives a probability of disease progression of 50%. We apply this framework to an external dataset containing biomarker data from 22 small cell lung cancer patients (four patients progressing during follow-up). Using only data up until the end of treatment (a total of 137 lactate dehydrogenase and 77 neuron-specific enolase observations), the statistical framework prospectively identified 75% of the individuals as having a predictable outcome in follow-up visits. This included two of the four patients who eventually progressed. In all identified individuals, the model-predicted outcomes matched the observed outcomes. This framework allows at risk patients to be identified early and therapeutic intervention/monitoring to be adjusted individually, which may improve overall patient survival.

Optimizing oncology therapeutics through quantitative translational and clinical pharmacology: challenges and opportunities

Despite advances in biomedical research that have deepened our understanding of cancer hallmarks, resulting in the discovery and development of targeted therapies, the success rates of oncology drug development remain low. Opportunities remain for objective dose selection informed by exposure-response understanding to optimize the benefit-risk balance of novel therapies for cancer patients. This review article discusses the principles and applications of modeling and simulation approaches across the lifecycle of development of oncology therapeutics. Illustrative examples are used to convey the value gained from integration of quantitative clinical pharmacology strategies from the preclinical-translational phase through confirmatory clinical evaluation of efficacy and safety.

A review of mixed-effects models of tumor growth and effects of anticancer drug treatment used in population analysis

Population modeling of tumor size dynamics has recently emerged as an important tool in pharmacometric research. A series of new mixed-effects models have been reported recently, and we present herein a synthetic view of models with published mathematical equations aimed at describing the dynamics of tumor size in cancer patients following anticancer drug treatment. This selection of models will constitute the basis for the Drug Disease Model Resources (DDMoRe) repository for models on oncology.

Association analysis of genetic polymorphisms in genes related to sunitinib pharmacokinetics, specifically clearance of sunitinib and SU12662

Interpatient variability in the pharmacokinetics (PK) of sunitinib is high. Single nucleotide polymorphisms (SNPs) in PK candidate genes have been associated with the efficacy and toxicity of sunitinib, but whether these SNPs truly affect the PK of sunitinib remains to be elucidated. This multicenter study involving 114 patients investigated whether these SNPs and haplotypes in genes encoding metabolizing enzymes or efflux transporters are associated with the clearance of sunitinib and its active metabolite SU12662. SNPs were tested as covariates in a population PK model. From univariate analysis, we found that the SNPs in CYP3A4, CYP3A5, and ABCB1 were associated with the clearance of both sunitinib and SU12662. In multivariate analysis, CYP3A4*22 was found to be eliminated last with an effect size of -22.5% on clearance. Observed effect sizes are below the interindividual variability in clearance and are therefore too limited to directly guide individual dosing of sunitinib.

PKPD Modeling of VEGF, sVEGFR-2, sVEGFR-3, and sKIT as Predictors of Tumor Dynamics and Overall Survival Following Sunitinib Treatment in GIST

The predictive value of longitudinal biomarker data (vascular endothelial growth factor (VEGF), soluble VEGF receptor (sVEGFR)-2, sVEGFR-3, and soluble stem cell factor receptor (sKIT)) for tumor response and survival was assessed based on data from 303 patients with imatinib-resistant gastrointestinal stromal tumors (GIST) receiving sunitinib and/or placebo treatment. The longitudinal tumor size data were well characterized by a tumor growth inhibition model, which included, as significant descriptors of tumor size change, the model-predicted relative changes from baseline over time for sKIT (most significant) and sVEGFR-3, in addition to sunitinib exposure. Survival time was best described by a parametric time-to-event model with baseline tumor size and relative change in sVEGFR-3 over time as predictive factors. Based on the proposed modeling framework to link longitudinal biomarker data with overall survival using pharmacokinetic-pharmacodynamic models, sVEGFR-3 demonstrated the greatest predictive potential for overall survival following sunitinib treatment in GIST.CPT: Pharmacometrics & Systems Pharmacology (2013) 2, e84; doi:10.1038/psp.2013.61; advance online publication 20 November 2013.