Assessing the impact of organ-specific lesion dynamics on survival in patients with recurrent urothelial carcinoma treated with atezolizumab or chemotherapy

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.

Tumor growth and overall survival modeling to support decision making in phase Ib/II trials: A comparison of the joint and two-stage approaches

Model-based tumor growth inhibition (TGI) metrics are increasingly used to predict overall survival (OS) data in Phase III immunotherapy clinical trials. However, there is still a lack of understanding regarding the differences between two-stage or joint modeling methods to leverage Phase I/II trial data and help early decision-making. A recent study showed that TGI metrics such as the tumor growth rate constant KG may have good operating characteristics as early endpoints. This previous study used a two-stage approach that is easy to implement and intuitive but prone to bias as it does not account for the relationship between the longitudinal and time-to-event processes. A relevant alternative is to use a joint modeling approach. In the present article, we evaluated the operating characteristics of TGI metrics using joint modeling, assuming an OS model previously developed using historical data. To that end, we used TGI and OS data from IMpower150-a study investigating atezolizumab in over 750 patients suffering from non-small cell lung cancer-to mimic randomized Phase Ib/II trials varying in terms of number of patients included (40 to 15 patients per arm) and follow-up duration (24 to 6 weeks after the last patient included). In this context, joint modeling did not outperform the two-stage approach and provided similar operating characteristics in all the investigated scenarios. Our results suggest that KG geometric mean ratio could be used to support early decision-making provided that 30 or more patients per arm are included and followed for at least 12 weeks.

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.

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.

Modeling tumor size dynamics based on real-world electronic health records and image data in advanced melanoma patients receiving immunotherapy

The development of immune checkpoint inhibitors (ICIs) has revolutionized cancer therapy but only a fraction of patients benefits from this therapy. Model-informed drug development can be used to assess prognostic and predictive clinical factors or biomarkers associated with treatment response. Most pharmacometric models have thus far been developed using data from randomized clinical trials, and further studies are needed to translate their findings into the real-world setting. We developed a tumor growth inhibition model based on real-world clinical and imaging data in a population of 91 advanced melanoma patients receiving ICIs (i.e., ipilimumab, nivolumab, and pembrolizumab). Drug effect was modeled as an ON/OFF treatment effect, with a tumor killing rate constant identical for the three drugs. Significant and clinically relevant covariate effects of albumin, neutrophil to lymphocyte ratio, and Eastern Cooperative Oncology Group (ECOG) performance status were identified on the baseline tumor volume parameter, as well as NRAS mutation on tumor growth rate constant using standard pharmacometric approaches. In a population subgroup (n = 38), we had the opportunity to conduct an exploratory analysis of image-based covariates (i.e., radiomics features), by combining machine learning and conventional pharmacometric covariate selection approaches. Overall, we demonstrated an innovative pipeline for longitudinal analyses of clinical and imaging RWD with a high-dimensional covariate selection method that enabled the identification of factors associated with tumor dynamics. This study also provides a proof of concept for using radiomics features as model covariates.

Tumor Dynamic Model-Based Decision Support for Phase Ib/II Combination Studies: A Retrospective Assessment Based on Resampling of the Phase III Study IMpower150

PURPOSE

Model-based tumor growth inhibition (TGI) metrics are increasingly incorporated into go/no-go decisions in early clinical studies. To apply this methodology to new investigational combinations requires independent evaluation of TGI metrics in recently completed Phase III trials of effective immunotherapy.

PATIENTS AND METHODS

Data were extracted from IMpower150, a positive, randomized, Phase III study of first-line therapy in 1,202 patients with non-small cell lung cancer. We resampled baseline characteristics and longitudinal sum of longest diameters of tumor lesions of patients from both arms, atezolizumab+ bevacizumab+chemotherapy (ABCP) versus BCP, to mimic Phase Ib/II studies of 15 to 40 patients/arm with 6 to 24 weeks follow-up. TGI metrics were estimated using a bi-exponential TGI model. Effect sizes were calculated as TGI metrics geometric mean ratio (GMR), objective response rate (ORR) difference (d), and progression-free survival (PFS), hazard ratio (HR) between arms. Correct and incorrect go decisions were evaluated as the probability to achieve desired effect sizes in ABCP versus BCP and BCP versus BCP, respectively, across 500 replicated subsamples for each design.

RESULTS

For 40 patients/24 weeks follow-up, correct go decisions based on probability tumor growth rate (KG) GMR <0.90, dORR >0.10, and PFS HR <0.70 were 83%, 69%, and 58% with incorrect go decision rates of 4%, 12%, and 11%, respectively. For other designs, the ranking did not change with TGI metrics consistently overperforming RECIST endpoints. The predicted overall survival (OS) HR was around 0.80 in most of the scenarios investigated.

CONCLUSIONS

Model-based estimate of KG GMR is an exploratory endpoint that informs early clinical decisions for combination studies.

Materials-driven approaches to understand extrinsic drug resistance in cancer

Metastatic cancer has a poor prognosis, because it is broadly disseminated and associated with both intrinsic and acquired drug resistance. Critical unmet needs in effectively killing drug resistant cancer cells include overcoming the drug desensitization characteristics of some metastatic cancers/lesions, and tailoring therapeutic regimens to both the tumor microenvironment and the genetic profiles of the resident cancer cells. Bioengineers and materials scientists are developing technologies to determine how metastatic sites exclude therapies, and how extracellular factors (including cells, proteins, metabolites, extracellular matrix, and abiotic factors) at metastatic sites significantly affect drug pharmacodynamics. Two looming challenges are determining which feature, or combination of features, from the tumor microenvironment drive drug resistance, and what the relative impact is of extracellular signals vs. intrinsic cell genetics in determining drug response. Sophisticated systems biology tools that can de-convolve a crowded network of signals and responses, as well as controllable microenvironments capable of providing discrete and tunable extracellular cues can help us begin to interrogate the high dimensional interactions governing drug resistance in patients.