Dynamic proteomic profiling of a unicellular cyanobacterium CyanotheceATCC51142 across light-dark diurnal cycles

Qualitative and quantitative profiling of the proteome

Cyanothece have been shown to store excess C and N in the form of glycogen and cyanophycin granules and to utilize these compounds during starvation [11]. Our experimental strategy involved growing cultures under the diel regime (12 h L: 12 h D) in ASP₂ medium without NaNO₃ (nitrogen starvation) for seven days prior to adding the labeled leucine. (Figure 1). Sampling of partially-labeled Cyanothece ATCC51142 cells was initiated thirteen hours after the addition of the labeled leucine, and samples were collected at 4 h intervals across 13 time points for a time-course analysis (Figure 1). Extractable proteins were divided into soluble and insoluble fractions for each collection [16] and subjected to high resolution LC-MS/MS analysis after tryptic digestion and SCX separation (Figure 1). Allowing < 5 ppm average mass deviations and 1% FDR on both peptide and protein levels, we were able to identify 19327 leucine-containing peptides irrespective of label information, which mapped to 2202 proteins with at least two peptides per protein, representing 47% of the entire predicted proteome (Figure 2A).

In cyanobacteria, the synchronization of cellular functions with appropriate phases of the diurnal cycle must originate, in part, from the dynamic changes in protein expression, which is the net result of the contribution from a new protein synthesis in competition with protein degradation at any given time [26]. However, under our experimental strategy, the rates of degradation were assumed to be identical due to the fact that these versions are chemically indistinguishable molecules, and are indeed inseparable for degradation. As a consequence, changes in the abundance of a specific protein can be attributed due to changes in new protein synthesis which can be determined by the measure of label incorporation. Accordingly, accumulation of the label into proteins would be proportional to the rate at which they are synthesized, and degraded. Under our experimental set-up, all newly synthesized proteins would increase in label through the time of the experiment, but at different rates. So, proteins with rapid turnover accumulate label much faster than the low turnover proteins. We used Relative Isotope Abundances (RIA) [21, 22, 24] to monitor the temporal changes in protein expression across time points. RIA provides a measure of the fraction of the labeled peptides in the total (labeled + unlabeled) peptide pool for a particular protein, and this value can be used to determine actual new protein synthesis at any given time point. Achieving a complete (saturated) labeling (RIA ~1) as assumed for pooled classical SILAC based measurements, was not the goal of this study, however, the concept of partial (incomplete) labeling with short labeling windows (dynamic metabolic labeling) can be a powerful approach to monitor and visualize the proteome dynamics. In incomplete labeling, each protein contains both labeled and unlabeled amino acid, and its tryptic digestion generates both labeled and unlabeled peptides. At each time point, proteins were identified, and the changes in the partition of peptide ion current between the labeled and unlabeled peptides in LC-MS measured by RIA provide information on the real-time protein synthesis. For example, with an RIA of 0.55 at D1, NifH (cce_0559) has at least 55% of this protein newly synthesized over the preceding 13 hours at that time point (Additional file 1, Table S2). However, as proteins can be synthesized from amino acids (i. e. leucine) that can be derived by de-novo synthesis, these measured values likely represent the lower limit of the true turnover [18].

Using this procedure, we identified 2837 individual peptides with the pair-wise information for both the heavy and corresponding natural light version of the peptides (Additional file 2, Table S1) which mapped to 721 proteins (13.8% of the predicted proteome) (Figure 2A and Additional file 1, Table S2). Figure 2B shows the distribution of these labeled proteins compared with the total protein observations (labeled + unlabeled) whereas the distribution of the labeled peptides and their corresponding proteins across all time points are shown in Additional File 3, Figure S1A. Even though the analysis was performed at a uniform peptide amount across all samples, we saw a distinct decreasing trend of total identifications across the entire time course (Figure 2B), which is more obvious during the light phases. This downward trend appeared to be predominately influenced by the decrease in unlabeled protein identifications over time.

