Proteomic analysis of mismatch repair-mediated alkylating agent-induced DNA damage response
© Chen et al.; licensee BioMed Central Ltd. 2013
Received: 9 July 2013
Accepted: 26 August 2013
Published: 19 September 2013
Mediating DNA damage-induced apoptosis is an important genome-maintenance function of the mismatch repair (MMR) system. Defects in MMR not only cause carcinogenesis, but also render cancer cells highly resistant to chemotherapeutics, including alkylating agents. To understand the mechanisms of MMR-mediated apoptosis and MMR-deficiency-caused drug resistance, we analyze a model alkylating agent (N-methyl-N’-nitro-N-nitrosoguanidine, MNNG)-induced changes in protein phosphorylation and abundance in two cell lines, the MMR-proficient TK6 and its derivative MMR-deficient MT1.
Under an experimental condition that MNNG-induced apoptosis was only observed in MutSα-proficient (TK6), but not in MutSα-deficient (MT1) cells, quantitative analysis of the proteomic data revealed differential expression and phosphorylation of numerous individual proteins and clusters of protein kinase substrates, as well differential activation of response pathways/networks in MNNG-treated TK6 and MT1 cells. Many alterations in TK6 cells are in favor of turning on the apoptotic machinery, while many of those in MT1 cells are to promote cell proliferation and anti-apoptosis.
Our work provides novel molecular insights into the mechanism of MMR-mediated DNA damage-induced apoptosis.
DNA mismatch repair (MMR) ensures genetic stability by correcting biosynthetic errors, suppressing non-homologous recombination, and mediating DNA damage-induced apoptosis [1–5]. While the first two functions of the MMR system prevent mismatch-derived mutations, the apoptotic function of the system is to initiate programmed cell death in cells with mutagenic and carcinogenic DNA lesions. The latter function has been recognized as an important factor in cancer chemotherapy, because tumor cells defective in MMR are highly resistant to chemotherapeutic drugs, including alkylating agents .
A key player in mammalian MMR is the major mismatch recognition protein MutSα, which is composed of MSH2 and MSH6 subunits. In a concerted action with other essential MMR factors that include MutLα, proliferating cell nuclear antigen (PCNA), replication protein A (RPA), exonuclease I, and DNA polymerases, MutSα targets the newly synthesized DNA strand for mismatch removal [2, 7]. Previous studies have revealed that besides mismatch recognition, MutSα also recognizes base pairs that contain a variety of modified or damaged bases, including alkylating adduct O6-methylguanine (O6-MeG) [6, 8, 9]. It has been suggested that recognition of the O6-MeG:T or O6-MeG:C pair by MutSα possibly in complex with MutLα, activates the ATR/ATM signaling network to induce apoptosis, via either the futile cycle mechanism or the direct interaction and signaling mechanism [1, 2, 6, 8, 10]. Many downstream targets, such as p53, Chk1, Chk2, CDC25A, and SMC1, have been shown to play important roles in the MMR-dependent alkylation-induced apoptosis [8, 10–12]. However, the molecular basis of the process, particularly how the DNA damage signal in the nucleus is transmitted to the cytoplasm/mitochondria to initiate apoptosis is unknown.
Using quantitative global proteomics and phosphoproteomics approaches, we have investigated in this study a pair of MMR-proficient and -deficient human lymphoblastoid cell lines, TK6 and MT1, for their responses to the treatment of N-methyl-N′-nitro-N-nitrosoguanidine (MNNG), a widely used model chemical for alkylating agent-induced mutagenesis, carcinogenesis, and cell killing [13, 14]. TK6 and MT1 are two closely-related cell lines with dramatically different susceptibility to MNNG-induced DNA damage. Both cells do not express methylguanine methyl transferase (MGMT), a direct reversal repair enzyme that protects cells from MMNG-induced cytotoxicity by transferring the methyl group from O6-MeG to an internal cysteine of the enzyme, restoring Watson-Crick base pairing at the site of the damage . MT1 cells were derived from MMR-proficient TK6 cells by mutagenesis , and are defective in the MSH6 subunit of MutSα, thereby defective in MMR [17, 18]. As a result, MT1 cells are 500-times more resistant to killing by MNNG cytotoxicity than their parental TK6 cells. The differential response of TK6 and MT1 cells to MNNG underscores the critical role played by MMR in alkylating agent-induced apoptosis/DNA damage response (DDR) and makes these cells the ideal subject for this study.
We demonstrate here that TK6 and MT1 cells display significant differences in protein phosphorylation species, stoichiometry and abundance, and in activation of protein kinases and signaling networks in response to MNNG treatment. Generally, MNNG-induced alterations in TK6 cells promote apoptosis, but changes in MT1 cells stimulate cell proliferation. Our data provide significant new insights into the mechanism of MMR-mediated DNA damage-induced apoptosis.
Results and discussion
Cell cycle progression in MNNG-treated TK6 and MT1 cells
Overview of changes in total and nuclear proteomes of MNNG-treated TK6 and MT1 cells
Number of peptides and proteins identified and quantified in MNNG-treated cells
Whole cell extract
1438 (FDR# = 1.9%)
1654 (FDR# = 1.6%)
1594 (FDR# = 2.7%)
1348 (FDR# = 2.3%)
1792 (FDR# = 4.9%)
1765 (FDR# = 4.5%)
1238 (FDR# = 4.5%)
1102 (FDR# = 4.8%)
Because significantly more up- or down- regulatory events were documented in the nuclear extract (Table 1), this directed our interest to the nuclear proteome, which is considered relevant to the MNNG-induced DDR. Therefore, we conducted detail analysis dissecting the key differences in the nuclear proteomes in MNNG-treated TK6 and MT1 cells.
System-level differences between the nuclear proteomes in MNNG-treated TK6 and MT1 cells
To identify key differences between the nuclear proteomes in MNNG-treated TK6 and MT1 cells, we performed quantitative proteomics analyses and obtained nuclear proteomics profiles of the two cells (see Additional file 1: Table S1 and Additional file 2: Table S2). Systematical comparison of the nuclear proteomes allowed us to extract out all the up- or down-regulated proteins and phosphoproteins from MNNG-treated TK6 and MT1 cells. We further analyzed for the extracted data for their biological function distribution using Gene Ontology (GO) annotation, and discovered that more than 90% of the up- and down-regulated proteins and phosphoproteins could be categorized into ten GO categories. These GO categories include: Transcription, DNA damage and repair, DNA replication, Cell cycle, RNA processing and splicing, Chromatin modification and nucleosome assembly, Translation, Transport, Apoptosis and Protein ubiquitination. The number of up- and down-regulated proteins in MNNG-treated TK6 and MT1 cells that fall into the GO categories is shown in Figure 2C and 2D (details are listed in Additional file 3: Table S3), which reveals the following interesting phenomena.
First, we found that there were essentially more protein expression changes, particularly up-regulated proteins (see red bars in Figure 2C), in MT1 cells than in TK6 cells (194 in MT1 vs. 74 in TK6), but there were more protein phosphorylation level changes, especially down-regulated phosphosites (see green bars in Figure 2D) in TK6 cells than in MT1 cells (203 in TK6 vs. 111 in MT1). Interestingly, the changes in protein expression (Figure 2C) and protein phosphorylation (Figure 2D) were independent of each other. For example, although more changes in protein expression were observed for MT1 Transcription category (Figure 2C), more protein phosphorylation level changes were observed in TK6 cells (Figure 2D).
Second, our data show strong up-regulation of “Translation proteins” in MNNG-treated MT1 cells (43 proteins), while significantly fewer (9 proteins) “Translation proteins” were up-regulated in MNNG-treated TK6 cells (Figure 2C). If this correlates with more active translation in MNNG-treated MT1 cells than in MNNG-treated TK6 cells, then it could explain the higher total number of up-regulated nuclear proteins in these cells (Table 1). We hypothesize that the preferentially up-regulated “Translation proteins” in MNNG-treated MT1 cells may be responsible for tolerance to killing by MNNG toxicity. Indeed, some of the up-regulated proteins that we identified in MT1 cells, including EF2 (Elongation factor 2), EIF4G1 (Eukaryotic translation initiation factor 4 gamma 1), RPL5 (60S ribosomal protein L5) and RPL23A (60S ribosomal protein L23a) (as listed in Additional file 1: Table S1), have been previously implicated in anti-apoptosis, drug resistance and cell survival [21–25].
Third, from the protein phosphorylation perspective, in nearly all the GO categories outlined, we identified more alteration in MNNG-treated TK6 cells than in MNNG-treated MT1 cells (Figure 2D). This is especially evident in the GO categories of “Transcription” (78 vs. 31), “DNA damage and repair” (24 vs. 8), “Cell cycle” (39 vs. 14), and “Transport” (20 vs. 9) (Figure 2D). Because activities in protein phosphorylation and dephosphorylation in general correlated with activities in cellular signaling, we hypothesized what we observed may reflected a functional MMR pathway in TK6 directly or indirectly regulated DNA damage-related protein subnetworks during MNNG-stimulation, suggesting that functional MMR or functional MutSα might play an important role in modulating a series of biological processes, and ultimately leading to apoptosis in TK6 cells (Figure 1A and 1B).
Unique proteomic changes in MNNG-treated TK6 and MT1 cells
In the heat map (Figure 3), the color indicated the level of phosphorylation at individual protein phosphorylation sites, with the green and red colors being the lowest and highest levels, respectively (see color diagram at the top panel). In each of the GO categories shown in Figure 3, protein phosphorylation sites are listed according to their phosphorylation levels (from low to high) in TK6. From this analysis, we deduced that: 1) the direction of change in protein phosphorylation was rarely opposite in TK6 and MT1 cells; 2) a few phosphorylation were up- or down-regulated in both MNNG-treated cell lines; but 3) most phosphorylation were up- or down-regulated in one of the cell lines and unaffected by MNNG-treatment in the other cell line.
Among the 133 phosphorylation sites listed in Figure 3, only one phosphorylation site, namely Ser720 of AKAP2 shown opposite change in protein phosphorylation (see the first entry in the Unknown Miscellaneous category). AKAP2 protein can bind to regulatory subunit (RII) of protein kinase A  and may be involved in establishing polarity in signaling systems. The biological significance of AKAP2 protein in differential responses of TK6 and MT1 toward MNNG-treatment remains to be investigated.
It is interesting to note that while we detected phosphorylation up-regulation of some important DDR factors in TK6 cells, we also identified more phosphorylation up-regulation of components involved in apoptosis in MNNG-tolerant MT1 cells than in MNNG-sensitive TK6 cells (see Figure 3, the GO category of “Apoptosis”). These proteins include anti-apoptotic ERCC5 (Ser384) and pro-apoptotic VDAC1 (Ser104), LAP4 (Ser1448), SH3K1 (Ser230), and SHIP1 (Ser971). This phenomenon appears to suggest that MNNG-induced phosphorylations have greatly altered the biological functions of these apoptosis-related proteins, i.e., inhibiting the pro-apoptotic activities of VDAC1, LAP4, SH3K1 and SHIP1, but enhancing the anti-apoptotic activity of ERCC5. However, further investigations are required to verify these possibilities.
Activation of nucleocytoplasmic shuttling in TK6 cells upon exposure to MNNG
The comparative mass spectrometry analysis of phosphoprotein between nuclear extract and whole cell extract allowed us to identify an important phenomenon associated with MNNG treated TK6 cells, i.e., activation of nucleocytoplasmic transport. As shown in Figure 4B, phosphorylation of Ser428 in Ran GTPase-activating protein 1 (RanGAP1), a critical regulator of the small GTPase Ran (RAN)-dependent nucleoplasmic transport [30–32], was up-regulated 1.8 fold in the nuclear proteome of MNNG-treated TK6 cells, but the up-regulation was not detected in MNNG-treated MT1 cells (upper spectra). Because RanGAP1 is located primarily within the cytoplasm  and because Ser428 phosphorylation was observed previously in an MNNG-independent manner , we suspected that the up-regulated RanGAP1 level in the TK6 nucleus is likely due to transport of the protein from the cytoplasm. Indeed, we observed 1.8-fold up-regulation of a reference peptide of RanGAP1 in the TK6 nuclear extract (Figure 4B, lower spectrum), but not in the TK6 whole cell extract (data not shown). These results strongly suggest that the up-regulated nuclear RanGAP1 in MNNG-treated TK6 cells is a result of translocation of the protein from cytoplasm. Similarly, we detected nucleocytoplasmic shuttling of RAN, an essential player of the RAN-dependent nucleocytopasmic transport [30, 31]. Since this phenomenon was only observed in MNNG-treated TK6 cells, the activation of the RAN/RanGAP1-dependent nucleocytoplasmic shuttling could be part of the MMR-mediated signaling network. Further studies are required to explore this possibility.
A DDR protein interaction network in TK6 cells
When comparing MNNG-treated nuclear phosphoproteome in TK6 and MT1 cells, we noticed significantly more protein phosphorylation level changes in TK6 cells than in MT1 cells (203 in TK6 vs. 111 in MT1). Are these regulations on protein phosphorylation related to a DDR protein interaction network? To address this question, all phospho- and nonphospho-proteins identified in TK6 (as listed in Additional file 1: Table S1 and Additional file 2: Table S2) were queried for protein interactions using the program STRING (version 8.2) (http://string-db.org/) . STRING integrates known and predicted interactions from multiple sources to generate a putative interaction network map for a given data set. A network of 1285 proteins (nodes) and 8427 connections (edges) were obtained (Additional file 6: Figure S2).
ATM (ataxia telangiectasia mutated) and ATR (ATM and Rad3-related), members of the phosphoinositide 3-kinase related kinase (PIKK) family, are well-recognized master regulators of DDR. Together with other regulator and mediator proteins, ATM and ATR stimulate phosphorylation of Chk1 and Chk2, which in turn activate signaling cascades via a large number of downstream effector proteins, such as MCMs, RFC, TopBP1, Rad51 and FANCD2 . Our proteomics data also showed increased expression or phosphorylation of BRCA1, NBN, RAD50, TP53BP1, SMC3, TOP2A, RBBP8, LIG3 and TERF2IP, all members of the DNA damage and repair subnetwork (Figure 5A). Previous studies suggested that NBN plays a role in the cellular response to MNNG , and that phosphorylation of SMC3, RAD50 and TP53BP1 increased in IR-treated cells . Together, these data help validate our proteomics approach and analysis. However, a number of proteins with known roles in DDR, such as Chk1, Chk2, p53, RAD17, XRCC2, ERCC1 and MDC1, were not identified as DDR network components in this study. This may indicate that these low abundance proteins fell below the detection limit associated with the instrumentation and analytical procedures used in this study. On the other hand, we identified a number of novel putative DDR proteins events, which have not previously been linked to DDR (see below).
Prediction of DDR protein kinase networks in MNNG-treated TK6 cells
Our quantitative phosphoproteomic analysis provided information regarding protein phosphorylation regulation upon MNNG-treatment. In order to use our data to explore protein kinases that may participate in DDR network, we used the NetworKIN algorithm to predict which kinases were responsible for MNNG-induced phosphorylation events within the DNA damage and repair, Apoptosis, Cell cycle and RNA processing and splicing subnetworks in TK6 cells (Figure 5B and Additional file 7: Figure S3).
We identified ATM or ATR as projected kinases for several proteins with up-regulation in phosphorylation. Because these protein substrates, NBN (Ser343), RAD50 (Ser635), TP53B (Ser580, Ser831) and SMC3 (Ser1083), are known players, and ATM/ATR are known master regulators in DDR, this finding further validated our phosphoproteomics approach in study cellular responses to MNNG treatment.
In addition, our analysis also project some other kinases, such as CDK2 (a cyclin-dependent protein kinase), Casein kinase II, and MAP kinase may be responsible for the phosphorylation of many protein substrates participating in “DNA damage and repair”, “Apoptosis”, “Cell cycle” and “RNA processing and splicing” in MNNG-treated TK6 cells (Figure 5B and Additional file 7: Figure S3), indicating that these kinases may also play important regulatory roles in MNNG-induced DDR.
Some proteins may be targeted by more than one protein kinase in MNNG-treated TK6 cells, suggesting that DDR involves coordinated regulatory activities and cross-talk between multiple kinases in TK6 cells. For example, S635 in RAD50, up-regulated 6.3-fold, is an ATM/ATR target, while T690 in RAD50, up-regulated 5.5-fold, is a target of Casein kinase II (Figure 5B).
Although unbiased phosphoproteomics analysis and system biology-based network investigation can provide insight into cellular signaling, in many cases, traditional molecular biology-based study, combined with functional assays are needed to identify and screen for the biological significance of predicted phosphorylation. Many protein phosphorylation sites we discovered are consistent with previous knowledge accumulated over the years using molecular biology method, giving credibility to our system biology level approach. For example, Ser343 phosphorylation in NBN (a ATM/ATR kinase substrate, Figure 5B) was identified previously as closely related with NBN-regulated apoptosis, and a mutation in the ATM/ATR phosphorylation site (Ser343A) leads to a severe inhibition of apoptosis . The CDK2 kinase substrate RBBP8 (Ser327) (Figure 5B), a putative tumor suppressor, was recently shown to be involved in regulating the G2/M checkpoint through BRCA1-dependent ubiquitination in response to Ser327 phosphorylation  and to counteract Rb-mediated G1 restraint .
The large amount of information about regulations on protein phosphorylation generated from our study may provide new hints into mechanism of protein interaction. For example, Thr690 of RAD was identified as being up-regulated in MNNG-treated TK6 nuclear, with Casein Kinase II as its potential kinase (Figure 5B). RAD50 has a Cys-X-X-Cys motif of the zinc-hook structure located in the middle of the coiled-coil domain, which functions as a dimerization domain between two RAD50 arms . The RAD50 hook is very important in forming a complex of MRE11/RAD50/NBN (MRN complex), which acts a double-strand break sensor for ATM and recruits ATM to broken DNA molecules and activates ATM to initiate phosphorylation-network-based apoptosis [42, 43]. Interestingly, the Thr690 residue is located in the zinc-hook of RAD50. Therefore, it is speculated that the up-regulated phosphorylation on Thr690 site affects RAD50 dimerization, further influences the formation of MRN complex, leading to a series of downstream apoptosis-related protein interactions.
Cell culture, flow cytometry and DNA fragmentation analysis
TK6 and MT1 cell lines were cultured in RPMI 1640 medium (GIBCO) supplemented with 10% heat-inactivated undialyzed fetal bovine serum (Hyclone), 100 U of penicillin and 100 μg/mL of streptomycin (Invitrogen) in a humidified incubator with 5% CO2 at 37°C.
For flow cytometry analysis, cells were treated or mock-treated with 0.5 μM MNNG and incubated for the indicated time periods. A total of 1.0 × 106 cells were washed with ice-cold PBS and fixed in ice-cold 70% ethanol overnight. The cells were then treated with 100 μg/mL RNase A for 1 h at 37°C and stained with 50 μg/mL propidium iodide for 30 min at 4°C in the dark. Cell cycle analysis was performed using a BECKMAN COULTER EPICS XL flow cytometer and the EXPO32ADC software.
For DNA fragmentation analysis, cells were treated or mock-treated with 0.5 μM MNNG and incubated for the indicated time periods. A total of 5 × 105 cells were washed with ice-cold PBS, centrifuged for 5 min at 400 × g, and resuspended in 40 μL cell lysis buffer (20 mM EDTA, 100 mM Tris, pH8.0 and 0.8% (w/v) SDS). RNase A (2 μL of 1 mg/mL, Fermentas) was added and incubated for 120 min at 37°C, then 10 μL Proteinase K (20 mg/mL, TaKaRa) was added and incubated for at least 90 min at 50°C. Sample was mixed with 10 × DNA loading buffer (TaKaRa), and loaded onto a 2% agarose gel. Electrophoresis was conducted at 2–4 V/cm for approximately 6 h. The DNA was stained with 1 μg/mL ethidium bromide (Sigma) in TAE buffer.
SILAC labeling, MNNG treatment and western blot
Cells were washed twice with RPMI 1640 medium without Arg and Lys (Thermo Scientific Pierce) and reconstituted in RPMI 1640 medium containing either 12C6,14 N4 Arg, 12C6,14 N2 Lys (Sigma) and normal L-proline or 13C6,15 N4 Arg, 13C6,15 N2 Lys (Sigma) and normal L-proline. The media also contain 10% heat inactivated dialyzed fetal bovine serum (Thermo Scientific Pierce), 100 U of penicillin and 100 μg/mL of streptomycin. Cells were cultured in a humidified incubator (5% CO2) at 37°C for ~8 cell doublings. The concentration of Arg and Lys used in SILAC labeling of TK6 and MT1 cell lines was 0.398 mM and 0.798 mM respectively. Adding additional L-proline (200 mg/liter) should reduce arginine to proline conversion .
Cells were treated with or without 0.5 μM MNNG for 4 hours, and harvested by centrifugation at 400 g. Cell pellets were washed three times with cold PBS, and lysed for 30 min at 4°C in RIPA buffer (25 mM Tris–HCl pH 7.6, 150 mM NaCl, 1% NP-40, 1% sodium deoxycholate, 0.1% SDS) containing phosphatase inhibitor cocktail (PhosSTOP, Roche Applied Science) and protease inhibitor cocktail (COMPLETE, Roche Applied Science). After centrifugation for 15 minutes at 16000 g, the supernatant was collected and saved as whole cell protein extract.
For the immunoblotting, a small aliquot of proteins from MNNG-treated TK6 and MT1 cells were separated by SDS-PAGE, transferred electrophoretically onto a PVDF membrane (Amersham) and blocked for 1 h with TBST containing 5% nonfat milk. The PVDF membranes were incubated with the primary antibody overnight at 4°C, washed three times with TBST, incubated with a horseradish peroxidase-conjugated secondary antibody for 1 h at room temperature, and then developed with a chemiluminescent reagent (SuperSignal® West Pico Chemiluminescent Substrate, PIERCE). Primary antibodies were rabbit anti-human Phospho-(Ser/Thr) ATM/ATR substrate antibody (Cell Signaling Technology) or rabbit anti-human β-Actin antibody (Cell Signaling Technology).
Nuclear protein extraction
All nuclear-protein-extraction steps were conducted on ice or at 4°C. Cell pellets were re-suspended in 10 volumes of ice-cold buffer A (10 mM Tris–HCl (pH7.4), 5 mM MgCl2, 10 mM NaCl, 1 mM DTT, plus proteinase inhibitor cocktail (Roche), phosphatase inhibitor cocktail (Roche)) and incubated on ice for 20 min. Then ice-cold buffer B (10 mM Tris–HCl (pH7.4), 5 mM MgCl2, 10 mM NaCl, 1 mM DTT, proteinase inhibitor cocktail, phosphatase inhibitor cocktail, 10% NP-40) was added, adjusting the final concentration of NP-40 to 0.5%. After vigorous vortexing on the highest setting for 5 s, the tube was incubated on ice for 1 min, then vortexed for another 5 s. The extract was centrifuged for 5 min at 500 × g, and the supernatant fraction (cytoplasmic extract) removed and saved. The pellet was washed with 5 volumes of ice-cold buffer A to minimize cytoplamic protein contamination, then re-centrifuged for 5 minutes at 500 × g to obtain the nuclei as pellet. The nuclear pellet was resuspended in 10 volumes of freshly prepared buffer C (20 mM HEPES-KOH (pH7.9), 1.5 mM MgCl2, 0.5 M NaCl, 1 mM DTT, 0.2 mM EDTA, 20% (V/V) glycerol, proteinase inhibitor cocktail, phosphatase inhibitor cocktail). The nuclear proteins were solubilized by vortexing 3 times for 15 s at 10 min intervals, after which the extract was centrifuged at 16000 × g for 30 min. The supernatant (nuclear proteins extract) was immediately transferred to a pre-chilled tube. The purity of the nuclear and cytoplasmic protein extracts was evaluated by Western Blot for Histone H3 (as a nuclear protein marker) and GAPDH (as a cytoplasmic protein marker) by using rabbit anti-human GAPDH Antibody (Cell Signaling Technology) and rabbit anti-human Histone H3 Antibody (Cell Signaling Technology) as primary antibodies.
Protein digestion and phosphopeptide enrichment
Proteins were precipitated with 3 volumes 50% acetone/50% ethanol/0.1% acetic acid on ice for 1 h. After centrifugation at 10000 × g for 15 min, the protein pellet was resuspended in 8 M urea/0.2 M Tris, pH 8/4 mM CaCl2. Proteins were reduced with 10 mM DTT for 1 h at 56°C, alkylated with 40 mM iodoacetamide for 30 min at room temperature in the dark. After diluting sample 7-fold, trypsin (Worthington) was added in a 1:50 (trypsin/protein) w/w ratio, and incubated overnight at 37°C.
Tryptic peptides from whole cell or nuclear protein extracts were loaded onto a 2 g Sep-Pak C18 column (Waters), washed twice with 10 mL 1% acetic acid, eluted with 7 mL 80% acetonitrile containing 0.1% acetic acid, lyophilized using a speed-vac, resuspended in 400 μL 1% acetic acid and loaded onto a mini column of 40 μL IMAC resin (prepared as previously described [46, 47]). The IMAC minicolumn was washed twice with 40 μL wash buffer containing 25% acetonitrile, 100 mM NaCl and 0.1% acetic acid, then washed once each with 40 μL 1% acetic acid and 20 μL deionized water, eluted with 120 μL 6% NH4OH, and dried under vacuum. The IMAC column flow-through, which contains mainly non-phosphopeptides, was also collected.
Enriched phosphopeptides and non-phosphopeptides were fractionated separately using a TSKgel Amide-80 column (2.0 mm × 150 mm, 5 μm particle size, 200 Å pore size) (TOSOH Bioscience) on a Agilent 1200 system (Agilent Technologies). A 60 min elution gradient was used for phosphopeptide separation with 90% ACN, 0.005% TFA as Buffer A and 0.005% TFA as Buffer B. The gradient elution profile composed of 0%-12% B for 5 min, 12-30% B for 25 min, 30-90% B for 5 min, then maintained at 90% B for 5 min, followed by 10-100% A for 5 min, ending at 100% A for 15 min. The flow rate was 0.15 mL/min. UV absorbance was monitored at 215 nm. A total of 26 0.3 mL fractions were collected. Fractions were dried via vacuum. Non-phosphopeptides were fractionated in essentially the same manner, except that the gradient elution profile was 0%-5% B for 5 min and then 5-30% B for 25 min.
A QSTAR ELITE mass spectrometer (Applied Biosystems) was coupled with an online Eksigent nano MDLC systems utilizing a nanospray ionization source. Peptides were first enriched with a CapTrap column (0.5 mm × 2 mm, MICHROM Bioresources, Inc.) followed by elution into an integrated nanoscale analytical column (MAGIC C18AQ, 100 μm × 150 mm, 3 μm particle size, 200 Å pore size, MICHROM Bioresources, Inc.). Mobile phase A (2% ACN, 0.1% formic acid) and mobile phase B (98% ACN, 0.1% formic acid) were used to establish a 130 min gradient comprised of 5 min 5% B, then 25 min 5-15% B, followed by 55 min 15-40% B, then 15 min 40-80% B, maintained at 80% B for 10 min, then 5 min 80-5% B, finally maintained at 5% B for 15 min. The flow rate was ~300 nL/min. We conducted MS from 400 to 1800 amu, with 1 s time spans. For MS/MS analysis, each scan cycle consisted of one full-scan mass spectrum (with m/z ranging from 400 to 1800 and charge states from 2 to 5) followed by five MS/MS events. Threshold count was set to 30, and the exclusion window was 90 s. Mass tolerance was 50 mDa. Automatic Collision Energy and Automatic MS/MS Accumulation were selected.
The raw data collected by QSTAR ELITE were presented by Mascot Daemon (version 2.2.2) (Matrix Science, London, UK) to an in-house MASCOT server (version 2.2) (Matrix Science, London, UK) and Distiller (version 220.127.116.11) (Matrix Science, London, UK). Briefly, peak lists were generated by Distiller and searched against a target/decoy SwissProt human protein database (version 57.7; 20405 sequences) by MASCOT server. Spectra match criteria were as follows: Fixed modification was carbamidomethyl at Cys residue, whereas variable modifications were oxidation at Met residue, and phosphorylation at Ser, Thr or Tyr residues, additionally Arg10 and Lys8 were set as exclusive modifications, which is a useful setting in MASCOT that can be thought of as a choice of fixed modifications to speed up the search and reduce the significance thresholds, and taxonomy was set to “human”. Peptide and MS/MS tolerances were 50 ppm and 0.2 Da, respectively. The peptide charges were 2+, 3+, 4+, or 5+, allowing up to two missed cleavages. The significance threshold was p < 0.05.
The rov data obtained from database search were opened by Mascot Distiller for quantitation. For the quantitation analysis, we set Fraction, Correlation and Std. Err at 0.5, 0.9 and 0.2, respectively. The peptide ratios were calculated as weighted average ratios (ion intensity versus ratio) if several spectra for the same peptide were available. Protein ratios were also calculated as the weighted average ratios (ion intensity versus ratio). In a few situations, outliers were removed for more accurate quantitation. The median of all quantitation data from non-phosphopeptides was used to normalize the peptide ratios.
STRING network analysis
Identified proteins (accession number) were presented to STRING 8.2 database (http://string-db.org/)  to construct networks. Then networks were extracted and loaded into Cytoscape (version 2.6.3) (http://www.cytoscape.org/)  for compilation and visualization. Protein networks identified by STRING consisted of protein IDs (nodes) and protein-protein interactions (edges). Only interactions with a score of 0.400, which represents the default medium confidence level in STRING, were kept.
In order to predict potential kinases for quantified phosphosites, we used the NetworKIN algorithm which augments motif-based predictions with the network context of kinases and phosphoproteins and can assign a specific kinase to an identified in vivo phosphosite with a 2.5-fold higher accuracy than previous methods such as Scansite and NetphosK . After loading the data to NetworKIN 2.0 (http://networkin.info/version_2_0/newPrediction.php), we obtained sif data which could be rendered with Cytoscape (version 2.6.3). Manual network layout was then performed to display interactions between kinases and phosphosites, using color-coding to indicate fold increase or decrease.
This work was supported by the National Basic Research Program of China (2007CB914200) and the 111 Project of China (B06018). GML is supported by the National Institute of Health grants CA167181 and GM089684.
- Jiricny J: The multifaceted mismatch-repair system. Nat Rev Mol Cell Biol. 2006, 7: 335-346.View ArticlePubMedGoogle Scholar
- Li GM: Mechanisms and functions of DNA mismatch repair. Cell Res. 2008, 18: 85-98. 10.1038/cr.2007.115View ArticlePubMedGoogle Scholar
- Kolodner RD, Marsischky GT: Eukaryotic DNA mismatch repair. Curr Opin Genet Dev. 1999, 9: 89-96. 10.1016/S0959-437X(99)80013-6View ArticlePubMedGoogle Scholar
- Kunkel TA, Erie DA: DNA mismatch repair. Annu Rev Biochem. 2005, 74: 681-710. 10.1146/annurev.biochem.74.082803.133243View ArticlePubMedGoogle Scholar
- Modrich P, Lahue R: Mismatch repair in replication fidelity, genetic recombination, and cancer biology. Annu Rev Biochem. 1996, 65: 101-133. 10.1146/annurev.bi.65.070196.000533View ArticlePubMedGoogle Scholar
- Li GM: The role of mismatch repair in DNA damage-induced apoptosis. Oncol Res. 1999, 11: 393-400.PubMedGoogle Scholar
- Modrich P: Mechanisms in eukaryotic mismatch repair. J Biol Chem. 2006, 281: 30305-30309. 10.1074/jbc.R600022200PubMed CentralView ArticlePubMedGoogle Scholar
- Yoshioka K, Yoshioka Y, Hsieh P: ATR kinase activation mediated by MutSalpha and MutLalpha in response to cytotoxic O6-methylguanine adducts. Mol Cell. 2006, 22: 501-510. 10.1016/j.molcel.2006.04.023PubMed CentralView ArticlePubMedGoogle Scholar
- Duckett DR, Drummond JT, Murchie AI, Reardon JT, Sancar A, Lilley DM, Modrich P: Human MutSalpha recognizes damaged DNA base pairs containing O6-methylguanine, O4-methylthymine, or the cisplatin-d(GpG) adduct. Proc Natl Acad Sci U S A. 1996, 93: 6443-6447. 10.1073/pnas.93.13.6443PubMed CentralView ArticlePubMedGoogle Scholar
- Stojic L, Mojas N, Cejka P, Di Pietro M, Ferrari S, Marra G, Jiricny J: Mismatch repair-dependent G2 checkpoint induced by low doses of SN1 type methylating agents requires the ATR kinase. Genes Dev. 2004, 18: 1331-1344. 10.1101/gad.294404PubMed CentralView ArticlePubMedGoogle Scholar
- Duckett DR, Bronstein SM, Taya Y, Modrich P: hMutSalpha- and hMutLalpha-dependent phosphorylation of p53 in response to DNA methylator damage. Proc Natl Acad Sci U S A. 1999, 96: 12384-12388. 10.1073/pnas.96.22.12384PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Qin J: MSH2 and ATR form a signaling module and regulate two branches of the damage response to DNA methylation. Proc Natl Acad Sci U S A. 2003, 100: 15387-15392. 10.1073/pnas.2536810100PubMed CentralView ArticlePubMedGoogle Scholar
- Bignami M, O'Driscoll M, Aquilina G, Karran P: Unmasking a killer: DNA O(6)-methylguanine and the cytotoxicity of methylating agents. Mutat Res. 2000, 462: 71-82. 10.1016/S1383-5742(00)00016-8View ArticlePubMedGoogle Scholar
- Karran P, Bignami M: Self-destruction and tolerance in resistance of mammalian cells to alkylation damage. Nucleic Acids Res. 1992, 20: 2933-2940. 10.1093/nar/20.12.2933PubMed CentralView ArticlePubMedGoogle Scholar
- Sedgwick B, Lindahl T: Recent progress on the Ada response for inducible repair of DNA alkylation damage. Oncogene. 2002, 21: 8886-8894. 10.1038/sj.onc.1205998View ArticlePubMedGoogle Scholar
- Goldmacher VS, Cuzick RA Jr, Thilly WG: Isolation and partial characterization of human cell mutants differing in sensitivity to killing and mutation by methylnitrosourea and N-methyl-N'-nitro-N-nitrosoguanidine. J Biol Chem. 1986, 261: 12462-12471.PubMedGoogle Scholar
- Kat A, Thilly WG, Fang WH, Longley MJ, Li GM, Modrich P: An alkylation-tolerant, mutator human cell line is deficient in strand-specific mismatch repair. Proc Natl Acad Sci USA. 1993, 90: 6424-6428. 10.1073/pnas.90.14.6424PubMed CentralView ArticlePubMedGoogle Scholar
- Papadopoulos N, Nicolaides NC, Liu B, Parsons R, Lengauer C, Palombo F, D'Arrigo A, Markowitz S, Willson JK, Kinzler KW, et al: Mutations of GTBP in genetically unstable cells. Science. 1995, 268: 1915-1917. 10.1126/science.7604266View ArticlePubMedGoogle Scholar
- di Pietro M, Marra G, Cejka P, Stojic L, Menigatti M, Cattaruzza MS, Jiricny J: Mismatch repair-dependent transcriptome changes in human cells treated with the methylating agent N-methyl-n'-nitro-N-nitrosoguanidine. Cancer Res. 2003, 63: 8158-8166.PubMedGoogle Scholar
- Shen J, Chen W, Yin X, Yu Y: Proteomic analysis of different temporal expression patterns induced by N-methyl-N'-nitro-N-nitrosoguanidine treatment. J Proteome Res. 2008, 7: 2999-3009. 10.1021/pr800133qView ArticlePubMedGoogle Scholar
- Zelivianski S, Liang D, Chen M, Mirkin BL, Zhao RY: Suppressive effect of elongation factor 2 on apoptosis induced by HIV-1 viral protein R. Apoptosis. 2006, 11: 377-388. 10.1007/s10495-006-4030-9View ArticlePubMedGoogle Scholar
- Tu L, Liu Z, He X, He Y, Yang H, Jiang Q, Xie S, Xiao G, Li X, Yao K, et al: Over-expression of eukaryotic translation initiation factor 4 gamma 1 correlates with tumor progression and poor prognosis in nasopharyngeal carcinoma. Mol Cancer. 2010, 9: 78. 10.1186/1476-4598-9-78PubMed CentralView ArticlePubMedGoogle Scholar
- Prevot D, Darlix JL, Ohlmann T: Conducting the initiation of protein synthesis: the role of eIF4G. Biol Cell. 2003, 95: 141-156. 10.1016/S0248-4900(03)00031-5View ArticlePubMedGoogle Scholar
- Bertram J, Palfner K, Hiddemann W, Kneba M: Overexpression of ribosomal proteins L4 and L5 and the putative alternative elongation factor PTI-1 in the doxorubicin resistant human colon cancer cell line LoVoDxR. Eur J Cancer. 1998, 34: 731-736. 10.1016/S0959-8049(97)10081-8View ArticlePubMedGoogle Scholar
- Jiang H, Lin JJ, Tao J, Fisher PB: Suppression of human ribosomal protein L23A expression during cell growth inhibition by interferon-beta. Oncogene. 1997, 14: 473-480. 10.1038/sj.onc.1200858View ArticlePubMedGoogle Scholar
- Wolf-Yadlin A, Hautaniemi S, Lauffenburger DA, White FM: Multiple reaction monitoring for robust quantitative proteomic analysis of cellular signaling networks. Proc Natl Acad Sci U S A. 2007, 104: 5860-5865. 10.1073/pnas.0608638104PubMed CentralView ArticlePubMedGoogle Scholar
- Huang LJ, Durick K, Weiner JA, Chun J, Taylor SS: D-AKAP2, a novel protein kinase A anchoring protein with a putative RGS domain. Proc Natl Acad Sci U S A. 1997, 94: 11184-11189. 10.1073/pnas.94.21.11184PubMed CentralView ArticlePubMedGoogle Scholar
- McGowan CH, Russell P: The DNA damage response: sensing and signaling. Curr Opin Cell Biol. 2004, 16: 629-633. 10.1016/j.ceb.2004.09.005View ArticlePubMedGoogle Scholar
- Stokes MP, Rush J, Macneill J, Ren JM, Sprott K, Nardone J, Yang V, Beausoleil SA, Gygi SP, Livingstone M, et al: Profiling of UV-induced ATM/ATR signaling pathways. Proc Natl Acad Sci U S A. 2007, 104: 19855-19860. 10.1073/pnas.0707579104PubMed CentralView ArticlePubMedGoogle Scholar
- Cole CN, Hammell CM: Nucleocytoplasmic transport: driving and directing transport. Curr Biol. 1998, 8: R368-R372. 10.1016/S0960-9822(98)70239-8View ArticlePubMedGoogle Scholar
- Sorokin AV, Kim ER, Ovchinnikov LP: Nucleocytoplasmic transport of proteins. Biochemistry (Mosc). 2007, 72: 1439-1457. 10.1134/S0006297907130032. 10.1134/S0006297907130032View ArticleGoogle Scholar
- Quimby BB, Dasso M: The small GTPase Ran: interpreting the signs. Curr Opin Cell Biol. 2003, 15: 338-344. 10.1016/S0955-0674(03)00046-2View ArticlePubMedGoogle Scholar
- Mahajan R, Delphin C, Guan T, Gerace L, Melchior F: A small ubiquitin-related polypeptide involved in targeting RanGAP1 to nuclear pore complex protein RanBP2. Cell. 1997, 88: 97-107. 10.1016/S0092-8674(00)81862-0View ArticlePubMedGoogle Scholar
- Beausoleil SA, Jedrychowski M, Schwartz D, Elias JE, Villen J, Li J, Cohn MA, Cantley LC, Gygi SP: Large-scale characterization of HeLa cell nuclear phosphoproteins. Proc Natl Acad Sci U S A. 2004, 101: 12130-12135. 10.1073/pnas.0404720101PubMed CentralView ArticlePubMedGoogle Scholar
- Jensen LJ, Kuhn M, Stark M, Chaffron S, Creevey C, Muller J, Doerks T, Julien P, Roth A, Simonovic M, et al: STRING 8-a global view on proteins and their functional interactions in 630 organisms. Nucleic Acids Res. 2009, 37: D412-D416. 10.1093/nar/gkn760PubMed CentralView ArticlePubMedGoogle Scholar
- Harper JW, Elledge SJ: The DNA damage response: ten years after. Mol Cell. 2007, 28: 739-745. 10.1016/j.molcel.2007.11.015View ArticlePubMedGoogle Scholar
- Schroering AG, Williams KJ: Rapid induction of chromatin-associated DNA mismatch repair proteins after MNNG treatment. DNA Repair (Amst). 2008, 7: 951-969. 10.1016/j.dnarep.2008.03.023. 10.1016/j.dnarep.2008.03.023View ArticleGoogle Scholar
- Matsuoka S, Ballif BA, Smogorzewska A, McDonald ER 3rd, Hurov KE, Luo J, Bakalarski CE, Zhao Z, Solimini N, Lerenthal Y, et al: ATM and ATR substrate analysis reveals extensive protein networks responsive to DNA damage. Science. 2007, 316: 1160-1166. 10.1126/science.1140321View ArticlePubMedGoogle Scholar
- Iijima K, Muranaka C, Kobayashi J, Sakamoto S, Komatsu K, Matsuura S, Kubota N, Tauchi H: NBS1 regulates a novel apoptotic pathway through Bax activation. DNA Repair (Amst). 2008, 7: 1705-1716. 10.1016/j.dnarep.2008.06.013. 10.1016/j.dnarep.2008.06.013View ArticleGoogle Scholar
- Yu X, Fu S, Lai M, Baer R, Chen J: BRCA1 ubiquitinates its phosphorylation-dependent binding partner CtIP. Genes Dev. 2006, 20: 1721-1726. 10.1101/gad.1431006PubMed CentralView ArticlePubMedGoogle Scholar
- Chen PL, Liu F, Cai S, Lin X, Li A, Chen Y, Gu B, Lee EY, Lee WH: Inactivation of CtIP leads to early embryonic lethality mediated by G1 restraint and to tumorigenesis by haploid insufficiency. Mol Cell Biol. 2005, 25: 3535-3542. 10.1128/MCB.25.9.3535-3542.2005PubMed CentralView ArticlePubMedGoogle Scholar
- Hopfner KP, Craig L, Moncalian G, Zinkel RA, Usui T, Owen BA, Karcher A, Henderson B, Bodmer JL, McMurray CT, et al: The Rad50 zinc-hook is a structure joining Mre11 complexes in DNA recombination and repair. Nature. 2002, 418: 562-566. 10.1038/nature00922View ArticlePubMedGoogle Scholar
- Lee JH, Paull TT: ATM activation by DNA double-strand breaks through the Mre11-Rad50-Nbs1 complex. Science. 2005, 308: 551-554. 10.1126/science.1108297View ArticlePubMedGoogle Scholar
- Gichner T, Veleminsky J: Genetic effects of N-methyl-N'-nitro-N-nitrosoguanidine and its homologs. Mutat Res. 1982, 99: 129-242. 10.1016/0165-1110(82)90057-4View ArticlePubMedGoogle Scholar
- Bendall SC, Hughes C, Stewart MH, Doble B, Bhatia M, Lajoie GA: Prevention of amino acid conversion in SILAC experiments with embryonic stem cells. Mol Cell Proteomics. 2008, 7: 1587-1597. 10.1074/mcp.M800113-MCP200PubMed CentralView ArticlePubMedGoogle Scholar
- Albuquerque CP, Smolka MB, Payne SH, Bafna V, Eng J, Zhou H: A multidimensional chromatography technology for in-depth phosphoproteome analysis. Mol Cell Proteomics. 2008, 7: 1389-1396. 10.1074/mcp.M700468-MCP200PubMed CentralView ArticlePubMedGoogle Scholar
- Chen X, Wu D, Zhao Y, Wong BH, Guo L: Increasing phosphoproteome coverage and identification of phosphorylation motifs through combination of different HPLC fractionation methods. J Chromatogr B Analyt Technol Biomed Life Sci. 2011, 879: 25-34. 10.1016/j.jchromb.2010.11.004View ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303PubMed CentralView ArticlePubMedGoogle Scholar
- Linding R, Jensen LJ, Ostheimer GJ, van Vugt MA, Jorgensen C, Miron IM, Diella F, Colwill K, Taylor L, Elder K, et al: Systematic discovery of in vivo phosphorylation networks. Cell. 2007, 129: 1415-1426. 10.1016/j.cell.2007.05.052PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.