Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Proteomics Clin Appl. Author manuscript; available in PMC 2010 November 22.
Published in final edited form as:
PMCID: PMC2989613

Quantitative Proteomic Analysis of Ovarian Cancer Cells Identified Mitochondrial Proteins Associated with Paclitaxel Resistance


Paclitaxel has been widely used as an anti-mitotic agent in chemotherapy for a variety of cancers and adds substantial efficacy as the first-line chemotherapeutic regimen for ovarian cancers. However, the frequent occurrence of paclitaxel resistance limits its function in long-term management. Despite abundant clinical and cellular demonstration of paclitaxel resistant tumors, the molecular mechanisms leading to paclitaxel resistance are poorly understood. Using genomic approaches, we have previously identified an association between a BTB/POZ gene, Nac1, and paclitaxel resistance in ovarian cancer. The experiments presented here have applied multiple quantitative proteomic methods to identify protein changes associated with paclitaxel resistance and Nac1 function. The SKOV-3 ovarian serous carcinoma cell line, which has inducible expression of dominant negative Nac1, was used to determine the paclitaxel treatment associated changes in the presence and absence of functional Nac1. Quantitative proteomic analyses were performed using iTRAQ labeling and mass spectrometry. Two label-free quantitative proteomic methods: LC-MS and spectral count were used to increase confidence of proteomic quantification. A total of 1371 proteins were quantified by at least one of the quantitative proteomic methods. Candidate proteins related to paclitaxel and NAC1 function were identified in this study. Go analysis of the protein changes identified upon paclitaxel resistance revealed that cell component enrichment related to mitochondria. Moreover, tubulin and mitochondrial proteins were the major cellular components with changes associated with paclitaxel treatment. This suggests that mitochondria may play a role in paclitaxel resistance.

Keywords: ovarian cancer, paclitaxel, Taxol, mass spectrometry, proteomics

1. Introduction

Paclitaxel (Taxol®) is a potent antimitotic agent which is currently employed for the treatment of many human cancers and as an inflammation deterrent in drug-eluting cardiovascular stents [1]. Paclitaxel is known to induce cytotoxicity by preventing tubulin depolymerization during the metaphase to anaphase transition of mitosis [2, 3] or by triggering apoptosis through regulating the expression of apoptosis-related proteins in both the caspase-dependent and caspase-independent pathways [1]. Unfortunately, while paclitaxel causes initial remission of ovarian cancer, the tumor often acquires resistance and recurs [4]. The molecular mechanisms underlying paclitaxel resistance remain unclear. The experiments presented here attempt to identify the proteins associated with paclitaxel resistance in ovarian cancer cells in order to facilitate the elucidation of molecular mechanisms of paclitaxel induced apoptosis and acquired resistance and discovery of potential drug targets for ovarian cancer with paclitaxel resistance.

Nac1, a member of BTB/POZ gene family, is a transcription repressor that is essential for the growth and survival of tumor cells [5]. We have previously associated Nac1 overexpression with tumor recurrence and paclitaxel resistance in ovarian serous carcinoma [5, 6]. However, the function of Nac1 for paclitaxel resistance is not well understood. To explore this function, we generated the SKOV-3 N130 cell line which is an ovarian serous carcinoma cell line (SKOV-3) with stable transfection of N130/EGFP controlled by tTA (tetracycline-controlled transactivator) [5]. This Tet-OFF inducible system can trigger the expression of N130 by removal of doxycycline which inhibits the function of Nac1.

The relationships between paclitaxel resistance and Nac1 were explored here using iTRAQ (isoabaric tags for relative and absolute quantitation) quantitation method, and also measured by label-free quantitation methods, LC-MS and spectral count. The three methods are among the several high throughput quantitative proteomic methods that have been developed in the past decade. Method development in this field has occurred in two directions: label-dependent and label-free. The label-dependent methods are widely used and include derivitizing methods [7] such as isotope-coded affinity tags (ICAT) [8] and isobaric tags for relative and absolute quantitation (iTRAQ) [9] and non-derivitizing methods such as stable isotope labeling with amino acids in cell culture (SILAC) [10] and 18O labeling [11]. When applied to proteomics, stable isotope labeling allows for the accurate measurement of the relative peptide abundance by direct comparison of light and heavy peptides in the same spectrum.

While label-dependent methods comprise the gold standard for quantitative techniques, only a limited number of samples can be quantified using isotopic derivitization in a single experiment due to the fixed number of channels from the labeling reagents. Thus an alternative direction of method development for quantitative proteomics, that of label-free quantitation methods, has evolved and includes the liquid chromatography-mass spectrometry (LC-MS) method [12, 13] and spectral count [14]. The LC-MS method determines the peptide abundance by comparing the intensity of the same peptide peak in multiple LC-MS runs. Quantitation of protein abundance by spectral count is based on the number of redundant spectra acquired for each protein from different samples in the LC-MS/MS analyses. Label-free quantitative methods are theoretically capable of quantifying an unlimited number of samples in a single study. The limitation of the label-free quantitative method is that the quantitation accuracy relies heavily on the reproducible analyses of different samples in multiple LC-MS and LC-MS/MS analyses [1518].

Therefore, it is clear that a high-throughput quantitative proteomic method has the potential to identify a large number of protein changes but that these measurements must be validated. The validation of these protein changes with traditional methods such as Western blots or immunohistochemistry is limited due to the availability, expense of the antibodies, and the low throughput of the assays. Therefore, the proteomics study should give the more confident ones for the further immune based validation. The study presented here identified protein changes related to paclitaxel resistance and Nac1 function using most reliable quantitation, iTRAQ labeling. In addition, the two label-free quantitative proteomic methods were used to increase the quantification confidence.

Using the iTRAQ quantitation, most of these changed proteins related to paclitaxel treatment were significantly overrepresented in mitochondria. Our results suggest a new role of mitochondria of ovarian cancer cells in paclitaxel resistance and define potential new targets for treatment of paclitaxel-resistant ovarian cancer. Most protein changes in mitochondria were also identified as up-regulated by the two label-free methods. In addition, we also list the candidate proteins related to NAC1 function.

2 Materials and methods

2.1 Materials

Sequencing grade trypsin was from Promega (Madison, WI); C18 Sep-Pak Vac columns were from Waters (Milford, MA); α-cyano-4-hydroxycinnasmic acid (CHCA) was from Agilent (Palo Alto, CA); iTRAQ reagent and mass calibration standards were from Applied Biosystems (Foster City, CA); BCA assay kit was from Pierce (Rockford, IL); SCX columns and C18 resin were from Sepax (Newark, DE).; Other chemicals were purchased from Sigma-Aldrich (St. Louis, MO).

2.2 Treatments of ovarian cancer cells

The N130-inducible SKOV-3 ovarian cancer cell line, stably expressed the inducible construct of N-terminal 130 amino acids (N130) of Nac1 gene and enhanced green fluorescent protein (EGFP), was generated and reported previously [5, 6]. N130-inducible SKOV-3 cells were cultured in G400D2 medium (RPMI medium contained 10% FBS, 1% penicillin/streptomycin, 400 μg/ml geneticin, and 2 μg/ml doxycycline). To induce the expression of N130, the cells were washed with PBS twice and cultured in G400 medium (G400D2 medium without doxycycline). Expression of N130-EGFP was confirmed by fluorescent microscope after 28-hour culture in G400 medium.

The N130-inducible SKOV-3 ovarian adenocarcinoma cells with and without N130 expression by culturing in two different medium as described above, were un-treated or treated with 20nM paclitaxel for 72-hour. The cells that remained alive after paclitaxel treatment were harvested at the end of paclitaxel treatment.

2.3 Peptide extraction

The cell pallets were collected and sonicated. Protein concentration was measured by BCA assay. The same amounts of proteins (1 mg) from each condition were denatured in 8M urea in 0.4M NH4HCO3, 0.1% (w/v) SDS solution (pH8.3), and 10mM TCEP (Tris (2-carboxythyl) phosphine) by incubation at 60°C for 1 hour. Proteins were alkylated with 16mM iodoacetamide by incubation at room temperature in the dark for 30 min. The sample was diluted 4-fold by trypsin digestion buffer (100mM NH4HCO3, pH8.3). Trypsin was added at a 1 to 50 part sample protein excess and allowed to digest at 37°C overnight. SDS-PAGE and silver staining was employed to ensure the completion of tryptic digest. The peptides were purified with C18 Sep-Pak Vac columns and resuspended in water with a final concentration of 10μg/μl.

2.4 iTRAQ labeling