The differentially labeled proteins with their individual patterns of RIA across the time-course were subjected to hierarchical clustering generating the dendrogram as represented in Figure 2C. These data show a clear increase in the incorporation of metabolic label per protein over time (green to red), suggesting that isotope incorporation is driven by the cell growth. However, the proportion of the labeled amino acid varied among the individual proteins (RIA values) suggesting that isotope incorporation is also dictated by intracellular protein turnover.

Functional categorization of the labeled proteome from KEGG pathway [27], and Gene Ontology (GO) analysis showed that the actively synthesized labeled proteins belong to diverse individual biochemical pathways such as CO₂ fixation and assimilation, N₂-fixation and N assimilation, glycogen synthesis, glycogen metabolism, glycolysis and oxidative PP pathways (Figure 2D), and that many of these were maximally expressed at different phases of the light-dark cycle. For example, the majority of proteins involved in glycolysis, TCA cycle and PP cycle usually showed higher expression in the dark whereas those involved in CO₂ fixation, heme, chlorophyll and phycobilins biosynthesis showed higher expression in the light cycle (Figure 2D). Again, proteins involved in N₂-fixation, protein degradation, glycogen synthesis and metabolism, heme and chlorophyll biosynthesis showed relatively higher levels of isotope incorporation (black or red on the heat map) as compared with other functional categories. Interestingly, a total of 182 (25%) proteins with unknown functions were differentially labeled in between the light-dark cycle suggesting their potential roles in cellular metabolism, and highlight the importance of systematic characterization of their functions for better understanding of the physiology and cellular metabolism of Cyanothece.

Comparison between transcriptional and proteomic data

Details of the transcriptional regulations of biological processes in Cyanothece and other unicellular cyanobacteria have been described in previous studies [12–15, 28, 29]. A robust approach to determine both cycling Cyanothece ATCC51142 transcripts and cycling proteins from near identical experiments has been previously reported [12, 16] which data was used in the current study for comparison. Indeed, a major focus of the current study was to investigate exactly how these cycling transcripts/proteins would be reflected when viewed through the real time metabolic synthesis of proteins. There was no additional cycling analysis performed on the current dataset, other than to map the detected proteins onto the previously published and characterized data. Although transcriptional and translational events are generally correlated, concordances and discordances exist [30]. Previously [16], we identified 250 proteins oscillating between light-dark cycles, of which 175 matched to the cycling genes [12]. The oscillating proteins were determined based on Fourier score [29], fold change (≥ 1.5) between spectral counts and computed P-value of < 5% using two sample t-test [16]. Of the 721 labeled proteins identified in this study, 425 proteins (58.9%) matched with the previously characterized cycling transcripts (Figure 3). Similarly, out of the 250 oscillating proteins identified from label free approach [16], 145 matched with the labeled proteins (Additional file 3, Figure S1B). This protein subset represents those under active synthesis and having correspondence to well-defined cyclic transcripts. Of potential interest are the 296 (41.1%) actively synthesized proteins having no correspondence to cycling transcripts, and possibly represent those that are likely regulated at a post-transcriptional (translational or post-translational) level. For example, circadian clock protein KaiC1 (cce_0422), which rhythmically phosphorylates and dephosphorylates in-vivo and serves as a timing cue for assembly/disassembly of other circadian clock proteins such as KaiA and KaiB [31], was actively synthesized in the dark cycle but had no distinct cycling transcript (Additional file 1, Table S2) [12]. Functional classifications showed that large proportions of the labeled proteins were of unknown function in both overlapping and non-overlapping sub-sets. Among the known functions, photosynthesis, respiration, CO₂-fixation, N₂- fixation, carbohydrate metabolism, cell division and protein translation showed higher representation in the sub-set that overlap with the cycling transcripts whereas DNA replication and repairs, transport and binding, nucleotides metabolisms, and tRNA synthesis and modification, and regulatory functions showed higher representation in the non-overlapped sub-set. Interestingly, all the 32 ribosomal proteins identified in this study matched with the corresponding cycling transcripts [12]. Such proteins are obviously integral to the ribosome and are involved in the translation of transcripts encoded in cellular genomes, and may exert a level of developmental control over cell/tissue fate and function [32]. Thus, major differences of proteins for various functional classes were apparent indicating that diverse groups of proteins are regulated differently.

Nitrogen fixation and assimilation

Expression of nitrogen fixation proteins NifH, D, and K were strongly cyclic across both diurnal cycles and limited to the dark cycle [16] (Figure 4A, first subset), with an actual loss of detection within the light compartment. NifB (cce_0554), a sensitive target for assessing the regulation of N₂-fixation [15], was also abundant in the dark whereas NifW (cce_0568), a small protein with 116 residues (13.3 kDa) located in the 35 nitrogenase gene clusters, was abundant near the end of the first dark transition, but not observed in the second dark phase. This was also true for NifX (cce_0565) and NifS (0557). Interestingly, out of the 8 labeled domain of unknown function (DUF) proteins, DUF269 (cce_0566) was distinct, showing cycling expression with higher incorporation in the dark (Additional file 1, Table S2), in agreement with transcript data [12]. Although, the biological function of DUF269 is not clear, its expression pattern and genetic location within the 35 nif gene cluster clearly suggests a role in N₂-fixation.

Under nitrogen limitation, cyanobacteria preferentially utilize alternative N sources [33]. GlnA (cce_4432) and GlsF (cce_4353) involved in nitrogen assimilation had increased synthesis in the dark (Figure 4A, second subset). Mechanistically, ammonium (NH₄⁺) in the cell is trapped by GlnA [34] and converted to glutamate via GS (GlnA)-GOGAT (GlsF) cycle, and incorporated into other amino acids through the action of aminotransferase or transaminase [35]. GlnB (cce_1775), which activates the transcription of many nitrogen regulated genes by regulating kinase and phosphatase activities [36] also showed increased synthesis in the dark (Figure 4A, second subset). On the other hand, proteins involved in amino acid biosynthesis such as chorismate synthase; AroC (cce_4200) and acetolactate synthase large subunit; IlvB (cce_4227) showed higher synthesis during the light cycle (Additional file 1, Table S2). AroC catalyzes chorismate production, which is a key intermediate of phenylalanine, tyrosine, and tryptophan biosynthesis whereas IlvB catalyzes leucine, isoleucine and valine biosynthesis. Differential labeling of these proteins provide evidences of their active synthesis and role in diurnal growth.

Photosynthesis and respiration

Carbon fixation machinery requires proper assembly and functioning of carboxysome [1]. Our data showed higher expression of IcfA1 (cce_2257) (Figure 4B, top panel) and bicarbontate transport system substrate-binding protein, CmpA (cce_0305) in the light cycle, and sodium dependent bicarbonate transporter; SbtA (cce_2939) in the dark cycle (Additional file 4, Figure S2). CmpA is highly homologous to the nitrate transport protein NrtA, and participates in the co-transport of bicarbonate, Ca²⁺ and Na⁺ or that Ca²⁺ and Na⁺ act as a cofactor in bicarbonate transport [37]. Glycolate oxidase (cce_3708), a key enzyme to convert 2-phosphoglycolate to PG and O₂ [38] showed higher amount of label incorporation in the light cycle (Additional file 1, Table S2). 2-Phosphoglycolate is produced by the oxidase activity of rubisco when the O₂-CO₂ ratio increases. Excitation of PSII during the light period is expected to increase the O₂-CO₂ ratio and expression of glycolate oxidase may constitute an important cellular adaptation under the light period [39]. So, higher synthesis of several carboxysome proteins (RbcSL, IcfA1, CcmM, CcmK1, CcmK2, and CcmK4) together with ParA (cce_2448) and RbcR (cce_3731) in the light cycle might be an indication of active carboxysome functioning under the light period (Additional file 4, Figure S2 and Additional file 1, Table S2). In Synechococcus elongatus PCC 7942, a parA-like gene (Synpcc7942_1833) is known to control carboxysome organization [1], and ParA protein was also actively synthesized in the light cycle for Cyanothece ATCC51142, providing clues of cellular processes that control carboxysome assembly and function.

