Modeling and experimental validation of the signal transduction via the Escherichia coli sucrose phospho transferase system

Self-replenishment cycles generate a threshold response

Many metabolic cycles, including the tricarboxylic acid cycle, glyoxylate cycle, Calvin cycle, urea cycle, coenzyme recycling, and substrate cycles, are well known to catabolize and anabolize different metabolites for efficient energy and mass conversion. In terms of stoichiometric structure, this study explicitly identifies two types of metabolic cycles. One is the well-known, elementary cycle that converts multiple substrates into different products and recycles one of the products as a substrate, where the recycled substrate is supplied from the outside to run the cycle. The other is the self-replenishment cycle that merges multiple substrates into two or multiple identical products and reuses one of the products as a substrate. The substrates are autonomously supplied within the cycle. This study first defines the self-replenishment cycles that many scientists have overlooked despite its functional importance. Theoretical analysis has revealed the design principle of the self-replenishment cycle that presents a threshold response without any bistability nor cooperativity. To verify the principle, three detailed kinetic models of self-replenishment cycles embedded in an E. coli metabolic system were simulated. They presented the threshold response or digital switch-like function that steeply shift metabolic status.

Sugar Influx Sensing by the Phosphotransferase System of Escherichia coli

The phosphotransferase system (PTS) plays a pivotal role in the uptake of multiple sugars in Escherichia coli and many other bacteria. In the cell, individual sugar-specific PTS branches are interconnected through a series of phosphotransfer reactions, thus creating a global network that not only phosphorylates incoming sugars but also regulates a number of cellular processes. Despite the apparent importance of the PTS network in bacterial physiology, the holistic function of the network in the cell remains unclear. Here we used Förster resonance energy transfer (FRET) to investigate the PTS network in E. coli, including the dynamics of protein interactions and the processing of different stimuli and their transmission to the chemotaxis pathway. Our results demonstrate that despite the seeming complexity of the cellular PTS network, its core part operates in a strikingly simple way, sensing the overall influx of PTS sugars irrespective of the sugar identity and distributing this information equally through all studied branches of the network. Moreover, it also integrates several other specific metabolic inputs. The integrated output of the PTS network is then transmitted linearly to the chemotaxis pathway, in stark contrast to the amplification of conventional chemotactic stimuli. Finally, we observe that default uptake through the uninduced PTS network correlates well with the quality of the carbon source, apparently representing an optimal regulatory strategy.

The bacterial phosphotransferase system (PTS) mediates uptake of multiple sugars from the environment and also controls cell physiology and swimming behavior in sugar gradients. In Escherichia coli and other bacteria, the PTS consists of a number of sugar-specific branches, interconnected via shared components through a series of phosphotransfer reactions. Whereas most previous studies have focused on understanding individual PTS branches, the holistic function of the entire PTS network in the cell remained elusive. In this study we address this question by investigating the dynamics of multiple protein interactions within the cellular PTS network upon stimulation with sugars and other metabolites. We demonstrate that despite its seeming complexity, the core part of the PTS network operates in a strikingly simple way, sensing the overall influx of PTS sugars and key metabolites into the cell and utilizing this information to control bacterial behavior. We further show that the default influx of the carbon source correlates with its quality, and we use computer simulations to demonstrate that this correlation apparently represents an optimal regulatory strategy.

Single-cell characterization of metabolic switching in the sugar phosphotransferase system of Escherichia coli

The utilization of several sugars in Escherichia coli is regulated by the Phosphotransferase System (PTS), in which diverse sugar utilization modules compete for phosphoryl flux from the general PTS proteins. Existing theoretical work predicts a winner-take-all outcome when this flux limits carbon uptake. To date, no experimental work has interrogated competing PTS uptake modules with single-cell resolution. Using time-lapse microscopy in perfused microchannels, we analyzed the competition between N-acetyl-glucosamine and sorbitol, as representative PTS sugars, by measuring both the expression of their utilization systems and the concomitant impact of sugar utilization on growth rates. We find two distinct regimes: hierarchical usage of the carbohydrates, and co-expression of the genes for both systems. Simulations of a mathematical model incorporating asymmetric sugar quality reproduce our metabolic phase diagram, indicating that under conditions of nonlimiting phosphate flux, co-expression is due to uncoupling of both sugar utilization systems. Our model reproduces hierarchical winner-take-all behaviour and stochastic co-expression, and predicts the switching between both strategies as a function of available phosphate flux. Hence, experiments and theory both suggest that PTS sugar utilization involves not only switching between the sugars utilized but also switching of utilization strategies to accommodate prevailing environmental conditions.

Understanding carbon catabolite repression in Escherichia coli using quantitative models

Carbon catabolite repression (CCR) controls the order in which different carbon sources are metabolized. Although this system is one of the paradigms of the regulation of gene expression in bacteria, the underlying mechanisms remain controversial. CCR involves the coordination of different subsystems of the cell that are responsible for the uptake of carbon sources, their breakdown for the production of energy and precursors, and the conversion of the latter to biomass. The complexity of this integrated system, with regulatory mechanisms cutting across metabolism, gene expression, and signaling, and that are subject to global physical and physiological constraints, has motivated important modeling efforts over the past four decades, especially in the enterobacterium Escherichia coli. Different hypotheses concerning the dynamic functioning of the system have been explored by a variety of modeling approaches. We review these studies and summarize their contributions to the quantitative understanding of CCR, focusing on diauxic growth in E. coli. Moreover, we propose a highly simplified representation of diauxic growth that makes it possible to bring out the salient features of the models proposed in the literature and confront and compare the explanations they provide.

EI of the Phosphotransferase System of Escherichia coli: Mathematical Modeling Approach to Analysis of Its Kinetic Properties

The mathematical model of the operation of the first enzyme of the Escherichia coli phosphotransferase system, EI, is proposed. Parameters of the kinetic model describing the operation of EI under different conditions are identified on the basis of a large amount of known experimental data. The verified model is employed to predict modes of operation of EI under both in vivo physiological conditions and in vitro nonphysiological conditions. The model predicts that under in vivo physiological conditions, the rate of phosphotransfer from EI to the second protein of the phosphotransferase system HPr by the dimer is much higher than by the monomer. A hypothesis is proposed on the basis of calculations that the transfer by a monomer plays a role in the regulation of chemotaxis. At submicromolar pyruvate concentration, the model predicts nonmonotonic dependence of the phosphotransfer rate on the substrate (PEP) concentration.

Influence of uncertainties in pH, pMg, activity coefficients, metabolite concentrations, and other factors on the analysis of the thermodynamic feasibility of metabolic pathways

Thermodynamic feasibility analysis (TFA) has been used as a tool capable of providing additional constraints to the mass balance-based methods of analysis of metabolic networks (e.g., flux balance analysis). Several publications have recently appeared in which TFA of different metabolic pathways from relatively simple to the genome-scale networks was described as a means of detecting the possible metabolic control steps. However, in order to perform TFA, many simplifying assumptions were necessary. On the other hand, it has been shown by applying TFA to the well-known pathway of glycolysis that erroneous simplifying assumptions may seriously bias the results of the analysis. A quantitative analysis of the influence of non-ideality of the biochemical system, pH, temperature, and complexation of the metabolites with Mg(2+) ions as well as a number of other factors on the TFA is reported. It is shown that the feasibility of glycolysis is very seriously limited by the reaction of oxidative phosphorylation of glyceraldehyde phosphate, and that the intracellular concentration of the main product of this reaction, biphosphoglycerate, must be anywhere from 10 to 100 times lower than published values. In addition, the driving force for this reaction, and consequently the feasibility of the entire pathway depend strongly on the intracellular pH and ionic strength and to a lesser extent on pMg and temperature. The analysis may also be influenced by uncertainties of the dissociation and magnesium complexation constants of glyceraldehyde phosphate. The analysis demonstrates the crucial importance of taking such factors into account when performing TFA. It also suggests an urgent need for experimental determinations of such factors as a prerequisite for sensible thermodynamic analysis of metabolism on a genome-wide scale.

HPrK regulates succinate-mediated catabolite repression in the gram-negative symbiont Sinorhizobium meliloti

The HPrK kinase/phosphatase is a common component of the phosphotransferase system (PTS) of gram-positive bacteria and regulates catabolite repression through phosphorylation/dephosphorylation of its substrate, the PTS protein HPr, at a conserved serine residue. Phosphorylation of HPr by HPrK also affects additional phosphorylation of HPr by the PTS enzyme EI at a conserved histidine residue. Sinorhizobium meliloti can live as symbionts inside legume root nodules or as free-living organisms and is one of the relatively rare gram-negative bacteria known to have a gene encoding HPrK. We have constructed S. meliloti mutants that lack HPrK or that lack key amino acids in HPr that are likely phosphorylated by HPrK and EI. Deletion of hprK in S. meliloti enhanced catabolite repression caused by succinate, as did an S53A substitution in HPr. Introduction of an H22A substitution into HPr alleviated the strong catabolite repression phenotypes of strains carrying Delta hprK or hpr(S53A) mutations, demonstrating that HPr-His22-P is needed for strong catabolite repression. Furthermore, strains with a hpr(H22A) allele exhibited relaxed catabolite repression. These results suggest that HPrK phosphorylates HPr at the serine-53 residue, that HPr-Ser53-P inhibits phosphorylation at the histidine-22 residue, and that HPr-His22-P enhances catabolite repression in the presence of succinate. Additional experiments show that Delta hprK mutants overproduce exopolysaccharides and form nodules that do not fix nitrogen.

Computer-aided rational design of the phosphotransferase system for enhanced glucose uptake in Escherichia coli

The phosphotransferase system (PTS) is the sugar transportation machinery that is widely distributed in prokaryotes and is critical for enhanced production of useful metabolites. To increase the glucose uptake rate, we propose a rational strategy for designing the molecular architecture of the Escherichia coli glucose PTS by using a computer-aided design (CAD) system and verified the simulated results with biological experiments. CAD supports construction of a biochemical map, mathematical modeling, simulation, and system analysis. Assuming that the PTS aims at controlling the glucose uptake rate, the PTS was decomposed into hierarchical modules, functional and flux modules, and the effect of changes in gene expression on the glucose uptake rate was simulated to make a rational strategy of how the gene regulatory network is engineered. Such design and analysis predicted that the mlc knockout mutant with ptsI gene overexpression would greatly increase the specific glucose uptake rate. By using biological experiments, we validated the prediction and the presented strategy, thereby enhancing the specific glucose uptake rate.

Correlation between growth rates, EIIACrr phosphorylation, and intracellular cyclic AMP levels in Escherichia coli K-12

In Escherichia coli K-12, components of the phosphoenolpyruvate-dependent phosphotransferase systems (PTSs) represent a signal transduction system involved in the global control of carbon catabolism through inducer exclusion mediated by phosphoenolpyruvate-dependent protein kinase enzyme IIA(Crr) (EIIA(Crr)) (= EIIA(Glc)) and catabolite repression mediated by the global regulator cyclic AMP (cAMP)-cAMP receptor protein (CRP). We measured in a systematic way the relation between cellular growth rates and the key parameters of catabolite repression, i.e., the phosphorylated EIIA(Crr) (EIIA(Crr) approximately P) level and the cAMP level, using in vitro and in vivo assays. Different growth rates were obtained by using either various carbon sources or by growing the cells with limited concentrations of glucose, sucrose, and mannitol in continuous bioreactor experiments. The ratio of EIIA(Crr) to EIIA(Crr) approximately P and the intracellular cAMP concentrations, deduced from the activity of a cAMP-CRP-dependent promoter, correlated well with specific growth rates between 0.3 h(-1) and 0.7 h(-1), corresponding to generation times of about 138 and 60 min, respectively. Below and above this range, these parameters were increasingly uncoupled from the growth rate, which perhaps indicates an increasing role executed by other global control systems, in particular the stringent-relaxed response system.

Analysis of global control of Escherichia coli carbohydrate uptake

Background

Global control influences the regulation of many individual subsystems by superimposed regulator proteins. A prominent example is the control of carbohydrate uptake systems by the transcription factor Crp in Escherichia coli. A detailed understanding of the coordination of the control of individual transporters offers possibilities to explore the potential of microorganisms e.g. in biotechnology.

Results

An o.d.e. based mathematical model is presented that maps a physiological parameter – the specific growth rate – to the sensor of the signal transduction unit, here a component of the bacterial phosphotransferase system (PTS), namely EIIA^Crr. The model describes the relation between the growth rate and the degree of phosphorylation of EIIA ^crr for a number of carbohydrates by a distinctive response curve, that differentiates between PTS transported carbohydrates and non-PTS carbohydrates. With only a small number of kinetic parameters, the model is able to describe a broad range of experimental steady-state and dynamical conditions.

Conclusion

The steady-state characteristic presented shows a relationship between the growth rate and the output of the sensor system PTS. The glycolytic flux that is measured by this sensor is a good indicator to represent the nutritional status of the cell.

How phosphotransferase system-related protein phosphorylation regulates carbohydrate metabolism in bacteria

The phosphoenolpyruvate(PEP):carbohydrate phosphotransferase system (PTS) is found only in bacteria, where it catalyzes the transport and phosphorylation of numerous monosaccharides, disaccharides, amino sugars, polyols, and other sugar derivatives. To carry out its catalytic function in sugar transport and phosphorylation, the PTS uses PEP as an energy source and phosphoryl donor. The phosphoryl group of PEP is usually transferred via four distinct proteins (domains) to the transported sugar bound to the respective membrane component(s) (EIIC and EIID) of the PTS. The organization of the PTS as a four-step phosphoryl transfer system, in which all P derivatives exhibit similar energy (phosphorylation occurs at histidyl or cysteyl residues), is surprising, as a single protein (or domain) coupling energy transfer and sugar phosphorylation would be sufficient for PTS function. A possible explanation for the complexity of the PTS was provided by the discovery that the PTS also carries out numerous regulatory functions. Depending on their phosphorylation state, the four proteins (domains) forming the PTS phosphorylation cascade (EI, HPr, EIIA, and EIIB) can phosphorylate or interact with numerous non-PTS proteins and thereby regulate their activity. In addition, in certain bacteria, one of the PTS components (HPr) is phosphorylated by ATP at a seryl residue, which increases the complexity of PTS-mediated regulation. In this review, we try to summarize the known protein phosphorylation-related regulatory functions of the PTS. As we shall see, the PTS regulation network not only controls carbohydrate uptake and metabolism but also interferes with the utilization of nitrogen and phosphorus and the virulence of certain pathogens.

How reliable are thermodynamic feasibility statements of biochemical pathways?

The driving force for organo- or lithotrophic growth as well as for each step in the metabolic network is the Gibbs reaction energy. For each enzymatic step it must be negative. Thermodynamics contributes therefore to the in-silico description of living systems. It may be used for assessing the feasibility of a given pathway because it provides a further constraint for those pathways which are feasible from the point of view of mass balance calculations (metabolic flux analysis) and the genetic potential of an organism. However, when this constraint was applied to lactic acid fermentation according to a method proposed by Mavrovouniotis (1993a, ISMB 93:273-283) it turned out that an unrealistically wide metabolite concentration range had to be assumed to make this well-known glycolytic pathway thermodynamically feasible. During a search for the reasons of this surprising result the insufficient consideration of the activity coefficients was identified as main cause. However, it is shown in the present contribution that the influence of the activity coefficients on Gibbs reaction energy can be easily taken into account based on the intracellular ionic strength. The uncertainty of the tabulated equilibrium constants and of the apparent standard Gibbs energies derived from them was found to be the second most important reason for the erroneous result of the feasibility analysis. Deviations of intracellular pH from the standard value and bad estimations of currency metabolites, e.g., NAD(+) and NADH, were found to be of lesser importance but not negligible. The pH dependency of Gibbs reaction enthalpy was proved to be easily taken into account. Therefore, the application of thermodynamics for a better in-silico prediction of the behavior of living cell factories calls predominantly for better equilibrium data determined under well defined conditions and also for a more detailed knowledge about the intracellular ionic strength and pH value.

A quantitative approach to catabolite repression in Escherichia coli

A dynamic mathematical model was developed to describe the uptake of various carbohydrates (glucose, lactose, glycerol, sucrose, and galactose) in Escherichia coli. For validation a number of isogenic strains with defined mutations were used. By considering metabolic reactions as well as signal transduction processes influencing the relevant pathways, we were able to describe quantitatively the phenomenon of catabolite repression in E. coli. We verified model predictions by measuring time courses of several extra- and intracellular components such as glycolytic intermediates, EII-ACrr phosphorylation level, both LacZ and PtsG concentrations, and total cAMP concentrations under various growth conditions. The entire data base consists of 18 experiments performed with nine different strains. The model describes the expression of 17 key enzymes, 38 enzymatic reactions, and the dynamic behavior of more than 50 metabolites. The different phenomena affecting the phosphorylation level of EIIACrr, the key regulation molecule for inducer exclusion and catabolite repression in enteric bacteria, can now be explained quantitatively.