Achieving optimal growth through product feedback inhibition in metabolism

Flux balance analysis of metabolic networks for efficient engineering of microbial cell factories

Metabolic engineering principles have long been applied to explore the metabolic networks of complex microbial cell factories under a variety of environmental constraints for effective deployment of the microorganisms in the optimal production of biochemicals like biofuels, polymers, amino acids, recombinant proteins. One of the methodologies used for analyzing microbial metabolic networks is the Flux Balance Analysis (FBA), which employs applications of optimization techniques for forecasting biomass growth and metabolic flux distribution of industrially important products under specified environmental conditions. The in silico flux simulations are instrumental for designing the production-specific microbial cell factories. However, FBA has some inherent limitations. The present review emphasizes how the incorporation of additional kinetic, thermodynamic, expression and regulatory constraints and integration of omics data into the classical FBA platform improve the prediction accuracy of FBA. A programmed comparison of the simulated data with the experimental observations is presented for supporting the claim. The review further accounts for the successful implementation of classical FBA in biotechnological applications and identifies areas in which classical FBA fails to make correct predictions. The analysis of the predictive capabilities of the different FBA strategies presented here is expected to help researchers in finding new avenues in engineering highly efficient microbial metabolic pathways and identify the key metabolic bottlenecks during the process. Based on the appropriate metabolic network design, fermentation engineers will be able to effectively design the bioreactors and optimize large-scale biochemical production through suitable pathway modifications.

How Does Escherichia coli Allocate Proteome?

Microorganisms are shown to actively partition their intracellular resources, such as proteins, for growth optimization. Recent experiments have begun to reveal molecular components unpinning the partition; however, quantitatively, it remains unclear how individual parts orchestrate to yield precise resource allocation that is both robust and dynamic. Here, we developed a coarse-grained mathematical framework that centers on guanosine pentaphosphate (ppGpp)-mediated regulation and used it to systematically uncover the design principles of proteome allocation in Escherichia coli. Our results showed that the cellular ability of resource partition lies in an ultrasensitive, negative feedback-controlling topology with the ultrasensitivity arising from zero-order amino acid kinetics and the negative feedback from ppGpp-controlled ribosome synthesis. In addition, together with the time-scale separation between slow ribosome kinetics and fast turnovers of ppGpp and amino acids, the network topology confers the organism an optimization mechanism that mimics sliding mode control, a nonlinear optimization strategy that is widely used in man-made systems. We further showed that such a controlling mechanism is robust against parameter variations and molecular fluctuations and is also efficient for biomass production over time. This work elucidates the fundamental controlling mechanism of E. coli proteome allocation, thereby providing insights into quantitative microbial physiology as well as the design of synthetic gene networks.

Spatially structured exchange of metabolites enhances bacterial survival and resilience in biofilms

Biofilm formation enhances bacterial survival and antibiotic tolerance, but the underlying mechanisms are incompletely understood. Here, we show that biofilm growth is accompanied by a reduction in bacterial energy metabolism and membrane potential, together with metabolic exchanges between the inner and outer regions in biofilms. More specifically, nutrient-starved cells in the interior supply amino acids to cells in the periphery, while peripheral cells experience a decrease in membrane potential and provide fatty acids to interior cells. Fatty acids facilitate the repair of starvation-induced membrane damage in inner cells and enhance their survival in the presence of antibiotics. Thus, metabolic exchanges between inner and outer cells contribute to survival of the nutrient-starved inner cells and contribute to antibiotic tolerance within the biofilm.

Revisiting kinetic Monte Carlo algorithms for time-dependent processes: From open-loop control to feedback control

Simulating stochastic systems with feedback control is challenging due to the complex interplay between the system's dynamics and the feedback-dependent control protocols. We present a single-step-trajectory probability analysis to time-dependent stochastic systems. Based on this analysis, we revisit several time-dependent kinetic Monte Carlo (KMC) algorithms designed for systems under open-loop-control protocols. Our analysis provides a unified alternative proof to these algorithms, summarized into a pedagogical tutorial. Moreover, with the trajectory probability analysis, we present a novel feedback-controlled KMC algorithm that accurately captures the dynamics systems controlled by an external signal based on the measurements of the system's state. Our method correctly captures the system dynamics and avoids the artificial Zeno effect that arises from incorrectly applying the direct Gillespie algorithm to feedback-controlled systems. This work provides a unified perspective on existing open-loop-control KMC algorithms and also offers a powerful and accurate tool for simulating stochastic systems with feedback control.

Metabolomics and Microbial Metabolism: Toward a Systematic Understanding

Over the past decades, our understanding of microbial metabolism has increased dramatically. Metabolomics, a family of techniques that are used to measure the quantities of small molecules in biological samples, has been central to these efforts. Advances in analytical chemistry have made it possible to measure the relative and absolute concentrations of more and more compounds with increasing levels of certainty. In this review, we highlight how metabolomics has contributed to understanding microbial metabolism and in what ways it can still be deployed to expand our systematic understanding of metabolism. To that end, we explain how metabolomics was used to (a) characterize network topologies of metabolism and its regulation networks, (b) elucidate the control of metabolic function, and (c) understand the molecular basis of higher-order phenomena. We also discuss areas of inquiry where technological advances should continue to increase the impact of metabolomics, as well as areas where our understanding is bottlenecked by other factors such as the availability of statistical and modeling frameworks that can extract biological meaning from metabolomics data.

Cell-cell metabolite exchange creates a pro-survival metabolic environment that extends lifespan

Metabolism is deeply intertwined with aging. Effects of metabolic interventions on aging have been explained with intracellular metabolism, growth control, and signaling. Studying chronological aging in yeast, we reveal a so far overlooked metabolic property that influences aging via the exchange of metabolites. We observed that metabolites exported by young cells are re-imported by chronologically aging cells, resulting in cross-generational metabolic interactions. Then, we used self-establishing metabolically cooperating communities (SeMeCo) as a tool to increase metabolite exchange and observed significant lifespan extensions. The longevity of the SeMeCo was attributable to metabolic reconfigurations in methionine consumer cells. These obtained a more glycolytic metabolism and increased the export of protective metabolites that in turn extended the lifespan of cells that supplied them with methionine. Our results establish metabolite exchange interactions as a determinant of cellular aging and show that metabolically cooperating cells can shape the metabolic environment to extend their lifespan.

Glutamate: A multifunctional amino acid in plants

Glutamate (Glu) is a versatile metabolite and a signaling molecule in plants. Glu biosynthesis is associated with the primary nitrogen assimilation pathway. The conversion between Glu and 2-oxoglutarate connects Glu metabolism to the tricarboxylic acid cycle, carbon metabolism, and energy production. Glu is the predominant amino donor for transamination reactions in the cell. In addition to protein synthesis, Glu is a building block for tetrapyrroles, glutathione, and folate. Glu is the precursor of γ-aminobutyric acid that plays an important role in balancing carbon/nitrogen metabolism and various cellular processes. Glu can conjugate to the major auxin indole 3-acetic acid (IAA), and IAA-Glu is destined for oxidative degradation. Glu also conjugates with isochorismate for the production of salicylic acid. Accumulating evidence indicates that Glu functions as a signaling molecule to regulate plant growth, development, and defense responses. The ligand-gated Glu receptor-like proteins (GLRs) mediate some of these responses. However, many of the Glu signaling events are GLR-independent. The receptor perceiving extracellular Glu as a danger signal is still unknown. In addition to GLRs, Glu may act on receptor-like kinases or receptor-like proteins to trigger immune responses. Glu metabolism and Glu signaling may entwine to regulate growth, development, and defense responses in plants.

Microbial communities form rich extracellular metabolomes that foster metabolic interactions and promote drug tolerance

Microbial communities are composed of cells of varying metabolic capacity, and regularly include auxotrophs that lack essential metabolic pathways. Through analysis of auxotrophs for amino acid biosynthesis pathways in microbiome data derived from >12,000 natural microbial communities obtained as part of the Earth Microbiome Project (EMP), and study of auxotrophic–prototrophic interactions in self-establishing metabolically cooperating yeast communities (SeMeCos), we reveal a metabolically imprinted mechanism that links the presence of auxotrophs to an increase in metabolic interactions and gains in antimicrobial drug tolerance. As a consequence of the metabolic adaptations necessary to uptake specific metabolites, auxotrophs obtain altered metabolic flux distributions, export more metabolites and, in this way, enrich community environments in metabolites. Moreover, increased efflux activities reduce intracellular drug concentrations, allowing cells to grow in the presence of drug levels above minimal inhibitory concentrations. For example, we show that the antifungal action of azoles is greatly diminished in yeast cells that uptake metabolites from a metabolically enriched environment. Our results hence provide a mechanism that explains why cells are more robust to drug exposure when they interact metabolically.

Using microbiome data analysis and a self-establishing metabolically cooperating yeast community model, the authors show that the presence of auxotrophs in a microbial community increases metabolic interactions between cells and fosters antimicrobial drug tolerance.

A simple regulatory architecture allows learning the statistical structure of a changing environment

Bacteria live in environments that are continuously fluctuating and changing. Exploiting any predictability of such fluctuations can lead to an increased fitness. On longer timescales, bacteria can 'learn' the structure of these fluctuations through evolution. However, on shorter timescales, inferring the statistics of the environment and acting upon this information would need to be accomplished by physiological mechanisms. Here, we use a model of metabolism to show that a simple generalization of a common regulatory motif (end-product inhibition) is sufficient both for learning continuous-valued features of the statistical structure of the environment and for translating this information into predictive behavior; moreover, it accomplishes these tasks near-optimally. We discuss plausible genetic circuits that could instantiate the mechanism we describe, including one similar to the architecture of two-component signaling, and argue that the key ingredients required for such predictive behavior are readily accessible to bacteria.

A unifying autocatalytic network-based framework for bacterial growth laws

