Quantitative Phosphoproteomic Analysis of Early Alterations in

The Role of AEG-1/MTDH/LYRIC in the Pathogenesis of Central Nervous System Disease. Evan K. Noch , Kamel Khalili. 2013,159-192 ...
2 downloads 0 Views 2MB Size
Article pubs.acs.org/jpr

Quantitative Phosphoproteomic Analysis of Early Alterations in Protein Phosphorylation by 2,3,7,8-Tetrachlorodibenzo‑p‑dioxin Melanie Schulz,†,‡ Stefanie Brandner,† Carola Eberhagen,† Friederike Eckardt-Schupp,§ Martin R. Larsen,‡ and Ulrich Andrae*,† †

Institute of Molecular Toxicology and Pharmacology, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany ‡ Department of Biochemistry and Molecular Biology, University of Southern Denmark, Odense, Denmark § Institute of Radiation Biology, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany S Supporting Information *

ABSTRACT: A comprehensive quantitative analysis of changes in protein phosphorylation preceding or accompanying transcriptional activation by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in 5L rat hepatoma cells was performed using the SILAC approach. Following exposure of the cells to DMSO or 1 nM TCDD for 0.5 to 2 h, 5648 phosphorylated peptides corresponding to 2156 phosphoproteins were identified. Eight peptides exhibited a statistically significantly altered phosphorylation because of TCDD exposure and 22 showed a regulation factor of ≥1.5 in one of the experiments per time point. The vast majority of the TCCD-induced phosphorylation changes had not been reported before. The transcription factor ARNT, the obligate partner for gene activation by the TCDD-bound Ah receptor, exhibited an up-regulation of its Ser77 phosphorylation, a modification known to control the differential binding of ARNT homodimers and heterodimers to different enhancers suggesting that this phosphorylation represents a novel mechanism contributing to the alteration of gene expression by TCDD. Other proteins with altered phosphorylation included, among others, various transcriptional coregulators previously unknown to participate in TCDD-induced gene activation, regulators of small GTPases of the Ras superfamily, UBX domain-containing proteins and the oncogenic protein LYRIC. The results open up new directions for research on the molecular mechanisms of dioxin action and toxicity. KEYWORDS: TCDD, dioxin, protein phosphorylation, phosphoproteomics, SILAC, SIMAC, ARNT, transcriptional regulation, GTPases, 5L cells



regulatory regions of aryl hydrocarbon-responsive genes.6−8 Binding to the DREs initiates chromatin remodeling and recruitment of coregulators9 which enables the transcription of a set of target genes comprising the so-called “AhR gene battery”. The induced gene products include several xenobioticmetabolizing enzymes such as cytochrome P4501A1 (CYP1A1). During the past few years, evidence has been accumulated that TCDD can also cause biochemical alterations that are easily detectable within a few minutes after addition of the compound to cells in vitro. These effects cannot be the consequence of AhR/ARNT-mediated gene activation as this is dependent on a stimulation of de novo transcription and translation which require a certain minimal time in the range of a few hours to become apparent. For these responses, which in the majority of the cases have been claimed to be dependent on a functional AhR, but not ARNT10 and which appear to be

INTRODUCTION 2,3,7,8-Tetrachlorodibeno-p-dioxin (TCDD) is a ubiquitous and persistent environmental contaminant with the potential to cause a broad spectrum of adverse health effects including a variety of developmental defects, cardiovascular disease, diabetes, porphyria, hormonal disturbances and cancer in humans.1,2 In spite of enormous efforts to uncover the mechanisms of action underlying the various deleterious effects of TCDD, the precise molecular events mediating toxicity are still largely unknown. TCDD is a ligand of the so-called aryl hydrocarbon receptor (AhR), a primarily cytosolic receptor which plays important roles in normal cell physiology3 and which mediates the toxic effects of TCDD.4,5 Binding of TCDD or a wide variety of other AhR agonists to the 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 specific DNA sequences known as xenobiotic response elements or dioxin response elements (DREs), enhancers located in © 2013 American Chemical Society

Received: October 8, 2012 Published: January 8, 2013 866

dx.doi.org/10.1021/pr3009429 | J. Proteome Res. 2013, 12, 866−882

Journal of Proteome Research

Article

kinase (MAPK) in an AhR-dependent, but, importantly, transcription-independent manner, suggesting a novel crosstalk between AhR and MAPK signaling. p38 activation was observed as early as 1 h after addition of TCDD and resulted in subsequent transcriptional induction of the proto-oncogene cjun. The signaling pathway resulting in p38 phosphorylation and the kinase involved remained unclear. In the present phosphoproteomic analysis of the effects of TCDD in 5L cells, short exposure periods of 0.5, 1, and 2 h were used, that is, treatments that were unlikely to result in significant alterations in the abundance of proteins with TCDD-regulated expression, which requires some time because of its dependence on de novo transcription and translation. The study was performed with the expectation to uncover proteins not previously known to exhibit rapid alterations in phosphorylation as a consequence of TCDD exposure and which might give clues to pathways involved in nongenomic signaling preceding alterations in gene expression. Moreover, it was anticipated to identify proteins with altered phosphorylation as a consequence of their involvement in the regulation of transcription induced by the DRE-bound liganded AhR/ ARNT heterodimer. Finally, the study was expected to provide novel information on the phosphoproteome of rat cells with respect to the identification of novel phosphorylation sites of known phosphoproteins and the discovery of proteins previously unknown to be phosphorylated.

largely mediated by rapid phosphorylation events, the term “protein phosphorylation pathway”11 or, more recently, “nongenomic” pathway was coined to discriminate them from the classical “genomic” pathway involving Ahr/ARNT-mediated induction of gene expression. Whereas the details of the phosphorylation-mediated pathways may vary between different types of cells, a salient feature in all cells examined is that they appear to be initiated by a rapid increase in intracellular Ca2+ concentration which has been consistently demonstrated in organ cultures of TCDD-exposed animals,12 primary cell cultures13−17 and mammalian cell lines.18−20 This increase in intracellular calcium, which takes place within a few minutes, is the earliest effect detectable in TCDD-exposed cells and has been suggested to constitute the trigger of the nongenomic action of TCDD.10 For certain cell lines, the elevation of the intracellular Ca2+ concentration by TCDD has been reported to result in the activation of the ubiquitous Ca2+-dependent cytosolic phospholipase A2 (cPLA2) resulting in the cleavage of arachidonic acid (AA) from glycerophospholipids.21−23 TCDDinduced cPLA2 activation was also detected via an increase in the activating Ser505 phosphorylation of the enzyme.21,24 The AA formed has been implicated in the rapid elevation of Src kinase activity that occurs in several cell types and which takes place in MCF10A cells within 30 min.21,25 Quite recently, a novel branch of the nongenomic pathway linking the TCDD-induced rapid elevation of intracellular Ca2+ concentration to the up-regulation of AhR target genes, that is, to the genomic response, via phosphorylation events has been identified.20,26,27 For MCF-7 cells it was shown that the increase in intracellular Ca2+ concentration results in a transient activation of the ubiquitous calmodulin (CaM)-dependent protein kinase Iα (CaMKI α) which turned out to be a necessary prerequisite for the TCDD-triggered nuclear import of the ligand-activated AhR and promoter activation of AhR target genes. Confirmatory data were presented for primary human macrophages.20 Evidence was obtained that exposure to TCDD was accompanied by an increased Thr177 phosphorylation of CaMKI α,20 which has been associated with full CaMKI α activation.28 The protein kinase mediating CaMKI α activation remained unclear.20 Whereas the genomic response to TCDD has been analyzed in numerous large-scale studies at the level of the transcriptome and the proteome both in vitro and in vivo, a systematic study of the proteins involved in the early nongenomic action of TCDD has not been performed previously. Therefore, the present study was conducted with the aim to identify TCDDinduced alterations in protein phosphorylation occurring early after the onset of TCDD exposure to achieve a better knowledge of the alterations in signal transduction potentially involved in the nongenomic effects of TCDD. The rat hepatoma cell line 5L was used as a model system. 5L cells are epithelial-like cells which express both the AhR and the ARNT proteins and which are highly responsive to TCDDinduced gene activation. They have been previously used in a variety of studies dealing with the role of AhR in cell cycle control and TCDD toxicity,29−37 the regulation of the expression of xenobiotic metabolizing enzymes,38−40 and the effects of AhR and TCDD on microRNA levels.41 In addition, they have been successfully employed in studies aiming at the detection of novel target genes of the AhR.42,43 With respect to nongenomic actions of TCDD in this cell line, Weiss et al.32 observed that TCDD exposure of the cells is associated with the rapid phosphorylation of the p38 Mitogen-activated protein



EXPERIMENTAL PROCEDURES

Cells

The 5L cell line is a descendant of the cell line H4IIEC3 established by Pitot et al.44 from the Reuber H35 rat hepatoma.45 Cells were grown as monolayers in 10-cm culture dishes containing RPMI 1640 medium supplemented with 10% FCS, 100 U/mL penicillin, and 100 μg/mL streptomycin at 37 °C in a humidified atmosphere of 95% air and 5% CO2. SILAC Labeling

For stable isotope labeling of the 5L cells using the SILAC (stable isotope labeling by amino acids in cell culture) approach, 13C6-lysine (Lys6) (Invitrogen, Carlsbad, U.S.A.), 15 N4-arginine (Arg4) (Cambridge Isotope Laboratories, Andover, U.S.A.) and 13C615N4-arginine (Arg10) (Invitrogen) were used according to the labeling scheme of van Hoof et al.46 in order to circumvent potential peptide quantitation problems associated with the conversion of isotopically labeled arginine into labeled proline. The amino acids were dissolved in RPMI medium without arginine and lysine supplemented with 10% dialyzed fetal bovine serum, 100 μg/mL streptomycin, and 100 U/ml penicillin. For the “light” condition, the medium was supplemented with Lys0 and Arg4, and for the “heavy” condition with Lys6 and Arg10. The concentrations of arginine and lysine were 75 mg/L. The cells were grown for about 7 cell cycles (7 days) with intermediate splitting after 3 days and replating in labeling media at a density of 7.5 × 105 cells/10-cm dish. Treatment of the Cells with TCDD

During the last 16 h before the initiation of the TCDD exposure, the cells were starved in fresh labeling medium without serum. The “heavy” cells were treated with 1 nM TCDD (purity ≥99%, Ö kometric, Bayreuth, Germany) dissolved in DMSO and the “light” cells with DMSO alone as a control for 0.5, 1, and 2 h. The DMSO concentration in the 867

dx.doi.org/10.1021/pr3009429 | J. Proteome Res. 2013, 12, 866−882

Journal of Proteome Research

Article

Hydrophilic Interaction Chromatography (HILIC)

medium was 0.1%. Safety considerations: Skin contact with TCDD-containing solutions and contamination of laboratory benches and laminar flow hoods with TCDD must be avoided.

The monophosphorylated and nonphosphorylated peptides from the SIMAC enrichment were fractionated by HILIC. The dried fraction was resuspended in B-solvent (80% ACN, 0.1% TFA) and loaded onto a 2 mm TSKGel Amide-80 HILIC column (Tosoh Bioscience, Stuttgart, Germany) connected to an Ä KTA Basic FPLC system (GE Healthcare, Hillerod, Denmark). Loading was performed in B-solvent at a flow rate of 200 μL/min for 12 min. The peptides were separated by an inverse gradient (100%−70% B-solvent within 28 min and 70%−0% B-solvent within 5 min) with an increasing amount of A-solvent (0.1% TFA) at a flow rate of 200 μL/min. Sixteen fractions were collected and dried down in a vacuum centrifuge.

Cell Lysis and Protein Precipitation

Following exposure of the cells the medium was removed and cells were washed twice with 10 mL ice-cold PBS. Then 1.5 mL of lysis buffer containing 1 μL/mL benzonase, 1 mM MgCl2, 25 mM Tris HCl, pH 7.4, 120 mM NaCl, 1% Triton X-100, one protease inhibitor tablet (Complete, Roche, Mannheim, Germany) per 50 mL buffer, 0.5 mL of phosphatase inhibitor cocktails 1 and 2 (Sigma-Aldrich) per 50 mL, and 1 mM sodium pervanadate were added to each 10-cm dish. The sodium pervanadate solution was prepared by treating a solution of 100 mM sodium orthovanadate, pH 10, with an equal volume of a 0.36% hydrogen peroxide solution immediately before use. The lysates were transferred to Eppendorf tubes, shaken for 30 min at 4 °C and centrifuged at 15000g for 10 min at 4 °C. Following determination of the protein concentrations with the Bradford assay, the supernatants from the “light” and the “heavy” cells were mixed 1:1. The proteins (800 μg) were precipitated with four volumes of ice-cold ethanol at −20 °C. After 15 min, four volumes of icecold acetone were added and the mixture was incubated at −20 °C overnight. The precipitated proteins were centrifuged at 15 000g for 20 min at 4 °C. The pellet was washed 3 times with ethanol/acetone/water 2:2:1.

Phosphopeptide Enrichment by Titanium Dioxide (TiO2) Chromatography

To enrich the mainly monophosphorylated peptides from the individual HILIC fractions, the TiO2 method was used essentially as described by Larsen et al.49 with slight modifications.50 TiO2 beads (GL Sciences, Tokyo, Japan) were washed with water, 100% ACN and loading buffer (1 M glycolic acid in 80% ACN, 5% TFA). Each HILIC fraction was dissolved in 5 μL 1% SDS, 200 μL loading buffer, and the appropriate amount of TiO2 was added. For 200 μg of peptides 0.6 mg TiO2 particles were used. The peptides and the TiO2 resin were incubated for 15 min under vigorous shaking. The beads were spun down by centrifugation at 2000g for 30 s. The supernatant was collected in a fresh tube and the TiO2 beads were washed with loading buffer (80% ACN/5% TFA), then 30% ACN/0.1% TFA and, finally, distilled water. The phosphopeptides were eluted from the beads with 50 μL 1% NH4OH using an incubation time of 10 min under vigorous shaking. The supernatant was acidified with 5 μL 100% FA. The phosphopeptides were desalted using Poros R3 P10 micro columns and analyzed by LC-MS/MS. The collected supernatant from the first TiO2 incubation with loading buffer was again treated with TiO2 as described above.

Proteolysis

The precipitated proteins were dissolved in 6 M urea, 2 M thiourea and 50 mM NH4HCO3. The proteins were reduced for 2 h at 30 °C with 20 mM DTT and alkylated with 40 mM iodoacetamide for 45 min in the dark at room temperature. Then they were digested with endoproteinase Lys-C (Wako, Richmond, U.S.A.) (1/100 w/w) for 3 h. For digestion with trypsin, the solutions were diluted in 50 mM NH4HCO3 to a urea concentration of less than 1 M and incubated with sequencing grade modified trypsin (Promega, Madison, U.S.A.) (1/50 w/w) at 30 °C for 16 h. The digestion was quenched by adding trifluoroacetic acid (TFA) to reach a pH of ≤3.

StageTip Peptide Purification

Peptides were purified using modified StageTips. StageTips were prepared by stamping out a small plug of a 3 M Empore C8 disk (3 M Bioanalytical Technologies, St. Paul, U.S.A.) and placing it at the end of a P10 tip. Poros R3 material in ACN was packed in the pipet tip onto the top of the C8 disk. The microcolumn was washed with 100% ACN and equilibrated with 0.1% TFA. Phosphopeptide samples were loaded onto the columns using a 1 mL disposable syringe. The columns were washed once with 0.1% TFA in water and the peptides were eluted with 50% ACN and 0.1% TFA. Samples were dried down in a vacuum centrifuge and redissolved in 0.5 μL 100% FA and 4.5 μL 0.1% FA prior to analysis by LC-MS/MS.

Phosphopeptide Enrichment by Sequential Elution from IMAC (SIMAC)

The phosphorylated peptides were enriched using SIMAC essentially as described.47 Briefly, PHOS-Select Iron Affinity Gel beads (Sigma-Aldrich) were washed twice with 0.1% TFA and once with loading buffer (50% acetonitrile (ACN), 0.1% TFA) to remove the glycerol. The beads were resuspended in loading buffer and the sample was added. For 200 μg of peptides, 50 μL gel beads were used. After incubation by endover-end rotation for 1 h, the peptide-bead mixture was transferred to a constricted P200 gel loading tip. The flowthrough and the first wash with loading buffer (100 μL) were collected and combined. The monophosphorylated peptides were eluted with 20% ACN, 1% TFA (100 μL) and combined with the flow-through. This solution was dried down in a vacuum centrifuge. The multiphosphorylated peptides were eluted off the IMAC material with 1% NH4OH (70 μL) and this solution was subsequently acidified with 100% formic acid (FA) to a final FA concentration of 1%. The multiphosphorylated peptides were desalted using Poros Oligo R3 reversed phase material (PerSeptive Biosystems, Framingham, U.S.A.) packed in P10 micropipet tips as previously described48 and analyzed by LC-MS/MS.

Phosphopeptide Identification and Quantitation

