
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Research Articles: Therapeutics, Targets, and Development
Prediction of drug combination chemosensitivity in human bladder cancer
Departments of 1 Molecular Physiology and 2 Public Health Sciences and 3 Mellon Prostate Cancer Institute, University of Virginia, Charlottesville, Virginia; 4 Department of Statistics, Korea University, Seoul, Korea; and 5 Genomics Institute of the Novartis Research Foundation, San Diego, California
Requests for reprints: Dan Theodorescu, Department of Molecular Physiology and Biological Physics, University of Virginia Health Sciences Center, Box 422, Charlottesville, VA 22908. Phone: 434-924-0042; Fax: 434-982-3652. E-mail: dt9d{at}virginia.edu
Abstract
The choice of therapy for metastatic cancer is largely empirical because of a lack of chemosensitivity prediction for available combination chemotherapeutic regimens. Here, we identify molecular models of bladder carcinoma chemosensitivity based on gene expression for three widely used chemotherapeutic agents: cisplatin, paclitaxel, and gemcitabine. We measured the growth inhibition elicited by these three agents in a series of 40 human urothelial cancer cell lines and correlated the GI50 (50% of growth inhibition) values with quantitative measures of global gene expression to derive models of chemosensitivity using a misclassification-penalized posterior approach. The misclassification-penalized posteriorderived models predicted the growth response of human bladder cancer cell lines to each of the three agents with sensitivities of between 0.93 and 0.96. We then developed an in silico approach to predict the cellular growth responses for each of these agents in the clinically relevant two-agent combinations. These predictions were prospectively evaluated on a series of 15 randomly chosen bladder carcinoma cell lines. Overall, 80% of the predicted combinations were correct (P = 0.0002). Together, our results suggest that chemosensitivity to drug combinations can be predicted based on molecular models and provide the framework for evaluation of such models in patients undergoing combination chemotherapy for cancer. If validated in vivo, such predictive models have the potential to guide therapeutic choice at the level of an individual's tumor. [Mol Cancer Ther 2007;6(2):57886]
Introduction
Most patients with epithelial cancers requiring systemic treatment undergo combination chemotherapy. However, a major challenge in these patients has been the prediction of chemotherapeutic efficacy based on the biological profile of the tumor. Recent studies have indicated that chemotherapeutic sensitivity of the NCI-60 cancer cell line panel can be, in part, predicted by gene expression (1). However, the predictive utility of these molecular correlates of chemosensitivity has not yet been validated for other cellular histologies or for drug combinations. Both issues have been partially addressed in bladder cancer where tumor expression profiling has been used to predict responsiveness to one combination chemotherapy (2). Unfortunately, the practical consequences of this result are limited because this combination is infrequently used, and similar studies are laking for more popular combinations. Furthermore, with the constant arrival of new agents, there is a steady evolution of combination chemotherapy regimens. Together, these two facts complicate the selection process of the best combination therapy for individual patients.
Here, we address this limitation by developing a novel approach that predicts chemotherapeutic responsiveness of human bladder cancer cells to combinations of gemcitabine with cisplatin (3), cisplatin with paclitaxel (4), and gemcitabine with paclitaxel (5) that have shown promising results in clinical trials. Furthermore, by virtue of using single-drug sensitivities to mathematically predict combination effects (rather than using effects of combination directly), this approach has the unique advantage of allowing the evaluation of any number of agents in combination and of allowing the integration of new agents into new combinations as needed.
Materials and Methods
Cell Lines, Cell Culture, Gene Expression Profiling, and Dose-Response Data Generation and Analyses
The human bladder cancer cell lines and the respective growth conditions used in this study have been previously described (6, 7). Cisplatin was purchased from Sigma (St. Louis, MO), dissolved in Dulbecco's PBS, and aliquoted in 1 mg/mL stocks. Paclitaxel was purchased from Sigma, dissolved in DMSO, and aliquoted in 1 mmol/L stocks. Gemcitabine was purchased from the University of Virginia Medical Center Pharmacy, dissolved in PBS, and aliquoted in 0.1 mol/L stocks. Cell lines were maintained in appropriate media, in a humidified atmosphere containing 5% CO2 in air, except CRL2169 (SW780), which requires no CO2 for its growth. Cell lines were subcultured in an aqueous solution of 0.05% trypsin (GibcoBRL, Grand Island, NY; 1:250) and 0.016% EDTA. Each cell line was used within 10 passages from its archival passage number to minimize any long-term cell culture effects. Gene expression analysis of bladder cell lines was carried out as previously described using the HG-U133A GeneChip array (Affymetrix, Santa Clara, CA; refs. 6, 7). The image file was analyzed with robust microarray analysis to obtain the expression intensity values of the microarray data (8).
Cell lines were seeded in 96-well cell culture plates (Costar, Corning, NY) at a density of 1,000 per well. Twenty-four hours later, cells were exposed to the drugs diluted in RPMI 1640, containing 10% fetal bovine serum, a concentration that is required by >75% of cell lines for their normal growth, at a total volume of 200 µL. Each drug dose was plated in triplicate, and the experiment was repeated four to seven times. The doses for cisplatin were 200, 400, 800, 1,600, 3,200, and 6,400 ng/mL; for paclitaxel were 0.0001, 0.001, 0.002, 0.005, 0.01, and 0.1 µmol/L; for gemcitabine were 0.001, 0.01, 0.1, 1, 10, and 100 µmol/L. Plates were incubated for 72 h with carrier or drug, and growth inhibition was assessed by Alamar Blue (BioSource International, Inc., Camarillo, CA; refs. 9, 10). Our doses for cisplatin, paclitaxel, and gemcitabine were chosen to be similar to the range of doses used by the National Cancer Institute (NCI) in their screening of the NCI-60 set of cell lines.6
Estimation of GI50 Values
From the dose-response data, log10(GI50) values (log base 10 of concentration required to inhibit cell growth by 50% in comparison with untreated control) were estimated for all the cell lines by deriving log(dose) concentration curves on cell count percents as described below. To estimate the GI50 values reliably, we computed Euclidean distances among all replicated experiments and excluded outlying experiments if they were in the top 20% among all measured distances. This percentage was determined heuristically based on the general observations in experimental quality control. Furthermore, we did not see significant changes in our results by slightly changing this proportion as several replicated experiments were averaged to estimate our GI50 values (data not shown). Subsequently, the data were fitted to the following nonlinear regression model for estimating dose-response curve of each cell line:
![]() |
and ß determine the shape of a fitted line. This sigmoidal regression function was used to capture the natural shapes of drug dose responses. Thus, the estimated ß is the predicted log10(GI50) value, the expected log concentration achieving the cell count reduction of 50%. Similarly, log10(GI30) and log10(GI70) values (i.e., the concentrations required to inhibit cell growth by 30% and 70%, respectively, in comparison with untreated control) were also calculated.
Determination of Sensitive and Resistant Cell Lines for Single Drug Sensitivity
Cell line drug sensitivity was classified using the GI estimates and application of a criterion dose concept. We defined the criterion dose as the minimum log10(drug dose) among the experimental dose concentrations of each compound at which at least 25% of the cell lines showed growth inhibition >50%. Criterion doses were determined as log10(400 ng/mL) for cisplatin, log10(0.005 µmol/L) for paclitaxel, and log10(0.1 µmol/L) for gemcitabine, which provided at least 10 drug-"sensitive" cell lines for each drug. Using these criterion dose concentrations, each cell line was defined as sensitive if log10(GI50)
criterion dose; strongly sensitive if log10(GI30) < criterion dose, or resistant if log10(GI70) > criterion dose; and intermediate if log10(GI50) > criterion dose and log10(GI70) < criterion dose.
Statistical Discovery of Molecular Chemosensitivity Prediction Models for Single Drugs
For statistical discovery of prediction models, all 22,215 genes on the HG-U133A array were first evaluated for their ability to differentiate sensitive and resistant cell lines; intermediate lines were excluded from the analysis. The most significant genes were selected both by local pooled error test (11) and significance analysis of microarray method (12). After candidate biomarker probes were identified for each tested compound on the basis of significant differential expression for drug sensitivity, we next searched among those candidate biomarkers for ones that would form optimal parsimonious models for prediction of the activity of the compound. For this, we used the "misclassification-penalized posterior" (MiPP) algorithm, which we introduced previously and is available at the open-source Bioconductor web site7 (13). MiPP is based on stepwise incremental classification modeling discovery for the optimal, most parsimonious prediction models and double cross-validated evaluation for each trained prediction model. Model training can be done from several different classification modeling techniques such as linear discriminant analysis, quadratic discriminant analysis, support vector machines, or logistic regression; linear discriminant analysis was used for most application in our current study. In the double cross-validation, the first cross-validation is based on random splitting of the whole data set into a training set and an independent test set for external model validation. The second is an n-fold cross-validation on the training set to avoid the pitfalls of a large-screening search and to obtain the most parsimonious optimal prediction models. Independent splits of the data result in multiple prediction models. MiPP generates multiple independent splits, which, in turn, results in multiple prediction models. The multiple models from different splits were re-evaluated on a large number of (e.g., 100) random splits of test and training sets to obtain their objective confidence bounds with the summary index, the so-called standardized MiPP (sMiPP score), which varies between 1 to 1, from the worst to the best. From this confidence interval evaluation, mean and lower 5% sMiPP scores were obtained for each of the candidate prediction models, together with mean misclassification rates (error rate). The final prediction of sensitive (or resistant) cell lines was done by averaging its (posterior) classification probabilities of the top three prediction models exceeding 5% sMiPP > 0.5. In performing MiPP analysis, we used the default values for many tuning variables of the MiPP Bioconductor R package. For example, n.fold, p.test, n.split, and n.seq were 5, 1/3, 20, and 3, respectively. However, we preselected the most significant top 1% genes by local pooled error and significance analysis of microarray and did not use the MiPP gene selection option by setting percent.cut = 0.
Statistical Chemosensitivity Prediction for Combination Drug Treatments
Prediction of combination drug efficacy was obtained based on the final single-drug prediction models, directly using classification probabilities of each cell line from these models. That is, assuming two different drug compounds acted independently, the combination chemosensitivity probability PAB of their combination treatment was derived as:
![]() |
0.75) for predicting chemosensitivity of each cell line to combination treatment. Results
Evaluation of In vitro Drug Sensitivity of Human Bladder Cell Lines to Single Agents
To approach the development of molecular models of chemotherapeutic sensitivity in human bladder cancer, we focused on a well-defined series of 40 urothelial cell lines for which we could measure sensitivity to relevant chemotherapeutic agents in vitro and correlate these responses with global measurements of gene expression (Fig. 1
). The in vitro sensitivity of these 40 bladder cancer cell lines to cisplatin, paclitaxel, and gemcitabine was carried out as described in Materials and Methods. Typical dose-response curves for representative sensitive and resistant cell lines are shown for each agent in Fig. 2A
. The cell lines were then divided into three groups (sensitive, intermediate, and resistant) based on GI estimates and the criterion dose (defined in Materials and Methods). Figure 2B to D shows the log10(GI30), log10(GI50), and log10(GI70) of the 40 cell lines for each of the agents. For cisplatin, we identified 16 sensitive and 11 resistant cell lines (Fig. 2B), 17 sensitive and 11 resistant cell lines for paclitaxel (Fig. 2C), and 8 sensitive and 11 resistant for gemcitabine (Fig. 2D). Cell lines that did not meet the "sensitive/resistant" criteria were excluded from further analyses. For some cell lines, log(GI) values could not be estimated due to flat response curves in nonlinear regression model fitting; thus, these log(GI) values of cell lines were thresholded at the maximum dose concentration and were classified as resistant.
|
|
|
|
We evaluated the performance of these in silico predictions by randomly selecting 15 of the 40 bladder carcinoma cell lines, attempting to roughly balance the numbers of predicted sensitive and resistant cell lines across the three drug combinations. We used the single-drug criteria dose and exposed cells to both drugs simultaneously. The growth of cell lines exposed to the drug combinations compared with control (no drug) was expected to be <55% for sensitive and >55% for resistant cell lines at these doses.
Overall, 35 of the 45 predictions were correct (binomial test, P = 0.0002; Table 2 ; Fig. 3D). Twelve of 15 cell lines (80%; binomial test, P = 0.03) were predicted correctly for the cisplatin-paclitaxel combination. Of the three misclassified cell lines, one sensitive line was predicted as resistant, and two resistant cell lines were predicted as sensitive. For the cisplatin-gemcitabine combination, 12 of 15 lines were also predicted correctly; three sensitive cell lines were incorrectly predicted as resistant (80% accuracy; binomial test, P = 0.03). Finally, for the combination of paclitaxel and gemcitabine, 11 of 15 lines were correctly classified; three sensitive and one resistant cell lines were misclassified as resistant and sensitive, respectively (73% accuracy; binomial test, P = 0.11).
|
Discussion
Here, we combined a novel mathematical approach (MiPP probabilities) with comprehensive gene expression profiles of 40 urothelial cell lines to discover high-performance molecular prediction models for single and combination chemotherapeutic sensitivity. The high-performance characteristics of the predictive models obtained in this study may be due to several factors. First, we used a panel of cancer cell lines derived from only one histologic type: urothelial cancer. In contrast to the NCI-60 cancer cell panel, which is comprised of cell lines from multiple anatomic origins, a single anatomic origin should eliminate confounding and biased gene expression signals that represent tissue-dependent sensitivity to different chemotherapy agents. Furthermore, the majority of the cell lines used in this study (14) are derived from invasive or metastatic human urothelial tumors that represent the typical patient population that would receive systemic chemotherapy. Hence, we anticipate that these prediction models may be applicable to clinical urothelial cancer. This conclusion is supported by the observation that cisplatin, a drug used in current clinical treatment of urothelial cancer (15), was highly effective in our assay (i.e., 16 of 40 cell lines meeting the chemosensitive criterion).
To identify gene prediction models for chemosensitivity, we used the MiPP method (13). Several studies have shown good predictive classification of cancer subtypes and prognosis using methods that require large numbers of (>50) genes (1618), whereas models that are dependent on only a small number of predictive genes has been limited despite the obvious practical advantages. The MiPP method combines the best of both approaches by maintaining excellent predictive accuracy with a small set of genes that are easy to evaluate in human tumors using currently available techniques such as real-time reverse transcription-PCR. This feature is a significant advantage as we begin to prospectively evaluate these genes for their ability to predict tumor response in patients treated with drug combinations.
The approach taken here led to the identification of predictive gene models for each of the three drugs. Cisplatin model 1 is comprised of TGM2, MOAP1, HIST2H2AA, and MRPS30; model 2 contains CAV2, LCP1, and MOAP1; and model 3 includes CCNG2, PEG10, and WNT5B. By examining the function of the genes encompassed by these models, a common functional theme was noted, that is, their direct (TGM2, MOAP1, and CAV2) or indirect (HIST2H2AA and LCP1) participation in apoptosis. Modulator of apoptosis 1 (MOAP1) is an important component of the pathway that links death receptors and the apoptotic machinery (19). Caveolin 2 (CAV2) is a major component of the inner surface of caveolae and is implicated in the control of cellular growth, signal transduction, lipid metabolism, and apoptosis (20). Lymphocyte cytosolic protein 1 (LCP1) is found in hemopoietic cell lineages and also in many types of malignant human cells of non-hemopoietic origin. Cyclin G2 (CCNG2) is a member of the cyclin family. Northern blot analysis revealed that cyclin G2 mRNA fluctuates throughout the cell cycle with peak expression in late S phase. Furthermore, cyclin G2 is induced by the DNA-damaging agent actinomycin D (21).
Models for paclitaxel included several genes involved in essential eukaryotic cell functions, such as protein modification (PLAT), spermatogenesis and cell differentiation (DZIP1), and negative autocrine growth factor regulation (LGALS1). However, perhaps the most interesting of this group is KIF14. This gene is responsible for microtubule motor activity and is expressed at very low levels in normal tissue samples, compared with significantly increased expression in the majority of tumor samples (22). Its overexpression may lead to rapid mitoses, potentially leading to aneuploidy. KIF14 overexpression is most striking in retinoblastoma, lung, breast, thymus, and tumors and associated with decreased survival in lung cancer (22). This relationship to paclitaxel sensitivity is intriguing because this drug promotes the assembly of microtubules from tubulin dimers and stabilizes microtubules by preventing depolymerization, thus inducing abnormal arrays of microtubules throughout the cell cycle. We are in the process of investigating this further.
Despite the robust and impressive predictive power of this approach and the lack of viable alternatives, our technique does have limitations. The most significant is its inability to predict drug synergy. For example, a mechanistic basis for gemcitabine/cisplatin (23), paclitaxel/cisplatin (24), and paclitaxel/gemcitabine (25) cytotoxic synergies have been described. Our model was less effective at predicting positive synergy (37%) than negative synergy or antagonism (12%). Fortunately, this differential performance is clinically useful because it implies that drug combinations predicted to be effective in a particular patient using our models are likely to be so and thus can be used as first line treatment. Conversely, those combinations that are predicted to be ineffective may have a low chance of being efficacious based on synergy and thus can be used as "salvage" therapy. Furthermore, it is possible that our approach could also be used to develop predictive models for drug synergy. However, this would require the generation of extensive in vitro drug response data for each of the different combinations. However, given the biological and clinical relevance of this issue, efforts in this direction are clearly warranted.
Our use of MiPP prediction model was mainly based on the linear discriminate analysis, one of its classification method options in this study. However, other classification options, such as quadratic discriminant analysis, support vector machines, and logistic regression methods, may be further explored in a future study. In addition, our statistical prediction approach to the combination treatments was relatively simple based on the independence assumption of any two compounds. Because our goal is to predict combination responses by two independent experiments of two or multiple compounds, this assumption is statistically inevitable. However, if such drug experiments are done on the same panel of cancer cell lines, we may be able to optimize such a prediction by considering their drug sensitivity through their gene models and then searching for the combination treatments with the highest combined response rates.
In summary, we have developed and validated a novel molecular chemosensitivity prediction model for commonly used combinations of cisplatin, paclitaxel, and gemcitabine, using only the results of their individual drug responses. We believe this prediction strategy warrants prospective validation in the clinical setting and, given the parsimonious nature of the predictions shown here, should be straightforward to implement.
Footnotes
Grant support: NIH grant CA075115.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received 8/15/06; revised 10/28/06; accepted 12/18/06.
References
This article has been cited by other articles:
![]() |
P. C. Black, G. A. Brown, T. Inamoto, M. Shrader, A. Arora, A. O. Siefker-Radtke, L. Adam, D. Theodorescu, X. Wu, M. F. Munsell, et al. Sensitivity to Epidermal Growth Factor Receptor Inhibitor Requires E-Cadherin Expression in Urothelial Carcinoma Cells Clin. Cancer Res., March 1, 2008; 14(5): 1478 - 1486. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |