Tumor growth inhibition modeling of individual lesion dynamics and interorgan variability in HER2-negative breast cancer patients treated with docetaxel

Predicting tumor dynamics in treated patients from patient-derived-xenograft mouse models: a translational model-based approach

This study presents a translational modeling framework designed to predict tumor size dynamics in cancer patients undergoing anticancer treatment, using data from patient-derived xenograft (PDX) mice. In the first step, a population tumor growth inhibition (TGI) model to estimate the distribution of exponential tumor growth rates and anticancer drug potency in PDX mice was built. Then, model parameters were allometrically scaled from mice to humans to inform a TGI model predicting tumor size dynamics in cancer patients. Longitudinal tumor dynamics predicted by the PDX-informed TGI model were expressed in terms of tumor progression events to allow validation against literature time-to-progression (TTP) data. The proposed approach was tested on two case studies: gemcitabine treatment for pancreatic cancer and sorafenib treatment for hepatocellular cancer. The framework successfully predicted median tumor size dynamics, closely aligned with clinical TTP curves for gemcitabine-pancreatic cancer case study. While predictions for extreme tumor size percentiles highlighted potential avenues for refinement, such as incorporating resistance mechanisms, the overall accuracy underscored the goodness of the approach. For the sorafenib-hepatocellular cancer case study, the framework provided plausible tumor size predictions, with TTP curves closely aligned with clinical observations, despite the limited availability of clinical data prevented a full validation. Overall, the translational modeling framework showed potential for predicting tumor dynamics in cancer patients, with results suggesting its applicability as a valid tool to support early decision-making in oncology.

Integrated modeling of biomarkers, survival and safety in clinical oncology drug development

Model-based approaches, including population pharmacokinetic-pharmacodynamic modeling, have become an essential component in the clinical phases of oncology drug development. Over the past two decades, models have evolved to describe the temporal dynamics of biomarkers and tumor size, treatment-related adverse events, and their links to survival. Integrated models, defined here as models that incorporate at least two pharmacodynamic/ outcome variables, are applied to answer drug development questions through simulations, e.g., to support the exploration of alternative dosing strategies and study designs in subgroups of patients or other tumor indications. It is expected that these pharmacometric approaches will be expanded as regulatory authorities place further emphasis on early and individualized dosage optimization and inclusive patient-focused development strategies. This review provides an overview of integrated models in the literature, examples of the considerations that need to be made when applying these advanced pharmacometric approaches, and an outlook on the expected further expansion of model-informed drug development of anticancer drugs.

Predicting Tumor Volume Doubling Time and Progression-Free Survival in Untreated Patients from Patient-Derived-Xenograft (PDX) Models: A Translational Model-Based Approach

Tumor volume doubling time (TVDT) has been shown to be a potential surrogate marker of biological tumor activity. However, its availability in clinics is strongly limited due to ethical and practical reasons, as its assessment requires at least two subsequent tumor volume measurements in untreated patients. Here, a translational modeling framework to predict TVDT distributions in untreated cancer patient populations from tumor growth data in patient-derived xenograft (PDX) mice is proposed. Eleven solid cancer types were considered. For each of them, a set of tumor growth studies in PDX mice was selected and analyzed through a mathematical model to characterize the distribution of the exponential tumor growth rate in mice. Then, assuming an exponential growth of the tumor mass in humans, the growth rates were scaled from PDX mice to humans through an allometric scaling approach and used to predict TVDTs in untreated patients. A very good agreement was found between model predicted and clinically observed TVDTs, with 91% of the predicted TVDT medians fell within 1.5-fold of observations. Further, exploiting the intrinsic relationship between tumor growth dynamics and progression free survival (PFS), the exponential growth rates in humans were used to generate the expected PFS curves in absence of anticancer treatment. Predicted curves were extremely close to published PFS data from studies involving patient cohorts treated with supportive care or low effective therapies. The proposed approach shows promise as a potential tool to increase knowledge about TVDT in humans without the need of directly measuring tumor dimensions in untreated patients, and to predict PFS curves in untreated patients, that could fill the absence of placebo-controlled arms against which to compare treaded arms during clinical trials. However, further validation and refinement are needed to fully assess its effectiveness in this regard.

Organ-specific tumor dynamics predict survival of patients with metastatic colorectal cancer

BACKGROUND

We aim to compare the prognostic value of organ-specific dynamics with the sum of the longest diameter (SLD) dynamics in patients with metastatic colorectal cancer (mCRC).

METHODS

All datasets are accessible in Project Data Sphere, an open-access platform. The tumor growth inhibition models developed based on organ-level SLD and SLD were used to estimate the organ-specific tumor growth rates (KGs) and SLD KG. The early tumor shrinkage (ETS) from baseline to the first measurement after treatment was also evaluated. The relationship between organ-specific dynamics, SLD dynamics, and survival outcomes (overall survival, OS; progression-free survival, PFS) was quantified using Kaplan-Meier analysis and Cox regression.

RESULTS

This study included 3687 patients from 6 phase III mCRC trials. The liver emerged as the most frequent metastatic site (2901, 78.7 %), with variable KGs across different organs in individual patients (liver 0.0243 > lung 0.0202 > lymph node 0.0127 > other 0.0118 [week-1]). Notably, the dynamics for different organs did not equally contribute to predicting survival outcomes. In liver metastasis cases, liver KG proved to be a superior prognostic indicator for OS and surpasses the predictive performance of SLD, (C-index, liver KG 0.610 vs SLD KG 0.606). A similar result can be found for PFS. Moreover, liver ETS also outperforms SLD ETS in predicting survival. Cox regression analysis confirmed liver KG is the most significant variable in survival prediction.

CONCLUSIONS

In mCRC patients with liver metastasis, liver dynamics is the primary prognostic indicator for both PFS and OS. In future drug development for mCRC, greater emphasis should be directed towards understanding the dynamics of liver metastasis development.

Nonlinear multilevel joint model for individual lesion kinetics and survival to characterize intra-individual heterogeneity in patients with advanced cancer

In advanced cancer patients, tumor burden is calculated using the sum of the longest diameters (SLD) of the target lesions, a measure that lumps all lesions together and ignores intra-patient heterogeneity. Here, we used a rich dataset of 342 metastatic bladder cancer patients treated with a novel immunotherapy agent to develop a Bayesian multilevel joint model that can quantify heterogeneity in lesion dynamics and measure their impact on survival. Using a nonlinear model of tumor growth inhibition, we estimated that dynamics differed greatly among lesions, and inter-lesion variability accounted for 21% and 28% of the total variance in tumor shrinkage and treatment effect duration, respectively. Next, we investigated the impact of individual lesion dynamics on survival. Lesions located in the liver and in the bladder had twice as much impact on the instantaneous risk of death compared to those located in the lung or the lymph nodes. Finally, we evaluated the utility of individual lesion follow-up for dynamic predictions. Consistent with results at the population level, the individual lesion model outperformed a model relying only on SLD, especially at early landmark times and in patients with liver or bladder target lesions. Our results show that an individual lesion model can characterize the heterogeneity in tumor dynamics and its impact on survival in advanced cancer patients.

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.

Dissecting sources of variability in patient response to targeted therapy: anti-HER2 therapies as a case study

BACKGROUND AND PURPOSE

Despite their use to treat cancers with specific genetic aberrations, targeted therapies elicit heterogeneous responses. Sources of variability are critical to targeted therapy drug development, yet there exists no method to discern their relative contribution to response heterogeneity.

EXPERIMENTAL APPROACH

We use HER2-amplified breast cancer and two agents, neratinib and lapatinib, to develop a platform for dissecting sources of variability in patient response. The platform comprises four components: pharmacokinetics, tumor burden and growth kinetics, clonal composition, and sensitivity to treatment. Pharmacokinetics are simulated using population models to capture variable systemic exposure. Tumor burden and growth kinetics are derived from clinical data comprising over 800,000 women. The fraction of sensitive and resistant tumor cells is informed by HER2 immunohistochemistry. Growth rate-corrected drug potency is used to predict response. We integrate these factors and simulate clinical outcomes for virtual patients. The relative contribution of these factors to response heterogeneity are compared.

KEY RESULTS

The platform was verified with clinical data, including response rate and progression-free survival (PFS). For both neratinib and lapatinib, the growth rate of resistant clones influenced PFS to a higher degree than systemic drug exposure. Variability in exposure at labeled doses did not significantly influence response. Sensitivity to drug strongly influenced responses to neratinib. Variability in patient HER2 immunohistochemistry scores influenced responses to lapatinib. Exploratory twice daily dosing improved PFS for neratinib but not lapatinib.

CONCLUSION AND IMPLICATIONS

The platform can dissect sources of variability in response to target therapy, which may facilitate decision-making during drug development.

Multistate Pharmacometric Model to Define the Impact of Second-Line Immunotherapies on the Survival Outcome of the IMpower131 Study

Overall survival is defined as the time since randomization into the clinical trial to event of death or censor (end of trial or follow-up), and is considered to be the most reliable cancer end point. However, the introduction of second-line treatment after disease progression could influence survival and be considered a confounding factor. The aim of the current study was to set up a multistate model framework, using data from the IMpower131 study, to investigate the influence of second-line immunotherapies on overall survival analysis. The model adequately described the transitions between different states in patients with advanced squamous non-small cell lung cancer treated with or without atezolizumab plus nab-paclitaxel and carboplatin, and characterized the survival data. High PD-L1 expression at baseline was associated with a decreased hazard of progression, while the presence of liver metastasis at baseline was indicative of a high risk of disease progression after initial response. The hazard of death after progression was lower for participants who had longer treatment response, i.e., longer time to progression. The simulations based on the final multistate model showed that the addition of atezolizumab to the nab-paclitaxel and carboplatin regimen had significant improvement in the patients' survival (hazard ratio = 0.75, 95% prediction interval: 0.61-0.90 favoring the atezolizumab + nab-paclitaxel and carboplatin arm). The developed modeling approach can be applied to other cancer types and therapies to provide a better understanding of efficacy of drug and characterizing different states, and investigate the benefit of primary therapy in survival while accounting for the switch to alternative treatment in the case of disease progression.