LC-MS/MS Analysis. The enriched and purified phosphopeptides were analyzed on a LTQ Orbitrap XL Mass Spectrometer (Thermo Fisher Scientific, Bremen, Germany). The samples were applied onto an EASY nano-HPLC system (Thermo Fisher Scientific, Odense, Denmark) comprising a self-packed 16-cm analytical column (100 μm I.D., 360 μm O.D., ReproSil-Pur C18 AQ 3 μm (Dr. Maisch, AmmerbuchEntringen, Germany). The peptides were loaded with a flow rate of 550 nL/min in 0.1% FA (100% phase A) as solvent and eluted at a flow rate of 250 nL/min with a gradient from 100% phase A to 34% phase B (95% ACN, 0.1% FA) for 100 min. 868

dx.doi.org/10.1021/pr3009429 | J. Proteome Res. 2013, 12, 866−882

Journal of Proteome Research

Article

peptide with the highest ion score was kept. If one peptide was detected with different IPI accessions, only the peptide belonging to the IPI identifier with the highest number of different peptides observed was retained.54 Then, a normalization of the calculated log2 peptide ratios was performed. A global normalization method involving the adjustment of the log median of the peptide ratios to zero was used. Statistical Analysis. Three different, independent biological replicates were analyzed for the data sets obtained for the 30-min and 1-h exposures to TCDD and two biological replicates were analyzed for the 2-h data set. Only peptides found at least twice for a given time point were used for the following statistical test. The normalized ratios were fitted to a linear model and an empirical Bayes moderated t test was conducted using the R/Bioconductor package “limma”. The resulting p-values were used to calculate the q-values36 with the R-package “q-value”. Peptides with a q-value ≤0.05 were accepted as statistically significantly regulated. Volcano plots were generated with the limma package. Annotation. Following statistical analysis, the IPI identifier of every unique phosphorylated peptide was mapped to a Uniprot, GO, KEGG, REFSEQ, and ENSEMBL identifier. Mapping was done using the software R in combination with the packages BioIDMapper and biomaRt. Phosphoset Analysis. The identification of known and hitherto unknown phosphorylation sites was performed with an in-house written Perl script. For this analysis, the Swiss-Prot accession numbers were used and searched against the UniProt database.

The nano-HPLC system was connected online to the mass spectrometer with a nano ES ion source and PicoTip emitter needles (New Objective, Berlin, Germany). The mass spectrometer was operated in positive ion mode and a datadependent acquisition with an automatic switch between MS and MS/MS was employed. Each MS scan was followed by 5 MS/MS scans. The full scan was acquired in the Orbitrap with an automatic gain control (AGC) target value of 1 × 106 ions and a maximum fill time of 500 ms. The top 5 most intense ions were selected for fragmentation in the linear trap quadrupole (LTQ). Fragmentation was induced by CID. For the LTQ, the AGC target value was set to 3 × 104 ions and a maximum fill time of 300 ms. For an improved fragmentation of phosphopeptides, the multistage algorithm51 was enabled for each MS/MS spectrum, which showed a loss of phosphoric acid from the parent ion (neutral loss of 97.97, 48.99, or 32.66 Thomson). Ions selected for MS/MS were dynamically excluded for a duration of 45 s. Database Searching. The acquired multistage activation spectra were converted into Mascot generic format (mgf) using the software DTASuperCharge, version 1.37 (http://msquant. sourceforge.net/). Deisotoping was performed using the default settings from DTASuperCharge. The processed spectra were searched against the rat IPI protein database version 3.60 which contains 39863 protein sequences. The search was performed using an in-house Mascot server (version 2.2) with carbamidomethyl (C) as a fixed modification and acetylation (protein N-terminal), oxidation (M), phosphorylation (ST and Y) and the SILAC labeling (K-6, R-4, R-10) as variable modifications. Trypsin was chosen as the enzyme and up to two miss-cleavages were allowed. For the parent ion a mass accuracy of 10 ppm was used and for the fragment ions a mass accuracy of 0.8 Da. A concatenated decoy database search was performed derived from the IPI rat database. Proteins with a false discovery rate (FDR) of ≤1% and with peptides having an ion score of ≥20 (MudPit scoring) were accepted. This ion score cutoff was slightly lower than that of ≥22 recently employed by Engholm-Keller et al.52 and selected because an evaluation of a representative part of low-scoring MS/MS spectra of the peptides identified in the present study had indicated that an ion score cutoff of ≥22 would have been too conservative as an ion score of ≥20 yielded highly reliable results (not shown). Quantitation and Phosphosite Localization. Relative quantitation of phosphopeptides and nonphosphorylated peptides was performed using the software MSQuant (version 1.0a70). To enable the processing of the data generated with the labeling scheme of van Hoof et al.,46 the source code of MSQuant was changed. MSQuant quantifies the peptides based on the first isotopic peak of the isotopically labeled peptides by averaging the ion intensities over multiple MS scans. As a measure of the reliability of the phosphorylation site assignment, the posttranslational modification (PTM) score for every phosphosite was calculated using MSQuant.53 Thereafter, every MS and MS/MS spectrum was manually validated. Peptides with a low signal-to-noise ratio, ion score 1.5. bThe asterisk (*) in the sequence indicates the position of the phosphorylation. cThe respective phosphorylation site of the protein had previously been observed for species other than the rat (known by similarity) or not been observed for any species (novel). dThe normalized SILAC ratio TCDD/DMSO is the normalized SILAC ratio “heavy/light” (H/L). eThe indicated two alternative positions of the phosphorylation site correspond to two splice isoforms that had been observed for the rat at the mRNA level.

871

dx.doi.org/10.1021/pr3009429 | J. Proteome Res. 2013, 12, 866−882

Journal of Proteome Research

Article

Table 2. Phosphopeptides Identified As Regulated (q-Value 6500 different phosphopeptides revealed that only a tiny fraction, less than 0.5%, exhibited evidence for a regulation of their phosphorylation status by TCDD. This finding supports the notion that early alterations in protein phosphorylation by dioxins are likely 878

dx.doi.org/10.1021/pr3009429 | J. Proteome Res. 2013, 12, 866−882

Journal of Proteome Research

Article

(4) Fernandez-Salguero, P. M.; Hilbert, D. M.; Rudikoff, S.; Ward, J. M.; Gonzalez, F. J. Aryl-hydrocarbon receptor-deficient mice are resistant to 2,3,7,8-tetrachlorodibenzo-p-dioxin-induced toxicity. Toxicol. Appl. Pharmacol. 1996, 140, 173−9. (5) Bunger, M. K.; Glover, E.; Moran, S. M.; Walisser, J. A.; Lahvis, G. P.; Hsu, E. L.; Bradfield, C. A. Abnormal liver development and resistance to 2,3,7,8-tetrachlorodibenzo-p-dioxin toxicity in mice carrying a mutation in the DNA-binding domain of the aryl hydrocarbon receptor. Toxicol. Sci. 2008, 106, 83−92. (6) Safe, S. Molecular biology of the Ah receptor and its role in carcinogenesis. Toxicol. Lett. 2001, 120, 1−7. (7) Sogawa, K.; Numayama-Tsuruta, K.; Takahashi, T.; Matsushita, N.; Miura, C.; Nikawa, J.-i.; Gotoh, O.; Kikuchi, Y.; Fujii-Kuriyama, Y. A novel induction mechanism of the rat CYP1A2 gene mediated by Ah receptor-Arnt heterodimer. Biochem. Biophys. Res. Commun. 2004, 318, 746−55. (8) Boutros, P. C.; Moffat, I. D.; Franc, M. A.; Tijet, N.; Tuomisto, J.; Pohjanvirta, R.; Okey, A. B. Dioxin-responsive AHRE-II gene battery: identification by phylogenetic footprinting. Biochem. Biophys. Res. Commun. 2004, 321, 707−15. (9) Beischlag, T. V.; Luis Morales, J.; Hollingshead, B. D.; Perdew, G. H. The aryl hydrocarbon receptor complex and the control of gene expression. Crit. Rev. Eukaryot. Gene Expr. 2008, 18, 207−50. (10) Matsumura, F. The significance of the nongenomic pathway in mediating inflammatory signaling of the dioxin-activated Ah receptor to cause toxic effects. Biochem. Pharmacol. 2009, 77, 608−26. (11) Enan, E.; Matsumura, F. Evidence for a second pathway in the action mechanism of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Significance of Ah-receptor mediated activation of protein kinase under cell-free conditions. Biochem. Pharmacol. 1995, 49, 249−61. (12) Canga, L.; Levi, R.; Rifkind, A. B. Heart as a target organ in 2,3,7,8-tetrachlorodibenzo-p-dioxin toxicity: decreased beta-adrenergic responsiveness and evidence of increased intracellular calcium. Proc. Natl. Acad. Sci. U.S.A. 1988, 85, 905−9. (13) McConkey, D. J.; Hartzell, P.; Duddy, S. K.; Håkansson, H.; Orrenius, S. 2,3,7,8-Tetrachlorodibenzo-p-dioxin kills immature thymocytes by Ca2+-mediated endonuclease activation. Science 1988, 242, 256−9. (14) Canga, L.; Paroli, L.; Blanck, T. J.; Silver, R. B.; Rifkind, A. B. 2,3,7,8-tetrachlorodibenzo-p-dioxin increases cardiac myocyte intracellular calcium and progressively impairs ventricular contractile responses to isoproterenol and to calcium in chick embryo hearts. Mol. Pharmacol. 1993, 44, 1142−51. (15) Hanneman, W. H.; Legare, M. E.; Barhoumi, R.; Burghardt, R. C.; Safe, S.; Tiffany-Castiglioni, E. Stimulation of calcium uptake in cultured rat hippocampal neurons by 2,3,7,8-tetrachlorodibenzo-pdioxin. Toxicology 1996, 112, 19−28. (16) Tannheimer, S. L.; Barton, S. L.; Ethier, S. P.; Burchiel, S. W. Carcinogenic polycyclic aromatic hydrocarbons increase intracellular Ca2+ and cell proliferation in primary human mammary epithelial cells. Carcinogenesis 1997, 18, 1177−82. (17) Lin, C.-H.; Juan, S.-H.; Wang, C. Y.; Sun, Y.-Y.; Chou, C.-M.; Chang, S.-F.; Hu, S.-Y.; Lee, W.-S.; Lee, Y.-H. Neuronal activity enhances aryl hydrocarbon receptor-mediated gene expression and dioxin neurotoxicity in cortical neurons. J. Neurochem. 2008, 104, 1415−29. (18) Puga, A.; Nebert, D. W.; Carrier, F. Dioxin induces expression of c-fos and c-jun proto-oncogenes and a large increase in transcription factor AP-1. DNA Cell Biol. 1992, 11, 269−81. (19) Puga, A.; Hoffer, A.; Zhou, S.; Bohm, J. M.; Leikauf, G. D.; Shertzer, H. G. Sustained increase in intracellular free calcium and activation of cyclooxygenase-2 expression in mouse hepatoma cells treated with dioxin. Biochem. Pharmacol. 1997, 54, 1287−96. (20) Monteiro, P.; Gilot, D.; Le Ferrec, E.; Rauch, C.; LagadicGossmann, D.; Fardel, O. Dioxin-mediated up-regulation of aryl hydrocarbon receptor target genes is dependent on the calcium/ calmodulin/CaMKIalpha pathway. Mol. Pharmacol. 2008, 73, 769−77. (21) Dong, B.; Matsumura, F. Roles of cytosolic phospholipase A2 and Src kinase in the early action of 2,3,7,8-tetrachlorodibenzo-p-

restricted to few signaling pathways. The observation that transcriptional regulators comprised the largest functional group of proteins with altered phosphorylation corroborates the suitability of the experimental approach applied as it is fully in line with the well-known activity of TCDD as a modulator of transcription. By contrast, the discovery that regulators of small GTPases constituted the second largest group was completely unexpected and sheds light on hitherto unknown activities of TCDD in cells. The importance of these alterations, like that of the phosphorylation changes of the other proteins identified in the present study, for the toxicity of TCDD can now be scrutinized in future investigations. It should also be noted that the identification of thousands of phosphorylation sites clearly unaffected by TCDD during short exposures has provided an invaluable source of information on regulatory pathways not subject to phosphorylation changes. These data could significantly aid in the selection of further experimental approaches dealing with the clarification of the mechanisms of TCDD toxicity. And, finally, the present study has identified phosphorylation changes which exhibited q-values just above the set threshold value and were therefore disregarded according to the stringent significance criteria applied. It appears very likely, however, that this group of phosphopeptides comprises further phosphorylation sites which are actually subject to regulation by TCDD. A closer look at these peptides might yield further stimuli for future targeted investigations on phosphorylation-mediated early alterations in cellular signaling by dioxins.



ASSOCIATED CONTENT

* Supporting Information S

Tables S1−S9 and Figures S1−S5. This material is available free of charge via the Internet at http://pubs.acs.org. Accession Codes

The MS/MS mass spectra of all individual phosphopeptides identified have been deposited with the ProteomeCommons.org Tranche network and may be accessed using the following URL: https://proteomecommons.org/data set.jsp?i= 7eOKZmaZM61cU58heX9amEydDZy7FjmVv Qi1JWnbHZznP8FRtlceK66QYGbrs7XBhg0As2ZJbp9JUZ2fP CyNogaozncAAAAAAABi2A%3D%3D



AUTHOR INFORMATION

Corresponding Author

*Tel.: +49 89 3187 2221. Fax: +49 89 3187 3449. E-mail: [email protected]. Notes

The authors declare no competing financial interest.



REFERENCES

(1) International Agency for Research on Cancer (IARC). IARC Monographs on the Evaluation of Carcinogenic Risks to Humans: Polychlorinated Dibenzo-para-dioxins and Polychlorinated Dibenzofurans, Vol. 69; International Agency for Research on Cancer: Lyon, France, 1997. (2) White, S. S.; Birnbaum, L. S. An overview of the effects of dioxins and dioxin-like compounds on vertebrates, as documented in human and ecological epidemiology. J. Environ. Sci. Health, Part C: Environ. Carcinog. Ecotoxicol. Rev. 2009, 27, 197−211. (3) Barouki, R.; Coumoul, X.; Fernandez-Salguero, P. M. The aryl hydrocarbon receptor, more than a xenobiotic-interacting protein. FEBS Lett. 2007, 581, 3608−15. 879

dx.doi.org/10.1021/pr3009429 | J. Proteome Res. 2013, 12, 866−882

Journal of Proteome Research

Article

dioxin through a nongenomic pathway in MCF10A cells. Mol. Pharmacol. 2008, 74, 255−63. (22) Li, W.; Matsumura, F. Significance of the nongenomic, inflammatory pathway in mediating the toxic action of TCDD to induce rapid and long-term cellular responses in 3T3-L1 adipocytes. Biochemistry 2008, 47, 13997−4008. (23) Dong, B.; Nishimura, N.; Vogel, C. F.; Tohyama, C.; Matsumura, F. TCDD-induced cyclooxygenase-2 expression is mediated by the nongenomic pathway in mouse MMDD1 macula densa cells and kidneys. Biochem. Pharmacol. 2010, 79, 487−97. (24) Sciullo, E. M.; Dong, B.; Vogel, C. F. A.; Matsumura, F. Characterization of the pattern of the nongenomic signaling pathway through which TCDD-induces early inflammatory responses in U937 human macrophages. Chemosphere 2009, 74, 1531−7. (25) Mazina, O.; Park, S.; Sano, H.; Wong, P.; Matsumura, F. Studies on the mechanism of rapid activation of protein tyrosine phosphorylation activities, particularly c-Src kinase, by TCDD in MCF10A. J. Biochem. Mol. Toxicol. 2004, 18, 313−21. (26) Le Ferrec, E.; Lagadic-Gossmann, D.; Rauch, C.; Bardiau, C.; Maheo, K.; Massiere, F.; Le Vee, M.; Guillouzo, A.; Morel, F. Transcriptional induction of CYP1A1 by oltipraz in human Caco-2 cells is aryl hydrocarbon receptor- and calcium-dependent. J. Biol. Chem. 2002, 277, 24780−7. (27) N’Diaye, M.; Le Ferrec, E.; Lagadic-Gossmann, D.; Corre, S.; Gilot, D.; Lecureur, V.; Monteiro, P.; Rauch, C.; Galibert, M.-D.; Fardel, O. Aryl hydrocarbon receptor- and calcium-dependent induction of the chemokine CCL1 by the environmental contaminant benzo[a]pyrene. J. Biol. Chem. 2006, 281, 19906−15. (28) Haribabu, B.; Hook, S. S.; Selbert, M. A.; Goldstein, E. G.; Tomhave, E. D.; Edelman, A. M.; Snyderman, R.; Means, A. R. Human calcium-calmodulin dependent protein kinase I: cDNA cloning, domain structure and activation by phosphorylation at threonine-177 by calcium-calmodulin dependent protein kinase I kinase. EMBO J. 1995, 14, 3679−86. (29) Göttlicher, M.; Cikryt, P.; Wiebel, F. J. Inhibition of growth by 2,3,7,8-tetrachlorodibenzo-p-dioxin in 5L rat hepatoma cells is associated with the presence of Ah receptor. Carcinogenesis 1990, 11, 2205−10. (30) Göttlicher, M.; Wiebel, F. J. 2,3,7,8-Tetrachlorodibenzo-p-dioxin causes unbalanced growth in 5L rat hepatoma cells. Toxicol. Appl. Pharmacol. 1991, 111, 496−503. (31) Wiebel, F. J.; Klose, U.; Kiefer, F. Toxicity of 2,3,7,8tetrachlorodibenzo-p-dioxin in vitro: H4IIEC3-derived 5L hepatoma cells as a model system. Toxicol. Lett. 1991, 55, 161−9. (32) Weiss, C.; Faust, D.; Dürk, H.; Kolluri, S. K.; Pelzer, A.; Schneider, S.; Dietrich, C.; Oesch, F.; Göttlicher, M. TCDD induces cjun expression via a novel Ah (dioxin) receptor-mediated p38-MAPKdependent pathway. Oncogene 2005, 24, 4975−83. (33) Weiss, C.; Kolluri, S. K.; Kiefer, F.; Gö t tlicher, M. Complementation of Ah receptor deficiency in hepatoma cells: negative feedback regulation and cell cycle control by the Ah receptor. Exp. Cell Res. 1996, 226, 154−63. (34) Ge, N. L.; Elferink, C. J. A direct interaction between the aryl hydrocarbon receptor and retinoblastoma protein. Linking dioxin signaling to the cell cycle. J. Biol. Chem. 1998, 273, 22708−13. (35) Reiners, J.; J J; Clift, R.; Mathieu, P. Suppression of cell cycle progression by flavonoids: Dependence on the aryl hydrocarbon receptor. Carcinogenesis 1999, 20, 1561−6. (36) Elferink, C. J.; Ge, N. L.; Levine, A. Maximal aryl hydrocarbon receptor activity depends on an interaction with the retinoblastoma protein. Mol. Pharmacol. 2001, 59, 664−73. (37) Levine-Fridman, A.; Chen, L.; Elferink, C. J. Cytochrome P4501A1 promotes G1 phase cell cycle progression by controlling aryl hydrocarbon receptor activity. Mol. Pharmacol. 2004, 65, 461−9. (38) Münzel, P. A.; Schmohl, S.; Buckler, F.; Jaehrling, J.; Raschko, F. T.; Köhle, C.; Bock, K. W. Contribution of the Ah receptor to the phenolic antioxidant-mediated expression of human and rat UDPglucuronosyltransferase UGT1A6 in Caco-2 and rat hepatoma 5L cells. Biochem. Pharmacol. 2003, 66, 841−7.

(39) Sawaya, R. M.; Riddick, D. S. Cytochrome P450 2C11 5′flanking region and promoter: Regulation by aromatic hydrocarbons in vitro. Toxicology 2008, 248, 104−12. (40) Vondrácek, J.; Krcmár, P.; Procházková, J.; Trilecová, L.; Gavelová, M.; Skálová, L.; Szotáková, B.; Buncek, M.; Radilová, H.; Kozubík, A.; Machala, M. The role of aryl hydrocarbon receptor in regulation of enzymes involved in metabolic activation of polycyclic aromatic hydrocarbons in a model of rat liver progenitor cells. Chem.Biol. Interact. 2009, 180, 226−37. (41) Moffat, I. D.; Boutros, P. C.; Celius, T.; Lindén, J.; Pohjanvirta, R.; Okey, A. B. microRNAs in adult rodent liver are refractory to dioxin treatment. Toxicol. Sci. 2007, 99, 470−87. (42) Kolluri, S. K.; Balduf, C.; Hofmann, M.; Göttlicher, M. Novel target genes of the Ah (dioxin) receptor: transcriptional induction of N-myristoyltransferase 2. Cancer Res. 2001, 61, 8534−9. (43) Kolluri, S. K.; Weiss, C.; Koff, A.; Göttlicher, M. p27(Kip1) induction and inhibition of proliferation by the intracellular Ah receptor in developing thymus and hepatoma cells. Genes Dev. 1999, 13, 1742−53. (44) Pitot, H. C.; Peraino, C.; Morse, P. A., Jr.; Potter, V. R. Hepatomas in tissue culture compared with adapting liver in vivo. Natl. Cancer Inst. Monogr. 1964, 13, 229−45. (45) Reuber, M. D. A transplantable bile-secreting hepatocellular carcinoma in the rat. J. Natl. Cancer Inst. 1961, 26, 891−9. (46) Van Hoof, D.; Pinkse, M. W. H.; Oostwaard, D. W.-V.; Mummery, C. L.; Heck, A. J. R.; Krijgsveld, J. An experimental correction for arginine-to-proline conversion artifacts in SILAC-based quantitative proteomics. Nat. Methods 2007, 4, 677−8. (47) Thingholm, T. E.; Jensen, O. N.; Robinson, P. J.; Larsen, M. R. SIMAC (sequential elution from IMAC), a phosphoproteomics strategy for the rapid separation of monophosphorylated from multiply phosphorylated peptides. Mol. Cell. Proteomics 2008, 7, 661−71. (48) Gobom, J.; Nordhoff, E.; Mirgorodskaya, E.; Ekman, R.; Roepstorff, P. Sample purification and preparation technique based on nano-scale reversed-phase columns for the sensitive analysis of complex peptide mixtures by matrix-assisted laser desorption/ ionization mass spectrometry. J. Mass Spectrom. 1999, 34, 105−16. (49) Larsen, M. R.; Thingholm, T. E.; Jensen, O. N.; Roepstorff, P.; Jorgensen, T. J. D. Highly selective enrichment of phosphorylated peptides from peptide mixtures using titanium dioxide microcolumns. Mol. Cell. Proteomics 2005, 4, 873−86. (50) Jensen, S. S.; Larsen, M. R. Evaluation of the impact of some experimental procedures on different phosphopeptide enrichment techniques. Rapid Commun. Mass Spectrom. 2007, 21, 3635−45. (51) Schroeder, M. J.; Shabanowitz, J.; Schwartz, J. C.; Hunt, D. F.; Coon, J. J. A neutral loss activation method for improved phosphopeptide sequence analysis by quadrupole ion trap mass spectrometry. Anal. Chem. 2004, 76, 3590−8. (52) Engholm-Keller, K.; Hansen, T. A.; Palmisano, G.; Larsen, M. R. Multidimensional strategy for sensitive phosphoproteomics incorporating protein prefractionation combined with SIMAC, HILIC, and TiO2 chromatography applied to proximal EGF signaling. J. Proteome Res. 2011, 10, 5383−97. (53) Olsen, J. V.; Mann, M. Improved peptide identification in proteomics by two consecutive stages of mass spectrometric fragmentation. Proc. Natl. Acad. Sci. U.S.A. 2004, 101, 13417−22. (54) Olsen, J. V.; Blagoev, B.; Gnad, F.; Macek, B.; Kumar, C.; Mortensen, P.; Mann, M. Global, in vivo, and site-specific phosphorylation dynamics in signaling networks. Cell 2006, 127, 635−48. (55) Pohl, R. J.; Fouts, J. R. A rapid method for assaying the metabolism of 7-ethoxyresorufin by microsomal subcellular fractions. Anal. Biochem. 1980, 107, 150−5. (56) Schwirzer, S. M.; Hofmaier, A. M.; Kettrup, A.; Nerdinger, P. E.; Schramm, K. W.; Thoma, H.; Wegenke, M.; Wiebel, F. J. Establishment of a simple cleanup procedure and bioassay for determining 2,3,7,8-tetrachlorodibenzo-p-dioxin toxicity equivalents of environmental samples. Ecotoxicol. Environ. Saf. 1998, 41, 77−82. 880

dx.doi.org/10.1021/pr3009429 | J. Proteome Res. 2013, 12, 866−882

Journal of Proteome Research

Article

dissociation tandem mass spectrometry. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 2199−204. (76) Dephoure, N.; Zhou, C.; Villén, J.; Beausoleil, S. A.; Bakalarski, C. E.; Elledge, S. J.; Gygi, S. P. A quantitative atlas of mitotic phosphorylation. Proc. Natl. Acad. Sci. U.S.A. 2008, 105, 10762−7. (77) Gauci, S.; Helbig, A. O.; Slijper, M.; Krijgsveld, J.; Heck, A. J. R.; Mohammed, S. Lys-N and trypsin cover complementary parts of the phosphoproteome in a refined SCX-based approach. Anal. Chem. 2009, 81, 4493−501. (78) Mayya, V.; Lundgren, D. H.; Hwang, S.-I.; Rezaul, K.; Wu, L.; Eng, J. K.; Rodionov, V.; Han, D. K. Quantitative phosphoproteomic analysis of T cell receptor signaling reveals system-wide modulation of protein−protein interactions. Sci. Signal. 2009, 2, ra46. (79) Abrink, M.; Ortiz, J. A.; Mark, C.; Sanchez, C.; Looman, C.; Hellman, L.; Chambon, P.; Losson, R. Conserved interaction between distinct Krüppel-associated box domains and the transcriptional intermediary factor 1 beta. Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 1422−6. (80) Chang, C.-W.; Chou, H.-Y.; Lin, Y.-S.; Huang, K.-H.; Chang, C.J.; Hsu, T.-C.; Lee, S.-C. Phosphorylation at Ser473 regulates heterochromatin protein 1 binding and corepressor function of TIF1beta/KAP1. BMC Mol. Biol. 2008, 9, 61. (81) Weiss, C.; Faust, D.; Schreck, I.; Ruff, A.; Farwerck, T.; Melenberg, A.; Schneider, S.; Oesch-Bartlomowicz, B.; Zatloukalová, J.; Vondrácek, J.; Oesch, F.; Dietrich, C. TCDD deregulates contact inhibition in rat liver oval cells via Ah receptor, JunD and cyclin A. Oncogene 2008, 27, 2198−207. (82) Dietrich, C.; Faust, D.; Budt, S.; Moskwa, M.; Kunz, A.; Bock, K.-W.; Oesch, F. 2,3,7,8-tetrachlorodibenzo-p-dioxin-dependent release from contact inhibition in WB-F344 cells: Involvement of cyclin A. Toxicol. Appl. Pharmacol. 2002, 183, 117−26. (83) Wiebel, F. J.; Klose, U.; Kiefer, F. Toxicity of 2,3,7,8tetrachlorodibenzo-p-dioxin in vitro: H4IIEC3-derived 5L hepatoma cells as a model system. Toxicol. Lett. 1991, 55, 161−9. (84) Sarioglu, H.; Brandner, S.; Jacobsen, C.; Meindl, T.; Schmidt, A.; Kellermann, J.; Lottspeich, F.; Andrae, U. Quantitative analysis of 2,3,7,8-tetrachlorodibenzo-p-dioxin-induced proteome alterations in 5L rat hepatoma cells using isotope-coded protein labels. Proteomics 2006, 6, 2407−21. (85) Underhill, C.; Qutob, M. S.; Yee, S. P.; Torchia, J. A novel nuclear receptor corepressor complex, N-CoR, contains components of the mammalian SWI/SNF complex and the corepressor KAP-1. J. Biol. Chem. 2000, 275, 40463−70. (86) Daub, H.; Olsen, J. V.; Bairlein, M.; Gnad, F.; Oppermann, F. S.; Korner, R.; Greff, Z.; Keri, G.; Stemmann, O.; Mann, M. Kinaseselective enrichment enables quantitative phosphoproteomics of the kinome across the cell cycle. Mol. Cell 2008, 31, 438−48. (87) Imami, K.; Sugiyama, N.; Kyono, Y.; Tomita, M.; Ishihama, Y. Automated phosphoproteome analysis for cultured cancer cells by two-dimensional nanoLC-MS using a calcined titania/C18 biphasic column. Anal. Sci. 2008, 24, 161−6. (88) Peterlin, B. M.; Price, D. H. Controlling the elongation phase of transcription with P-TEFb. Mol. Cell 2006, 23, 297−305. (89) Narita, T.; Yung, T. M. C.; Yamamoto, J.; Tsuboi, Y.; Tanabe, H.; Tanaka, K.; Yamaguchi, Y.; Handa, H. NELF interacts with CBC and participates in 3′ end processing of replication-dependent histone mRNAs. Mol. Cell 2007, 26, 349−65. (90) Tian, Y.; Ke, S.; Chen, M.; Sheng, T. Interactions between the aryl hydrocarbon receptor and P-TEFb. Sequential recruitment of transcription factors and differential phosphorylation of C-terminal domain of RNA polymerase II at cyp1a1 promoter. J. Biol. Chem. 2003, 278, 44041−8. (91) Geyer, M.; Wittinghofer, A. GEFs, GAPs, GDIs, and effectors: Taking a closer (3D) look at the regulation of Ras-related GTPbinding proteins. Curr. Opin. Struct. Biol. 1997, 7, 786−92. (92) Sarioglu, H.; Brandner, S.; Haberger, M.; Jacobsen, C.; Lichtmannegger, J.; Wormke, M.; Andrae, U. Analysis of 2,3,7,8tetrachlorodibenzo-p-dioxin-induced proteome changes in 5L rat hepatoma cells reveals novel targets of dioxin action including the

(57) Van Hoof, D.; Munoz, J.; Braam, S. R.; Pinkse, M. W. H.; Linding, R.; Heck, A. J. R.; Mummery, C. L.; Krijgsveld, J. Phosphorylation dynamics during early differentiation of human embryonic stem cells. Cell Stem Cell 2009, 5, 214−26. (58) Storey, J. D.; Tibshirani, R. Statistical significance for genomewide studies. Proc. Natl. Acad. Sci. U.S.A. 2003, 100, 9440−5. (59) Haase, H.; Podzuweit, T.; Lutsch, G.; Hohaus, A.; Kostka, S.; Lindschau, C.; Kott, M.; Kraft, R.; Morano, I. Signaling from betaadrenoceptor to L-type calcium channel: identification of a novel cardiac protein kinase A target possessing similarities to AHNAK. FASEB J. 1999, 13, 2161−72. (60) Hohaus, A.; Person, V.; Behlke, J.; Schaper, J.; Morano, I.; Haase, H. The carboxyl-terminal region of ahnak provides a link between cardiac L-type Ca 2+ channels and the actin-based cytoskeleton. FASEB J. 2002, 16, 1205−16. (61) Alvarez, J.; Hamplova, J.; Hohaus, A.; Morano, I.; Haase, H.; Vassort, G. Calcium current in rat cardiomyocytes is modulated by the carboxyl-terminal ahnak domain. J. Biol. Chem. 2004, 279, 12456−61. (62) Matza, D.; Badou, A.; Kobayashi, K. S.; Goldsmith-Pestana, K.; Masuda, Y.; Komuro, A.; McMahon-Pratt, D.; Marchesi, V. T.; Flavell, R. A. A scaffold protein, AHNAK1, is required for calcium signaling during T cell activation. Immunity 2008, 28, 64−74. (63) Bootman, M. D.; Collins, T. J.; Mackenzie, L.; Roderick, H. L.; Berridge, M. J.; Peppiatt, C. M. 2-Aminoethoxydiphenyl borate (2APB) is a reliable blocker of store-operated Ca2+ entry but an inconsistent inhibitor of InsP3-induced Ca2+ release. FASEB J. 2002, 16, 1145−50. (64) Kohle, C.; Gschaidmeier, H.; Lauth, D.; Topell, S.; Zitzer, H.; Bock, K. W. 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD)-mediated membrane translocation of c-Src protein kinase in liver WB-F344 cells. Arch. Toxicol. 1999, 73, 152−8. (65) Thomas, S. M.; Brugge, J. S. Cellular functions regulated by Src family kinases. Annu. Rev. Cell. Dev. Biol. 1997, 13, 513−609. (66) Ruff, S. J.; Chen, K.; Cohen, S. Peroxovanadate induces tyrosine phosphorylation of multiple signaling proteins in mouse liver and kidney. J. Biol. Chem. 1997, 272, 1263−7. (67) Kewley, R. J.; Whitelaw, M. L. Phosphorylation inhibits DNAbinding of alternatively spliced aryl hydrocarbon receptor nuclear translocator. Biochem. Biophys. Res. Commun. 2005, 338, 660−7. (68) Sogawa, K.; Nakano, R.; Kobayashi, A.; Kikuchi, Y.; Ohe, N.; Matsushita, N.; Fujii-Kuriyama, Y. Possible function of Ah receptor nuclear translocator (Arnt) homodimer in transcriptional regulation. Proc. Natl. Acad. Sci. U.S.A. 1995, 92, 1936−40. (69) Huffman, J. L.; Mokashi, A.; Bächinger, H. P.; Brennan, R. G. The basic helix-loop-helix domain of the aryl hydrocarbon receptor nuclear transporter (ARNT) can oligomerize and bind E-box DNA specifically. J. Biol. Chem. 2001, 276, 40537−44. (70) Whitelaw, M.; Pongratz, I.; Wilhelmsson, A.; Gustafsson, J. A.; Poellinger, L. Ligand-dependent recruitment of the Arnt coregulator determines DNA recognition by the dioxin receptor. Mol. Cell. Biol. 1993, 13, 2504−14. (71) Xue, Y.; Ren, J.; Gao, X.; Jin, C.; Wen, L.; Yao, X. GPS 2.0, a tool to predict kinase-specific phosphorylation sites in hierarchy. Mol. Cell. Proteomics 2008, 7, 1598−1608. (72) Beausoleil, S. A.; Jedrychowski, M.; Schwartz, D.; Elias, J. E.; Villén, J.; Li, J.; Cohn, M. A.; Cantley, L. C.; Gygi, S. P. Large-scale characterization of HeLa cell nuclear phosphoproteins. Proc. Natl. Acad. Sci. U.S.A. 2004, 101, 12130−5. (73) Kim, J.-E.; Tannenbaum, S. R.; White, F. M. Global phosphoproteome of HT-29 human colon adenocarcinoma cells. J. Proteome Res. 2005, 4, 1339−46. (74) Yu, L.-R.; Zhu, Z.; Chan, K. C.; Issaq, H. J.; Dimitrov, D. S.; Veenstra, T. D. Improved titanium dioxide enrichment of phosphopeptides from HeLa cells and high confident phosphopeptide identification by cross-validation of MS/MS and MS/MS/MS spectra. J. Proteome Res. 2007, 6, 4150−62. (75) Molina, H.; Horn, D. M.; Tang, N.; Mathivanan, S.; Pandey, A. Global proteomic profiling of phosphopeptides using electron transfer 881

dx.doi.org/10.1021/pr3009429 | J. Proteome Res. 2013, 12, 866−882

Journal of Proteome Research

Article

mitochondrial apoptosis regulator VDAC2. Mol. Cell. Proteomics 2008, 7, 394−410. (93) Qi, H.; Fournier, A.; Grenier, J.; Fillion, C.; Labrie, Y.; Labrie, C. Isolation of the novel human guanine nucleotide exchange factor Src homology 3 domain-containing guanine nucleotide exchange factor (SGEF) and of C-terminal SGEF, an N-terminally truncated form of SGEF, the expression of which is regulated by androgen in prostate cancer cells. Endocrinology 2003, 144, 1742−52. (94) Ellerbroek, S. M.; Wennerberg, K.; Arthur, W. T.; Dunty, J. M.; Bowman, D. R.; DeMali, K. A.; Der, C.; Burridge, K. SGEF, a RhoG guanine nucleotide exchange factor that stimulates macropinocytosis. Mol. Biol. Cell 2004, 15, 3309−19. (95) Zerial, M.; McBride, H. Rab proteins as membrane organizers. Nat. Rev. Mol. Cell. Biol. 2001, 2, 107−17. (96) Zougman, A.; Wiśniewski, J. R. Beyond linker histones and high mobility group proteins: global profiling of perchloric acid soluble proteins. J. Proteome Res. 2006, 5, 925−34. (97) Laughlin, K. M.; Luo, D.; Liu, C.; Shaw, G.; Warrington, J.; Kenneth, H.; Law, B. K.; Harrison, J. K. Hematopoietic- and neurologic-expressed sequence 1 (Hn1) depletion in B16.F10 melanoma cells promotes a differentiated phenotype that includes increased melanogenesis and cell cycle arrest. Differentiation 2009, 78, 35−44. (98) Lu, K. H.; Patterson, A. P.; Wang, L.; Marquez, R. T.; Atkinson, E. N.; Baggerly, K. A.; Ramoth, L. R.; Rosen, D. G.; Liu, J.; Hellstrom, I.; Smith, D.; Hartmann, L.; Fishman, D.; Berchuck, A.; Schmandt, R.; Whitaker, R.; Gershenson, D. M.; Mills, G. B.; Bast, J.; Robert, C. Selection of potential markers for epithelial ovarian cancer with gene expression arrays and recursive descent partition analysis. Clin. Cancer Res. 2004, 10, 3291−300. (99) Miserey-Lenkei, S.; Waharte, F.; Boulet, A.; Cuif, M.-H.; Tenza, D.; El Marjou, A.; Raposo, G.; Salamero, J.; Héliot, L.; Goud, B.; Monier, S. Rab6-interacting protein 1 links Rab6 and Rab11 function. Traffic 2007, 8, 1385−403. (100) Echard, A.; Opdam, F. J.; de Leeuw, H. J.; Jollivet, F.; Savelkoul, P.; Hendriks, W.; Voorberg, J.; Goud, B.; Fransen, J. A. Alternative splicing of the human Rab6A gene generates two close but functionally different isoforms. Mol. Biol. Cell 2000, 11, 3819−33. (101) Cuif, M. H.; Possmayer, F.; Zander, H.; Bordes, N.; Jollivet, F.; Couedel-Courteille, A.; Janoueix-Lerosey, I.; Langsley, G.; Bornens, M.; Goud, B. Characterization of GAPCenA, a GTPase activating protein for Rab6, part of which associates with the centrosome. EMBO J. 1999, 18, 1772−82. (102) Alexandru, G.; Graumann, J.; Smith, G. T.; Kolawa, N. J.; Fang, R.; Deshaies, R. J. UBXD7 binds multiple ubiquitin ligases and implicates p97 in HIF1alpha turnover. Cell 2008, 134, 804−16. (103) Ye, Y. Diverse functions with a common regulator: ubiquitin takes command of an AAA ATPase. J. Struct. Biol. 2006, 156, 29−40. (104) Lim, P. J.; Danner, R.; Liang, J.; Doong, H.; Harman, C.; Srinivasan, D.; Rothenberg, C.; Wang, H.; Ye, Y.; Fang, S.; Monteiro, M. J. Ubiquilin and p97/VCP bind erasin, forming a complex involved in ERAD. J. Cell Biol. 2009, 187, 201−17. (105) Seifert, A.; Katschinski, D. M.; Tonack, S.; Fischer, B.; Navarrete Santos, A. Significance of prolyl hydroxylase 2 in the interference of aryl hydrocarbon receptor and hypoxia-inducible factor1 alpha signaling. Chem. Res. Toxicol. 2008, 21, 341−8. (106) Moritz, A.; Li, Y.; Guo, A.; Villén, J.; Wang, Y.; MacNeill, J.; Kornhauser, J.; Sprott, K.; Zhou, J.; Possemato, A.; Ren, J. M.; Hornbeck, P.; Cantley, L. C.; Gygi, S. P.; Rush, J.; Comb, M. J. AktRSK-S6 kinase signaling networks activated by oncogenic receptor tyrosine kinases. Sci. Signal. 2010, 3, ra64. (107) Matza, D.; Flavell, R. A. Roles of Ca(v) channels and AHNAK1 in T cells: the beauty and the beast. Immunol. Rev. 2009, 231, 257−64. (108) Shao, Y.; Czymmek, K. J.; Jones, P. A.; Fomin, V. P.; Akanbi, K.; Duncan, R. L.; Farach-Carson, M. C. Dynamic interactions between L-type voltage-sensitive calcium channel Cav1.2 subunits and ahnak in osteoblastic cells. Am. J. Physiol. Cell. Physiol. 2009, 296, C1067−78.

(109) Benaud, C.; Gentil, B. J.; Assard, N.; Court, M.; Garin, J.; Delphin, C.; Baudier, J. AHNAK interaction with the annexin 2/ S100A10 complex regulates cell membrane cytoarchitecture. J. Cell Biol. 2004, 164, 133−44. (110) Komuro, A.; Masuda, Y.; Kobayashi, K.; Babbitt, R.; Gunel, M.; Flavell, R. A.; Marchesi, V. T. The AHNAKs are a class of giant propeller-like proteins that associate with calcium channel proteins of cardiomyocytes and other cells. Proc. Natl. Acad. Sci. U.S.A. 2004, 101, 4053−8. (111) Kim, E.; Sheng, M. PDZ domain proteins of synapses. Nat. Rev. Neurosci. 2004, 5, 771−81. (112) Emdad, L.; Lee, S.-G.; Su, Z. Z.; Jeon, H. Y.; Boukerche, H.; Sarkar, D.; Fisher, P. B. Astrocyte elevated gene-1 (AEG-1) functions as an oncogene and regulates angiogenesis. Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 21300−5. (113) Sarkar, D.; Emdad, L.; Lee, S.-G.; Yoo, B. K.; Su, Z.-Z.; Fisher, P. B. Astrocyte elevated gene-1: far more than just a gene regulated in astrocytes. Cancer Res. 2009, 69, 8529−35. (114) Villen, J.; Beausoleil, S. A.; Gerber, S. A.; Gygi, S. P. Large-scale phosphorylation analysis of mouse liver. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 1488−93. (115) Zhou, H.; Ye, M.; Dong, J.; Han, G.; Jiang, X.; Wu, R.; Zou, H. Specific phosphopeptide enrichment with immobilized titanium ion affinity chromatography adsorbent for phosphoproteome analysis. J. Proteome Res. 2008, 7, 3957−67. (116) Zahedi, R. P.; Lewandrowski, U.; Wiesner, J.; Wortelkamp, S.; Moebius, J.; Schütz, C.; Walter, U.; Gambaryan, S.; Sickmann, A. Phosphoproteome of resting human platelets. J. Proteome Res. 2008, 7, 526−34. (117) Hu, G.; Wei, Y.; Kang, Y. The multifaceted role of MTDH/ AEG-1 in cancer progression. Clin. Cancer Res. 2009, 15, 5615−20.

882

dx.doi.org/10.1021/pr3009429 | J. Proteome Res. 2013, 12, 866−882