Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

Open Access Research article

Modeling and analysis of flux distributions in the two branches of the phosphotransferase system in Pseudomonas putida

Andreas Kremling1*, Katharina Pflüger-Grau1, Max Chavarría2, Jacek Puchalka35, Vitor Martins dos Santos4 and Víctor de Lorenzo2

Author affiliations

1 Fachgebiet Systembiotechnologie, Technische Universität München, Garching b. München, Germany

2 Systems Biology Program, Centro Nacional de Biotecnología-CSIC, Campus Cantoblanco, Madrid, Spain

3 , Helmholtz Center for Infection Research, Braunschweig, Germany

4 , Wageningen University & Research centre, Agrotechnology & Food Sciences, Wageningen, Netherlands

5 , Present address: University Children’s Hospital

For all author emails, please log on.

Citation and License

BMC Systems Biology 2012, 6:149  doi:10.1186/1752-0509-6-149

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/6/149


Received:18 October 2011
Accepted:31 October 2012
Published:6 December 2012

© 2012 Kremling et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Signal transduction plays a fundamental role in the understanding of cellular physiology. The bacterial phosphotransferase system (PTS) together with the PEP/pyruvate node in central metabolism represents a signaling unit that acts as a sensory element and measures the activity of the central metabolism. Pseudomonas putida possesses two PTS branches, the C-branch (PTSFru) and a second branch (PTSNtr), which communicate with each other by phosphate exchange. Recent experimental results showed a cross talk between the two branches. However, the functional role of the crosstalk remains open.

Results

A mathematical model was set up to describe the available data of the state of phosphorylation of PtsN, one of the PTS proteins, for different environmental conditions and different strain variants. Additionally, data from flux balance analysis was used to determine some of the kinetic parameters of the involved reactions. Based on the calculated and estimated parameters, the flux distribution during growth of the wild type strain on fructose could be determined.

Conclusion

Our calculations show that during growth of the wild type strain on the PTS substrate fructose, the major part of the phosphoryl groups is provided by the second branch of the PTS. This theoretical finding indicates a new role of the second branch of the PTS and will serve as a basis for further experimental studies.

Keywords:
Phosphotransferase System (PTS); Flux Balance Analysis (FBA); Kinetic modelling; Metabolic Control Analysis (MCA)

Background

Mathematical modelling of biological processes is a powerful tool towards the thorough understanding of a biological system. In the mathematical simulation, in the first step, experimental data is reproduced and subsequently, the model can be used to predict the behaviour of the system. This type of iterative model-based analysis is a hallmark of systems biology research that, in the future, is expected to be very helpful in enhancing the understanding of cellular systems in a better way. Here, we chose to analyse the PTSNtr of P. putida with the help of mathematical tools, in order to obtain a functional model that describes the system in detail and may provide new ideas on the physiological role of the PTSNtr and its cross-talk with the sugar PTS.

Pseudomonas putida KT2440 is a ubiquitous Gram-negative, saprophytic soil bacterium, important in biotechnological and systems biological research. This organism is characterized by its metabolic versatility, which enables the strain to use a variety of natural and man-made compounds as carbon and energy source [1], making it an ideal organism to be used in a broad range of biotechnological applications, such as bioremediation or biotransformation processes [2-5].

On one hand, it can degrade degrade a variety of toxic compounds including methoxylated or hydroxylated aromatic acids. On the other hand, P. putida is able to degrade fructose [6] and glucose [6-8], two of the most abundant sugars present in plant root exudates [9]. Glucose catabolism occurs in this strain by the simultaneous operation of three peripheral pathways that converge at the level of 6-phosphogluconate, which is then further metabolized via the Entner-Doudoroff pathway [8]. Fructose degradation additionally occurs via the Embden-Meyerhof-Parnas pathway [6].

Central metabolism of P. putida is different in comparison to e.g. E. coli. However, a common structural metabolic unit that is found in these two organisms, as well as in other bacteria, is the phosphoenolpyruvate-dependent: phosphotransferase system (PTS). In E. coli the PTS is the main uptake system for glucose and for other carbohydrates like mannitol, mannose, and N-acetyl-glucosamine. Besides transport properties, the PTS together with the reactions of the central metabolism represents an important signaling unit that measures the glycolytic and gluconeogenetic flux and maps it to the degree of phosphorylation of the PTS proteins [10]. These proteins are involved in a number of signal transduction processes coordinating metabolism and motility of the bacteria [11].

Like many other eubacteria, P. putida KT2440 possesses also an alternative PTS, which is encoded by the genes ptsP, ptsN and ptsO. This unusual PTS preserves the phosphotransfer domains but lacks the permease component, ruling out its participation in sugar intake. Because ptsN and ptsO map adjacently to the rpoN gene that encodes nitrogen σ factor 54, it was generally believed that this alternative PTS was related to N metabolism. More recent data suggest that this PTS branch predominantly controls the influx of potassium ions [12,13], which can indirectly affect other processes [14]. Furthermore, it could well happen that most reported nitrogen-related phenotypes of PTS-Ntr mutants of E. coli could be artifactual [15], thereby raising doubts on the correct designation of the system [16]. While we take notice of current controversy on the name of this PTS branch -and notwithstanding a possible future change in the nomenclature, we will use the name PTS-Ntr [17] for the rest of the article to designate the system composed of the proteins PtsP (EI-Ntr), PtsO (NPr) and PtsN (EIIA-Ntr). Here, we applied a mathematical approach to estimate the flow of the high energy phosphate through the PTSNtr. This was performed in order to get an idea of how the phosphate flow is distributed between the two PTS (PTSFru and PTSNtr) and to shed some light on the role of both systems in the process of fructose uptake.

Mathematical modeling of the E. coli PTS has progressed in the last decade and several models are currently available (reviewed in [10,11]). The models provide evidence that - depending on the type of transport - a clear relationship between the specific growth rate and the degree of phosphorylation of the PTS proteins is valid for a broad range of specific growth rates. Moreover, it could be shown, that the PTS together with glycolytic reactions represents a robust structure in E. coli[18]. These results suggest a general principle of cellular control by the PTS that should also be found in other bacteria. Therefore, we set out to analyze in detail the PTS of P. putida by mathematically modeling the degree of phosphorylation of its proteins. Since steady-state conditions for the intracellular metabolites (e.g. during a (more or less) constant specific growth rate) contain little information on kinetic parameters, experimental data from a perturbation experiment was used. In these experiments exponential growth of P. putida on casamino acids (CAA) was compared with growth on CAA plus fructose. Since uptake of fructose takes place in the late exponential phase [19], this can be seen as a small perturbation and results in a different intracellular steady-state. In a first step - based on data for the wild type strain and mutant strains - a mathematical model was developed that links two input parameters, namely the PEP/pyruvate ratio and the fructose uptake rate to the output, the degree of phosphorylation of one PTS protein, namely PtsN. With the data available it was possible to estimate all equilibrium constants for the biochemical reactions. Based on these initial findings the behavior of the system was predicted for different conditions and the results could be verified in new experiments. In a second step, for both environmental conditions (CAA and CAA plus fructose), flux balance analysis (FBA) was performed based on a stoichiometric model that was introduced previously [20]. Combining results from the PTS model and FBA, new insights into the kinetic properties of the involved enzymes were found. Specifically, we found that the metabolic PEP/pyruvate ratio is faithfully reflected in the phosphorylation state of PtsN, thereby suggesting a mechanism by which the gross physiological state is translated into multiple downstream consequences.

Biological background

Pseudomonas putida possesses two distinct phosphoenolpyruvate (PEP) phosphotransferase systems. One is responsible for fructose uptake (carbohydrate branch, PTSFru) and the other is suggested to be involved in signal transduction. In both branches, a high-energy phosphate from PEP is transferred by a number of proteins either to the incoming fructose or to PtsN the last protein in the second branch. As depicted in Figure 1, cross talk may occur and the phosphoryl group from the C branch is transferred to the second branch. Growth of the wild type and three isogenic mutants in the corresponding pts genes of P. putida in minimal media with different C sources (casamino acids (CAA), CAA plus fructose and CAA plus glucose) revealed different degree of phosphorylation of PtsN. The figure also shows the reactions considered in the mathematical model.

thumbnailFigure 1. The PTS reactions in P. putida. The C branch of the PTS is shown on the left hand side from PEP (reaction r4) while the second branch is shown on the right hand side (reactions r1r3). The phosphryl group from the final protein in the second branch, PtsN is not metabolized. However, there is evidence that cross talk occurs between the two branches. The C branch transfers the phosphryl group to the second branch directly to PtsN (reaction r5).

By in vivo phosphorylation studies using dielectric breakdown and subsequent Western blotting, the direction of the flow of the high-energy phosphate within the PTSNtr was recently confirmed to succeed in the direction PEP → PtsP → PtsO → PtsN in analogy to the sugar PTS [19,21]. Furthermore, it was shown that both PTS systems do not proceed strictly separated from each other, but that the PTSNtr and the PTSFru communicate with each other by phosphate exchange under specific metabolic conditions [22]. Phosphorylation of PtsN by PTSFru was restricted to conditions in which the PTSFru was active [22], i.e. it only occurred when cells were grown on fructose, but not on carbon sources that enter the central carbon degradation pathway below the level of glycolysis, like succinate. It is assumed that the FruB protein delivers the phosphate moiety directly to PtsN via its HPr domain, as phosphorylation could also be shown in a ptsO mutant strain [22]. Thus, via this cross talk information about the nature of the degraded sugar or the degradation pathway is integrated directly into the phosphorylation state of PtsN. On the other hand it was shown recently that not only the carbon source, but also other physiological conditions, as the nitrogen source, or growth stage influences phosphorylation of PtsN [19]. Monitoring the phosphorylation state of PtsN in the presence of different carbon or nitrogen sources at different growth stages revealed that phosphorylated PtsN was detectable at all growth stages and strongly accumulated in the stationary growth phase, while the non-phosphorylated form exclusively appeared during rapid growth. However, non-phosphorylated PtsN did not disappear at the same growth stage in every case. Instead, the decline of non-phosphorylated PtsN happened earlier in the presence of fructose than in the presence of glucose and earlier in the presence of nitrate than with ammonia. Therefore, it is assumed that the phosphorylation state of PtsN reflects the physiological status of the cell [19]. The catabolic reaction network in P. putida is slightly different from E. coli, e.g. glucose is metabolized a different way. Figure 2 shows the main reactions including several uptake systems. The central nodes are PEP and pyruvate that serve as a “turning device” that distributes all incoming fluxes, irrespective from their source to their specific metabolic modules. Carbohydrates as well as aromatic carboxylic acids like benzoate feed in at various nodes. The results from FBA (see below) suggest that the main input node from casamino acid uptake is oxaloacetate.

thumbnailFigure 2. Central metabolic reactions of P. putida. PEP and pyruvate are hubs that distribute the fluxes to other parts of the network.

Methods

Model for the phosphotransferase system

For the involved PTS proteins shown in Figure 1 the respective differential equations are as follows (only the phosphorylated components are considered):

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M1">View MathML</a>

(1)

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M2">View MathML</a>

(2)

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M3">View MathML</a>

(3)

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M4">View MathML</a>

(4)

Reactions r1 to r3 describe the phosphoryl transfer from PEP via PtsP and PtsO to PtsN, reaction r4 describes the phosphoryl transfer needed for fructose uptake, and reaction r5 describes the cross talk between the two branches. A general form of mass action kinetics for the reaction rate is used. The ri read:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M5">View MathML</a>

(5)

where X (e.g. PtsP) is the phosphoryl donor, XP the phosphorylated form (e.g. PtsPP) and Y , YP are the unphosphorylated and the phosphorylated recipient (e.g. PtsO, PtsOP). Conservation relations can be taken into account for the four proteins:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M6">View MathML</a>

(6)

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M7">View MathML</a>

(7)

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M8">View MathML</a>

(8)

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M9">View MathML</a>

(9)

The system can be described with five equilibrium constants Ki, five velocity reactions ki and constant entities representing the overall concentration of the proteins.

Scaling is a well known tool to reduce the number of parameters. Scaling on the overall concentration Xi,0 of the respective compound and on a chosen time constant leads to a set of equations with seven velocity constants, five equilibrium constants Ki and a scaled uptake rate <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M10">View MathML</a>. System inputs are the PEP/pyruvate ratio and the fructose uptake rate. The respective equations are given in the Additional file 1. Matlab-files to simulate the system are available on request.

Additional file 1. The files describes the model equations and the kinetic parameters.

Format: PDF Size: 137KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Model for flux balance analysis

The flux distributions were computed by applying Flux Balance Analysis (FBA) on the genome-scale metabolic reconstruction of P. putida metabolism iJP815 [20]. Cellular networks are characterized by the observation that there are more unknown reaction rates than equations (for every compound a mass balance equation can be set up). To circumvent this problem an objective function is defined and from all solutions that are admissible the one that maximizes the objective function is chosen. In many cases the maximization of the yield or the maximization of ATP generating fluxes are a good choice, however, different objective functions are possible as well and are analyzed in a systematic way [23]. The following procedure was applied here: The experimentally measured specific growth rates and uptake rate of fructose constituted the constraints of the reconstruction. First, the flux distribution for bacteria growing solely on fructose was identified by setting the specific growth rate to the experimentally measured value and minimizing the uptake rate of fructose. In the first round, however, the experimentally measured uptake rate could not be reached, due to the values of the maintenance parameters. In order to cope with this problem the Non-Growth Associated Maintenance (NGAM) parameter was decreased, so as to obtain the experimentally measured value. The same NGAM value was used for the computations for the remaining conditions (CAA+fructose, CAA). In both cases the flux distributions were ascertained by limiting the specific growth rate and the fructose uptake (where applicable) to the experimentally measured values and minimizing the uptake of CAAs. The latter was achieved by introducing a virtual compound “carbon” along with a set of reactions, each converting “carbon” into a particular amino acid with the stoichiometry defined by the number of carbon atoms possessed by the particular amino acid, and, subsequently, minimizing the amount of carbon needed to be fed into the system. Flux variability analysis (FVA) was performed to show that the optimal flux distribution has narrow ranges and that alternate optima with a complete different pattern can be excluded.

Kinetics

Based on the experimental data [19,22], flux distributions for different growth conditions are available that fulfill the steady state conditions. Information on basic kinetic parameters can be obtained with the following approach: Metabolic control analysis is applied to identify rate parameters that, in principle, can be estimated. This is illustrated with the following simple network with two metabolites and four reactions given in Figure 3. To determine all fluxes at least two fluxes have to be measured. This results in two sub-matrices with unknown rates (Nu) and known rates (Nkn). Introducing kinetic rate laws allows the calculation of the respective elasticities ∈;ij=ri/cj. If it is assumed that the two input fluxes are not independent, but are related by a factor (an assumption that is reasonable while taking into account that for higher uptake rates (e.g. ra) the demands of the cell are changing (e.g. rb)), the following condition holds true for the concentration control coefficients dcj/dri:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M11">View MathML</a>

(10)

thumbnailFigure 3. Network with two metabolites and four reactions. Given is also the stoichiometric matrix.

with <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M12">View MathML</a> the derivative of the known rates with respect to the chosen input flux u, and ∈; is the matrix with all elasticities. To determine the concentration control coefficients matrix, matrix ∈; has to be invertible. This requires at least two entries different from zero. Since specific entries of matrix ∈; are usually not known, an inversion of the matrix ∈; is only possible if it has the maximal structural rank. To determine the structural rank, we refer to the following definition [24]: “The structural rank (s-Rank) of a structural matrix is the maximum rank a linear matrix of this structure can have.” Therefore, it can be determined by finding the maximum size set of nonzero entries that do not share a row or column of the matrix. If ra and rd are the measured rates, and assuming that all remaining reactions rb(Pyr), rc(Pyr) are irreversible, matrix ∈; is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M13">View MathML</a>

(11)

where stars mark the respective dependency. The matrix has structural rank one. In case of ra and rb are measured, the respective matrix can be inverted in (nearly) all cases since it is of the form

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M14">View MathML</a>

(12)

Having the concentration control coefficients at hand, they are used to determine the kinetic properties of the system.

Strains and experimental conditions

Strains, media, and growth conditions

All Pseudomonas strains used in this work were derived from strain Pseudomonas putida MAD2. P. putida MAD2 variants bear directed chromosomal insertions of each gene (ptsP, ptsO, ptsN, fruB), with either a kanamycin (Km) resistance gene or the xylE marker and all mutants used (ptsP, ptsO, ptsN, fruB) in this work have been described previously [25,26]. For determination of physiological parameters, bacteria were grown at 30°C in 250 ml baffled shake flasks with 50 ml synthetic M9 medium [27] with 0.2% CAA and 0.2% glucose or fructose as C-source.

Analytical procedures and physiological parameters

Cell growth was monitored spectrophotometrically at 600 nm (OD600) and fructose/ glucose concentrations were determined enzymatically with the fructose/ glucose assay kit (Sigma-Aldrich) according to the supplier’s manual. The following physiological parameters were determined by regression analysis during the exponential growth phase in batch culture, as described elsewhere [28,29]: maximal specific growth rate, biomass yield on glucose (or fructose), specific glucose (or fructose) consumption. The correlations factors (rc) between cellular dry weight (CDW) and OD600 were determined from batch cultures of each mutant. Therefore, CDW was measured from at least three parallel 10 ml cell suspensions by harvesting the cells by fast filtration through pre-weight nitrocellulose filters (0.45 μm), which were subsequently washed with 0.9% NaCl and dried at 105°C for 24 h to a constant weight. We have followed the procedures described in [30-32].

Determination of the phosphorylation state of PtsN

The phosphorylation state of PtsN in P. putida Mad2 and pts mutants grown on CAA and glucose was determined by dielectric breakdown and Western blotting as described elsewhere [22].

Results and discussion

The PTS model describes the available data points

Since metabolic reactions are faster than the growth rate, a pseudo-steady-state can be considered. For the complete model the number of data points (see section Methods) are not sufficient to determine all kinetic parameters directly. However, using the data for the mutant strains allows to decompose the system in a smaller set of equations. Table 1 summarizes restrictions for all reaction rates ri for the different conditions. The following procedure was applied: The PEP/pyruvate ratio for growth only on CAA was set to 1. As can be seen in Figure 4, mutation in the respective genes of the system prevents the flow of the phosphoryl groups and allows to formulate conditions for some of the kinetic parameters (details for the calculation of the parameters are given in the Additional file 1). Furthermore, the complete system contains a cycle. This means that the second branch and the C branch of the PTS have the same overall equilibrium constant:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M15">View MathML</a>

(13)

thumbnailFigure 4. Network fluxes under different conditions (black lines indicate active fluxes). A The network for the wild type strain growing on CAA as well as FruB mutant growing on CAA plus fructose. B The network of PtsP and PtsO mutants growing on CAA and fructose.

Table 1. Summary of the restrictions for the determination of the kinetic parameters

providing an additional constraint for the parameters.

Figure 5 shows an comparison of the measured values and the simulation results. Shown are two bars for every condition tested: growth on CAA (4 strains), growth on CAA plus fructose (4 strains), and growth on CAA plus glucose (only wild type). For growth on CAA it assumed that the fructose operon is not induced, that is, the value for the overall concentration of FruB is set to a low value.

thumbnailFigure 5. Comparison of the experimental data with the simulated data. Growth conditions as indicated on the x-axis. Black bars represent the simulation results while the grey ones the experimental data. Given is the degree of phosphorylation of PtsN. Kinetic parameters are summarized in the Additional file 1 .

The degree of phosphorylation of PtsN of the mutant strains PtsP and PtsO on fructose is characterized by the cross talk between the C and the second branch of the PTS. Although the phosphoryl flux in these mutants is interrupted in the second branch, a certain degree of phosphorylation of PtsN is detected. This can only be explained with a transfer of the phosphoryl group from FruB to PtsN (reaction r5). In contrast, for growth on CAA, it is assumed that the fru operon is not induced and the levels of phosphorylation for both mutants are low. The PEP pyruvate ratio was estimated to:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M16">View MathML</a>

(14)

For growth on CAA plus glucose, the values for ppCAA + Glc and the concentration of FruB represent degrees of freedom and have to be adjusted to describe the data. From the available data on the degree of phosphorylation, FruB is assumed to be present also on growth on glucose, however, it has to be further analyzed whether this is a specific effect of glucose or an artefact of residual presence of fructose in the medium composition.

Model prediciton and verification

To check the performance of the model, the behavior of the model for growth of P. putida under different environmental conditions was predicted. Here, growth of mutant strains PtsP, PtsO, and FruB on glucose was predicted and finally verified experimentally. Table 2 summarizes the respective values. Figure 6 shows the comparison of the simulation results with the new experimental data. For the PtsP and PtsO mutant, again, the cross talk between the C and the second branch of the PTS can be seen. The model predicts the values for all conditions very well.

thumbnailFigure 6. Comparison of the experimental data with the predicted simulated data for growth in glucose. Black bars represent the simulation results while grey bars are the experimental data. Given is the degree of phosphorylation of PtsN.

