Time-resolved Analysis of Proteome Dynamics by Tandem Mass Tags and Stable Isotope Labeling in Cell Culture (TMT-SILAC) Hyperplexing*

Recent advances in mass spectrometry have enabled system-wide analyses of protein turnover. By globally quantifying the kinetics of protein clearance and synthesis, these methodologies can provide important insights into the regulation of the proteome under varying cellular and environmental conditions. To facilitate such analyses, we have employed a methodology that combines metabolic isotopic labeling (Stable Isotope Labeling in Cell Culture - SILAC) with isobaric tagging (Tandem Mass Tags - TMT) for analysis of multiplexed samples. The fractional labeling of multiple time-points can be measured in a single mass spectrometry run, providing temporally resolved measurements of protein turnover kinetics. To demonstrate the feasibility of the approach, we simultaneously measured the kinetics of protein clearance and accumulation for more than 3000 proteins in dividing and quiescent human fibroblasts and verified the accuracy of the measurements by comparison to established non-multiplexed approaches. The results indicate that upon reaching quiescence, fibroblasts compensate for lack of cellular growth by globally downregulating protein synthesis and upregulating protein degradation. The described methodology significantly reduces the cost and complexity of temporally-resolved dynamic proteomic experiments and improves the precision of proteome-wide turnover data.

Recent advances in mass spectrometry have enabled system-wide analyses of protein turnover. By globally quantifying the kinetics of protein clearance and synthesis, these methodologies can provide important insights into the regulation of the proteome under varying cellular and environmental conditions. To facilitate such analyses, we have employed a methodology that combines metabolic isotopic labeling (Stable Isotope Labeling in Cell Culture -SILAC) with isobaric tagging (Tandem Mass Tags -TMT) for analysis of multiplexed samples. The fractional labeling of multiple time-points can be measured in a single mass spectrometry run, providing temporally resolved measurements of protein turnover kinetics. To demonstrate the feasibility of the approach, we simultaneously measured the kinetics of protein clearance and accumulation for more than 3000 proteins in dividing and quiescent human fibroblasts and verified the accuracy of the measurements by comparison to established nonmultiplexed approaches. The results indicate that upon reaching quiescence, fibroblasts compensate for lack of cellular growth by globally downregulating protein synthesis and upregulating protein degradation. The described methodology significantly reduces the cost and complexity of temporally-resolved dynamic proteomic experiments and improves the precision of proteomewide turnover data. Analysis of protein expression is central to our understanding of cellular physiology and regulation. Recent advances in the field of proteomics have provided powerful experimental tools for global analyses of protein expression. Tandem mass spectrometry (LC-MS/MS) 1 can now readily capture information on composition, expression levels and modifications of very complex protein mixtures (1). Though immensely powerful, these techniques typically provide temporally static snapshots of the proteome. Within a cell, proteins are in a state of continual flux, and kinetics of synthesis and degradation (turnover) are intimately linked to the function and regulation of proteins (2,3). As such, there is great interest in the development of proteomic methodologies capable of analyzing protein turnover dynamics on a global scale (4).
Proteome-wide analyses of turnover have been greatly propelled by the development of stable isotope labeling in cell culture (SILAC)-a set of standardized protocols, reagents and analysis tools for quantifying the incorporation of 15 N/ 13 C labeled amino acids into proteins by tandem mass spectrometry (LC-MS/MS) (5,6). Dynamic (i.e. time-resolved or pulsed) SILAC experiments (Fig. 1) can effectively measure the labeling kinetics of proteins and provide rates of protein degradation and synthesis on a global scale. Such analyses have successfully obtained proteome-wide surveys of protein dynamics for a number of in vitro and in vivo model systems and have provided important insights into the regulation of the proteome under varying environmental and cellular conditions (2,(7)(8)(9)(10).
A practical limitation of SILAC-based dynamic proteomic analyses is that each labeling time point produces an experimental sample that must be analyzed in a distinct LC-MS/MS run. Thus, conducting kinetic analyses composed of large numbers of time-points is expensive and time consuming. Furthermore, SILAC-based dynamic proteomic analyses of-ten suffer from the problem of "missing values": proteins that can be identified and quantified in some time-points but not others. To circumvent these limitations, we sought to develop a more efficient procedure for obtaining multi-time-point kinetic data for proteome-wide analyses of protein turnover. The methodology combines and builds upon three previously reported techniques: integration of SILAC and tandem mass tags (TMT) to produce highly multiplexed samples ("hyperplexing" (11,12)), synchronous precursor selection (SPS) to minimize co-isolation of interfering ions and enhance accuracy of MS n quantitation of TMT reporters (13,14) and incorporation of isobaric and metabolic labeling for analysis of protein dynamics (15,16). Using this combined approach, we analyzed the global impact of cellular quiescence on the kinetics of protein degradation and synthesis in human dermal fibroblasts.

EXPERIMENTAL PROCEDURES
Experimental Design and Statistical Rationale-The theoretical rationale for the assay workflow ( Fig. 2A) is discussed in the Results section. Briefly, two replicate experiments were conducted. In the first experiment, 10 cultures of human dermal fibroblasts were grown to confluency and labeled for varying times in SILAC media. The timepoints were then labeled with different TMT tags, multiplexed and the fractional labeling of the peptides were analyzed at the MS1, MS2 and MS3 levels using the procedures described below. In a second biologically replicate experiment, the quiescence experiment was repeated with four time-points and the samples were multiplexed with similarly labeled dividing cells for comparison. The determination of the rate constants for fractional labeling was conducted by least squares regression analysis of the time-resolved data fitted to equations derived from the kinetic model described below and the goodness of fit was assessed by R 2 values. The spread of rate constants for peptides encompassing the same protein was assessed by measuring the coefficient of variation, as the mean rates for peptides for different proteins had variable magnitudes. The correlation of measured rates at the protein level using different procedures was assessed by conducting pairwise comparisons. For this analysis, Spearman rank correlation was used, as the measured degradation rates were not normally distributed.
Cell Culture and Stable Isotope Labeling-Human dermal fibroblasts (HCA2-hTert) (17,18) were maintained in Eagle's Minimum Essential Medium (ATCC) supplemented with 15% fetal bovine serum (Invitrogen, Carlsbad, CA), 100 U/ml penicillin, 100 U/ml streptomycin at 37°C with 5% CO 2 . The media utilized for isotopic labeling was Eagle's minimum essential medium (ATCC) supplemented with 15% dialyzed fetal bovine serum (Thermo Scientific, Boston, MA), 100 U/ml penicillin, and 100 U/ml streptomycin. Cells were gradually adapted to the labeling media and were then plated at a density of 500,000 cells per 10 cm plate. For dividing cells, 1 day after plating, the cultures were switched to MEM labeling media for SILAC (Thermo Scientific) supplemented with L-arginine:HCl ( 13 C6, 99%) and L-lysine: 2HCl ( 13 C6, 99%; Cambridge Isotope Laboratories, Tewksbury, MA) at concentrations of 0.1264 g/l and 0.087 g/l and 15% dialyzed fetal bovine serum (Thermo Scientific) and were subsequently collected after 0, 24, 48, 72 h of labeling. Cell numbers were determined by a hemocytometer. The number of cells increased in an approximate exponential fashion during this period. Least-square fitting of a plot of cell numbers versus time to an exponential equation was used to determine the growth rate (k div ) of the cells. For quiescent cells, the cultures were grown for 8 days to achieve a state of quiescence by contact inhibition (see (19) for additional details). The confluent quiescent cultures were switched to the labeling media and cells were collected after 0, 6,12,24,36,48,72,96,144,192, 336 h of labeling. All cells were washed with PBS and frozen as cell pellets prior to further analysis.
Mass Spectrometry Sample Preparation-Cell Lysis-Lysis buffer was composed of 8 M urea, 75 mM NaCl, 50 mM Tris, pH 8.5. 50 l of lysis buffer was used for one million cells. Cell lysis was conducted by adding lysis buffer to pelleted cells, followed by vortexing and sonication using a QSonica sonicator. Sonication cycles consisted of 5 ϫ 10 s sonications with 1 min incubations on ice between each cycle. Samples were then centrifuged for 10 min at 15,000 ϫ g, and the supernatant was collected. Protein concentration was determined by BCA assays.
Trypsinization and TMT Labeling-One hundred micrograms of each extract was removed, and brought up to 100 l in 100 mM ammonium bicarbonate. This dilution lowered the urea concentration to less than 1.5 M such that trypsinization efficiency was not inhibited. 2 g of trypsin (Thermo Scientific) was added to each extract and the samples were allowed to incubate overnight at 37°C. Disulfide bonds were reduced with the addition of Bond-Breaker TCEP Solution (Thermo Scientific) to a final concentration of 5 mM with incubation at 55°C for 1 h. Reduced cysteines were alkylated by the addition of iodoacetamide at a final concentration of 10 mM by incubating in the dark at room temperature for 30 min. Subsequently, an additional 2 g of trypsin was added to the samples and incubated for 3 h at 37°C. Formic acid was added to a final concentration of 1% to stop trypsin activity. Samples were then desalted with 100 l C18 Tips (Thermo Scientific) to remove primary amines that would interfere with TMT labeling. After drying down, the samples were re-suspended in 100 l of 100 mM triethylammonium bicarbonate (TEAB). TMT tenplex reagents (0.8 mg) (Thermo Scientific) were removed from Ϫ20°C and allowed to reach room temperature before being re-suspended in 40 l of acetonitrile. Ten microliters of each reagent was added to the corresponding sample. For analysis of degradation kinetics in quiescent cells, the ten samples derived from different SILAC labeling time-points were tagged with TMT-126, TMT-127N, TMT-127C, TMT-128N, TMT-128C, TMT-129N, TMT-129C, TMT-130N, TMT-130C, and TMT-131 in order of increasing time. For comparison of degradation and synthesis kinetics between dividing and quiescent cells, the dividing cells SILAC labeled for 0, 24, 48 and 72 h were tagged with TMT-126, TMT-127N, TMT-127C, and TMT-128N, respectively, and the quiescent cells were labeled SILAC labeled for 0, 48, 96 and 144 h were tagged by TMT-129C, TMT-130N, TMT-130C, and TMT-131, respectively. After the addition of TMT tags, the samples were vortexed and incubated at room temperature for 1 h before being quenched with 8 l of 5% hydroxylamine for 15 min. 10 g of each TMT-labeled samples were combined according to the experimental design described in Results and dried down by SpeedVac (Labconco, Kansas City, MO). The remaining labeled peptides were left un-mixed for quantitation using the precursor MS1 SILAC peaks. Un-mixed labeled peptides were desalted with 100 l C18 tips, and dried down by SpeedVac.
High pH Fractionation-Digested samples were reconstituted in 300 l of 0.1% TFA and were fractionated using the Pierce High pH Reversed-Phase Peptide Fractionation Kit (Thermo Scientific). Briefly, the columns were prepared by adding 300 l of acetonitrile and centrifuging at 5000 ϫ g for 2 min (twice), then washed with 300 l 0.1% TFA (twice). Samples were then loaded onto the column and centrifuged at 3000 ϫ g for 2 min, then washed once with HPLC grade water (JT Baker). Peptides were eluted in 8 fractions, with 10%, 12.5%, 15%, 17.5%, 20%, 22.5%, 25%, 50% acetonitrile in 0.1% triethylamine. Fractions were dried down by SpeedVac, resuspended in 50 l of 0.1% TFA, and analyzed as described below.

LC-MS/MS Analysis-LC-MS/MS for MS2 TMT and SILAC Analysis Using Q Exactive Plus
Mass Spectrometer-Six microliters of each high pH fraction was analyzed on a Q Exactive Plus mass spectrometer (Thermo Scientific) coupled with an Easy nLC-1000 pump (Thermo Scientific). Columns were hand-pulled using 100 m fused silica, which was packed with 30 cm of 1.8 m, 120 Angstrom C18 beads (Sepax). Mobile phase A was 0.1% formic acid in water, and mobile phase B was 0.1% formic acid in acetonitrile. Peptides were separated using a gradient of 8 -30% B over 145 min, 30 -50% B over 10 min, and 50 -70% B over 4 min, holding at 70% B for 5 min. The gradient returned to 0% B over 4 min, and held there for 10 min to prepare for the next injection. The flow rate throughout the run was held at 300 nL/minute. A data-dependent top 10 MS2 method was used for peptide detection and fragmentation. For each cycle, one full scan from m/z 400 -1700 was acquired in the Orbitrap at a resolution of 70,000 at m/z ϭ 200, with an AGC target of 3e6, and a maximum injection time of 50 ms. After each full scan, the top 10 most intense ions were selected for fragmentation in the HCD cell, followed by detection in the Orbitrap at a resolution of 35,000 at m/z ϭ 200. The AGC target was set to 1e5, with a maximum injection time of 150 ms, an isolation width of 0.5 Da, and a normalized collision energy of 35. A dynamic exclusion window of 25 s was used, while excluding unassigned, ϩ1, and greater than ϩ5 charge states. The polysiloxane 445.12003 m/z was used for lockmass. For SILAC LC-MS/MS analysis on the TMT-labeled, unmixed samples, all settings remained the same except for the isolation width, which was adjusted up to 1.0 Da, with an offset of 0.3 Da.

LC-MS/MS for SPS-MS3 Analysis Using Orbitrap Fusion Tribrid Mass
Spectrometer-Six microliters of each fraction was analyzed on a Fusion mass spectrometer (Thermo Scientific). Chromatographic separation of derived peptides were accomplished on a SCIEX ekspert TM nanoLC 425 system (SCIEX, Redwood City, CA). Mobile phase A and B were 0.1% formic acid in 2% acetonitrile and 0.1% formic acid in 88% acetonitrile, respectively. Peptides were loaded onto a reversed-phase trap (300 m ID ϫ1 cm, packed with Zorbax 3 m C18 material), with 1% mobile phase B at a flow rate of 10 l/min, and the trap was washed for 3 min. A series of nanoflow gradients (flow rate at 250 nL/min) was used to back-flush the trapped samples onto the nano-LC column (75-m ID ϫ 100 cm, packed with Pepmap 3 m C18 material) for separation. The nano-LC column was heated at 52°C to improve both chromatographic resolution and reproducibility. A 3 h gradient was used to achieve sufficient peptide separation. The optimized gradient profile was as follows: 4-10% B over 3-8 min; 10 -32% B over 8 -125 min; 32-44% B over 125-130 min; 44 -60% B over 130 -135 min; 60 -97% B in 135-136 min, and isocratic at 97% B over 136 -153 min. A data-dependent multi-notch synchronous precursor selection (SPS) method was employed (13). For each cycle, one MS1 scan was done in the Orbitrap with a mass range of 400 -1500 m/z at a resolution of 120,000 at 200 m/z. The AGC was set to 2e5 with a maximum ion injection time of 50 ms. Dynamic exclusion was used with a 10-ppm tolerance, and a 40 s exclusion duration. MS2 scans were done in the ion trap with CID, fragmenting the most intense precursors in the Top Speed mode, with an isolation width of 0.5 Da. The "Ion Trap Scan Rate" was set to Rapid, with an AGC target of 1e4, and a maximum injection time of 35 ms. Only ions with a charge state between 2 and 5 were selected for fragmentation. MS3 scans were done by isolating the 10 most intense fragment ions from CID fragmentation over a 400 -2000 m/z range, excluding 40 Da below and 15 Da above the precursor ion. These fragment ions were sent to the HCD cell using the SPS waveform with an isolation width of 2.0 Da and an AGC of 1e5, and subsequently fragmented with an HDC collision energy of 65%, then detected in the Orbitrap at a resolution of 60,000 at 200 m/z. The maximum injection time was set to 120 ms.
Data Analysis-Peptide Identification-All LC-MS/MS runs were analyzed using Mascot version 2.3.2 (Matrix Science) within Proteome Discoverer 2.1 (Thermo Scientific) against the SwissProt human database (2015/12/ 29, 42,085 entries). The eight fractions that were run on each instrument were searched together to create one report. A 10 ppm MS1 error tolerance was used, with a 25 mmu error used for MS2 product ions on the Q Exactive Plus, and 0.8 Da on the Fusion, because the MS2s were done in the ion trap. Trypsin was set as the enzyme, allowing for 1 missed cleavage. Methionine oxidation (15.994915 Da) was set as a variable modification, and carbamidomethyl (57.021464 Da) on cysteines was set as a static modification for all searches. Because we could not search the heavy lysine as a dynamic modification with the TMT6 plex on lysine as a static modification, we had to conduct separate searches to get data for both the unlabeled and labeled SILAC peptides. When searching for the unlabeled SILAC peak, TMT6 plex modifications (229.162932 Da) on both lysine and peptide N termini were used. When searching for the labeled SILAC peak, a new modification was created that added the heavy lysine mass with the TMT6 plex mass, for a total mass of 235.183061 Da. This was set as a static modification, along with heavy labeled arginine (6.020129 Da), TMT6 plex on peptide N termini, and carbamidomethyl on cysteine. Percolator was used as the false discovery rate calculator, and all peptides were filtered at the "Strict" Target FDR level of 0.01.
SILAC Quantitation-SILAC peptides and protein quantification was performed with Proteome Discoverer 1.4 using the Precursor Ions Quantifier node. For each peptide, the H/L ratio was determined by a regression model fitted to all isotopic peaks within all scans that the peptide eluted in. H/L ratio for each protein was determined as the median of all the peptides assigned to the protein. Analysis of the isotopic distribution of incompletely cleaved peptides that contained two lysines or arginines indicated that Ͼ98% of the internal amino acid pool was labeled at all time-points (8). For a given protein, fractional labeling was calculated as the H/(HϩL) ratio based on the H/L ratio from Proteome Discoverer outputs. Peptide fractional labeling measurements at all time-points were combined to obtain an aggregated plot for the kinetics of labeling. The aggregated plots were fitted to a single exponential function by least-square fitting to obtain first order degradation rate constant (k deg ) for each protein.
TMT Quantitation-For peak integration of the TMT product ions, an error tolerance of 15 ppm with the integration method set to "Centroid with the Smallest Delta Mass" within the "Reporter Ions Quantifier" node was used. The MS Order was set to MS2 for the Q Exactive Plus data, and MS3 when searching the Fusion data. Reporter abundances were based on Signal/Noise ratio (S/N) if all spectrum files had a S/N value. Otherwise, reporter abundances were based on signal intensities. Peptides with an isolation interference of greater than 30% were not used for quantitation. Fractional labeling of PSMs at all time-points were combined to obtain an aggregated plot for the kinetics of labeling. The aggregated plots were fitted to a single exponential function by least-square fitting to obtain first order degradation rate constant (k deg ) for each protein. When aggregating PSM data, the kinetics of labeling for individual PSMs were initially fit to a single exponential function by least-square fitting. The data were then filtered to exclude PSMs that had R 2 Ͻ0.8 prior to aggregation and protein-level determination of k deg . When fitting exponential curves for decay of signal from TMT reporters derived from SILAC-unlabeled peptides, TMT measurements were normalized with respect to the t0 reporter (TMT-126) and a constant fractional baseline at infinite time was incorporated in the fitting equation. The equation used for fitting the curves was therefore: intensity ϭ baseline ϩ ͑1 Ϫ baseline͒ ϫ e Ϫkdegϫtime . When fitting exponential curves for increase of signal from TMT reporters derived from SILAC-labeled peptides, TMT measure-ments were normalized with respect to the t0 reporter (TMT-126) and the amplitude of the curve was treated as a fitted parameter in the regression analysis. The equation used for fitting the curves was therefore: intensity ϭ 1 ϩ amplitude ϫ ͑1 Ϫ e Ϫkdegϫtime ). For measurement of k deg in proliferating cells, the rate constant for cell division (0.40 d Ϫ1 ) was subtracted from the measured rate constants (see below). This subtraction should theoretically yield a positive rate constant. However, for a fraction of proteins, a negative calculated k deg value was measured, likely because of experimental error. These proteins were not considered in subsequent analyses.
Kinetic Model-The kinetic model applied in this study is based on the following assumptions: 1 Protein synthesis is a zero-order process with respect to protein concentration.
2 Protein degradation occurs at a constant fractional rate that is uniform for the entire protein pool. Thus, protein degradation can be modeled as a first order process with respect to protein concentration.
3 The total protein concentration of each cell does not change during the experimental time-course and the system is at steady-state.
Based on these assumptions, we can devise the following rate equations for clearance of unlabeled proteins and accumulation of labeled proteins: Where k syn is the zero-order rate constant for protein synthesis, k deg is the first order rate constant for protein degradation and k div is the first order rate constant for cell division. We

RESULTS
Assay Workflow-Our approach takes advantage of the multiplexing capabilities of isobaric tandem mass tags (TMT) (20 -22) to measure fractional SILAC labeling of multiple timepoints in a single LC-MS/MS run (Fig. 1B). In this procedure, isotopically labeled amino acids are introduced to cultured cells. Over time, as cellular proteins are degraded and replaced by newly synthesized proteins, the fractional labeling of the protein pool increases. At each time point, protein extracts are collected, digested with trypsin and labeled with distinct TMT tags. Tagged peptides from different time-points are subsequently combined into a single multiplexed sample and analyzed by LC-MS/MS with higher-energy collisional dissociation (HCD). Isobaric TMT tags have identical intact masses but produce reporter ions of varying masses following HCD. Hence, tryptic peptides of a given sequence obtained from different time-points will have identical masses in MS1 spectra but will produce reporter ions of varying masses in MS2 spectra. As protein turnover leads to the gradual disappearance of SILAC-unlabeled and appearance of SILAC-labeled MS1 peaks, the kinetics of this process are revealed by the relative intensities of the TMT reporter ions in the MS2 spectra. Specifically, the relative intensities of TMT reporter ions generated by SILAC-unlabeled precursor peptides monitors their clearance and the relative intensities of TMT reporter ions generated by SILAC-labeled precursor peptides monitors their accumulation.
As a proof of principle, we used the above approach to globally measure the kinetics of protein turnover in primary human dermal fibroblasts (HCA2-hTert) (17,18) (Fig. 2A). Cells were grown to confluency, resulting in the arrest of cell division by contact inhibition. After reaching quiescence, the cultures were switched to a media containing 13 C labeled lysines and arginines. The fact that cells were labeled in a non-dividing state ensured that dilution of cellular proteins by cell division did not contribute to the kinetics of fractional labeling, and protein turnover was the sole determinant of labeling kinetics (2,23). Cells were collected at 10 different time-points ranging from 0 to 336 h. For each time point, cell extracts were digested with trypsin, covalently tagged with ten different TMT tags before being combined in a single multiplexed sample. Following offline fractionation, the multiplexed sample was analyzed by LC-MS/MS as described under Experimental Procedures. Additionally, in order to com-pare the results with a standard dynamic SILAC approach, TMT-tagged peptides at five selected time-points were individually analyzed at the MS1 level without multiplexing.
Representative TMT-SILAC Hyperplexed Data-To illustrate a typical datapoint, Figs. 2B-2F show the labeling kinetics of a peptide belonging to neuroblast differentiation-associated protein (AHNAK) (24). Protein turnover was quantified by five different parameters: 1 Fractional labeling of intact precursor peptides in MS1 spectra obtained from individual non-multiplexed time-points (Dynamic SILAC; Fig. 2B).
3 Relative intensities of TMT reporter ions obtained from MS2 spectra of SILAC-labeled precursor peptides (Fig. 2D). 4 Relative intensities of TMT reporter ions obtained in MS3 spectra of SILAC-unlabeled precursor peptides by isolation of multiple fragment ions in the MS2 spectra using isolation waveform with multiple notches (synchronous precursor selection; SPS-MS3 (13, 14)) (Fig. 2E). The SPS-MS3 method has been reported to diminish TMT ratio distortion by partially alleviating the co-isolation and fragmentation of interfering ions (13,14). 5 Relative intensities of TMT reporter ions obtained by SPS-MS3 spectra of SILAC-labeled precursor peptides. (Fig. 2F).
The kinetic model used to interpret labeling kinetics is presented under Experimental Procedures. The salient assumptions of the model are as follows: (1) the system is at steadystate and population-averaged cellular protein concentrations are not changing during the course of the experiment, (2) protein degradation is a first order process with respect to protein concentration, (3) protein synthesis is a zero order process with respect to protein concentration, and (4) given Resulting curves for either the decay of unlabeled peptides or appearance of labeled peptides are fitted to exponential equations to obtain the turnover rate constants or half-lives of proteins. B, In TMT-SILAC hyperplexing, experimental samples from different SILAC labeling time-points are reacted with different TMT tags and combined in a single multiplexed sample prior to LC-MS/MS. At different time-points, a peptide generated from the same protein will have identical (overlapped) masses in MS1 spectra but will generate reporter ions with distinct masses in MS2 spectra. Thus, fragmentation of SILAC-unlabeled peaks will provide the kinetics of decay for unlabeled peptides and fragmentation of SILAC-labeled peaks will provide the kinetics of appearance for labeled peptides.