Bacterial cells contain various autocatalytic cycles, e.g., the ribosome cycle, where ribosomes translate ribosomal proteins that subsequently self-assemble to form new ribosomes. Here, we show that the transcription–translation machinery couples all cellular autocatalytic cycles, resulting in balanced exponential growth. Each autocatalytic cycle generates two types of growth laws. We derive the RNA polymerase (RNAP) growth law based on the RNAP autocatalytic cycle, where RNAPs transcribe messenger RNAs (mRNAs) of its constituent Rpo protein subunits. Before degrading, these mRNAs catalyze Rpo proteins employing ribosomes. The Rpo proteins subsequently self-assemble, forming new RNAPs, thus completing the cycle. Contrary to ribosome growth law, a reduction in growth rate due to shortage in RNAPs occurs without affecting the ribosomal protein mass fraction.

Recently discovered simple quantitative relations, known as bacterial growth laws, hint at the existence of simple underlying principles at the heart of bacterial growth. In this work, we provide a unifying picture of how these known relations, as well as relations that we derive, stem from a universal autocatalytic network common to all bacteria, facilitating balanced exponential growth of individual cells. We show that the core of the cellular autocatalytic network is the transcription–translation machinery—in itself an autocatalytic network comprising several coupled autocatalytic cycles, including the ribosome, RNA polymerase, and transfer RNA (tRNA) charging cycles. We derive two types of growth laws per autocatalytic cycle, one relating growth rate to the relative fraction of the catalyst and its catalysis rate and the other relating growth rate to all the time scales in the cycle. The structure of the autocatalytic network generates numerous regimes in state space, determined by the limiting components, while the number of growth laws can be much smaller. We also derive a growth law that accounts for the RNA polymerase autocatalytic cycle, which we use to explain how growth rate depends on the inducible expression of the rpoB and rpoC genes, which code for the RpoB and C protein subunits of RNA polymerase, and how the concentration of rifampicin, which targets RNA polymerase, affects growth rate without changing the RNA-to-protein ratio. We derive growth laws for tRNA synthesis and charging and predict how growth rate depends on temperature, perturbation to ribosome assembly, and membrane synthesis.

Modeling microbial metabolic trade-offs in a chemostat

Microbes face intense competition in the natural world, and so need to wisely allocate their resources to multiple functions, in particular to metabolism. Understanding competition among metabolic strategies that are subject to trade-offs is therefore crucial for deeper insight into the competition, cooperation, and community assembly of microorganisms. In this work, we evaluate competing metabolic strategies within an ecological context by considering not only how the environment influences cell growth, but also how microbes shape their chemical environment. Utilizing chemostat-based resource-competition models, we exhibit a set of intuitive and general procedures for assessing metabolic strategies. Using this framework, we are able to relate and unify multiple metabolic models, and to demonstrate how the fitness landscape of strategies becomes intrinsically dynamic due to species-environment feedback. Such dynamic fitness landscapes produce rich behaviors, and prove to be crucial for ecological and evolutionarily stable coexistence in all the models we examined.

Allosteric Feedback Inhibition Enables Robust Amino Acid Biosynthesis in E. coli by Enforcing Enzyme Overabundance

Microbes must ensure robust amino acid metabolism in the face of external and internal perturbations. This robustness is thought to emerge from regulatory interactions in metabolic and genetic networks. Here, we explored the consequences of removing allosteric feedback inhibition in seven amino acid biosynthesis pathways in Escherichia coli (arginine, histidine, tryptophan, leucine, isoleucine, threonine, and proline). Proteome data revealed that enzyme levels decreased in five of the seven dysregulated pathways. Despite that, flux through the dysregulated pathways was not limited, indicating that enzyme levels are higher than absolutely needed in wild-type cells. We showed that such enzyme overabundance renders the arginine, histidine, and tryptophan pathways robust against perturbations of gene expression, using a metabolic model and CRISPR interference experiments. The results suggested a sensitive interaction between allosteric feedback inhibition and enzyme-level regulation that ensures robust yet efficient biosynthesis of histidine, arginine, and tryptophan in E. coli.

•

Amino acid biosynthesis enzymes do not normally operate at maximum capacity

•

Allosteric feedback inhibition ensures that enzymes are overabundant

•

Enzyme overabundance provides robustness against decreases in gene expression

Discovery and Functional Characterization of a Yeast Sugar Alcohol Phosphatase

Sugar alcohols (polyols) exist widely in nature. While some specific sugar alcohol phosphatases are known, there is no known phosphatase for some important sugar alcohols (e.g., sorbitol-6-phosphate). Using liquid chromatography-mass spectrometry-based metabolomics, we screened yeast strains with putative phosphatases of unknown function deleted. We show that the yeast gene YNL010W, which has close homologues in all fungi species and some plants, encodes a sugar alcohol phosphatase. We term this enzyme, which hydrolyzes sorbitol-6-phosphate, ribitol-5-phosphate, and (d)-glycerol-3-phosphate, polyol phosphatase 1 or PYP1. Polyol phosphates are structural analogs of the enediol intermediate of phosphoglucose isomerase (Pgi). We find that sorbitol-6-phosphate and ribitol-5-phosphate inhibit Pgi and that Pyp1 activity is important for yeast to maintain Pgi activity in the presence of environmental sugar alcohols. Pyp1 expression is strongly positively correlated with yeast growth rate, presumably because faster growth requires greater glycolytic and accordingly Pgi flux. Thus, yeast express the previously uncharacterized enzyme Pyp1 to prevent inhibition of glycolysis by sugar alcohol phosphates. Pyp1 may be useful for engineering sugar alcohol production.