Longitudinal and time-to-event modeling for prognostic implications of radical surgery in retroperitoneal sarcoma

Retroperitoneal sarcoma (RPS) is a rare malignancy which can be difficult to manage due to the variety of clinical behaviors. In this study, we aimed to develop a parametric modeling framework to quantify the relationship between postoperative dynamics of several biomarkers and overall/progression-free survival of RPS. One hundred seventy-four patients with RPS who received surgical resection with curative intent at the Peking University Cancer Hospital Sarcoma Center were retrospectively included. Potential prognostic factors were preliminarily identified. Longitudinal analyses of body mass index (BMI), serum total protein (TP), and white blood cells (WBCs) were performed using nonlinear mixed effects models. The impacts of time-varying and time-invariant predictors on survival were investigated by parametric time-to-event (TTE) models. The majority of patients experienced decline in BMI, recovery of TP, as well as transient elevation in WBC counts after surgery, which significantly correlated with survival. An indirect-response model incorporating surgery effect described the fluctuation in percentage BMI. The recovery of TP was captured by a modified Gompertz model, and a semimechanistic model was selected for WBCs. TTE models estimated that the daily cumulative average of predicted BMI and WBC, the seventh-day TP, as well as certain baseline variables, were significant predictors of survival. Model-based simulations were performed to examine the clinical significance of prognostic factors. The current work quantified the individual trajectories of prognostic biomarkers in response to surgery and predicted clinical outcomes, which would constitute an additional strategy for disease monitoring and intervention in postoperative RPS.

A longitudinal model for the Mayo Clinical Score and its sub-components in patients with ulcerative colitis

Clinical trials in patients with ulcerative colitis (UC) face the challenge of high and variable placebo response rates. The Mayo Clinical Score (MCS) is used widely as the primary endpoint in clinical trials to describe the clinical status of patients with UC. The MCS is comprised of four subscores, each scored 0, 1, 2 and 3: rectal bleeding (RB), stool frequency (SF), physician's global assessment (PGA), and endoscopy (ENDO) subscore. Excluding the PGA subscore gives the modified MCS. Quantitative insight on the placebo response, and its impact on the components of the MCS over time, can better inform clinical trial design and interpretation. Longitudinal modeling of the MCS, and the modified MCS, can be challenging due to complex clinical trial design, population heterogeneity, and limited assessments for the ENDO subscore. The current study pooled patient-level placebo/standard of care (SoC) arm data from five clinical trials in the TransCelerate database to develop a longitudinal placebo response model that describes the MCS over time in patients with UC. MCS subscores were modeled using proportional odds models, and the removal of patients from the placebo/SoC arm, or "dropout", was modeled using logistic regression models. The subscore and dropout models were linked to allow for the prediction of the MCS and the modified MCS. Stepwise covariate modeling identified prior exposure to TNF-α antagonists as a statistically significant predictor on the RB + SF subscore. Patients with prior exposure to TNF-α antagonists had higher post-baseline RB + SF subscores than naive patients.

Semi-Mechanistic Model for the Antitumor Response of a Combination Cocktail of Immuno-Modulators in Non-Inflamed (Cold) Tumors

Simple Summary

The clinical efficacy of immunotherapies when treating cold tumors is still low, and different treatment combinations are needed when dealing with this challenging scenario. In this work, a middle-out strategy was followed to develop a model describing the antitumor efficacy of different immune-modulator combinations, including an antigen, a toll-like receptor-3 agonist, and an immune checkpoint inhibitor in mice treated with non-inflamed tumor cells. Our results support that clinical response requires antigen-presenting cell activation and also relies on the amount of CD8 T cells and tumor resistance mechanisms present. This mathematical model is a very useful platform to evaluate different immuno-oncology combinations in both preclinical and clinical settings.

Abstract

Immune checkpoint inhibitors, administered as single agents, have demonstrated clinical efficacy. However, when treating cold tumors, different combination strategies are needed. This work aims to develop a semi-mechanistic model describing the antitumor efficacy of immunotherapy combinations in cold tumors. Tumor size of mice treated with TC-1/A9 non-inflamed tumors and the drug effects of an antigen, a toll-like receptor-3 agonist (PIC), and an immune checkpoint inhibitor (anti-programmed cell death 1 antibody) were modeled using Monolix and following a middle-out strategy. Tumor growth was best characterized by an exponential model with an estimated initial tumor size of 19.5 mm³ and a doubling time of 3.6 days. In the treatment groups, contrary to the lack of response observed in monotherapy, combinations including the antigen were able to induce an antitumor response. The final model successfully captured the 23% increase in the probability of cure from bi-therapy to triple-therapy. Moreover, our work supports that CD8⁺ T lymphocytes and resistance mechanisms are strongly related to the clinical outcome. The activation of antigen-presenting cells might be needed to achieve an antitumor response in reduced immunogenic tumors when combined with other immunotherapies. These models can be used as a platform to evaluate different immuno-oncology combinations in preclinical and clinical scenarios.