A distinct p53 target gene set predicts for response to the selective p53–HDM2 inhibitor NVP-CGM097
Figures
![](https://iiif.elifesciences.org/lax/06498%2Felife-06498-fig1-v3.tif/full/617,/0/default.jpg)
TP53 wild-type status is necessary but not sufficient to predict sensitivity to NVP-CGM097 and NVP-CFC218.
(A) Chemical structure of NVP-CGM097 and NVP-CFC218. (B) In vitro activity of NVP-CFC218 and NVP-CGM097 in TR-FRET binding assay (a) and cellular proliferation assay in human cancer cell lines (b). Data are expressed as concentration causing 50% inhibition and shown as mean ± SD from multiple (n ≥ 8) independent experiments. (c) Selectivity is determined by the ratio of GI50 obtained using the HCT-116 p53-null and the HCT-116 p53WT isogenic pair of cell lines. (d) Selectivity is determined by the ratio of GI50 obtained using SAOS-2 (p53-null) and SJSA-1 (p53WT and HDM2-amplified) osteosarcoma pair of cell lines. (C and D) Scatter plot showing IC50 values expressed in μM of NVP-CFC218 in cell viability assays of p53 wild-type cell lines (C) and p53 mutated cell lines (D), colored by their response to NVP-CFC218. The data used to generate these plots, as well as cell line identity is available in Figure 1—source data 1. (E) Contingency table indicating the total number of sensitive and insensitive cell lines to NVP-CFC218. The p-value of 5.1 × 10−23 shows a significant association between sensitivity to NVP-CFC218 and TP53 wild-type status. (F) Main enriched compound target p-values from the Global Compound Selectivity Analysis. p-values are minus log10 transformed. The red color refers to compounds that are more selective in the wild-type p53 (p53WT) strata than in the mutated p53 (p53MUT) strata. The blue color indicates the reverse profile. Brighter colors indicate which target classes pass the 0.25 FDR cut-off. The length of each red/blue segment corresponds to the proportion of p53WT selective/p53MUT selective compounds in each target class.
-
Figure 1—source data 1
List of cell lines tested for their sensitivity to NVP-CFC218 (n = 356).
- https://doi.org/10.7554/eLife.06498.004
![](https://iiif.elifesciences.org/lax/06498%2Felife-06498-fig2-v3.tif/full/617,/0/default.jpg)
Predictive modeling of NVP-CFC218 chemical sensitivity from CCLE expression data.
(A) Prediction accuracy estimation by bootstrapping analysis upon increasing feature set size. Estimates are averaged over 20 bootstrapping repeats. The maximum accuracy is observed for 17 features (red dot). Sensitivity and specificity are shown in Figure 2—figure supplement 1. The averaged data used to generate the accuracy, sensitivity and specificity plots is available in Figure 2—source data 1. (B) Feature occurrence frequencies of 13 feature selections in 100 bootstrapping repeats. Fifty five features were selected at least once. The 13 most frequently occurring features (red dots) were selected more than 30 times. (C) ROC curves for the three predictive models under comparison, i.e., the p53 mutation status, the 215-feature and the 13-gene signature models. (D) Precision-Recall plot for the 13-gene signature. Five curves are typically shown since cross-validation was repeated five times. (E) Performance estimates of the three compared predictive models: AUC (Area Under the Curve, from the ROC curve shown in C), Sensitivity (fraction of correctly predicted sensitive cell lines), Specificity (fraction of correctly predicted insensitive cell lines), PPV (positive predicted value, fraction of sensitive cell lines predicted as such), and NPV (negative predictive value, fraction of insensitive cell lines predicted as such). Measures are averaged over the 5 iterations of 5-fold cross-validations.
-
Figure 2—source data 1
Classifier performance with increasing feature set size.
- https://doi.org/10.7554/eLife.06498.006
![](https://iiif.elifesciences.org/lax/06498%2Felife-06498-fig2-figsupp1-v3.tif/full/617,/0/default.jpg)
Prediction performance estimation by bootstrapping analysis upon increasing feature set size.
Estimates are averaged over 20 bootstrapping repeats. Maximum sensitivity is observed for 13 features (sensitivity plot, red dot). Maximum specificity is observed for 17 features (specificity plot, red dot).
![](https://iiif.elifesciences.org/lax/06498%2Felife-06498-fig2-figsupp2-v3.tif/full/617,/0/default.jpg)
Prediction performance of NVP-CFC218 in p53WT CCLE cell lines.
The measures of prediction accuracy are calculated restricted to the wild-type TP53 subset of CCLE cell lines (n=68). (A) ROC curves for the 215-feature set and the 13-gene signature predictive models. (B) Precision-Recall plot for the 13-gene signature. Five curves are typically shown since cross-validation was repeated 5 times. (C) Performance estimates of the 13-gene signature and TP53 mutation status as predictive models.
![](https://iiif.elifesciences.org/lax/06498%2Felife-06498-fig3-v3.tif/full/617,/0/default.jpg)
Validation of the predictive 13-gene signature in an external set of cell lines.
A total of 52 or 38 cell lines were tested against NVP-CGM097 (A) or NVP-CFC218 (B), respectively, in a 3-day proliferation assay. Sensitivity call for each compound was applied according to a cut-off of 3 µM. Cells that were predicted as sensitive are shown in pink, while the ones predicted as insensitive are shown in blue. The data used to generate the waterfall plots and cell line details are available in Figure 3—source data 1. Performance values of the 13-gene signature were derived from this experiment, and are shown in (C).
-
Figure 3—source data 1
Sensitivity prediction and sensitivity to NVP-CFC218 and NVP-CGM097 of an external set of cell lines (n = 52).
- https://doi.org/10.7554/eLife.06498.011
![](https://iiif.elifesciences.org/lax/06498%2Felife-06498-fig3-figsupp1-v3.tif/full/617,/0/default.jpg)
Validation of the predictive 13-gene signature in p53WT cell lines.
The same visualization as in Figure 3 is shown restricted to the wild-type p53 cell lines. A total of 42 or 30 p53WT cell lines were tested against NVP-CGM097 (A) or NVP-CFC218 (B), respectively, in a 3-day proliferation assay. The data used to generate the waterfall plots and cell line details are available in Figure 3—source data 1. Performance values of the 13-gene signature were derived from this experiment, and are shown in (C).
![](https://iiif.elifesciences.org/lax/06498%2Felife-06498-fig4-v3.tif/full/617,/0/default.jpg)
Prevalence of the 13-gene signature in human primary tumors and patient-derived tumor xenograft models.
(A) Fraction of human primary tumors predicted sensitive by the 13-gene signature. Samples are grouped by primary site. Expression data from 21,300 human primary tumors were compiled and only primary sites consisting of more than 50 samples were included in the analysis. The primary site nomenclature is based on the COSMIC Classification System. (B) Correlation between the fraction of human primary tumor samples predicted sensitive and the observed sensitivity ratio in the CCLE cell line data. (C) Correlation between the sensitivity predictions from the human primary tumor sample collection and the patient-derived tumor xenograft model collection. Human primary tumor samples, patient-derived tumor xenografts, and CCLE cell lines are organized by lineage. The dashed line corresponds to the identity line.
![](https://iiif.elifesciences.org/lax/06498%2Felife-06498-fig5-v3.tif/full/617,/0/default.jpg)
Activity of NVP-CGM097 in SJSA-1 tumor-bearing mice.
(A) PK/PD relationship. Mice (n = 3/time-point) bearing established subcutaneous SJSA-1 xenografts received vehicle or a single oral dose of 100 mg/kg NVP-CGM097. Levels of NVP-CGM097 were measured in plasma and in tumor at 7 different time-points within 48 hr following the dose. As a pharmacodynamic readout, the p53 target gene p21 mRNA (left panel) and protein (right panel) levels were assessed in the tumor samples by qRT-PCR and reverse phase protein array at 7 time-points within 48 hr following NVP-CGM097 single dose. Data are plotted as mean ± SEM. (B) In vivo anti-tumor activity of NVP-CGM097. Mice (n = 6/dosing group) bearing established subcutaneous SJSA-1 xenografts received vehicle or 25, 50, or 100 mg/kg of an oral suspension of NVP-CGM097 daily. Tumor volumes were calipered throughout the study, and data are plotted as mean ± SEM (left panel). The body-weight (BW) is expressed as the percentage change in BW relative to day 0 (initiation of treatment). Data are plotted as mean ± SEM (right panel). *, p < 0.05 vs vehicle control. (C and D) Tumor pharmacodynamics effects of NVP-CGM097. Tumor samples were retrieved 3 hr post-last dose at the end of the efficacy study described in (B). Tumor p21 mRNA levels were assessed by qRT-PCR and data are plotted as mean ± SEM (C). Tumor samples were retrieved and paraffin-embedded 3 hr post-last dose at the end of the efficacy experiment described in (B). p21 and HDM2 protein levels were assessed by immunohistochemistry (D).
![](https://iiif.elifesciences.org/lax/06498%2Felife-06498-fig6-v3.tif/full/617,/0/default.jpg)
Validation of the predictive 13-gene signature in patient-derived tumor xenografts (PDX).
PDX models were established by subcutaneous implantation in nude mice of surgical tumor tissues from treatment-naive cancer patients. 55 PDX models of multiple cancer types were used for the study: melanoma (n = 19), colorectal cancer (n = 17), liposarcoma (n = 3), renal cell carcinoma (n = 3), hepatocellular carcinoma (n = 7), breast cancer (n = 1), pancreatic cancer (n = 2), and lung cancer (n = 3). Tumor-bearing mice (n = 4/dosing group/model) received vehicle or 100 mg/kg of NVP-CGM097 daily for 4 weeks. The response is reported as percentage change in tumor volume at a given day of treatment relative to day 0 (start of treatment). The cut-off used for sensitivity to NVP-CGM097 was based on RECIST adapted for full tumor volume measurement. The sensitivity call was made at the maximum effect time-point during the treatment period. Sensitivity predictions for each PDX model were generated using the 13-gene signature (Figure 6—source data 2). (A) Waterfall plot showing the tumor response to NVP-CGM097 for all the PDX models in the study (n = 55), color-coded by the prediction output. (B) The same visualization as in (A) is shown restricted to the wild-type p53 PDX models (n = 33). The data used to generate the waterfall plots and PDX model details are available in Figure 6—source data 1. Performance values of the 13-gene signature were derived from this experiment and are shown in (C).
-
Figure 6—source data 1
Sensitivity prediction and sensitivity to NVP-CGM097 of a set of in vivo PDX models (n = 55).
- https://doi.org/10.7554/eLife.06498.017
-
Figure 6—source data 2
Expression data of each of the 13 genes included in the gene signature in the set of in vivo PDX models (n = 55).
- https://doi.org/10.7554/eLife.06498.018
![](https://iiif.elifesciences.org/lax/06498%2Felife-06498-resp-fig1-v3.tif/full/617,/0/default.jpg)
Tables
13-gene signature selected for the sensitivity prediction to p53–HDM2 inhibitors
Features | Fold change | Wilcoxon p-value | Proposed function in p53 pathway |
---|---|---|---|
Expr. HDM2 | 2.18 | 3.09 × 10−14 | Negative feedback loop |
Expr. CDKN1A | 3.88 | 3.76 × 10−11 | Cell cycle/senescence |
Expr. ZMAT3 | 2.81 | 1.20 × 10−10 | Positive feedback loop |
Expr. DDB2 | 2.55 | 1.22 × 10−10 | DNA repair |
Expr. FDXR | 2.42 | 1.08 × 10−09 | Apoptosis |
Expr. RPS27L | 1.97 | 1.23 × 10−09 | Positive feedback loop |
Expr. BAX | 2.12 | 1.84 × 10−09 | Apoptosis |
Expr. RRM2B | 2.06 | 4.25 × 10−09 | DNA repair |
Expr. SESN1 | 2.27 | 1.04 × 10−08 | Oxidative stress |
Expr. CCNG1 | 1.69 | 4.81 × 10−08 | Cell cycle |
Expr. XPC | 1.62 | 1.56 × 10−07 | DNA repair |
Expr. TNFRSF10B | 1.91 | 3.30 × 10−06 | Apoptosis |
Expr. AEN | 1.46 | 3.30 × 10−06 | Apoptosis |
-
Expr: gene-level expression values generated by the Affymetrix GeneChip technology with the HG-U133 plus 2 arrays, summarized according to the RMA normalization method.
PK parameters of NVP-CFC218 and NVP-CGM097 in preclinical species
PK parameters | Mouse | Rat | Dog | Monkey |
---|---|---|---|---|
NVP-CFC218 | ||||
CL (mL/min.kg) | 9 | 12 ± 1 | 3 ± 1 | 11 ± 2 |
Vss (L/kg) | 3.1 | 8.0 ± 1.8 | 3.6 ± 0.5 | 2.1 ± 0.3 |
t1/2app. (h) | 3.6 | 9.3 ± 1.9 | 14.4 ± 0.7 | 2.6 ± 0.3 |
AUCinf (µM.h) i.v.* | 2.80 | 2.21 ± 0.21 | 8.08 ± 1.68 | 2.37 ± 0.31 |
AUCinf (µM.h) p.o.* | 1.20 | 1.07 ± 0.34 | 5.74 ± 0.68 | 0.34 ± 0.05 |
Cmax (µM) p.o.* | 0.094 | 0.075 ± 0.023 | 0.240 ± 0.043 | 0.077 ± 0.014 |
Tmax p.o. (h) | 2.0 | 4.0 ± 1.4 | 3.3 ± 1.2 | 1.0 ± 0.1 |
Oral BA (%F) | 43 | 49 ± 15 | 71 ± 8 | 14 ± 2 |
NVP-CGM097 | ||||
CL (mL/min.kg) | 5 | 7 ± 1 | 3 ± 1 | 4 ± 1 |
Vss (L/kg) | 2.6 | 6.4 ± 0.4 | 3.8 ± 0.4 | 2.0 ± 0.4 |
t1/2app. (h) | 6.4 | 12.1 ± 1.1 | 14.2 ± 1.8 | 8.3 ± 0.9 |
AUCinf (µM.h) i.v.* | 4.70 | 3.66 ± 0.27 | 9.44 ± 1.43 | 6.59 ± 0.74 |
AUCinf (µM.h) p.o.* | 3.34 | 2.97 ± 0.40 | 7.08 ± 1.05 | 3.73 ± 0.78 |
Cmax (µM) p.o.* | 0.207 | 0.134 ± 0.016 | 0.250 ± 0.030 | 0.363 ± 0.020 |
Tmax p.o. (h) | 4.0 | 4.5 ± 1.9 | 2.7 ± 1.2 | 1.3 ± 0.6 |
Oral BA (%F) | 71 | 81 ± 11 | 75 ± 11 | 57 ± 12 |
-
Mouse (n = 1 per time point) and rat (n = 3 per time point) were treated either with a single dose of the indicated compound at 1 mg/kg i.v. or 3 mg/kg p.o. Dog and monkey (n = 3 per time-point) were treated either with a single dose of the indicated compound at 0.1 mg/kg i.v. or 0.3 mg/kg p.o.
-
*AUC and Cmax values are shown dose-normalized to 1 mg/kg.