Pseudomonas aeruginosa reverse diauxie is a multidimensional, optimized, resource utilization strategy

Pseudomonas aeruginosa is a globally-distributed bacterium often found in medical infections. The opportunistic pathogen uses a different, carbon catabolite repression (CCR) strategy than many, model microorganisms. It does not utilize a classic diauxie phenotype, nor does it follow common systems biology assumptions including preferential consumption of glucose with an ‘overflow’ metabolism. Despite these contradictions, P. aeruginosa is competitive in many, disparate environments underscoring knowledge gaps in microbial ecology and systems biology. Physiological, omics, and in silico analyses were used to quantify the P. aeruginosa CCR strategy known as ‘reverse diauxie’. An ecological basis of reverse diauxie was identified using a genome-scale, metabolic model interrogated with in vitro omics data. Reverse diauxie preference for lower energy, nonfermentable carbon sources, such as acetate or succinate over glucose, was predicted using a multidimensional strategy which minimized resource investment into central metabolism while completely oxidizing substrates. Application of a common, in silico optimization criterion, which maximizes growth rate, did not predict the reverse diauxie phenotypes. This study quantifies P. aeruginosa metabolic strategies foundational to its wide distribution and virulence including its potentially, mutualistic interactions with microorganisms found commonly in the environment and in medical infections.


Results
Growth physiology and substrate preference of rCCR . Pseudomonas aeruginosa strain 215 (Pa 215) is a medical isolate from a chronic wound 33,34 . Pa 215 was grown in chemically-defined, glucose containing, CSP G medium (materials and methods, supplementary material S1). Cultures exhibited two distinct exponential growth phases followed by stationary phase (Fig. 1). A subset of amino acids was consumed preferentially during the first exponential growth phase, which had the highest specific growth rate (Fig. 1, supplementary material  S2). The second exponential growth phase corresponded with the catabolism of second and third tier amino acids and glucose. CSP G media contained a small concentration (3 mM) of citrate, which was added as an ion chelator; however, the citrate was readily catabolized as a preferred substrate during the first growth phase. The cultures did not exhibit an overflow metabolism defined by the secretion of reduced metabolic byproducts like acetate, as is typical of microorganisms expressing cCCR phenotypes 17,30 . Trace amounts of gluconate were secreted during glucose metabolism but were quickly depleted (supplementary material S3). Amino acid deamination products like α-ketoglutarate and pyruvate were not observed in spent medium or found in only trace amounts (< 1 mm), respectively.
Substrate utilization order for Pa 215 was quantified using five different formulations of CSP G medium supplemented with permutations of additional carbon sources: lactate (L), acetate (A), and succinate (S) (Fig. 2). Pa 215 grown on CSP GL medium, preferentially consumed the top tier amino acids, represented by aspartate in Fig. 2, followed by lower tier amino acids (data for each measured amino acid can be found in supplementary material S3-S8) and lactate before finally catabolizing glucose. Glucose catabolism was not observed while lactate was present. Pa 215 grown on CSP GA consumed the top tier amino acids followed by lower tier amino acids and acetate and finally glucose after the acetate was exhausted. Pa 215 grown on CSP GLA preferentially consumed top tier amino acids, then lower tier amino acids and lactate followed by acetate and glucose. Finally, Pa 215 grown on CSP GLAS preferentially catabolized the top tier amino acids, followed by succinate, lactate, acetate, and ultimately glucose. Ion chelator, citrate, was readily catabolized as a preferred substrate in all media formulations. No or minimal overflow metabolism (< 4 mM acetate, ~ 3% of lactate and glucose carbon moles in CSP GL medium) was observed. An exception was CSP GLAS grown cultures which accumulated acetate (~ 10 mM) above the initial medium concentrations. Upon exhaustion of succinate and lactate, the acetate was catabolized prior to glucose catabolism. The glucose was not completely catabolized in CSP GLAS medium because the medium was nitrogen limited (supplementary material S1). Culture parameters are summarized in supplementary material S2 and data is available in supplementary material S3-S7.
The order of amino acid catabolism was assessed for all five CSP medium formulations by binning the amino acids into three categories based on their time of exhaustion (Table 1). Binning was used, as opposed to using an absolute time metric, because each medium formulation had a different number of substrates leading to different total growth times. Other substrate usage metrics were considered including the initial time of substrate catabolism and nonlinear fitting of the temporal metabolite profiles to calculate a substrate 'half-life' value 35 ; these metrics were sensitive to experimental variability during the initial growth phase where small fluctuations in substrate concentration, based likely on analytical techniques, influenced predictions. All amino acid data, The preferred amino acids, referred to here as top tier, included aspartate, asparagine, glutamine, glutamate, and alanine. Most top tier amino acids were binned consistently across the five CSP media formulations as quantitated by their small standard deviations. The middle tier amino acids had more variability which may have been CCR-related or based on the temporal granularity of the experimental sampling schedule. The temporal trends in medium pH reflected the metabolism of different substrates. Catabolism of amino acids increased medium pH based on nitrogen chemistry. Catabolism of organic acids also raised the culture pH because the bacterium imports the protonated base, removing a proton from the medium. Approximately 90% of the anabolic nitrogen in CSP G medium was in the form of amino acids (supplementary data S1). P. aeruginosa can use ammonium as the sole nitrogen source 36 . CSP G medium formulations were modified with the addition of 2 g/L ammonium chloride to test the effect of nitrogen form. The rCCR phenotype was not changed by the presence of ammonium. The cultures consumed the amino acids as preferred substrates followed by lactate and then glucose (supplementary material S9). The amino acid utilization order remained largely unchanged (supplementary material S10).
The common laboratory strain of P. aeruginosa, PAO1, was also grown on CSP GLAS medium. The PAO1 substrate utilization order of organic acids and glucose was the same as Pa 215 and the order of amino acid consumption was very similar to Pa 215 (supplementary material S11). Cultures demonstrated two exponential growth phases highlighted with different background shading. Amino acids were binned into three categories based on their time of exhaustion. Top tier amino acids were consumed during the first exponential growth phase while lower tier amino acids and glucose were consumed during the second exponential growth phase. All values are averages of three biological replicates, and metabolite values are also averaged from two technical replicates. Additional data can be found in supplementary material S3. www.nature.com/scientificreports/ Proteomics quantifies a constitutive, respiration-centric metabolism. Proteomic data were collected from CSP G and CSP GL grown cultures. Proteomic data are more predictive of cell function than transcriptomic or genomic data alone because they represent an actual allocation of resources into relatively stable, macromolecular pools 37,38 . Phenotypes were analyzed using label-free proteomics with mass-spectrometry (MS) of whole-cell lysates collected mid-first, exponential growth phase (4 h), early-second, exponential growth phase (7 h), and late-second, exponential growth phase (11 h). The proteomics data were analyzed with focus on central metabolism proteins associated with catabolizing the available substrates and with producing cellular energy. Enzymes from the tricarboxylic acid (TCA) cycle and associated auxiliary enzymes had largely, constitutive abundances regardless of the medium formulation and the growth phase (Fig. 3). All TCA cycle enzymes except the membrane-associated succinate dehydrogenase were detected and quantified. Additionally, the enzymes oxaloacetate decarboxylase (PA4872) and PEP synthase which process metabolic intermediates from the TCA cycle for gluconeogenesis, were expressed constitutively. The abundance of ATP synthase subunits was also constitutive. Membrane-associated, electron transport chain (ETC) enzymes were not detected. It was assumed that the ETC enzymes were also constitutively expressed based on the TCA cycle and the ATP synthase protein abundances and the lack of an overflow metabolism.
Enzymes associated with the processing of specific substrates did change in abundance based on presence and concentration of substrates, contrary to most TCA cycle enzymes (Fig. 4). Aspartate was plotted as a representative top tier amino acid (Table 1). Protein abundance for aspartate ammonia-lyse (AspA), responsible for the   36,39 . Abundance of ED phosphogluconate dehydratase (Edd) increased while glucose was being metabolized, after the exhaustion of top tier amino acids. Acetate kinase (AckA) abundance increased for the CSP GL culture at the exhaustion of lactate and glucose. AckA quickly metabolized the small amount of acetate (< 4 mM) secreted at the exhaustion of lactate and glucose. Proteomic analysis measured additional proteins that displayed changes in expression during exponential growth and stationary phases. Data can be found at ftp://massi ve.ucsd.edu/MSV00 00855 90/.
In silico analysis of rCCR phenotypes. CCR is a regulation scheme that contributes to metabolic plasticity. CCR regulation schemes have evolved to control expression of metabolic strategies that favor cellular fitness. The order of substrate utilization is hypothesized to reflect the ecological strategy used by P. aeruginosa to thrive in environmental and medical niches. Computational systems biology was used to test hypotheses regarding what fitness properties were being optimized in the Pa 215 cultures, with predictions compared to experimental data. In silico analyses used flux balance analysis (FBA) of a published, genome-scale, metabolic model of P. aeruginosa updated here with genome-supported, amino acid catabolism reactions 40,41 (supplementary material S12). Stoichiometric modeling methods, such as FBA, can be utilized with a minimum number of a priori fitting parameters. The applied FBA considered only steady state simulations. Temporally-resolved simulations require enzyme kinetic parameters for every considered substrate, which are not available in the literature for P. aeruginosa.
In silico testing of ecological strategies was applied first to amino acid utilization order and included all experimentally measured amino acids except for aromatic and sulfur containing amino acids due to their specialty chemistries. The experimental amino acid utilization order did not correlate with the amino acid frequency in genome open reading frames (Fig. 5a) indicating the amino acids were not consumed solely for protein assembly; amino acids were also used as anabolic building blocks for other macromolecules and catabolized for cellular energy (supplementary material S8). Therefore, simulations considered either the production of cellular energy (e.g. ATP) or cellular growth which was quantified as carbon moles (Cmol) of biomass.
The first round of in silico analyses considered six separate, single dimension, optimization criteria which were informed by previous studies that examined numerous optimization criteria 22,30,42 . The criteria included (1) maximizing biomass or energy production rates based on electron donor, (2) maximizing biomass or energy production rates based on electron acceptor (O 2 ), or (3) minimizing nutrient investment into the proteome required for either biomass or energy synthesis. The results of these simulations are presented in the next two sections and a summary of the results and analyses can be found in supplementary material S13.
Amino acid utilization order did not correlate with in silico maximization of rates. Computational approaches for studying metabolism often assume cells utilize metabolic potential to maximize growth rate [43][44][45] . The experimental amino acid utilization order, as quantified across five media formulations (Table 1), was analyzed using this maximum rate theory. Separate, steady state simulations were run for each individual amino acid and for either biomass production or cellular energy production. Simulations identified the optimal phenotypes for the conversion of each substrate into product. The in silico phenotypes which maximized elec- Table 1. Amino acid utilization order for P. aeruginosa 215 cultures grown on five different, chemicallydefined media supplemented with various additional carbon sources. Amino acids were binned into three categories (1, 2, 3) based on the time of exhaustion, averaged between three biological replicates for each of five medium conditions, n = 15. Data can be found in supplementary material S8. www.nature.com/scientificreports/ tron donor yields also mineralized the substrates secreting only CO 2 , consistent with in vitro cultures which did not utilize overflow metabolisms. Maximizing cellular energy yield per amino acid, on a mole substrate or Cmol substrate basis, did not correlate with amino acid utilization order, having r 2 values of 0.06 and 0.15 respectively (Fig. 5b). The in silico phenotypes which maximized yields were used to calculate maximum product rates using enzyme parameters from survey studies 46,47 and experimental medium composition (supplementary material S1). This common theoretical treatment linked product yields to product rates such that maximizing one maximizes the other, this was an assumption of convenience and the advantages and disadvantages of its application have been discussed in the literature 30,43,48 . Maximizing the rate of cellular energy production (or growth) did not predict amino acid utilization order, as the correlation was r 2 = 0.17 (Fig. 5c) (supplementary material S13-S18). Proteomics data for P. aeruginosa 215 grown on chemically-defined, CSP G and CSP GL medium which contained glucose (G) or glucose and lactate (GL), respectively. CSP G culture data are represented by dark blue bars and CSP GL culture data are represented by light blue bars. Bars quantify the abundance of the enzyme during the first exponential growth phase (4 h), early second exponential growth phase (7 h), and late second exponential growth phase (11 h). Presented enzymes are involved in the tricarboxylic acid (TCA) cycle, anaplerotic reactions, and ATP synthesis. All values are averaged from three biological replicates. *p value < 0.05; **p value < 0.01. www.nature.com/scientificreports/ P. aeruginosa has a respiration-centric metabolism. The rate maximization criterion was also applied to O 2 which is required to mineralize the amino acids under the experimental conditions. This alternative, rate maximization criterion predicted substantial overflow metabolisms for most of the amino acids. This predicted phenotypic trait was not consistent with the experimental data indicating the criterion was not relevant for Pa 215 metabolism (supplementary material S13).

Scientific Reports
Amino acids with high cellular energy yields (mol ATP (mol amino acid) −1 ) also had high biomass yields (Cmol biomass (mol amino acid) −1 ); the two yields correlated with an r 2 value of 0.99 (supplementary data S18). Therefore, maximizing rates for cellular energy production or biomass production had similar trends and neither predicted rCCR phenotype (supplementary material S13).  Analysis considered a base metabolism scenario where all enzyme-catalyzed reactions were considered and a refined, no core metabolism scenario where only enzymecatalyzed reactions extraneous to the experimentally-measured, constitutive proteome core where considered. A subset of the core enzymes is shown in Fig. 3 while an explicit list can be found in supplementary material S15.
(e) Proteome investment analysis using the minimal proteome investment proxy to predict experimental amino acid utilization order. Analysis considered base metabolism scenario and refined, no core metabolism scenario. www.nature.com/scientificreports/ resource investment into catabolic pathways. Explicit investment models were not possible in non-model organism, P. aeruginosa. Therefore, two previously developed resource investment proxies were applied to estimate relative, proteome investment into central metabolism 19,22,49,50 . The flux minimization proxy assumes the total network flux is proportional to the enzymatic resources required to synthesize the necessary proteome 19,22,42,51 .
Another proxy for protein investment minimizes the number of enzyme catalyzed reactions (a.k.a. minimal proteome investment) identifying the smallest proteome required to realize an in silico phenotype 19 . These hypotheses assumed central metabolism enzymes could be approximated as having the same molecular weight with the same amino acid distribution. Both proxies for resource investment, when applied as the single optimization criterion, predicted overflow metabolisms for most amino acids (supplementary material S13,S15,S16). The in vitro experimental cultures did not demonstrate substantial overflow metabolisms indicating these single dimension criteria were not relevant for Pa 215 phenotypes.
Substrate utilization order was consistent with a resource utilization strategy optimizing substrate oxidation and proteome investment. Life occurs in multifactorial environments with multiple stressors influencing phenotypes 20,42,52 . Two dimensional, optimizations of in silico phenotypes were performed where the first dimension considered the optimal conversion of substrate into cellular energy which completely oxidized the substrate. The second dimension approximated the nutrient investment into the enzymes required to realize the in silico phenotype, for example the amount of anabolic nitrogen required to synthesize the proteome or the amount of ATP required to form the associated peptide bonds 19,22 (supplementary material S15,S16). Both the flux minimization and minimal proteome investment proxies were tested. Two-dimensional optimization (2-DO) using the flux minimization proxy had poor correlations with the observed amino acid utilization order (Fig. 5d). Alternatively, 2-DO using complete substrate oxidation and the minimal proteome proxy predicted the experimental utilization order for amino acids (r 2 = 0.67) (Fig. 5e).
2-DO was further refined using experimentally measured proteomics data. The constitutively expressed TCA cycle, anaplerotic enzymes, ATP synthase, and electron transport chain (Fig. 3) were considered part of a core, constitutive proteome, independent of substrate. The refined, 2-DO theory considered only the resource investment extraneous to the conserved, core proteome. This theory lead to improved predictions of amino acid utilization order with the minimal proteome investment theory but not the flux minimization theory (Fig. 5d, e). The outlier amino acid in Fig. 5e was serine. Serine is catabolized via the L-serine dehydratase enzyme which is O 2 -labile suggesting higher cell densities and lower O 2 concentrations were necessary for its functionality 53 . The predictive accuracy of the analysis improved to a correlation of r 2 = 0.88 if serine data were excluded.
2-DO, considering complete oxidation of substrate and minimal proteome investment, was extended to the other CSP media substrates including organic acids and glucose. Analysis applied the minimal proteome investment with conserved core proteome assumption and considered both cellular energy production as well as the more complex biomass production (Fig. 6a). The experimental substrate utilization hierarchy, which was determined using culturing data from the five CSP media formulations, was used to assess the accuracy of the predictions (Fig. 2). The in silico analysis accurately predicted, substrate utilization order with r 2 correlations of 0.94 and 0.73 for cellular energy and biomass production, respectively. The cellular energy simulations had a noteworthy correlation with experimental data suggesting ATP production was a superior in silico optimization criterion for Pa 215. The biomass simulations considered an aggregate amino acid substrate pool containing all 20 metabolites, which was not considered for cellular energy simulations. As anticipated, the aggregate amino acid pool greatly reduced the requirement for enzymatic steps by negating de novo amino acid synthesis reactions (Fig. 6a, supplementary material S15,S16). The correlation between the predicted and experimental substrate utilization order for the biomass simulations was not as strong due largely to the predicted order of lactate and citrate utilization. The discrepancy could be due to a couple factors. First, the computational approach approximated the resource investment necessary to synthesize the in silico proteome by assuming all enzyme-catalyzed reactions required the same amount of anabolic resource. This was a necessary simplification due to the lack of detailed data for P. aeruginosa that could be improved as more data becomes available. Additionally, the experimental utilization order was based on five separate media formulations which resulted in five separate, dynamic, batch growth profiles each with their own intricacies including cometabolism of substrates.
The maximization of rate criterion was also tested with the additional substrates. The analysis assumed optimal product yields on substrate were proportional to the optimal product rates 43 . The maximum rate criterion did not predict the experimental utilization order for organic acids over glucose. In fact, the predicted utilization order had negative correlations with the experimental data (Fig. 6b). Additional optimizations and aggregate substrate simulations were considered (supplementary material S17-S19). None outperformed the presented approach in terms of accuracy and simplicity.

Discussion
P. aeruginosa preferentially consumes nonfermentable, lower energy substrates, such as succinate over glucose in a strategy known as reverse diauxie or rCCR. The term has been defined in terms of substrate preference relative to cCCR organisms E. coli and B. subtilis. The rCCR preference for nonfermentable substrates is associated with minimal overflow metabolism and, under certain conditions, can result in cultures preferentially catabolizing substrates that do not maximize cellular growth rates 10,36,[54][55][56][57] . The CSP media studied here did not result in this property. The term 'inverse diauxie' has been proposed to describe microorganisms that prefer substrates that sustain lower growth rates 58 . The rCCR strategy has enabled the broad, global distribution of P. aeruginosa in both environmental and medical niches including chronic, diabetic ulcers. The hierarchy of substrate preferences for Pa 215 was: amino acids such as aspartate, followed by citrate, succinate, lactate, acetate, and finally glucose. www.nature.com/scientificreports/ These preferences were also observed with P. aeruginosa PAOI grown on CSP GLAS medium (supplementary material S11). Pa 215 maintained, constitutively, core TCA cycle enzymes and regulated the abundance of the proteins required for specific substrates as needed to convert the substrates into central metabolism intermediates. Analysis using an in silico metabolic model and FBA determined the rCCR phenotype was consistent with a multidimensional, resource utilization strategy where substrate utilization order was based on minimizing the proteome investment required to mineralize the metabolite. Optimization of multiple cellular functions simultaneously has been reported previously for cCCR model organism E. coli, albeit with optimization criteria 22,52 . The rCCR phenotypes were not consistent with the commonly applied, systems biology criterion, which maximizes growth rate 43,44 . Pseudomonads are commonly found in consortia and the rCCR metabolism is proposed to provide fitness advantages in these competitive environments 59,60 . Consortia with populations expressing rCCR and cCCR phenotypes have the metabolic basis for an effective division of labor, thus avoiding overlapping substrate preferences which can result in species competition and exclusion [61][62][63][64][65] (Fig. 7). cCCR organisms prefer primary substrates like glucose where fast growth rate is likely a strong fitness determinant. Glucose can be catabolized by cCCR microorganisms via respiration, fermentation, or a combination of the two strategies. This flexibility enables tradeoffs between high yields during a fully respiratory catabolism and fast rates, with the associated overflow metabolism, during (partially) fermentative catabolism 30,48 . rCCR organisms avoid competition for  (5) glucose. (a) Substrate utilization order predictions for cellular energy production and biomass production using a two-dimensional optimization including complete substrate oxidation and minimal proteome investment. Simulations used a refined, core proteome theory where only enzyme-catalyzed reactions extraneous to the experimentally-measured, constitutive, core proteome where considered. A subset of the core enzymes is shown in Fig. 3 while the explicit list can be found in supplementary material S15. (b) Predicted substrate preference based on the 'maximization of rate' criterion for cellular energy production and biomass production. In silico product yield on substrate was assumed proportional to rate. www.nature.com/scientificreports/ fermentable carbon sources by preferring to catabolize secondary, byproducts of cCCR metabolism. Fast growth may not be as central to fitness as the efficient extraction of energy from the lower energy substrates. Many rCCR preferred organic acids are nonfermentable which precludes the rCCR metabolism from utilizing the tradeoffs inherent to either high yield or high rate strategies, as commonly observed in cCCR organisms 48 . Both rCCR and cCCR phenotypes are found in microorganisms described as generalists, e.g. P. aeruginosa and E. coli. Both rCCR and cCCR phenotypes can be predicted using resource investment theories, albeit with the rCCR organism expressing a respiration-centric phenotype and the cCCR organism expressing a glycolysiscentric, overflow phenotype. Natural environments are often limited by anabolic nutrients including nitrogen 18 . Division of labor, thru rCCR and cCCR based phenotypes, can theoretically enable higher consortia fluxes, from a scarce nitrogen supply, based on the nonlinear relationship between enzyme flux and resource investment 17,66 (Fig. 7c). This kinetic effect can translate into consortia having a better metabolic return on limiting nutrients, leading to higher biomass accumulation and higher host bioburden 17,[67][68][69][70] . Additionally, the CCR-based, division of labor could create a positive feedback mechanism by removing inhibitory organic acids and preventing environmental acidification via both organic acid consumption and amino acid catabolism, ultimately increasing consortia productivity by permitting a more complete depletion of substrates 71 . Amino acid catabolism and the release of ammonia can also function as an intercellular communication strategy where the small metabolite influences phenotype, like antibiotic susceptibility, in distant populations [72][73][74] .
Most virulence mechanisms are nutrient acquisition strategies that are also effective in medical niches 8 . CCR regulates a wide range of social behaviors and likely modulates division of labor which would facilitate substrate acquisition by rCCR microorganisms 17 . P. aeruginosa preference for non-fermentable substrates like succinate makes it a secondary resource specialist that requires terminal electron acceptors like O 2 or nitrate 75 . However, O 2 is often limiting in biofilms where cellular O 2 consumption rates are faster than diffusion rates 34,76 . P. aeruginosa possesses effective mechanisms to acquire scarce resources like O 2 70,77-80 . For example, P. aeruginosa secretes a cocktail of moieties such as pyocyanin, quinolones, and cyanide 59,[81][82][83][84] . Exposure to this cocktail can manipulate the S. aureus cCCR phenotype, driving it toward overflow and fermentative metabolisms 70,85 . Collectively, the compounds enable a secondary consumer to influence the metabolism of neighboring cells directing their phenotypes toward secreting preferred substrates including organic acids while reserving the O 2 for P. aeruginosa 69,70,84 .
Lactate has remarkable connections to P. aeruginosa substrate preference and medical niches including diabetic wounds. Elevated lactate levels found in diabetic wounds come from two sources. First, diabetic patients can have higher levels of serum lactate due to diabetic ketoacidosis, and secondly, lactate is associated with wound bed colonization by bacteria which produce it as a byproduct [86][87][88] . > 80% of chronic wounds are colonized by P. aeruginosa 89 while 90% of chronic leg ulcers are colonized by S. aureus 78,[90][91][92][93] which displays a cCCR phenotype 94,95 . Not surprisingly, these bacteria are often co-isolated 55,85,96 . Wounds colonized by multispecies can be more difficult to treat and can have more negative outcomes than wounds colonized by a single species 69,91,93 . Mutualistic interactions in consortia, based on complementary rCCR and cCCR metabolisms, could lead to www.nature.com/scientificreports/ emergent properties such as enhanced biomass productivity based on enhanced resource acquisition and better metabolic return on investment of scarce nutrients, ultimately leading to greater virulence. Mitigating these consortia, through rational countermeasures, will require quantitative knowledge of the metabolic organization which forms the bases of all virulence mechanisms.

