|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular & Cellular Proteomics 7:394-410, 2008.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
TCDD is a ligand of the aryl hydrocarbon receptor (AhR), and it is assumed that most, if not all, toxic effects of TCDD are mediated by the AhR (4). Binding of TCDD to the cytosolic receptor results in the translocation of the ligand receptor complex into the cell nucleus and the dimerization with a related protein, the AhR nuclear translocator (ARNT). This heterodimer can act as a ligand-activated transcription factor by binding to a specific DNA sequence known as aryl hydrocarbon response element I (AHRE-I), xenobiotic-responsive element, or dioxin-responsive element, an enhancer located in regulatory regions of aryl hydrocarbon-responsive genes. Binding to the AHRE-I activates the transcription of target genes and thus results in an alteration of gene expression (5). In addition, the liganded AhR/ARNT heterodimer can function as a ligand-activated coactivator of the transcription of a different set of responsive genes via the interaction with another response element called AHRE-II (6, 7). TCDD-induced changes in gene expression may also occur indirectly as a consequence of an AhR/ARNT-dependent induction of oxidative stress and the transcriptional activation of still another set of genes via a different enhancer, the so-called antioxidant response element (ARE) (8). Moreover several AhR-dependent mechanisms not involving heterodimer binding to an AHRE have been suggested to contribute to the toxicity of TCDD (9–11). Thus, the mechanisms by which TCDD may interfere with signaling pathways appear to be extremely complex, and the elucidation of the resulting alterations in gene expression and their importance for the various aspects of TCDD toxicity clearly remains a continuing challenge.
The alterations in gene expression at the level of the transcriptome in hepatic cells exposed to TCDD have been explored in a number of studies on hepatoma cell lines from different species (12–14) and on rodents (7, 15–20). These studies have provided a large body of information on the effects of TCDD on gene expression at the mRNA level. There are, however, well documented, strong limitations regarding the correlation between mRNA and protein levels, and for the vast majority of the observed changes in mRNA abundance it has remained open to question whether they are of biological relevance because it is not known whether they are actually translated into changes of the abundance of their cognate proteins, i.e. at the functional level.
Studies on the effects of TCDD on hepatic cells of mammals at the level of the proteome are still scarce. Three in vivo investigations have dealt with the effects of dioxin on the protein profile of rats (21, 22) and marmosets (23). Although a comparison of the obtained results is difficult or even impossible because of the use of, in part, different species, rat strains, application routes, doses, and exposure periods, it is evident that relatively low numbers (up to a maximum of 10) of up- or down-regulated proteins were detected in the individual studies. This might reflect a limited biological activity of TCDD in vivo under the experimental conditions used, but it appears more likely that the low numbers were due to technical limitations of the in vivo approach. A major limiting factor of the studies, which all used two-dimensional gel electrophoresis (2-DE) for protein separation, may have been the use of IPG strips covering the broad pH range from 3 to 10 in all of these investigations. Unfortunately the use of these low resolution strips is almost mandatory for in vivo studies because the necessity to use reasonably large groups of animals to account for the unavoidable variability between individual animals practically precludes the use of combinations of the much better resolving IPG strips covering overlapping, narrow pH ranges (24, 25).
With the objective to obtain detailed insight into the proteome changes induced by TCDD in hepatic cells and to obtain new clues to the signaling pathways mediating its cellular activity, we are currently performing a series of in vitro studies aimed at characterizing the changes in protein abundance and posttranslational protein modifications in TCDD-treated 5L rat hepatoma cells. 5L cells are epithelial-like cells that express both the AhR and the ARNT protein. They are highly responsive to the toxic effects of TCDD, and they have been used in investigations on the role of AhR and ARNT in TCDD toxicity, cell cycle regulation, and cell signaling (11, 26–33) and on the induction of xenobiotic-metabolizing enzymes by TCDD (34). Moreover they have been successfully used in projects aiming at the identification of novel target genes of the AhR (14, 35). We recently reported results of a proteomics study in which the effects of TCDD on 5L cells were investigated using an LC-based mass spectrometric approach, the so-called isotope-coded protein label (ICPL) method (36, 37). In this study, we identified 89 protein species of various functional categories as up- or down-regulated by TCDD. It was likely, however, that the identified proteins represented only a fraction of the proteins actually altered by TCDD because several studies dealing with the detection of differentially expressed proteins by diverse techniques for quantitative proteome analysis reported a very limited overlap of the proteins identified by the different approaches (38, 39). It thus appeared expedient to complement our findings with results obtained by an experimental system for proteome analysis unrelated to the ICPL approach.
We have therefore conducted a quantitative proteomics analysis of the proteome alterations induced by TCDD in 5L cells using a gel-based approach, the classical 2-DE for protein separation and quantitation in combination with mass spectrometry for protein identification. To optimize protein resolution, separation of the proteins in the first dimension, i.e. by IEF, was performed on four individual gel strips with different narrow, overlapping IPGs covering a total pH range of 3–11. This approach has been shown to considerably increase the number of proteins detectable in comparison with the IPG strips of pH 3–10 or 4–7 most frequently used in 2-DE analyses and, thus, to yield a much more detailed picture of the proteome under study (24, 25).
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Cells
The 5L cell line is a variant of the Fao line (40), a descendant of the cell line H4IIEC3 established by Pitot et al. (41) from the Reuber H35 rat hepatoma (42). BP8 cells represent a variant clone of 5L cells that lacks the AhR (26, 29). Cells were grown as monolayers in 10-cm culture dishes containing RPMI 1640 medium supplemented with 10% FCS, 100 units/ml penicillin, and 100 µg/ml streptomycin ("complete RPMI 1640 medium") at 37 °C in a humidified atmosphere of 95% air and 5% CO2.
Treatment of the Cells with TCDD
Cells were seeded at a density of 1.5 x 106 cells/10-cm dish in complete RPMI 1640 medium. 24 h later, the medium was replaced by 5 ml of fresh complete RPMI 1640 medium, and 48 h after seeding 1 nM TCCD dissolved in Me2SO or solvent alone was added for the times indicated. The Me2SO concentration in the medium was 0.1%.
Harvest and Solubilization of Cells for 2-DE
After incubation, the medium was removed, and cells were washed three times with ice-cold PBS without Ca2+/Mg2+. To detach the cells from the plates, the monolayers were exposed to 3 ml of trypsin/EDTA at room temperature, the trypsin was removed completely, and the cells were detached by vigorous tapping of the plates on the table. The cells were rinsed off the plates with 3 ml of PBS-Ca2+/Mg2+ (three times) and pelleted by centrifugation (0 °C for 5 min at 1000 x g). The cell pellet was carefully resuspended in 1 ml of 0.5x PBS without Ca2+/Mg2+ (0 °C) and centrifuged again (0 °C for 5 min at 1000 rpm). The cell pellet was then covered with 500 µl of solubilization buffer at room temperature and immediately vortexed until the cells had dissolved completely. Solubilized cells from four identically treated dishes were each combined. The solubilization buffer contained 7 M urea, 2 M thiourea, 2% (w/v) CHAPS, 65 mM DTE for IPG strips pH 3.0–5.6 and 5.3–6.5 or 15 mM DTE for pH 6.2–7.5 and 7–11 (see below), 0.4% (w/v) Biolyte pH 3–10 (Bio-Rad), 0.4% (w/v) Pharmalyte pH 3–10 (Amersham Biosciences), and one protease inhibitor tablet (Complete, Roche Applied Science)/50 ml of buffer. Before use, a stock solution of urea and thiourea was treated with a mixed-bed ion exchanger (1% (w/v) Serdolit MB-1; Serva, Heidelberg, Germany) to remove ionic impurities.
Two-dimensional Gel Electrophoresis
2-DE was performed essentially as described by Görg et al. (25). In every experiment four samples for each pH range and treatment condition (TCDD or Me2SO) were processed in parallel. Proteins were separated in the first dimension by IEF in a Multiphor II or IPGphor apparatus (Amersham Biosciences) using four different types of 24-cm IPG strips (Amersham Biosciences) covering pH ranges of 3.0–5.6, 5.3–6.5, 6.2–7.5, and 7–11. The strips were rehydrated in rehydration cassettes (Amersham Biosciences) with buffer containing 6 M urea and 2 M thiourea (pretreated with 1% (w/v) Serdolit MB-1), 1% (v/v) CHAPS, the appropriate IPG buffer (0.2%, w/v), a trace of bromphenol blue, and, depending on the pH range of the strips, 15 mM DTE for pH ranges 3.0–5.6 and 5.3–6.5 or 15 mg/ml DeStreak (Amersham Biosciences) for pH ranges pH 6.2–7.5 and 7–11. The IPG buffers used were 3.5–5.0 for pH range 3.0–5.6, 5.5–6.7 for pH range 5.3–6.6, 6–11 for pH range 6.2–7.5, and 7–11 for pH range 7–11. Following overnight rehydration of the strips, they were placed in the Multiphor II or IPGphor apparatus, and protein samples (100 µl) were applied to the IPG strips by cup loading at the anodic end (IPG strips pH 5.3–6.5, 6.2–7.5, and 7–11) or the cathodic end (pH 3.0–5.6). 250 µg of protein were applied with the exception of the pH 7–11 strips, which were run with 100 µg of protein. IPG strips were subjected to IEF at 20 °C for 95 kV-h (pH 3–5.6), 150 kV-h (pH 5.3–6.5 and pH 6.2–7.5), and 100 kV-h (pH 7–11). The focusing protocols included a change of the electrode strips at low voltage to remove most of the salts in the samples.
For reduction and alkylation of the proteins, the IPG strips were equilibrated for 15 min in a solution of 50 mM Tris-HCl (pH 6.8), 6.0 M urea, 30% (v/v) glycerol, and 2% (w/v) SDS containing 1% (w/v) DTE and then for an additional 15 min in the same solution except that DTE was replaced by 4% (w/v) iodoacetamide. After equilibration, proteins were separated on SDS-polyacrylamide gels (12% T, 2.6% C) at 20 °C using the Ettan DALT II separation unit (Amersham Biosciences) (3 watts/gel for the first 30 min and 18 watts/gel until the end). The running buffer consisted of 50 mM Tris-HCl (pH 8.8), 386 mM glycine, and 2% (w/v) SDS. The Precision Plus Protein Standard (Bio-Rad) was used as molecular weight marker.
Protein Staining and 2-DE Gel Image Analysis
Gels were stained with the fluorescent dye ruthenium(II)-tris(bathophenanthrolinedisulfonate) as described by Lamanda et al. (43). Stained gels were scanned with the Fuji FLA-3000 phosphorimaging system (Raytest) using the software BASReader Version 3.01 with the following parameters: resolution, 100 µm; 16-bit image (65,536 gray levels); sensitivity, 100; excitation wavelength, 473 nm; detection filter, O580.
For detection and quantitation of protein spots and comparison of gels, Proteomweaver software (release 4.0.0.5 beta, Bio-Rad) was used. The following parameters were used for spot detection: minimum spot radius limit of 4, minimum spot intensity (volume above base level) of 2000, and minimum contrast limit (height above base level) of 40. Gels from each experiment were processed by pair-match-based normalization. Subsequently protein spots of every gel of both groups were automatically matched with each other, and the matching was carefully checked manually and corrected in case of obvious mismatches to avoid assigning false positives. To identify significant differences in spot intensities between the two groups the following rules were established: only spots matched in at least 50% of the gel images in a group were considered, and to get reasonable protein amounts for mass spectrometric identification only spots exceeding an average intensity threshold of 0.2 were taken for further analysis. The remaining spots were sorted according to their regulation factors, and for each experimental set the threshold regulation factor for the significance level (p < 0.05) was determined by the software using Student's t test. Differences were considered significant at p < 0.05 corresponding to a thresholding for 1.54-fold changes in spot intensities (i.e. 1.54 for up-regulation and 0.65 for down-regulation). Only spots exhibiting intensities significantly changed after TCDD treatment were further considered as candidate spots and subsequently subjected to manual verification.
Protein Identification by Mass Spectrometry
In-gel Proteolysis—
Protein spots with altered intensities were selected for identification and, following silver staining of the 2-DE gels using a modification of the mass spectrometry-compatible version of Shevchenko et al. (44), manually excised with a scalpel. The gel pieces were washed with buffer containing 50 mM NH4HCO3 in 30% ACN, destained with 30 mM potassium hexacyanoferrate(III) and 100 mM sodium thiosulfate, and washed again with 50 mM NH4HCO3 in 30% ACN. The pieces were shrunk with 100% (v/v) ACN prior to proteolytic digestion. Following addition of 0.2–0.4 mg of modified trypsin (Promega, Madison, WI), proteolysis was performed overnight at 37 °C. The supernatant was collected and combined with the eluates of two subsequent elution steps with 80% (v/v) ACN and 1% (v/v) TFA. The combined eluates were dried in a SpeedVac centrifuge and dissolved in 5–10 µl of 70% (v/v) ACN and 0.1% (v/v) TFA. Samples that did not yield successful protein identifications in the subsequent MS analysis were reanalyzed after desalting and concentrating peptides with ZipTip® pipette tips containing C18 reversed-phase medium (Millipore, Billerica, MA).
MS Analysis—
Mass spectra were acquired using a Proteomics Analyzer 4700 (MALDI-TOF/TOF) mass spectrometer (Applied Biosystems, Framingham, MA). Measurements were performed with a 355 nm neodymium-doped yttrium aluminium garnet (Nd:YAG) laser in positive reflector mode with a 20-kV acceleration voltage. For each MS and MS/MS spectrum, 3000 shots in a m/z window of 800–4000 were accumulated. To obtain a mass accuracy better than 50 ppm, internal calibration using peptides resulting from autoproteolysis of trypsin (2163.05 and 2211.10 Da) was utilized. In addition, external calibration in the mass range from 800 to 3000 with a peptide mass standard consisting of eight peptides (Bruker, Billerica, MA) was used. For sample preparation, 0.5-µl aliquots of the digests were spotted on a target, dried, and mixed 1:1 with matrix consisting of 2 mg/ml
-cyano-4-hydroxycinnamic acid in 70% (v/v) ACN and 0.1% (v/v) TFA. Peak labeling was automatically done by 4000 series Explorer software Version 3.6 (Applied Biosystems) without any kind of smoothing of the peaks or base line. In addition to peptide mass fingerprint spectra, the five most abundant precursor ions having a signal-to-noise ratio higher than 70 were chosen for MS/MS fragmentation, which was performed using medium collision energy. An exclusion list of known contaminant ion masses of keratin and trypsin was used.
Data Processing—
The processing of spectra obtained with the 4700 Proteomics Analyzer was carried out with the software GPS ExplorerTM Version 3.6 (Applied Biosystems). The combined MS and MS/MS spectra were searched against the Swiss-Prot (230,133 sequences, 84,471,903 residues; date July, 25, 2006) and the National Center for Biotechnology Information non-redundant (NCBInr) databases (4,878,246 sequences, 1,686,729,293 residues; date April 26, 2007) using an in-house version of the Mascot (Version 2.0) search engine (Matrix Science Ltd., London, UK) with the following parameters. One miscleavage was allowed. As taxon we chose Rattus, and as enzyme we chose trypsin. Carbamidomethylation of cysteine was set as fixed modification, and oxidation of methionine was set as variable modification. For MS spectra, mass errors of 65 ppm for precursor ions and 0.3 Da for fragment ions were allowed. For MS/MS of peptides, further filtering criteria were applied: the maximum peptide rank was set at 3, and the minimum ion score was set at 30% confidence level, which corresponds to a Mascot score of 10–12.
Criteria for Protein Identification—
The GPS Explorer 3.6 software reports three different scores: The Mascot best ion score, i.e. the highest score of a single peptide; a total ion score, i.e. the sum of all peptide scores of one protein; and a protein score that includes a score for the peptide mass fingerprint in addition to the total ion score. The significance level is usually >20 for a peptide score and >40–50 for a protein score. Because different database searches have different Mascot significance levels due to different databases sizes and different numbers of masses submitted for a search, scores cannot be compared directly. For this reason, the software calculates a confidence interval from Mascot protein scores or ion scores, and the Mascot significance level for each search is defined as the 95% confidence level. Therefore, the total protein score confidence level is a reliable and comparable parameter for the significance of a database search.
A protein was regarded as identified if the following three criteria were fulfilled. (i) The probability-based MOWSE (molecular weight search) protein score was above the 5% significance threshold for the respective database (95% confidence level). (ii) The matched peptide masses were abundant in the spectrum. (iii) The theoretical pI and the molecular mass of the protein suggested by the search result could be correlated with the position of the corresponding spot in the 2-DE gel. If the latter was not the case, posttranslational modifications (PTMs) were taken into consideration. In most cases only rank 1 proteins were identified for each 2-DE spot as the peptides of the proteins of the following rank positions were carefully searched for peptides shared with the identified protein. Additional protein hits were only included in the protein list when a significant score of the non-shared peptides was obtained. In the case where all peptides detected were shared by individual proteins, as was the case for e.g. myosin regulatory light chain 2-B and 2-A isoforms, this ambiguity was clearly marked in the tables.
A summary table is available as supplemental Table S1. MS/MS data for each peptide match including precursor masses with observed mass errors, peptide sequences, ion scores, ion score confidence levels, and modifications are listed. Only MS/MS data of peptides with an ion score confidence level higher than 95% are shown.
Western Blot Analysis
For Western blot analysis of VDAC2 abundance, cells were lysed in Cell Extraction Buffer (BIOSOURCE, Camarillo, CA), and 10 µg of total protein were resolved on 12% polyacrylamide gels (Protean II XL, Bio-Rad) and blotted onto PVDF membranes. VDAC2 was detected using polyclonal anti-VDAC2 antibodies (ab22170 or ab37985, Abcam, Cambridge, UK), a horseradish peroxidase-labeled secondary antibody, and ECL detection kit (Amersham Biosciences). The actin band, which served as loading control, was detected using an anti-actin antibody (clone AC-40, Sigma).
Real Time Quantitative PCR
The expression of VDAC2 and cytochrome P4501A1 (CYP1A1) at the mRNA level was quantitated by real time quantitative PCR with a LightCycler instrument (Roche Applied Science) using the following primers: VDAC2: forward, 5'-TTGGAGTGGGCTATACTCAGACT-3'; reverse, 5'-GCATTAAAGCTCTTCCCGTCT-3'; CYP1A1: forward, 5'-CCTTGGAGCTGGGTTTG-3'; reverse, 5'-GCTGTGGGGGATGGTGAA-3'.
Measurement of the Generation of Reactive Oxygen Species
The formation of reactive oxygen species (ROS) in 5L cells was determined by measuring the intracellular oxidation of dihydrorhodamine 123 (DHR123) (Invitrogen-Molecular Probes) to the fluorescent rhodamine 123 (R123) by flow cytometry. Cells were exposed to 1 nM TCDD for 8 h as described above or to 50 µM menadione for 30 min, and 15 min before the end of the exposure DHR123 was added to a final concentration of 5 µM. Cells were rinsed with PBS, harvested with trypsin, and immediately placed on ice for flow cytometry. 5000 cells were analyzed on an LSR II flow cytometer (BD Biosciences) using excitation at 488 nm and a 530/30-nm band pass filter for detection. Autofluorescence was determined on 5L cells incubated without the dye and used for background subtraction. Populations of viable single cells (based on forward sideward scatter) were gated identically for ROS determinations for each treatment group. Mean fluorescence values were acquired using BD FACSDiva software and used for quantitative analysis. Histogram overlays were performed using WinMDI (Scripps Research Institute). All quantitative data were analyzed by single factor analysis of variance (Microsoft Excel) test for significance (p < 0.05).
| RESULTS |
|---|
|
|
|---|
TCDD-induced Proteome Alterations
Lysates of 5L cells exposed to 1 nM TCDD in Me2SO or to solvent alone were separated by 2-DE using four different narrow, overlapping 24-cm IPG strips covering a total pH range from 3 to 11 for the IEF step. Following fluorescence staining and quantitative comparison of the gels obtained from treated and untreated cells, gels were silver-stained, and protein spots that had exhibited a statistically significant, reproducible alteration by TCDD were cut out, digested with trypsin, analyzed by MALDI-TOF/TOF mass spectrometry, and identified by database search. Fig. 1 shows representative examples of the fluorescence-stained gels obtained for the four different pH ranges.
|
4200 spots that resulted from about 830, 1780, 1070, and 520 spots for the pH ranges 3–5.6, 5.3–6.5, 6.2–7.5, and 7–11, respectively (mean values of eight gels each). Using previously reported Proteomweaver software parameters for spot detection (45), these numbers were increased to
8800 spots that resulted from about 1830, 2790, 2620, and 1570 spots for the respective pH ranges. The protein species in 78 of the spots with intensities significantly altered by TCDD were identified by mass spectrometry. These protein species, which corresponded to 73 individual proteins, five of which apparently occurred in different posttranslationally modified forms, are listed in Table I.
|
In Table I, proteins were tentatively categorized according to common functional features. This classification is, however, to a large extent arbitrary because many of these proteins are multitasking proteins playing different roles in different cellular contexts. For the ease of orientation, proteins were consecutively numbered according to the order of their appearance in the table.
Known and Putative New Members of the AhR and ARE Gene Batteries—
Among the proteins identified as up-regulated by TCDD were several well known members of the AhR and ARE gene batteries, namely aldehyde dehydrogenase 3A1 (ALDH3A1, number 42), NAD(P)H:quinone oxidoreductase 1 (NQO1, number 1), and, most likely, glutathione S-transferase A2 (GSTA2, formerly called Ya2) (number 2). GSTA2 was identified in an up-regulated protein spot that may have also contained another class Alpha GST isoenzyme, GSTA3 (formerly Yc1) because no peptides unique to GSTA2 as compared with GSTA3 (and vice versa) were identified (see supplemental Table S1). Although it appears likely that GSTA2 induction was involved in the up-regulation of the GST protein spot, a contribution of GSTA3 cannot be ruled out as in our previous study (37) we had observed a significant up-regulation of peptides characteristic for both GSTA3 and GSTA5 (formerly Yc2). As the rat GSTA5 gene contains an AHRE-I within the first intron that is active in regulating transcription (46), an induction of this GST isoform, which has also been detected as up-regulated in a recent microarray analysis on gene expression in TCDD-treated rats (19), by TCDD would also appear plausible.
In addition to the induction of GSTA2 (or, possibly, GSTA3), we detected an induction of GSTP1 (number 3). GSTP1 was identified in two individual up-regulated protein spots (numbers 3a and 3b) of apparently the same molecular weight but with different pI values that probably represented distinct individual forms of the enzyme carrying different PTMs.
The proteins with increased abundance following TCDD treatment also included several for which there is evidence from other experimental systems that their induction may be controlled via the ARE. These proteins were insulin-degrading enzyme (number 32), transaldolase (number 38), lamin A (number 62), and tubulin β-5 chain (number 63), which were recently identified as regulated in an ARE-dependent manner in mice (47, 48). To the best of our knowledge, they have not been known previously to be affected by TCDD. Moreover we observed that two subunits (β and
, numbers 8 and 9) of the chaperonin containing t-complex polypeptide 1 (chaperonin containing TCP1 (CCT)) were up-regulated by TCDD. As the expression of the
subunit of CCT has been reported to be regulated via the ARE in mouse liver (47) and as it appears likely that the expression of the multiple chaperonin subunits occurs coordinately, it seems possible that the observed up-regulation of the β and
subunits also depends on an ARE-mediated mechanism.
We also observed that TCDD affected the enzyme aflatoxin B1-aldehyde reductase (AFAR1 or AKR7A1) (number 4) in 5L cells. The AFAR1 gene contains numerous ARE-related sequences and also one AHRE-I in its promoter region (49) and thus might be expected to be up-regulated by TCDD. However, we did not observe an up-regulation but a significant down-regulation of an AFAR1 species. Thus, further work is clearly required to characterize the effects of TCDD on the level and PTMs of AFAR1 protein in 5L cells.
Cytoskeletal Proteins—
A large group of proteins with altered abundance or PTMs was comprised of structural proteins. Thus, in addition to several proteins of or associated with the cytoskeleton that we recently identified as altered by TCDD (37) and that recurred in the present study (numbers 50 and 57–60), we observed quantitative changes of protein spots corresponding to several actin-interacting proteins (50) not previously known to be affected by TCDD. They include capping protein (number 51), calponin-3 (number 52), the GTPase septin 2 (number 53), drebrin (number 56), and tubulin β-5 chain (number 62). These alterations coincide with the up-regulation of two subunits of CCT (numbers 8 and 9), a cytosolic molecular chaperone composed of eight subunits that assists in the folding of actin and tubulin and other cytosolic proteins (51). The data support the notion (37) that TCDD exposure results in a highly complex, coordinated reorganization of the cytoskeleton of 5L cells.
Proteins Involved in the Maintenance of Mitochondrial Homeostasis—
Several findings of the present study strongly support our previous observation of a marked effect of TCDD on mitochondrial homeostasis (37). Thus, we observed significant alterations of two components of the respiratory chain, namely a down-regulation of ubiquinol-cytochrome c reductase core protein I (UQCRC1, number 40), a core subunit of mitochondrial respiratory complex III, and an up-regulation of the β chain of the mitochondrial ATP synthase (number 41). Up-regulation of the latter is a characteristic cellular response to certain toxicants that has been attributed to an adaptation of cells to mitochondrial damage induced by ROS (52, 53).
An intriguing finding with regard to the response of mitochondria to dioxin exposure was the observation that TCDD clearly affected the appearance of the voltage-dependent anion channel-selective protein 2 (VDAC2) in 5L cells. VDAC2 functions as a pore in the outer mitochondrial membrane by mediating metabolic exchange between mitochondria and the cytoplasm and acts as a key regulator of mitochondrial BAK-mediated apoptosis (54, 55).
TCDD treatment for 8 h altered the intensities of two closely spaced protein spots, both of which were subsequently identified as VDAC2 (Table I, number 39, and Fig. 2A). Whereas the more acidic spot VDAC2b disappeared almost completely in the TCDD-treated cells, the intensity of its more basic counterpart, VDAC2a, increased by a factor of
2.1 (Fig. 2B). Because a TCDD-induced increase in the abundance of VDAC2 would potentially have important consequences for the regulation of apoptosis in 5L cells, we monitored the abundance of VDAC2 protein by Western blotting over treatment periods of up to 24 h. TCDD treatment clearly caused a marked and time-dependent induction of total VDAC2 protein (Fig. 2C), indicating that the observed increase in VDAC2a was not merely the consequence of a conversion of the more acidic to the more basic VDAC2 form. Thus, it appears that TCDD exposure results in both an altered PTM of VDAC2 and an increase in its abundance.
|
TCDD-induced Oxidative Stress
Transcription of ARE-regulated genes can be activated by both electrophiles and oxidants (56–58). As TCDD does not yield electrophilic metabolites, it appeared possible that the observed increases in the abundance of proteins with ARE-containing gene promoters occurred in response to TCDD-induced oxidative stress. To determine whether TCDD treatment actually induces oxidative stress in 5L cells, we used the redox-active probe DHR123. Cells were treated with 1 nM TCDD in Me2SO or only Me2SO for 8 h and then loaded with DHR123. This compound, an uncharged and nonfluorescent dye, passively diffuses across cell membranes and reacts with ROS, resulting in the formation of cationic fluorescent R123 that can be quantitated by flow cytometry.
Fig. 3A shows that the mean R123 fluorescence in TCDD-treated cells increased by about 165% as compared with that of the Me2SO-treated control. Treatment of the cells with 50 µM menadione (2-methyl-1,4-naphthoquinone), which produces reactive oxygen species by redox cycling and which was used as a positive control, for 30 min increased R123 fluorescence by about 265% (Fig. 3B). These results show that TCDD treatment in fact results in the production of ROS in 5L cells.
|
| DISCUSSION |
|---|
|
|
|---|
In a quantitative proteomics study on the effects of TCDD on the proteins of 5L cells using the ICPL method we have recently identified a large number of proteins with altered abundance or PTMs upon exposure of the cells to the dioxin (37). The majority of these alterations had not been known before. In view of the large number of up- or down-regulated protein species detected we reckoned that an independent, technically unrelated approach for proteome analysis would presumably identify an additional, different set of TCCD-responsive proteins that might give further clues to the largely enigmatic mechanisms of dioxin toxicity. Based on this expectation we performed the present study in which a 2-DE approach was used for proteome analysis. To maximize the resolving power of the gels and, thus, the sensitivity of the analysis, we made use of a set of four overlapping, narrow immobilized pH gradients for the IEF step instead of the single standard IPG strip covering a pH range of 3–10 or 4–7 commonly utilized in 2-DE studies.
In conformance with our previous ICPL study (37), we used an exposure period of 8 h to allow the detection of TCDD-induced early stimulatory effects on gene expression and, at the same time, provide sufficient time to detect a decrease in the abundance of specific proteins due to an inhibition of gene expression or an acceleration of protein degradation. This treatment scheme was also chosen to minimize possible secondary changes, such as those associated with altered cell proliferation or the induction of apoptosis, which have been reported at the transcriptome level after longer incubation periods (59).
The 2-DE approach outlined above resulted in the identification of 73 protein species with TCDD-altered abundance or PTMs. Most of the detected alterations, some of which will be discussed in the following sections, had not been described before and, as we had expected, complemented the results of the ICPL study.
Known and Putative New Members of the AhR and ARE Gene Batteries—
The responsiveness of the 5L cells used in the present study to characteristic dioxin actions was clearly demonstrated by the strong induction of CYP1A1 activity and the up-regulation of several classical members of the AhR gene battery, i.e. ALDH3A1, NQO1, and, most likely, GSTA2. The genes of these enzymes contain AHRE-I in their promoter regions (60–62), and the increased expression of their cognate proteins, like that of CYP1A1, is a hallmark of the action of dioxins and other AhR agonists (63). Despite the massive enhancement of CYP1A1 activity, we detected neither CYP1A1 nor other CYP forms known to be responsive to TCDD, such as CYP1A2 or CYP1B1 (64), among the up-regulated proteins. This was not unexpected, however, because CYPs are integral proteins of the endoplasmic reticulum membrane, and their highly hydrophobic nature apparently impedes their analysis by 2-DE (65, 66). Surprisingly, however, we found CYP3A2 (which is not regarded as a member of the AhR and ARE gene batteries) among the protein spots subjected to protein identification (not shown). Its abundance was reduced by 31% and thus not considered as significantly affected by TCDD. It should yet be noted that the CYP3A2 species, which exhibited the expected molecular mass of about 60 kDa, had a pI of 5.4–5.5, which was much lower than the calculated theoretical pI of 8.9 (Scansite, Massachusetts Institute of Technology). It appears possible that the CYP3A2 spot represented a species with multiple acidic PTMs, such as phosphorylations (67, 68), which may have caused a higher hydrophilicity and its appearance on the two-dimensional gel.
Even though ALDH3A1, NQO1, and several of the TCDD-inducible GST isoforms have been shown or suggested to be regulated in an AHRE-I-dependent manner, it appears unlikely that their induction in 5L cells is mediated solely via AHRE-I. The regulatory regions of the genes of rat ALDH3A1 (22), rat NQO1 (69), and the Alpha class GSTs of the rat, GSTA2 (70) and GSTA5 (71), also contain an enhancer termed ARE that is regulatable by electrophiles and oxidative stress, suggesting that ARE-dependent gene activation may also be involved in the induction of these enzymes in 5L cells. This conjecture is supported by the observed induction of GSTP1. The promoter region of rat GSTP1 does not contain AHRE-I but an enhancer designated glutathione transferase P enhancer I that is closely related to the ARE (72), and the enzyme has recently been shown to be induced via the same mechanism as ARE-dependent proteins in rat liver (73).
The expression of ARE-driven genes is mediated by the transcription factor "nuclear factor E2 p45-related factor 2" (Nrf2) (56, 74). In a microarray study on the effect of the chemopreventive agent 3H-1,2-dithiole-3-thione, which is known to activate the Nrf2 pathway, on hepatic gene expression profiles of nrf2 wild-type and nrf2-deficient mice, Kwak et al. (47) recently identified a large number of genes regulated in an Nrf2-dependent manner. The genes included those of several well known ARE-regulated xenobiotic-metabolizing enzymes like NQO1 and GSTA2 but also a large number of others from functional classes that had not been shown previously to be inducible in an Nrf2/ARE-dependent way. Thimmulappa et al. (48) recently identified several additional genes as regulated by Nrf2 in the small intestine of mice. Among the genes identified by the two groups were several coding for proteins found to be up-regulated by TCDD in our present study, namely lamin A, tubulin β-5 chain, insulin-degrading enzyme, and transaldolase. Indirect evidence (47) suggests that the up-regulation of the β and
subunits of the chaperonin CCT detected in the present study may also occur via the Nrf2/ARE pathway. Thus, substantially more proteins than previously thought are probably induced indirectly by TCDD in an Nrf2/ARE-dependent manner, and in view of the data of Kwak et al. (47) and Thimmulappa et al. (48) it appears highly likely that these proteins actually vastly outnumber those regulated directly by TCDD via the AHREs.
Nrf2 is normally sequestered in the cytoplasm by Keap1, a cysteine-rich protein bound to the actin cytoskeleton. Electrophiles, which alkylate specific cysteine residues of Keap1, or oxidants, which cause the formation of intersubunit disulfide bonds, dissociate Keap1 from Nrf2 resulting in nuclear translocation of Nrf2 and induction of ARE-regulated genes (56–58). The Keap1-Nrf2 complex thus acts as a sensor for electrophilic and oxidative stress. TCDD is not metabolized to electrophilic products but clearly induces the formation of reactive oxygen species in 5L cells as indicated by the intracellular oxidation of the redox-active probe DHR123. We therefore conclude that the observed increases in the abundance of proteins with ARE-containing gene promoters occurred in response to TCDD-induced oxidative stress.
It thus appears likely that the observed induction of the Nrf2/ARE-regulated proteins represents a protective mechanism against the deleterious effects of TCDD-induced ROS and their oxidizing reaction products. GSTs of the Alpha family, for example, exhibit GSH peroxidase activity (75) toward lipid hydroperoxides and use the reactive and toxic lipid peroxidation product 4-hydroxynonenal as a preferred substrate (76). The expression of GSTP1 has long been known to be regulated by the cellular redox status (77) and recently been shown to be inducible in an Nrf2-dependent manner in rat liver (73). In addition to its enzymatic function, the enzyme constitutes a sensor able to transmit redox variations to the apoptosis machinery by modulating the stress kinases pathway by inhibiting c-Jun N-terminal kinase (JNK) (78, 79). The expression of GSTP is inversely correlated with JNK activity, and increased expression of GSTP and inhibition of JNK have been linked to the inhibition of apoptosis (78). On the other hand, oxidative stress has been shown to decrease JNK inhibition by GSTP via an oligomerization of the latter (78), to activate JNK (80), and to promote apoptosis (81). The observation that TCDD-treated 5L cells undergo apoptosis that becomes detectable after about 30 h of treatment2 indicates that the observed induction of GSTP1 by TCCD is not sufficient to prevent apoptosis.
The observed up-regulation of transaldolase, another protein not previously recognized to be affected by TCDD, probably also serves to increase the antioxidant capacity of the cells in response to TCDD-induced oxidative stress. This enzyme is a key player in the regulation of the balance between the reversible non-oxidative branch of the pentose phosphate pathway and the irreversible oxidative branch. The latter is responsible for the generation of NADPH required for the maintenance of GSH in its reduced state and, thus, the protection of cells from reactive oxygen intermediates. Thimmulappa et al. (48) recently identified transaldolase as an Nrf2-regulated gene in the small intestine of mice and showed that it was induced coordinately with two genes of the irreversible pathway, glucose 6-phosphate dehydrogenase and 6-phosphogluconate dehydrogenase, by the chemoprotectant sulforaphene. As both enzymes generate NADPH, the induction of transaldolase would be expected to facilitate the formation of reducing equivalents that, in turn, would maintain the antioxidative power of GSH and also directly contribute to an increased antioxidative capacity of the cells (82). Our observation that TCDD markedly stimulates glucose oxidation via the pentose phosphate pathway in 5L cells3 suggests that a mechanism similar to that described above (48) also operates in 5L cells.
Whereas the induction of ARE-regulated proteins by TCDD is generally expected to result in enhanced protection against oxidative stress, we also observed alterations of protein abundances that would be expected to enhance the susceptibility to oxidants. Thus, a protein spot identified as the aldehyde reductase, aldose reductase-related protein 2 (number 43), was down-regulated by
50% following exposure of the cells to TCDD. The protein, also termed Akr1b8 protein or FR-1, rapidly reduces various hydrophobic aldehydes, such as 4-hydroxynonenal, as well as aromatic aldehydes (83). The expression of this enzyme, which shares some substrates with ALDH3A1 and GSTs of the Alpha family, has been shown to be positively regulated by fibroblast growth factor 1 and related peptides (84). Whether the down-regulation of aldose reductase-related protein 2 by TCDD actually reflects an interference of the dioxin with growth factor signaling remains to be determined.
Cytoskeletal Proteins—
The present study corroborates our recent observation that TCDD strongly impinges on the organization of the cytoskeleton in 5L cells (37), and it markedly expands the knowledge of the proteins involved. A particularly interesting member of the group of actin-interacting proteins newly discovered to be affected by TCDD treatment was drebrin, which was down-regulated by
50%. Drebrin has been recently identified as a connexin-43 binding partner that links gap junctions to the submembrane actin cytoskeleton and that is required for maintaining connexin-43-containing gap junctions in their functional state at the plasma membrane (85). Depletion of drebrin protein with small interfering RNA was shown to result in impaired cell-cell coupling, internalization of gap junctions, and targeting of connexin-43 for degradation (85). Thus, it appears that the down-regulation of drebrin may be conducive to the inhibition of gap junctional intercellular communication that occurs in TCDD-treated rat hepatocytes in primary culture and that has been implicated in tumor promotion by TCDD (86).
Diry et al. (87) recently reported that in human MCF7 breast cancer cells TCDD exposure resulted in morphological changes associated with a prominent cytoskeletal remodeling, increases in cell plasticity and motility, and a down-regulation of E-cadherin. Human HepG2 hepatoma cells responded similarly. These changes, which were dependent on an AhR-dependent activation of JNK and observed after a treatment period of 48 h, were tentatively linked to a potential action of dioxin on tumor progression. At present it is unclear whether the alterations of the numerous cytoskeletal proteins detected in our proteomics studies after 8 h of TCDD exposure are mechanistically related to the much later effects observed by Diry et al. (87).
Proteins Involved in the Maintenance of Mitochondrial Homeostasis—
Several observations of the present and our previous studies show that TCDD exposure of 5L cells results in functionally important alterations of the mitochondrial proteome, and they suggest that some of these alterations may occur in response to the induction of mitochondrial oxidative stress. One clue is the observation of an up-regulation of the β chain of the mitochondrial ATP synthase as up-regulation of ATP synthase by toxicants has been attributed to an adaptation of the cells to mitochondrial damage induced by ROS (52, 53). For liver mitochondria of rats treated with lipopolysaccharide, up-regulation of the ATP synthase
chain was shown to be accompanied by an up-regulation of mitochondrial manganese-dependent superoxide dismutase (SOD2), and the expression of both enzymes was correlated with ROS generation suggesting that it represented a coordinated response to compensate for low levels of ATP and to increase mitochondrial antioxidant capacity (53). In fact, for rat and mouse liver, the induction of mitochondrial oxidative stress by TCDD is well established (88–91). Recent studies reported that the increased mitochondrial reactive oxygen production observed during treatment of rats with TCDD was accompanied by decreased ATP levels in the liver (90, 91). Our observations that TCDD induces oxidative stress and up-regulates at least one of the ATP synthase subunits as well as mitochondrial superoxide dismutase (37) in 5L cells suggest that a similar ROS-dependent adaptive mechanism is operative in mitochondria of TCDD-exposed 5L cells.
The notion that TCDD causes an imbalance in ATP demand and supply in 5L cells is indirectly supported by results from our previous ICPL study (37), which had revealed an interference by TCDD with the abundance or PTM of ADP/ATP translocase 1 (ANT1), a component of the inner mitochondrial membrane, and the glycolytic enzyme hexokinase I. ANT1 and the mitochondrial hexokinases I and II (HKI and HKII) interact with the VDACs of the outer mitochondrial membrane (92–94) that constitute the major pores for metabolic exchange of ATP, ADP, and inorganic phosphate and that have been suggested to act as dynamic regulators of global mitochondrial function by ultimately controlling the activity of the respiratory chain (95). The observed TCDD-induced alterations of ANT1 and HKI may thus indicate an interference of TCDD with the function of VDAC, although the precise nature of the alterations remains to be defined.
A particularly important outcome of the present study was the discovery that TCDD markedly affects the VDAC isoform VDAC2 in 5L cells. VDAC2, in addition to its function as a pore, mediates the import of the proapoptotic protein BAK into the outer mitochondrial membrane (54), inhibits the activity of BAK by preventing its oligomerization (55), and, thus, provides a connection between mitochondrial physiology and the core apoptotic pathway. Our data show that TCDD exposure results in both an altered PTM of VDAC2 and a marked increase in its abundance. The nature of the difference of the two forms of VDAC2, which may be related to the two VDAC2 species detected by 2-DE analysis of various rat tissues by Yamamoto et al. (96), is still unclear. VDAC2 is known to be subject to tyrosine phosphorylation (97) and lysine acetylation (98), but the very small pI difference of the two observed spots of <0.1 pH unit suggests that the two forms do not differ by one of these PTMs but, presumably, some other modification.
Up-regulation of VDAC2 protein was dependent on the AhR in 5L cells as indicated by the absence of any effect of TCDD on the amount of VDAC2 in the AhR-deficient 5L subline BP8. In contrast to the proteins of the AhR and ARE gene batteries that are regulated at the level of transcription, quantitative PCR analyses showed that up-regulation of VDAC2 was not mediated by an increase in the steady-state level of VDAC2 mRNA. We therefore suppose that the increase in the amount of VDAC2 protein during exposure to TCDD occurs at the translational level or via a prolongation of the half-life of the protein. We are not aware of any other report on an increase in VDAC2 protein in cells in response to exposure to chemicals such as drugs or xenobiotics. VDAC2 was one of a large number of genes of mouse cerebral cortical neurons with increased expression in response to hypoxia (99), but as this increase occurred at the transcriptional level it must have been mediated by a different mechanism. Because retroviral overexpression of VDAC2 has been shown to result in an inhibition of the mitochondrial apoptotic pathway via an inhibition of BAK activation (55), the up-regulation of VDAC2 protein observed in our study may constitute a general, hitherto unrecognized component of a concerted endeavor of cells to prevent cell death due to an imbalance between cellular energy demand and mitochondrial ATP supply. As the function of VDAC2 as a pore in the outer mitochondrial membrane is shared with that of the more abundant (96) VDAC1, the increased capability of the cells to keep BAK in check may be more momentous than the increase in permeability for metabolic exchange provided by the up-regulation of VDAC2 protein.
Taken together, our findings indicate that within 8 h of treatment TCDD induces concerted changes in the mitochondrial proteome that probably occur to counteract oxidative stress, to prevent ROS-induced energy depletion, and to maintain cell viability by sustaining ATP supply and inhibiting the triggering of the intrinsic, mitochondria-mediated apoptotic pathway. The observation that the cells nevertheless undergo apoptosis that becomes detectable in TCCD-exposed 5L cells after about 30 h of treatment2 shows, however, that these alterations cannot prevent TCDD-induced cell death in the long run.
Comparison of the Results with Those from the ICPL Study (37)—
In our previous ICPL study on the effects of TCDD on the proteome of 5L cells that used cell culture and treatment conditions identical to those of the present investigation, 89 different proteins were identified as up- or down-regulated by TCDD (37). This figure is similar to the number of 78 protein species, which corresponded to 73 different proteins, detected as "regulated" in the present study. Interestingly an inspection of the two protein lists reveals that the overlap between the two studies, i.e. the number of proteins detected by both methods, amounted to only nine (numbers 6, 9, 35, 42, 50, and 57–60). Of the 80 proteins detected only by the ICPL method, 27 proteins, largely histones and ribosomal proteins, had pI values above 9.5 and would thus be expected to cause problems in 2-DE analyses. Five had a molecular mass of >200 kDa and seven had a molecular mass of <15 kDa and were, therefore, also outside the range readily amenable to standard 2-DE analysis. Of the remaining 41 proteins, several may have been too hydrophobic for an analysis by 2-DE, but for the majority of them there is no straightforward explanation why they were not detected in the present study. Conversely it is not readily apparent why 64 altered proteins were only detected using the 2-DE but not the ICPL method. However, this limited concordance is actually not too surprising as several studies have shown that the choice of the methodology used for global proteome analysis strongly affects the spectrum of proteins detected (compare e.g. Refs. 38, 39, and 100). The results of our investigations on the impact of TCDD on the proteome of 5L cells demonstrate that the concurrent use of the ICPL and the 2-DE approach for quantitative proteome analysis can in fact yield highly complementary information on a biological system under study.
| FOOTNOTES |
|---|
Published, MCP Papers in Press, November 12, 2007, DOI 10.1074/mcp.M700258-MCP200
1 The abbreviations used are: TCDD, 2,3,7,8-tetrachlorodibenzo-p-dioxin; AFAR1, aflatoxin B1-aldehyde reductase; AhR, aryl hydrocarbon receptor; AHRE, aryl hydrocarbon response element; ALDH3A1, aldehyde dehydrogenase 3A1; ANT1, ADP/ATP translocase 1; ARE, antioxidant response element; ARNT, AhR nuclear translocator; CCT, chaperonin containing t-complex polypeptide 1; CYP, cytochrome P450; 2-DE, two-dimensional gel electrophoresis; DHR123, dihydrorhodamine 123; DTE, dithioerythritol; EROD, 7-ethoxyresorufin O-deethylase; HK, hexokinase; ICPL, isotope-coded protein label; JNK, c-Jun N-terminal kinase; NQO1, NAD(P)H:quinone oxidoreductase 1; Nrf2, nuclear factor E2 p45-related factor 2; PTM, posttranslational modification; R123, rhodamine 123; ROS, reactive oxygen species; TCP1, t-complex polypeptide 1; VDAC, voltage-dependent anion channel-selective protein. ![]()
2 U. Andrae, unpublished observation. ![]()
3 F. Kiefer, personal communication. ![]()
* The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. ![]()
S The on-line version of this article (available at http://www.mcponline.org) contains supplemental material. ![]()
¶ To whom correspondence should be addressed: GSF-Forschungszentrum für Umwelt und Gesundheit, Institut für Toxikologie, Ingolstädter Landstr. 1, D-85764 Neuherberg, Germany. Tel.: 49-89-3187-2221; Fax: 49-89-3187-3449; E-mail: andrae{at}gsf.de
| REFERENCES |
|---|
|
|
|---|
interacts with and inhibits the permeability transition pore in cardiac mitochondria.
Circ. Res. 92, 873
–880This article has been cited by other articles:
![]() |
V. Gillardin, F. Silvestre, M. Dieu, E. Delaive, M. Raes, J.-P. Thome, and P. Kestemont Protein Expression Profiling in the African Clawed Frog Xenopus laevis Tadpoles Exposed to the Polychlorinated Biphenyl Mixture Aroclor 1254 Mol. Cell. Proteomics, April 1, 2009; 8(4): 596 - 611. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Carpi, M. Korkalainen, L. Airoldi, R. Fanelli, H. Hakansson, V. Muhonen, J. Tuukkanen, M. Viluksela, and R. Pastorelli Dioxin-Sensitive Proteins in Differentiating Osteoblasts: Effects on Bone Formation In Vitro Toxicol. Sci., April 1, 2009; 108(2): 330 - 343. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Journal of Biological Chemistry |
| Journal of Lipid Research | ASBMB Today |