Evaluation of Tumor Size Response Metrics to Predict Survival in Oncology Clinical Trials

A bootstrapping method to optimize go/no-go decisions from single-arm, signal-finding studies in oncology

Phase Ib trials are common in oncology development but often are not powered for statistical significance. Go/no-go decisions are largely driven by observed trends in response data. We applied a bootstrapping method to systematically compare tumor dynamic end points to historical control data to identify drugs with clinically meaningful efficacy. A proprietary mathematical model calibrated to phase Ib anti-PD-1 therapy trial data (KEYNOTE-001) was used to simulate thousands of phase Ib trials (n = 30) with a combination of anti-PD-1 therapy and four novel agents with varying efficacy. A redacted bootstrapping method compared these results to a simulated phase III control arm (N = 511) while adjusting for differences in trial duration and cohort size to determine the probability that the novel agent provides clinically meaningful efficacy. Receiver operating characteristic (ROC) analysis showed strong ability to separate drugs with modest (area under ROC [AUROC] = 83%), moderate (AUROC = 96%), and considerable efficacy (AUROC = 99%) from placebo in early-phase trials (n = 30). The method was shown to effectively move drugs with a range of efficacy through an in silico pipeline with an overall success rate of 93% and false-positive rate of 7.5% from phase I to phase III. This model allows for effective comparisons of tumor dynamics from early clinical trials with more mature historical control data and provides a framework to predict drug efficacy in early-phase trials. We suggest this method should be employed to improve decision making in early oncology trials.

Bridging the gap between two-stage and joint models: The case of tumor growth inhibition and overall survival models

Many clinical trials generate both longitudinal biomarker and time-to-event data. We might be interested in their relationship, as in the case of tumor size and overall survival in oncology drug development. Many well-established methods exist for analyzing such data either sequentially (two-stage models) or simultaneously (joint models). Two-stage modeling (2stgM) has been challenged (i) for not acknowledging that biomarkers are endogenous covariable to the survival submodel and (ii) for not propagating the uncertainty of the longitudinal biomarker submodel to the survival submodel. On the other hand, joint modeling (JM), which properly circumvents both problems, has been criticized for being time-consuming, and difficult to use in practice. In this paper, we explore a third approach, referred to as a novel two-stage modeling (N2stgM). This strategy reduces the model complexity without compromising the parameter estimate accuracy. The three approaches (2stgM, JM, and N2stgM) are formulated, and a Bayesian framework is considered for their implementation. Both real and simulated data were used to analyze the performance of such approaches. In all scenarios, our proposal estimated the parameters approximately as JM but without being computationally expensive, while 2stgM produced biased results.

Utility and impact of quantitative pharmacology on dose selection and clinical development of immuno-oncology therapy

Immuno-oncology (IO) therapies have changed the cancer treatment landscape. Immune checkpoint inhibitors (ICIs) have improved overall survival in 20-40% of patients with malignancies that were previously refractory. Due to the uniqueness in biology, modalities and patient responses, drug development strategies for IO differed from that traditionally used for cytotoxic and target therapies in oncology, and quantitative pharmacology utilizing modeling approach can be applied in all phases of the development process. In this review, we used case studies to showcase how various modeling methodologies were applied from translational science and dose selection through to label change, using examples that included anti-programmed-death-1 (anti-PD-1), anti-programmed-death ligand-1 (anti-PD-L1), anti-cytotoxic T-lymphocyte-associated protein 4 (anti-CTLA-4), and anti-glucocorticoid-induced tumor necrosis factor receptor-related protein (anti-GITR) antibodies. How these approaches were utilized to support phase I-III dose selection, the design of phase III trials, and regulatory decisions on label change are discussed to illustrate development strategies. Model-based quantitative approaches have positively impacted IO drug development, and a better understanding of the biology and exposure-response relationship may benefit the development and optimization of new IO therapies.

Comparison of two-stage and joint TGI-OS modeling using data from six atezolizumab clinical studies in patients with metastatic non-small cell lung cancer

Two-stage and joint modeling approaches are the two main approaches to investigate the link between longitudinal tumor size data and overall survival (OS) and anticipate clinical trial outcome. We here used a large database composed of one phase II and five phase III clinical trials evaluating atezolizumab (an immunotherapy) in monotherapy or in combination with chemotherapies in 3699 patients with non-small cell lung cancer to evaluate the differences between both approaches in terms of parameter estimates, magnitude of covariate effects, and ability to predict OS. Although the two-stage approach may underestimate the magnitude of the impact of tumor growth rate (KG ) on OS compared to joint modeling approach (hazard ratios [HRs] of 0.42-2.52 vs. 0.25-2.85, respectively, for individual KG varying from the 5th and 95th percentiles), this difference did not lead into poorer performance of the two-stage approach to describe the OS distribution in the six clinical studies. Overall, two-stage and joint modeling approaches accurately predicted OS HR with a median (range) difference with the observed OS HR of 0.02 (0.01-0.18) and 0.03 (0.00-0.19), in all cases considered, respectively (e.g., for IMpower150: 0.80 [0.66-0.95] vs. 0.82 [0.70-0.95], respectively, whereas the observed OS HR was 0.80). In our setting, the two-stage approach accurately predicted the benefit of atezolizumab on OS. Further work is needed to verify if similar results are achieved using phase Ib or phase II clinical trials where the number of patients and measurements is limited as well as in other cancer indications.

Support to early clinical decisions in drug development and personalised medicine with checkpoint inhibitors using dynamic biomarker-overall survival models

Longitudinal models of biomarkers such as tumour size dynamics capture treatment efficacy and predict treatment outcome (overall survival) of a variety of anticancer therapies, including chemotherapies, targeted therapies, immunotherapies and their combinations. These pharmacological endpoints like tumour dynamic (tumour growth inhibition) metrics have been proposed as alternative endpoints to complement the classical RECIST endpoints (objective response rate, progression-free survival) to support early decisions both at the study level in drug development as well as at the patients level in personalised therapy with checkpoint inhibitors. This perspective paper presents recent developments and future directions to enable wider and robust use of model-based decision frameworks based on pharmacological endpoints.