Expressions of PSI and PSII proteins were not strongly coupled with the light-dark cycles (Additional file 4, Figure S2). In particular, in PSII, Psb28 (cce_1599), PsbE (cce_1307) and PsbA3 (cce_0267) were more abundant in the dark while other PSII proteins (PsbD1, PsbC, Psb28-2, PsbO, Psaf, PsbP) were more abundant in the light cycle. This disparity might be a result of continuous assembly and disassembly of PSI and PSII proteins due to frequent light-mediated damage under normal cellular function. Phycobiliproteins, which play an important role for light harvesting and distributing energy to PSI and PSII, usually had higher levels of synthesis in the light phase (Additional file 4, Figure S2). In contrast, CoxB1 (cce_1977), a subunit of the terminal oxidase in the respiratory transport chain was predominant in the dark cycle (Figure 4B, second panel). Respiration has an important role in diazotrophy in cyanobacteria [7, 40], and CoxB1 is likely to play a role in scavenging free O₂ in the cytoplasm to prevent nitrogenase deactivation [7]. The transcripts of coxA were shown to increase in the dark cycle both under N₂-suffcient and N₂-deficient conditions [41]. Interestingly, soluble electron carrier ferredoxin-NADP oxidoreductase, PetH (cce_0994) showed increased synthesis in the dark cycle whereas plastocyanin PetE (cce_0590) showed increased synthesis in the light cycle (Figure 4B, third panel) suggesting that PetH might be involved in the respiratory and PetE in the photosynthesis electron transport.

Central intermediary metabolism

Efficient cellular metabolism relies on the compartmentalization of enzymatic reactions [1]. In cyanobacteria, formation of glycogen during photosynthesis and their metabolism during N₂-fixation is well known [7]. Consistent with this, synthesis of glycogen synthase; GlgA1 (cce_3396) peaked at the first light cycle (L1-3) with no detection at later time points (Figure 4A, third subset). Similarly, synthesis of 1, 4-glucan debranching enzyme; GlgB2 (cce_4595), which is involved in converting glucose to glycogen also peaked in the first light cycle and remained consistently high in the second light cycle (Figure 4B, 4th panel). In contrast, glycogen phosphorylase; GlgP1 (cce_1629) usually had higher levels of synthesis in the dark compared to the light cycle, suggesting higher rate of glycogen metabolism during N₂-fixation (Figure 4B, 4th panel) [7, 16]. Overall, we observed active synthesis of the majority of glycolysis and PP pathway enzymes at varied time points during the dark cycle (Figure 4A). Glucose 6-phosphate dehydrogenase, Zwf (cce_2536) redirects C flow and initiates the oxidation of glucose concomitantly to sugars and reducing equivalents in the form of NADPH, and represents the essential catabolic route for providing reductant for N₂-fixation and respiration in Cyanothece ATCC51142 [12, 16]. There was a slight trend in the synthesis rate of Zwf (Figure 4B, 5th panel) being increased in the dark, however there is no significant different in abundances between light-dark cycles. In cyanobacteria, the TCA cycle is incomplete [42]. Isocitrate dehydrogenase (Icd) converts isocitrate to 2-OG, which is utilized for the biosynthesis of glutamate via GS-GOGAT cycle [43]. Concurrently, we observed a slight trend in the synthesis of Icd (cce_3202) being increased in the dark and GabD, another enzyme known to participate in glutamate metabolism for the conversion of 2-OG to succinate [44] in the light cycle (Figure 4B, 6th panel), and might suggest that cells utilize different substrates as C and N sources. Increased expression of icd transcripts under N₂-fixing conditions has been observed in other cyanobacteria [43].

Synthesis of cyanophycin synthetase; CphA (cce_2237) (Figure 4B; bottom panel) was not coupled with the light-dark cycle, but cyanophycinase, CphB (cce_2236), which breaks down cyanophycin into arginine and aspartic acid, showed higher synthesis in the light cycle (Figure 4B, bottom panel). Arginine and aspartic acid can be utilized for the synthesis glutamate by combining with 2-OG and ammonia [45] (Figure 5) and might be utilized via the urea cycles as reported in Synechocystis PCC6803 [11] to obtain both alternative C and N sources..

Based on the quantitative information of new protein synthesis, we provided an overview of the temporal separation of cellular processes and metabolism in Cyanothece ATCC51142 during the light-dark cycle (Figure 5). Pathways were divided according to the diurnal period in which most of the proteins in a group are up-regulated. The timing was set based on the time course when the mean RIA of all members in a functional group (for example N₂- fixation) for a given time-point was higher than the mean of all the time points. Different enzymes involved in multiple pathways linking carbon, organic acid and amino acid metabolism showed increased synthesis at different phases of the light-dark cycles suggesting a robust cellular programming in Cyanothece ATCC51142 to undergo diurnal cycles. Enzymes involved in carbohydrate metabolism (glycolysis, TCA cycle and oxidative PP cycle) usually showed higher synthesis in the dark and enzyme involved in CO₂ fixation, stress and redox regulation, and co-factors synthesis showed higher synthesis in the light cycle (Figure 5). Importantly, enzymes involved in fatty acid and lipid metabolism including GlpK (cce_0973), FadD (cce_1133), FabD (cce_1583), PlsX (cce_1585) and GlpD (cce_2610) had a slight trend for increased synthesis either in the light cycle or during the transition from light to dark cycle. This might suggest that fatty acids and lipids are important sources of energy for many cellular functions. Additionally, photo-oxidative stress during active photosynthesis generates reactive oxygen species (ROS) and triggers redox signaling, apparently expressing enzymes of cellular protection and detoxification including Gor, Trx, Gst and SodB.

Crosstalk between carbon metabolism and other intermediary pathways

The ability to link and quantify enzymes that link carbon, organic acid and amino acid metabolism is crucial to understanding metabolic regulation at multiple levels. Our analysis enables the linking of carbon metabolism with other central metabolism including lipids, fatty acids, amino acids (leucine and isoleucine), nucleic acids (purines and pyrimidines), heme and chlorophyll metabolism (Figure 5).

The heme biosynthetic pathway is fundamental for phycobilin and chlorophyll biosynthesis [46]. In parallel with phycobiliproteins, those involved in pigment biosynthesis also showed increased synthesis in the light cycle (Additional file 4, Figure S2) suggesting their role for chlorophyll biosynthesis and for the formation of PSI, PSII and cytochrome b6f complex [47]. Integration of this data suggests a novel pathway-crosstalk between central metabolism and the heme and chlorophyll biosynthetic pathway (Figure 6). Significant expression of GabD as mentioned before (Figure 4B, 6th panel) might possibly suggest the existence of an alternative route for the conversion of 2-OG to succinate (Figure 6). So, a portion of 2-OG might be channeled to produce succinate, and might be utilized as an alternative carbon source. Although, this pathway has been described in plants [48], to our knowledge, a similar short-cut shunt has not been described in cyanobacteria yet. Furthermore, significant up-regulation of PHK (cce_3607) during the dark cycle is striking (RIA = 0.62 at D2-1) (Figure 6), as this enzyme catalyzes the first step reaction to convert xylulose- 5P to pyruvate. In glycolytic-PP pathway, conversion of xylulose-5P to pyruvate generates an average of 1.67 moles of ATP per mole of xylulose-5P [49]. In contrast, conversion of xylulose-5P to pyruvate and acetyl-P by PHK pathway generates an average of 2 moles of ATP, and additional 1 mole of ATP is further produced by the conversion of acetyl-P to acetate. Therefore, the existence of a more energy efficient PHK pathway in cyanobacteria might be a cellular strategy and play an important role for improved cellular metabolism under limited condition.