Engineering Microbial Metabolite Dynamics and Heterogeneity

As yields for biological chemical production in microorganisms approach their theoretical maximum, metabolic engineering requires new tools, and approaches for improvements beyond what traditional strategies can achieve. Engineering metabolite dynamics and metabolite heterogeneity is necessary to achieve further improvements in product titers, productivities, and yields. Metabolite dynamics, the ensemble change in metabolite concentration over time, arise from the need for microbes to adapt their metabolism in response to the extracellular environment and are important for controlling growth and productivity in industrial fermentations. Metabolite heterogeneity, the cell-to-cell variation in a metabolite concentration in an isoclonal population, has a significant impact on ensemble productivity. Recent advances in single cell analysis enable a more complete understanding of the processes driving metabolite heterogeneity and reveal metabolic engineering targets. The authors present an overview of the mechanistic origins of metabolite dynamics and heterogeneity, why they are important, their potential effects in chemical production processes, and tools and strategies for engineering metabolite dynamics and heterogeneity. The authors emphasize that the ability to control metabolite dynamics and heterogeneity will bring new avenues of engineering to increase productivity of microbial strains.

Glucose becomes one of the worst carbon sources for E.coli on poor nitrogen sources due to suboptimal levels of cAMP

In most conditions, glucose is the best carbon source for E. coli: it provides faster growth than other sugars, and is consumed first in sugar mixtures. Here we identify conditions in which E. coli strains grow slower on glucose than on other sugars, namely when a single amino acid (arginine, glutamate, or proline) is the sole nitrogen source. In sugar mixtures with these nitrogen sources, E. coli still consumes glucose first, but grows faster rather than slower after exhausting glucose, generating a reversed diauxic shift. We trace this counterintuitive behavior to a metabolic imbalance: levels of TCA-cycle metabolites including α-ketoglutarate are high, and levels of the key regulatory molecule cAMP are low. Growth rates were increased by experimentally increasing cAMP levels, either by adding external cAMP, by genetically perturbing the cAMP circuit or by inhibition of glucose uptake. Thus, the cAMP control circuitry seems to have a ‘bug’ that leads to slow growth under what may be an environmentally rare condition.

Cellular Growth Arrest and Persistence from Enzyme Saturation

Metabolic efficiency depends on the balance between supply and demand of metabolites, which is sensitive to environmental and physiological fluctuations, or noise, causing shortages or surpluses in the metabolic pipeline. How cells can reliably optimize biomass production in the presence of metabolic fluctuations is a fundamental question that has not been fully answered. Here we use mathematical models to predict that enzyme saturation creates distinct regimes of cellular growth, including a phase of growth arrest resulting from toxicity of the metabolic process. Noise can drive entry of single cells into growth arrest while a fast-growing majority sustains the population. We confirmed these predictions by measuring the growth dynamics of Escherichia coli utilizing lactose as a sole carbon source. The predicted heterogeneous growth emerged at high lactose concentrations, and was associated with cell death and production of antibiotic-tolerant persister cells. These results suggest how metabolic networks may balance costs and benefits, with important implications for drug tolerance.

In bacteria, changes in gene expression, with resulting changes in protein concentration, can drastically change how fast cells and cellular populations grow. This fact has big implications for how we treat infectious disease, which types of organisms make up our microbiomes, and what patterns of gene regulation have undergone evolutionary selection. Here, we show how, in principle, the expression level of a single enzyme can affect bacterial population growth by creating a threshold where cells grow optimally fast just below it, but rapidly reach a state of no growth just above it because metabolic byproducts build up and halt growth. The narrow margin between these two states makes entering either of them possible for the same bacterium because of intrinsic uncertainty, or "noise", in gene expression. The predicted result is a variety of growth rates in a single population of genetically identical cells, manifested as a mix of fast- and slow-growing cells. We created laboratory conditions that reproduce the effect in the model organism E. coli, and showed that there may be a benefit to having slower growing cells, because they can survive antibiotic exposure for longer.

Real-time metabolome profiling of the metabolic switch between starvation and growth

Metabolic systems are often the first networks to respond to environmental changes, and the ability to monitor metabolite dynamics is key for understanding these cellular responses. Because monitoring metabolome changes is experimentally tedious and demanding, dynamic data on time scales from seconds to hours are scarce. Here we describe real-time metabolome profiling by direct injection of living bacteria, yeast or mammalian cells into a high-resolution mass spectrometer, which enables automated monitoring of about 300 compounds in 15-30-s cycles over several hours. We observed accumulation of energetically costly biomass metabolites in Escherichia coli in carbon starvation-induced stationary phase, as well as the rapid use of these metabolites upon growth resumption. By combining real-time metabolome profiling with modeling and inhibitor experiments, we obtained evidence for switch-like feedback inhibition in amino acid biosynthesis and for control of substrate availability through the preferential use of the metabolically cheaper one-step salvaging pathway over costly ten-step de novo purine biosynthesis during growth resumption.