Materials and methods
Bacterial strain and cultivation. All experiments used P. aeruginosa str. 215, a clinical isolate obtained from a chronic wound 33,97 . Frozen stocks of P. aeruginosa 215 were prepared by growing cultures in 10 mL of 1/10 strength tryptone soy broth (TSB) at 37 °C with shaking (150 rpm), mixed with 3 mL of 20% glycerol, and stored at − 80 °C. Frozen stocks were plated on tryptic soy agar (TSA) at 37 °C for 12 h, five colonies were picked to inoculate 10 mL of Clostridium, Staphylococcus, Pseudomonas (CSP) medium in culture tubes. CSP is a chemically defined medium developed to support the growth of P. aeruginosa, Staphylococcus aureus, and Clostridium perfringens as monocultures or consortia 34 (supplementary material S1). For CSP supplemented with one or more organic acids, 22 mM of each organic acid specified was added.
A total of three culturing tubes containing 10 mL of CSP were each inoculated with about five colonies from the TSA plates, incubated at 37 °C with shaking at 150 rpm (tubes were placed at a 45° angle in the shaker to increase mixing) and grown until the cultures reach an OD 600 of 0.5. 1 mL of each culture was then added to 49 mL of fresh CSP medium in 250 mL baffled flasks and an OD 600 of 0.010. The baffled flasks were capped with gas permeable foam lids and incubated at 37 °C with shaking at 150 rpm. Sampling occurred about every hour during the first 12 h and less frequently afterwards.
Culture sampling. Samples were drawn from each flask for OD 600 , pH, amino acid, and carbon metabolite measurements. An aliquot of 1.5 mL of culture was collected at each sampling, cells were separated from the supernatant using centrifugation at 7000 rpm for 10 min (Eppendorf 5415D microcentrifuge). Supernatants were then filtered using 0.22 µm syringe filters prior to being stored at − 20 °C.
At each sampling, a volume of culture was collected for OD 600 measurement. OD 600 readings were blanked with fresh CSP and samples were diluted, if necessary, to keep OD 600 measurements ≤ 0.30.
Organic acid and sugar analyses. HPLC analysis of select carbon metabolites including glucose and organic acids was performed with an Agilent 1200 series HPLC equipped with a refractive index detector (RID) and an Aminex HPX-87H ion exclusion column, 300 mm × 7.8 mm. A mobile phase of 5 mM H 2 SO 4 was run at a flow rate of 0.6 mL/min for 25 min/injection. A volume of 200 µL of sample was added to an HPLC vial with 200 µL of an internal standard of 1 g/L fucose dissolved in 10 mM H 2 SO 4 . Each sample was injected twice for a total of two technical replicates for each of the three biological replicates for each time point. HPLC analysis of culture supernatants were compared to NMR analysis (Chenomx library and an internal standard of DSS) as a verification of metabolite identities and to ensure no major metabolites were being missed 34 . Supplementary material S3 includes a comparison of HPLC-and NMR-based analyses of select metabolites during batch growth.
Amino acid analysis. HPLC analysis of amino acids was performed with an Agilent 1100 series equipped with a diode array detector (DAD) and a ZORBAX Eclipse XDB-C18 column, 4.6 mm ID × 250 mm (5 µm) 80 Å. This setup was used with the Agilent protocol for HPLC analysis of amino acids 98 .
Cell dry weight measurement. A correlation curve between OD 600 and grams of cell dry weight (g CDW) per liter was constructed. 5 mL aliquots of P. aeruginosa culture harvested in mid-exponential growth phase diluted to a range of densities were dried at 80 °C for 24 h in aluminum drying pans and weighed. Correlation equation was: (g CDW/L) = (OD 600 )*2.23 −1 .
Proteomics analyses. P. aeruginosa cultures (n = 3) were collected via centrifugation at 3600×g and washed three times with phosphate buffered solution to remove residual media. Cells were resuspended in 1.25 mL radioimmunoprecipitation assay buffer consisting of 12.5 μL of protease inhibitor (Halt Protease Cocktail Inhibitor, Thermo Fisher Scientific, Rockford, IL) to prevent enzymatic degradation upon cell lysis, 0.1 mg/ml lysozyme to solubilize the cell peptidoglycan layer, and 5 mM dithiothreitol (DTT) to cleave protein disulfide bonds. Cells were lysed mechanically in a beadbeater (Mini-beadbeater-1, Biospec Products, Inc, Bartlesville, OK) at 4800 oscillations/min with the remainder of the vial filled with 0.1 mm diameter zirconia/silica beads for a total time of 2.5 min (five cycles at a duration of 30 s each with chilling in an ice water bath between cycles).
Protein concentrations were determined by protein assay kit (DC Bradford Reagent, Thermo Fisher Scientific). 15 µg of proteins were taken from all samples and transferred to centrifugal filter units (Microcon-30 kDa Centrifugal Filter Unit with Ultracel-30 membrane, Millipore Sigma, Billerica, MA). The samples were then processed following the filter aided sample preparation method 99 , in which proteins were reduced with DTT, alkylated with iodoacetamide to prevent disulfide bond reformation, then enzymatically digested overnight into peptides with trypsin at 1:50 enzyme:substrate (w:w) at 37 °C. The resulting peptides were desalted using a C18 column (Macro SpinColumn, Harvard Apparatus, Holliston, MA), dried in a centrifugal evaporator, then resuspended in 5% acetonitrile with 0.1% formic acid to a concentration of 0.2 μg/μL. 200 ng samples of peptides were separated in a HPLC system (1260 Infinity LC System, Agilent, Santa Clara, CA) and a C18 column (3.5 μm particle size, 150 mm length × 75 μm internal diam, Zorbax 300SB, Agilent) using a 60 min mobile phase gradient ranging from 5 to 85% organic (0.1% formic acid in water to 0.1% formic acid www.nature.com/scientificreports/ in acetonitrile) at 250 nL/min flow rate for a total run time of 75 min. Following LC separation, peptides were ionized by nanoelectrospray with a spray voltage of 1.90 kV and 275 °C capillary temperature, then analyzed in a high resolution Orbitrap mass spectrometer (Orbitrap Velos Pro purchased in 2007, Thermo Scientific, Waltham, MA) with automatic gain control set at 10 6 ions and injection times of 1-200 ms. Full scan mass spectra from m/z 400 to 2000 at 30,000 mass resolution were collected in data-dependent acquisition mode. Ten precursor ions were selected from each full mass scan for analysis by tandem mass spectrometry (MS/MS, using 30% energy in HCD mode for fragmentation). Raw mass spectra data files were processed for protein identification using the MaxQuant software (v. 1.5.3.30) 100 with main search parameters of 4.5 ppm peptide tolerance, 20 ppm MS/MS match tolerance, 10 ppm MS/MS de novo tolerance, seven minimum peptide length, carbamidomethyl as fixed modification, 0.01 FDR, oxidation and acetylation variable modification, and enabled search for contaminants. Protein abundances were further data processed and normalized with log transformation (base 2) for data visualization and statistical analyses (ANOVA, p value < 0.05) using the Perseus software (v.1.5.4.0) 101 . The Search Tool for Retrieval of Interacting Genes (STRING) database (v. 10.5) 102 was used for protein-protein interactions amongst the statistically significant proteins, set at medium confidence of 0.4 and protein annotation (functional enrichment analyses).
In silico analysis of metabolism and resource allocation. A genome-scale, stoichiometric model of P. aeruginosa (iMO1086) 40,41 was analyzed using flux balance analysis (FBA) via the COBRA Toolbox (https ://openc obra.githu b.io/cobra toolb ox/stabl e/cite.html) in MATLAB using the Gurobi optimization program (http://www.gurob i.com) (supplementary material S12,S20). Carbon and O 2 limitations were modeled by setting the carbon (5 mmol/g/h) or O 2 (20 mmol/g/h) uptake rates, respectively, for each of the examined carbon sources and maximizing the production of biomass or cellular energy (i.e. quantified as the number of ATP bonds hydrolyzed). Enzyme limitation was modeled by minimizing the number of participating reactions for specified substrate uptake rates while producing biomass or cellular energy. Suboptimal solutions between the minimal total flux and the maximum product yield (biomass or cellular energy) or the minimum proteome and the maximum product yield (biomass or cellular energy) were identified by minimizing an aggregate objective function. The aggregate objective function was the sum of either total flux or total proteome and a weighted flux through the substrate transport reaction of interest. Optimization between proteome investment and product yield was achieved by changing the weight of the flux through the carbon or oxygen transport reaction of interest. The algorithms can be found in supplementary material S20.