Table 2. Degree of phosphorylation of PtsN grown on CAA plus glucose

Flux distributions at nodes PEP and pyruvate

FBA is widely used to explore the capabilities of genome-scale networks (for a review, see [33]). However, it was also shown that FBA can be used to estimate flux distributions in cellular networks. In particular, for the genome scale model available for P. putida it could be shown, that there is a good agreement between the calculated values and measured values [20].

FBA can be used to estimate flux distributions given a stoichiometric network, uptake and production rates for compounds in the medium and an objective function. Here, the network was analyzed for growth on CAA and growth on CAA plus fructose. The complete flux distribution can be found in the Additional file 2. Figure 7 highlights the flux distributions at nodes PEP and pyruvate. Shown are the values for growth on CAA plus fructose (upper value) and for CAA only (lower value). The number of reactions where either PEP or pyruvate or both are involved for these conditions is 15. The main part of the fluxes is coming via oxaloacetate to the node pyruvate and is then further distributed to acetyl CoA and PEP. Minor drain fluxes are directed to alanine and other compounds. Note that in the stoichiometric network only the metabolic part is considered (C branch) but not the detailed reactions for the second branch of the PTS. To check the results for alternate optima, FVA was performed and the respective minimal and maximal values are computed. Table 3 summarizes the values for the fluxes where PEP and pyruvate are involved. As can bee seen, the minimal and maximal values show only a small range. The deviation between minimal and maximal values based on the nominal value is smaller than 3% except for the drain from pyruvate to acetyl CoA. This indicates that alternate optima with a complete different flux distribution could be excluded.

Additional file 2. The xls-sheet gives all values for the flux distributions based on a FBA analysis for different growth conditions. The xls-sheet is subdivided into two parts.In the first part, the values for the growth conditions are summarized; afterwards the flux for each reaction for the different conditions is provided. For all conditions, three columns are shown: given are the nominal value, and minimal/maximal values based on FVA. Reactions with PEP and pyruvate as substrate (S) or product (P) are highlighted.

Format: XLS Size: 248KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

thumbnailFigure 7. Flux map at node PEP pyruvate for growth conditions of the wild type on CAA plus fructose (upper values) and CAA only (lower values). The ratio of PEP and pyruvate determines the degree of phosphorylation of PtsN.

Table 3. Summary of FVA

Kinetic properties at PEP and pyruvate node

Having the results from the previous sections at hand, a mathematical model taking into account the kinetics of reactions where PEP and pyruvate are involved was set up. The focus is on the steady-state values of intracellular metabolites PEP and pyruvate in dependence on the input flux. Therefore, fluxes that produce PEP or pyruvate or have PEP or pyruvate as substrate were summed up. This resulted in a network with 5 fluxes and two nodes as shown in Figure 8 on the left side. Based on the structure of the network, the approach introduced above (Figure 3), and the FBA results, kinetic expressions were set up for reactions rc and rd, since in this case matrix ∈; has full rank. To perform the calculation for different input fluxes ra, the values for rb had to be interpolated as shown in Figure 8 to estimate a factor k (the relation between the two fluxes, see Methods). A linear interpolation was used and the values for rb are given as follows:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M17">View MathML</a>

(15)

thumbnailFigure 8. Left: Reduced flux map at node PEP/pyruvate. The network is represented with 5 reactions rates and 2 nodes. Right: An increasing input flux ra leads to increasing other fluxes. The increase in flux rb is calculated with the slope <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M18">View MathML</a> .

with <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M19">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M20">View MathML</a> are the respective values for growth in CAA plus fructose and k is the slope as shown in Figure 8. For the following considerations, the very low PTS flux is not considered; this results in the following steady-state conditions:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M21">View MathML</a>

(16)

The elasticity matrix is in the form:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M22">View MathML</a>

(17)

and the respective concentration control coefficient are:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M23">View MathML</a>

(18)

For the unknown reaction rates rc and rd power law kinetics are applied:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M24">View MathML</a>

(19)

and the respective elasticities are

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M25">View MathML</a>

(20)

From the observation that for increasing input flux ra the PEP/pyruvate ratio is also increasing, a necessary condition for the positive slope can be calculated from the following condition:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M26">View MathML</a>

(21)

and therefore

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M27">View MathML</a>

(22)

Inserting Equations (18) and (20) result in:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M28">View MathML</a>

(23)

which give first insights for the kinetic parameters. Since for the (relative) simple model, the steady-state values for PEP and pyruvate can be calculated, an estimation of the kinetic parameters is possible. With rc=rd, αkd=kc, and n=1/nd−1/nc the PEP/pyruvate ratio is:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/149/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/149/mathml/M29">View MathML</a>

(24)

Having two measurements available, two values (α and n, corresponding to kd and nd) could be calculated when kc=1 and nc=1 (this choice is reasonable since we are mainly interested in a comparison of the two enzymes rather than in absolute values).

Figure 9 shows the relation between the PEP/pyrvuate ratio and the input flux.The curve for growth on CAA is shown. In case of growth on CAA plus fructose, the flux through the PTS is very small. Therefore, the curve is nearly identical and was not shown in the figure. Based on the flux values reported in Table 3, value k (see Figure 8) might change also which results in a different behavior of the kinetics. In Figure 9 two extreme cases were also simulated (dashed lines) representing the largest and the smalles value of k.