Chaperones, heat shock and stress responsive proteins

Quantitative profiles of chaperons and stress response proteins peaked at different stages of light-dark cycle (Additional file 4, Figure S2). While larger proportions showed relatively higher abundance in the light, others such as superoxide dismutase; SOD (cce_1620), glutathione S-transferase; Gst1 (cce_0174) and chaperonin 2; GroEL2 (cce_3330) showed higher abundance in the dark cycle. So, regulation of these proteins is presumably controlled by the metabolic state of the cell rather than the direct response to the light-dark cycle. ATP dependent Clp proteases also peaked at different time points (Additional file 4, Figure S2). Notably, regulatory subunits ClpC1 (cce_4247) and ClpB1 (cce_1941) showed higher synthesis in the dark but proteolytic subunits such as ClpP1 (cce_0641), ClpP2 (cce_4439), ClpP3 (cce_1911) and ClpP4 (cce_1910) showed higher synthesis in the light. Protein degradation plays a vital role in regulating cellular enzymes and regulatory proteins in cyanobacteria [50]. It is also responsible for targeting and efficient removal of damaged or otherwise functionally abnormal proteins such as oxygen sensitive nitrogenase before they reach levels toxic to metabolic process. This selective proteolysis in cyanobacteria is performed by different ATP dependent proteases, and their active synthesis suggests a robust targeted protein degradation mechanism, an agreement with the increased protein aggregation during the transition from dark to light cycle or the vice-versa and removal of aggregated proteins by different proteases. Additionally, carboxy terminal protease (cce_1372) and amino terminal peptidase, PepN (ce_0144) were also actively synthesized (Additional file 4, Figure S2). These peptidases are involved in proteolytic processing of proteins into mature forms [51], where carboxy terminal protease is responsible for cleavage of C-terminus residues, and PepN is involved in ATP-independent steps in the cleavage of amino terminus residues. As photosynthesis inherently produces strong oxidants and reductants, and utilizes multiple redox active intermediate electron transport chain, maintenance of redox balance is particularly challenging in photosynthetic microbes [52]. So, active synthesis of Clp proteases, chaperones and antioxidants suggest an extensive cellular response mechanism in Cyanothece ATCC51142 to deal with production of reactive oxygen species (ROS) and protein degradation during photosynthesis and N₂-fixation.

Culture conditions and labeling experiment

Stock cultures of Cyanothece ATCC51142 were maintained in ASP₂ medium [3] with NaNO₃ under conditions of continuous light at 50 μmol photons.m^-2.s^-1 [12]. For the experiment, Cyanothece ATCC51142 cells were grown in ASP₂ medium [3] without NaNO₃ at 30°C under 12 h light/12 h dark cycles and 50 μmol photons.m^-2.s^-1 of light for 7 days (chlorophyll concentration of 0.97 μg/ml of cells) before adding 20 mg/l of [¹³C₆ ¹⁵N₇]-leucine (U-¹³C₆, 98%; ¹⁵N₇, 98%; Cambridge Isotope Laboratories Inc., MA). Cultures were grown for a 12 h period in the light to allow for sufficient isotope incorporation before first sampling at 1 h after the onset of the dark period. Samples were collected at 4 h intervals. A total of 13 samples were collected for analysis with the last sample collected at 1 h into the third diurnal cycle (D3-1; Figure 1).

Sample preparation

