- Research article
- Open Access
Dynamic proteomics reveals bimodal protein dynamics of cancer cells in response to HSP90 inhibitor
BMC Systems Biologyvolume 11, Article number: 33 (2017)
Drugs often kill some cancer cells while others survive. This stochastic outcome is seen even in clonal cells grown under the same conditions. Understanding the molecular reasons for this stochastic outcome is a current challenge, which requires studying the proteome at the single cell level over time. In a previous study we used dynamic proteomics to study the response of cancer cells to a DNA damaging drug, camptothecin. Several proteins showed bimodal dynamics: they rose in some cells and decreased in others, in a way that correlated with eventual cell fate: death or survival. Here we ask whether bimodality is a special case for camptothecin, or whether it occurs for other drugs as well. To address this, we tested a second drug with a different mechanism of action, an HSP90 inhibitor. We used dynamic proteomics to follow 100 proteins in space and time, endogenously tagged in their native chromosomal location in individual living human lung-cancer cells, following drug administration.
We find bimodal dynamics for a quarter of the proteins. In some cells these proteins strongly rise in level about 12 h after treatment, but in other cells their level drops or remains constant. The proteins which rise in surviving cells included anti-apoptotic factors such as DDX5, and cell cycle regulators such as RFC1. The proteins that rise in cells that eventually die include pro-apoptotic factors such as APAF1. The two drugs shared some aspects in their single-cell response, including 7 of the bimodal proteins and translocation of oxidative response proteins to the nucleus, but differed in other aspects, with HSP90i showing more bimodal proteins. Moreover, the cell cycle phase at drug administration impacted the probability to die from HSP90i but not camptothecin.
Single-cell dynamic proteomics reveals sub-populations of cells within a clonal cell line with different protein dynamics in response to a drug. These different dynamics correlate with cell survival or death. Bimodal proteins which correlate with cell fate may be potential drug targets to enhance the effects of therapy.
Cancer drugs often kill some cells while other cells survive [1–5]. This stochastic outcome occurs even in clonal cells that are under identical conditions such as sister cells on the same plate. This stochastic resistance is non-genetic: The surviving cells, when re-plated, often give rise to populations that again show the same fraction of death versus survival in response to the drug [4, 6–8]. Inherited resistance evolves much slower, and usually occurs only after many such passages [3, 6, 9, 10]. The stochastic survival of cells may be one reason that cancer drugs do not always succeed in eliminating tumors, and understanding how some cells survive is therefore a pressing need.
In order to understand the molecular basis for the stochastic outcome of a drug, one needs to view the proteome in individual cells over time. Most existing proteomic methods average over millions of cells and therefore mask single-cell effects [1, 11]. Techniques for single-cell analysis based on immunostaining [12, 13] or transcriptomics  require fixing the cells and thus preclude studying the dynamics and eventual fate of each cell.
We have previously established a dynamic proteomics approach that addresses these issues and is able to follow proteins in single living human cancer cells over time. Dynamic proteomics is based on a library of cancer cell clones. Each clone expresses a full length tagged protein from its endogenous chromosomal locus [14–16]. We used this method to study the response of cells to the chemotherapy drug camptothecin (CPT) . CPT is a topoisomerase poison which causes DNA damage  in dividing cells. Survival and death of different cells was found not to be due to cell-cycle differences. Instead, several proteins were found with different dynamics in individual cells, which correlated with cell fate. These proteins were called bimodal proteins: their level rose 20 h after CPT treatment in some cells, but decreased in other cells. Two proteins rose primarily in cells that survived, DDX5 and RFC1. Knocking down these proteins enhanced killing by CPT, suggesting a causal effect .
Here we ask whether bimodality of protein dynamics is specific to CPT, or whether it occurs also for another drug. For this purpose we used dynamic proteomics to analyze the response to a drug with a different mechanism of action, an HSP90 inhibitor (HSP90i). The HSP90i class of drugs inhibits the chaperone HSP90 and therefore interferes with protein degradation, which is thought to affect cancer cells more strongly than other cells [18–22].
We used dynamic proteomics to study how cells respond to an HSP90 inhibitor. We followed 100 proteins in space and time in living cancer cells following HSP90i treatment. We find 16 bimodal proteins whose protein dynamics are correlated with cell fate (survival/death). Seven of these were bimodal also in the previous CPT study. This indicates the formation of subpopulations of cells 12 h post treatments within a clonal cell line, long before death begins to at about 20 h. Unlike CPT, we find that HSP90i efficacy is correlated with the cell cycle. We further find localization changes of oxidative-stress response proteins which indicate that oxidative stress response is activated about 20 h after treatment. The dynamic proteomic approach thus suggests that bimodality may be a widespread response of cells to anti-cancer drugs, indicating that subpopulations with differential survival are induced by the drug many hours before death occurs.
Fluorescently-tagged protein clone library
The clones in this study were taken from the LARC library, in which proteins were fluorescently-tagged in their native locus and under their native regulation by CD-tagging as previously described [2, 23]. Briefly, H1299 (non-small cell lung carcinoma commercial cell line) cells were double tagged using the CD-tagging method, such that all clones contains a red fluorescent mCherry protein fused to a protein that is brighter in the nucleus and dimmer in the cytoplasm and is used for automated image analysis, and each clone expresses a yellow tag (eYFP or Venus) fused to the protein of interest. Each clone in the library originated from a single cell obtained by sorting for YFP fluorescence.
The H1299 cells that used to establish the LARC library originated from a commercial cell line, as described in Cohen et al. .
For this study we chose a subset of clones from the LARC library. Clones were selected to provide coverage of different cellular pathways (such as: cell division, apoptosis, cell signaling), protein functions (such as: kinase, phosphatase, protease) and localization (such as: nuclei, nucleoli, cytoplasm, membrane). We chose 14 proteins which were found to be bimodal in the previous study of CPT (Table 1).
Tissue Culture Media
Cells were grown in RPMI 1640 supplied with (+) L-Glutamine (GIBCO, catalog number 21875) medium supplemented with 10% Fetal Calf Serum (certified fetal bovine serum, membrane filtered, Biological Industries, catalog number 04-001-1A) and 0.05% Penicillin-Streptomycin antibiotics (Biological Industries, catalog number 03-031-1B), in incubators at 37 °C and 8% CO2.
HSP90 inhibitor (MSC1972516 (EMD-614684)) , was dissolved in DMSO (hybri-max, D2650 Sigma) giving a stock solution of 1 mM. The drug was diluted to a final concentration of 1 μM by 1:1000 dilution in transparent growth medium (RPMI 1640, 0.05% Penicillin-Streptomycin antibiotics, 10% FCS, with L-Glutamine, lacking riboflavin and phenol red, Bet Haemek, Biological Industries Cat. No. 06-1100-26-1A). Cells were grown in 96-well plates (7 × 103 cells per well) with normal (red) RPMI medium. Before each experiment, normal (red) RPMI medium was replaced by transparent RPMI 1640 (lacking riboflavin and phenol red to decrease medium fluorescence) with the desired drug (volume of 200ul per well). HGF (Hepatocyte Growth factor, Recombinant Human HGF, Catalog Number: 294-HGN, reconstituted at 50 μg/mL in sterile PBS containing 0.1% human or bovine serum albumin) was added at a final concentration of 100 ng/ml.
Time-lapse movies were obtained at 20× magnification, in IDEA Bio-Medical Ltd. WiScan system which includes Olympus inverted fluorescence microscope (modified IX71), live cell environmental incubators maintaining 37 °C (37-2 digital and Heating unit, PeCon, Germany, Leica #15531719), humidity and 8% CO2 (PeCon, GmbH, Germany #0506.000-230, Leica #11521733) and automated stage movement control (Corvus, ITK, GmbH, Germany); the stage was surrounded by a custom enclosure to maintain constant temperature, CO2 concentration, and humidity. Transmitted and fluorescent light paths were controlled by electronic shutters (Uniblitz, model VMM-D1, Rochester, NY); Fluorescent light sources were Short ARC Lamp HXP R 120w/45C VIS (OSRAM, Germany). iXon (Andor technology) cooled 14 bit CCD camera was used. The filters used were from Chroma Technology: Olympus Dual channel dichroic with excitation and emission filters mounted on fast filter wheels (Chroma set 490/20x, 577/25x, 535/30 m, 632/60 m).
Cells were visualized in 96-well optical glass-bottom plates (96 well optical CVG sterile w/lid Black, catalog number 164588 Thermo scientific nunc) coated with 10 μM fibronectin 0.1% (solution from bovine plasma, Sigma, Cat. No. F1141) diluted 1:100 in Dulbecco’s Phosphate Buffered Saline, PBS (Sigma, Cat. No. D8537). For each well, time-lapse movies were obtained at four fields of view. Each movie was taken for 2 days with 30 min resolution (100 time points). Each time point included three frames: a transmitted light image, and two fluorescent channels (red and yellow). In the present study we use the first 24 h of each movie.
Image Analysis of Time-Lapse Movies
The image analysis software described in Cohen et al.  was used to analyze the time lapse movies in this study. The main steps in this software include background correction (flat field and background subtraction), segmentation, cell tracking, and automated identification of cell phenotypes (mitosis and cell death). Global image threshold (Otsu, 1979) followed by watershed segmentation was used to segment the nuclei and cytoplasm based on the red fluorescent images of the red tagged protein found in all clones. Tracking was performed by analysis of the movie from end to start and linking of each segmented cell to the cell in the previous image with the closest centroid, and other attributes like protein levels. Automated cell death identification algorithm utilized the morphological changes such as envelope breakdown, and cell roundness correlated with dying cells, and cells that undergo mitosis.
Definition of clones with bimodal dynamics
For every clone, we calculated the slopes between every two time-points starting from t0 = 12 h after drug administration until t1 = 20 h after drug administration (total of 17 time points). We calculated for every cell the mean and median of the slopes between t0 and t1. Mean and median calculation showed similar results. Next we divided the cells in the clone into two groups: cells that had a positive slope on average–meaning that the protein level increased, and cells that had negative slope on average–meaning a decrease in protein level. We removed clones in which there were fewer than 10 cells in at least one of the two groups, or in which the movie quality was judged to be insufficient for analysis, remaining with 75 clones. We calculated the distribution of mean slopes in each group (positive and negative mean slope), and applied three tests for the null hypothesis H0 that these two groups are drawn from the same distribution: t-test, KS test and the nonparametric Mann Whitney test. To control for multiple hypothesis testing we applied the Benjaminy-Hochberg procedure with FDR = 0.1 to the results from all three tests. This result in a list of 25 proteins in which H0 can be rejected (p < 0.05 in all three tests).
Bimodal dynamics and cell fate
We analyzed the 25 bimodal clones in terms of cell survival or death at 24 h. We divided the cells into two groups according to the fate of the cells, compared the slope distribution in these two groups, as described in the section above. We find that for 16 clones, cell fate significantly differed in the two groups (t-test, KS and MW tests, p < 0.05, Benjamini Hochberg corrected with FDR = 0.1).
To study how cells respond to an HSP90 inhibitor we used dynamic proteomics. We measured the dynamics of 100 proteins involved in a wide range of cellular functions including signaling, growth and death, in response to 1 uM of an Hsp90 inhibitor (Merc Serrono MSC1972516) (proteins listed in Table 1). Cell death occurred 15 h+/−3 after drug addition. Protein dynamics were assayed using clones from the LARC library, described in . In these clones of the H1299 lung-cancer cell line, proteins were tagged with YFP at their endogenous chromosomal location (Fig. 1) using a CD-tagging approach. Previous work showed that the clones express full length tagged proteins from their native locus with intact regulatory sequences. All clones also express red fluorescent proteins that aid fully-automated image analysis. Dynamics of protein level at a resolution of 30 min were examined for 24 h after drug addition in 96-well plates in a dedicated fluorescence microscopy system with controlled CO2 and temperature. Protein YFP level and localization, as well as cell mitosis and death events were automatically tracked and quantitated for hundreds of individual cells per clone as described  (Fig. 1). Protein YFP levels averaged over all cells showed diverse dynamics with most proteins dropping in levels on average and several other rising at late times (Fig. 2). Such dynamics that are averaged over all cells mask single cell-effects, to which we now turn.
A quarter of the proteins show bimodal behavior, decaying in some cells and increasing in others, in a way that correlates with cell fate
We studied the dynamics of the proteins at the individual cell level. We used the criteria of  to define bimodal dynamics: we computed the rate of accumulation of YFP protein level for each cell defined by the slope of the YFP curve over time. We separated the cells into two groups, with positive and negative slopes. Bimodal proteins were proteins from which the slope distribution of the two groups was significantly different as judged by using statistical tests (see Methods, t-test results shown in Table 2). We find that 25 out of 100 proteins (listed in Table 2) displayed a bi-modal response to HSP90i with expression increasing in some cells but decreasing in other cells 12 h after drug addition (Fig. 3).
For 16 out of the 25 bimodal proteins the protein dynamics correlated with cell fate (death or survival at the end of the experiment). These proteins and their function are listed in Table 2. For 12 proteins, such as STMN1, a protein important for mitosis spindle organization, cells that increase the protein level are preferentially killed (p = 7e−15). Many of these proteins are not naturally associated with cell death, including proteins involved in mitosis and in anti-apoptotic functions. The increase of 4 proteins correlated with cell survival, including DDX5 and RCF1. For 9 of the bimodal proteins we found no significant correlation with cell death or cell survival.
Six of the bimodal proteins, DDX5, RFC1, BAG1, BAG2, SPCS1, CALM2 and PCMTD2 were previously identified as bimodal also for CPT. The other 18 bimodal proteins found here were not bimodal in CPT, suggesting that the stochastic mechanisms in the response to the two drugs have only partial overlap.
Cells that attempt mitosis 12 h or more after drug addition are preferentially killed
We next asked whether the cell-cycle phase at the moment of drug administration affects cell survival. We automatically detected mitosis events as described , using morphological changes of cells, primarily cell rounding (see Methods) [2, 14, 16, 24]. We find that cells that attempt mitosis 12 h or more after drug administration were not able to complete cell division and were preferentially killed (82%), whereas cells that show the morphological correlate of mitosis in the first 12 h were able to complete mitosis and were less frequently killed (35%), Fig. 4. This finding suggests that HSP90i effect depends on the cell-cycle stage of the cells. In contrast CPT showed no dependence on cell cycle stage, highlighting another difference between the drugs at the level of individual cell responses.
Oxidative stress proteins increase in levels and enter the nucleus after drug addition
The dynamic proteomic approach also allows detection of localization changes in response to drugs. We observe translocations of thioredoxin (TXN) and thioredoxin reductase (TrxR), which show nuclear entry after 20 h (Fig. 5). This translocation occurs upon oxidative stress response [25–27]. A similar translocation of these two proteins occurred also for CPT . Oxidative stress response is an off-target effect for these two drugs.
We asked whether bimodal protein dynamics, which were found for CPT, also occur for another drug with a different mechanism of action. We therefore used the dynamic proteomics approach to study how individual cancer cells respond to an HSP90 inhibitor. We find that a quarter of the proteins show bimodal dynamics: they rise in some cells and drop in others. This bimodality begins about 12 h post-drug, and indicates the existence of cell subpopulations. Most of the bimodal proteins have dynamics that correlate with eventual cell fate: death or survival.
The existence of cell subpopulations with different dynamics may explain why some cells survive and others die in response to the drug. The proteins that rise in correlation with cell survival may be potential drug targets to enhance HSP90i effects. In previous work two of the bimodal protein that rose preferentially in surviving cells, DDX5 and RFC1, were shown to increase effectiveness of CPT upon knockdown . Two additional proteins were found in the present study that rose in correlation with survival: a MAPK that lies upstream of ERK1,2 which signals for proliferation, and a tubulin stabilizing protein TBCA. Future work can test the causal link of these proteins to survival under HSP90i.
Proteins that rose in cells that eventually die include mitosis and anti-apoptotic proteins. This result is unexpected because one might expect death-related proteins to rise in cells that are eventually killed rather than proteins associated with cell division and survival. One way to understand this finding is our observation that the cells that attempt mitosis more than 12 h after drug addition are preferentially killed. These cells may attempt to progress through the cell cycle and to complete mitosis in the face of the drug but fail to do so.
The present single cell approach also allowed us to identify localization changes indicative of oxidative stress. Thioredoxin and thioredoxin reductase both move to the nucleus 20 h post drug administration . Such off-target effects may be relevant in understanding and addressing the toxicity of HSP90i drug candidates. Similar localization changes were observed also after camptothecin treatment .
Differences in cell cycle phase play a major role in cell-cell variability. We find that cells that attempted mitosis in the first 12 h of drug treatment were preferentially killed. This is in contrast with camptothecin that showed no measurable cell-cycle dependence. An advantage of the present approach is that it can directly view mitosis events. This non-invasive measurement avoids some of the concerns associated with chemical synchronization of cells to study cell-cycle effects. This advantage allowed dynamic proteomics to accurately detect cell-cycle dependent proteins in previous studies [28–30].
Future work can address the origin of the subpopulations found here. One possibility is that the subpopulations are formed dynamically after the drug is added. Another possibility is that subpopulations exist before the drug treatment due to genetic mutation or a slowly reversible epigenetic change . Previous work on the H1299 cell line, the parental cell line for the present clones, suggests that it may harbor preexisting subpopulations. For H1299, cisplatin treatment of unsorted or stem cell marker CD44+ cells resulted in no significant increase of apoptotic cells compared to the untreated control, while CD44- cells showed significant increase in apoptosis after cisplatin treatment . H1299 telomere-lengthening correlates with telomerase activity, but exhibits clonal variability in telomere-lengthening and telomerase activity . Suboptimal doses of cisplatin on H1299 cells and re-plating select for an enhanced malignant phenotype . These findings are in line with the idea that H1299 includes at least two different subclones with different drug-sensitivities.
It would be intriguing to further study the nature of the cell subpopulations in more detail. One can ask whether bimodal proteins go up and down in a coordinated way in each cell, or whether they are uncoordinated. This may require multi-color tagging of these proteins . It would be important to test other drugs and cell lines in order to explore the generality of bimodality that correlates with cell fate.
This study used single-cell dynamic proteomics to reveal sub-populations of cells that have different protein dynamics in response to a drug within a clonal cell line, in a way that correlates with cell survival or death. Bimodal proteins which correlate with cell fate may be potential drug targets to enhance the effects of therapy.
Apoptotic Peptidase Activating Factor 1 (gene)
BCL2 Associated Athanogene 1 (gene)
BCL2 Associated Athanogene 2 (gene)
Calmodulin 2 (gene)
Central Dogma tagging (a method to label genes by integrating a DNA sequence coding for a fluorescent tag into the genome)
Camptothecin (an anti-cancer drug)
DEAD (Asp-Glu-Ala-Asp)-Box Helicase 5 (gene)
Extracellular signal-regulated kinases 1 or 2 (genes)
Human non-small cell lung carcinoma cell line
Heat shock protein 90
Heat shock protein 90 inhibitor
- LARC library:
Library of annotated reporter cell clones, a collection of proteins which were fluorescently tagged in their endogenous chromosomal locus
Mitogen-activated protein kinas (a type of protein kinase)
Protein-L-Isoaspartate (D-Aspartate) O-Methyltransferase Domain Containing 2 (gene)
Replication Factor C Subunit 1 (gene)
Signal Peptidase Complex Subunit 1 (gene)
Stathmin 1 (gene)
Tubulin Folding Cofactor A (gene)
Thioredoxin reductase (gene)
Altschuler SJ, Wu LF. Cellular Heterogeneity: Do Differences Make a Difference? Cell. 2010;141:559–63.
Cohen AA, Geva-Zatorsky N, Eden E, Frenkel-Morgenstern M, Issaeva I, Sigal A, et al. Dynamic proteomics of individual cancer cells in response to a drug. Science. 2008;322:1511–6.
Marusyk A, Almendro V, Polyak K. Intra-tumour heterogeneity: a looking glass for cancer? Nat Rev Cancer. 2012;12:323–34.
Sharma SV, Lee DY, Li B, Quinlan MP, Takahashi F, Maheswaran S, et al. A chromatin-mediated reversible drug-tolerant state in cancer cell subpopulations. Cell. 2010;141:69–80.
Tirosh I, Izar B, Prakadan SM, Wadsworth MH, Treacy D, Trombetta JJ, et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science. 2016;352:189–96.
Flusberg DA, Sorger PK. Surviving apoptosis: life–death signaling in single cells. Trends Cell Biol. 2015;25:446–58.
Snijder B, Pelkmans L. Origins of regulated cell-to-cell variability. Nat Rev Mol Cell Biol. 2011;12:119–25.
Spencer SL, Gaudet S, Albeck JG, Burke JM, Sorger PK. Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis. Nature. 2009;459:428–32.
Brock A, Chang H, Huang S. Non-genetic heterogeneity—a mutation-independent driving force for the somatic evolution of tumours. Nat Rev Genet. 2009;10:336–42.
Meacham CE, Morrison SJ. Tumour heterogeneity and cancer cell plasticity. Nature. 2013;501:328–37.
Wang D, Bodovitz S. Single cell analysis: the new frontier in “omics.”. Trends Biotechnol. 2010;28:281–90.
Yuan X, Curtin J, Xiong Y, Liu G, Waschsmann-Hogiu S, Farkas DL, et al. Isolation of cancer stem cells from adult glioblastoma multiforme. Oncogene. 2004;23:9392–400.
Zhang S, Balch C, Chan MW, Lai H-C, Matei D, Schilder JM, et al. Identification and Characterization of Ovarian Cancer-Initiating Cells from Primary Human Tumors. Cancer Res. 2008;68:4311–20.
Cohen AA, Kalisky T, Mayo A, Geva-Zatorsky N, Danon T, Issaeva I, et al. Protein Dynamics in Individual Human Cells: Experiment and Theory. PLoS ONE. 2009;4:e4901.
Frenkel-Morgenstern M, Cohen AA, Geva-Zatorsky N, Eden E, Prilusky J, Issaeva I, et al. Dynamic Proteomics: a database for dynamics and localizations of endogenous fluorescently-tagged proteins in living human cells. Nucleic Acids Res. 2010;38:D508–12.
Issaeva I, Cohen AA, Eden E, Cohen-Saidon C, Danon T, Cohen L, et al. Generation of Double-Labeled Reporter Cell Lines for Studying Co-Dynamics of Endogenous Proteins in Individual Human Cells. PLoS ONE [Internet]. 2010;5. Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2958823/. cited 25 Feb 2016.
Hsiang YH, Hertzberg R, Hecht S, Liu LF. Camptothecin induces protein-linked DNA breaks via mammalian DNA topoisomerase I. J Biol Chem. 1985;260:14873–8.
Sidera K, Patsavoudi E. HSP90 Inhibitors: Current Development and Potential in Cancer Therapy. Recent Pat Anticancer Drug Discov. 2014;9:1–20.
Proia DA, Bates RC. Ganetespib and HSP90: translating preclinical hypotheses into clinical promise. Cancer Res. 2014;74:1294–300.
Porter JR, Fritz CC, Depew KM. Discovery and development of Hsp90 inhibitors: a promising pathway for cancer therapy. Curr Opin Chem Biol. 2010;14:412–20.
Jhaveri K, Modi S. Ganetespib: research and clinical development. OncoTargets Ther. 2015;8:1849–58.
Whitesell L, Lindquist SL. HSP90 and the chaperoning of cancer. Nat Rev Cancer. 2005;5:761–72.
Sigal A, Danon T, Cohen A, Milo R, Geva-Zatorsky N, Lustig G, et al. Generation of a fluorescently labeled endogenous protein library in living human cells. Nat Protoc. 2007;2:1515–27.
Issaeva I, Cohen AA, Eden E, Cohen-Saidon C, Danon T, Cohen L, et al. Generation of Double-Labeled Reporter Cell Lines for Studying Co-Dynamics of Endogenous Proteins in Individual Human Cells. PLoS ONE. 2010;5:e13524.
Lee K-H, Jang A-H, Yoo C-G. 17-Allylamino-17-Demethoxygeldanamycin and the Enhancement of PS-341–Induced Lung Cancer Cell Death by Blocking the NF-κB and PI3K/Akt Pathways. Am J Respir Cell Mol Biol. 2015;53:412–21.
Byrne BM, Welsh J. Altered thioredoxin subcellular localization and redox status in MCF-7 cells following 1,25-dihydroxyvitamin D3 treatment. J Steroid Biochem Mol Biol. 2005;97:57–64.
Karimpour S, Lou J, Lin LL, Rene LM, Lagunas L, Ma X, et al. Thioredoxin reductase regulates AP-1 activity as well as thioredoxin nuclear localization via active cysteines in response to ionizing radiation. Oncogene. 2002;21:6317–27.
Geva-Zatorsky N, Dekel E, Cohen AA, Danon T, Cohen L, Alon U. Protein Dynamics in Drug Combinations: a Linear Superposition of Individual-Drug Responses. Cell. 2010;140:643–51.
Farkash-Amar S, Eden E, Cohen A, Geva-Zatorsky N, Cohen L, Milo R, et al. Dynamic proteomics of human protein level and localization across the cell cycle. PLoS One. 2012;7:e48722.
Sigal A, Milo R, Cohen A, Geva-Zatorsky N, Klein Y, Alaluf I, et al. Dynamic proteomics in individual human cells uncovers widespread cell-cycle dependence of nuclear proteins. Nat Methods. 2006;3:525–31.
Leung EL-H, Fiscus RR, Tung JW, Tin VP-C, Cheng LC, Sihoe AD-L, et al. Non-Small Cell Lung Cancer Cells Expressing CD44 Are Enriched for Stem Cell-Like Properties. PLoS One. 2010;5:e14062.
Savre-Train I, Gollahon LS, Holt SE. Clonal heterogeneity in telomerase activity and telomere length in tumor-derived cell lines. Proc Soc Exp Biol Med Soc Exp Biol Med N Y N. 2000;223:379–88.
Hsieh J-L, Lu C-S, Huang C-L, Shieh G-S, Su B-H, Su Y-C, et al. Acquisition of an enhanced aggressive phenotype in human lung cancer cells selected by suboptimal doses of cisplatin following cell deattachment and reattachment. Cancer Lett. 2012;321:36–44.
This work was supported by the Israel Science Foundation and the Clore-Katz Weizmann internal grant.
Availability of data and materials
The datasets generated during and/or analyzed during the current study available from the corresponding author on reasonable request.
AZ designed and preformed the experiments and wrote the paper. SA participated in the design of the study and the image analysis. TD participated and assisted in the experimental procedures. UA designed the study and wrote the paper. All authors read and approved the final manuscript.
UA is the incumbent of the Abisch-Frenkel Professorial Chair.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate