Modeling and analysis of flux distributions in the two branches of the phosphotransferase system in Pseudomonas putida
© Kremling et al.; licensee BioMed Central Ltd. 2012
Received: 18 October 2011
Accepted: 31 October 2012
Published: 6 December 2012
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.
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.
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.
KeywordsPhosphotransferase System (PTS) Flux Balance Analysis (FBA) Kinetic modelling Metabolic Control Analysis (MCA)
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 , 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  and glucose [6–8], two of the most abundant sugars present in plant root exudates . 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 . Fructose degradation additionally occurs via the Embden-Meyerhof-Parnas pathway .
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 . These proteins are involved in a number of signal transduction processes coordinating metabolism and motility of the bacteria .
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 . Furthermore, it could well happen that most reported nitrogen-related phenotypes of PTS-Ntr mutants of E. coli could be artifactual , thereby raising doubts on the correct designation of the system . 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  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. 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 , 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 . 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.
Model for the phosphotransferase system
The system can be described with five equilibrium constants K i , five velocity reactions k i 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 K i and a scaled uptake rate . 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.
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 . 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 . 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.
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  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 .
Results and discussion
The PTS model describes the available data points
Summary of the restrictions for the determination of the kinetic parameters
(r up =0)
CAA + Fru
providing an additional constraint for the parameters.
For growth on CAA plus glucose, the values for p pCAA + Glcand 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
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 ). 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 .
Summary of FVA
PEP to 2PG
CAA - FBA
CAA + Fru - FBA
Kinetic properties at PEP and pyruvate node
Having two measurements available, two values (αand n ′ , corresponding to k d and n d ) could be calculated when k c =1and n c =1(this choice is reasonable since we are mainly interested in a comparison of the two enzymes rather than in absolute values).
Interrelation between the two branches of the PTS
The PTS is a key player in the coordination of catabolic reactions in E. coli. The second branch of this signaling unit recently gained attention  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 . 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 K eq =K1K2K3=K4K5=0.02. This value is a factor 35 smaller than reported for the E. coli PTS . However, the value used in  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 p pCAA + 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 r c and r d 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 (r c ) to 1 the ratio of the parameters of the enzyme for r d (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 , the accumulation of polyhydroxyalcanoates , or the expression of the toluene degradation pathway , 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 , 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 . 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 . 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.
AK was funded in part by the FORSYS initiative from the German Federal Ministry of Education and Research (BMBF).
- Nelson KE, et al.: Complete genome sequence and comparative analysis of the metabolically versatile Pseudomonas putida KT2440. Environ Microbiol. 2002, 4: 799-808. 10.1046/j.1462-2920.2002.00366.x.View ArticleGoogle Scholar
- Jimenez JI, et al.: Genomic insights in the metabolism of aromatic compounds in Pseudomonas. Pseudomonas. Edited by: Ramos JL. 2004, Plenum Publishers, 425-462.View ArticleGoogle Scholar
- Rojo F, et al.: Assemblage of ortho cleavage route for simultaneous degradation of chloro- and methylaromatics. Science. 1987, 238: 1395-1398. 10.1126/science.3479842.View ArticleGoogle Scholar
- Timmis KN, et al.: Designing microorganisms for the treatment of toxic wastes. Annu Rev Microbiol. 1994, 48: 525-557. 10.1146/annurev.mi.48.100194.002521.View ArticleGoogle Scholar
- Walsh UF, et al.: Pseudomonas for biocontrol of phytopathogens: from functional genomics to commercial exploitation. Curr Opin Biotechnol. 2001, 12: 289-295. 10.1016/S0958-1669(00)00212-3.View ArticleGoogle Scholar
- 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. 10.1128/JB.186.24.8267-8275.2004.View ArticleGoogle Scholar
- 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. 10.1128/JB.00679-07.View ArticleGoogle Scholar
- 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. 10.1128/JB.00203-07.View ArticleGoogle Scholar
- 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. 10.1094/MPMI-19-0250.View ArticleGoogle Scholar
- Kremling A, et al.: Catabolite repression in E. coli: A comparison of modeling approaches. FEBS J. 2009, 276: 594-602. 10.1111/j.1742-4658.2008.06810.x.View ArticleGoogle Scholar
- Deutscher J, et al.: How Phosphotransferase system-related protein Phosphorylation regulates carbohydrate metabolism in bacteria. Microbiol Mol Biol Rev. 2006, 70 (4): 939-1031. 10.1128/MMBR.00024-06.View ArticleGoogle Scholar
- Lee CR, et al.: Escherichia coli enzyme IIANtr regulates the K+ transporter TrkA. PNAS. 2007, 104: 4124-4129. 10.1073/pnas.0609897104.View ArticleGoogle Scholar
- 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. 10.1111/j.1365-2958.2009.06704.x.View ArticleGoogle Scholar
- Lee CR, et al.: Potassium mediates Escherichia coli enzyme IIA(Ntr) - dependent regulation of sigma factor selectivity. Mol Microbiol. 2010, 7: 1468-1483.View ArticleGoogle Scholar
- 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. 10.1128/JB.00325-11.View ArticleGoogle Scholar
- Ninfa A: Unnecessary signaling: Poorly named?. J Bacteriol. 2011, 193 (18): 4571-4573. 10.1128/JB.05682-11.View ArticleGoogle Scholar
- 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. 10.1074/jbc.270.9.4822.View ArticleGoogle Scholar
- Kremling A: Analysis of global control of Escherichia coli carbohydrate uptake. BMC Syst Biol. 1: 42-
- Pflüger K, de Lorenzo V: Growth-dependent phosphorylation of the PtsN (EIINtr) protein of Pseudomonas putida. J Biol Chem. 2007, 282: 18206-18211. 10.1074/jbc.M611110200.View ArticleGoogle Scholar
- 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-10.1371/journal.pcbi.1000210.View ArticleGoogle Scholar
- 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. 10.1016/j.mimet.2007.08.004.View ArticleGoogle Scholar
- 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. 10.1128/JB.02002-07.View ArticleGoogle Scholar
- Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli. Mol Syst Biol. 2007, 3: 119-Google Scholar
- Steffen T: Control reconfiguration of dynamical systems. 2005, SpringerGoogle Scholar
- 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. 10.1074/jbc.274.22.15562.View ArticleGoogle Scholar
- Cases I, de Lorenzo V: The black cat/white cat principle of signal integration in bacterial promoters. Embo J. 2001, 20: 1-11. 10.1093/emboj/20.1.1.View ArticleGoogle Scholar
- Miller JH: Experiments in Molecular Genetics. 1972, New York: Cold Spring Harbor LaboratoryGoogle Scholar
- 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.Google Scholar
- 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. 10.1128/JB.00033-07.View ArticleGoogle Scholar
- Ostergaard S: Increasing galactose consumption by Saccharomyces cerevisiae through metabolic engineering of the GAL gene regulatory network. Nat Biotechnol. 2000, 18: 1283-1286. 10.1038/82400.View ArticleGoogle Scholar
- Rapid quantifiable assessment of nutritional parameters influencing pediocin production by, Pediococcus acidilactici NRRL B5627 Bioresource technology. 2008, 99: 6646-6650. 10.1016/j.biortech.2007.11.068.Google Scholar
- Kylma AK, et al.: Microbial Biotech. 2004, 63: 584-591. 10.1007/s00253-003-1406-4.View ArticleGoogle Scholar
- Feist AM, Palsson BO: The growing scope of applications of genome-scale metabolic reconstructions using Escherichia coli. Nat Biotechnol. 2008, 26: 659-67. 10.1038/nbt1401.View ArticleGoogle Scholar
- Pflüger-Grau K, Görke B: Regulatory roles of the bacterial nitrogen-related phosphotransferase system. Trends Microbiol. 2010, 18: 205-214. 10.1016/j.tim.2010.02.003.View ArticleGoogle Scholar
- 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. 10.1016/j.bbagen.2011.01.002.View ArticleGoogle Scholar
- Sawyer MH, et al.: Pathways of D-fructose catabolism in species of Pseudomonas. Arch Microbiol. 1977, 112: 49-55. 10.1007/BF00446653.View ArticleGoogle Scholar
- Van Dijken JP, Quayle J: R: Fructose metabolism in four Pseudomonas species. Arch Microbiol. 1977, 114: 281-286. 10.1007/BF00446874.View ArticleGoogle Scholar
This article is published under license to 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.