Early Decision Making in a Randomized Phase II Trial of Atezolizumab in Biliary Tract Cancer Using a Tumor Growth Inhibition-Survival Modeling Framework

We assess the longitudinal tumor growth inhibition (TGI) metrics and overall survival (OS) predictions applied to patients with advanced biliary tract cancer (BTC) enrolled in IMbrave151 a multicenter randomized phase II, double-blind, placebo-controlled trial evaluating the efficacy and safety of atezolizumab with or without bevacizumab in combination with cisplatin plus gemcitabine. Tumor growth rate (KG) was estimated for patients in IMbrave151. A pre-existing TGI-OS model for patients with hepatocellular carcinoma in IMbrave150 was modified to include available IMbrave151 study covariates and KG estimates and used to simulate IMbrave151 study outcomes. At the interim progression-free survival (PFS) analysis (98 patients, 27 weeks follow-up), clear separation in tumor dynamic profiles with a faster shrinkage rate and slower KG (0.0103 vs. 0.0117 week-1 ; tumor doubling time 67 vs. 59 weeks; KG geometric mean ratio of 0.84) favoring the bevacizumab containing arm was observed. At the first interim analysis for PFS, the simulated OS hazard ratio (HR) 95% prediction interval (PI) of 0.74 (95% PI: 0.58-0.94) offered an early prediction of treatment benefit later confirmed at the final analysis, observed HR of 0.76 based on 159 treated patients and 34 weeks of follow-up. This is the first prospective application of a TGI-OS modeling framework supporting gating of a phase III trial. The findings demonstrate the utility for longitudinal TGI and KG geometric mean ratio as relevant end points in oncology studies to support go/no-go decision making and facilitate interpretation of the IMbrave151 results to support future development efforts for novel therapeutics for patients with advanced BTC.

Model-based prediction of progression-free survival for combination therapies in oncology

Progression-free survival (PFS) is an important clinical metric for comparing and evaluating similar treatments for the same disease within oncology. After the completion of a clinical trial, a descriptive analysis of the patients' PFS is often performed post hoc using the Kaplan-Meier estimator. However, to perform predictions, more sophisticated quantitative methods are needed. Tumor growth inhibition models are commonly used to describe and predict the dynamics of preclinical and clinical tumor size data. Moreover, frameworks also exist for describing the probability of different types of events, such as tumor metastasis or patient dropout. Combining these two types of models into a so-called joint model enables model-based prediction of PFS. In this paper, we have constructed a joint model from clinical data comparing the efficacy of FOLFOX against FOLFOX + panitumumab in patients with metastatic colorectal cancer. The nonlinear mixed effects framework was used to quantify interindividual variability (IIV). The model describes tumor size and PFS data well, and showed good predictive capabilities using truncated as well as external data. A machine-learning guided analysis was performed to reduce unexplained IIV by incorporating patient covariates. The model-based approach illustrated in this paper could be useful to help design clinical trials or to determine new promising drug candidates for combination therapy trials.

Comparison of sequential and joint nonlinear mixed effects modeling of tumor kinetics and survival following Durvalumab treatment in patients with metastatic urothelial carcinoma

Standard endpoints such as objective response rate are usually poorly correlated with overall survival (OS) for treatment with immune checkpoint inhibitors. Longitudinal tumor size may serve as a more useful predictor of OS, and establishing a quantitative relationship between tumor kinetics (TK) and OS is a crucial step for successfully predicting OS based on limited tumor size measurements. This study aims to develop a population TK model in combination with a parametric survival model by sequential and joint modeling approaches to characterize durvalumab phase I/II data from patients with metastatic urothelial cancer, and to evaluate and compare the performance of the two modeling approaches in terms of parameter estimates, TK and survival predictions, and covariate identification. The tumor growth rate constant was estimated to be greater for patients with OS ≤ 16 weeks as compared to that for patients with OS > 16 weeks with the joint modeling approach (k_g= 0.130 vs. 0.0551 week^-1, p-value < 0.0001), but similar for both groups (k_g = 0.0624 vs.0.0563 week^-1, p-value = 0.37) with the sequential modeling approach. The predicted TK profiles by joint modeling appeared better aligned with clinical observations. Joint modeling also predicted OS more accurately than the sequential approach according to concordance index and Brier score. The sequential and joint modeling approaches were also compared using additional simulated datasets, and survival was predicted better by joint modeling in the case of a strong association between TK and OS. In conclusion, joint modeling enabled the establishment of a robust association between TK and OS and may represent a better choice for parametric survival analyses over the sequential approach.

Bayesian forecasting of tumor size metrics and overall survival

The tumor size ratio (TSR), time-to-tumor growth (TTG), and tumor growth rate (kG) are frequently suggested as model-based predictors of overall survival (OS) for different types of tumors. When the tumor metrics are applied in forecasting of the outcome for individual patients at an early stage, the tumor data might be sparse resulting in imprecise prediction. This simulation study aimed to investigate how the tumor follow-up data and estimation approaches influence the accuracy in the tumor size metrics and the predicted hazard of death for individual patients. Longitudinal tumor size and OS data were simulated using tumor growth inhibition and Weibull distribution models, respectively. Based on the model and increasing measurement durations, the accuracy (defined as 80-125% of the simulated "true" value) in individual metrics and hazard was computed. TSR week 6 (TSRw6) accuracy was adequate for 91% of the patients when tumor size was measured up to 12 weeks. For TTG and kG metrics, the highest accuracy observed was lower (43 and 77%, respectively) and occurred later (42 and 60 weeks, respectively). The simultaneous (joint) and sequential estimation approaches resulted in similar accuracies, however, in general, the sequential approach where individual tumor size parameters are fixed, demonstrated inferior estimation properties. The TSRw6 and the model-predicted tumor time course (absolute or relative change) had better forecasting properties than TTG or kG. The population pharmacokinetic (PK) parameters and data approach performed similarly or better than the simultaneous approach and had a better accuracy in estimating individuals' hazard of death than the individual PK parameters method.

