Fluxome study of Pseudomonas fluorescens reveals major reorganisation of carbon flux through central metabolic pathways in response to inactivation of the anti-sigma factor MucA
- Stina K Lien†1,
- Sebastian Niedenführ†2,
- Håvard Sletta3,
- Katharina Nöh2Email author and
- Per Bruheim1Email author
© Lien et al.; licensee BioMed Central. 2015
Received: 28 August 2014
Accepted: 27 January 2015
Published: 18 February 2015
The bacterium Pseudomonas fluorescens switches to an alginate-producing phenotype when the pleiotropic anti-sigma factor MucA is inactivated. The inactivation is accompanied by an increased biomass yield on carbon sources when grown under nitrogen-limited chemostat conditions. A previous metabolome study showed significant changes in the intracellular metabolite concentrations, especially of the nucleotides, in mucA deletion mutants compared to the wild-type. In this study, the P. fluorescens SBW25 wild-type and an alginate non-producing mucA- ΔalgC double-knockout mutant are investigated through model-based 13C-metabolic flux analysis (13C-MFA) to explore the physiological consequences of MucA inactivation at the metabolic flux level. Intracellular metabolite extracts from three carbon labelling experiments using fructose as the sole carbon source are analysed for 13C-label incorporation in primary metabolites by gas and liquid chromatography tandem mass spectrometry.
From mass isotopomer distribution datasets, absolute intracellular metabolic reaction rates for the wild type and the mutant are determined, revealing extensive reorganisation of carbon flux through central metabolic pathways in response to MucA inactivation. The carbon flux through the Entner-Doudoroff pathway was reduced in the mucA- ΔalgC mutant, while flux through the pentose phosphate pathway was increased. Our findings also indicated flexibility of the anaplerotic reactions through down-regulation of the pyruvate shunt in the mucA- ΔalgC mutant and up-regulation of the glyoxylate shunt.
Absolute metabolic fluxes and metabolite levels give detailed, integrated insight into the physiology of this industrially, medically and agriculturally important bacterial species and suggest that the most efficient way of using a mucA- mutant as a cell factory for alginate production would be to use non-growing conditions and nitrogen deprivation.
Pseudomonas is a bacterial genus containing species of industrial, agricultural and medical relevance due to their metabolic diversity and ability to colonise a wide variety of ecological niches, including not only soil and water but also insects, plants, and animals [1,2]. The biosynthetic diversity of this genus includes production of the polysaccharide alginate, which has many commercial applications, especially in food and medicine , but alginate is a complicating factor to the human host during Pseudomonas infections .
Alginate is not produced by Pseudomonas fluorescens SBW25 wild type, but its production can be induced by inactivation of the anti-sigma factor MucA . MucA inactivation is a frequent mechanism behind the conversion of P. aeruginosa to a mucoid phenotype infecting the lungs of cystic fibrosis patients . In its active form, MucA prevents alginate biosynthesis by binding and sequestering AlgU, the sigma factor necessary for transcription of the alg operon. The alg operon codes for all of the structural genes needed for alginate synthesis except AlgC , the role of which is to convert mannose 6-phosphate (M6P) to mannose 1-phosphate (M1P). M6P is produced from fructose-6-phosphate (F6P), implying that alginate synthesis draws from the hexose-phosphate pool of central metabolism. In both P. aeruginosa and P. fluorescens, MucA inactivation has been shown to affect numerous genes other than those involved in alginate biosynthesis [5,7,8].
In our previous metabolome study , P. fluorescens SBW25 wild type and mucA-mutants were cultivated on fructose and glycerol as sole carbon sources. Intracellular concentrations of a wide range of metabolic intermediates, including phosphometabolites, amino acids, organic acids, and nucleotides, were measured. The study revealed significant carbon-source-dependent differences in primary metabolite and nucleotide levels for both the wild type and the mucA mutants. A dramatic decrease in energy charge (EC) was found in MucA mutant strains, which was accompanied by an increased growth yield (YX/S). Interestingly, no counter response action of the energy-deprived mutants’ levels of central intermediates and metabolic genes was observed. This observation led to the hypothesis that mutants do not perceive the low EC as a critical energy limitation. To understand this puzzling robustness, the reorganisation of carbon flows within the biological network of P. fluorescens strains upon MucA inactivation needs to be elucidated. This is the aim of the present study.
In contrast to other “-omics” technologies, metabolic fluxes cannot be directly measured but must be inferred from carbon labelling experiments (CLEs) in combination with mathematical modelling . In CLEs, 13C-enriched carbon substrates are fed to the cells and, driven by the cells’ metabolic and enzymatic activities, isotopically labelled carbon atoms are incorporated into the metabolite pools . The incorporation of labelled material causes all central metabolic intermediates to exist as a set of different isotopomers, which are a collection of isomers that have different labelling patterns of 12C and 13C carbon atoms. Labelling data along with physiological data on observed fluxes (uptake and production fluxes) then enable intracellular fluxes to be quantified through 13C-metabolic flux analysis (13C-MFA) or, analogous to other “omics”-type analyses, fluxomics . During the last few decades, 13C-MFA has become a mature tool to quantify the actual in vivo physiological state of microorganisms, as proven by numerous application studies from the fields of metabolic engineering and systems biology [13-15].
Published fluxome studies for Pseudomonas species are scarce. A flux ratio analysis comparing seven bacterial species (P. fluorescens wild type, P. putida, Agrobacterium tumefaciens, Sinorhizobium meliloti, Rhodobacter sphaeroides, Zymomonas mobilis, and Paracoccus versutus) with the two model organisms E. coli and Bacillus subtilis has been performed with glucose as the carbon source . CLEs were carried out in shake flasks, and samples were taken during exponential growth to derive the relative net flux contribution of converging pathways in selected branch points of the metabolic network from the labelling patterns of protein-bound amino acids. The authors showed that compared to the model organisms, the seven investigated strains, including P. fluorescens, behaved significantly different, as they all mainly relied on the Entner-Doudoroff pathway (EDP) and used the pentose phosphate pathway (PPP) only for biosynthetic functions. In addition, the aerobic bacteria had no by-product formation and to a larger extent used the pyruvate shunt (malate (MAL) to pyruvate (PYR) to oxaloacetate (OAA)) for conversion of MAL to OAA instead of direct conversion by malate dehydrogenase . In another study, intracellular fluxes for P. aeruginosa elucidated alginate production using glucose as the carbon source [17-19]. Significant involvement of the EDP in alginate biosynthesis was shown; glucose was shuttled via triose phosphate glyceraldehyde 3-phosphate (GAP) and dihydroxyacetone phosphate (DHAP) to alginate.
The objective of the current study is to quantify the effect of MucA inactivation on the primary metabolism of P. fluorescens through determination of the levels of absolute fluxes. For comparability with experiments promoting alginate production  on the one hand and the conditions of our previous metabolomics study  on the other hand, nitrogen limitation was imposed and fructose was used as the carbon source. To investigate the effect of a mucA mutant, we focus on the flux reorganisation of the primary metabolism. However, alginate production starts upon MucA inactivation , which affects the metabolic fluxes due to the removal of the alginate precursor metabolite M6P, concurrently with the effects of MucA inactivation itself. To separate these two effects, we use a mucA- ΔalgC double-knockout strain in which the direct effect of MucA inactivation can be observed without indirect “distortion” by alginate production. It was previously shown that a single deletion of algC does not alter gene expression  or metabolite profiles  compared to the wild-type strain. The determined absolute flux distributions are compared with data from the previous metabolome study conducted with the same strains. The comparison of both flux maps and metabolite levels generates complementary views of the primary metabolism of P. fluorescens. In particular, flux maps reveal new insights into the mechanisms by which primary metabolism is controlled by MucA, and how MucA may influence alginate production through the re-routing of carbon flow. In summary, this study provides the first integrative perspective on the metabolism of P. fluorescens and lays the groundwork for the future metabolic engineering of this organism.
Chemostat cultivations and carbon labelling experiments (CLE)
To quantify the changes in carbon flows in P. fluorescens upon MucA inactivation, the wild-type strain and an alginate non-producing mucA- ΔalgC double-knockout mutant were grown with fructose as the sole carbon source under nitrogen limitation. For 13C-MFA, it is essential to ensure that cell physiology and metabolism are held in a pseudo-steady state throughout the experiment. Therefore, CLEs are conducted in chemostats with growth rates of 0.04 h-1 under defined environmental conditions. CLEs were performed in two rounds: results derived from one wild-type cultivation (a “pre-CLE”) were used for an experimental design study to determine a statistically maximally informative fructose isotopomer composition, as detailed in the Materials and Methods section. In turn, the calculated mixture was administered in the main CLEs to the two strains.
Measured and calculated cultivation data for P. fluorescens CLEs
Specific fructose uptake [mmolC/gDW h]
Specific alginate production [mmolC/ gDW h]
Specific biomass production [mmolC/gDW h] 3)
Specific CO 2 excretion [mmol/gDW h] 4)
Y xs [%]
Mass spectrometry measurements
Metabolite extracts from the CLEs were analysed by LC-MS/MS and GC-MS/MS to generate mass isotopomer datasets. Both chromatographic techniques were used to increase the information content and the coverage of labelling patterns for flux estimation. GC-MS was mainly used to analyse soluble free amino acids, while LC-MS was used for central carbon metabolites (organic acids and sugar phosphates). Based on the MS/MS datasets, intracellular fluxes were estimated for both strains.
Quantification of metabolic fluxes for the wild-type and mucA- ΔalgC strains
A model of P. fluorescens central carbon metabolism was formulated that incorporated reactions for fructose uptake (Carbon Uptake), glycolysis/Embden – Meyerhof – Parnas pathway (EMP), the pentose phosphate pathway (PPP), the Entner-Doudoroff pathway (EDP), the tricarboxylic acid cycle (TCA) and glyoxylate shunt, the anaplerotic reactions (ANA), as well as amino acid biosynthesis pathways (BM/BS), covering the measured metabolite pools (cf. Additional file 1 for specification of the metabolic network model including carbon atom transitions). Details of the 13C-MFA modelling workflow are given in the Materials and Methods section.
Although the wild-type and the mucA- ΔalgC strains both utilise the same main route for initial fructose catabolism, there are differences between the two strains at important branch points. One of the most striking distinctions occurs at the 6-phosphogluconate (6PGN) node, where the EDP and the PPP split. In the wild-type strain, the major carbon flux from the 6PGN node is directed to the EDP (61 ± 8%, EC 188.8.131.52), whereas the minor proportion goes to the PPP (39 ± 7%, EC 184.108.40.206) (Figures 1 and 3B). However, the situation is reversed for the mucA- ΔalgC strain. As shown in Figures 2 and 3B, the major carbon flux goes to the PPP (71 ± 17%) and the minor efflux for this strain is directed to the EDP (29 ± 5%). Notably, even with an almost 40% decrease in the fructose uptake rate, the mucA- ΔalgC strain still has a higher absolute flux through the PPP (1.6 ± 0.3 mmol/(gDW h) versus 1.18 ± 0.20 mmol/(gDW h) for the wild-type strain). It is thus clear that inactivating the anti-sigma factor MucA leads to a reduced flux through the EDP and an increased flux through the PPP.
The second important branch point that is found to be distinct for the two strains is isocitrate (ICIT) in the TCA. The wild-type strain does not utilise the glyoxylate (GLX) shunt to any significant extent (flux value of 0% from ICIT to GLX with an asymmetric confidence interval [0, 4]% corresponding to an absolute confidence interval of [0.00,0.10] mmol/(gDW h)), while the mucA- ΔalgC strain shuttles 22 ± 7% (0.33 ± 0.16 mmol//(gDW h)) of the ICIT influx to GLX (cf. Figures 1, 2 and 3C). The fructose uptake rate and EDP flux are lower for the mucA- ΔalgC strain than for the wild-type, leading to a significantly lower influx to ICIT compared to the wild type (1.51 ± 0.09 mmol/(gDW h) vs. 3.10 ± 0.17 mmol/(gDW h)). Thus, the utilisation of the GLX shunt in response to MucA inactivation in turn leads to an even further reduction in flux through the remainder of the TCA.
Finally, comparison of the flux maps reveals differences between the two strains, although the current dataset does not allow for the precise estimation of ANA flux values (see broad confidence intervals in Additional file 2). Nevertheless, the results indicate that the mucA- ΔalgC strain might cycle phosphoenolpyruvate (PEP) to OAA and back, with the two net fluxes almost cancelling each other out (to replace the drain of TCA intermediates to biomass, a small net flux to TCA remains). In contrast, for the wild-type strain, the flux from PEP to OAA might be significantly higher than the reverse flux (cf. Figures 1, 2 and inset of Figure 3D portraying the magnitude of the net fluxes). The results also show that although the PYR shunt (MAL → PYR → OAA) is utilised to a significant extent by the wild-type strain, shuttling 28 ± 9% (0.84 ± 0.27 mmol/(gDW h)) of the MAL influx to OAA via PYR, this shunt is significantly less active in the mucA- ΔalgC strain (a flux of 13% or 0.20 mmol/(gDW h) within the confidence range of [0, 0.45] mmol/(gDW h) of the MAL influx goes to PYR) (Figures 1, 2 and main chart of Figure 3D). Thus, down-regulation of the pyruvate shunt must be a consequence of MucA inactivation.
NADH and NADPH production is significantly higher in wild type
In contrast, the net production of NADH and NADPH differs for the two strains, both being higher for the wild type (Figure 4A). However, the relative proportion of NADH to NADPH production is approximately the same, even though there are large differences in major intracellular carbon flow between the two strains (Figure 4B, with data normalised to the fructose uptake rate). In the wild-type strain, the production of NADH is significantly higher because of higher uptake of fructose, more of which is shuttled through the EDP and to the TCA cycle than in the mucA- ΔalgC strain. Likewise, NADPH production is higher for the wildtype, again predominantly because of the wildtype’s higher TCA flux (isocitrate dehydrogenase is NADPH-dependent, EC. 220.127.116.11). On a relative scale, compared to the influx of the network (Table 1), NADPH production is similar in both strains. Because biomass production compared to substrate uptake is higher for the mucA- ΔalgC strain, the wild-type strain generates a surplus of NADPH that is not needed for anabolic purposes.
Pseudomonas fluorescens is an attractive organism in industrial, medical and agricultural biotechnology. The first step towards optimising alginate production with this microorganism is to accurately characterise the influencing factors and underlying metabolic regulation mechanisms that reorganise the carbon flow of alginate-producing mutant strains in comparison to the wildtype. This study investigated how the effect of MucA modulates the carbon flow under conditions favouring alginate production. To this end, information-optimised CLEs were designed and, in turn, conducted in well-controlled chemostats under steady-state conditions. Mass isotopomer data were collected from LC- as well as GC-coupled MS devices, generating comprehensive labelling datasets with high coverage. Labelling data, in combination with measured uptake and biomass formation rates, enable the precise determination of absolute intracellular flux rates by 13C-MFA. Finally, the fluxome of P. fluorescens is related to previously determined transcriptome and metabolome data.
Our fluxome study revealed significant reorganisation of major intracellular metabolic fluxes due to inactivation of the pleiotropic anti-sigma factor MucA. The most significant findings from the mucA- ΔalgC double-knockout mutant are as follows: 1. a reduction in EDP flux, 2. an increase in PPP flux, and 3. an activation of the GLX shunt. The data for the wildtype pre-experiment and the two cultivations of the main experiment (Table 1) coincide well with previously published data for P. fluorescens cultivations under similar conditions [5,9], thus confirming that the results from these down-scaled CLE cultivations (120 mL cultivations) can be correlated with transcriptome and metabolome results generated at standard lab-scale operating conditions (750 mL). The experimental set-up with nitrogen-limited chemostats provides the sole carbon source, fructose, in excess. The cells can therefore adjust their carbon source consumption to their metabolic requirements.
Reduced fructose uptake rate and EP pathway flux in the mucA- ΔalgC mutant with lowest energy charge
Interestingly, the wild-type strain has a 60% higher fructose uptake rate than the mucA- ΔalgC strain (Figures 1, 2 and Table 1). Because the growth rate of both strains is the same, the biomass yield on substrate (YXS) is 60% higher in the mucA- ΔalgC strain. Thus, the wildtype either uses more of the fructose for maintenance processes, has a less energy-efficient metabolism (i.e., lower P/O ratio), and/or disposes of surplus energy through futile cycling, which was frequently described for other organisms such as Escherichia coli, Bacillus subtilis, and Corynebacterium glutamicum [16,21-24].
The reduction in fructose uptake in the mucA- ΔalgC strain cannot be a direct effect of MucA inactivation alone, as an alginate-producing mucA- single-deletion mutant had a two-fold increased fructose uptake rate in the metabolome study (18.3 mmol C/(gDW h) in the mucA- mutant vs. 9.4 mmol C/(gDW h) in the wild type) . MucA inactivation thus seems to exert its effect on the EDP, the PPP and the TCA rather than on the uptake system and the initial reactions of F6P and G6P alone.
The previous metabolome study  showed metabolite pool changes by MucA inactivation in the presence and absence of alginate production. The most striking result was a significantly decreased EC ((ATP + 0.5ADP)/(ATP + ADP + AMP)) in both the alginate-non-producing mucA- ΔalgC mutant used in this fluxome study and an alginate-producing mucA- single-knockout mutant . The reduction in EC was caused both by a decrease in ATP concentration and an increase in ADP and AMP concentrations, the latter almost 10 times higher. The higher EC of the wild type correlates with its larger NADH production (Figure 4). However, the mucA- ΔalgC strain does not exploit the full capacity to increase its EC in the way the wild type does by, for example, increasing fructose uptake. The likely reason for this is that mucA- mutants do not perceive their relatively low EC as a stress situation, and therefore, no cellular adaptation response is stimulated.
In theory, P. fluorescens could utilise glycolysis alone for fructose assimilation, converting FBP to PYR and generating two net ATP. By converting FBP to PYR via the EDP, no net ATP is generated, as no ATP is formed via fructose-1,6-bisphosphatase (EC 18.104.22.168). From this perspective, both P. fluorescens strains preferably use the less energy-efficient EDP route. Flamholz and co-workers analysed glycolysis and the EDP in terms of thermodynamic and kinetic constraints and concluded that glycolysis incurs a higher protein cost (i.e., more enzymatic proteins are needed to sustain the same flux) than the EDP . The lower protein cost of the EDP might especially favour this pathway under the studied conditions because the chemostats were nitrogen limited, raising the pressure on protein synthesis. Thus, nitrogen limitation seems advantageous for alginate production on fructose.
Concerning the preferred route for fructose uptake, it is also interesting to note that FBPase in P. fluorescens is negatively regulated by AMP, as it is in other bacteria, e.g., E. coli . Because alginate synthesis starts with the conversion of F6P to M6P, the induction of alginate synthesis in the mucA- mutant must alleviate the possible rate-limitation of the FBPase. The eight-fold increase in the mucA- mutants’ AMP concentration could have led to the direct conversion of FBP to DHAP and GAP by FBP aldolase (EC 22.214.171.124) in the mucA- ΔalgC strain, but that was not observed in the current fluxome study.
Disposal of surplus energy in the wild-type strain
According to our fluxome study, a large proportion of the excess fructose uptake of the wild type is shuttled via EDP to TCA and complete oxidation to CO2, with energy conservation taking place through oxidative phosphorylation. This occurs with no obvious use of this extra energy compared to the mucA- ΔalgC strain. Based on current knowledge of co-factor dependencies , the NADPH as well as NADH production are higher in the wild type (Figure 4). Strikingly, the NADPH production is higher than needed for anabolic processes in the wild type because the biomass production rate is the same for both strains. The cyclic operation of anaplerotic reactions in the wild type points to futile cycling dissipating energy by hydrolysing energy equivalents, as is also hypothesised for Corynebacterium species [27-29]. The wild type uses the PYR shunt, which generates one NADPH and consumes one ATP, which is a less energy-conserving solution than the direct oxidation of MAL to OAA with the NADH-producing MAL dehydrogenase reaction. It has been previously demonstrated that in batch cultivations growing on glucose, P. fluorescens primarily relies on the PYR shunt for conversion of MAL to OAA in contrast to other Gram-negative (E. coli, P. putida, A. tumefaciens, R. sphaeroides and Z. mobilis) and Gram-positive (B. subtilis and P. versutus) bacteria . Thus, the pyruvate shunt might be active in situations where NADH production is not required.
Yet there remain open questions: how does the cell regenerate or allocate NADP if needed? Activation of a transhydrogenase (EC 126.96.36.199) may provide an explanation to this question . Due to the lower fructose uptake rate of the mucA- ΔalgC strain, it is unlikely that this mutant is energy limited. It is puzzling that the mutant activates the GLX shunt, which is mostly associated with growth on two-carbon substrates, and generates one net OAA from acetyl-CoA without any energy conservation. Additionally, and similarly to the wild type, the ANA cycling of the mucA- ΔalgC strain through PEP carboxylase (EC 188.8.131.52) and PEP carboxykinase (EC 184.108.40.206) could be regarded as a futile cycle because one net ATP is consumed. However, the extent of futile cycling could not be quantified because, as mentioned, ANA reactions were not resolved in the forward and backward directions with statistical confidence. To this end, further CLEs with specially designed 13C-label mixtures are needed [28,31].
Comparison with other -omics data
The comparative investigation of changes in the transcriptome  and the fluxome of the investigated strains reveals that the strong reorganisation of central metabolic carbon flows is not tightly coupled to expression levels of the corresponding metabolic genes, which remained remarkably constant.
In this work, a mucA- ΔalgC double-knockout P. fluorescens mutant was compared to the wild-type strain at the metabolic flux level to illuminate the global metabolic response to MucA inactivation on central carbon metabolism without a biasing carbon re-routing by alginate production. Corroborating transcriptome and metabolome results from former studies, a strong influence of the mucA- ΔalgC double-knockout mutant was shown, involving fine-tuned and significant reorganisation of carbon utilisation. Salient alterations lead to a decrease of the overall fructose uptake and an activation of the PPP and the GLX shunt that is accompanied by a deactivation of the EDP. Although both strains produce a similar amount of energy (NADH and ATP) per mole of substrate, the mutant seems to exhibit more energy-efficient behaviour than the wildtype.
Clearly, to mechanistically elucidate the metabolic control mechanisms of alginate production on fluxes, mucA- mutants need to be further investigated. These experiments bear the potential to elucidate the complex control mechanisms of MucA on alginate production and central carbon metabolism. These control points may become important targets for drug development to prevent the formation of the mucoid phenotype, which leads to pulmonary infection in cystic fibrosis patients, or to increase industrial alginate production. This requires extending the multi-omics characterisation to further mutants, including catabolite repressor/activator cra mutants, and experimental conditions, i.e., growth rates and carbon sources, that are clinically and industrially relevant. From the modelling point of view, the application of non-stationary 13C-MFA in combination with advanced experimental design strategies would be advantageous to advance the integrative analysis of metabolite levels and fluxes and increase the statistical significance of the results, particularly with regard to anaplerosis.
The strains used in this study were P. fluorescens SBW25 and P. fluorescens SBW25 mucA-ΔalgC. Construction of the mutant strain was previously presented in detail . Flux distribution results for these strains can be compared to metabolome results  for the same strains grown under identical conditions. A confounding factor when using the mucA-ΔalgC strain to look for effects of MucA inactivation other than alginate production is that it produces all biosynthetic enzymes required for alginate production other than AlgC.
Cultivations and sample preparation
Chemostat cultivations were performed as described earlier , except that the working volume was scaled down from 750 mL to 120 mL in a custom-made 0.5-L reactor vessel and the carbon source concentration was scaled down from 40 g/L to be less in excess, using the equivalent of 25 g/L and 20 g/L of naturally labelled fructose for the wild type and the mucA- ΔalgC strain, respectively. These carbon source concentrations lead to a fructose concentration of 2.5 g/L or above in the fermenters at steady state. In brief, the chemostat cultivations were nitrogen limited using fructose as a carbon source with a dilution rate of 0.04 h-1. Temperature was maintained at 25°C, pH was maintained at 7.8 by automated sodium hydroxide addition, air flow was 1.5 VVM, and oxygen saturation was maintained above 20% by adjustment of the stirring speed. Cultivations were operated as batches for 24 h before initiating continuous operation with a feed containing naturally labelled fructose. After 1.5-2.0 volume exchanges, the feed stock was changed. Afterwards, the feedstock contained 13C-labelled fructose (two different lots of 1-13C-fructose with chemical purities of 99.1 and 98.0% (w/w) and isotope enrichments of 99% (mole/mole); and one lot of U-13C-fructose with chemical purity of 100% (w/w) and isotopic enrichment of 99.5% (mole/mole), Cambridge Isotope Laboratories, Inc, Andover). Sampling was performed after an additional 3 volume exchanges. Process data of double knock out mutant is given in Additional file 3.
Three cultivations were performed; the only difference between the first and the two last was the fructose isotopomer composition used in the 13C-labelled feed. In the first cultivation, the wild-type strain was grown using a 60.04% 1-13C-, 20.31% U-13C- and 19.65% 12C-fructose isotopomer composition (values based on experimental weight of dry substances adjusted for chemical impurities). In the two main cultivations, the wild type and the mucA-ΔalgC strain were grown with the designed fructose isotopomer mixture of 60.08% 1-13C-, 39.92% U-13C-fructose.
Immediately prior to sampling for metabolite extracts, an aliquot of 5 mL was taken from the culture for OD660 and residual fructose measurement. The CO2 concentration in the off gas was measured with a Rosemount Binos 100 CO2 analyser. These data were used to calculate the carbon balance as described previously . For preparation of metabolite extracts, the cultivations were sampled two (wild-type experimental design cultivation samples for LC-MS/MS analysis) or three (for all other conditions) times for each of the two MS-methods (biological replicates) to be employed, using fast vacuum filtration and quenching of metabolism in 37.5% cold methanol (-20°C). Details of the sampling method are as described in our previous metabolome study , except that only one filter with pore size 0.8 μm (Millipore cat#AAWP04700) was used instead of a stack of three filters and the culture volume filtered was 6 mL instead of 2 mL, leading to metabolites from 3 mL of culture being present in the final freeze-dried metabolite extracts.
Sample and evaluation of mass isotopomer datasets
Metabolite extracts were analysed by an LC-MS/MS method and a GC-MS/MS method to generate mass isotopomer distribution datasets for the P. fluorescens SBW25 wild-type and the mucA-ΔalgC strains. Prior to LC-MS/MS analysis, samples were reconstituted in 500 μL 60% (v/v) methanol. Analysis was performed on an Agilent 1200 series LC connected via an electrospray ion source to an Agilent 6410 triple-quadrupole MS using an adaptation of the reverse phase tributylamine ion exchange method introduced by Luo et al. . The same LC-MS/MS method was used in the metabolome study of P. fluorescens . A total of eight central carbon metabolism metabolites were analysed for mass isotopomer distributions by LC-MS/MS and used in the study. For each compound, the number of precursor-to-fragment transitions (multiple reaction monitoring (MRM)-transitions) was expanded to account for all mass isotopomers of the compound, i.e., (n-m + 1)(m + 1) MRM-transitions for each compound, where n is the number of carbons in the precursor ion and m is the number of carbons in the product ion, respectively . A list of all MRM-transitions and accompanying MS-settings for the LC-MS/MS method can be found in Additional file 4.
Samples for GC-MS/MS analysis were derivatised prior to analysis as described previously  using an adaptation of the methyl chloroformate derivatisation protocol [38,39]. Analysis was performed on an Agilent 7890A GC – 7000 triple-quadrupole MS interfaced with a chemical ionisation (CI) ion source based on a previously described MS/MS method . The GC method was as described previously , with the exception that the injected volume was increased from 1 μL to 3 μL. The MS method used positive CI with methane as the reagent gas and nitrogen as the collision gas. MRM-transitions to account for all mass isotopomers of the 14 detected compounds (central carbon metabolism metabolites and amino acids) were included in the fluxomics models. A list of all MRM-transitions and accompanying GC-MS/MS-settings can be found in Additional file 4.
The resulting mass isotopomer fractions are listed in Additional file 5. Average relative standard deviations (STD) for the mass isotopomers are calculated from the biological replicates for the estimation of intracellular fluxes. The obtained LC-MS/MS STD values were 13% for the wild-type experimental design cultivation, and 10% and 8% for the main experiment cultivations on wild type and mucA-ΔalgC, respectively. For mass isotopomers detected by the GC-MS/MS method, these STD values were 10%, 12% and 15%, respectively.
Computational 13C-MFA: modelling, flux estimation and design of informative CLEs
A biochemical reaction network of central carbon metabolism of P. fluorescens SBW25 was derived from the genome-scale metabolic model iSB1139 reconstructed by Borgos et al.  as well as from the KEGG database (http://www.genome.jp/kegg). Construction of our model was guided by 1) representing degradation pathways by effluxes from the system and excluding low-abundance biomass components (see Additional file 1: Tables S.1.2 and S.1.4); 2) pooling metabolites due to limited measurement information, e.g., 1,3-bisphosphoglycerate, 2-phosphoglycerate and 3-phosphoglycerate, which are represented by one lumped pool 23PG while preserving the carbon atom transitions; and 3) aggregating parallel fluxes, such as isoenzymes. Like other Pseudomonas strains, P. fluorescens does not possess a phosphofructokinase enzyme for conversion of F6P to FBP; thus, only the gluconeogenetic conversion of FBP to F6P is included in the model (EC 220.127.116.11).
The model of P. fluorescens central carbon metabolism consists of routes for fructose uptake (Carbon Uptake), EMP, EDP, PPP, TCA cycle including the GLX shunt, the ANA section, as well as the biomass equation, which defines the net fluxes of metabolites into amino acid and fatty acid synthesis. To formulate the latter effluxes, the biomass composition for P. putida devised by Nogales et al.  was adopted and integrated as a set of constraints for the fluxes into building block synthesis (Additional file 1). In total, the model contains 66 metabolites and 116 reactions (18 thereof reversible), each supplemented with carbon atom transitions (Additional file 1). The model has 29 degrees of freedom (13 net and 15 exchange fluxes) to be estimated from the measured datasets generated for the P. fluorescens wild type and the mucA-ΔalgC mutant strain.
The generation of isotopomer balances, the simulation of measurements, the experimental mixture design, the flux estimation (minimisation of a weighted least squares objective repeated 100 times with randomly sampled starting values combined with a globalised optimisation strategy to detect multiple equally good but essentially different flux solutions), and the statistical assessment of resulting flux confidence intervals were performed using the software tool 13CFLUX2 . For further details of the computational 13C-MFA protocols, we refer to classical review papers, e.g., . Finally, the resulting flux estimations are represented in the context of the metabolic network model using the visualisation software Omix .
An optimal experimental design (OED) study to determine an information-optimised fructose label mixture was performed based on the fluxes estimated from the pre-CLE with P. fluorescens SBW25 wild type. The main CLEs for the wild type and the mucA-ΔalgC strain were performed using the calculated OED mixture consisting of 60.08% 1-13C-, 39.92% U-13C-fructose isotopomers. When inspecting the data thoroughly with the isotopomer network at hand, some peak measurements have to be excluded due to low peak quality. Specifically, we used the following procedure to achieve an acceptable fit: 1) we assumed a minimum STD of 1% to account for the noise of measurements close to the detection limit of the MS instruments, 2) we obtained the best fit by minimising the difference between simulated and measured data, and 3) we excluded those measurements with low abundance that have extraordinary large discrepancies (values for experimental and simulated measurements are found in Additional file 6). Optimisation runs provided a unique optimum for each dataset and flux estimation. For the best-fitting flux maps, the statistical analyses were performed with constraining several main exchange fluxes (as indicated by the absence of uncertainties for fluxes in Additional file 2.
This study was financed by an internal grant at NTNU (SKL) and the Norwegian Research Council via the Era-Net SysMO program. The authors thank R. Aune for technical assistance during the 13C-fructose cultivations.
- Mahajan-Miklos S, Rahme LG, Ausubel FM. Elucidating the molecular mechanisms of bacterial virulence using non-mammalian hosts. Mol Microbiol. 2000;37(5):981–8.View ArticlePubMedGoogle Scholar
- Poblete-Castro I, Becker J, Dohnt K, dos Santos VM, Wittmann C. Industrial biotechnology of Pseudomonas putida and related species. Appl Microbiol Biotechnol. 2012;93(6):2279–90.View ArticlePubMedGoogle Scholar
- Skjåk-Bræk G, Draget KI. Alginates: properties and applications. In: Matyjaszewski K, Möller M, editors. Polymer science: a comprehensive reference, vol. 10. Amsterdam: Elsevier BV; 2012. p. 213–20.View ArticleGoogle Scholar
- Mathee K, Ciofu O, Sternberg C, Lindum PW, Campbell JIA, Jensen P, et al. Mucoid conversion of Pseudomonas aeruginosa by hydrogen peroxide: a mechanism for virulence activation in the cystic fibrosis lung. Microbiology-Sgm. 1999;145:1349–57.Google Scholar
- Borgos SE, Bordel S, Sletta H, Ertesvag H, Jakobsen O, Bruheim P, et al. Mapping global effects of the anti-sigma factor MucA in Pseudomonas fluorescens SBW25 through genome-scale metabolic modeling. BMC Syst Biol. 2013;7:19.View ArticlePubMed CentralPubMedGoogle Scholar
- Morea A, Mathee K, Franklin MJ, Giacomini A, O’Regan M, Ohman DE. Characterization of algG encoding C5-epimerase in the alginate biosynthetic gene cluster of Pseudomonas fluorescens. Gene. 2001;278(1–2):107–14.View ArticlePubMedGoogle Scholar
- Firoved AM, Deretic V. Microarray analysis of global gene expression in mucoid Pseudomonas aeruginosa. J Bacteriol. 2003;185(3):1071–81.View ArticlePubMed CentralPubMedGoogle Scholar
- Jones AK, Fulcher NB, Balzer GJ, Urbanowski ML, Pritchett CL, Schurr MJ, et al. Activation of the Pseudomonas aeruginosa AlgU Regulon through mucA Mutation Inhibits Cyclic AMP/Vfr Signaling. J Bacteriol. 2010;192(21):5709–17.View ArticlePubMed CentralPubMedGoogle Scholar
- Lien SK, Sletta H, Ellingsen TE, Valla S, Correa E, Goodacre R, et al. Investigating alginate production and carbon utilization in Pseudomonas fluorescens SBW25 using mass spectrometry-based metabolic profiling. Metabolomics. 2013;9(2):403–17.View ArticleGoogle Scholar
- Sauer U. Metabolic networks in motion: 13C-based flux analysis. Mol Syst Biol. 2006;2:62.View ArticlePubMed CentralPubMedGoogle Scholar
- Wiechert W. 13C metabolic flux analysis. Metab Eng. 2001;3(3):195–206.View ArticlePubMedGoogle Scholar
- Weitzel M, Nöh K, Dalman T, Niedenführ S, Stute B, Wiechert W. 13CFLUX2 - high-performance software suite for 13C-metabolic flux analysis. Bioinformatics (Oxford, England). 2013;29(1):143–5.View ArticleGoogle Scholar
- Woolston BM, Edgar S, Stephanopoulos G. Metabolic engineering: past and future. In: Prausnitz JM, editor. Annual review of chemical and biomolecular engineering, Vol 4. vol. 4th ed. Palo Alto: Annual Reviews; 2013. p. 259–88.Google Scholar
- Nöh K, Wiechert W. The benefits of being transient: isotope-based metabolic flux analysis at the short time scale. Appl Microbiol Biotechnol. 2011;91(5):1247–65.View ArticlePubMedGoogle Scholar
- Zamboni N, Sauer U. Novel biological insights through metabolomics and 13C-flux analysis. Curr Opin Microbiol. 2009;12(5):553–8.View ArticlePubMedGoogle Scholar
- Fuhrer T, Fischer E, Sauer U. Experimental identification and quantification of glucose metabolism in seven bacterial species. J Bacteriol. 2005;187(5):1581–90.View ArticlePubMed CentralPubMedGoogle Scholar
- Lynn AR, Sokatch JR. Incorporation of isotope from specifically labeled glucose into alginates of Pseudomonas aeruginosa and Azotobacter vinelandii. J Bacteriol. 1984;158(3):1161–2.PubMed CentralPubMedGoogle Scholar
- Narbad A, Hewlins MJE, Gacesa P, Russell NJ. The use of 13C-NMR spectroscopy to monitor alginate biosynthesis in mucoid Pseudomonas aeruginosa. Biochem J. 1990;267(3):579–84.PubMed CentralPubMedGoogle Scholar
- Narbad A, Russell NJ, Gacesa P. Radiolabelling patterns in alginate of Pseudomonas aeruginosa synthesized from specifically-labelled 14C-monosaccharide precursors. Microbios. 1988;54(220–21):171–9.PubMedGoogle Scholar
- Wiechert W. The thermodynamic meaning of metabolic exchange fluxes. Biophys J. 2007;93(6):2255–64.View ArticlePubMed CentralPubMedGoogle Scholar
- Emmerling M, Dauner M, Ponti A, Fiaux J, Hochuli M, Szyperski T, et al. Metabolic flux responses to pyruvate kinase knockout in Escherichia coli. J Bacteriol. 2002;184(1):152–64.Google Scholar
- Petersen S, de Graaf AA, Eggeling L, Möllney M, Wiechert W, Sahm H. In vivo quantification of parallel and bidirectional fluxes in the anaplerosis of Corynebacterium glutamicum. J Biol Chem. 2000;275(46):35932–41.Google Scholar
- Sauer U, Hatzimanikatis V, Bailey JE, Hochuli M, Szyperski T, Wüthrich K. Metabolic fluxes in riboflavin-producing Bacillus subtilis. Nat Biotechnol. 1997;15(5):448–52.Google Scholar
- Yang C, Hua Q, Baba T, Mori H, Shimizu K. Analysis of Escherichia coli anaplerotic metabolism and its regulation mechanisms from the metabolic responses to altered dilution rates and phosphoenolpyruvate carboxykinase knockout. Biotechnol Bioeng. 2003;84(2):129–44.Google Scholar
- Flamholz A, Noor E, Bar-Even A, Liebermeister W, Milo R. Glycolytic strategy as a tradeoff between energy yield and protein cost. Proc Natl Acad Sci U S A. 2013;110(24):10039–44.View ArticlePubMed CentralPubMedGoogle Scholar
- Hines JK, Kruesel CE, Fromm HJ, Honzatko RB. Structure of inhibited fructose-1,6-bisphosphatase from Escherichia coli - Distinct allosteric inhibition sites for AMP and glucose 6-phosphate and the characterization of a gluconeogenic switch. J Biol Chem. 2007;282(34):24697–706.Google Scholar
- Blombach B, Riester T, Wieschalka S, Ziert C, Youn JW, Wendisch VF, et al. Corynebacterium glutamicum tailored for efficient isobutanol production. Appl Environ Microbiol. 2011;77(10):3300–10.Google Scholar
- Petersen S, Mack C, de Graaf AA, Riedel C, Eikmanns BJ, Sahm H. Metabolic consequences of altered phosphoenolpyruvate carboxykinase activity in Corynebacterium glutamicum reveal anaplerotic regulation mechanisms in vivo. Metab Eng. 2001;3(4):344–61.Google Scholar
- Sauer U, Eikmanns BJ. The PEP-pyruvate-oxaloacetate node as the switch point for carbon flux distribution in bacteria. Fems Microbiol Rev. 2005;29(4):765–94.View ArticlePubMedGoogle Scholar
- Bruheim P, Butler M, Ellingsen TE. A theoretical analysis of the biosynthesis of actinorhodin in a hyper-producing Streptomyces lividans strain cultivated on various carbon sources. Appl Microbiol Biotechnol. 2002;58(6):735–42.Google Scholar
- Möllney M, Wiechert W, Kownatzki D, de Graaf AA. Bidirectional reaction steps in metabolic networks: IV, Optimal design of isotopomer labeling experiments. Biotechnol Bioeng. 1999;66(2):86–103.View ArticlePubMedGoogle Scholar
- Saier MH, Ramseier TM. The catabolite repressor/activator (Cra) protein of enteric bacteria. J Bacteriol. 1996;178(12):3411–7.PubMed CentralPubMedGoogle Scholar
- Kochanowski K, Volkmer B, Gerosa L, van Rijsewijk BRH, Schmidt A, Heinemann M. Functioning of a metabolic flux sensor in Escherichia coli. Proc Natl Acad Sci U S A. 2013;110(3):1130–5.Google Scholar
- Sarkar D, Siddiquee KAZ, Arauzo-Bravo MJ, Oba T, Shimizu K. Effect of cra gene knockout together with edd and iclR genes knockout on the metabolism in Escherichia coli. Arch Microbiol. 2008;190(5):559–71.Google Scholar
- Chavarria M, Santiago C, Platero R, Krell T, Casasnovas JM, de Lorenzo V. Fructose 1-phosphate is the preferred effector of the metabolic regulator Cra of Pseudomonas putida. J Biol Chem. 2011;286(11):9351–9.Google Scholar
- Luo B, Groenke K, Takors R, Wandrey C, Oldiges M. Simultaneous determination of multiple intracellular metabolites in glycolysis, pentose phosphate pathway and tricarboxylic acid cycle by liquid chromatography-mass spectrometry. J Chromatogr A. 2007;1147:153–64.View ArticlePubMedGoogle Scholar
- Choi J, Antoniewicz MR. Tandem mass spectrometry: a novel approach for metabolic flux analysis. Metab Eng. 2011;13(2):225–33.View ArticlePubMedGoogle Scholar
- Smart KF, Aggio RBM, Van Houtte JR, Villas-Boas SG. Analytical platform for metabolome analysis of microbial cells using methyl chloroformate derivatization followed by gas chromatography-mass spectrometry. Nat Protoc. 2010;5(10):1709–29.View ArticlePubMedGoogle Scholar
- Villas-Bôas SG, Delicado DG, Åkesson M, Nielsen J. Simultaneous analysis of amino and nonamino organic acids as methyl chloroformate derivatives using gas chromatography-mass spectrometry. Anal Biochem. 2003;322:134–8.View ArticlePubMedGoogle Scholar
- Kvitvang HFN, Andreassen T, Adam T, Villas-Boas SG, Bruheim P. Highly sensitive GC/MS/MS method for quantitation of amino and nonamino organic acids. Anal Chem. 2011;83(7):2705–11.View ArticlePubMedGoogle Scholar
- Nogales J, Palsson BO, Thiele I. A genome-scale metabolic reconstruction of Pseudomonas putida KT2440: iJN746 as a cell factory. BMC Syst Biol. 2008;2:79.Google Scholar
- Zamboni N, Fendt SM, Ruhl M, Sauer U. 13C-based metabolic flux analysis. Nat Protoc. 2009;4(6):878–92.View ArticlePubMedGoogle Scholar
- Droste P, Nöh K, Wiechert W. Omix - a visualization tool for metabolic networks with highest usability and customizability in focus. Chemie Ingenieur Technik. 2013;85(6):849–62.View ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.