Regulation of indole-3-acetic acid biosynthesis by branched-chain amino acids in Enterobacter cloacae UW5

The soil bacterium Enterobacter cloacae UW5 produces the rhizosphere signaling molecule indole-3-acetic acid (IAA) via the indolepyruvate pathway. Expression of indolepyruvate decarboxylase, a key pathway enzyme encoded by ipdC, is upregulated by the transcription factor TyrR in response to aromatic amino acids. Some members of the TyrR regulon may also be controlled by branched-chain amino acids and here we show that expression from the ipdC promoter and production of IAA are downregulated by valine, leucine and isoleucine. Regulation of the IAA synthesis pathway by both aromatic and branched-chain amino acids suggests a broader role for this pathway in bacterial physiology, beyond plant interactions.

Krebs cycle metabolon formation: metabolite concentration gradient enhanced compartmentation of sequential enzymes

Dynamics of metabolon formation in mitochondria was probed by studying diffusional motion of two sequential Krebs cycle enzymes in a microfluidic channel. Enhanced directional co-diffusion of both enzymes against a substrate concentration gradient was observed in the presence of intermediate generation. This reveals a metabolite directed compartmentation of metabolic pathways.

Computing the functional proteome: recent progress and future prospects for genome-scale models

Constraint-based models enable the computation of feasible, optimal, and realized biological phenotypes from reaction network reconstructions and constraints on their operation. To date, stoichiometric reconstructions have largely focused on metabolism, resulting in genome-scale metabolic models (M-Models). Recent expansions in network content to encompass proteome synthesis have resulted in models of metabolism and protein expression (ME-Models). ME-Models advance the predictions possible with constraint-based models from network flux states to the spatially resolved molecular composition of a cell. Specifically, ME-Models enable the prediction of transcriptome and proteome allocation and limitations, and basal expression states and regulatory needs. Continued expansion in reconstruction content and constraints will result in an increasingly refined representation of cellular composition and behavior.

Coordination of microbial metabolism

Beyond fuelling cellular activities with building blocks and energy, metabolism also integrates environmental conditions into intracellular signals. The underlying regulatory network is complex and multifaceted: it ranges from slow interactions, such as changing gene expression, to rapid ones, such as the modulation of protein activity via post-translational modification or the allosteric binding of small molecules. In this Review, we outline the coordination of common metabolic tasks, including nutrient uptake, central metabolism, the generation of energy, the supply of amino acids and protein synthesis. Increasingly, a set of key metabolites is recognized to control individual regulatory circuits, which carry out specific functions of information input and regulatory output. Such a modular view of microbial metabolism facilitates an intuitive understanding of the molecular mechanisms that underlie cellular decision making.

Optimal control of gene expression for fast proteome adaptation to environmental change

Bacterial populations growing in a changing world must adjust their proteome composition in response to alterations in the environment. Rapid proteome responses to growth medium changes are expected to increase the average growth rate and fitness value of these populations. Little is known about the dynamics of proteome change, e.g., whether bacteria use optimal strategies of gene expression for rapid proteome adjustments and if there are lower bounds to the time of proteome adaptation in response to growth medium changes. To begin answering these types of questions, we modeled growing bacteria as stoichiometrically coupled networks of metabolic pathways. These are balanced during steady-state growth in a constant environment but are initially unbalanced after rapid medium shifts due to a shortage of enzymes required at higher concentrations in the new environment. We identified an optimal strategy for rapid proteome adjustment in the absence of protein degradation and found a lower bound to the time of proteome adaptation after medium shifts. This minimal time is determined by the ratio between the Kullback-Leibler distance from the pre- to the postshift proteome and the postshift steady-state growth rate. The dynamics of optimally controlled proteome adaptation has a simple analytical solution. We used detailed numerical modeling to demonstrate that realistic bacterial control systems can emulate this optimal strategy for rapid proteome adaptation. Our results may provide a conceptual link between the physiology and population genetics of growing bacteria.

Nitrogen assimilation in Escherichia coli: putting molecular data into a systems perspective

We present a comprehensive overview of the hierarchical network of intracellular processes revolving around central nitrogen metabolism in Escherichia coli. The hierarchy intertwines transport, metabolism, signaling leading to posttranslational modification, and transcription. The protein components of the network include an ammonium transporter (AmtB), a glutamine transporter (GlnHPQ), two ammonium assimilation pathways (glutamine synthetase [GS]-glutamate synthase [glutamine 2-oxoglutarate amidotransferase {GOGAT}] and glutamate dehydrogenase [GDH]), the two bifunctional enzymes adenylyl transferase/adenylyl-removing enzyme (ATase) and uridylyl transferase/uridylyl-removing enzyme (UTase), the two trimeric signal transduction proteins (GlnB and GlnK), the two-component regulatory system composed of the histidine protein kinase nitrogen regulator II (NRII) and the response nitrogen regulator I (NRI), three global transcriptional regulators called nitrogen assimilation control (Nac) protein, leucine-responsive regulatory protein (Lrp), and cyclic AMP (cAMP) receptor protein (Crp), the glutaminases, and the nitrogen-phosphotransferase system. First, the structural and molecular knowledge on these proteins is reviewed. Thereafter, the activities of the components as they engage together in transport, metabolism, signal transduction, and transcription and their regulation are discussed. Next, old and new molecular data and physiological data are put into a common perspective on integral cellular functioning, especially with the aim of resolving counterintuitive or paradoxical processes featured in nitrogen assimilation. Finally, we articulate what still remains to be discovered and what general lessons can be learned from the vast amounts of data that are available now.

