Skip to main content

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

Abstract

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.

Peer Review reports

Background

There is still an unmet medical need in the treatment of non-small cell lung cancer (NSCLC) in both the first-line and recurrent settings. Maintenance treatment has been investigated in patients with disease control (i.e. without progressive disease) during first-line therapy in a number of trials with the goal to prolong time to disease progression (progression-free survival, PFS), improve quality of life and ultimately prolong overall survival (OS) [14]. However, the risk-benefit ratio of maintenance therapy in NSCLC is still unclear, and several aspects of this strategy have raised considerable debate [2]. Therefore models that could predict the clinical outcomes of maintenance therapy may be of great importance to practitioners and drug developers.

Modeling and simulation may provide quantitative support for interpretation of interim data and development decisions in oncology [5, 6]. Tumor response of patients can be characterized using tumor growth inhibition (TGI) metrics, which are estimated based on modeling of longitudinal tumor size data. TGI metrics have been shown to predict treatment effect on OS in solid tumors and in multiple myeloma [5]. These TGI metrics include model-based estimates of change in tumor size from baseline at end of cycle 2 (e.g. week 6 or 8), tumor growth rate and time to tumor regrowth [5]. TGI metrics could be used as alternative endpoints [7] in early clinical studies to optimize drug dosing, support clinical trial design for investigational anti-cancer treatments [5, 6].

Although a few models linking OS with TGI metrics and prognostic factors have been published for NSCLC first-line [810] and second-line [8] therapies, there has been no investigation of TGI metrics and of their link to OS in the context of maintenance therapy to date. Furthermore, there is insufficient published external validation of such models. External validation is critical for assessing treatment independence of the models and favour acceptance [5]. Finally, the OS models are assumed to be disease-specific but treatment-independent. However, to date, there has been insufficient validation of the treatment-independence assumption.

Accumulation of valuable clinical data has made it possible to build and externally validate a TGI-OS model for maintenance therapy in NSCLC patients whose disease did not progress during first-line therapy. Erlotinib maintenance prolonged both PFS [11] and OS [12] in the SATURN trial. The addition of erlotinib to bevacizumab during maintenance therapy significantly prolonged PFS but not OS compared to the bevacizumab-only maintenance in the ATLAS trial [13]. The AVAPERL trial compared maintenance bevacizumab plus pemetrexed vs. bevacizumab alone and showed a significant prolongation of PFS [14] but not of OS [15] following bevacizumab plus pemetrexed compared to bevacizumab alone.

The objectives of this work were 1) to develop a model for OS after maintenance therapy in NSCLC based on erlotinib data from SATURN and ATLAS, 2) to prospectively predict the probability to success of AVAPERL study and perform an external validation by simulating the OS outcomes of AVAPERL study (pemetrexed data) based on interim tumor size data (up to PFS database lock).

Methods

Trials and data

Data were collected from all patients enrolled in three studies evaluating maintenance treatment after first-line NSCLC therapy. In all studies, patients whose disease did not progress after four cycles of first-line treatment were randomized to maintenance treatment. Details of the studies can be found in the respective papers, in the introduction section and in Table 1. The studies complied with the Declaration of Helsinki and Good Clinical Practice guidelines, and were approved at all investigating centers by local ethics committees. All patients provided written informed consent for participation and publication of the data [1115]. An ethics statement was not required for this analysis as they have been provided in each of the three individual studies [1115].

Table 1 Study summary

The SATURN trial compared maintenance erlotinib vs. placebo in patients whose disease did not progress after four cycles of platinum-based first-line chemotherapy [11, 12]. The ATLAS trial compared maintenance erlotinib plus bevacizumab vs. bevacizumab alone in patients whose disease did not progress after four cycles of platinum-doublet chemotherapy in combination with bevacizumab [13]. The AVAPERL trial compared maintenance bevacizumab plus pemetrexed vs. bevacizumab alone in patients whose disease did not progress after four cycles of first-line chemotherapy of cisplatin plus pemetrexed in combination with bevacizumab [14, 15].