Tryptic peptides (50 μg) from each sample were mixed with 20 μl of dissolution buffer provided with iTRAQ kit. The iTRAQ 4-plex reagents were dissolved in 70μl of methanol respectively and strongly vortexed. Each iTRAQ labeling reagent was then added to the sample and mixed. The mixture was incubated at room temperature for 1 hour followed by cleaning up by SCX column.

2.5 Mass spectrometry analysis

For protein quantification by spectral count, each peptide mixture was analyzed twice by the LTQ ion trap mass spectrometer (Thermo Finnigan, San Jose, CA). For protein identification and quantitative analysis using LC-MS, an ESI-QSTAR mass spectrometer (Applied Biosystems, Foster City, CA) was used. In both systems, 2 μl (2 μg) peptides were injected into a peptide cartridge packed with C18 resin, and then passed through a 10 cm × 75 μm i.d. microcapillary HPLC (μLC) column packed with C18 resin. The effluent from the μLC column entered an electrospray ionization source in which peptides were ionized and passed directly into the mass spectrometers. A linear gradient of acetonitrile from 5%–32% over 100 min at flow rate of ~300 nL/min was applied. During the LC-MS mode, data was acquired in the m/z range of 400 and 2000. The MS/MS was also turned on to collect CID using data dependent mode. Each sample was analyzed three times by QSTAR to increase the accuracy of quantification.

iTRAQ labeled peptide was analyzed by both QSTAR and 2-D LC (Nano, Eksigent, Dublin, CA) MALDI TOF/TOF (ABI 4800, Applied Biosystems, Foster City, CA). The analysis on the QSTAR was performed in the same setting as described above. For the analysis by 2-D Nano LC and MALDI 4800-TOF/TOF, on-line integration of 15-cm-long 300μm strong cation exchange column (SCX) with 15-cm-long 300 μm of C18-reverse phase liquid chromatography (RPLC) was employed. Four SCX fractions of 0, 5, 50 and 500mM KCl and 3–45% linear acetonitrile gradient (containing 0.1% TFA and acetonitrile) of RPLC for each fraction were applied before analysis by MALDI-TOF/TOF. Peptides eluted from columns were directly mixed with CHCA and spotted on a MALDI target plate with 768 spots followed by analysis with MS and MS/MS using the ABI 4800 MALDI-TOF/TOF.

2.6 Peptide identifications

The iTRAQ data analyzed either by QSTAR or MALDI and label-free data from LTQ were searched by ProteinPilot software 2.0 [19] against the human International Protein Index database (IPI, version 2.28) using the cut-off probability score of 0.9.

Tandem MS spectra of label-free peptides from the QSTAR were searched with SEQUEST [20] against the same human IPI protein database (version 2.28). The peptide mass tolerance is 2.0Da. Other parameters of database searching are modified as following: cysteine modification (add cysteine 57) and oxidized methionine (add methionine with 16 Da). The output files were evaluated by INTERACT [21] and ProteinProphet [22]. The cutoff of ProteinProphet analysis is the probability score ≥ 0.9 so that low probability protein identifications can be filtered out. For each identified peptide, peptide sequence, protein name, precursor m/z value, peptide mass, charge state, retention time where the MS/MS was acquired, and probability of the peptide identification being correct were recorded and outputted using INTERACT [21].

2.7 Quantitative proteomic analyses

The ratio of the four channels of iTRAQ labeling was determined by the ProteinPilot software.

A suite of software tools of SpecArray were used to analyze the LC-MS data as described previously [23]. For each peptide peak, an abundance ratio of matched peptides in different samples was determined for each peptide peak. An in-house Perl script was then used to link the peptide identification from MS/MS spectra to their corresponding MS peaks by matching precursor mass within 1 Da, retention time within 10 min, and charge state of the peptides.

The identified peptides from LTQ with a probability score ≥ 0.9 were used for the spectral count. To determine the number of MS/MS spectra used for identification of each protein in different conditions using our in-house developed software tool. For peptide sequence that could come from multiple proteins, the spectral count is equally distributed to all proteins with the identified peptide. Due to random sampling of mass spectrometer in collecting MS/MS spectra used for spectral count, we only quantified proteins with at least 4 spectral counts in total from the four cell states.

2.8 Evaluation of the cut-off of protein abundance ratio for proteins changes

To correct for any systematic errors of protein ratio introduced by sample handling and to determine the appropriate cut-off for protein changes, the distribution of abundance ratios in different cell states was generated for each quantitative method. Since the majority of proteins were not expressed differently in two cell states, we normalized the ratio based on the distribution of the protein abundance ratios from two cell states. Proteins fell out of the normal distribution from the abundance ratio of two cell states were considered as altered proteins. The threshold to select protein changes was based on the ratio distribution of two cell states. The mean and standard deviation of ratio from two cell states were calculated, and the abundance of proteins with an abundance ratio outside of one standard deviation from the mean were flagged as altered.

2.9 Cellular component classification of changed genes

To classify the changed proteins into cellular component, GO (Gene Ontology) [24] analysis ( was performed. All the identified and quantified proteins by iTRAQ quantitation were used as background. Protein changes due to paclitaxel quantified by iTRAQ were used as changed proteins. P value was calculated using one-side Fisher exact test. To correct for multiple testing errors, p value was adjusted the minimum P method of Westfall and Young[25].

3 Results and discussion

3.1 Inducible expression of N130 in SKOV-3 cells

To identify proteins related to paclitaxel treatment and resistance, SKOV-3 cells with inducible expression of N130 protein [6, 26] was used in this study. The quantitative proteomic analyses of paclitaxel treatment for SKOV-3 N130 cells with and without expression of N130 are schematically illustrated in Figure 1 and consist of four steps: 1) two dishes were treated to induce the expression of N130 and two dishes were untreated as controls; 2) One dish with expression of N130 and one dish without expression of N130 were treated with paclitaxel, and the other two dishes were not treated with paclitaxel act as controls; 3) the peptides were extracted from cell lysate of the four cell states by sonication and trypsin digestion; 4) the tryptic peptides were identified and quantified by iTRAQ labeling, LC-MS, and spectral count.

Figure 1
Flowchart of the quantitative study of Nac1 function and paclitaxel treatment

To evaluate the expression of N130, the florescence of EGFP was monitored as an indicator to determine whether N130 was induced after removing doxycycline from culture medium. SKOV-3 N130 cells were observed after 28-hour culture in medium with and without doxycycline according to our previous study [6, 26]. The cells were observed under florescence microscope (Figure 2A and B). The induced expression of N130/EGFP by doxycycline withdrawal was indicated by the green fluorescent (Figure 2D), which was not observed in cells cultured with doxycycline (shown in Figure 2C). Thus the expression of N130 can be robustly induced in SKOV-3 N130 cells.

Figure 2
Expression of fusion protein N130-EGFP in SKOV-3 ovarian cell line: A) SKOV-3 N130 cultured with doxycycline, the N130-EGFP expression was off; B) SKOV-3 N130 cultured without doxycycline, the N130-EGFP expression was on; C) no fluorescent came from SKOV-3 ...

3.2 Quantitative proteomic analyses to identify protein changes

To determine the protein changes related to paclitaxel treatment and Nac1 function, the ovarian cancer cells with and without Nac1 function were treated with paclitaxel. After treated with 20nM paclitaxel for 72-hour, around 60 % cells was alive, which were considered as cells resistant to paclitaxel. The cells that remained alive after paclitaxel treatment were harvested at the end of paclitaxel treatment. The cell pallets were sonicated and 1 mg proteins from each cell states were digested by trypsin, followed by quantitative proteomic analysis using iTRAQ. A portion of tryptic peptides (50 μg) from N130-ON without paclitaxel treatment (ON−T), N130-ON with paclitaxel treatment (ON+T), N130-OFF without paclitaxel treatment (OFF−T) and N130-OFF with paclitaxel treatment (OFF+T) cells were labeled with 114, 115, 116 and 117 of iTRAQ reagents and analyzed by LC-MALDI TOF/TOF and LC-QSTAR for quantitative proteomic analysis. We were able to identify and quantify 850 proteins using iTRAQ labeling.

We then determined the protein changes in two cell states quantified by iTRAQ. When cells are induced of N130 expression or treated of paclitaxel, majority of cellular proteins are expected to be not affected and stay in the same level [27]. However, due to errors introduced by analytical procedures, such as sample handling and quantification process, the protein ratio from majority of proteins may be shifted. To determine the proteins with abundance changes in two cell states, histogram was used to generate the number of proteins in different abundance ratio (Figure 3). The threshold was set as < 0.7 and >1.3 for iTRAQ labeling. Majority of proteins (554 proteins, 65%) were distributed within one standard deviation (0.25) from the mean (1.046) and were considered as unchanged. Proteins that fell out of one standard deviation of the normal distribution curve were considered as with changed. A total of 296 proteins were changed due to paclitaxel treatment (160 proteins) or Nac1 inactivation (93 proteins) or both of paclitaxel treatment and Nac1 inactivation (181 proteins).