Tumor growth inhibition modeling to support the starting dose for dacomitinib

Dacomitinib is a second-generation, irreversible EGFR tyrosine kinase inhibitor for first-line treatment of patients with metastatic non-small cell lung cancer and EGFR-activating mutations. A high rate of dose reductions in the pivotal trial led to an observed inverse exposure-response (ER) relationship with the primary end points. Three ER models were developed to determine if the starting dose from the pivotal trial, 45 mg once daily (q.d.) dose, is appropriate: a longitudinal logistic regression model for adverse event-related dose changes, a Claret tumor growth inhibition (TGI) model, and a Cox model for progression-free survival (PFS) based on the TGI model predictions. This analysis included 266 patients taking dacomitinib with a starting dose of 45 mg (N = 250) or 30 mg (N = 16) q.d. The ER relationships with the time-varying exposure metrics, most recent maximum plasma concentration (C_max ) and average concentration (C_avg ) from the first dose, were established for the dose reduction and TGI models, respectively. The TGI model characterized the tumor inhibition over time with constant growth rate (k_L  = 0.0012 years^-1 ) and highly variable kill rate (k_D  = 1.002 years^-1 /[μg/L]^θcavg , coefficient of variation [CV] = 89%) and drug resistance (λ = 14.47 years^-1 , CV = 96%) leading to prolonged tumor shrinkage. The ER relationship was characterized using an exposure parameter with a power parameterization (θcavg = 0.454, p < 0.0001). The Cox model found that baseline tumor size (p = 0.0166) and week 8 tumor shrinkage rate (p = 0.0726) were the best predictors of PFS. Simulations of dose reductions and drug interruptions on tumor shrinkage over time showed greater and more prolonged tumor shrinkage with a starting dose of 45 mg q.d.

Aiding early clinical drug development by elucidation of the relationship between tumor growth inhibition and survival in relapsed/refractory multiple myeloma patients

Early prognosis of clinical efficacy is an urgent need for oncology drug development. Herein, we systemically examined the quantitative approach of tumor growth inhibition (TGI) and survival modeling in the space of relapsed and refractory multiple myeloma (MM), aiming to provide insights into clinical drug development. Longitudinal serum M-protein and progression-free survival (PFS) data from three phase III studies (N = 1367) across six treatment regimens and different patient populations were leveraged. The TGI model successfully described the longitudinal M-protein data in patients with MM. The tumor inhibition and growth parameters were found to vary as per each study, likely due to the patient population and treatment regimen difference. Based on a parametric time-to-event model for PFS, M-protein reduction at week 4 was identified as a significant prognostic factor for PFS across the three studies. Other factors, including Eastern Cooperative Oncology Group performance status, prior anti-myeloma therapeutics, and baseline serum ß2-microglobulin level, were correlated with PFS as well. In conclusion, patient disease characteristics (i.e., baseline tumor burden and treatment lines) were important determinants of tumor inhibition and PFS in MM patients. M-protein change at week 4 was an early prognostic biomarker for PFS.

Deep learning for the prediction of early on-treatment response in metastatic colorectal cancer from serial medical imaging

In current clinical practice, tumor response assessment is usually based on tumor size change on serial computerized tomography (CT) scan images. However, evaluation of tumor response to anti-vascular endothelial growth factor therapies in metastatic colorectal cancer (mCRC) is limited because morphological change in tumor may occur earlier than tumor size change. Here we present an analysis utilizing a deep learning (DL) network to characterize tumor morphological change for response assessment in mCRC patients. We retrospectively analyzed 1,028 mCRC patients who were prospectively included in the VELOUR trial (NCT00561470). We found that DL network was able to predict early on-treatment response in mCRC and showed better performance than its size-based counterpart with C-Index: 0.649 (95% CI: 0.619,0.679) vs. 0.627 (95% CI: 0.567,0.638), p = 0.009, z-test. The integration of DL network with size-based methodology could further improve the prediction performance to C-Index: 0.694 (95% CI: 0.661,0.720), which was superior to size/DL-based-only models (all p < 0.001, z-test). Our study suggests that DL network could provide a noninvasive mean for quantitative and comprehensive characterization of tumor morphological change, which may potentially benefit personalized early on-treatment decision making.

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.

Prediction of overall survival in patients across solid tumors following atezolizumab treatments: A tumor growth inhibition-overall survival modeling framework

The objectives of the study were to use tumor size data from 10 phase II/III atezolizumab studies across five solid tumor types to estimate tumor growth inhibition (TGI) metrics and assess the impact of TGI metrics and baseline prognostic factors on overall survival (OS) for each tumor type. TGI metrics were estimated from biexponential models and posttreatment longitudinal data of 6699 patients. TGI‐OS full models were built using parametric survival regression by including all significant baseline covariates from the Cox univariate analysis followed by a backward elimination step. The model performance was evaluated for each trial by 1000 simulations of the OS distributions and hazard ratios (HR) of the atezolizumab‐containing arms versus the respective controls. The tumor growth rate estimate was the most significant predictor of OS across all tumor types. Several baseline prognostic factors, such as inflammatory status (C‐reactive protein, albumin, and/or neutrophil‐to‐lymphocyte ratio), tumor burden (sum of longest diameters, number of metastatic sites, and/or presence of liver metastases), Eastern Cooperative Oncology Group performance status, and lactate dehydrogenase were also highly significant across multiple studies in the final multivariate models. TGI‐OS models adequately described the OS distribution. The model‐predicted HRs indicated good model performance across the 10 studies, with observed HRs within the 95% prediction intervals for all study arms versus controls. Multivariate TGI‐OS models developed for different solid tumor types were able to predict treatment effect with various atezolizumab monotherapy or combination regimens and could be used to support design and analysis of future studies.

A meta-analysis of tumour response and relapse kinetics based on 34,881 patients: A question of cancer type, treatment and line of treatment

PURPOSE

Cancer disease burden is commonly assessed radiologically in solid tumours in support of response assessment via the RECIST criteria. These longitudinal data are amenable to mathematical modelling and these models characterise the initial tumour size, initial tumour shrinkage in responding patients and rate of regrowth as patient's disease progresses. Knowing how these parameters vary between patient populations and treatments would inform translational modelling approaches from non-clinical data as well as clinical trial design.

EXPERIMENTAL DESIGN

Here a meta-analysis of reported model parameter values is reported. Appropriate literature was identified via a PubMed search and the application of text-based clustering approaches. The resulting parameter estimates are examined graphically and with ANOVA.

RESULTS

Parameter values from a total of 80 treatment arms were identified based on 80 trial arms containing a total of 34,881 patients. Parameter estimates are generally consistent. It is found that a significant proportion of the variation in rates of tumour shrinkage and regrowth are explained by differing cancer and treatment: cancer type accounts for 66% of the variation in shrinkage rate and 71% of the variation in reported regrowth rates. Mean average parameter values by cancer and treatment are also reported.

CONCLUSIONS

Mathematical modelling of longitudinal data is most often reported on a per clinical trial basis. However, the results reported here suggest that a more integrative approach would benefit the development of new treatments as well as the further optimisation of those currently used.

Application of Machine Learning for Tumor Growth Inhibition - Overall Survival Modeling Platform

Machine learning (ML) was used to leverage tumor growth inhibition (TGI) metrics to characterize the relationship with overall survival (OS) as a novel approach and to compare with traditional TGI‐OS modeling methods. Historical dataset from a phase III non‐small cell lung cancer study (OAK, atezolizumab vs. docetaxel, N = 668) was used. ML methods support the validity of TGI metrics in predicting OS. With lasso, the best model with TGI metrics outperforms the best model without TGI metrics. Boosting was the best linear ML method for this dataset with reduced estimation bias and lowest Brier score, suggesting better prediction accuracy. Random forest did not outperform linear ML methods despite hyperparameter optimization. Kernel machine was marginally the best nonlinear ML method for this dataset and uncovered nonlinear and interaction effects. Nonlinear ML may improve prediction by capturing nonlinear effects and covariate interactions, but its predictive performance and value need further evaluation with larger datasets.

Clone Wars: Quantitatively Understanding Cancer Drug Resistance

A key aim of early clinical development for new cancer treatments is to detect the potential for efficacy early and to identify a safe therapeutic dose to take forward to phase II. Because of this need, researchers have sought to build mathematical models linking initial radiologic tumor response, often assessed after 6 to 8 weeks of treatment, with overall survival. However, there has been mixed success of this approach in the literature. We argue that evolutionary selection pressure should be considered to interpret these early efficacy signals and so optimize cancer therapy.

Characterizing Exposure-Response Relationship for Therapeutic Monoclonal Antibodies in Immuno-Oncology and Beyond: Challenges, Perspectives, and Prospects

Recent data from immuno-oncology clinical studies have shown the exposure-response (E-R) relationship for therapeutic monoclonal antibodies (mAbs) was often confounded by various factors due to the complex interplay of patient characteristics, disease, drug exposure, clearance, and treatment response and presented challenges in characterization and interpretation of E-R analysis. To tackle the challenges, exposure relationships for therapeutic mAbs in immuno-oncology and oncology are reviewed, and a general framework for an integrative understanding of E-R relationship is proposed. In this framework, baseline factors, drug exposure, and treatment response are envisioned to form an interconnected triangle, driving the E-R relationship and underlying three components that compose the apparent relationship: exposure-driven E-R, baseline-driven E-R, and response-driven E-R. Various strategies in data analysis and study design to decouple those components and mitigate the confounding effect are reviewed for their merits and limitations, and a potential roadmap for selection of these strategies is proposed. Specifically, exposure metrics based on a single-dose pharmacokinetic model can be used to mitigate response-driven E-R, while multivariable analysis and/or case control analysis of data obtained from multiple dose levels in a randomized study may be used to account for the baseline-driven E-R. In this context, the importance of collecting data from multiple dose levels, the role of prognostic factors and predictive factors, the potential utility of clearance at baseline and its change over time, and future directions are discussed.

A New Method to Model and Predict Progression Free Survival Based on Tumor Growth Dynamics

Progression-free survival (PFS) has been increasingly used as a primary endpoint for early clinical development. The aim of the present work was to develop a model where target lesion dynamics and risk for nontarget progression are jointly modeled for predicting PFS. The model was developed based on a pooled platinum-resistant ovarian cancer dataset comprising four different treatments and a wide range of dose levels. The target lesion progression was derived from tumor growth dynamics based on the Response Evaluation Criteria in Solid Tumors (RECIST) criteria. The nontarget progression hazard was correlated to the first derivative of target lesion tumor size with respect to time. The PFS time was determined by the first occurring event, target lesion progression, or nontarget progression. The final joint model not only captured target lesion tumor growth dynamics but also predicted PFS well. A similar approach can potentially be used to predict PFS in future oncology studies.

Progress and Opportunities to Advance Clinical Cancer Therapeutics Using Tumor Dynamic Models

There is a need for new approaches and endpoints in oncology drug development, particularly with the advent of immunotherapies and the multiple drug combinations under investigation. Tumor dynamics modeling, a key component to oncology "model-informed drug development," has shown a growing number of applications and a broader adoption by drug developers and regulatory agencies in the past years to support drug development and approval in a variety of ways. Tumor dynamics modeling is also being investigated in personalized cancer therapy approaches. These models and applications are reviewed and discussed, as well as the limitations and issues open for further investigations. A close collaboration between stakeholders like clinical investigators, statisticians, and pharmacometricians is warranted to advance clinical cancer therapeutics.

Alternative dosing regimens for atezolizumab: an example of model-informed drug development in the postmarketing setting

PURPOSE

To determine the exposure-response (ER) relationships between atezolizumab exposure and efficacy or safety in patients with advanced non-small cell lung cancer (NSCLC) or urothelial carcinoma (UC) and to identify alternative dosing regimens.

METHODS

ER analyses were conducted using pooled NSCLC and UC data from phase 1 and 3 studies (PCD4989g, OAK, IMvigor211; ClinicalTrials.gov IDs, NCT01375842, NCT02008227, and NCT02302807, respectively). Objective response rate, overall survival, and adverse events were evaluated vs pharmacokinetic (PK) metrics. Population PK-simulated exposures for regimens of 840 mg every 2 weeks (q2w) and 1680 mg every 4 weeks (q4w) were compared with the approved regimen of 1200 mg every 3 weeks (q3w) and the maximum assessed dose (MAD; 20 mg/kg q3w). Phase 3 IMpassion130 (NCT02425891) data were used to validate the PK simulations for 840 mg q2w. Observed safety data were evaluated by exposure and body weight subgroups.

RESULTS

No significant ER relationships were observed for safety or efficacy. Predicted exposures for 840 mg q2w and 1680 mg q4w were comparable to 1200 mg q3w and the MAD and consistent with observed PK data from IMpassion130. Observed safety was similar between patients with a C_max above and below the predicted C_max for 1680 mg q4w and between patients in the lowest and upper 3 body weight quartiles.

CONCLUSION

Atezolizumab regimens of 840 mg q2w and 1680 mg q4w are expected to have comparable efficacy and safety as the approved regimen of 1200 mg q3w, supporting their interchangeable use and offering patients greater flexibility.

Optimizing drug development in oncology by clinical trial simulation: Why and how?

In therapeutic research, the safety and efficacy of pharmaceutical products are necessarily tested on humans via clinical trials after an extensive and expensive preclinical development period. Methodologies such as computer modeling and clinical trial simulation (CTS) might represent a valuable option to reduce animal and human assays. The relevance of these methods is well recognized in pharmacokinetics and pharmacodynamics from the preclinical phase to postmarketing. However, they are barely used and are poorly regarded for drug approval, despite Food and Drug Administration and European Medicines Agency recommendations. The generalization of CTS could be greatly facilitated by the availability of software for modeling biological systems, by clinical trial studies and hospital databases. Data sharing and data merging raise legal, policy and technical issues that will need to be addressed. Development of future molecules will have to use CTS for faster development and thus enable better patient management. Drug activity modeling coupled with disease modeling, optimal use of medical data and increased computing speed should allow this leap forward. The realization of CTS requires not only bioinformatics tools to allow interconnection and global integration of all clinical data but also a universal legal framework to protect the privacy of every patient. While recognizing that CTS can never replace 'real-life' trials, they should be implemented in future drug development schemes to provide quantitative support for decision-making. This in silico medicine opens the way to the P4 medicine: predictive, preventive, personalized and participatory.

Association Between Tumor Size Kinetics and Survival in Patients With Urothelial Carcinoma Treated With Atezolizumab: Implication for Patient Follow-Up

We characterized the association between tumor size kinetics and survival in patients with advanced urothelial carcinoma treated with atezolizumab (anti-programmed death-ligand 1, Tecentriq) using a joint model. The model, developed on data from 309 patients of a phase II clinical trial, identified the time-to-tumor growth and the instantaneous changes in tumor size as the best on-treatment predictors of survival. On the validation dataset containing data from 457 patients from a phase III study, the model predicted individual survival probability using 3-month or 6-month tumor size follow-up data with an area under the receptor-occupancy curve between 0.75 and 0.84, as compared with values comprised between 0.62 and 0.75 when the model included only information available at treatment initiation. Including tumor size kinetics in a relevant statistical framework improves the prediction of survival probability during immunotherapy treatment and may be useful to identify most-at-risk patients in "real-time."

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.

Seamless Designs: Current Practice and Considerations for Early-Phase Drug Development in Oncology

Traditionally, drug development has evaluated dose, safety, activity, and comparative benefit in a sequence of phases using trial designs and endpoints specifically devised for each phase. Innovations in drug development seek to consolidate the phases and rapidly expand accrual with "seamless" trial designs. Although consolidation and rapid accrual may yield efficiencies, widespread use of seamless first-in-human (FiH) trials without careful consideration of objectives, statistical analysis plans, or trial oversight raises concerns. A working group formed by the National Cancer Institute convened to consider and discuss opportunities and challenges for such trials as well as encourage responsible use of these designs. We reviewed all abstracts presented at American Society of Clinical Oncology annual meetings from 2010 to 2017 for FiH trials enrolling at least 100 patients. We identified 1786 early-phase trials enrolling 57 559 adult patients. Fifty-one of the trials (2.9%) investigated 50 investigational new drugs, were seamless, and accounted for 14.6% of the total patients. The seamless trials included a median of 3 (range = 1-13) expansion cohorts. The overall risk of clinically significant treatment-related adverse events (grade 3-4) was 49.1% (range = 0.0-100%), and seven studies reported at least one toxic death. Rapid expansion of FiH trials may lead to earlier drug approval and corresponding widespread patient access to active therapeutics. Nevertheless, seamless designs must adhere to established ethical, scientific, and statistical standards. Protocols should include prospectively planned analyses of efficacy in disease- or biomarker-defined cohorts of sufficient rigor to support accelerated approval.

Comparison of tumor size assessments in tumor growth inhibition-overall survival models with second-line colorectal cancer data from the VELOUR study

PURPOSE

To compare lesion-level and volumetric measures of tumor burden with sum of the longest dimensions (SLD) of target lesions on overall survival (OS) predictions using time-to-growth (TTG) as predictor.

METHODS

Tumor burden and OS data from a phase 3 randomized study of second-line FOLFIRI ± aflibercept in metastatic colorectal cancer were available for 918 patients out of 1216 treated (75%). A TGI model that estimates TTG was fit to the longitudinal tumor size data (nonlinear mixed effect modeling) to estimate TTG with: SLD, sum of the measured lesion volumes (SV), individual lesion diameters (ILD), or individual lesion volumes (ILV). A parametric OS model was built with TTG estimates and assessed for prediction of the hazard ratio (HR) for survival.

RESULTS

Individual lesions had consistent dynamics within individuals. Between-lesion variability in rate constants was lower (typically < 27% CV) than inter-patient variability (typically > 50% CV). Estimates of TTG were consistent (around 12 weeks) across tumor size assessments. TTG was highly significant in a log-logistic parametric model of OS (median over 12 months). When individual lesions were considered, TTG of the fastest progressing lesions best predicted OS. TTG obtained from the lesion-level analyses were slightly better predictors of OS than estimates from the sums, with ILV marginally better than ILD. All models predicted VELOUR HR equally well and all predicted study success.

CONCLUSION

This analysis revealed consistent TGI profiles across all tumor size assessments considered. TTG predicted VELOUR HR when based on any of the tumor size measures.

Rendering the 3 + 3 Design to Rest: More Efficient Approaches to Oncology Dose-Finding Trials in the Era of Targeted Therapy

Selection of the maximum tolerated dose (MTD) as the recommended dose for registration trials based on a dose-escalation trial using variations of an MTD/3 + 3 design often occurs in the development of oncology products. The MTD/3 + 3 approach is not optimal and may result in recommended doses that are unacceptably toxic for many patients and in dose reduction/interruptions that might have an impact on effectiveness. Instead of the MTD/3 + 3 approach, the authors recommend an integrated approach. In this approach, typically an adaptive/Bayesian model provides a general framework to incorporate and make decisions for dose escalation based on nonclinical data, such as animal efficacy and toxicity data; clinical data, including pharmacokinetics/pharmacodynamics data; and dose/exposure-response data for efficacy and safety. To improve dose-ranging trials, model-based estimation, rather than hypothesis testing, should be used to maximize and integrate the information gathered across trials and doses. This approach may improve identification of optimal recommended doses, which can then be confirmed in registration trials. Clin Cancer Res; 22(11); 2623-9. ©2016 AACR SEE ALL ARTICLES IN THIS CCR FOCUS SECTION, "NEW APPROACHES FOR OPTIMIZING DOSING OF ANTICANCER AGENTS".

Population Modeling of Tumor Kinetics and Overall Survival to Identify Prognostic and Predictive Biomarkers of Efficacy for Durvalumab in Patients With Urothelial Carcinoma

Durvalumab is an anti‐PD‐L1 monoclonal antibody approved for patients with locally advanced or metastatic urothelial carcinoma (UC) that has progressed after platinum‐containing chemotherapy. A population tumor kinetic model, coupled with dropout and survival models, was developed to describe longitudinal tumor size data and predict overall survival in UC patients treated with durvalumab (NCT01693562) and to identify prognostic and predictive biomarkers of clinical outcomes. Model‐based covariate analysis identified liver metastasis as the most influential factor for tumor growth and immune‐cell PD‐L1 expression and baseline tumor burden as predictive factors for tumor killing. Tumor or immune‐cell PD‐L1 expression, liver metastasis, baseline hemoglobin, and albumin levels were identified as significant covariates for overall survival. These model simulations provided further insights into the impact of PD‐L1 cutoff values on treatment outcomes. The modeling framework can be a useful tool to guide patient selection and enrichment strategies for immunotherapies across various cancer indications.

On the relationship between tumour growth rate and survival in non-small cell lung cancer

A recurrent question within oncology drug development is predicting phase III outcome for a new treatment using early clinical data. One approach to tackle this problem has been to derive metrics from mathematical models that describe tumour size dynamics termed re-growth rate and time to tumour re-growth. They have shown to be strong predictors of overall survival in numerous studies but there is debate about how these metrics are derived and if they are more predictive than empirical end-points. This work explores the issues raised in using model-derived metric as predictors for survival analyses. Re-growth rate and time to tumour re-growth were calculated for three large clinical studies by forward and reverse alignment. The latter involves re-aligning patients to their time of progression. Hence, it accounts for the time taken to estimate re-growth rate and time to tumour re-growth but also assesses if these predictors correlate to survival from the time of progression. I found that neither re-growth rate nor time to tumour re-growth correlated to survival using reverse alignment. This suggests that the dynamics of tumours up until disease progression has no relationship to survival post progression. For prediction of a phase III trial I found the metrics performed no better than empirical end-points. These results highlight that care must be taken when relating dynamics of tumour imaging to survival and that bench-marking new approaches to existing ones is essential.

Model-Based Meta-Analysis for Multiple Myeloma: A Quantitative Drug-Independent Framework for Efficient Decisions in Oncology Drug Development

The failure rate for phase III trials in oncology is high; quantitative predictive approaches are needed. We developed a model‐based meta‐analysis (MBMA) framework to predict progression‐free survival (PFS) from overall response rates (ORR) in relapsed/refractory multiple myeloma (RRMM), using data from seven phase III trials. A Bayesian analysis was used to predict the probability of technical success (PTS) for achieving desired phase III PFS targets based on phase II ORR data. The model demonstrated a strongly correlated (R² = 0.84) linear relationship between ORR and median PFS. As a representative application of the framework, MBMA predicted that an ORR of ∼66% would be needed in a phase II study of 50 patients to achieve a target median PFS of 13.5 months in a phase III study. This model can be used to help estimate PTS to achieve gold‐standard targets in a target product profile, thereby enabling objectively informed decision‐making.