The following baseline patient characteristics were tested as prognostic factors for OS based on SATURN and ATLAS data: age, gender, ethnicity, Eastern Cooperative Oncology Group (ECOG) score, smoking status, tumor size, and histology. In addition, study effects and response to first-line therapy were investigated. Interim AVAPERL data consisted in longitudinal tumor size collected by the time of PFS database lock (data cutoff: July 2011) and baseline patient characteristics only.

Tumor growth inhibition metrics

The full TGI profile was modeled using equations adapted from previously published simplified TGI models [16] (Fig. 1) that were fit to data from evaluable patients using a nonlinear mixed-effect modeling (population) approach (NONMEM, version 7, FOCE algorithm with interaction) [17]. To be evaluable in this analysis, patients had to have at least one tumor size measurement after randomization to maintenance treatment. Tumor size was assessed as the sum of longest diameters of target lesions by Response Evaluation Criteria In Solid Tumors (RECIST) [18, 19]. Shrinkage in model-parameter estimates was estimated as previously described [20]. Model fitting was assessed using standard goodness-of-fit plots.

Fig. 1
figure 1

Theoretical tumor size profile over first-line treatment run-in phase and maintenance treatment phase. t1: time of randomization to maintenance treatment; KL: growth rate (assumed to be same during the two treatment phases); TS, KDE and λ: tumor size, drug effect and drug effect decay over time, respectively, for the first-line (TS1, KDE1 and λ1) and maintenance (TS2, KDE2 and λ2) phase; RND: randomization

Two patient-level TGI metrics were calculated based on individual posthoc parameter estimates: the time to tumor regrowth (TTG) [16], and the week 8 ECTS (early change in tumor size) that represented early tumor shrinkage and was calculated as the ratio of model-predicted tumor size at week 8 to baseline estimated by the model. Equations are displayed in Fig. 1. Only the TGI metrics during the maintenance phase were of interest and were calculated.

Overall survival model development

Data from SATURN and ATLAS were used to build the OS model. The impact of individual factors on OS was assessed using Kaplan-Meier and Cox regression analyses using survfit and coxph functions, respectively in R (version 2.15.0) [21]. The baseline patient prognostic factors together with the TGI metrics were tested to explain variability in OS.

A parametric survival regression model (using the survreg function in R version 2.15.0) was developed that describes OS distribution. The probability density function that best describes the observed survival time was selected among normal, lognormal, Weibull, logistic, log-logistic, and exponential by using difference in Akaike information criterion (AIC) [22] of the alternative models.

A “full” model was built by including all significant covariates (baseline prognostic factors, TGI metrics) from the Cox univariate analysis with a significance level of p < 0.05 per the log-likelihood ratio test where the difference in −2*log-likelihood (score) between alternative models follows a χ2 distribution. The score indicates the level of significance for the association between this covariate and OS: the higher the score, the more significantly this covariate is associated with OS. Then a backward stepwise elimination was carried out. At each elimination step, one covariate was removed from the model. If the reduced model (without this removed covariate) became significantly worse (p < 0.01), the removed covariate stayed in the model. The relative influence of each remaining covariate on the model was re-evaluated by deleting it from the reduced model on an individual basis with a significance level of p < 0.01. The backward elimination resulted in the final model, in which all covariates were significant.

The model simulation performances were evaluated using a posterior predictive check. OS distributions and hazard ratios (HR) in SATURN and ATLAS were simulated 1000 times. Model parameters were sampled from the estimated mean values and uncertainty in parameter estimates for each of the simulated study replicate. Censoring was assumed to be 30 % as in the original data.

Simulations

OS of AVAPERL study were simulated based on TGI metrics estimated using interim tumor size data to predict the likelihood of a successful OS outcome for AVAPERL and further assess performance of the OS model (external validation). In order to calculate the prediction interval and make statistical inferences, the study was simulated multiple times (20,000) by sampling survival model parameters from their estimated uncertainty distribution. Patient survival times were drawn from the appropriate survival distribution defined by model parameters, baseline prognostic factors and TGI metric of AVAPERL patients. Censoring was simulated in sampling patient study duration, assumed to be independent of death. Patient survival times were censored assuming a uniform distribution of patient study duration from 50 to 140 weeks, which was consistent with the minimum and the maximum time period the patient stayed in the SATURN study without a death event. For each of the replicates, simulated data were analyzed by Kaplan-Meier estimation and Cox regression. Kaplan-Meier estimates of OS distributions and HR used to compare both arms were summarized by median and 95 % prediction interval (PI) across the replicates.

Results

Data

Patients with at least one post-randomization tumor size measurement were included in this analysis. Overall 1534 patients were evaluable to estimate TGI metrics used for building the OS model: 837 (94 %) out of 889 patients from SATURN, and 697 (94 %) out of 743 patients from ATLAS. Interim AVAPERL data used as the external validation dataset were collected by the time of PFS database lock (data cutoff: July 2011) and included 231 evaluable patients out of 245 (94 %) randomized to maintenance treatment.

Tumor size model

The simplified TGI model adequately described the observed tumor size data, as shown by goodness-of-fit plots and individual fits (Additional file 1: Figure S1 and Additional file 2: Figures S2). Parameters were adequately estimated with small standard errors and shrinkage (Table 2) except that inter-individual variability could not be estimated on λ1 due to the limited number of observations during first-line treatment phase. TGI metric estimates (TTG and week 8 ECTS) that were calculated from the TGI model parameters (Table 2) using equations displayed in Fig. 1 were highly variable: the range from 5th to 95th percentile was 0.721 (i.e. decrease in tumor size from baseline) to 1.24 (i.e. increase in tumor size from baseline) for week 8 ECTS, and −23 weeks to 70 weeks for TTG after randomization. TTG may take negative values when KL > KDE2, i.e. in patients with progression at the first assessment after randomization to maintenance phase (Additional file 2: Figure S2).

Table 2 Parameter estimates of the simplified TGI model

Overall survival model

In univariate Cox analysis (Table 3), TTG was the most significant covariate associated with OS (score 151.7) and much better than week 8 ECTS (score 45.1). The most significant baseline prognostic factors and patient characteristics were tumor size, gender, smoking status, Asian ethnicity and ECOG score (scores 8 to 50, p < 0.0001). Also OS tended to be longer in erlotinib treated patients and in ATLAS trial compared to SATURN (p < 0.01). OS distribution by quartiles of TTG is shown in Fig. 2.

Table 3 Screening of the potential covariates for overall survival using the Cox model
Fig. 2
figure 2

Overall survival by quartiles of TTG. Each group represents 25 % of the patients. TTG: time to tumor regrowth (week). OS: overall survival (week). Median estimates are reported in the insert

A lognormal distribution had the best likelihood to describe the OS distribution (lower AIC than other distributions). All covariates that were significant in the Cox univariate analysis were included in the “full” model, and underwent backward stepwise elimination. The final model included TTG and the following baseline prognostic factors: baseline tumor size, ECOG score (0 vs. >0), Asian ethnicity, age and gender. All parameters in the final OS model were estimated with good precision (Table 4). According to the model, good prognostic is predicted for patients with longer TTG (treatment effect), small baseline tumor size, age below 55 years, Asian ethnicity, ECOG score 0 and for female patients.

Table 4 Parameter estimates of the final overall survival model

The model was evaluated by simulating OS distributions in each of the study arms (Fig. 3) and the HR of treatment vs. control arm in SATURN and ATLAS (Fig. 4a and b). The observed HR (0.79 for SATURN and 0.93 for ATLAS) was within the 95 % PI by the model (0.74–0.97 for SATURN and 0.70–1.00 for ATLAS).

Fig. 3
figure 3

Posterior predictive check of the final OS model by studies: a) SATURN and b) ATLAS. Solid line: observed OS. Band: 95 % prediction interval of OS. OS: overall survival

Fig. 4
figure 4