FIG. 2. TMT-SILAC hyperplexing experimental design and sample data.
A, Experimental design. HCA2-hTert cells were grown to confluency and SILAC labeling was conducted over the course of 336 h as cells were in a contact-inhibited quiescent state. Extracts were collected from 10 time-points and, following trypsinization, each was labeled with a different TMT tag. The TMT-tagged peptides were combined in a single multiplexed sample and analyzed by LC-MS/MS as described under Experimental Procedures. Additionally, a number of TMT-tagged time-points (boxed) were individually analyzed by LC-MS/MS without multiplexing using a standard dynamic SILAC approach. B, The fractional labeling of a tryptic peptide ion (GEYDVTVPK, z ϭ 2) belonging to neuroblast differentiation-associated protein analyzed at the MS1 level by dynamic SILAC. The green and orange peaks indicate the relative intensities of unlabeled and labeled peptides, respectively. The plot indicates the decay curve for the fractional labeling of the unlabeled peaks for GEYDVTVPK (green circles) and all other peptides mapped to the protein (transparent blue circles). The green line indicates a fit to a single exponential equation. The fractional labeling of labeled peaks (calculated as one minus the fractional labeling of unlabeled peaks) is indicted by a dotted orange line. C, The kinetics for clearance of GEYDVTVPK analyzed at the MS2 level by TMT-SILAC hyperplexing. The MS2 spectrum indicates the relative intensities of TMT reporter ions for the SILAC-unlabeled peptide ion. The plot indicates intensities of peaks relative to the t0 sample (TMT-126) for GEYDVTVPK (green circles) that the cells are in a quiescent state, protein dilution by cell division is not impacting the fractional labeling of proteins. Based on this model, both protein clearance and accumulation should follow single exponential kinetics with the apparent rate constant equaling the rate constant for protein degradation (k deg ). In support of this model, Fig. 2 indicates that peptide labeling measured by all five methods fit well to single exponential equations and have similar apparent first order rate constants (0.013-0.014 h Ϫ1 ).
It is interesting to note that the kinetic curves measured by TMT-SILAC hyperplexing methods have non-zero baseline intensities (Fig. 2C-2F, dashed lines). In MS2/MS3 analyses of SILAC-unlabeled peptides, the kinetic curve would ordinarily be expected to level off at zero intensity as the protein population becomes fully labeled at progressively longer timepoints. Conversely, SILAC-labeled peptides would be expected to have zero intensities at time zero. However, at the MS2 level, both curves have non-zero baselines that have a magnitude of ϳ10 -15% of the total amplitude of the curve. This suggests that at the MS2 level, a fraction of TMT reporter ions are not generated by the targeted peptide precursors but are instead produced by co-isolation of interfering ions. This interference phenomenon has been described by a number of previous publications (13,14,(25)(26)(27)(28)(29)(30). We observed that for most peptides, the relative contribution of interfering precursor ions to total TMT signal is roughly constant at varying time-points. For example, for MS2 TMT analyses of SILAC-labeled precursors peptides mapped to neuroblast differentiation-associated protein (Fig. 2C), ϳ12% of the total intensity of the reporter ions remained at long time-points where Ͼ99% of the protein was expected to have turned over (Fig. 2B). Nonetheless, the fractional labeling curve was well fit by a single exponential equation with a baseline of 0.12. Thus, although co-isolation of interfering precursors are likely distorting the measurements of fractional labeling, their influence on the measured rate constant is minimized by the fact that our regression analysis accounts for the presence of a constant baseline intensity that is independently determined as a fitted parameter.
To further mitigate the contribution of interfering species, we used synchronous precursor selection (SPS) to co-isolate and fragment multiple MS2 fragment ions generated by the targeted precursor as described by McAlister et al. (13). The reporter ions generated in the MS3 spectra were quantified as described above for MS2 spectra (Fig. 2E, 2F). The reduction in co-isolation of interfering ions was evident by decreases in baseline intensities (compare dotted baselines between Figs. 2C, 2D and Figs. 2E, 2F).

Global Validation of Precision and Accuracy of TMT-SILAC
Hyperplexing-Coverage statistics for all proteome-wide analyses are summarized in supplemental Table S1. By combining the four TMT-SILAC hyperplexing analyses, we were able to measure the degradation kinetics of peptides belonging to 3009 unique protein groups in one or more analysis (supplemental Table S2). To obtain a protein-level data set, we aggregated PSMs whose labeling kinetics were able to be fit with a first order rate equation having R 2 values greater than 0.8, and then limited our analysis to 1307 protein groups with 2 or more unique quantified peptides that passed this criterion in at least one of the analyses (supplemental Table S3). We next compared the TMT-hyperplexing data to dynamic SILAC experiments based on quantitation of the MS1 spectra (supplemental Tables S4, S5). For these comparisons, we limited our analyses to 229 proteins that were quantified by all five approaches. First, we analyzed the correspondence of degradation rates of different peptides mapped to the same protein (Fig.  3). For most proteins, it would be expected that all peptide ions matched to a given protein should have the same apparent k deg . Indeed, we observe that the range of values for coefficients of variation (CV) for filtered peptide k deg s within a protein is very low in comparison to the CV for all analyzed peptide k deg s within the proteome. The range of CVs for intra-protein peptides measured by all approaches had similar magnitudes. However, the mean CV for dynamic SILAC analyses were somewhat lower than the TMT-SILAC approaches, suggesting that it has somewhat better experimental precision. The comparison also indicates that SPS-MS3 improves the precision of the TMThyperplexing approach.
Next, we analyzed the correlation between protein k deg measurements obtained from the different analyses in pairwise comparisons. (Fig. 4). The data indicate that there is significant correlation between the measurements obtained from the five methods, validating the general accuracy of k deg measurements. However, there is considerable experimental noise in the pairwise comparisons. The hyperplexing data obtained at the MS2 level appeared to have a weaker correlation with the MS1 data in comparison to the SPS-MS3 experiments. This is likely because of the higher presence of interfering co-isolated ions that distort the fractional labeling measurements. In general, there also appeared to be more experimental noise in the TMT-hyperplexing data obtained from SILAC-labeled precursor peptides in comparison to SILAC-unlabeled data. This is likely because of increased error in the regression analysis. Unlike unlabeled precursors that have a uniform final baseline signal near zero that can be and all other peptides mapped to the protein (transparent blue circles). The green line indicates a fit to a single exponential equation with a non-zero baseline indicated by a dashed line and double-headed arrow. D, The kinetics for appearance of GEYDVTVPK analyzed at the MS2 level by TMT-SILAC hyperplexing. The MS2 spectrum indicates the relative intensities of TMT reporter ions for the SILAC-labeled peptide ion. The plot indicates intensities of peaks relative to the fitted intensity of the exponential curve (orange line) at infinite time (t ϱ ). The dashed line indicates the intensity of the t0 sample (TMT-126) and represents the initial baseline of the curve. E and F, Data analyzed as in C and D, respectively, except at the MS3 level using SPS-MS3. assigned as a constant in the regression analysis, the signal of labeled peaks at infinite time is undetermined and must be treated as a variable fitted parameter. Our extended labeling time-course allowed many proteins to reach near complete labeling during the experiment and thus the intensity of fully labeled peptides could be determined with high confidence. However, for proteins whose labeling had not reached near completion after 14 days of labeling, a considerable error could have potentially been introduced in determination of their labeling rate constants.
Global Impact of Cellular Quiescence on Protein Synthesis and Degradation Rates-In the above experiments, the apparent rate constants for clearance and accumulation are measured by analysis of fractional labeling. In other words, the intensities of unlabeled and labeled reporters were internally normalized with respect to fully unlabeled or labeled The solid line indicates a fit to a single exponential equation for measurement of clearance rate constants. For quiescent cells, the clearance rate constant equals the time-points, respectively. Because the kinetics of protein degradation typically follows first order kinetics whereas protein synthesis is a zero-order process with respect to protein concentrations, the apparent rate constants for fractional labeling for both the disappearance of unlabeled reporters and appearance of labeled reporters correspond to the first-order rate constant for degradation (k deg ). However, the TMT-SILAC hyperplexing methodology can be readily modified to additionally measure relative changes in synthesis rates (k syn ) between two experimental conditions under steady-state conditions. By conducting TMT-SILAC hyperplexing analyses of labeling time-points obtained from two conditions, ratios of steady-state protein levels, as well as protein degradation rates can be measured in a single experiment. In turn, these two parameters can be used to calculate the ratio of k syn values in the two experimental conditions (Fig. 5, Experimental Procedures). We employed this approach to compare global changes in both protein degradation and synthesis as cells transition from a proliferating to a quiescent state.
Proliferating and quiescent cells have widely different metabolic demands. Much of the energy expenditure in dividing cells is devoted to accumulation of biomolecules necessary for formation of new daughter cells (31). Exempt from this anabolic requirement, quiescent cells typically downregulate the synthesis of proteins at the level of transcription and translation (32)(33)(34). This homeostatic adjustment ensures that the overall rate of protein production does not exceed the rate of clearance and protein concentrations are maintained at steady-state levels in quiescent cells.
Moreover, a newly quiescent cell can potentially counteract the reduced requirement for protein accumulation not only by reducing rates of protein synthesis, but also by increasing rates of protein degradation. The relationship between cellular growth and protein degradation has previously been investigated by a number of classic radioactive tracer studies. By analyzing the degradation of the bulk cellular protein population, some of these studies have reported significant changes in degradations kinetics in quiescent cells (35)(36)(37)(38)(39), whereas others have failed to observe such an effect (40 -45). Here, we addressed this question by simultaneously measuring proteomewide changes in synthesis and degradation rates between dividing and quiescent cells using TMT-SILAC hyperplexing.
We combined four SILAC-labeled time-points for both proliferating and quiescent cells in a single hyperplexed analysis (Fig. 5). The experiment allowed us to simultaneously measure relative ratios of steady-state protein levels and degradation rates between quiescent and proliferating cells. The relative ratios of steady-state protein levels were obtained by com-paring reporter tag intensities for SILAC-unlabeled (t0) peptides between dividing and quiescent cells. Protein degradation rate constants were determined by measuring the first order rate constant for clearance of SILAC-unlabeled peptides and correcting for the contribution of protein dilution by cell division in dividing cells. In turn, the experimentally measured relative ratios of steady-state protein levels and degradation ratios were used to calculate the relative ratio of synthesis rates according to the model presented under Experimental Procedures. The relative log2 ratios of steadystate protein levels ( ͓ protein stead yϪstate ͔ quiescent ͓ protein stead yϪstate ͔ proliferating ), degradation rates ( k deg quiescent k deg proliferating ) and synthesis rates ( k syn quiescent k syn proliferating ) were measured for 1277 protein and are tabulated in supplemental Table S5.
As expected, our data indicates that the synthesis rate of proteins are globally decreased in quiescent cells. However, somewhat unexpectedly, our data also provides clear evidence that for a large cross-section of the proteome, protein degradation rates are enhanced in quiescent cells (Fig. 5D). Interestingly, the enhancement of degradation rates, but not reduction of synthesis rates, is particularly pronounced in longlived proteins (Fig. 5D, see Discussion). Together, our data indicate that quiescent fibroblasts compensate for the reduced requirement for protein accumulation by both decreasing global rates of synthesis and enhancing rates of protein of degradation. The results reported here provide the first quantitative global census of these two effects in quiescent cells.

DISCUSSION
The proof-of-principle experiments described above validate a potentially powerful approach for global analyses of proteome dynamics. By combining dynamic SILAC and TMT technologies, multi-time-point kinetics of protein turnover can be analyzed in a single LC-MS/MS run, significantly reducing the cost and complexity of dynamic proteomic experiments. For example, here, we were able to analyze the kinetics of protein clearance at ten distinct time-points in a single 3-hour LC-MS/MS run with nearly no missing time-points for each analyzed protein. If each time point were to be analyzed independently in equivalent dynamic SILAC experiments, more than 24 h of instrument time would be needed and most proteins would not have been quantifiable at all time-points. However, our analyses also indicate that in comparison to conventional dynamic SILAC experiments, this improvement in time-resolution and instrument time comes at the cost of quantitative precision and proteome coverage. These limitations are caused by a number of factors. First, the co-isolation degradation rate constant. For proliferating cells, the rate of cell division (0.40 d Ϫ1 ) is subtracted from the clearance rate in order to determine the degradation rate (dashed blue line). D, The distribution of log2 ratios of synthesis and degradation rate constants between quiescent and dividing cells. The histograms are shown for all proteins, long-lived proteins (half-lives greater than 3 days in proliferating cells) and short-lived proteins (half-lives less than 3 days in proliferating cells). The dotted lines on the rightmost plot indicate the distribution of long-lived proteins for comparison. and fragmentation of interfering ions introduce error into the quantitation of TMT reporter ions. Although this effect can be somewhat mitigated by accounting for a constant baseline signal in the kinetic analysis and utilization of the SPS-MS3 protocol, in general we found more experimental noise in the labeling kinetics of reporter MSn ions in comparison to MS1 ions. Second, TMT labeling leads to a reduction in observed peptides because of a general increase in charge state on the labeled peptides. Finally, because of the small mass differences between some TMT 10plex labels, a relatively high resolution is needed to ensure peak separation, increasing the scan time and reducing the number of MS2 scans in the run. Future improvements to the TMT-SILAC hyperplexing methodology will require the development of further strategies for maintaining proteome coverage while minimizing the co-isolation of interfering ions with targeted precursor ions. Until then, the decision to employ conventional dynamic SILAC or TMT-SILAC hyperplexing for analysis of proteome dynamics will likely hinge on the relative importance of temporal resolution and experimental cost versus quantitative precision and maximal proteome coverage.
When conducted under steady-state conditions, apparent rate constants of fractional isotopic labeling typically provide information on protein degradation kinetics (e.g. half-lives). However, by co-analyzing samples collected under two or more environmental conditions, TMT-SILAC hyperplexing has the potential ability to provide proteome-wide measurements of relative ratios of steady-state protein levels in addition to degradation kinetics. Together, these two measured parameters can be used to measure relative changes in protein synthesis rates under steady-state conditions. To demonstrate this capability, we measured differences in protein synthesis and degradation rates between proliferating and dividing cells. Our data indicate that quiescent cells globally decrease rates of protein synthesis and increase rates of protein degradation.
Decreased rates of protein synthesis are a well-described feature of the quiescent state and reductions in overall rates of transcription and translation have been described in bacteria, yeast and higher eukaryotes (31)(32)(33)46). In stationaryphase S. cerevisiae, where this phenomenon has been best described, the inhibition of the mTOR pathway appears to be a key player in inhibiting protein synthesis, although other contributing pathways, such as protein kinase A and C are also known to play a role (47). The inhibition of the mTOR pathway results in a number of downstream effects that promotes the quiescent state, including the shutdown of ribosomal biogenesis required for protein synthesis (48).
Our data also demonstrate that rates of protein degradation are enhanced in quiescent cells and that this effect is most pronounced for long-lived proteins. In dividing cells, two general processes contribute to protein clearance: cytoplasmic dilution by cell division and hydrolysis by cellular degradation pathways (23,49). Indeed, if we consider a simplified kinetic model of protein expression (see Experimental Procedures), the rate constant of protein clearance is the mathematical sum of the rate constants for these two flux pathways. As has been previously noted (23), degradation is the predominant route of clearance for short-lived proteins, whereas long-lived proteins are predominantly cleared from a cell through the process of proliferation. Thus, as cells transition from a dividing to a non-dividing state, it is the clearance kinetics of long-lived proteins that are particularly impacted. If no compensatory response is initiated, this phenomenon would result in a global increase in cellular concentrations of long-lived proteins relative to short-lived proteins in quiescent cells. By selectively increasing rates of degradation for long-lived proteins, quiescent cells can counter this potential proteostatic disruption and ensure that their production does not exceed cellular demand.
It has previously been shown that basal macroautophagy is up-regulated in quiescent fibroblasts (50). Given that macroautophagy is known to preferentially degrade long-lived proteins (19,51), its activation may form the molecular basis of the proteome-wide trends observed in our study. Further studies are required to assess this hypothesis and investigate the role of other potential pathways in quiescence-induced protein degradation. It is also important to note that the proteostatic response to quiescence may vary in different cell types. Indeed, the quiescent state of dermal fibroblasts is somewhat unique in that fibroblasts maintain a rather active metabolic state in order to synthesize and secrete extracellular matrix proteins such as collagen fibers that form the basement membrane of connective tissues (31). Conversely, many other cell types enter a quiescent state primarily as a response to nutritional deprivation. Whether different cell types differ in their global proteostatic response to quiescence remains to be determined. More generally, the methodology described in this study may provide a useful proteomic approach for globally investigating biological processes in which the kinetics of protein synthesis and degradation are jointly modulated.
Associated Content-Spectral search results, SILAC and TMT quantitations, and rate measurements at peptide and protein levels are provided in tabular format in the supplementary information.
Accession Numbers-All raw and processed data are available at ProteomeXchange Consortium via the PRIDE database (accession number: PXD004725).