Modeling and simulation of maintenance treatment in first-line non-small cell lung cancer with external validation

Background

Maintenance treatment (MTx) in responders following first-line treatment has been investigated and practiced for many cancers. Modeling and simulation may support interpretation of interim data and development decisions. We aimed to develop a modeling framework to simulate overall survival (OS) for MTx in NSCLC using tumor growth inhibition (TGI) data.

Methods

TGI metrics were estimated using longitudinal tumor size data from two Phase III first-line NSCLC studies evaluating bevacizumab and erlotinib as MTx in 1632 patients. Baseline prognostic factors and TGI metric estimates were assessed in multivariate parametric models to predict OS. The OS model was externally validated by simulating a third independent NSCLC study (n = 253) based on interim TGI data (up to progression-free survival database lock). The third study evaluated pemetrexed + bevacizumab vs. bevacizumab alone as MTx.

Results

Time-to-tumor-growth (TTG) was the best TGI metric to predict OS. TTG, baseline tumor size, ECOG score, Asian ethnicity, age, and gender were significant covariates in the final OS model. The OS model was qualified by simulating OS distributions and hazard ratios (HR) in the two studies used for model-building. Simulations of the third independent study based on interim TGI data showed that pemetrexed + bevacizumab MTx was unlikely to significantly prolong OS vs. bevacizumab alone given the current sample size (predicted HR: 0.81; 95 % prediction interval: 0.59–1.09). Predicted median OS was 17.3 months and 14.7 months in both arms, respectively. These simulations are consistent with the results of the final OS analysis published 2 years later (observed HR: 0.87; 95 % confidence interval: 0.63–1.21). Final observed median OS was 17.1 months and 13.2 months in both arms, respectively, consistent with our predictions.

Conclusions

A robust TGI-OS model was developed for MTx in NSCLC. TTG captures treatment effect. The model successfully predicted the OS outcomes of an independent study based on interim TGI data and thus may facilitate trial simulation and interpretation of interim data. The model was built based on erlotinib data and externally validated using pemetrexed data, suggesting that TGI-OS models may be treatment-independent. The results supported the use of longitudinal tumor size and TTG as endpoints in early clinical oncology studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2455-2) contains supplementary material, which is available to authorized users.

Simulations to Predict Clinical Trial Outcome of Bevacizumab Plus Chemotherapy vs. Chemotherapy Alone in Patients With First-Line Gastric Cancer and Elevated Plasma VEGF-A

To simulate clinical trials to assess overall survival (OS) benefit of bevacizumab in combination with chemotherapy in selected patients with gastric cancer (GC), a modeling framework linking OS with tumor growth inhibition (TGI) metrics and baseline patient characteristics was developed. Various TGI metrics were estimated using TGI models and data from two phase III studies comparing bevacizumab plus chemotherapy vs. chemotherapy as first‐line therapy in 976 GC patients. Time‐to‐tumor‐growth (TTG) was the best TGI metric to predict OS. TTG, Eastern Cooperative Oncology Group (ECOG) score, albumin level, and Asian ethnicity were significant covariates in the final OS model. The model correctly predicted a decreased hazard ratio favorable to bevacizumab in patients with high baseline plasma VEGF‐A above the median of 113.4 ng/L. Based on trial simulations, in trials enrolling patients with elevated baseline plasma VEGF‐A (500 patients per arm), the expected hazard ratio was 0.82 (95% prediction interval: 0.70–0.95), independent of ethnicity.

Predicting survival of pancreatic cancer patients treated with gemcitabine using longitudinal tumour size data

Purpose

Measures derived from longitudinal tumour size data have been increasingly utilised to predict survival of patients with solid tumours. The aim of this study was to examine the prognostic value of such measures for patients with metastatic pancreatic cancer undergoing gemcitabine therapy.

Methods

The control data from two Phase III studies were retrospectively used to develop (271 patients) and validate (398 patients) survival models. Firstly, 31 baseline variables were screened from the training set using penalised Cox regression. Secondly, tumour shrinkage metrics were interpolated for each patient by hierarchical modelling of the tumour size time-series. Subsequently, survival models were built by applying two approaches: the first aimed at incorporating model-derived tumour size metrics in a parametric model, and the second simply aimed at identifying empirical factors using Cox regression. Finally, the performance of the models in predicting patient survival was evaluated on the validation set.

Results

Depending on the modelling approach applied, albumin, body surface area, neutrophil, baseline tumour size and tumour shrinkage measures were identified as potential prognostic factors. The distributional assumption on survival times appeared to affect the identification of risk factors but not the ability to describe the training data. The two survival modelling approaches performed similarly in predicting the validation data.

Conclusions

A parametric model that incorporates model-derived tumour shrinkage metrics in addition to other baseline variables could predict reasonably well survival of patients with metastatic pancreatic cancer. However, the predictive performance was not significantly better than a simple Cox model that incorporates only baseline characteristics.

Electronic supplementary material

The online version of this article (doi:10.1007/s00280-016-2994-x) contains supplementary material, which is available to authorized users.

Increasing the Time Interval between PCV Chemotherapy Cycles as a Strategy to Improve Duration of Response in Low-Grade Gliomas: Results from a Model-Based Clinical Trial Simulation

Background. We previously developed a mathematical model capturing tumor size dynamics of adult low-grade gliomas (LGGs) before and after treatment either with PCV (Procarbazine, CCNU, and Vincristine) chemotherapy alone or with radiotherapy (RT) alone. Objective. The aim of the present study was to present how the model could be used as a simulation tool to suggest more effective therapeutic strategies in LGGs. Simulations were performed to identify schedule modifications that might improve PCV chemotherapy efficacy. Methods. Virtual populations of LGG patients were generated on the basis of previously evaluated parameter distributions. Monte Carlo simulations were performed to compare treatment efficacy across in silico clinical trials. Results. Simulations predicted that RT plus PCV would be more effective in terms of duration of response than RT alone. Additional simulations suggested that, in patients treated with PCV chemotherapy, increasing the interval between treatment cycles up to 6 months from the standard 6 weeks can increase treatment efficacy. The predicted median duration of response was 4.3 years in LGGs treated with PCV cycles given every 6 months versus 3.1 years in patients treated with the classical regimen. Conclusion. The present study suggests that, in LGGs, mathematical modeling could facilitate clinical research by helping to identify, in silico, potentially more effective therapeutic strategies.

Prediction of Response to Temozolomide in Low-Grade Glioma Patients Based on Tumor Size Dynamics and Genetic Characteristics

Both molecular profiling of tumors and longitudinal tumor size data modeling are relevant strategies to predict cancer patients' response to treatment. Herein we propose a model of tumor growth inhibition integrating a tumor's genetic characteristics (p53 mutation and 1p/19q codeletion) that successfully describes the time course of tumor size in patients with low-grade gliomas treated with first-line temozolomide chemotherapy. The model captures potential tumor progression under chemotherapy by accounting for the emergence of tissue resistance to treatment following prolonged exposure to temozolomide. Using information on individual tumors' genetic characteristics, in addition to early tumor size measurements, the model was able to predict the duration and magnitude of response, especially in those patients in whom repeated assessment of tumor response was obtained during the first 3 months of treatment. Combining longitudinal tumor size quantitative modeling with a tumor''s genetic characterization appears as a promising strategy to personalize treatments in patients with low-grade gliomas.

A Tumor Growth Inhibition Model Based on M-Protein Levels in Subjects With Relapsed/Refractory Multiple Myeloma Following Single-Agent Carfilzomib Use

Change in tumor size estimated using longitudinal tumor growth inhibition (TGI) modeling is an early predictive biomarker of clinical outcomes for multiple cancer types. We present the application of TGI modeling for subjects with multiple myeloma (MM). Longitudinal time course changes in M‐protein data from relapsed and/or refractory MM subjects who received single‐agent carfilzomib in phase II studies (n = 456) were fit to a TGI model. The tumor growth rate estimate was similar to that of other anti‐myeloma agents, indicating that the model is robust and treatment‐independent. An overall survival model was subsequently developed, which showed that early change in tumor size (ECTS) at week 4, Eastern Cooperative Oncology Group performance status (ECOG PS), hemoglobin, sex, percent bone marrow cell involvement, and number of prior regimens were significant independent predictors for overall survival (P < 0.001). ECTS based on M‐protein modeling could be an early biomarker for survival in MM following exposure to single‐agent carfilzomib.

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.

Early tumour shrinkage as a prognostic factor and surrogate end-point in colorectal cancer: a systematic review and pooled-analysis

PURPOSE

Early tumour shrinkage (ETS), defined as a reduction of at least 20% in tumour size at first reassessment, has been recently investigated retrospectively in first-line trials of metastatic colorectal cancer (CRC), and appears to be associated with better outcomes. We have performed a systematic review and meta-analysis of published trials to evaluate the prognostic value of ETS in CRC in terms of overall survival (OS) and progression-free survival (PFS).

MATERIAL AND METHODS

An electronic search of the PubMed, SCOPUS, EMBASE, the Web of Science, and the Cochrane Central Register of Controlled Trial databases identified trials that compared outcomes of patients with or without ETS during first-line chemotherapy for metastatic CRC. The OS, reported as a hazard ratio (HR) with a 95% confidence interval (CI), was the primary outcome measure; the correlation coefficient (R) between ETS with median OS was also estimated.

RESULTS

Twenty-one trials from 10 publications were analysed. Overall, patients with ETS were associated with a better OS (HR, 0.58; 95% CI, 0.53 to 0.64; P<0.00001) and PFS (HR, 0.57; 95% CI, 0.47-0.69; P<0.00001) compared with patients who were early non-responders. However, ETS was poorly correlated with OS in terms of surrogacy (R=0.37; 95% CI - 0.31-0.78; P=0.28).

CONCLUSIONS

ETS is a good prognostic factor but an inappropriate surrogate for predicting outcome in CRC patients. These findings support ETS as prognostic tool in ascertaining earlier non-responders; however, its role as a surrogate end-point deserves further evaluation.

Model-Based Design of a Decision Tree for Treating HER2+ Cancers Based on Genetic and Protein Biomarkers

Human cancers are incredibly diverse with regard to molecular aberrations, dependence on oncogenic signaling pathways, and responses to pharmacological intervention. We wished to assess how cellular dependence on the canonical PI3K vs. MAPK pathways within HER2+ cancers affects responses to combinations of targeted therapies, and biomarkers predictive of their activity. Through an integrative analysis of mechanistic model simulations and in vitro cell line profiling, we designed a six-arm decision tree to stratify treatment of HER2+ cancers using combinations of targeted agents. Activating mutations in the PI3K and MAPK pathways (PIK3CA and KRAS), and expression of the HER3 ligand heregulin determined sensitivity to combinations of inhibitors against HER2 (lapatinib), HER3 (MM-111), AKT (MK-2206), and MEK (GSK-1120212; trametinib), in addition to the standard of care trastuzumab (Herceptin). The strategy used to identify effective combinations and predictive biomarkers in HER2-expressing tumors may be more broadly extendable to other human cancers.

Developing Exposure/Response Models for Anticancer Drug Treatment: Special Considerations

Anticancer agents often have a narrow therapeutic index (TI), requiring precise dosing to ensure sufficient exposure for clinical activity while minimizing toxicity. These agents frequently have complex pharmacology, and combination therapy may cause schedule-specific effects and interactions. We review anticancer drug development, showing how integration of modeling and simulation throughout development can inform anticancer dose selection, potentially improving the late-phase success rate. This article has a companion article in Clinical Pharmacology & Therapeutics with practical examples.