Soluble and insoluble fractions from total cell extracts were prepared as described before [11, 16]. Protein concentrations were determined using the bicinchoninic acid (BCA) assay (Pierce Chemical Co., Rockford, IL). Soluble proteins were denatured and reduced in 8 M urea and 5 mM tributylphosphine (Sigma-Aldrich) at 37°C for 60 min. The insoluble fraction was treated identically, except for the addition of 1% CHAPS for 45 min prior to digestion. All samples were then diluted 5 times in 25 mM ammonium bicarbonate prior to tryptic digestion using sequencing-grade modified porcine trypsin (Promega, Madison, WI) at a 1:50 (w/w) trypsin-to-protein ratio for 5 h at 37°C. Peptides were desalted using a 1 ml C-18 SPE column (Supelco, St. Louis, MO) [53]. Samples were separated using strong cation exchange chromatography (SCX) with a PolySulfoethyl A, 200 mm × 2.1 mm, 5 μM, 300-Å column and a 10 mm × 2.1 mm guard column (PolyLC, Inc., Columbia, MD) at a flow rate of 0.2 ml/min as previously described [54, 55]. The peptides were re-suspended in 900 μl of mobile phase A (0.2% FA in water), and separated on an Agilent 1100 system (Agilent, Palo Alto, CA) equipped with a quaternary pump, degasser, diode array detector, Peltier-cooled auto-sampler and fraction collector (both set at 4°C). A total of 24 fractions were collected for each sample (soluble or insoluble protein) and pooled in such a way as to make 6 fractions for soluble protein and 4 fractions for insoluble protein for subsequent MS analysis.

Reversed Phase Liquid Chromatography (LC) Separation and Tandem MS (MS/MS) Analysis

Samples were analyzed using high throughput hybrid LC-MS/MS of the trypsin-digested proteome samples [16]. Six (6) μl of the tryptic peptides containing 3 μg of the total peptides from each SCX fraction out of 60 μl per fraction were analyzed by capillary LC (150 μm i.d. × 360 μm o.d. × 65 cm capillary; Polymicro Technologies Inc., Phoenix, AZ) under constant pressure of 5000 psi coupled with a Finnigan LTQ Iontrap Mass Spectrometer (ThermoFinnigan, San Jose, CA) and an electrospray ionization source manufactured in-house as described in [56]. The LC column was prepared by slurry-packing 3-μm Jupiter C₁₈ bonded particles (Phenomenex, Torrence, CA) into a 65-cm long, 75-μm inner diameter fused silica capillary (Polymicro Technologies, Phoenix, AZ). After loading peptides onto the column, peptides were separated for a 100 min run at a flow rate of ~1.8 μl/min using the following settings: 100% mobile phase solvent A (0.2% formic acid in water) for 20 min, followed by a linear gradient from 0 to 70% of solvent B (0.2% formic acid in 100% acetonitrile) over 80 min before reverting to 100% solvent A. Each full MS scan (m/z 400-2000) was followed by collision-induced MS/MS spectra (normalized collision energy setting of 35%) for the 10 most abundant ions per duty cycle selected for fragmentation. The dynamic exclusion window was set to 30 s, the heated capillary was maintained at 200°C, and the ESI voltage was held at 2.2 kV.

Data analysis and relative isotope abundance (RIA) measurement

The MS/MS data were processed using MaxQuant software [57] version 1.1.1.2. The MaxQuant package integrates a suite of software packages for peak list generation, labelled peptide pairs, SILAC- and XIC-based quantitation, false discovery rates [58], and peptide identification based on Mascot [42] search [18]. The step by step details of this procedure can be found in reference [57]. Briefly, MaxQuant executes data analysis in 4 sequential steps. First step performs 'pre-processing and quantitation' to identify all the peptide pairs without knowing the identity, and then performs Mascot search in the second step. Mascot search was performed using version 2.2.04 against the concatenated forward and reverse versions of the Cyanothece ATCC51142 protein sequence databases (updated on June 21, 2010) downloaded from NCBI http://www.ncbi.nlm.nih.gov/protein that contained 5304 protein sequences. Mascot search parameters included: enzyme trypsin, one missed cleavage, oxidation of methionine and ¹³C¹⁵N-leucine (7 Da) as variable modifications. The maximum allowed mass error tolerances for MS and MS/MS scans were 5 ppm and 0.5 Da, respectively. Peptides with ≥6 amino acid residues were only allowed, and proteins that had at least two peptides with ion scores > 20 were considered for identification and quantitation. A target decoy database approach using forward-reverse sequences and contaminants was used to determine false discovery rate (FDR) and FDR threshold was set to the default values of 1% for both the peptides and the proteins and posterior error probability (PEP) of 0.5%. The identified and quantified peptide pairs are then visualized into two-dimensional (2D) and three-dimensional (3D) peak formats. 2D peaks are generated based on local maxima of the intensity as a function of m/z, centroid and assembled into 3D peak hills over the m/z-retention time plane. Two peaks in the neighbouring scans are connected if their centroid m/z positions differ by < 7 ppm. At least two centroids should match to form a 3D peak and the centroids that cannot be matched to centroids in the two previous or the two next scans are discarded. This removes a major fraction of the noise in form of peaks not reproducible in time. Additionally, intensity profiles of the pairs should have a sufficient overlap in retention time by requiring them to have a uncentered Pearson correlation of > 0.6. Then, the intensity profiles are smoothed over retention time applying a window mean filter of ± 1 scan width and checked for local minima. Intensities of heavy and light peptides thus generated were used to determine the Relative Isotope Abundances (RIA) described previously [21, 23, 24] as:

RIA_t = I_H/(I_H + I_L), where I_H and I_L represent the intensities of heavy and light peptides, respectively at time t. RIA measures the relative abundance of the label peptides (I_H) to the total peptide pool (I_H+I_L) in the cell. So, rapidly synthesized proteins accumulate label at a faster rate than the low turnover proteins, and proportionally have higher RIA values. The protein RIA at any particular time point was calculated based on the averaged peptides RIAs per protein. As labelled peptides contain the mixture of both fully and partially labelled peptides (in case of multiple Leu containing peptides), both peptide mixtures were used to calculate the RIA. In each time point, a labelled protein was identified by the pair of stable isotope labelled and the corresponding unlabeled variant, and accordingly, under the experimental set-up, most of the proteins would increase in the label through the time of the experiment, but at different rates. In LC-MS, the partition of the peptide ion current between the heavy and light pairs provided the rate at which that protein was newly synthesized. The Mascot output files were combined with time-course information, metabolic functions and other relevant information using Microsoft Office Access 2007 and then loaded into Data Analysis Tool Extension (DAnTE) version 1.2 [59] for clustering and downstream analysis. DAnTE uses abundance (expression ratios) data or the spectral count data, and allows users to apply many analysis algorithms for data normalization, peptides rollup into proteins, data plotting and statistical analysis. The RIA values of individual peptides were log2 transformed, normalized to correct variations in labelling intensities of peptides, rolled up into proteins (Reference peptide based roll up; RRollup) and performed hierarchical clustering and visualized as heat maps. Labelled proteins were excluded from the list if they were identified by a single peptide pair and observed only at one time-point out of the 13 time points. But, proteins identified by a single peptide pair were accepted if that peptide was unique to that protein and identified in at least 2 out of the 13 time points. Proteins identified only at one time point were accepted if they contain ≥ 2 labelled peptides pairs. Proteins identified by shared peptides were accepted if they included at least one unique peptide pair.

A robust approach to determine both cycling Cyanothece transcripts and cycling proteins from near identical experiments have been previously reported [12, 16], which was used in the current study for a baseline comparison of the current results.