Figure 3
The histogram analysis of peptide ratios of different cell states quantified by iTRAQ

Nac1 was determined as changed upon the induced expression of N130. A total of seven peptides were identified and quantified from Nac1, and all the identified peptides were located in the N-terminus 1-130 amino acids (with 71% sequence coverage of N 130), indicating that the identified peptides were likely from overexpressed N130 instead of endogenous Nac1 protein. The amount of N130 in N130 ON cells was measured as about 10-fold higher than it in N130 OFF cells (Table 1). The quantitative results of N130 expression confirmed that: 1) the inducible Tet-OFF system was efficient in inducing N130 expression; 2) the iTRAQ quantitative methods was able to determine the relative abundance of proteins and could be used for identification of other protein changes.

Table 1
Overexpression of N130 determined by three quantitative proteomic methods. For iTRAQ quantitation, all ratios are normalized to the reporter ion at m/z 116 (OFF-T). For LC-MS quantitation, the number showed the peak intensities.

3.3 Altered proteins related to Nac1 function and paclitaxel treatment

The 296 unique protein alterations determined by iTRAQ labeling were listed in Table 2 and grouped into 3 classes: 1) Protein changes upon paclitaxel treatment (OFF−T vs. OFF+T). These include proteins elevated upon paclitaxel treatment such as tubulin beta-5 chain, tubulin alpha-4 chain, mitochondrial proteins such as cytochrome c and ATP synthase, mitochondrial inner membrane protein, acute-phase proteins such as hemoglobin, cell surface antigens such as CD44 and 4F2 cell surface antigen, etc., and proteins with decreased abundance upon paclitaxel treatment such as seven subunits of ribosomal proteins, proteins regulating cell meiosis, mitosis and postmitotic functions such as mitogen-activated protein kinase3, etc. 2) Changed protein expression upon inactivation of Nac1 (OFF−T vs. ON−T). Induced expression of N130 (ON−T) inhibits the function of Nac1. Since Nac1 is a potential transcriptional repressor [26], the proteins with altered expression after Nac1 inhibition could be controlled by Nac1. These proteins include Ras-related protein Rab-8, transcription repressor, eukaryotic translation initiation factor 3, etc. 3) Changes of protein abundance upon paclitaxel treatment and induced N130 expression (OFF−T vs. ON+T), and proteins in this class might associate with the function of Nac1 gene in the response to paclitaxel treatment. Proteins in this class include Ras GTPase-activating-like protein IQGAP1, polyadenylate-binding protein 1, etc.

Table 2
Proteins associated with paclitaxel treatment with and without N130 overexpression determined by iTRAQ

Although the proteins identified in this proteomic study need further investigation to facilitate the understanding of the biological mechanism of Nac1 function or paclitaxel treatment, the results provide a list of proteins and cellular machinery, including ribosomal complexes, cell surface antigens, and stress response proteins, such as heat shock proteins and acute-phase proteins. These protein changes associated with Nac1 or paclitaxel resistance can be exploited as targets for treatment of paclitaxel resistance.

To further analyze the relationship of the protein changes upon paclitaxel treatment, the GO categories of protein changes were classified. Cellular components analysis revealed that the protein changes are significantly overrepresented in mitochondrion in the set of all proteins identified by iTRAQ (p value of Fisher’s product: 8.8 e-5, p value corrected by multiple testing: 0.024).

Furthermore, we found some interesting co-regulation of tubulin and mitochondrial proteins after paclitaxel treatment. Tubulin is a well-known target for paclitaxel function and responsible for paclitaxel induced cell death [33]. One of the mechanisms of paclitaxel function is believed to induce cell death by altering microtubule assembly through the binding to the microtubule polymer so as to stabilize microtubules [34], as a result, it disrupts the normal re-assembling of microtubule network which is required by mitosis and cell proliferation [3]. Another protein, cytochrome c, was reported previously of release from mitochondrion thus inducing cell apoptosis upon paclitaxel treatment [35]. However, how cytochrome c was released upon paclitaxel treatment is not clear. Interestingly, in this study, both α-4 and β-5 subunit of tubulin were observed of up-regulated after paclitaxel treatment (Table 1), so were many mitochondrial proteins including mitochondrial trifunctional enzyme, mitochondrial ATP synthase, cytochrome c, Serine hydroymethyltransferase, GrpE protein homolog 1, Mitochondrial inner membrane protein, Complement component 1 Q subcomponent binding protein, Thioredoin-dependent peroide reductase, and mitochondrial malate dehydrogenase, etc.

This observation suggests a regulation of mitochondrial function associated with paclitaxel treatment and tubulins. The regulation of mitochondrial function by tubulins was also reported by several studies recently. The regulation might be the result of direct interaction of the voltage-dependent anion channel (VDAC) on mitochondrial outer membrane with tubulin [3638]. Taken together, a hypothesis is that mitochondria may be involved in the response to paclitaxel treatment. Mitochondrial function is the key player for cell apoptosis, and the mechanism of paclitaxel treatment might be to induce apoptosis through tubulin polymerization and regulation of mitochondrial function.

3.4 protein changes determined by label-free quantitation

Quantitative analysis using different quantitative proteomic methods may increase the confidence of the protein changes if the proteins could be identified and quantified by multiple methods consistently. To this end, label free quantitation methods were also employed in this study. The tryptic peptides from the four cell states without iTRAQ labeling were analyzed three times with the QSTAR for the LC-MS quantitative analysis and two times with the LTQ for spectral count (Figure 1). A total of 383 proteins were quantified by the LC-MS method, and 757 were quantified by spectral count (Figure 4).

Figure 4
Venn diagram depicts the number of proteins quantified by each quantitation method

We then determined the proteins that were changed in two cell states quantified by LC-MS and spectral count. Similar as iTRAQ quantitation, proteins that fell out of one standard deviation of the normal distribution curve were considered as with changed expression. The thresholds were determined as <0.75 and >1.15 for both spectral count and LC-MS.

Since N130 was induced expressed in cells, N130 was the perfect internal control for quantitation. N130 should be over-expressed in N130 ON cells compared to the N130 OFF cells. All the three quantitation showed the higher abundance of N130 in N130 ON cells (Table 1). However, the detection limitation varied among the three methods. The N130 overexpression were undetectable in N130 OFF cells for LC-MS and spectral count (Table 1), which may come from different instrumentations with various dynamic range and background. Further work and experiments will help define the most accurate quantifications along with better standards for calibrating the ratio of protein abundance from different methods. This is needed in proteomics to improve quantitative accuracy [28]. Nevertheless, the quantitative results of N130 expression confirmed that the three quantitative proteomic methods could be used to increase the confidence of quantitation.

For the mitochondria protein changes upon piclitaxel treatment, 7 out of 14 proteins determined by iTRAQ were also measured by the label free methods in the same track, e.g. ATP synthase, cytochrome c, Trifunctional enzyme, and Enoyl-CoA hydratase, etc (Table 3). Those proteins consistently determined by the two label-free methods confirmed the real changes of mitochondrial proteins upon piclitaxel treatment. This study represents the first proteomic study to discover the association of paclitaxel treatment and mitochondria protein changes in ovarian cancer cells, which may offer a new direction for studying the mechanism of drug resistance of cancer cells.

Table 3
Mitochondrial protein changes related to paclitaxel quantified by iTRAQ were also measured by label-free quantitation methods

4 Concluding remarks

In this study, 1371 proteins were identified and quantified from Nac1 dominant negative model, SKOV-3 N130 cell line, associated with paclitaxel resistance and Nac1 function using iTRAQ quantitation, LC-MS method and spectral count. Candidate proteins related to paclitaxel resistance and NAC1 function were determined. Go analysis of the protein changes upon paclitaxel resistance revealed that protein changes significantly overrepresented in mitochondria. The co-regulation of tubulins and mitochondrial proteins was found, which suggests the roles of mitochondria in response to paclitaxel treatment. The identified proteins will be useful for further study of biological functions of Nac1 and elucidation of the molecular mechanism of paclitaxel treatment and resistance.

Supplementary Material





This work was supported by federal funds from the National Cancer Institute, National Institutes of Health, by grant R21-CA-114852 and RO1-CA-103937 (IMS) and Early Detection and Research Network (EDRN). We gratefully acknowledge the support from the Mass Spectrometry Facility at the Johns Hopkins University and the support from Applied Biosystems.