Posterior predictive check of HR in SATURN (a) and ATLAS (b) for maintenance erlotinib compared to placebo and simulation of HR in AVAPERL (c) for maintenance pemetrexed vs. placebo and comparison to observed HR. CI: confidence interval. HR: hazard ratio. PI: prediction interval

Simulation

The final OS model was applied to prospectively predict the expected OS outcome of AVAPERL study (external validation). The goal was to predict the likelihood of a successful OS outcome using interim tumor size data collected by the time of PFS database lock (data cutoff: July 2011). This dataset was not used for model-building (Table 1). Median OS was not yet reached at the time of data cutoff, and the immature OS data that were observed by the time of data cutoff were not used. Patients in AVAPERL study had more favorable prognostic factors than those from SATURN and ATLAS with a smaller proportion of ECOG score >0 (52 % vs. 66–69 %) and smaller baseline tumor size (5.2 cm vs. >6 cm) (Table 1). Simulations indicated that pemetrexed plus bevacizumab as maintenance treatment in AVAPERL was unlikely to demonstrate a significant OS prolongation vs. bevacizumab alone. The expected HR was 0.81 with a 95 % PI of 0.59–1.09 (62 % of events), which contained 1 (Fig. 4c). Predicted median OS was 17.3 and 14.7 months in both arms, respectively. These prospective simulations were consistent with the results of the final OS analysis published recently [15]: the final observed HR was 0.87 with a 95 % confidence interval of 0.63–1.21 (58 % of events). The final observed median OS was 17.1 and 13.2 months in both arms, respectively.

Discussion

Maintenance treatment in responders after induction first-line treatment, without waiting for disease progression and start of a new line of therapy, is a therapeutic strategy investigated and used in several tumor types including adult and pediatric acute lymphocytic leukemia [23, 24], follicular non-Hodgkin lymphoma [25, 26], multiple myeloma [27], breast cancer [28], metastatic colorectal cancer [29, 30], and advanced ovarian cancer [3133]. Although well established for certain hematologic cancers, maintenance therapy has only recently become a treatment option for NSCLC [13]. The risk-benefit ratio of maintenance therapy in NSCLC is still unclear, and the thoracic oncology community has seen considerable debate over several aspects of this strategy [2]. Even when maintenance treatment allows prolonging PFS and possibly OS, it is unclear whether OS is prolonged compared to classical first-line followed by second-line paradigm. The selection of patients likely to benefit warrants further research [13].

Model-based approaches are gaining momentum to optimize anti-cancer drug usage and development [6]. Estimates of TGI metrics from modeling of longitudinal tumor size data have been used to predict clinical outcomes and simulate clinical trials [5] in variety of settings including first- and second-line treatment of NSCLC [810]. We present here an adaptation of the modeling framework for maintenance treatment in NSCLC. The framework is developed based on two erlotinib maintenance studies and assessed in simulating outcome of an independent pemetrexed study. As observed in first-line treatment [9, 10], an estimate of time to tumor regrowth (TTG) after start of maintenance treatment captured drug effect, i.e. an OS model incorporating TTG and baseline prognostic factors was able to simulate erlotinib HR in SATURN and ATLAS. Baseline prognostic factors in the model are well known prognostic factors for OS: good prognostic for patients with small baseline tumor size, age below 55 years, Asian ethnicity, ECOG score 0 and for female patients. Smoking status and histology (squamous vs. non-squamous) that were of significant prognosis in the univariate analysis were not retained in the final multivariate model.

As previously discussed [16], the TGI model does not account for exposure to the treatment drugs and is not subjected to any simulation-based assessment (e.g. visual predictive check) because it is not meant to be used for simulation but only to estimate the TGI metrics to be used in the OS model. The TGI model could be in other forms as well, such as a combination of exponential and/or linear models [8, 34] or a simple spline function. Therefore the fundamental assumption of constant exposure over time that was previously used [35] to derive this TGI model from the more complex exposure-driven model is irrelevant here as the model is not used in simulations of response for alternative exposure. There is also no need to assess covariate effects on the TGI model parameters because the model is not used to simulate tumor sizes in new patients.

We performed a two-stage analysis, meaning that we first estimated TGI metrics and then developed the OS model, and we thereby ignored time-dependent hazard driven by time-dependent tumor size. In a typical clinical trial setting, tumor size is only observed until disease progression when treatment is stopped. Median time of last tumor size observation was 11–18 weeks while median OS was four times as long (45–63 weeks) in our model-building dataset. As a result TTG is much shorter than time to death as shown with the median estimates of TTG and OS in Fig. 2. Accounting for tumor size-dependent hazard would have implied an extrapolation substantially beyond last tumor size observation, leading to unrealistically large tumor sizes as the model assumes exponential growth after end of treatment. Information about subsequent treatments is usually unavailable, while tumor size-dependent hazard could only be implemented and evaluated with richer data that could be obtained during routine care of the patients across several lines of treatments when tumor size data could be observed and hazard defined up to patient death. This approach has been explored with PFS, which does not suffer this problem [36]. Additionally, simulations have shown that TTG was not confounded with OS [37, 38].

In the OS model, the censoring model is meant to mimic the duration (treatment plus follow-up period) a patient stays in the study if no death event occurs. The distribution of this duration is defined per protocol by the maximum duration of the study and the patient inclusion rate. If a patient is predicted to die after his predicted duration in the study, this patient is censored. The distribution of study duration is independent of OS and TGI data and doesn’t require simultaneous modeling.

Another limitation of our analysis is that patients needed to have at least two tumor size measurements in the maintenance phase to be evaluable in the TGI model because the TGI parameters were unidentifiable with only one tumor size measurement. These excluded patients who died or dropped out of the study early before the first tumor size measurement may have rapidly growing tumors. However, this may not have a significant impact on this analysis because 94 % of the patients were evaluable.

The model successfully simulated the OS outcomes of the pemetrexed maintenance study AVAPERL based on interim tumor size data collected by the time of PFS database lock before median OS was even reached. This is the first modeling framework for maintenance treatment and one of the few such frameworks validated in simulating an independent study with a drug with a different mechanism of action (pemetrexed) compared to the one used to develop the model (erlotinib), providing support to the hypothesis that TGI metrics capture drug effect independent of treatment [5]. This framework may be used to support design and interim analysis of upcoming maintenance studies and to help in the selection of patients most likely to benefit from maintenance treatment.

Conclusion

In conclusion, a robust TGI-OS model linking OS with TGI metrics and prognostic factors was developed for maintenance therapy following first-line NSCLC treatment. The model successfully predicted the OS outcomes of an independent study (AVAPERL) based on interim tumor size data (up to PFS database lock), indicating that the model may be used for trial simulation and facilitate interpretation of interim data and development decisions. 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 also supported the use of longitudinal tumor size and TTG as endpoints in early clinical oncology studies.

Abbreviations

AIC, Akaike information criterion; ECOG, Eastern Cooperative Oncology Group; ECTS, early change in tumor size; HR, hazard ratio; MTx, maintenance treatment; NSCLC, non-small cell lung cancer; OS, overall survival; PFS, progression-free survival; PI, prediction interval; TGI, tumor growth inhibition; TTG, time to tumor regrowth

References

  1. Sculier JP, Berghmans T, Castaigne C, et al. Maintenance chemotherapy for small cell lung cancer: a critical review of the literature. Lung Cancer. 1998;19:141–51.

    Article  CAS  PubMed  Google Scholar 

  2. Gerber DE, Schiller JH. Maintenance chemotherapy for advanced non-small-cell lung cancer: new life for an old idea. J Clin Oncol. 2013;31(8):1009–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Rosell R, Karachaliou N. Maintenance therapy and precision medicine in NSCLC. Nat Rev Clin Oncol. 2013;10:549–50.

    Article  CAS  PubMed  Google Scholar 

  4. Polo V, Besse B. Maintenance strategies in stage IV non-small-cell lung cancer (NSCLC): in which patients, with which drugs? Ann Oncol. 2014;25(7):1283–93.

    Article  CAS  PubMed  Google Scholar 

  5. Bruno R, Mercier F, Claret L. Evaluation of tumor-size response metrics to predict survival in oncology clinical trials. Clin Pharmacol Ther. 2014;95:386–93.

    Article  CAS  PubMed  Google Scholar 

  6. Venkatakrishnan K, Friberg LE, Ouellet D, Mettetal JT, Stein A, Trocóniz IF, et al. Optimizing oncology therapeutics through quantitative clinical pharmacology: challenges and opportunities. Clin Pharmacol Ther. 2015;97:37–54.

    Article  CAS  PubMed  Google Scholar 

  7. Sharma MR, Gray E, Goldberg RM, et al. Resampling the N9741 trial to compare tumor dynamic versus conventional end points in randomized phase II trials. J Clin Oncol. 2015;33:36–41.

    Article  PubMed  Google Scholar 

  8. Wang Y et al. Elucidation of relationship between tumor size and survival in non-small-cell lung cancer patients can aid early decision making in clinical drug development. Clin Pharmacol Ther. 2009;86:167–74.

    Article  CAS  PubMed  Google Scholar 

  9. Claret L, Mancini P, Sebastien B, Veyrat-Follet C, Bruno R. Model-based estimates of tumor growth inhibition metrics to predict for overall survival in first-line non-small cell lung cance. J Clin Oncol. 2013;31:e19049.

    Google Scholar 

  10. Claret L, Bruno R, Lu JF, Sun YN, Hsu CP. Exploratory modeling and simulation to support development of motesanib in Asian patients with non-small cell lung cancer based on MONET1 study results. Clin Pharmacol Ther. 2014;95:446–51.

    Article  CAS  PubMed  Google Scholar 

  11. Cappuzzo F, Ciuleanu T, Stelmakh L, et al. Erlotinib as maintenance treatment in advanced non-small-cell lung cancer: a multicentre, randomised, placebo controlled phase 3 study. Lancet Oncol. 2010;11:521–9.

    Article  CAS  PubMed  Google Scholar 

  12. Coudert B, Ciuleanu T, Park K, Wu YL, Giaccone G, Brugger W, et al. Survival benefit with erlotinib maintenance therapy in patients with advanced non-small-cell lung cancer (NSCLC) according to response to first-line chemotherapy. Ann Oncol. 2012;23:388–94.

    Article  CAS  PubMed  Google Scholar 

  13. Johnson BE, Kabbinavar F, Fehrenbacher L, Hainsworth J, Kasubhai S, Kressel B, et al. ATLAS: Randomized, double-blind, placebo controlled, phase IIIb trial comparing bevacizumab therapy with or without erlotinib after completion of chemotherapy with bevacizumab for first-line treatment of locally advanced, recurrent, or metastatic non-small cell lung cancer (NSCLC). J Clin Oncol. 2013;31(31):3926–34.

    Article  CAS  PubMed  Google Scholar 

  14. Barlesi F, Scherpereel A, Rittmeyer A, et al. Randomized phase III trial of maintenance bevacizumab with or without pemetrexed after first-line induction with bevacizumab, cisplatin, and pemetrexed in advanced nonsquamous non-small-cell lung cancer: AVAPERL (MO22089). J Clin Oncol. 2013;31:3004–11.

    Article  CAS  PubMed  Google Scholar 

  15. Barlesi F, Scherpereel A, Gorbunova V, Gervais R, Vikström A, Chouaid C, et al. Maintenance bevacizumab-pemetrexed after first-line cisplatin-pemetrexed-bevacizumab for advanced nonsquamous nonsmall-cell lung cancer: updated survival analysis of the AVAPERL (MO22089) randomized phase III trial. Ann Oncol. 2014;25(5):1044–52.

    Article  CAS  PubMed  Google Scholar 

  16. Claret L, Gupta M, Han K, et al. Evaluation of tumor size response metrics to predict overall survival in western and Chinese patients with first line metastatic colorectal cancer. J Clin Oncol. 2013;31:2110–4.

    Article  CAS  PubMed  Google Scholar 

  17. Beal SL, Sheiner LB. NONMEM users guides. San Francisco: NONMEM Project Group, University of California at San Francisco; 1992.

    Google Scholar 

  18. Therasse P, Arbuck SG, Eisenhauer EA, et al. New guidelines to evaluate the response to treatment in solid tumors: European organization for research and treatment of cancer, National cancer Institute of the United States, National cancer Institute of Canada. J Natl Cancer Inst. 2000;92:205–16.

    Article  CAS  PubMed  Google Scholar 

  19. Eisenhauer EA, Therasse P, Bogaerts J, et al. New response evaluation criteria in solid tumours: revised RECIST guideline (version 1.1). Eur J Cancer. 2009;45:228–47.

    Article  CAS  PubMed  Google Scholar 

  20. Savic RM, Karlsson MO. Importance of shrinkage in empirical bayes estimates for diagnostics: problems and solutions. AAPS J. 2009;11(3):558–69.

    Article  PubMed  PubMed Central  Google Scholar 

  21. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014. URL: http://www.R-project.org/.

    Google Scholar 

  22. Akaike H. A new look at the statistical model identification. IEEE Trans Autom Control. 1974;19:716–23.

    Article  Google Scholar 

  23. Ellison RR, Mick R, Cuttner J, et al. The effects of postinduction intensification treatment with cytarabine and daunorubicin in adult acute lymphocytic leukemia: a prospective randomized clinical trial by cancer and leukemia group B. J Clin Oncol. 1991;9:2002–15.

    CAS  PubMed  Google Scholar 

  24. Childhood ALL Collaborative Group. Duration and intensity of maintenance chemotherapy in acute lymphoblastic leukaemia: overview of 42 trials involving 12 000 randomised children. Lancet. 1996;347:1783–8.

    Article  Google Scholar 

  25. Salles G, Seymour JF, Offner F, et al. Rituximab maintenance for 2 years in patients with high tumour burden follicular lymphoma responding to rituximab plus chemotherapy (PRIMA): a phase 3, randomised controlled trial. Lancet. 2011;377:42–51.

    Article  CAS  PubMed  Google Scholar 

  26. Vidal L, Gafter-Gvili A, Salles G, et al. Rituximab maintenance for the treatment of patients with follicular lymphoma: an updated systematic review and meta-analysis of randomized trials. J Natl Cancer Inst. 2011;103:1799–806.

    Article  CAS  PubMed  Google Scholar 

  27. McCarthy PL, Owzar K, Hofmeister CC, et al. Lenalidomide after stem-cell transplantation for multiple myeloma. N Engl J Med. 2012;366:1770–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Muss HB, Case LD, Richards 2nd F, et al. Interrupted versus continuous chemotherapy in patients with metastatic breast cancer: the piedmont oncology association. N Engl J Med. 1991;325:1342–8.

    Article  CAS  PubMed  Google Scholar 

  29. Maughan TS, James RD, Kerr DJ, et al. Comparison of intermittent and continuous palliative chemotherapy for advanced colorectal cancer: a multicentre randomised trial. Lancet. 2003;361:457–64.

    Article  CAS  PubMed  Google Scholar 

  30. Giuliani F, De Vita F, Colucci G, et al. Maintenance therapy in colon cancer. Cancer Treat Rev. 2010;36 suppl 3:S42–5.

    Article  CAS  PubMed  Google Scholar 

  31. Markman M, Liu PY, Wilczynski S, et al. Phase III randomized trial of 12 versus 3 months of maintenance paclitaxel in patients with advanced ovarian cancer after complete response to platinum and paclitaxel-based chemotherapy: a southwest oncology group and gynecologic oncology group trial. J Clin Oncol. 2003;21:2460–5.

    Article  CAS  PubMed  Google Scholar 

  32. Burger RA, Brady MF, Bookman MA, et al. Incorporation of bevacizumab in the primary treatment of ovarian cancer. N Engl J Med. 2011;365:2473–83.

    Article  CAS  PubMed  Google Scholar 

  33. Perren TJ, Swart AM, Pfisterer J, et al. A phase 3 trial of bevacizumab in ovarian cancer. N Engl J Med. 2011;365:2484–96.

    Article  CAS  PubMed  Google Scholar 

  34. Stein WD, Gulley JL, Schlom J, Madan RA, Dahut W, Figg WD, et al. Tumor regression and growth rates determined in five intramural NCI prostate cancer trials: the growth rate constant as an indicator of therapeutic efficacy. Clin Cancer Res. 2011;17(4):907–17.

    Article  CAS  PubMed  Google Scholar 

  35. Claret L, Girard P, Hoff PM, Van Cutsem E, Zuideveld KP, Jorga K, et al. Model-based prediction of phase III overall survival in colorectal cancer on the basis of phase II tumor dynamics. J Clin Oncol. 2009;27(25):4103–8.

    Article  CAS  PubMed  Google Scholar 

  36. Claret L, Bruno R. Wang B, Marchand M, Li C, Girish S, et al. PAGE 24 (2015) Abstr 3391 www.page-meeting.org/?abstract=3391

  37. Claret L, Bruno R. PAGE 23 (2014) Abstr 3109 www.page-meeting.org/?abstract=3109

  38. Claret L, Bruno R. Assessment of tumor growth inhibition metrics to predict overall survival. Clin Pharmacol Ther. 2014;96(2):135–7.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

We acknowledge Bert Lum from Genentech Inc. for assistance.

Funding

Genentech Inc.

Availability of Data and Materials

The data supporting the findings in this article can be found in respective publications [1115].

Authors’ contributions

KH has contributed substantially to conception, design, acquisition of data, analysis, interpretation of data, drafting or revising the paper as well as giving final approval for submission. LC has contributed substantially to conception, design, acquisition of data, analysis, interpretation of data, drafting or revising the paper as well as giving final approval for submission. AS contributed substantially to conception, design, acquisition of data, analysis, interpretation of data, drafting or revising the paper as well as giving final approval for submission. AD contributed substantially to conception, design, acquisition of data, analysis, interpretation of data, drafting or revising the paper as well as giving final approval for submission. JJ contributed substantially to conception, design, acquisition of data, analysis, interpretation of data, drafting or revising the paper as well as giving final approval for submission. RB has contributed substantially to conception, design, acquisition of data, analysis, interpretation of data, drafting or revising the paper as well as giving final approval for submission.

Competing interest

Kelong Han, Alan Sandler, Asha Das, and Jin Jin are employees of Genentech and hold stocks of Roche Pharmaceuticals; In addition, Jin Jin also holds stocks in Eli Lilly; Laurent Claret and Rene Bruno are employees of Pharsight Consulting Services, Pharsight, a CertaraTM Company, and were paid consultants to Pfizer in connection with the development of this manuscript.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The studies complied with the Declaration of Helsinki and Good Clinical Practice guidelines, and were approved at all investigating centers by local ethics committees. All patients provided written informed consent for participation and publication of the data [1115]. An ethics statement was not required for this analysis as they have been provided in each of the three individual studies [1115].

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Kelong Han or Rene Bruno.

Additional files

Additional file 1: Figure S1.

Goodness of fit plot of the simplified tumor growth inhibition (TGI) model. (TIFF 3241 kb)

Additional file 2: Figure S2.

Model predicted vs. observed tumor size for 16 individuals taken at random. Black solid points: observed tumor size. Orange solid line: model predicted tumor size. Grey dash line: time of randomization (start of maintenance treatment). Grey solid line: estimated time to tumor regrowth (TTG). (TIFF 3041 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Han, K., Claret, L., Sandler, A. et al. Modeling and simulation of maintenance treatment in first-line non-small cell lung cancer with external validation. BMC Cancer 16, 473 (2016). https://doi.org/10.1186/s12885-016-2455-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-016-2455-2

Keywords