thumbnailFigure 9. Characteristic curve for the dependency of the PEP/pyruvate ratio on the input flux ra. The two measured data points are indicated with a line.

Interrelation between the two branches of the PTS

Having both parts of the model available, an estimation of the complete flux distribution, considering fluxes through both parts of the PTS branches, is possible. Table 1 shows all equilibrium conditions for the individual reaction rates. As can be seen, for the wild type strain and growth on CAA plus fructose, none of the individual reaction is in equilibrium (column WT, row CAA + Fru). Therefore this experiment could not used for the determination of the kinetic parameters. The equation system in this case (representing a steady-state) given in the Additional file 1 (Equation (4)) is valid to determine the individual fluxes given all kinetic parameters. Taking into account that the equations are scaled (as described in the Additional file 1), a recalculation of the fluxes, based on literature values for the intracellular components, leads to the flux distribution shown in Figure 10. The most interesting result is that – under the investigated condition of a small fructose uptake rate – approx. 78% of the required phosphoryl groups for the fructose uptake are provided by the second branch of the PTS. The strength of this coupling is unexpected since the C branch of the PTS was thought to be an autonomous module. Furthermore, there are no theoretical studies with E. coli available that show that the second branch of the PTS is involved in carbohydrate transport. The explanation is as follows: The overall equilibrium constant requires that both branches (PEP to PtsN via the C branch, and PEP to PtsN via the second branch) have the same value. From the experimental data, cross talk occurs with a rather low value of the equilibrium constant K5. This leads to a rather small value for K4 and corresponding to a high affinity of PEP to FruB. The high affinity leads to a high phosphorylation of FruB, 99%, hence, the substrate, unphosphorylated FruB, is rare and consequently, the flux r4 is low. This result is rather unexpected and poses new questions on the function of cross talk.

thumbnailFigure 10. Complete flux distribution at node PEP and pyruvate for growth on CAA plus fructose. Note, that the major part for the phosphorylation of incomming fructose comes from the second branch of the PTS.

Conclusion

The PTS is a key player in the coordination of catabolic reactions in E. coli[11]. The second branch of this signaling unit recently gained attention [34] since it is involved in the regulation of nitrogen metabolism and carbon assimilation. Moreover, this unit plays also a role in potassium homoeostasis by regulating a potassium transport system [12]. Several mathematical models for the sugar phosphotransferase system in E. coli are published [10,18] that reveal a clear relationship between the input flux of carbon and the degree of phosphorylation of the PTS protein EIIA. This was shown by a comparison of experimental data for growth on different PTS and non-PTS sugars. However, kinetic data for the N branch of the system can hardly be found, since measurement of the degree of phosphorylation of the respective proteins “in vivo” are missing. In contrast, for Pseudomonas putida experimental data on the phsophorylation state of the PTS protein PtsN is available from a perturbation experiment and this data was used to set up and validate a mathematical model.

The model comprises two parts: The first one relates the PEP/pyruvate ratio for different growth conditions of a wild type strain and mutant strains (PtsP, PtsO, FruB) to the degree of phosphorylation of the PTS protein PtsN. Based on the information for the different mutant strains, the complete set of equations could be simplified for every growth condition/strain and the equilibrium constants were calculated or chosen to describe the experimental data at best. Moreover, it was possible to estimate the PEP/pyruvate ratio for growth on CAA plus fructose if the PEP/pyruvate ratio for growth on CAA alone was given. The experimental results are in good agreement with the simulated data; this is reflected by the fact that the mean of the residuals is in the range of 10% of the measured values (the errors in the experimental data are between 5% and 30%). The experimental data revealed differences in the phosphorylation state of PtsN between the wild type and the FruB mutant strain for growth on CAA as well as for the PtsP and PtsO mutant for growth on fructose. These differences could not be reproduced by the model, since in both cases the mathematical equations are the same. However, it was assumed in the model that the PEP/pyruvate ratio is the same for one substrate. So, the differences could be explained with slightly different values of the PEP/pyruvate ratio for growth on the same substrate. The overall equilibrium constant for the PTS was estimated to be Keq=K1K2K3=K4K5=0.02. This value is a factor 35 smaller than reported for the E. coli PTS [18]. However, the value used in [18] is based on different studies with enzymes in vitro and therefore a fair comparison is not possible. The equilibrium constant for the cross talk between the two branches, parameter K5, is 654.6, indicating a very weak connection between the two branches.

Based on the kinetic parameters obtained from the experiments with the wild type P. putida strain and the mutant strains, the PEP/pyruvate ratio for growth in CAA plus glucose was adjusted to describe the data. Since the degree of phosphorylation of PtsN under this condition is much lower than for growth on CAA or CAA plus fructose, the PEP/pyruvate ratio is an order of magnitude smaller ppCAA + Glc=0.05. Based on the data, a prediction of the degree of phosphorylation of PtsN was performed for different mutant strains (PtsP, PtsO, and FruB). The experimental data could be reproduced very well.

Having PEP/pyruvate ratios for growth on CAA and CAA plus fructose available, in a second step, kinetic properties of the system were analyzed based on the respective flux distributions. As expected the flux pattern was different for the two conditions. The main fluxes from or to the both nodes PEP and pyruvate are shown in Figure 7. The uptake of fructose results in a decrease of the main flux to pyruvate and simultaneously in the flux from PEP to the upper part of glycolysis, while other fluxes to alanine and others remain more or less constant. FVA was performed and indicates only small ranges for the minimal and maximal values. This excludes alternate optima with a complete different flux distribution. Taking power law kinetics into consideration for the two rates rc and rd and using the other fluxes as input fluxes, conditions for the hill coefficient for both enzymes could be calculated. Furthermore, setting the two parameters for enzyme PEP synthase (rc) to 1 the ratio of the parameters of the enzyme for rd (representing gluconeogenetic reactions) is calculated. As a result, the hill coefficient for the gluconeogenetic reactions is smaller then the one for PEP synthase (factor 0.30) while the reaction constant k is higher (factor 2.73). Taking into account the values from flux variability analysis a confidence region for the kinetics shown in Figure 9 can be given.

An interesting observation was seen when calculating all reaction rates in case of growth on CAA plus fructose. In this case, the individual reactions are not in equilibrium since fructose as additional input enhances the PTS reaction. Approx. 78% of the required phosphoryl groups for the fructose uptake are provided by the second branch of the PTS by cross talk.

This result allows to speculate on a complete new and unexpected role for the PTSNtr, at least in P. putida. It might not only be involved in the regulation of various processes, as the activity of central enzymes [35], the accumulation of polyhydroxyalcanoates [29], or the expression of the toluene degradation pathway [25], but serve also as the main provider of phosphoryl groups for fructose uptake. One could speculate that by having the phosphoryl groups cycling in the PTSNtr, they provide a storage system of rapidly available phosphate. This system comes into action when fructose is provided to the cells. Fructose in P. putida is degraded by both pathways, the EMP and the ED pathway [36], whereas glucose is degraded almost exclusively via the ED pathway in P. putida as the fructose-6-phosphate kinase, able to catalyze the conversion of fructose-6-phsophate to fructose-1,6-biphosphate and thereby making the connection to the EMP pathway, is missing [6]. Degradation of fructose via the EMP pathway is thermodynamically more favourable than degradation via the ED pathway. It is estimated that per molecule fructose degraded through the EMP pathway one molecule ATP is gained compared to degradation following the ED pathway [37]. Thus, the PTSNtr might serve as a “pre”-adaptation to the potential presence of fructose, enabling the cell to rapidly and efficiently metabolize fructose, when it is available. This hypothesis is currently under investigation in our laboratory.

Modeling of signal transduction units together with genome scale stoichiometric models will help for a better understanding of the cellular system. Especially the PTS is an important system that is involved in the coordination of catabolic reactions. The proposed model is good starting point to extend research in direction of the coordination between the carbon and other networks.

Competing interest

The authors declare that they have no competing interests.

Authors’ contributions

AK performed the modeling of the PTS, KPG and MC performed the experiments, JP and VMdS performed the flux balance analysis, AK and VdL designed the study, AK and KPG wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

AK was funded in part by the FORSYS initiative from the German Federal Ministry of Education and Research (BMBF).

References

  1. Nelson KE, et al.: Complete genome sequence and comparative analysis of the metabolically versatile Pseudomonas putida KT2440.

    Environ Microbiol 2002, 4:799-808. PubMed Abstract | Publisher Full Text OpenURL

  2. Jimenez JI, et al.: Genomic insights in the metabolism of aromatic compounds in Pseudomonas. In Pseudomonas. Edited by Ramos JL. Plenum Publishers; 2004:425-462. OpenURL

  3. Rojo F, et al.: Assemblage of ortho cleavage route for simultaneous degradation of chloro- and methylaromatics.

    Science 1987, 238:1395-1398. PubMed Abstract | Publisher Full Text OpenURL

  4. Timmis KN, et al.: Designing microorganisms for the treatment of toxic wastes.

    Annu Rev Microbiol 1994, 48:525-557. PubMed Abstract | Publisher Full Text OpenURL

  5. Walsh UF, et al.: Pseudomonas for biocontrol of phytopathogens: from functional genomics to commercial exploitation.

    Curr Opin Biotechnol 2001, 12:289-295. PubMed Abstract | Publisher Full Text OpenURL

  6. Velazquez F, et al.: Genetic evidence that catabolites of the Entner-Doudoroff pathway signal C source repression of the sigma54 Pu promoter of Pseudomonas putida.

    J Bacteriol 2004, 186:8267-8275. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. del Castillo T, Ramos JL: Simultaneous catabolite repression between glucose and toluene metabolism in Pseudomonas putida is channeled through different signaling pathways.

    J Bacteriol 2007, 189:6602-6610. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. del Castillo T, et al.: Convergent peripheral pathways catalyze initial glucose catabolism in Pseudomonas putida: genomic and flux analysis.

    J Bacteriol 2007, 189:5142-5152. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Kamilova F, et al.: Organic acids, sugars, and L-tryptophane in exudates of vegetables growing on stonewool and their effects on activities of rhizosphere bacteria.

    Mol Plant Microbe Interact 2006, 19:250-256. PubMed Abstract | Publisher Full Text OpenURL

  10. Kremling A, et al.: Catabolite repression in E. coli: A comparison of modeling approaches.

    FEBS J 2009, 276:594-602. PubMed Abstract | Publisher Full Text OpenURL

  11. Deutscher J, et al.: How Phosphotransferase system-related protein Phosphorylation regulates carbohydrate metabolism in bacteria.

    Microbiol Mol Biol Rev 2006, 70(4):939-1031. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Lee CR, et al.: Escherichia coli enzyme IIANtr regulates the K+ transporter TrkA.

    PNAS 2007, 104:4124-4129. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Lüttmann D, et al.: Stimulation of the potassium sensor KdpD kinase activity by interaction with the phosphotransferase protein IIANtr in Escherichia coli.

    Mol Microbiol 2009, 72:978-994. PubMed Abstract | Publisher Full Text OpenURL

  14. Lee CR, et al.: Potassium mediates Escherichia coli enzyme IIA(Ntr) - dependent regulation of sigma factor selectivity.

    Mol Microbiol 2010, 7:1468-1483. OpenURL

  15. Reaves ML, Rabinowitz JD: Characteristic phenotypes associated with ptsN-null mutants in Escherichia coli K-12 are absent in strains with functional ilvG.

    J Bacteriol 2011, 193:4576-4581. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Ninfa A: Unnecessary signaling: Poorly named?

    J Bacteriol 2011, 193(18):4571-4573. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Powell BS, et al.: Novel proteins of the phosphotransferase system encoded within the rpoN operon of Escherichia coli: Enzyme IIANtr affects growth on organic nitrogen and the conditional lethality of an erats mutant.

    J Biol Chem 1995, 270:4822-4839. PubMed Abstract | Publisher Full Text OpenURL

  18. Kremling A, et al.: Analysis of global control of Escherichia coli carbohydrate uptake.

    BMC Syst Biol 1:42. OpenURL

  19. Pflüger K, de Lorenzo V: Growth-dependent phosphorylation of the PtsN (EIINtr) protein of Pseudomonas putida.

    J Biol Chem 2007, 282:18206-18211. PubMed Abstract | Publisher Full Text OpenURL

  20. Puchalka J, et al.: Genome-scale reconstruction and analysis of the Pseudomonas putida KT2440 metabolic network facilitates applications in biotechnology.

    PLOS Comput Biol 2008, 4:e1000210. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Pflüger K, et al.: Non-disruptive release of Pseudomonas putida proteins by in situ electric breakdown of intact cells.

    J Microbiol Methods 2007, 71:179-185. PubMed Abstract | Publisher Full Text OpenURL

  22. Pflüger K, de Lorenzo V: Evidence of in vivo cross talk between the nitrogen-related and fructose-related branches of the carbohydrate phosphotransferase system of Pseudomonas putida.

    J Bacteriol 2008, 190:3374-3380. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli

    Mol Syst Biol 2007, 3:119. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Steffen T: Control reconfiguration of dynamical systems. Springer; 2005. OpenURL

  25. Cases I, et al.: The IIANtr (PtsN) protein of Pseudomonas putida mediates the C source inhibition of the sigma54-dependent Pu promoter of the TOL plasmid.

    J Biol Chem 1999, 274:15562-15568. PubMed Abstract | Publisher Full Text OpenURL

  26. Cases I, de Lorenzo V: The black cat/white cat principle of signal integration in bacterial promoters.

    Embo J 2001, 20:1-11. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Miller JH: Experiments in Molecular Genetics. New York: Cold Spring Harbor Laboratory; 1972. OpenURL

  28. Sauer U, et al.: Metabolic flux ratio analysis of genetic and environmental modulations of Escherichia coli central carbon metabolism.

    J Bacteriol 1999, 181:6679-6688. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Velazquez F, et al.: The phosphotransferase system formed by PtsP, PtsO, and PtsN proteins controls production of polyhydroxyalkanoates in Pseudomonas putida.

    J Bacteriol 2007, 189:4529-4533. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Ostergaard S: Increasing galactose consumption by Saccharomyces cerevisiae through metabolic engineering of the GAL gene regulatory network.

    Nat Biotechnol 2000, 18:1283-1286. PubMed Abstract | Publisher Full Text OpenURL

  31. Rapid quantifiable assessment of nutritional parameters influencing pediocin production by Pediococcus acidilactici NRRL B5627

    Bioresource technology 2008, 99:6646-6650. PubMed Abstract | Publisher Full Text OpenURL

  32. Kylma AK, et al.:

    Microbial Biotech. 2004, 63:584-591. Publisher Full Text OpenURL

  33. Feist AM, Palsson BO: The growing scope of applications of genome-scale metabolic reconstructions using Escherichia coli.

    Nat Biotechnol 2008, 26:659-67. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Pflüger-Grau K, Görke B: Regulatory roles of the bacterial nitrogen-related phosphotransferase system.

    Trends Microbiol 2010, 18:205-214. PubMed Abstract | Publisher Full Text OpenURL

  35. Pflüger-Grau K, et al.: The interplay of the EIIA(Ntr) component of the nitrogen-related phosphotransferase system (PTSNtr) of Pseudomonas putida with pyruvate dehydrogenase.

    Biochim Biophys Acta 2011, 1810:995-1005. PubMed Abstract | Publisher Full Text OpenURL

  36. Sawyer MH, et al.: Pathways of D-fructose catabolism in species of Pseudomonas.

    Arch Microbiol 1977, 112:49-55. PubMed Abstract | Publisher Full Text OpenURL

  37. Van Dijken JP, Quayle J: R: Fructose metabolism in four Pseudomonas species.

    Arch Microbiol 1977, 114:281-286. PubMed Abstract | Publisher Full Text OpenURL