The authors have declared no conflict of interest.


1. Khayat D, Antoine EC, Coeffic D. Taxol in the management of cancers of the breast and the ovary. Cancer Invest. 2000;18:242–260. [PubMed]
2. Jordan MA, Wendell K, Gardiner S, Derry WB, et al. Mitotic block induced in HeLa cells by low concentrations of paclitaxel (Taxol) results in abnormal mitotic exit and apoptotic cell death. Cancer Res. 1996;56:816–825. [PubMed]
3. Amos LA, Lowe J. How Taxol stabilises microtubule structure. Chem Biol. 1999;6:R65–69. [PubMed]
4. Sangrajrang S, Fellous A. Taxol resistance. Chemotherapy. 2000;46:327–334. [PubMed]
5. Ishibashi M, Nakayama K, Yeasmin S, Katagiri A, et al. A BTB/POZ gene, NAC-1, a tumor recurrence-associated gene, as a potential target for Taxol resistance in ovarian cancer. Clin Cancer Res. 2008;14:3149–3155. [PubMed]
6. Nakayama K, Nakayama N, Davidson B, Sheu JJ, et al. A BTB/POZ protein, NAC-1, is related to tumor recurrence and is essential for tumor growth and survival. Proc Natl Acad Sci U S A. 2006;103:18739–18744. [PubMed]
7. Zhang H, Li XJ, Martin DB, Aebersold R. Identification and quantification of N-linked glycoproteins using hydrazide chemistry, stable isotope labeling and mass spectrometry. Nat Biotechnol. 2003;21:660–666. [PubMed]
8. Gygi SP, Rist B, Gerber SA, Turecek F, et al. Quantitative analysis of complex protein mixtures using isotope-coded affinity tags. Nat Biotechnol. 1999;17:994–999. [PubMed]
9. Ross PL, Huang YN, Marchese JN, Williamson B, et al. Multiplexed protein quantitation in Saccharomyces cerevisiae using amine-reactive isobaric tagging reagents. Mol Cell Proteomics. 2004;3:1154–1169. [PubMed]
10. Ong SE, Blagoev B, Kratchmarova I, Kristensen DB, et al. Stable isotope labeling by amino acids in cell culture, SILAC, as a simple and accurate approach to expression proteomics. Mol Cell Proteomics. 2002;1:376–386. [PubMed]
11. Yao X, Freas A, Ramirez J, Demirev PA, Fenselau C. Proteolytic 18O labeling for comparative proteomics: model studies with two serotypes of adenovirus. Anal Chem. 2001;73:2836–2842. [PubMed]
12. Bondarenko PV, Chelius D, Shaler TA. Identification and relative quantitation of protein mixtures by enzymatic digestion followed by capillary reversed-phase liquid chromatography-tandem mass spectrometry. Anal Chem. 2002;74:4741–4749. [PubMed]
13. Chelius D, Bondarenko PV. Quantitative profiling of proteins in complex mixtures using liquid chromatography and mass spectrometry. J Proteome Res. 2002;1:317–323. [PubMed]
14. Liu H, Sadygov RG, Yates JR., 3rd A model for random sampling and estimation of relative protein abundance in shotgun proteomics. Anal Chem. 2004;76:4193–4201. [PubMed]
15. Zhang H, Yi EC, Li XJ, Mallick P, et al. High throughput quantitative analysis of serum proteins using glycopeptide capture and liquid chromatography mass spectrometry. Mol Cell Proteomics. 2005;4:144–155. [PubMed]
16. Olson MT, Blank PS, Sackett DL, Yergey AL. Evaluating reproducibility and similarity of mass and intensity data in complex spectra--applications to tubulin. J Am Soc Mass Spectrom. 2008;19:367–374. [PubMed]
17. Frewen BE, Merrihew GE, Wu CC, Noble WS, MacCoss MJ. Analysis of peptide MS/MS spectra from large-scale proteomics experiments using spectrum libraries. Anal Chem. 2006;78:5678–5684. [PubMed]
18. Lam H, Deutsch EW, Eddes JS, Eng JK, et al. Development and validation of a spectral library searching method for peptide identification from MS/MS. Proteomics. 2007;7:655–667. [PubMed]
19. Shilov IV, Seymour SL, Patel AA, Loboda A, et al. The Paragon Algorithm, a next generation search engine that uses sequence temperature values and feature probabilities to identify peptides from tandem mass spectra. Mol Cell Proteomics. 2007;6:1638–1655. [PubMed]
20. Eng JM, AL, Yates JR., 3rd An approach to correlate tandem mass spectral data of peptides with amino acid sequences in a protein database. J Am Soc Mass Spectrom. 1994;5:13. [PubMed]
21. Han DK, Eng J, Zhou H, Aebersold R. Quantitative profiling of differentiation-induced microsomal proteins using isotope-coded affinity tags and mass spectrometry. Nat Biotechnol. 2001;19:946–951. [PMC free article] [PubMed]
22. Keller A, Nesvizhskii AI, Kolker E, Aebersold R. Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal Chem. 2002;74:5383–5392. [PubMed]
23. Li XJ, Yi EC, Kemp CJ, Zhang H, Aebersold R. A software suite for the generation and comparison of peptide arrays from sets of data collected by liquid chromatography-mass spectrometry. Mol Cell Proteomics. 2005;4:1328–1340. [PubMed]
24. Ashburner M, Ball CA, Blake JA, Botstein D, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–29. [PMC free article] [PubMed]
25. Westfall PaYS. Resampling-Based Multiple Testing: Examples and Methods for p-Value Adjustment. 1993.
26. Nakayama K, Nakayama N, Wang TL, Shih Ie M. NAC-1 controls cell growth and survival by repressing transcription of Gadd45GIP1, a candidate tumor suppressor. Cancer Res. 2007;67:8058–8064. [PubMed]
27. Li XJ, Zhang H, Ranish JA, Aebersold R. Automated statistical analysis of protein abundance ratios from data generated by stable-isotope dilution and tandem mass spectrometry. Anal Chem. 2003;75:6648–6657. [PubMed]
28. Lau KW, Jones AR, Swainston N, Siepen JA, Hubbard SJ. Capture and analysis of quantitative proteomic data. Proteomics. 2007;7:2787–2799. [PMC free article] [PubMed]
29. Li Y, Sokoll LJ, Rush J, Zou N, Chan DW. Targeted detection of prostate cancer proteins in serum using heavy peptide standards and MALDI-TOF/TOF. Proteomics-Clinical Applications. 2009 In press.
30. Stahl-Zeng J, Lange V, Ossola R, Eckhardt K, et al. High sensitivity detection of plasma proteins by multiple reaction monitoring of N-glycosites. Mol Cell Proteomics. 2007;6:1809–1817. [PubMed]
31. Keshishian H, Addona T, Burgess M, Kuhn E, Carr SA. Quantitative, multiplexed assays for low abundance proteins in plasma by targeted mass spectrometry and stable isotope dilution. Mol Cell Proteomics. 2007;6:2212–2229. [PMC free article] [PubMed]
32. Anderson L, Hunter CL. Quantitative mass spectrometric multiple reaction monitoring assays for major plasma proteins. Mol Cell Proteomics. 2006;5:573–588. [PubMed]
33. Umezu T, Shibata K, Kajiyama H, Terauchi M, et al. Taxol resistance among the different histological subtypes of ovarian cancer may be associated with the expression of class III beta-tubulin. Int J Gynecol Pathol. 2008;27:207–212. [PubMed]
34. Orr GA, Verdier-Pinard P, McDaid H, Horwitz SB. Mechanisms of Taxol resistance related to microtubules. Oncogene. 2003;22:7280–7295. [PMC free article] [PubMed]
35. Kuo HC, Lee HJ, Hu CC, Shun HI, Tseng TH. Enhancement of esculetin on Taxol-induced apoptosis in human hepatoma HepG2 cells. Toxicol Appl Pharmacol. 2006;210:55–62. [PubMed]
36. Rostovtseva TK, Sheldon KL, Hassanzadeh E, Monge C, et al. Tubulin binding blocks mitochondrial voltage-dependent anion channel and regulates respiration. Proc Natl Acad Sci U S A. 2008;105:18746–18751. [PubMed]
37. Carre M, Andre N, Carles G, Borghi H, et al. Tubulin is an inherent component of mitochondrial membranes that interacts with the voltage-dependent anion channel. J Biol Chem. 2002;277:33664–33669. [PubMed]
38. Rostovtseva TK, Bezrukov SM. VDAC regulation: role of cytosolic proteins and mitochondrial lipids. J Bioenerg Biomembr. 2008;40:163–170. [PMC free article] [PubMed]