Abstract
Aim: To investigate the prognostic value of modeled CA-125 kinetic parameters regarding surgery outcomes and time to second-line chemotherapy in a population of ovarian cancer patients treated with neoadjuvant chemotherapy followed by interval cytoreduction. Patients and Methods: This retrospective study included consecutive patients with FIGO stage IIIc/IV ovarian cancer treated between 2006 and 2014. We characterized CA-125 kinetics and identified modeled kinetic parameters. Results: Fifty four patients were included. KELIM (modeled CA-125 elimination rate constant) was an independent predictive parameter of optimal cytoreduction (OR=0.18; 95% CI=0.04-0.69; p=0.02). In the optimal cytoreduction population (40 patients, 74.1%), E50 (concentration producing 50% of the maximum chemotherapy effect) was a significant prognostic parameter regarding time to second-line chemotherapy (HR=0.38; 95% CI=0.173-0.854; p=0.019). Conclusion: We identified CA-125 modeled kinetic parameters during neoadjuvant chemotherapy harboring potential predictive values regarding the likelihood of optimal cytoreduction, along with time to second-line chemotherapy in optimally-cytoreduced patients.
- Epithelial ovarian cancer
- neoadjuvant chemotherapy
- CA-125 kinetics
- mathematical modeling
- cytoreduction prediction
- second-line chemotherapy
Epithelial ovarian cancer (EOC) is the seventh most common cancer in the world among women. Moreover it is the eighth cause of cancer-related deaths (1). Sixty to 75% of patients have advanced disease with peritoneal involvement at diagnosis (2, 3), and only 20-25% can be expected to be long-term survivors.
The standard therapy of advanced epithelial ovarian cancer relies on primary debulking surgery followed by adjuvant chemotherapy when optimal resection is possible (4). Post-operative residual disease is one of the most important prognostic factors regarding survival (5-7). Unfortunately, optimal debulking surgery is not feasible in patients with too extensive involvement of peritoneum and gastro-intestinal structures. The use of NAC (neoadjuvant chemotherapy) before IDS (interval debulking surgery) has been developed to reduce the disease bulk and resections, as well as to select patients with chemosensitive tumors. Based on the data of recent trials suggesting the efficacy of NAC before IDS regarding progression free survival and overall survival (6, 8), American Society of Clinical Oncology has recognized this strategy as a standard strategy in patients with a high perioperative risk profile or a low likelihood of achieving cytoreduction to ≤1 cm of residual disease (5).
However, following an apparently effective NAC, it is common to discover during the surgery procedure that optimal resection is finally not possible (6, 8). Among multiple reasons for such finding, the current imagining tests lack of accuracy for predicting tumor response to treatments, due to the small sizes of peritoneal carcinomatosis lesions (9). As a consequence, there is a need for non-invasive tools able to predict the chances of resectability after NAC.
Some authors have developed predictive models based on clinical characteristics, computed tomography imaging, or CA-125 levels (10, 11). Most of these approaches present strong limitations, and none of them have been yet validated. Mathematical modeling of CA-125 kinetics during neoadjuvant chemotherapy may be a promising strategy. Indeed, this modern approach presents several advantages, such as low dependence on selected timepoints, longitudinal analysis of tumor marker kinetics, along with consideration of inter- and intra-individual variability (12). The reproducible predictive values of modeled kinetic parameters of other serum tumor markers, such as CA-125, prostate-specific antigen (PSA), human chorionic gonadotropin (hCG), alpha-fetoprotein (AFP) and circulating tumor cells (CTCs) have been reported in previous works (13-17). A modeled kinetic parameter related to CA-125 elimination, KELIM (modeled CA-125 elimination rate constant), was reported with a strong prognostic value regarding progression free survival in patients with recurrent ovarian cancer treated with carboplatin-based chemotherapy regimens in CALYPSO phase III trial (17). The prognostic value of KELIM was then validated in 3 independent phase III databases in patients treated with 1st line chemotherapy (13, 14, 16, 17).
We hypothesize that mathematical modeling of CA-125 kinetics during NAC may help predict the chance of optimal resection and time to second-line chemotherapy in ovarian cancers patients. The aim of the present exploratory retrospective study was to investigate the prognostic value of modeled CA-125 kinetic parameters regarding surgery outcomes.
Patients and Methods
The primary objective of the study was to characterize CA-125 kinetics with a kinetic semi-mechanistic model. The secondary objective was to identify modeled kinetic parameters which might harbor predictive values regarding the chance of optimal cytoreductive surgery, and longer time to second-line chemotherapy.
The data of patients consecutively treated in a gynecology department for pathologically proven stage IIIc or IV FIGO ovarian cancers with NAC followed by IDS between 2006 and 2014 were collected. To build a semi-mechanistic model of CA-125 kinetics, we included patients with at least 1 abnormal CA-125 level (>35 IU/ml) measured within 8 days before the first chemotherapy cycle, and at least 3 CA-125 level determinations during NAC. We excluded all patients with unavailable CA-125 level data, or incomplete clinicopathological data. Age at diagnosis, tumor histology, chemotherapy type, number of cycles, cycles dates, CA-125 level monitoring, surgery date and outcomes (cytoreductive surgery: optimal cytoreduction or not), were collected. Surgical procedures were all performed by different skilled gynecology oncology surgeons, considered as experts in radical surgical procedures, with the intent to achieve optimal cytoreduction. Surgical indications were validated in multidisciplinary meetings. All tumors were considered resectable based on pre-operative imaging. Optimal cytoreductive surgery was defined according to the 2009 Gynecologic Oncology group (GOG) criteria (18), based on the size of residual tumor nodules, as it was the reference during this period. Surgery was considered optimal if the residual tumor nodules did not exceed 1 cm in maximum diameter. The procedure included at least total abdominal hysterectomy and bilateral salpingo-oophorectomy, supra-mesocolic omentectomy, pelvic and para-aortic lymph node dissection. The neoadjuvant chemotherapy regimens were determined by medical oncologists based on international recommendations.
Individual CA-125 data were analyzed using a population kinetic approach with a non-linear mixed effect model (19, 20), and were Log transformed. Basic details of population kinetic approach have been presented elsewhere (12). A semi-mechanistic model based on kinetic-pharmacodynamic (K-PD) approach was used to fit serum CA-125 values available before IDS (21). Indeed, K-PD models are commonly used in the absence of pharmacokinetic (PK) data for pharmacokinetic studies. As shown in Figure 1, treatment kinetic was described using a 2-virtual-compartment model, where A1 represents chemotherapy in the central compartment, A2 chemotherapy in the transit compartment and A3 the CA-125 rate. The inhibition of CA-125 production induced by treatment was expressed by an indirect effect model using an Emax relationship (22). K is the treatment kinetic rate constant (years−1); KPROD is the CA-125 basal tumor production rate (UI ml−1 per year); E50 is the concentration producing 50% of the maximum effect (arbitrary unit; AU); KELIM is the CA-125 elimination rate (years−1); and CA-1250 is CA-125 baseline determined by the model (UI ml−1).
To validate the best model, we assessed standard errors of estimated parameters, goodness-of-fit plots and distribution of the normalized prediction distribution errors (NPDE). In addition CA-125 decline profiles were simulated in 43,200 virtual patients (800 simulations) using the final individual parameters estimated from the best model to perform a visual predictive check (VPC, i.e., 10th and 90th percentiles of the observed CA-125 data values were compared to a 95% confidence interval [CI] of CA-125 titers from the 800 simulated replicates). We used VPC to assess the internal validity of the model. The predictive values of modeled kinetics parameters (K; KPROD; E50; KELIM; CA-1250) categorized by their medians were tested regarding outcomes using univariate and multivariate analyses. Moreover, predictive CA-125 cut-offs previously reported in the literature were also investigated: CA125 value ≤ or >75 IU/ml after the third cycle of chemotherapy (CA-125Post C3 and CA-125Post C3MOD); and CA125 value ≤ or >20 IU/ml before IDS (CA-125Pre IDS and CA-125PreIDSMOD) (10, 11). Multivariate logistic regressions were used to assess their predictive values regarding surgery outcomes (i.e. optimal cytoreduction: yes or no).
Kaplan-Meier curves and Log-rank tests were used to assess the predictive values of these parameters regarding time to second-line chemotherapy in patients with optimal cytoreduction using univariate analyses. All variables significant with p≤0.3 were included in the multivariate Cox model. A landmark time to second-line chemotherapy in the optimally cytoreduced population was used with landmark time point at IDS.
The NONNEM 7.3.0 software (ICON Development Solutions, Ellicott city, MD, USA) was used to fit CA-125 kinetics to the semi-mechanistic model using the SAEM method. Data handling, survival analyses and graphical representations were performed in R © 3.2.1 software. All tests were implemented using a two-sided 0.05 alpha risk.
Results
The data from 54 patients treated with NAC between 2006 and 2014 were collected and analyzed for modeling purpose (Figure 2). The median follow-up was 28.1 months (range=9.0-117.5 months). The patient clinical characteristics are summarized in Table I. Optimal IDS was achieved in 40 patients (74.1%). However 14 patients were not optimally cytoreduced (25.9%) because 5 had massive diffuse peritoneal carcinomatosis, 4 had miliary carcinomatosis, 3 had iliac fossa or parametrium deep lymph node involvement, 1 had pleural carcinomatosis and 1 had hepatic pedicle infiltration. Mean age was 63 years for optimal IDS, and 67 years for non-optimal IDS. Thirty three (82.5%) and 10 patients (71.4%) had stage IIIc diseases in optimal and in non-optimal surgery groups, respectively. The median CA-125 levels before the first cycle of chemotherapy was 1,283 IU/ml for optimal IDS (range=37-5,966 IU/ml) against 3,370 IU/ml for non-optimal IDS (range=73-17,621 IU/ml). All patients who underwent optimal surgery had received taxane and platinum-based NAC. Five optimally cytoreduced patients had been in CHIVA trial assessing the efficacy of an anti-angiogenic drug (nintedanib) in addition to first line pre-operative chemotherapy. The median number of neoadjuvant cycles was 5.0 (range=3-7) and 5.4 (range=3-7) in optimal and non-optimal IDS groups respectively.
Typical parameters and between subject variability (BSV) estimated from the semi-mechanistic model for CA-125 kinetics are presented in Table II. Relative standard errors (RSE) of typical parameters were all less than 40%, excepted for KPROD. The goodness-of-fit plots demonstrated that CA-125 kinetics were properly fitted by the model before IDS. Observations versus individual predictions were close to the identity line (Figure 3a). NPDE did not deviate from a standard normal distribution and suggested that the model and population parameters distributions were correct (Figure 3b). VPC showed that most of the CA-125 values were included within the 90% confidence interval boundaries of simulated CA-125 predictions (Figure 3c of Figure 3), thereby suggesting the good predictive performance of the model.
Using univariate logistic analysis, CA-125 elimination rate parameter KELIM exhibited a strong predictive value regarding optimal IRS (Table III) (odds ratio (OR)=0.18, 95% confidence interval [0.04 to 0.69]; p=0.02). The other modeled kinetic parameters were not significant, including treatment kinetic rate constant K (p=0.21), CA-125 tumor production rate KPROD (p=0.07), concentration producing 50% of the maximum effect E50 (p=0.53) and CA-125 baseline CA-1250 (p=0.21). Among covariates based on CA-125 thresholds, only CA-125Post C3, was found significant (OR=4.15; 95% C=1.15-15.96; p=0.03). Using multivariate analysis performed including parameters with p≤30% significance in univariate analyses (KELIM, KPROD, CA-125Post C3), KELIM was the only independent parameter regarding the likelihood of optimal cytoreduction (OR=0.18; 95% CI=0.04-0.69; p=0.02).
The data of 40 patients who underwent optimal cytoreduction surgery were analyzed. Using Kaplan-Meier curves, 2 modeled kinetic parameters had strong predictive values: E50 (11.7 months if E50 is < median vs. 25.0 months if E50 ≥ median, p=0.015) (Figure 4 and Table III); baseline CA-125 (29.5 months if CA-1250 < median vs. 13.2 months if CA-1250 ≥median, p=0.03). In a multivariate Cox regression model including these significant prognostic factors, E50 was the only significant covariate (HR=0.38; 95% CI=0.173-0.854; p=0.019).
Discussion
The residual tumor size after cytoreductive surgery is one of the most important prognostic factors regarding survival in EOC (5-7). The GOG defined optimal surgery as residual tumor nodules measuring 1 cm or less in maximum diameter (18). In our population, 74.1% of patients who underwent IDS had residual nodules tumors <1 cm after surgery, with the intent to achieve optimal cytoreduction. Data from clinical trials subsequently suggested that resection to microscopic tumor burden was associated with improved survival time (23) and gynecology oncologists now consider that the aim of cytoreductive surgery should ideally be no residual disease.
Although pre-operative prediction of optimal surgery could guide treatment decisions, few studies reported approaches able to forecast the amount of residual disease after cytoreductive surgery and NAC. Some authors have attempted to use computed tomography imaging to predict the amount of residual tumor after cytoreductive surgery. However, the radiologic models had poor reproducibility for predicting any residual disease after surgery, due to the small and heterogeneous size of peritoneal lesions (9). Other authors tried to define CA-125 thresholds to predict optimal cytoreduction. Furukawa et al. reported a 20 IU/ml CA-125 as a specific and sensitive cut-off before IDS (sensitivity 74.4%, and specificity 82.1%) (10). Pelissier et al. determined that a CA-125 value of less than 75 IU/ml after 3 cycles of chemotherapy was an independent predictor of complete IDS (sensitivity: 52% [40%-63%], specificity: 68% [62%-74%]; p=0.02) in a retrospective study with 148 patients (11). These covariates were tested in our study, and only a CA-125 value of less than 75 IU/ml after 3 cycles of chemotherapy was significant using univariate analysis (OR=4.15; 95% C=1.15-15.96; p=0.03). The limited reproducibility of CA-125 thresholds is probably due to the fact that they were arbitrarily set by authors after comparing the predictive values of different cut-offs. These simple static approaches rely on one time-point, without consideration of intra and inter-individual variabilities related to tumor biomarker assays (12). None of them have been shown to be reproducible yet.
Integration of several CA125 measurements in mathematical models reduces the effects of intra-individual variability by assessing the equations describing the whole decline curves, on an individual basis. In previous studies involving recurrent ovarian cancer patients, KELIM, a kinetic parameter describing CA125 elimination rate, was a strong independent predictive value of tumor lesion size changes (13), progression-free survival (17) and overall survival in independent phase III trial databases (24). Consistently, we analyzed the roles of modeled CA-125 kinetic parameters as potential predictors of optimal cytoreduction after NAC in the present study (17). KELIM, already reported previously, was found to be an independent predictive factor of optimal cytoreduction. Moreover, another modeled kinetic parameter E50 was a strong predictive factor of time to second-line chemotherapy in optimally cytoreduced patients. E50 represents the concentration producing 50% of the maximum effect of chemotherapy. As a consequence, it somehow represents the strength of chemotherapy, and suggests that initial chemotherapy efficiency might explain a part of the recurrence risk in such patients.
There are some limitations in our study. It is a retrospective study with a limited number of patients treated in a single institute. IDS and chemotherapy were validated in multidisciplinary meetings based on international recommendations that evolved during the study period. The PCI and Fagotti modified score are now routinely used to describe the extent of peritoneal spread, and the resectability of ovarian cancer (25-27). It was not the case before 2010, thereby explaining why data about peritoneal cancer index (PCI) could not be extracted for the present work. In the present real-life study, cytoreductive surgery was usually performed after 3 cycles of chemotherapy, as assessed in the CHORUS trial (8). However, patients were sometimes operated after 6 cycles, when tumor response was insufficient for optimal surgery, thereby inducing some heterogeneity in the cohort data.
Despite these limitations, the outcomes of the present exploratory study suggest that modeled kinetic parameters may help identify patients who might have optimal cytoreduction after neoadjuvant chemotherapy, and shorter time to second-line chemotherapy in optimally cytoreduced patients. Additional confirmatory studies are needed, and we have planned to i) assess the predictive value of KELIM and E50 in the CHIVA randomized phase II trial dababase, and ii) to perform a prospective study to validate the value of these modeled kinetic parameters.
Acknowledgements
The Authors thank the physicians and staffs who were involved in patient management and who contributed to generate the data.
- Received August 31, 2017.
- Revision received October 18, 2017.
- Accepted October 25, 2017.
- Copyright© 2017, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved