Temporal Dynamics of the Saccharopolyspora erythraea Phosphoproteome*

Actinomycetes undergo a dramatic reorganization of metabolic and cellular machinery during a brief period of growth arrest (“metabolic switch”) preceding mycelia differentiation and the onset of secondary metabolite biosynthesis. This study explores the role of phosphorylation in coordinating the metabolic switch in the industrial actinomycete Saccharopolyspora erythraea. A total of 109 phosphopeptides from 88 proteins were detected across a 150-h fermentation using open-profile two-dimensional LC-MS proteomics and TiO2 enrichment. Quantitative analysis of the phosphopeptides and their unphosphorylated cognates was possible for 20 pairs that also displayed constant total protein expression. Enzymes from central carbon metabolism such as putative acetyl-coenzyme A carboxylase, isocitrate lyase, and 2-oxoglutarate dehydrogenase changed dramatically in the degree of phosphorylation during the stationary phase, suggesting metabolic rearrangement for the reutilization of substrates and the production of polyketide precursors. In addition, an enzyme involved in cellular response to environmental stress, trypsin-like serine protease (SACE_6340/NC_009142_6216), decreased in phosphorylation during the growth arrest stage. More important, enzymes related to the regulation of protein synthesis underwent rapid phosphorylation changes during this stage. Whereas the degree of phosphorylation of ribonuclease Rne/Rng (SACE_1406/NC_009142_1388) increased during the metabolic switch, that of two ribosomal proteins, S6 (SACE_7351/NC_009142_7233) and S32 (SACE_6101/NC_009142_5981), dramatically decreased during this stage of the fermentation, supporting the hypothesis that ribosome subpopulations differentially regulate translation before and after the metabolic switch. Overall, we show the great potential of phosphoproteomic studies to explain microbial physiology and specifically provide evidence of dynamic protein phosphorylation events across the developmental cycle of actinomycetes.

Actinomycetes undergo a dramatic reorganization of metabolic and cellular machinery during a brief period of growth arrest ("metabolic switch") preceding mycelia differentiation and the onset of secondary metabolite biosynthesis. This study explores the role of phosphorylation in coordinating the metabolic switch in the industrial actinomycete Saccharopolyspora erythraea. A total of 109 phosphopeptides from 88 proteins were detected across a 150-h fermentation using open-profile two-dimensional LC-MS proteomics and TiO 2 enrichment. Quantitative analysis of the phosphopeptides and their unphosphorylated cognates was possible for 20 pairs that also displayed constant total protein expression. Enzymes from central carbon metabolism such as putative acetyl-coenzyme A carboxylase, isocitrate lyase, and 2-oxoglutarate dehydrogenase changed dramatically in the degree of phosphorylation during the stationary phase, suggesting metabolic rearrangement for the reutilization of substrates and the production of polyketide precursors. In addition, an enzyme involved in cellular response to environmental stress, trypsin-like serine protease (SACE_6340/NC_009142_6216), decreased in phosphorylation during the growth arrest stage. More important, enzymes related to the regulation of protein synthesis underwent rapid phosphorylation changes during this stage. Whereas the degree of phosphorylation of ribonuclease Rne/Rng (SACE_1406/NC_009142_1388) increased during the metabolic switch, that of two ribosomal proteins, S6 (SACE_7351/NC_009142_7233) and S32 (SACE_ 6101/NC_009142_5981), dramatically decreased during this stage of the fermentation, supporting the hypothesis that ribosome subpopulations differentially regulate translation before and after the metabolic switch. Overall, we show the great potential of phosphoproteomic studies to explain microbial physiology and specifically provide evidence of dynamic protein phosphorylation events across the developmental cycle of actinomycetes. Microorganisms have evolved mechanisms that enable them to grow and rapidly adapt to changing environmental conditions. Regulation of protein activity can occur at transcriptional, translational, and/or post-translational levels. Transcriptional and translational control are slow and have high energy costs due to de novo synthesis of proteins (i.e. transcription, translation, and protein-folding processes). Conversely, protein post-translational modifications drive adaptive cellular responses more efficiently by adding or removing functional groups from specific protein residues (1). Among the posttranslational modifications that regulate protein functionality, phosphorylation is by far the most studied in bacteria (2,3).
Two-component systems involve the phosphorylation of histidine and aspartate residues and were the first studied bacterial signal transduction mechanisms (3). Pioneering studies in Escherichia coli and Bacillus subtilis demonstrated extensive serine, threonine, and tyrosine phosphorylation, predominantly at sites with eukaryotic-like phosphorylation signatures (4,5). Since then, numerous studies have extended the repertoire of serine, threonine, and tyrosine kinases and eukaryotic-like phosphorylated proteins present in different bacteria (6 -15). Apart from some studies on specific enzymes of interest in E. coli and other model organisms, these studies have focused on mapping phosphorylation sites rather than identifying the biological role of phosphorylation. In order for phosphorylation to play a role in adaptive responses, it must display quantitative and dynamic variation. Recently, in vivo phosphorylation controlling enzyme functionality was studied in E. coli and Streptomyces coelicolor developmental cycles (16,17), showing the yet poorly explored effect of dynamic protein phosphorylation in microbial physiology.
Actinomycetes produce a large variety of secondary metabolites, including approximately half the antibiotics in current use (18). The production of secondary metabolites occurs after a critical culture transition known as the "metabolic switch" (19), believed to be triggered by nutrient limitation or oxidative stress. Understanding the regulatory mechanisms underpinning the reorganization of the metabolic and cellular machinery at the metabolic switch (20) is of both fundamental and practical importance, as efficient induction is essential for high-level production. Although some cellular regulatory mechanisms have been explained with regard to transcription (19,20), in vivo protein phosphorylation has not been yet explored at the global level and promises to fill an important gap in the understanding of the biology of actinomycetes.
Saccharopolyspora erythraea is a soil-dwelling actinomycete from the Pseudonocaridaceae family. This soil bacterium contains within its 8.29-Mb genome the machinery required for the synthesis of more than 25 different secondary metabolites, including erythromycin, the first clinically used macrolide antibiotic (21). Although highly exploited in industry, the S. erythraea secondary metabolism remains mostly unexplored; in fact, more than 17 secondary metabolites produced by this bacterium have unknown function and chemical structure (22). In addition, even though the S. erythraea genome was completed more than half a decade ago, industrial titers of erythromycin are obtained mostly via classical methods of random mutagenesis and fermentation media optimization using complex carbon and nitrogen sources (23). Several genomic and transcriptomic studies have compared genome sequences and gene transcription between wild-type and industrial erythromycin overproducing strains (24 -26). These investigations show that regulation of the erythromycin gene cluster is complex and may be regulated at the post-translational level.
Here, we present a dynamic phosphoproteomic study of the erythromycin-producing actinomycete S. erythraea. Using samples taken across a fermentation time course, a total of 109 phosphorylation sites were identified in discovery mode, before the degree of phosphorylation at each site was monitored using multiple reaction monitoring (MRM). 1 Quantitatively significant changes in phosphorylation were observed for many proteins during the metabolic switch and the stationary phase, and several of these events can be directly linked to known metabolic effects. We thus present a timeresolved dynamic study of protein phosphorylation in S. erythraea that specifically provides new insights into the physiology of actinomycetes.

EXPERIMENTAL PROCEDURES
Strain and Culture Conditions-S. erythraea strain NRRL23338 was purchased from the American Type Culture Collection (ATCC number 11635 TM ). Unless otherwise specified, all chemicals were purchased from Sigma. Medium ISP 2 (yeast extract, 4 g/l; malt extract, 10 g/l; dextrose, 4 g/l; agar, 20 g/l) was used for spore germination and seed cultures. Medium MM-101 used in the bioreactors contained (per liter) 7 g of NH 4  Samples were extracted from two different fermentations in 2-l Applikon reactors (Applikon Biosciences, Schiedam, The Netherlands) operated at working volumes of 1.4 l. The temperature and pH remained constant at 30°C and 7, respectively, throughout the fermentation; the pH was controlled by the addition of 20% sodium hydroxide (NaOH) or 10.9% hydrochloric acid (HCl). Dissolved oxy-gen was maintained at 45% to 60% saturation by increasing the air flow and reactor mixing. Carbon dioxide production was measured using an HPR20 QIC mass spectrometer (Hiden Analytical Ltd., Warrington, UK) attached to the bioreactor's condensers.
Protein Extraction and Trypsin Digestion-To characterize and quantitate the phosphoproteome of S. erythraea, we used monitoringinitiated detection and sequencing with information-dependant acquisition. Samples were extracted at six time points across the fermentation: one time point during the exponential phase (34 h), three around the metabolic switch (48,49, and 51 h), and finally two during the stationary phase (78 and 98 h). Proteins were extracted according to Ref. 4, incorporating some modifications. Cell pellets (100 ml of culture broth) were harvested from the bioreactor, pelleted (4°C at 5000 rpm) using an Allegra X-15R centrifuge (Beckman Coulter), washed in 100 mM NaCl, 25 mM Tris-HCl (pH 7.5), and resuspended in lysis buffer (50 mM Tris-HCl, 5 mg/ml lysozyme, and 5 mM of each of the following phosphatase inhibitors: sodium fluoride, 2-glycerol phosphate, sodium vanadate, and sodium pyrophosphate). After 10 min of incubation at room temperature, N-octylglucoside was added at a final concentration of 1% for the solubilization of membrane proteins. Enhanced cell disruption was achieved by homogenizing cell lysates with glass beads for 5 min at 4800 rpm using a Mini-bead beater (Extech Equipment Pty Ltd, Wantirna South, Australia). Cellular debris was removed by centrifuging (10 min, 4°C, 13,000 rpm) using a Microfuge (22R, Beckman Coulter). To remove nucleic acids, samples were incubated with DNase I (100 g/ml; Fermentas, Vilnius, Lithuania) and RNase A (100 g/ml; Fermentas) for 10 min at 37°C. Finally, protein extracts were dialyzed for 16 h in deionized water using 3.5 molecular weight cutoff Side-A-Lyser dialysis cassettes (Thermo Scientific). The total protein content was quantified using a 2D Quant assay kit (GE Healthcare). Protein extract aliquots (10, 5, and 1 mg and 500, 200, 100, and 50 g) were freeze-dried using Alpha 1-4 LSC (John Morris Scientific Pty Ltd, Willoughby, Australia) for 16 h and stored at Ϫ20°C. Trypsin digestion was performed as described in Ref. 4. Briefly, 10 mg of dried protein was resuspended in 6 M urea, 2 M thiourea, and 2% CHAPS for protein denaturation and incubated for 45 min at room temperature with 1 mM DTT. Iodoacetamide (2.5 mM final concentration) was added, and samples were incubated for 45 min at room temperature in the dark. Samples were then diluted with 25 mM ammonium bicarbonate to obtain a final concentration of urea of Ͻ800 mM. Trypsin (Promega Gold, MS grade) was added at a ratio of 1/100 (trypsin/protein), and samples were incubated for 16 h at 37°C. Reversed-phase chromatography was used to clean the digested peptide mixtures using C-18 cartridges (Sep-Pak tC18, Waters, Milford, MA) eluting with 60% acetonitrile (v/v). Residual acetonitrile was removed by vacuum centrifugation (Eppendorf, Hamburg, Germany) and resuspended in 0.1% formic acid prior to analysis. All samples were processed in parallel.
Phosphopeptide Identification-Peptide identification was performed in an LC MS/MS QSTAR Elite (AB Sciex, Ontario, Canada).
The LC system was equipped with a Vydac MS C18 300-Å, 150 mm ϫ 0.3 mm column (Grace Davison Discovery Sciences, Deerfield, IL) operated at 30°C with a 0%-80% acetonitrile gradient (in 0.1% formic acid) for 105 min at a flow rate of 3 l/min. All MS/MS raw data are available online through the University of Queensland website (S. erythraea phosphoproteomic data (385 MB)). Proteins were identified via advanced information-dependent acquisition of the fragmentation spectra of one to five charged peptides with a precursor selection window of m/z 100 -1800 using enhanced pulsed extraction of fragments (using Analyst 1.5.2; AB Sciex), employing specific features such as "Smart Collision" and "Smart Exit" (fragment intensity multiplier set to 2.0 and maximum accumulation time of 1.5 s) to obtain MS/MS spectra. Tandem mass spectra were acquired for 1 s, and fragmented peptides were selected for sequencing for 12 s in positive mode. The Paragon search algorithm (27) from Protein Pilot 4.0 software (Applied Biosystems, Foster City, CA) was used to identify all phosphoproteins. The mass tolerance values for precursor ions and fragment ions were set to the default values of the Paragon search algorithm. Trypsin was specified as the digesting protease, the iodoacetamide derivative of cysteine (carboxyamidomethylcysteine) was specified as the fixed modification, and urea denaturation and phosphorylation were specified as special factors. The sequence database used was taken from the S. erythraea Genome Project website (release version 15/03/2007). A false discovery rate analysis was performed for all searches. Hits were considered positive when at least two peptides with more than six residues and 90% confidence were detected. A monitoring-initiated detection and sequencing workflow was used for the relative quantification of peptide phosphorylation (28).
MRM Development-MRM assays for all the identified phosphopeptides and their unphosphorylated counterparts were developed and optimized from TiO 2 -enriched samples with MRMPilot TM 2.0 software (AB Sciex) according to the software manual. All quantitative experiments were performed from total protein extracts (without TiO 2 enrichment or chromatographic fractionation) using a triple quadrupole mass spectrometer (QTRAP TM 4000, Applied Biosystems) with an electrospray ion source configured in positive mode. LC was performed using a 100 mm ϫ 2.1 mm 2.6-m, 100-Å Kinetix C18 column (Phenomenex, Torrance, CA) running a gradient of 2%-80% acetonitrile (in 0.1% formic acid) for 100 min. Mobile phase A was an aqueous solution of 2% acetonitrile and 0.1% formic acid, and solvent B was 2% MilliQ water and 0.1% formic acid in 100% acetonitrile. Every peptide and phosphorylated cognate was monitored by at least two transition ions overlapping in retention time in the LC chromatogram (supplemental Table S7). The MS scan was performed with the following parameters: ion source voltage, 5400 V; temperature, 350°C; curtain gas, 20 psi; collisionally activated dissociation gas, high; TurboIonSpray nebulizer gas or Atmospheric-pressure chemical ionization nebulizer gas (Gas 1), 60 psi; TurboIonSpray heater gas (Gas 2), 60 psi; orifice (differential pressure), 80. To minimize technical variations, all injections were performed in triplicate.
Data Processing-NetPhos (29) and NetPhos-Bac (30) were used as primary sources for identifying phosphorylation sites. PhosCalc (31) was used to assign probabilities to the potential phosphorylation sites on the identified peptides. All possible phosphorylation positions, namely, tyrosine, threonine, serine, histidine, and aspartic acid, were evaluated. DAVID Bioinformatics Resources 6.7 (32) was used to perform gene-enrichment and functional annotation analysis. For all MRM experiments, automatic quantification methods were built in Analyst 1.5.2 (AB Sciex). All peaks were manually verified by visualizing their shape and correct elution time. This procedure was performed by two different people for two biological replicates in order to avoid bias in the analysis. The abundance for each peptide and its phosphorylated cognate was estimated as the average peak area of three technical replicates. The phosphorylation ratio was estimated from the area of the phosphorylated and unphosphorylated peptides. The phosphorylation ratios were log transformed to satisfy the linear model assumption of residual normality and to stabilize variances (phosphorylation ratio range of Ϫ3 to ϩ3). The total protein amount was estimated from the logged sum of peak areas for the phosphorylated and nonphosphorylated peptide cognates. Cluster 3.0 was used for all hierarchical cluster analyses. Complete linkage clustering was used to compute distances between groups of phosphoproteins. Dendrograms were made using Java TreeView 1.1.5r2 (33). Similarities between biological replicates were evaluated according to the Spearman rank correlation coefficient using a coefficient cutoff of 0.7 and p values Ͻ 0.2. The statistical significance (p value) for the correlation coefficients was estimated using R.
Analytical Procedures-Cell growth was monitored by measuring A 450 and by determining the cell dry weight during the fermentation. The cell dry weight was quantified by filtering 15 ml of broth using 22-m nylon membranes (Millipore, Bedford, MA) and drying the samples at 60°C for 8 h. The concentration of glucose was determined via high-performance liquid chromatography (Agilent 1200 HPLC system) (as described in Ref. 34) using a Phenomenex Rezex RHM Monosaccharide H ϩ column (300 ϫ 7.8 mm; 8 m) and detected by Refractive index. Erythromycin production was determined via LC-MS using a Dionex UltiMate 3000 liquid chromatography system (Dionex, Sunnyvale, CA) coupled with an AB Sciex 4000 QTRAP mass spectrometer (AB Sciex, Ontario, Canada) as described in Ref. 35.

S. erythraea Displays Biphasic Growth in a Bioreactor-S.
erythraea fermentation displayed a two-stage growth curve common to secondary metabolite producers. For the first 48 h, biomass was produced exponentially at a constant rate of 0.11 h Ϫ1 . A transient 30-to 60-min growth interruption (known as the metabolic switch) was observed around the 48th hour of culture. We monitored the metabolic switch online by connecting a mass spectrometer to the outlet of the bioreactor condenser. A drop of CO 2 production was clearly observed in the middle of the exponential phase, which is an indicator of growth arrest (Fig. 1). After the metabolic switch, the culture recommenced growth at the initial rate for ϳ10 h before entering the stationary phase. The main antibiotic produced by S. erythraea, erythromycin, was produced in a growth-independent manner for the first 90 h at a rate of 1.01 mg/g dry weight per hour. As the culture entered stationary phase, a significant increase in erythromycin production/release was observed (Fig. 1).
The S. erythraea Phosphoproteome: General Features-Openprofile two-dimensional LC-MS across six culture time points (see "Experimental Procedures") identified a total of 109 phosphopeptides from 88 proteins ( Fig. 2 and supplemental Table S1). Genes encoding the phosphorylated peptides were evenly distributed across the core region (53%), which encodes essential genes and primary metabolic functions, and the noncore region of the chromosome (47%), which codes for secondary metabolic functions ( Fig. 2A). Enrichment scores (ESs) were estimated to systematically assess whether some subset of the phosphoproteins showed significant overrepresentation of certain cellular functions (the greater the ES, the more represented a certain cellular function is in the dataset). Functional annotation analysis revealed that the S. erythraea phosphopro- Erythromycin production (red), dry cell weight (gDCW) (blue), and CO 2 (black) were measured during the fermentation. The fermentation is characterized by a growth arrest stage known as the metabolic switch. To characterize post-transcriptional modifications involved in the metabolic switch, samples were harvested from the bioreactor at several time points, indicated in the graph with black arrows.  (Fig. 2A).
Only half of the phosphorylation sites identified were detected by NetPhos (29) or NetPhos-Bac (30) (Fig. 2B). This low identification rate is likely because these prediction tools were constructed based on model organisms such as E. coli, B. subtilis, and Saccharomyces cerevisiae. Using PhosCalc (31), 60% of the phosphosites were assigned to a single residue (supplemental Table S3). The most phosphorylated residue was serine (26.5%), followed by threonine (23%), histidine (5.3%), and tyrosine (4.4%) (Fig. 2B). Annotated spectra of all detected phosphopeptides and their corresponding NetPhos and NetPhos-Bac scores are presented in supplementary File S1.
The S. erythraea Phosphoproteome Displays Dynamic Behavior across the Fermentation Time Course-We estimated phosphorylation dynamics by developing MRM assays for all the identified phosphopeptides and their unphosphorylated cognates. To this end, data from discovery proteomics were used to build an MS/MS library of spectra and facilitate the design process. Discovery data were loaded into MRMPilot and used to design and select the best fragment ions. In order to rectify the identity of the analytes, all MRMs were verified by ensuring complete co-elution/overlap of at least four fragment ions for each analyte (supplementary File S2). Furthermore, we verified the quantifiable signal-to-noise ratio in total protein extracts sampled across six time points of the fermen-tation. Out of 109 phosphopeptides discovered from enriched samples, we restricted our analysis to 60 phosphopeptides (supplemental Table S4) that satisfied the following criteria: (i) at least two independent consistent MRMs each for the phosphorylated and unphosphorylated peptides, (ii) consistency in retention time, (iii) identified in both biological replicates, and (iv) signal intensity based on the limit of quantification and signal-to-noise ratio.
The degree of phosphorylation was approximated based on the ratio of the signal intensities of the phosphopeptide ion and its unmodified cognate, as previously done to assess protein functionality regulated by phosphorylation in yeast central carbon metabolism (36). Reproducible phosphorylation events among two independent experiments were selected by estimating the Spearman rank correlation coefficients (see "Experimental Procedures"). A total of 44 phosphopeptides with correlation coefficients of Ͼ0.5 were clustered and classified according to their degree of phosphorylation and functional annotation (Fig. 3).
Profile analysis identified seven distinct phosphorylation clusters (supplemental Fig. S1 and supplemental Table S5). Ignoring positional information, these clusters can be mapped into three general profiles (Fig. 3). Profile A includes 15 phosphopeptides with constant degrees of phosphorylation across the fermentation (see clusters III and V in the supplementary figures). Among this group we found proteins related to general cellular processes, such as elongation factor TufA (SACE_6838/NC_009142_6720) and cell division control protein Cdc48 (SACE_0519/NC_009142_05191), and enzymes related to metabolic functions, such as peptide ABC transporter (SACE_6324/NC_009142_6201) and type I polyketide synthase (SACE_4140/NC_009142_4094).
Proteins with an increased degree of phosphorylation during the metabolic switch and stationary phase were grouped in Profile B (see clusters I, IV, and VII in the supplementary figures). These proteins are presumably involved in the activation of cellular processes during the metabolic switch. Among this group, we found a ribonuclease of the Rne/Rng family (SACE_1406/NC_009142_1388) and ArsR-family transcriptional regulator (SACE_4355/NC_009142_4304). Central carbon metabolic enzymes found with this profile (isocitrate lyase (SACE_1449/NC_009142_1431), adenosine kinase (SACE_ 3899/NC_009142_3856), pyruvate carboxylase (SACE_6118/ NC_009142_5998), 2-oxoglutarate dehydrogenase (SACE_6385/ NC_009142_6262), and phosphoglycerate mutase (SACE_ 3448/NC_009142_3404)) suggest that repression of anabolic enzymes and activation of systems for nutrient recycling might occur during this stage of fermentation. Additionally, this profile also includes four hypothetical proteins located in the core region of the chromosome (SACE_3064/NC_009142_3026, SACE_3268/ NC_009142_3226, SACE_3282/NC_009142_3240, and SACE_ 4206/NC_009142_4158) and three poorly characterized integral membrane proteins (SACE_6667/NC_009142_6544, SACE_ 6793/NC_009142_6674, and SACE_7119/NC_009142_7002).
Lastly, Profile C includes proteins with decreased phosphorylation during the metabolic switch and stationary phase (see clusters II and VI in the supplementary figures). Three enzymes with metabolic functions-acetyl/propionyl-CoA carboxylase (SACE_4237/NC_009142_4188), modular polyketide synthase (SACE_2875/NC_009142_2839), and FMNH 2 -dependent monooxygenase (SACE_4061/NC_009142_ 4018)-were found within this group. We also found a trypsinlike serine protease (SACE_6340/NC_009142_6216), a protein involved in cell dismantling and oxidative stress mechanisms in related actinomycetes (37,38). Interestingly, ribosomal proteins S32 and S6 (SACE_6101/NC_009142_5981 and SACE_7023/ NC_009142_6904, respectively) showed similar dephosphorylation behavior during the growth arrest stage. Lastly, SACE_ 6254/NC_009142_6132 and SACE_7023/NC_009142_6904, two hypothetical proteins located in the non-core region of the chromosome, were found in this group of proteins.
The changes observed in the phosphorylation ratios in Profiles B and C might reflect kinase-or phosphatase-mediated regulation events controlling metabolic and developmental pro-cesses. However, an altered ratio of protein phosphorylation may also result from dilution due to increased expression or preferential degradation of phospho-or nonphosphorylated proteins. The logged sum of peak areas for the phosphorylated and nonphosphorylated peptide cognates provides an indication of the total protein amount. For 20 out of 29 peptides in Profiles B and C, the difference across the profile was less than 2% (Fig. 3). Given that biological interpretation is more difficult when the total amount of protein changes across the developmental cycle, we limited our discussion to phosphoproteins that are differentially phosphorylated with less than 2% variation in the total amount of protein (Fig. 3, blue shaded boxes). DISCUSSION Actinomycetes undergo a dramatic reorganization of metabolic and cellular machinery during the metabolic switch leading to mycelia differentiation and the onset of secondary metabolite biosynthesis (38 -40). We speculated that protein phosphorylation might play an important role in the regulation of this process. The S. erythraea phosphoproteome reported  3. Selection criteria employed for the identification of proteins potentially regulated by phosphorylation. The signal-to-noise ratio was first verified independently for two biological replicates. A total of 60 quantifiable phosphopeptides were obtained. The Spearman rank correlation coefficient was estimated using a cutoff threshold of 0.5 (horizontal axis). The resulting 44 phosphopeptides were clustered according to their phosphorylation ratio profiles and classified into three function categories ("Metabolism," "Signaling and Cellular Processes," and "Hypothetical Proteins") according to the genome annotation. A total of 29 phosphopeptides showed dynamic behavior (table, second and third rows). The total amount of peptide quantified was estimated by adding the peak area of peptide ϩ phosphopeptide across the time course, and a maximum average variation of 2% was used as a threshold (blue-shaded boxes). For simplicity, the recent genome annotation has been omitted from the figure. here consists of 109 phosphopeptides from 88 proteins (Fig.  2), which is similar in size to most phosphoproteomes reported for other bacteria (6 -14). Sixty percent of phosphorylation events could be reliably assigned to a single site, and, as in other phosphoproteomes, serine and threonine accounted for the vast majority of events and were roughly equal in frequency (47% on serine and 45% on threonine) (41). Only half the phosphorylation events could have been predicted with NetPhos or NetPhos-Bac software.
Out of the 109 phosphopeptides, quantitative-dynamic MRM analysis was achieved for 60 phosphopeptides and their cognate non-phosphopeptides across six time points of the fermentation (Fig. 1). The failure to detect the remaining 49 peptides may be attributed to the absence of TiO 2 enrichment for paired quantification. A total of 44 peptides showed reproducible phosphorylation profiles across the two experiments (Spearman rank correlation coefficient Ͼ 0.5).
Clustering analysis revealed three phosphorylation profiles (Fig. 3). Fifteen peptides showed constant phosphorylation throughout the fermentation (Profile A) and thus are unlikely to be involved in regulation. Of the 29 peptides displaying dynamic phosphorylation (Profiles B and C), 20 also had near-constant total protein contents, indicating kinase/phosphatase activity potentially regulating the developmental cycle. The phosphoproteins were functionally classified according to the genome annotation in three groups, namely, unknown and hypothetical proteins, metabolism, and cellular processes/signaling.

Hypothetical proteins
SACE_3064-hypothetical protein SACE_6667-surface protein SACE_6793-integral membrane protein SACE_3268-hypothetical protein SACE_3282-selenium binding protein SACE_7023-hypothetical protein SACE_6254-hypothetical portein SACE_6793/NC_009142_6674 as the only protein with a homolog among related actinomycetes. SACE_6793/NC_ 009142_6674 putatively encodes an integral membrane protein part of a secretion system highly conserved among actinomycetes (ESX-1) (42). In S. coelicolor, knock-out of ESX-1 generates deficient coordination of cell division with segregation of nucleoids (43). Although the functionality of the ESX-1 system depends on a seven-amino-acid peptide signal in Mycobacterium tuberculosis (44), it remains unknown how this system is targeted for secretion in S. coelicolor. It is unclear what role increased phosphorylation over the stationary phase may play, but its presence suggests that SACE_6793/NC_009142_6674 could be involved in the regulation of this important secretion system.
Metabolic Processes-Eight proteins related to either anabolic or catabolic pathways were detected with dynamic phosphorylation profiles. A putative acetyl-CoA carboxylase (SACE_4237/NC_009142_4188) was found to be phosphorylated only during the metabolic switch (Fig. 5). Acetyl-CoA carboxylases catalyze the conversion of acetyl-CoA to malonyl-CoA in the first step of the biosynthesis of fatty acids and the formation of polyketide precursors (45). Isocitrate lyase (ICL) and the E1 subunit of 2-oxoglutarate dehydrogenase (ODH-E1) had opposite phosphorylation profiles (Fig. 6). This result correlates with the natural competition for carbon flux between the glyoxylate shunt and the TCA cycle. Sequence analyses demonstrate that isocitrate lyase is highly conserved; however, their regulation seems to differ significantly. In E. coli, decreased isocitrate lyase activity is observed when the enzyme is dephosphorylated (46), whereas the opposite effect has been demonstrated in yeast (47). 2-Oxoglutarate dehydrogenase was dephosphorylated as cells entered the stationary phase, but the role of this phosphorylation event remains unknown. Phosphorylation of 2-oxoglutarate dehydrogenase has been identified in other bacterial phosphoproteomes (11). Beyond their well-known metabolic functions, both enzymes are involved in the biosynthesis of TCA intermediates, specifically succinate and succinyl-CoA. We consider this relevant, as succinyl-CoA is the precursor of the polyketide backbone molecule methyl malonyl-CoA (48). Several metabolic engineering studies have identified the availability of acyl-CoA derivatives as the main factor for improving polyketide production (49,50). Correlating these phosphorylation events during the stationary phase (when secondary metabolism is active) suggests a regulatory mechanism that facilitates the biosynthesis of methyl malonyl-CoA. This result not only shows a link between enzyme phosphorylation and the transition between primary and secondary metabolism, but also provides new insights into the use of regulatory mechanisms as strategies to improve secondary metabolite producing strains.
Another enzyme found with an increasing phosphorylation ratio during the late stage of the culture was adenosine kinase (SACE_3899/NC_009142_3856), which catalyzes the biosynthesis of adenosine monophosphate from adenosine. The phosphorylation site is not in the putative catalytic site, and increased phosphorylation may be linked to increased enzyme regulation. Adenosine kinase is essential for sporulation (51); thus the increase in phosphorylation as cells enter the stationary phase may reflect increased enzyme activity. Finally, several enzymes involved in sulfur and redox metabolism were identified: adenylylsulfate kinase (SACE_1473/ NC_009142_1456), cysteine desulfurase (SACE_2177/NC_ 009142_2150), FAD-linked oxidoreductase (SACE_0195/ NC_009142_0193), and FMNH 2 -dependent monooxygenase (SACE_4061/NC_009142_4018).
Signaling and Global Cellular Processes-Five proteins involved in global cellular processes with dynamic phosphorylation profiles were identified (Fig. 6). The first protein, SACE_4355/NC_009142_4304, encodes for a transcriptional regulator of the ArsR family. The phosphorylation degree of this regulator was constant at the early stages of the culture and increased during the metabolic switch and the stationary phase. Although this class of transcriptional regulators is conserved among actinomycetes, their function has not been deeply investigated. In several studies, this regulator has been related to the cellular response to hypoxic adaptation and the regulation of the phoP regulon (52,53).
Ribonuclease E/G (SACE_1406/NC_009142_1388) phosphorylation dramatically increased during the metabolic switch (Fig. 6). Ribonucleases (RNases) of the Rne/Rng family are specific endonucleases involved in processing 16S and FIG. 6. Two-dimensional plots of proteins related to cellular processes displaying dynamic phosphorylation profiles across the fermentation. The phosphorylation ratio (log P-ratio) and total peptide (log peak area (PA)) are plotted as solid red and dotted blue lines, respectively. Graphs represent Rne/Rng family ribonuclease, SACE_1406 -2/NC_009142_1388 -2 (upper left panel); trypsin-like serine protease, SACE_6340/NC_009142_6216 (upper middle panel); ribosomal protein S6, SACE_7351/NC_009142_7233 (upper right panel); ArsR family transcriptional regulator, SACE_4355/NC_009142_4304 (lower left panel); and ribosomal protein S32, SACE_6101/NC_009142_5981 (lower middle panel). Total peptide measurements were estimated as log (peak area of phosphopeptide ϩ log peak area of peptide). For simplicity, the recent genome annotation has been omitted from the figure.
5S RNA precursors (54). Targeted RNA degradation has been suggested as a global post-transcriptional mechanism for regulating the developmental cycle and antibiotic production in S. erythraea (20) and S. coelicolor (55). RNase catalysis regulated by phosphorylation has been reported in eukaryotic systems (56) and E. coli for the RNase catalytic domain of the polynucleotide phosphorylase (57). This study provides evidence that dynamic phosphorylation of the ribonuclease E/G is correlated with the metabolic switch, suggesting cross-talk between post-transcriptional and post-translational regulatory mechanisms.
Another protein found with dynamic phosphorylation across the developmental cycle was a trypsin-like serine protease (SACE_6340/NC_009142_6216). This enzyme, which was phosphorylated during early stages of the fermentation, rapidly decreased in phosphorylation degree at the metabolic switch (Fig. 6). Sequence alignment showed that this enzyme is homologous to high temperature requirement A (HtrA) protease, widely conserved in single-celled and multicellular organisms (58). HtrA protease is essentially involved in ATPindependent protein quality control. For example, bacterial HtrA proteases are essential for the stress response to protein denaturation (58). In addition, HtrA functions as a nuclear serine protease essential for apoptosis in S. cerevisiae (59). Despite being widely characterized, phosphorylation-mediated activation has been reported only for mitochondrial HtrA2 (60). In this study, the phosphorylation profile observed for HtrA suggested that this enzyme could be regulated at the post-translational level, and this provides new insights into the regulation of protein degradation processes during the metabolic switch in actinomycetes.
Finally, two ribosomal proteins, S32 and S6 (SACE_6101/ NC_009142_5981 and SACE_7351/NC_009142_7233, respectively), were found to have notable phosphorylation dynamics. While constitutively transcribed, both proteins decreased dramatically in phosphorylation degree at the metabolic switch (Fig.  6). This result is consistent with the "ribosome filter hypothesis" (61), which establishes that the phosphorylation of some ribosomal proteins, particularly S6, generates variations in the stoichiometry of proteins associated with the ribosome. These ribosomal subpopulations selectively translate specific subsets of mRNAs, thus regulating translation (62). Evidence of ribosomal subpopulations in actinomycetes has been reported for S. erythraea (20) and S. coelicolor (63)(64)(65), in which rRNA and ribosomal protein operons are differentially expressed and undergo targeted degradation. The results presented here suggest that ribosomal subpopulations might exist in S. erythraea, regulating gene expression at the translational level.
In summary, this work represents the only in vivo dynamic phosphoproteomic study of S. erythraea to date and highlights the importance of studying phosphorylation dynamics to understand physiological processes in actinomycetes. In addition, the results presented provide new hypotheses for the regulation of central carbon metabolic enzymes and developmental programs during the metabolic switch and stationary phase.