Transcriptional regulation is insufficient to explain substrate-induced flux changes in Bacillus subtilis

Regulation of enzyme expression is one key mechanism by which cells control their metabolic programs. In this work, a quantitative analysis of metabolism in a model bacterium under different conditions shows that expression alone cannot explain the majority of the observed metabolic changes.

Most enzymes are indeed highly expressed in conditions where they are more active.

Quantitatively, however, the observed changes in expression between conditions do not match the changes in activity for most enzymes.

A good quantitative match is only observed for enzymes involved in the TCA cycle.

Metabolomics reveals that increased substrate availability explains only a few instances of changes in activity.

One of the key ways in which microbes are thought to regulate their metabolism is by modulating the availability of enzymes through transcriptional regulation. However, the limited success of efforts to manipulate metabolic fluxes by rewiring the transcriptional network has cast doubt on the idea that transcript abundance controls metabolic fluxes. In this study, we investigate control of metabolic flux in the model bacterium Bacillus subtilis by quantifying fluxes, transcripts, and metabolites in eight metabolic states enforced by different environmental conditions. We find that most enzymes whose flux switches between on and off states, such as those involved in substrate uptake, exhibit large corresponding transcriptional changes. However, for the majority of enzymes in central metabolism, enzyme concentrations were insufficient to explain the observed fluxes—only for a number of reactions in the tricarboxylic acid cycle were enzyme changes approximately proportional to flux changes. Surprisingly, substrate changes revealed by metabolomics were also insufficient to explain observed fluxes, leaving a large role for allosteric regulation and enzyme modification in the control of metabolic fluxes.

Coordination of bacterial proteome with metabolism by cyclic AMP signalling

The cyclic AMP (cAMP)-dependent catabolite repression effect in Escherichia coli is among the most intensely studied regulatory processes in biology. However, the physiological function(s) of cAMP signalling and its molecular triggers remain elusive. Here we use a quantitative physiological approach to show that cAMP signalling tightly coordinates the expression of catabolic proteins with biosynthetic and ribosomal proteins, in accordance with the cellular metabolic needs during exponential growth. The expression of carbon catabolic genes increased linearly with decreasing growth rates upon limitation of carbon influx, but decreased linearly with decreasing growth rate upon limitation of nitrogen or sulphur influx. In contrast, the expression of biosynthetic genes showed the opposite linear growth-rate dependence as the catabolic genes. A coarse-grained mathematical model provides a quantitative framework for understanding and predicting gene expression responses to catabolic and anabolic limitations. A scheme of integral feedback control featuring the inhibition of cAMP signalling by metabolic precursors is proposed and validated. These results reveal a key physiological role of cAMP-dependent catabolite repression: to ensure that proteomic resources are spent on distinct metabolic sectors as needed in different nutrient environments. Our findings underscore the power of quantitative physiology in unravelling the underlying functions of complex molecular signalling networks.

Protein-DNA binding dynamics predict transcriptional response to nutrients in archaea

Organisms across all three domains of life use gene regulatory networks (GRNs) to integrate varied stimuli into coherent transcriptional responses to environmental pressures. However, inferring GRN topology and regulatory causality remains a central challenge in systems biology. Previous work characterized TrmB as a global metabolic transcription factor in archaeal extremophiles. However, it remains unclear how TrmB dynamically regulates its ∼100 metabolic enzyme-coding gene targets. Using a dynamic perturbation approach, we elucidate the topology of the TrmB metabolic GRN in the model archaeon Halobacterium salinarum. Clustering of dynamic gene expression patterns reveals that TrmB functions alone to regulate central metabolic enzyme-coding genes but cooperates with various regulators to control peripheral metabolic pathways. Using a dynamical model, we predict gene expression patterns for some TrmB-dependent promoters and infer secondary regulators for others. Our data suggest feed-forward gene regulatory topology for cobalamin biosynthesis. In contrast, purine biosynthesis appears to require TrmB-independent regulators. We conclude that TrmB is an important component for mediating metabolic modularity, integrating nutrient status and regulating gene expression dynamics alone and in concert with secondary regulators.

Catabolite regulation analysis of Escherichia coli for acetate overflow mechanism and co-consumption of multiple sugars based on systems biology approach using computer simulation

It is quite important to understand the basic principle embedded in the main metabolism for the interpretation of the fermentation data. For this, it may be useful to understand the regulation mechanism based on systems biology approach. In the present study, we considered the perturbation analysis together with computer simulation based on the models which include the effects of global regulators on the pathway activation for the main metabolism of Escherichia coli. Main focus is the acetate overflow metabolism and the co-fermentation of multiple carbon sources. The perturbation analysis was first made to understand the nature of the feed-forward loop formed by the activation of Pyk by FDP (F1,6BP), and the feed-back loop formed by the inhibition of Pfk by PEP in the glycolysis. Those together with the effect of transcription factor Cra caused by FDP level affected the glycolysis activity. The PTS (phosphotransferase system) acts as the feed-back system by repressing the glucose uptake rate for the increase in the glucose uptake rate. It was also shown that the increased PTS flux (or glucose consumption rate) causes PEP/PYR ratio to be decreased, and EIIA-P, Cya, cAMP-Crp decreased, where cAMP-Crp in turn repressed TCA cycle and more acetate is formed. This was further verified by the detailed computer simulation. In the case of multiple carbon sources such as glucose and xylose, it was shown that the sequential utilization of carbon sources was observed for wild type, while the co-consumption of multiple carbon sources with slow consumption rates were observed for the ptsG mutant by computer simulation, and this was verified by experiments. Moreover, the effect of a specific gene knockout such as Δpyk on the metabolic characteristics was also investigated based on the computer simulation.

Genome-scale genetic engineering in Escherichia coli

Genome engineering has been developed to create useful strains for biological studies and industrial uses. However, a continuous challenge remained in the field: technical limitations in high-throughput screening and precise manipulation of strains. Today, technical improvements have made genome engineering more rapid and efficient. This review introduces recent advances in genome engineering technologies applied to Escherichia coli as well as multiplex automated genome engineering (MAGE), a recent technique proposed as a powerful toolkit due to its straightforward process, rapid experimental procedures, and highly efficient properties.

Functional modules, structural topology, and optimal activity in metabolic networks

Modular organization in biological networks has been suggested as a natural mechanism by which a cell coordinates its metabolic strategies for evolving and responding to environmental perturbations. To understand how this occurs, there is a need for developing computational schemes that contribute to integration of genomic-scale information and assist investigators in formulating biological hypotheses in a quantitative and systematic fashion. In this work, we combined metabolome data and constraint-based modeling to elucidate the relationships among structural modules, functional organization, and the optimal metabolic phenotype of Rhizobium etli, a bacterium that fixes nitrogen in symbiosis with Phaseolus vulgaris. To experimentally characterize the metabolic phenotype of this microorganism, we obtained the metabolic profile of 220 metabolites at two physiological stages: under free-living conditions, and during nitrogen fixation with P. vulgaris. By integrating these data into a constraint-based model, we built a refined computational platform with the capability to survey the metabolic activity underlying nitrogen fixation in R. etli. Topological analysis of the metabolic reconstruction led us to identify modular structures with functional activities. Consistent with modular activity in metabolism, we found that most of the metabolites experimentally detected in each module simultaneously increased their relative abundances during nitrogen fixation. In this work, we explore the relationships among topology, biological function, and optimal activity in the metabolism of R. etli through an integrative analysis based on modeling and metabolome data. Our findings suggest that the metabolic activity during nitrogen fixation is supported by interacting structural modules that correlate with three functional classifications: nucleic acids, peptides, and lipids. More fundamentally, we supply evidence that such modular organization during functional nitrogen fixation is a robust property under different environmental conditions.

Biological networks are an inherent concept in systems biology that is useful in elucidating how biological entities—as metabolites or proteins—work together in supporting specific phenotypes in microorganisms. Notably, topological analyses carried out over these networks have shown that modular organization is a ubiquitous property at different levels of biological organization, in such a way that modular organization may serve as an organizing principle governing the metabolic activity in microorganisms. With the aim of elucidating the relationship among functional modules, network topology, and optimal metabolic activity, here we present an integrative study that combines computational modeling and metabolome data for evaluation of the metabolic activity of the soil bacterium Rhizobium etli during symbiotic nitrogen fixation with Phaseolus vulgaris. As a result, we supply experimental and computational evidence supporting the concept that the optimal metabolic activity during this biological process is guided by modular structures in the metabolic network of R. etli. Even more fundamentally, we suggest that these biochemical modules interact among each other to ensure an optimal phenotype during nitrogen fixation. Finally, through the in silico analysis on the genome scale metabolic reconstruction for R.etli, we give some examples that suggest that these modular structures supporting nitrogen fixation are robust to external physiological conditions.

Transcriptional changes underlying elemental stoichiometry shifts in a marine heterotrophic bacterium

Marine bacteria drive the biogeochemical processing of oceanic dissolved organic carbon (DOC), a 750-Tg C reservoir that is a critical component of the global C cycle. Catabolism of DOC is thought to be regulated by the biomass composition of heterotrophic bacteria, as cells maintain a C:N:P ratio of ∼50:10:1 during DOC processing. Yet a complicating factor in stoichiometry-based analyses is that bacteria can change the C:N:P ratio of their biomass in response to resource composition. We investigated the physiological mechanisms of resource-driven shifts in biomass stoichiometry in continuous cultures of the marine heterotrophic bacterium Ruegeria pomeroyi (a member of the Roseobacter clade) under four element limitation regimes (C, N, P, and S). Microarray analysis indicated that the bacterium scavenged for alternate sources of the scarce element when cells were C-, N-, or P-limited; reworked the ratios of biomolecules when C- and P- limited; and exerted tighter control over import/export and cytoplasmic pools when N-limited. Under S limitation, a scenario not existing naturally for surface ocean microbes, stress responses dominated transcriptional changes. Resource-driven changes in C:N ratios of up to 2.5-fold and in C:P ratios of up to sixfold were measured in R. pomeroyi biomass. These changes were best explained if the C and P content of the cells was flexible in the face of shifting resources but N content was not, achieved through the net balance of different transcriptional strategies. The cellular-level metabolic trade-offs that govern biomass stoichiometry in R. pomeroyi may have implications for global carbon cycling if extendable to other heterotrophic bacteria. Strong homeostatic responses to N limitation by marine bacteria would intensify competition with autotrophs. Modification of cellular inventories in C- and P-limited heterotrophs would vary the elemental ratio of particulate organic matter sequestered in the deep ocean.

α-Ketoglutarate coordinates carbon and nitrogen utilization via enzyme I inhibition

Microbes survive in a variety of nutrient environments by modulating their intracellular metabolism. Balanced growth requires coordinated uptake of carbon and nitrogen, the primary substrates for biomass production. The mechanisms that balance carbon and nitrogen uptake are, however, poorly understood. We find in Escherichia coli that a sudden increase in nitrogen availability results in an almost immediate increase in glucose uptake. The concentrations of known glycolytic intermediates and regulators, however, remain homeostatic. Instead, we find that α-ketoglutarate, which accumulates in nitrogen limitation, directly blocks glucose uptake by inhibiting Enzyme I, the first step of the phosphotransferase system (PTS). This enables rapid modulation of glycolytic flux without marked concentration changes in glycolytic intermediates by simultaneously accelerating glucose import and consumption of the terminal glycolytic intermediate phosphoenolpyruvate. Quantitative modeling shows that this previously unidentified regulatory connection is in principle sufficient to coordinate carbon and nitrogen utilization.

Modeling metabolic adjustment in Mycobacterium tuberculosis upon treatment with isoniazid

Complex biological systems exhibit a property of robustness at all levels of organization. Through different mechanisms, the system tries to sustain stress such as due to starvation or drug exposure. To explore whether reconfiguration of the metabolic networks is used as a means to achieve robustness, we have studied possible metabolic adjustments in Mtb upon exposure to isoniazid (INH), a front-line clinical drug. The redundancy in the genome of M. tuberculosis (Mtb) makes it an attractive system to explore if alternate routes of metabolism exist in the bacterium. While the mechanism of action of INH is well studied, its effect on the overall metabolism is not well characterized. Using flux balance analysis, inhibiting the fluxes flowing through the reactions catalyzed by Rv1484, the target of INH, significantly changes the overall flux profiles. At the pathway level, activation or inactivation of certain pathways distant from the target pathway, are seen. Metabolites such as NADPH are shown to reduce drastically, while fatty acids tend to accumulate. The overall biomass also decreases with increasing inhibition levels. Inhibition studies, pathway level clustering and comparison of the flux profiles with the gene expression data indicate the activation of folate metabolism, ubiquinone metabolism, and metabolism of certain amino acids. This analysis provides insights useful for target identification and designing strategies for combination therapy. Insights gained about the role of individual components of a system and their interactions will also provide a basis for reconstruction of whole systems through synthetic biology approaches.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (doi:10.1007/s11693-011-9075-6) contains supplementary material, which is available to authorized users.

Dispensability of Escherichia coli's latent pathways

Gene-knockout experiments on single-cell organisms have established that expression of a substantial fraction of genes is not needed for optimal growth. This problem acquired a new dimension with the recent discovery that environmental and genetic perturbations of the bacterium Escherichia coli are followed by the temporary activation of a large number of latent metabolic pathways, which suggests the hypothesis that temporarily activated reactions impact growth and hence facilitate adaptation in the presence of perturbations. Here, we test this hypothesis computationally and find, surprisingly, that the availability of latent pathways consistently offers no growth advantage and tends, in fact, to inhibit growth after genetic perturbations. This is shown to be true even for latent pathways with a known function in alternate conditions, thus extending the significance of this adverse effect beyond apparently nonessential genes. These findings raise the possibility that latent pathway activation is in fact derivative of another, potentially suboptimal, adaptive response.

Exploring and exploiting allostery: Models, evolution, and drug targeting

The concept of allostery was elaborated almost 50years ago by Monod and coworkers to provide a framework for interpreting experimental studies on the regulation of protein function. In essence, binding of a ligand at an allosteric site affects the function at a distant site exploiting protein flexibility and reshaping protein energy landscape. Both monomeric and oligomeric proteins can be allosteric. In the past decades, the behavior of allosteric systems has been analyzed in many investigations while general theoretical models and variations thereof have been steadily proposed to interpret the experimental data. Allostery has been established as a fundamental mechanism of regulation in all organisms, governing a variety of processes that range from metabolic control to receptor function and from ligand transport to cell motility. A number of studies have shed light on how evolutionary pressures have favored and molded the development of allosteric features in specific macromolecular systems. The widespread occurrence of allostery has been recently exploited for the development and design of allosteric drugs that bind to either physiological or non-physiological allosteric sites leading to gain of function or loss of function. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.