| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |




From the Centre dOncologia Molecular,* Institut de Recerca Oncològica-Institut de Investigació Biomèdica de Bellvitge, Hospital Duran i Reynals, Barcelona; and the Grup de Bioinformàtica Estructural (Unitat de Recerca dInformàtica Biomèdica-Institut Municipal dInvestigacions Mediques)
and the Departament de Ciències Experimentals i de la Salut,
Unitat de Proteòmica, Universitat Pompeu Fabra, Barcelona, Spain
| Abstract |
|---|
|
|
|---|
Because the disseminating tumor cells are freely and ubiquitously distributed by the hematogeneous or lymphatic system, the arrival to many organs is not sufficient for the development of secondary tumors. Unless a growing tumor colony is established at a new site, the metastatic process is not fulfilled because the local organ microenvironment at the site where the tumor cells are lodged determines whether metastases emerge or not.5
Tumor dormancy describes a prolonged quiescent state during which metastasis progression is not clinically detected.6 It has been suggested that dormancy might be related to anti-angiogenic factors that indirectly promote apoptosis and thus oppose proliferative tendencies.7 Indeed, metastatic cells may be present after surgery but remain dormant due to an inability to induce angiogenesis, or to change the balance between other growth-inducing/inhibiting factors in the tumor microenvironment, a circumstance that may also determine the length of the period between dissemination and the appearance of clinical metastases.8,9 How cells that have been selected at the primary site for acquisition of self-sufficiency can suppress these activities and remain dormant for years is an enigma. It has been suggested that tumor cells may disseminate in a far less progressed genomic state than previously thought, acquiring genomic aberrations typical of metastatic cells thereafter.10 Other groups have found that solid tumors carrying a gene-expression signature were mostly associated with metastasis and poor clinical outcome, suggesting that the metastatic potential of human tumors is encoded in the bulk of a primary tumor and that tumors likely to metastasize are fundamentally different.11-13
The molecular and cellular mechanisms responsible for the metastatic phenotype in breast cancer involve, among other factors, a number of gene products that participate in apoptosis.14-17 Anti-apoptotic genes have a role in displacing the balance between death and proliferation factors toward growth, possibly shortening the period between dissemination and the appearance of clinical metastasis.8
Bcl-xL expression in breast cancer cells increases metastatic activity. This might result from resistance to apoptosis against cytokines, increasing cell survival in circulation, and enhancing anchorage-independent growth.18 We have described that overexpression of anti-apoptotic Bcl-xL plays a role in breast cancer dormancy, promoting survival of cells in metastatic foci by counteracting the proapoptotic signals in the microenvironment and favoring the successful development of metastases in specific organs,19 preferentially lodging in peripheral lymph nodes.20 Thus, by selecting the organ-specific most adaptive phenotype, anti-apoptotic proteins might be a hallmark of metastasis and resistance to therapies.21 The aim of this work is to provide insight into the metastasis phenotype of breast cancer cells overexpressing Bcl-xL, which results in a short dormancy period in several organs.
The systematic large scale profiling of protein expression by proteomic methods allows the molecular characterization of cellular events associated with cancer progression.22 We have used two-dimensional gel electrophoresis (2DE), with subsequent identification of specific proteins of interest by mass spectrometry, to define phenotypes that in vivo enhance the ubiquity of clinical metastasis of Bcl-xL-overexpressing cells. On the assumption that the progenitor of the later-arising metastasis must be present among those cells that have disseminated to distant sites, our analysis included the orthotopic tumor, lymph node, lung and bone metastasis, and the previously originated blood stream survival cells (BSCs).18
We have identified 17 proteins that are either overexpressed (more than twofold) or underexpressed (less than twofold) in metastases with regard to the tumor. Furthermore, using recently developed protein-protein interaction prediction algorithms we have proposed interactors for proteins for which no reliable interactions could be inferred from experimental data. We have used protein-protein interaction databases and in silico prediction approaches to find cell functions relevant to the 17 identified proteins and which could guide additional research into expression patterns and metabolic and signal pathways.
| Materials and Methods |
|---|
|
|
|---|
MDA-MB 435 cell cultures (435) were maintained in 1:1 (v/v) mixture of Dulbeccos modified Eagles medium and Hams F12 medium supplemented with 10% fetal bovine serum, 1 mmol/L piruvate, and 2 mmol/L L-glutamine in 5% CO2-95% air at 37°C in a humidified incubator. Transfections of pSFFV-Neo Bcl-xL or pSFFV-Neo were performed with Lipofectin (Life Technologies, Inc., Gaithersburg, MD). Selection of 435/Bcl-xL and 435/Neo cells started 48 hours after transfection using 500 µg/ml of neomycin G418 (Life Technologies). The pool of transfectants was used to avoid the possibility of clonal variations.
We generated orthotopic primary tumors in 7-week-old athymic nude BALB/c female mice by inoculation of 1 x 106 cells in 0.05 ml of medium without serum in the right inguinal mammary gland (i.m.f.p.), following previous protocols.18
Tumors grew at day 45 from cell injection, and animals were put down and organs removed, weighed, and examined for metastasis at day 110. Metastatic involvement was explored in paraffin sections by hematoxylin and eosin (H&E) stain. Lung metastasis from tumors overexpressing Bcl-xL reached more than 3 mm in diameter, whereas only small metastases of
3 mm were found in the lungs of 435/Neo-injected mice. Monolayers of metastatic cells were obtained from trypsin-treated histocultures (metastasis fragments of
1 mm3) maintained until growth in medium supplemented with 20% fetal bovine serum, in the presence of 500 µg/ml of neomycin G418. Because lymph node and bone metastasis appeared in mice from Bcl-xL clones but not in mice from 435/Neo control cells, metastatic variants were established from bone, lung, and lymph node metastasis from mice carrying 435/Bcl-xL tumors, and lung metastasis from mice carrying 435/Neo tumors.
Previously, we had obtained and characterized survival cells in the blood stream (BSCs).18 Briefly, a mixture of cells 3:1, 435/Neo and 435/Bcl-xL, respectively, was injected into the tail vein of nude mice. We assessed differences in survival rates between the two clones by comparing the original mixture and cell populations recovered from the lungs. Mice were killed at 15 minutes, 2, 4, 8, and 16 hours, and recovered cells were expanded in culture and the DNA was hybridized with a radiolabeled 0.8-kb fragment of Bcl-xL to analyze the Bcl-xL cDNA in survival clones. We used for experiments as a BSC the cell variant recovered from lungs 4 hours after injection into the tail vein of nude mice.
Sample Preparation
Cultured cells derived from metastasis and primary tumors were used for proteomic analyses. Cells (1 x 106) were plated onto a 75-mm flask in complete medium for 48 hours, starved for 24 hours in serum-free medium, and replaced by complete medium for 24 hours more. Cells were rinsed in ice-cold phosphate-buffered saline (PBS) before trypsinization with 0.05% trypsin and 0.5 mmol/L ethylenediamine tetraacetic acid. Variants in which a floating subpopulation was growing were preserved at each pass, by recovering the PBS wash before trypsin treatment. The cell pellets were frozen in N2 and stored at 80°C until electrophoresis. For this purpose, 1 x 106 cells were solubilized in 200 µl of lysis buffer containing 8 mol/L urea, 4% CHAPS (w/v), and 40 mmol/L Tris, and then homogenized by passing through a 25-gauge needle. Insoluble material was removed by centrifugation at 13,000 rpm for 10 minutes at 4°C. Protein concentration was determined by BCA protein assay reagent (Pierce, Rockford, IL).
Two-Dimensional Gel Electrophoresis (2DE)
Immobilized pH 4 to 7 gradient (IPG) 18-cm strips, (Immobiline DryStrip 4-7 NL; Amersham Pharmacia Biotech AB, Uppsala, Sweden) were rehydrated with 350 ml of rehydration solution (8 mol/L urea, 4% CHAPS, 1% IPG buffer, 13 mmol/L dithiothreitol), in which 200 µg of protein were applied. Isoelectric focusing was performed using a IPGphor apparatus (Amersham Pharmacia Biotech) for a total of 80 kV/hour, 50 mA/strip at 20°C, according to the manufacturers instructions, in a first-dimension isoelectric focusing system (Amersham Pharmacia Biotech). Immediately after being focused, IPG strips were equilibrated for 15 minutes in 6.5 mmol/L dithiothreitol, 50 mmol/L Tris, 6 mol/L urea, 30% glycerol, and 2% sodium dodecyl sulfate, followed by 13.5 mmol/L iodoacetamide, 50 mmol/L Tris, 6 mol/L urea, 30% glycerol, and 2% sodium dodecyl sulfate, 15 minutes before running in the second dimension, which was performed on 12.5% polyacrylamide gels of 18 x 20 cm in an electrophoresis tank Protean II (Bio-Rad, Hercules, CA). For protein identification (see below) preparative 2DE gels loaded with 300 µg of protein were run following the same protocol. After electrophoresis, gels were fixed and stained using a silver-staining low-fixation protocol compatible with mass spectrometry, dried between cellophane sheets, scanned in a densitometer, and processed with the Molecular Analyst Software (Bio-Rad, Richmond, CA) into Adobe Photoshop.
Analysis of Gel Images
The digitalized 2DE gel images were normalized and comparatively analyzed with the PDQuest program (version 7.0; Bio-Rad, Hercules, CA). Four gels from each specimen, in two independent experiments, were analyzed and the mean gel was used for comparisons. Normalization was assessed using the formula: normalized spot quantity = raw spot quantity x scaling factor (106 ppm) x normalization factor (total density in gel image). Significant changes were ratios <0.5, underexpressed, or >2, overexpressed, with regard to the referent sample. Molecular weight (MW) and isoelectric point (pI) were determined after comparison with reference gels in the SWISS-2DPAGE (Expasy) database.
Protein Fingerprinting by Matrix-Assisted Laser Desorption Ionization/Time of Flight (MALDI-TOF) Mass Spectrometry
Protein spots were selected according to the silver stain signal without overlap. The selected spots were excised from 2DE silver-stained gels and submitted to trypsin in-gel digestion as described before23
with minor modifications. After reduction and alkylation the gel pieces were swollen in 10 µl of a digestion buffer containing 50 mmol/L NH4HCO3 and 5 ng/µl of trypsin (modified porcine, sequencing grade; Promega, Madison, WI) in an ice bath. After 45 minutes 30 µl of 50 mmol/L NH4HCO3 were added and the digestion allowed to proceed overnight at 37°C. The samples were acidified using 5% formic acid and 10 to 20 µl of the supernatant containing tryptic peptides purified in a Poros R2 column24
(Applied Biosystems, Framingham, MA) packed in a constricted GelLoader Tip (Eppendorf, Hamburg, Germany) previously washed with 10 µl of acetonitrile and equilibrated with 10 µl of 5% formic acid. The sample was washed with 10 µl of 5% formic acid and eluted directly to the MALDI plate with 0.8 µl of matrix solution (20 mg/ml
-cyano-4-hydroxycinnamic acid in 70% in 0.1% formic acid). Peptide mass fingerprinting spectra were recorded in a Voyager STR MALDI-TOF (Applied Biosystems) in positive reflector mode with delayed extraction. The spectra were analyzed using the program m/z (Proteometrics, New York, NY). Protein identification was performed against a nonredundant database (NCBI) using the MASCOT program (http://www.matrixscience.com).
Experimental Validation of Specific Proteins Expressed in Metastasis
Western Blot Analysis
Cells from exponential cultures were lysed in 200 µl of RIPA buffer (50 mmol/L Tris, 150 mmol/L NaCl, 0.1% sodium dodecyl sulfate, 1% Nonidet P-40, 0.5% sodium deoxycholate). Sample volumes were adjusted to contain 50 µg of protein as determined by the BCA protein assay reagent (Pierce), and electrophoresed in a 12% polyacrylamide gel. The separated proteins were transferred to polyvinylidene difluoride membranes (Immobilon-p; Millipore Corp., Bedford, MA) and nonspecific protein-binding sites blocked by a 5% solution of nonfat dried milk with 0.01% Tween 20 in PBS. The following antibodies (Abs) were used: polyclonal rabbit Ab specific for human Bcl-x (S-18) at 1/1000 dilution (Santa Cruz Biotechnology, Santa Cruz, CA); polyclonal rabbit Ab anti-stathmin 1 (OP18) at 1/1500 dilution (Sigma, St. Louis, MO); monoclonal Ab anti-peroxiredoxin 3 (AOP-1) at 1/8000 dilution (Sigma); polyclonal rabbit Ab anti-peroxiredoxin 2 at 1/1000 dilution (LabFrontier, Seoul, Korea); polyclonal goat Ab anti-enolase 1 (eno
) at 1/500 dilution (Santa Cruz Biotechnology); polyclonal rabbit Ab anti-glutathione transferase (GST P1-1) at 1/1000 dilution (Calbiochem, CN Biosciences Ltd., Nottingham, UK); a monoclonal anti-tyrosine 3/tryptophan 5-monooxygenase, 14-3-3-
(clone 12) at 1/500 dilution (BD Biosciences, Pharmingen, Franklin Lakes, NJ). Peroxidase-conjugated goat anti-rabbit secondary antibody 1/2000 (Amersham), or anti-mouse secondary antibody 1/2000 (Pierce, Perbio Science Ltd., Cheshire, UK), or anti-goat secondary antibody 1/3000 (Santa Cruz Biotechnology) were used as corresponded in each case. Immunoreactive bands were visualized on CL-x posure Film (Pierce) using the Super Signal West-Pico (Pierce). MWs were established using Benchmark prestained MW markers (Invitrogen, San Diego, CA). An anti-human actin monoclonal antibody 1/2000 (Sigma) was also used as an internal standard for densitometric analysis of X-ray film. Densitometry analysis with the Quantity One program was performed using the quantity of a band that it is the sum of the intensities of all of the pixels within the band boundary multiplied by the area of each pixel.
Immunohistochemistry
Paraffin-embedded tissue was cut into 6-µm sections, deparaffinized, and rehydrated. On antigen retrieval with citrate buffer (pH, 6), we incubated with the following primary Abs: anti-stathmin 1 (OP18) at 1/500 dilution (Sigma); anti-peroxiredoxin 3 (AOP-1) at 1/500 dilution (Sigma), and anti-enolase 1 (eno
) at 1/100 dilution (Santa Cruz Biotechnology). Bound antibody was visualized with biotinylated anti-IgG (Vector Laboratories, Burlingame, CA) after diaminobenzidine staining and counterstaining with hematoxylin in an Olympus BX60 (Olympus Optical Co., Ltd., Tokyo, Japan) and digitalized in a digital camera Spot using the Spot 4.2 software (Diagnostic Instruments, Inc., Sterling Heights, MI).
Bioinformatics
We used the STRING (search tool for the retrieval of interacting genes/proteins25 ) to find interactors of the proteins identified by mass spectrometry. STRING contains protein-protein interactions that have been predicted using different methods such as gene fusion,26 phylogenetic profiles,27 or gene neighborhood.28 For each of these predictions, STRING gives a confidence score, which measures how likely it is that the prediction is a true interaction. Furthermore, we have developed an algorithm that predicts interactions between proteins by propagating known interactions in one species to ortholog proteins in another species (ie, assuming protein-protein interologs).29 For this algorithm we have used the protein-protein interactions contained in the data base of interacting proteins.30
Using interactions from STRING (at a confidence level of 400) and the protein-protein interologs found by our algorithm, we built a protein-protein interaction network for each of the metastatic organs studied in this work. Taking as root proteins of the network those that in the metastatic organ were under- or overexpressed with respect to the tumor, we obtained networks that attempt to describe the relationship between proteins relevant to metastasis in that organ. We then grouped proteins by their cellular function and extracted those functions that were connecting at least two root proteins. These cellular functions and the root proteins that they connected provided us with clues about processes that are potentially involved in metastasis.
| Results |
|---|
|
|
|---|
For these experiments we used MDA-MB 435 cells transfected with the anti-apoptotic gene Bcl-xL or the control 435/Neo cells. The parental cells expressed low levels of Bcl-xL protein, and the transgene protein expression was twofold higher than in the control cells. Comparative analyses of Bcl-xL protein expression by Western blot showed stable overexpression of these genes in 435/Neo and 435/Bcl-xL tumor cells (Figure 1)
.
|
|
Identification of Bcl-xL-Mediated Protein Expression Changes in Metastasis
Seventeen proteins displaying differential spots in 2DE gels from either tumor cells or metastases (Figure 3)
could be identified by MALDI-TOF peptide fingerprint analysis, from a total of 29 picked out. These proteins are listed in Table 1
, with their corresponding MW, pI, and recognized function (according to SwissProt database), as well as references to their relevance to breast cancer, primarily derived from experiments using array technologies.
|
|
-chain (twofold), tyrosine 3/tryptophan 5-monoxygenase activation protein (48-fold), ATP synthase ß-chain (twofold), My 032 protein (twofold), stathmin 1 (fourfold), familial ALS mutant G37r (fourfold), nonmetastatic cells 1 protein nm-23 (fivefold), peroxiredoxin 2 (threefold), glutathione transferase (twofold), aminoacylase 1 (fivefold), peroxiredoxin 3 (threefold), and enolase 1 (sixfold).
|
-chain (threefold), tyrosine 3/tryptophan 5-monoxygenase activation protein (33-fold), ATP synthase ß-chain (threefold), familial ALS mutant G37r (fourfold), aminoacylase 1 (sixfold), peroxiredoxin 3 (twofold), and enolase 1 (twofold). Only three proteins were overexpressed in bone metastasis: multicatalytic endopeptidase complex
-chain (threefold), tyrosine 3/tryptophan 5-monoxygenase activation protein (20-fold), nonmetastatic cells 1 protein nm-23 (threefold), and enolase 1 (twofold).
Protein underexpression with regard to the tumor was more frequent in bone metastasis. We found less ß-tubulin co-factor A (0.01-fold), F-actin capping protein
-1 subunit (0.3-fold), hypothetical protein MGC15429 (0.1-fold), and carbonyl reductase 3 (0.3-fold). In contrast, only Ran BP1 protein was underexpressed (0.4-fold) in 435/Bcl-xL lung metastasis. In lymph node metastasis we did not find underexpression of these proteins relative to the tumor. Thus, the phenotype evolving in lung metastasis was quite different from that of bone metastasis.
In contrast to 435/Bcl-xL lung metastasis, neither of the above described changes were found in lung metastasis from 435/Neo cells, where the only observed modifications, with regard to its tumor, were overexpression of F-actin capping protein
-1 subunit (twofold), and underexpression of tyrosine 3/tryptophan 5-monoxygenase activation protein (0.4-fold). Furthermore, neither of the 435/Bcl-xL metastasis-related proteome alterations were found in 435/Bcl-xL tumor cells with regard to 435/Neo tumor cells, thus expression changes induced by Bcl-xL must be emerging during tumor development and modulated throughout metastatic steps.
Survival in the Blood Stream (BSC) versus Organ-Specific Changes
Because protein expression profiles varied between metastatic sites, we analyzed the expression of proteins in cells living in bloodstream and rescued from lung 4 hours after intravenous injection.18 This allowed us to discern whether changes in protein expression resulted from cross-talk between metastatic cells and the different environments that finally select the phenotype, or belonged to subpopulations living in tumors selected during progression to metastasis.
The expression of aminoacylase 1 and enolase 1 in BSCs with regard to 435/Bcl-xL tumor cells increased by 9- and 11-fold, respectively. Moreover, no over- or underexpression was found in the rest of proteins. Thus, the enrichment of both proteins is a phenotype selected in the bloodstream and maintained in the following steps of the metastatic process.
Some organ-specific changes were already present in BSCs, such as underexpression of ß-tubulin co-factor A (0.03-fold), hypothetical protein MGC15429 (0.4-fold), and carbonyl reductase 3 (0.3-fold), present in bone metastasis; or Ran BP1 (0.04-fold) present only in lung metastasis. Conversely, other organ-specific proteins were not enough represented in BSCs but did occur in metastatic foci, such as My 032 protein and stathmin 1 (twofold and fourfold, respectively) in lung cells, and F-actin capping protein
-1 subunit (0.3-fold) in bone cells.
Finally, changes common to all metastases but not evident in BSCs were multicatalytic endopeptidase complex
-chain and tyrosine 3/tryptophan 5-monoxygenase activation protein (14-3-3-
). We used Western blot analyses with specific antibodies against some of these proteins to assess with an alternative technique the different expression between tumor, BSCs, and metastases (Figure 4A)
. The analyses confirmed the overexpression (more than twofold) of enolase 1 and glutathione transferase in BSCs, and in all metastases from 435/Bcl-xL cells with regard to the tumor; the increased expression of peroxiredoxin 2 and stathmin 1 in 435/Bcl-xL lung metastasis with regard to tumor and BSCs; and the overexpression of 14-3-3-
in all metastasis with regard to the tumor and BSCs (Table 3)
.
|
|
An in Silico Proposal of Interacting Proteins
The bioinformatic analysis revealed several potential human-protein interactors within the 17 proteins (root proteins hereafter) for which an organ-specific change in expression was experimentally found. These protein-connectors were clustered according to functional criteria: 1) cell division; 2) DNA-repair, synthesis, and translation; 3) redox; 4) structure and motion (ie, dynein, tubulin, actin, etc); 5) ATP synthase; 6) kinase-phosphatase; 7) calcium metabolism; 8) protein folding stability; 9) Krebs cycle; and 10) notch signaling. Proteins for which no function or human homologue was found were ignored.
The largest set of putative interactors was found linking proteins involved in lung metastasis, which could be expected from the fact that most root proteins exhibited differential expression in lung metastasis. Table 4
lists proteins bridging two or more proteins with organ-specific changes in expression. For all organs a central connection with a proteasome subunit, P28066, was found. It is also noteworthy that calmodulin (P62158), involved in calcium metabolism, connects P28066 and NDP kinase A, a metastasis-associated protein, P15531.
|
In lung metastasis there are more proteins linking the two types of peroxiredoxin P32119 and P30048. Among them it is worth to point out: 1) O15160, another peroxiredoxin; 2) Q8N9M6, assigned to cell division function (see above); and 3) P40818 (ubiquitin thiolesterase 8) and Q86YA5 (SMARCA3), two proteins involved in ubiquitination and containing domains for DNA-binding, helicase C activity, and a ring zinc finger. Among the proteins found to link proteins with altered expression in lung metastasis, many are kinases, usually related to signal cascades: Q9UBV2, a homologue of Sel-1; HIPK2 (Q9H2X6), which phosphorylates homeodomain transcription factors; Rab22 (Q9UL26), in the GTPase family; and Cdc14A3 phosphatase (O60728), which is a Tyr-phosphatase that may play a role in cell division.
The functions exerted by metastatic cells have been validated in silico. It has to be understood that these links are proposed through predictions based on homology because no real experiments have been performed with these proposed human proteins. Therefore, only the implication on a cluster of related proteins should be considered as potential real connections. This limits the prediction to few clusters: for proteins involved in a redox system (putatively affecting other functions involved in DNA helicase and cell division and ubiquitination), the set of kinases (which putatively could implicate some signal transduction), and the ATP synthase complex. Taking into account known proteins and their interactions, we suggest a model of breast cancer metastasis progression of Bcl-xL-overexpressing tumors, in which Bcl-xL could mediate changes in metabolic pathways of breast cancer cells, with the result of ubiquity and organ specificity of a progressively selected phenotype from the survival in the blood stream to organ-specific clinical metastasis (Figure 5)
.
|
| Discussion |
|---|
|
|
|---|
Proteomic profiling of metastatic human breast cancer cells by 2DE and mass spectrometry identified a Bcl-xL-mediated phenotype that harbors cells from the primary tumor to the metastatic state and appears to be useful for shortening dormancy with the result of clinical metastasis in many organs. Several proteins that function in redox system are overexpressed (familial Als mutant G37r, peroxiredoxin 2, peroxiredoxin 3) or underexpressed (carbonyl reductase 3), particularly in lung metastasis from 435/Bcl-xL cells. The presence of these proteins has been reported in some tumor cells and tissues,34-36 validating their implications in cancer. Furthermore, the protein-protein interaction program shows a common core related to the redox system, involving peroxiredoxins and carbonyl reductase, that connects several structural proteins, proteins implicated in ubiquitination, in DNA binding, cell division, kinase-phosphatase activity, and ATP synthase complex. These results suggest that cellular response to oxidative stress is a mechanism implicated in the selective growth of cells at the metastatic foci.
Cells living under aerobic conditions, as might be the case of lung metastasis, need to protect themselves against the damage caused by reactive oxygen species, arising from either the incomplete reduction of oxygen or exposure to external agents.37 Peroxiredoxins are mitochondrial or cytoplasmic enzymes that catalyze the destruction of peroxides and use cellular thiol to regenerate their active state.38 It has been reported that peroxiredoxin 3 protects cancer cells against hypoxia and drug-induced hydrogen peroxide-dependent apoptosis.39 Experiments using mitochondria-specific fluorescent probes demonstrate that peroxiredoxin 3 is essential for maintaining mitochondrial mass and membrane potential in transformed rat and human cells.40
The protein-protein interaction program showed a relationship between enolase 1 and aminoacylase 1 proteins, both cytoplasmic, with a high reliability score, suggesting that glucose and amino group metabolism could be associated to metastatic activity in terms of functional gain. Furthermore, the interaction of both with redox system proteins suggests that glycolysis metabolism and antioxidant pathways might be connected to the achievement of metastatic activity.
Enolase 1 is an enzyme that converts 2-phospho-D-glycerate to phosphoenol-pyruvate in the glycolysis pathway. The gene
enolase encodes a nuclear shorter monomeric structural protein by an alternative translation start, which has been found to bind to an element in the c-myc promoter. It is expressed in mammary gland and in breast carcinoma with endocrinal cells.41
Aminoacylase 1 is a hydrolase involved in the urea cycle and metabolism of amino groups. Mutational inactivation has been found in small cell lung cancer.42
The presence of both genes has also been established in array studies in which tumor cell lines and tissues, including breast, were analyzed and different expression patterns among samples were found.34-36
Expression of Bcl-xL influences a signaling network that links changes in mitochondria metabolism to alterations in nuclear gene expression, even under nonapoptotic conditions.43 Thus, genetic predisposition at primary tumor must exist because tumor cells stably overexpressed Bcl-xL, but a higher restriction takes place in the blood stream, ie, enolase 1 and aminoacylase 1 are overexpressed in BSCs with regard to breast tumor, and they are further represented in metastases. Because Bcl-xL inhibits anoikis and increases the survival of BSCs,18 overexpression of these proteins might be a result of a selection process. Indeed, the role of Bcl-xL in metastases might be associated to its ability to prevent apoptosis by coupling oxidative phosphorylation and glycolysis.44
Some phenotypic changes appeared only in metastatic foci, such as, tyrosine 3/tryptophan 5-monoxygenase activation protein, which activate tyrosine and tryptophan hydroxylases in the presence of calmodulin protein kinases; and multicatalytic endopeptidase complex
-chain, a proteasome subunit, characterized by its ability to cleave peptides in an ATP/ubiquitin-dependent nonlysosomal proteolytic pathways.45
Because both are overexpressed in all metastases, these proteins could have a role in common metabolic pathways involved in shortening dormancy.
Tyrosine 3/tryptophan 5-monoxygenase activation protein belongs to the 14-3-3 family of cytoplasmic proteins that mediate signal transduction binding to phosphoserine-containing proteins. A down-regulation has been reported in breast cancer cell lines MCF-7 and MDA-MB-231, and in primary breast carcinomas as compared with normal breast epithelial cells.46 In agreement with these results, we found underexpression of this protein in 435/Bcl-xL tumor cells with regard to control tumor. Moreover, a higher expression occurred in metastases, thus emphasizing the importance of the cross-talk between metastatic cells and the microenvironment in the development of clinical metastasis. Indeed, 14-3-3 proteins prevent apoptosis by directly binding to proapoptotic proteins as Bax and Bad.47 Moreover, it interacts with CDC25 phosphatases, RAF1 and IRS1 proteins, suggesting its role in diverse biochemical activities related to cell division.45
The organ-specific phenotype of tumor cells might be as a result of the selection in both the blood stream and the lodging site: differential expression of hypothetical protein MGC15429 and ß-tubulin co-factor A underexpressed in bone metastasis, and Ran BP1, a GTP-binding nuclear protein, underexpressed in lung metastasis, are still represented in blood cells with regard to the tumor. Moreover, My 032 protein, glutathione transferase, and stathmin 1 might be a phenotype positively selected in the lodging site.
Stathmin 1, an ubiquitous cytosolic phosphoprotein that binds to two
/ß-tubulin heterodimers, functions as an intracellular relay integrating regulatory signals of the cellular environment.48,49
Its overexpression has been found in highly proliferative breast carcinomas.50,51
Furthermore, glutathione transferase belongs to an enzyme family that plays an important role in detoxification by catalyzing the conjugation of many hydrophobic and electrophilic compounds with reduced glutathione, with a prognostic significance on invasive breast cancer.52
Increased expression of stathmin, enolase, peroxiredoxin 2, and some 14-3-3 protein isoforms in human breast ductal carcinoma in situ with regard to the normal breast has been described.53
Moreover, a serological proteomic study has identified a 14-3-3 protein
as a breast cancer biomarker.54
These proteins are not included in the protein expression map of the normal human breast cells, suggesting that they belong to the earliest stage of breast cancer progression.55
The prevailing paradigm of carcinogenesis suggests that epithelial cells sequentially accumulate multiple genetic and epigenetic changes. It has been shown that cells isolated from metastases are frequently more highly metastatic than the bulk population of cells from primary tumors.56 Progenitors of the later-arising metastasis must be present among those cells that have disseminated to distant sites before the primary tumor is removed, and metastatic variants can arise by genetic or epigenetic mechanisms. Therefore, a signature of metastatic genes can exist.11-13 Moreover, tumor cells may disseminate in a far less progressed phenotype and acquire the metabolic properties of cells that live under the stress of no adhesion conditions, and the metabolic properties of cells with organ-specific growth, during the eventual development of metastatic disease.10 Because disseminate cells progress independently from the primary tumor, a simple extrapolation from primary tumor data to disseminated cancer cells is impossible.
Up-regulation of Bcl-xL protein may be a marker of tumor progression. High levels of Bcl-xL expression in breast tumors correlate with poor patient survival and increased metastasis.57 We suggest that overexpression of Bcl-xL in breast cancer cells act as an adaptor protein mediating breast cancer progression by connecting several metabolic pathways with the result of a phenotype that includes several functional pathways useful to metastasis development. Because metastasis is a state that emerges from the tumor-host interface, in which there is an interactive dialog between the neoplastic cells and the local activated host cells,58 the cross-talk results in different facilitated functions in each microenvironment (metastatic foci, cells that survive in blood, and tumor cell populations), arguing for the importance of anti-apoptotic proteins in paracrine and autocrine signaling as a crucial feature to therapeutic strategies to prevent metastases.
| Footnotes |
|---|
Supported by the Ministerio de Sanidad y Consumo (grants FIS 01/1469 and FIS/PI041937), the EC MetaBre (contract no. LSHC-CT-2004-506049), the Spanish Ministerio de Ciencia y Tecnologia (MCyT BIO02-0369 to R.A. and B.O.), and the Fundación Ramón Areces (to R.A. and B.O.).
L.E. and B.M. contributed equally in this project.
Accepted for publication June 23, 2005.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
O. Mendez, B. Martin, R. Sanz, R. Aragues, V. Moreno, B. Oliva, V. Stresing, and A. Sierra Underexpression of transcriptional regulators is common in metastatic breast cancer cells overexpressing Bcl-xL Carcinogenesis, June 1, 2006; 27(6): 1169 - 1179. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Aragues, D. Jaeggi, and B. Oliva PIANA: protein interactions and network analysis Bioinformatics, April 15, 2006; 22(8): 1015 - 1017. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE |