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

Information on individual lesion dynamics and organ location are often ignored in pharmacometric modeling analyses of tumor response. Typically, the sum of their longest diameters is utilized. Herein, a tumor growth inhibition model was developed for describing the individual lesion time-course data from 183 patients with metastatic HER2-negative breast cancer receiving docetaxel. The interindividual variability (IIV), interlesion variability (ILV), and interorgan variability of parameters describing the lesion time-courses were evaluated. Additionally, a model describing the probability of new lesion appearance and a time-to-event model for overall survival (OS), were developed. Before treatment initiation, the lesions were largest in the soft tissues and smallest in the lungs, and associated with a significant IIV and ILV. The tumor growth rate was 2.6 times higher in the breasts and liver, compared with other metastatic sites. The docetaxel drug effect in the liver, breasts, and soft tissues was greater than or equal to 1.2 times higher compared with other organs. The time-course of the largest lesion, the presence of at least 3 liver lesions, and the time since study enrollment, increased the probability of new lesion appearance. New lesion appearance, along with the time to growth and time-course of the largest lesion at baseline, were identified as the best predictors of OS. This tumor modeling approach, incorporating individual lesion dynamics, provided a more complete understanding of heterogeneity in tumor growth and drug effect in different organs. Thus, there may be potential to tailor treatments based on lesion location, lesion size, and early lesion response to provide better clinical outcomes.

WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC?

The individual tumor lesion data, which comprise the sum of the longest diameters (SLDs), contain lesion dynamic data and organ location that are often ignored in the investigation of predictors of long‐term clinical end points, such as overall survival (OS). The lesion data may better describe disease progression and OS more so than a composite SLD.

WHAT QUESTION DID THIS STUDY ADDRESS?

The time‐course of individual metastatic lesions was characterized and quantified the interlesion and interorgan variability in patients with HER2‐negative breast cancer treated with docetaxel. Further, the ability of various lesion‐related metrics to predict appearance of new lesions, dropout from tumor measurements, and OS were evaluated.

WHAT DOES THIS STUDY ADD TO OUR KNOWLEDGE?

The developed tumor modeling approach, incorporating individual lesion dynamics, provided a more complete understanding of heterogeneity in tumor dynamics in different organs. The tumor growth rate and drug‐induced shrinkage were different between metastatic sites, whereas the rate of appearance of resistance toward docetaxel was the same in all organs. Metastasis of the liver was associated with a poor prognosis as it was associated with increased hazard of both the appearance of new lesions as well as death. Appearance of new lesions along with metrics driven from tumor dynamics were predictive of OS.

HOW MIGHT THIS CHANGE DRUG DISCOVERY, DEVELOPMENT, AND/OR THERAPEUTICS?

The developed modeling framework provided information about the treatment response/tumor growth in various metastatic sites. This modeling approach can be applied to other cancer types and therapies to provide a better understanding of the relationship among drug exposure, individual lesion‐organ response, and new lesion appearance/OS and thereby facilitate early clinical interventions to improve anticancer therapy.

INTRODUCTION

Pharmacometric models are increasingly applied to data collected from oncology trials toward a better understanding of the drug response in patients over time. The time‐course of the sum of the longest diameters (SLDs), or metrics thereof, are investigated as predictors of long‐term clinical end points, such as overall survival (OS). , , , The individual tumor lesion data, which comprise the SLD, contain lesion dynamic data, and organ location that are often ignored. These data may better describe disease progression and OS more so than a composite SLD measurement.

Tumors that originate from the same organ can behave and respond to treatment differently at different metastatic sites even though they are histologically similar, indicating that the growth and drug‐induced shrinkage of individual tumor lesions may be highly dependent on the microenvironment. , , , , It is noteworthy that the metastasis of vital organs is one of the main reasons for death in patients with cancer, and metastasis‐associated death was reported as high as 90% in some cancer types. , , Breast cancer, the most common cancer in women, has a 5‐year survival rate of 99% if the patient had cancer only in the breasts; the survival rate drops to 26% for patients with metastasis, and ~ 75% of deaths in breast cancer were associated with metastasis. Liver, lungs, and brain are the most common metastasis sites having a strong relation to cancer relapse and OS. ,

Traditionally the tumor (SLD) changes are usually assessed based on the response evaluation criteria in solid tumors (RECIST). The initial version of RECIST (version 1.0) was based on the sum of 10 lesions and a later version (version 1.1) is based on the sum of 5 lesions. When SLD is evaluated, the progression of an individual lesion may be hidden by shrinkage in other lesions, or progression may not be apparent if the increasing lesion has a relatively small absolute size compared with the other target lesions. Consequently, SLD‐based response categorization ignores individual lesions’ and the metastatic organs’ contribution to clinical end points. ,

Pharmacometric models have an advantage in that different levels of variability can be quantified; however, there are only a few examples where pharmacometric modeling has been applied to characterize lesion dynamics. , , , , For example, analyses by Mercier et al. Schindler et al., , and Claret et al. found significant interlesion variability (ILV) in model parameters, suggesting tumor heterogeneity. Additionally, in a nonparametric analysis using a novel methodology (i.e., “classification clustering of individual lesions”), similarities between lesion dynamics within the same anatomic site classification were observed; however, many patients demonstrated very diverse lesion dynamics when comparing across different tumor site classifications. A model‐based approach, incorporating the understanding of the time‐course of the disease in distinct metastatic sites, and the contribution of lesions to overall response, would help in predicting the expected responses at different timepoints.

In the present study, a pharmacometric modeling approach was developed to describe the time‐course of individual lesions. Detailed survival analysis was performed where, in addition to baseline covariates, tumor dynamics in an organ or of a lesion, and appearance of new lesions were investigated as predictors of the OS. The approach was developed based on data from patients with HER2‐negative metastatic breast cancer receiving docetaxel therapy. In addition to survival analysis, we explored the ability of various lesion‐related metrics to predict new lesions and dropout from tumor measurements.

METHODS

Study population

The tumor data, including metastatic sites, target lesion diameters (LDs), and appearance of new lesions were available from the docetaxel arm of the phase III AVADO trial where the efficacy and safety of combining bevacizumab with docetaxel were investigated in patients with HER2‐negative metastatic breast cancer (ClinicalTrials.gov Identifier: NCT00333775). From the docetaxel arm (n = 241), patients who had no measurable target lesions at baseline (n = 34) or had received bevacizumab before disease progression (n = 24) were excluded. Thus, the current study population comprised of 183 subjects with a median age of 55 years (range 29–83 years). Patients received docetaxel 100 mg/m2 (i.e., 164 mg for a typical woman) infused over 1 h on day 1 of each 3‐week cycle. The tumor size response was evaluated according to RECIST version 1.0 (i.e., up to 10 lesions/patient). The lesions were in the population distributed in seven different metastatic sites, namely breasts, liver, lungs, lymph nodes, mediastinum, soft tissues, and other sites (including adrenal, skin, pelvis, and kidneys). Lesion size and locations were evaluated from computed tomography scans and the median tumor size follow‐up was 32 weeks (range 6–145 weeks). The treatment continued until disease progression, as assessed by RECIST version 1.0, or until intolerable toxicity occurred. The main characteristics of the study population are summarized in Table 1. The AVADO trial was conducted according to the Declaration of Helsinki, the Good Clinical Practice guidelines of the International Conference on Harmonization, and the laws and regulations of the countries involved. The protocol was approved by local ethics committees and written informed consent was obtained from all patients before the screening.

TABLE 1

Summary of clinical characteristics of the study population

Characteristics	Median	Range
Total number of patients, n	183	‐
Age, years	55	29–83
Tumor baseline size, mm	34	10–140
Sum of longest diameters at baseline	69	10–308
Tumor follow‐up, a weeks	32	5–145
Number of organs with metastasis, per patient	2	1–6
Number of lesions, per patient	3	1–10
New lesion appearance (yes), n	121	66%
Time of new lesion appearance, weeks	34.4	5.43–111
Death events, n	93	51%
Survival follow‐up, weeks	108	12–160
Time to death, weeks	50	12–145

Lesion size and locations were evaluated from computed tomography scans every 9 weeks during the first 36 weeks, and thereafter every 12 weeks.

Lesion model

A tumor growth inhibition (TGI) model, without modification, was applied to describe the change in individual lesions over time. In this model, the tumor growth rate is described by an exponential growth function with a first‐order rate constant (k GROW). The tumor size shrinkage rate during treatment is explained by drug exposure, the drug‐specific cell kill rate constant (k SHR), and the emergence of resistance to the treatment (λ; Equation 1). As docetaxel concentrations were not available, a population pharmacokinetic‐pharmacodynamic (PK‐PD) modeling approach (Equation 2) was used along with the TGI model.

where LD ,1 is the diameter of lesion 1 in organ “a,” is the growth rate constant for organ a, is the docetaxel‐related cell kill rate constant for organ a of the i subject. is the elimination rate constant in the PK‐PD model and Docetaxel(t) describes the docetaxel exposure driven by dose and elimination rate.

In addition to interindividual (IIV) variability, ILV, , and interorgan variabilities (IORGVs) in lesion baseline, k GROW, λ, and k SHR parameters, were initially considered but it led to model instability and estimation errors. Therefore, ILV in the lesion at baseline and IORGV in k GROW, λ, and k SHR parameters were explored. The typical baseline value for the lesions was evaluated to be organ‐specific. For the parameter for the i th subject (), and an organ a, the value for lesion, j = 1, 2…5 can be written as:

where, is the typical parameter value for organ a, is the random IIV that is common to all organs and lesions within individual i, and K is the random ILV specific to the lesion, organ, and individual.

The IORGV in k GROW, λ, and k SHR was evaluated as fixed effects. The IORVG in k GROW parameter for different organs a–e for ith subject () can be written as:

where, are the typical parameter value for k GROW for organs a, b, c, d, and e; and is the random IIV. The parametrization for λ and k SHR were similar to k GROW parametrization (Equation 4). Both η and κ are assumed to be normally distributed with mean 0 and variance of ω 2 and π 2, respectively. The ILV for a specific organ was assumed to have a common variance magnitude regardless of the lesion number (i.e., ). The organ‐specific k GROW, λ, and k SHR were allowed to be shared between organs (e.g., ) if they do not differ statistically.

New lesion appearance and dropout from tumor size measurements

Patient dropout was defined as the stopping of docetaxel treatment due to disease progression or “other” event, or at the scheduled end of the trial. The probability of a new lesion or dropout from tumor size measurements can be high for those individuals who have a poor treatment response and/or develop resistance to docetaxel. In simulations, to account for the dropping out from tumor measurement follow‐up, a logistic regression model was used (Equation 5). The probability of identifying new lesion(s) at scheduled measurement time points was described using a logistic regression model. A sequential approach similar to population PK parameters and data (PPP&D) was used in lesion‐new lesion appearance‐dropout modeling.where,

The different baseline and time‐varying factors that were tested as predictors on the probability of developing new lesions are summarized in Table S1. The predictors evaluated for dropout from tumor size measurements were a 20% increase (in combination with at least a 5 mm absolute increase) in SLD from the tumor nadir, time in the study, baseline SLD, time‐course of SLD, and the appearance of new lesions. ,

Overall survival model

A parametric time to event (TTE) model with a Weibull function was used for describing the baseline hazard (h 0(t)) of OS events. The effect of potential predictors on h 0(t) was explored using a sequential analysis method similar to the individual PK parameter approach. The long computational times of the lesion model did not allow for using other estimation approaches like the PPP&D or simultaneous approaches when linking the lesion model predictions to the OS model. Individual empirical Bayes estimates (EBEs) from the final lesion model were consequently used to explore the lesion time‐courses, or model‐derived metrics, as predictors of OS. Predictors were investigated alone in an initial univariate analysis and then evaluated in combination (multivariate). Four separate analyses were performed with different types of predictors; (i) baseline characteristics only, (ii) appearance of new lesions, (iii) baseline LD and metrics derived from LD, and (iv) baseline SLD and metrics derived from SLD (sum of either up to 10 lesions [RECIST version 1.0] or up to 5 lesions [RECIST version 1.1]). New lesion was evaluated as a time‐varying categorical covariate (NewLesion = 0 until the new lesion appeared, thereafter NewLesion = 1).

The tumor dataset based on RECIST version 1.1 was created by selecting one (largest) lesion per organ because none of the patients had more than 5 metastatic sites at baseline. For individuals who had a total number of lesions less than five and had more than one lesion/organ, a second lesion/organ was selected from (in the order of priority) liver, lymph nodes, lungs, or other available sites to have a total number of lesions equal to five.

Model development and evaluation

Population models were developed using the nonlinear mixed‐effect modeling software (NONMEM version 7.4). The first‐order conditional estimation method with interaction was used in lesion model development and the Laplace method was used in logistic regression and TTE models. The EBEs from the final lesion model were obtained using the Laplace method with interaction because the estimation method used in regression/TTE models was the Laplace method. Logistic regression models for new lesion appearance and dropout from tumor size measurements were developed in combination and estimated simultaneously. Model development was assisted by Pirana (version 2.9.9) for run management, the Perl‐speaks‐NONMEM (PsN) toolkit for handling NONMEM run commands, R (version 3.6), and Xpose (version 4.1) for model diagnostics and graphical analysis.

The objective function value (OFV; −2∙log‐likelihood) and graphical diagnostics were used in the evaluation of model performance. An OFV decrease of 6.64 (p < 0.01) was considered as significant for the addition of one parameter (1 degree of freedom) in the lesion model and an OFV decrease of 3.84 (p < 0.05) was applied in TTE and logistic regression models. The actual decrease in OFV (dOFV) required for p less than 0.05 and 1 degree of freedom in the TTE model was determined to 3.36 using the randomization test (randtest) in the PsN toolkit. While investigating shared k GROW, λ, or k SHR parameter between two metastatic organs (i.e., reduction in one parameter [1 less degree of freedom], or an increase in OFV of 10.83 [p < 0.001]) was considered as significant. Parameter uncertainties were derived using Sampling Importance Resampling or the NONMEM‐provided R matrix. Visual predictive checks (VPCs) for the lesion model and Kaplan–Meier VPCs for the logistic and TTE models were used for evaluating their predictive performance. For VPCs of lesion, new lesion, and dropout models, the final joint model (lesion‐new lesion appearance‐dropout model) was used in the simulations to account for the response‐related dropout observed in the study. Patients were recruited into the study at different dates and they censored from OS follow‐up at the study closure date or at a random time. A competing TTE model was developed in OS simulations to describe the censoring from OS follow‐up.

RESULTS

Data

The lesion data consisted of 2045 observations from 586 lesions distributed among 7 metastatic sites in 183 patients with metastatic HER2‐negative breast cancer (Figure S1). At baseline, a median of three (range 1–10) lesions per patient, from two (range 1–4) metastatic sites were present. More than one metastatic site was present in 44% (n = 81) of the patients. Liver and lymph nodes were the most common metastatic sites, where 49% (n = 90) of the patients in the population had at least one lesion in the liver and 46% (n = 84) had at least one lesion in the lymph nodes.

A TGI model described the individual lesion data and the schematic representation of the structural lesion model are given in Figure S2. The baseline lesion size was characterized by an organ‐specific value along with IIV and ILV (Equation 4). At study enrollment, the lesions were typically largest in soft tissues (44 mm), and smallest in the lungs (19 mm). The most common metastatic organs— liver (24 mm) and lymph nodes (20 mm)—had about half of the size of the typical soft tissue lesions. The ILV was not significantly different between different sites (41% coefficient of variation [CV]) except for the ‘other sites” (67% CV). Notably, ILV was higher than IIV (29% CV) for the baseline lesion diameter. The growth rate (k GROW) and docetaxel drug effect (k SHR) were significantly different between organs, whereas λ was not affected by metastatic site (i.e., the rate of appearance of resistance was the same in all metastatic sites). The fastest growing lesions were from the liver and breasts, the k GROW, Liver was not significantly different from k GROW, Breast, and the estimate was 0.00917 week−1 (a doubling time () of ~ 1.5 years). The slowest growth rate was estimated in k GROW, Lymph nodes, k GROW, Mediastinum, and k GROW, Other, and was 0.00341 week−1 (doubling time of ~ 4 years). The drug effect in the liver, breasts, and soft tissues did not differ statistically and were estimated to be 1.7 times higher than the lungs, mediastinum, and other sites, and 1.23 times higher than the lymph nodes. There was significant IIV in growth rate (135% CV), λ (63% CV), and drug effect (43% CV).

VPCs, illustrated in Figure 1, indicate that the final lesion model can adequately describe the data from liver and lymph nodes (most common organs where lesions were observed) and for SLD (=lesion1+lesion2+…+lesion10). VPCs of all organs and a few individual fits are provided in Figure S3 and Figure S4, respectively. The final lesion model parameters along with their uncertainties are presented in Table 2.

FIGURE 1

Visual predictive checks (VPCs) of the final lesion model. Sum of their longest diameter (SLD; left) was defined as lesion1+lesion2+…+lesion10, most common metastatic sites (i.e., liver [middle] and lymph nodes [right]). The red solid line represents the median of the observed tumor measurements (circles) and the blue solid lines represent the 5th and 95th percentiles of the observed tumor measurements. The inner shaded region represents the 95% confidence interval of the model simulated median (red dashed line). The outer shaded regions represent the 95% confidence intervals of the model simulated 5th and 95th percentiles and black dashed lines are the simulated median of the corresponding percentiles. Vertical lines indicate binning intervals for VPCs

TABLE 2

Parameter estimates and their uncertainty in the final lesion model

Parameter	Description (unit)	Typical value (RSE b )	IIV, CV% (RSE b )	ILV, CV% (RSE b )
k GROW, Lung	Tumor growth rate constant (week−1)	0.00453 (37)	135 (15)	‐
k GROW, Liver		0.00917 (32)		‐
k GROW, Breast
k GROW, Lymph nodes		0.00341 (37)		‐
k GROW, Mediastinum
k GROW, Other
k GROW, Soft tissue		0.00676 (48)		‐
k DRUG, Lung	Tumor kill rate constant (week−1)	0.000742 (27)	43 (19)
k DRUG, Mediastinum
k DRUG, Other
k DRUG, Liver		0.00123 (27)
k DRUG, Breast
k DRUG, Soft tissue
k DRUG, Lymph nodes		0.00101 (28)
Λ	Resistance parameter (week−1)	0.126 (26)	63 (16)	–
LD0, Lung	Baseline lesion diameter (mm)	19.1 (14)	31 (14)	40 (9)
LD0, Liver		24.2 (11)
LD0, Lymph nodes		20.4 (11)
LD0, Soft tissue		44.0 (33)
LD0, Breast		41.9 (24)
LD0, Mediastinum		23.0 (25)
LD0, Other		28.8 (31)		67 (16)
k KPD	Docetaxel elimination rate constant in PK‐PD model (week−1)	0.576 (34)	62 (16)	–
RUV a	Residual unexplained variability (%)	20.5 (4)	–	–

Abbreviations: CV%, percent coefficient of variation; IIV, interindividual variability; ILV, interlesion variability; PD, pharmacodynamic; PK, pharmacokinetic; RSE, relative standard error.

Additive residual error model on log transformed data.

Obtained from Sampling Importance Resampling.

New lesion appearance and dropout from tumor size follow‐up

The probability of new lesion appearance was dependent on the time‐course of the largest lesion at baseline, three or more liver lesions at enrollment, and time from the start of treatment (on logarithmic scale). The most significant predictor in the univariate analysis, the time‐course of the largest lesion (dOFV = 13), was associated with a hazard ratio (HR) of 1.71 for every 10% increase in lesion diameter. The patients who had 3 or more target lesions in the liver at baseline had a 62% higher risk (HR = 1.62) of developing a new lesion compared with patients who had less than 3 lesions in the liver. Every week in the study added about 4% increase in risk (HR = 1.04) of developing a new lesion.

As expected, both 20% increase in SLDs (HR = 8.92) and the appearance of new lesions (HR = 5.47) were significantly related to dropout. No other predictors retained significance level after inclusion of 20% increase in SLDs and the appearance of new lesions as predictors of dropout from tumor follow‐up. The final parameter estimates and their uncertainty are given in Table 3 and predictive performances as illustrated in VPCs are shown in Figure 2.

TABLE 3

Parameter estimates and their uncertainty in the final new lesion appearance and dropout from tumor size model

Parameter	Description	Typical value (RSE b )
New lesion appearance model
InterceptNew_lesion	Parameter relating to baseline probability of developing new lesion	−7.07 (8.0)
β_Largest_lesion(t)	Coefficient of the effect of the time course of largest lesion	0.488 (23)
β_>3Liver lesion	Coefficient of the effect of ≥3 liver lesions (y/n)	0.481 (47)
β_Time	Coefficient of the effect of time since study start	0.502 (24)
Dropout from tumor size measurement model
InterceptDropout	Parameter relating to baseline probability of dropping out	−4.59 (3.0)
β_New_lesion	Coefficient of the effect of appearance of new lesions on dropout	2.88 (7.0)
β_Disease_progression	Coefficient of the effect of disease progression a on dropout	1.99 (14)

Abbreviations: RSE, relative standard error; SLD, sum of their longest diameter.

Disease progression defined as 20% increase in SLD from tumor nadir and an absolute increase of 5 mm.

Obtained from NONMEM R‐matrix.

FIGURE 2

Kaplan–Meier visual predictive checks for the final new lesion appearance model (left) and dropout from tumor size (right). The observed Kaplan–Meier curve (black line) is compared to the 95% confidence interval (shaded area) derived from model simulations (100 samples)

Final overall survival model

A parametric model with Weibull function adequately described the increase in the hazard of death over time. The final OS model showed that the appearance of new lesions, increase in tumor time‐course, and shorter time‐to‐tumor growth (TTG) of SLD (RECIST version 1.0) were associated with a higher hazard of death. In univariate analysis, the most significant predictor of OS was the appearance of a new lesion. However, the hazard of death was not significantly related to the number of new lesions (i.e., regardless of how many new lesions that were present at the time point of new lesion identification), the patient would have a worse prognosis (HR = 2.83) after the identification of new lesions. The OS model based on time‐course and TTG of the largest lesion at baseline (OS model iii) had similar level of statistical significance (dOFV =28, p < 0.001) as that of the OS model based on time‐course and TTG of SLD‐RECIST version 1.0 (dOFV = 34, p < 0.001; OS model iv). Addition of time‐course or TTG of the largest lesion at baseline or RECIST version 1.1 based SLD, on top of SLD‐RECIST version 1.0 based predictors did not show any improvements. Furthermore, the presence of liver lesions, dynamics of liver lesions, more than two target lesions at baseline, and baseline SLD were significant in univariate analysis, but did not retain level of significance in multivariate analysis. Detailed results from OS model i–iv analysis are given in Supplementary OS results. The parameter estimates and VPCs of the final OS model are given in Table S2 and Figure 3, respectively.

FIGURE 3

Kaplan–Meier visual predictive checks for the final overall survival model. The observed Kaplan–Meier curve (black line) is compared to the 95% confidence interval (shaded area) derived from model simulations (200 samples)

DISCUSSION

In this study, we successfully applied a modeling approach to characterize the longitudinal lesion time‐course and quantified different types of variability in the lesion dynamics in patients with HER2‐negative breast cancer receiving docetaxel treatment. ILV was found to be significant to include in the baseline lesion diameter, but not for the tumor growth rate constant or the drug‐induced shrinkage rate. The results show a significant difference in drug‐induced shrinkage rate and growth rate between different metastatic sites within an individual. Besides, baseline and tumor‐related metrics were explored for their association with new lesion appearance, dropout from tumor size measurements, and OS.

Interestingly, the presence of liver lesions increased the hazard of both the appearance of new lesions as well as death. To our knowledge, there are few other pharmacometric analyses on the appearance of new lesions despite that it is one important criterion to define progressive disease and indicator of poor prognosis. , Liver metastasis was a significant predictor of new lesions, and the appearance of a new lesion was identified as a predictor of OS in a pharmacometric analysis by Zecchin et al. The liver was the most common metastatic site in the current study data and thus the observed shrinkage and regrowth as monitored in SLDs might have been primarily driven by liver lesion dynamics. A high incidence of liver metastasis in breast cancer and the analysis using traditional statistics demonstrated liver metastasis as an independent predictor of worse OS in different subtypes of breast cancer. , The current modeling framework could be applied to other cancers as well where metastasis site data are collected and to identify the influential metastasis on OS.

Our study results suggest that OS can be predicted by the appearance of new lesions, model‐predicted lesion time‐course, along with TTG of the SLD. The tumor time‐course , and TTG , , , , of the SLD or lesion level data, as well as combination of time‐course and TTG, have earlier been suggested as predictors of OS. However, in the present study, we have compared predictors of survival based on both SLD and lesion level information. A study by Claret et al. determined that the TTG of the fastest progressing lesion was the best predictor of OS and predictions based on the lesion model performed slightly better than SLD. In contrast, the current study results have shown TTG of SLD had performed marginally better than the single, largest lesion‐based metric. In both univariate and multivariate analysis (OS model iii), the dynamics of the largest lesion at baseline has valuable information on assessing the risk of death of an individual patient. Moreover, the same metrics (i.e., TS(t) along with TTG) were included in the final OS models based on lesion (OS model iii) and SLD (OS model iv) as predictors, but the OFV drop associated with the addition of predictors based on SLD was slightly better than those of lesion based. In addition, we could conclude that following 5 lesions, as in RECIST version 1.1, would in this study population perform just as good as the information on data on 10 lesions (RECIST version 1.0) in assessing the risk of death. A couple of earlier studies , have compared the tumor response categorization in RECIST version 1.0 versus RECIST version 1.1 and found an excellent agreement in assessing tumor response between versions. Our study results support these findings and add valuable information that regardless of the version of RECIST applied, the tumor time‐course and model predicted TTG predict survival. Instead of the sum of lesions, even the time‐course of the largest lesion provided an equal level of significance (p < 0.001) to predict OS in the current analysis. The selection of the lesions in a patient with multiple lesions could be subjective, and collecting a maximum of 2 lesions per organ, as in RECIST version 1.1, is likely enough to account for the variability in the response. Interestingly, baseline covariates were either not significant or did not retain level of significance (SLD0 and LES10) once the time‐course (which contain baseline tumor information) was included as a predictor of OS.

A model based on SLD combined with risk for new lesion would be enough for investigating predictors of the survival. Thus, the modeling of individual lesions is not intended to replace the traditional SLD based analysis—instead, it provides more detailed information about the treatment response/tumor growth in various metastatic sites. Individual lesion modeling provides a quantitative measure of various variability for an anticancer drug and indication and may find further application in combination chemotherapies, when one drug may directly target a specific receptor expression organ, and the other drug has nonspecific targeting. Furthermore, this modeling framework allows evaluation of individual lesion dynamics and organ‐specific dynamics as a predictor of OS. An extension to this work could be to explore if information of individual lesion dynamics would improve the prediction of the individual hazard of death compared with SLD‐based metrics. The number of differential equations (one differential equation for every lesion‐organ combination) and the estimation of different levels of variability lead to long (~ 9 h at 15 nodes) computational run times, hence a parsimonious model was applied in the current study. The TGI model, applied in the current study, has been used for different anticancer drug classes in various indications. , , , , , Alternative approaches that could be considered include tumor models requiring multiple compartments, such as biexponential, sensitive‐resistant, , , , or proliferative‐quiescent. , The loss of further tumor measurements at the time of progression restricted the evaluation of more complex growth curves and the developed model may not necessarily be valid for extrapolation beyond disease progression. However, the individual fits and VPCs based on TGI model did a reasonable job of describing the profiles at both lesion and SLD level data.

In conclusion, this model framework successfully described the lesion time‐course and quantification of IIV, ILV, and IORGV in different model parameters. ILV was here significant in the organ‐specific baseline lesion diameters, whereas significant IORGV was estimated in growth and shrinkage parameters. The probability of developing a new lesion and the hazard of death were increased by the presence of liver metastasis, indicating that the liver is the major metastatic site in the treatment prognosis for this group of patients. The appearance of new lesions along with tumor time‐course and TTG of SLD were identified as the best predictors of the clinical end point OS. This modeling approach can be applied to other cancer types and therapies to provide a better understanding of the relationship among drug exposure, individual lesion‐organ response, and new lesion appearance/OS and thereby facilitate early clinical interventions to improve anticancer therapy.

CONFLICT OF INTEREST

B.C.B. and A.L.Q. are employees of Genentech Inc. L.E.F. has acted as a paid consultant to Genentech Inc. All other authors declared no competing interests for this work.

AUTHOR CONTRIBUTIONS

S.M.K., S.S.L., B.C.B., A.L.Q., and L.E.F. wrote the manuscript. S.M.K., S.S.L., B.C.B., A.L.Q., and L.E.F. designed the research. S.M.K. and S.S.L. performed the research. S.M.K., S.S.L., B.C.B., A.L.Q., and L.E.F. analyzed the data.

DISCLAIMER

As Deputy Editor‐in‐Chief of CPT: Pharmacometrics & Systems Pharmacology, Lena E. Friberg was not involved in the review or decision process for this paper.

Supplementary Material

Click here for additional data file.

Supplementary Material

Click here for additional data file.

Supplementary Material

Click here for additional data file.

Supplementary Material

Click here for additional data file.