A generalized computational framework to streamline thermodynamics and kinetics analysis of metabolic pathways

Metabolic engineering is a critical biotechnological approach in addressing global energy and environment challenges. Most engineering efforts, however, consist of laborious and inefficient trial-and-error of target pathways, due in part to the lack of methodologies that can comprehensively assess pathway properties in thermodynamics and kinetics. Metabolic engineering can benefit from computational tools that evaluate feasibility, expense and stability of non-natural metabolic pathways. Such tools can also help us understand natural pathways and their regulation at systems level. Here we introduce a computational toolbox, PathParser, which, for the first time, integrates multiple important functions for pathway analysis including thermodynamics analysis, kinetics-based protein cost optimization and robustness analysis. Specifically, PathParser enables optimization of the driving force of a pathway by minimizing the Gibbs free energy of least thermodynamically favorable reaction. In addition, based on reaction thermodynamics and enzyme kinetics, it can compute the minimal enzyme protein cost that supports metabolic flux, and evaluate pathway stability and flux in response to enzyme concentration perturbations. In a demo analysis of the Calvin–Benson–Bassham cycle and photorespiration pathway in the model cyanobacterium Synechocystis PCC 6803, the computation results are corroborated by experimental proteomics data as well as metabolic engineering outcomes. This toolbox may have broad application in metabolic engineering and systems biology in other microbial systems.     

Dynamic control in metabolic engineering: Theories,tools, and applications

Metabolic engineering has allowed the production of a diverse number of valuable chemicals using microbial organisms. Many biological challenges for improving bio-production exist which limit performance and slow the commercialization of metabolically engineered systems. Dynamic metabolic engineering is a rapidly developing field that seeks to address these challenges through the design of genetically encoded metabolic control systems which allow cells to autonomously adjust their flux in response to their external and internal metabolic state. This review first discusses theoretical works which provide mechanistic insights and design choices for dynamic control systems including two-stage, continuous, and population behavior control strategies. Next, we summarize molecular mechanisms for various sensors and actuators which enable dynamic metabolic control in microbial systems. Finally, important applications of dynamic control to the production of several metabolite products are highlighted, including fatty acids, aromatics, and terpene compounds. Altogether, this review provides a comprehensive overview of the progress, advances, and prospects in the design of dynamic control systems for improved titer, rate, and yield metrics in metabolic engineering.     

Monte-Carlo Modeling of the Central Carbon Metabolism of Lactococcus lactis: Insights into Metabolic Regulation

Metabolic pathways are complex dynamic systems whose response to perturbations and environmental challenges are governed by multiple interdependencies between enzyme properties, reactions rates, and substrate levels. Understanding the dynamics arising from such a network can be greatly enhanced by the construction of a computational model that embodies the properties of the respective system. Such models aim to incorporate mechanistic details of cellular interactions to mimic the temporal behavior of the biochemical reaction system and usually require substantial knowledge of kinetic parameters to allow meaningful conclusions. Several approaches have been suggested to overcome the severe data requirements of kinetic modeling, including the use of approximative kinetics and Monte-Carlo sampling of reaction parameters. In this work, we employ a probabilistic approach to study the response of a complex metabolic system, the central metabolism of the lactic acid bacterium Lactococcus lactis, subject to perturbations and brief periods of starvation. Supplementing existing methodologies, we show that it is possible to acquire a detailed understanding of the control properties of a corresponding metabolic pathway model that is directly based on experimental observations. In particular, we delineate the role of enzymatic regulation to maintain metabolic stability and metabolic recovery after periods of starvation. It is shown that the feedforward activation of the pyruvate kinase by fructose-1,6-bisphosphate qualitatively alters the bifurcation structure of the corresponding pathway model, indicating a crucial role of enzymatic regulation to prevent metabolic collapse for low external concentrations of glucose. We argue that similar probabilistic methodologies will help our understanding of dynamic properties of small-, medium- and large-scale metabolic networks models.     

Metabolic engineering of Bacillus subtilis fueled by systems biology: Recent advances and future directions

By combining advanced omics technology and computational modeling, systems biologists have identified and inferred thousands of regulatory events and system-wide interactions of the bacterium Bacillus subtilis, which is commonly used both in the laboratory and in industry. This dissection of the multiple layers of regulatory networks and their interactions has provided invaluable information for unraveling regulatory mechanisms and guiding metabolic engineering. In this review, we discuss recent advances in the systems biology and metabolic engineering of B. subtilis and highlight current gaps in our understanding of global metabolism and global pathway engineering in this organism. We also propose future perspectives in the systems biology of B. subtilis and suggest ways that this approach can be used to guide metabolic engineering. Specifically, although hundreds of regulatory events have been identified or inferred via systems biology approaches, systematic investigation of the functionality of these events in vivo has lagged, thereby preventing the elucidation of regulatory mechanisms and further rational pathway engineering. In metabolic engineering, ignoring the engineering of multilayer regulation hinders metabolic flux redistribution. Post-translational engineering, allosteric engineering, and dynamic pathway analyses and control will also contribute to the modulation and control of the metabolism of engineered B. subtilis, ultimately producing the desired cellular traits. We hope this review will aid metabolic engineers in making full use of available systems biology datasets and approaches for the design and perfection of microbial cell factories through global metabolism optimization.     

Trace elements in diabetes mellitus. Peculiarities and clinical validity of determinations in blood cells

Significantly more information about trace element status can be obtained by investigating concentrations in blood cells instead of only evaluating the concentrations in plasma. This can be explained by the fact that essential trace elements such as zinc, copper, chromium and selenium take part in a variety of enzymatic processes on a molecular cellular level. Ignoring these important biochemical roles, trace element concentrations determined in whole blood or plasma very often lead to conclusions contrary to the actual intracellular concentration. Especially in metabolic diseases like diabetes mellitus, conclusions drawn from trace element concentrations in blood cells usually offer more valuable clinical information about the metabolic state than trace element concentrations in plasma or whole blood. In the present investigation copper and zinc concentrations were increased in all blood fractions of diabetic patients (IDDM). In insulin-dependent diabetic children significantly higher values of zinc in erythrocytes were also found, and they were higher in patients with poor metabolic control (HbA1c>9%). When different blood fractions in diabetic patients (NIDDM) were compared with a control group, chromium was significantly increased in plasma and polymorphonuclear cells. Patients with IDDM had pronounced decreased selenium concentrations in erythrocytes as compared to controls.     

Metabolic engineering of Escherichia coli for the production of benzoic acid from glucose

Benzoic acid (BA) is an important platform aromatic compound in chemical industry and is widely used as food preservatives in its salt forms. Yet, current manufacture of BA is dependent on petrochemical processes under harsh conditions. Here we report the de novo production of BA from glucose using metabolically engineered Escherichia coli strains harboring a plant-like β-oxidation pathway or a newly designed synthetic pathway. First, three different natural BA biosynthetic pathways originated from plants and one synthetically designed pathway were systemically assessed for BA production from glucose by in silico flux response analyses. The selected plant-like β-oxidation pathway and the synthetic pathway were separately established in E. coli by expressing the genes encoding the necessary enzymes and screened heterologous enzymes under optimal plasmid configurations. BA production was further optimized by applying several metabolic engineering strategies to the engineered E. coli strains harboring each metabolic pathway, which included enhancement of the precursor availability, removal of competitive reactions, transporter engineering, and reduction of byproduct formation. Lastly, fed-batch fermentations of the final engineered strain harboring the β-oxidation pathway and the strain harboring the synthetic pathway were conducted, which resulted in the production of 2.37 ± 0.02 g/L and 181.0 ± 5.8 mg/L of BA from glucose, respectively; the former being the highest titer reported by microbial fermentation. The metabolic engineering strategies developed here will be useful for the production of related aromatics of high industrial interest.     

Development of a growth coupled and multi-layered dynamic regulation network balancing malonyl-CoA node to enhance (2S)-naringenin biosynthesis in Escherichia coli

Metabolic heterogeneity and dynamic changes in metabolic fluxes are two inherent characteristics of microbial fermentation that limit the precise control of metabolisms, often leading to impaired cell growth and low productivity. Dynamic metabolic engineering addresses these challenges through the design of multi-layered and multi-genetic dynamic regulation network (DRN) that allow a single cell to autonomously adjust metabolic flux in response to its growth and metabolite accumulation conditions. Here, we developed a growth coupled NCOMB (Naringenin-Coumaric acid-Malonyl-CoA-Balanced) DRN with systematic optimization of (2S)-naringenin and p-coumaric acid-responsive regulation pathways for real-time control of intracellular supply of malonyl-CoA. In this scenario, the acyl carrier protein was used as a novel critical node for fine-tuning malonyl-CoA consumption instead of direct repression of fatty acid synthase commonly employed in previous studies. To do so, we first engineered a multi-layered DRN enabling single cells to concurrently regulate acpH, acpS, acpT, acs, and ACC in malonyl-CoA catabolic and anabolic pathways. Next, the NCOMB DRN was optimized to enhance the synergies between different dynamic regulation layers via a biosensor-based directed evolution strategy. Finally, a high producer obtained from NCOMB DRN approach yielded a 8.7-fold improvement in (2S)-naringenin production (523.7 ± 51.8 mg/L) with a concomitant 20% increase in cell growth compared to the base strain using static strain engineering approach, thus demonstrating the high efficiency of this system for improving pathway production.     

Applications of genetically-encoded biosensors for the construction and control of biosynthetic pathways

Cells are filled with biosensors, molecular systems that measure the state of the cell and respond by regulating host processes. In much the same way that an engineer would monitor a chemical reactor, the cell uses these sensors to monitor changing intracellular environments and produce consistent behavior despite the variable environment. While natural systems derive a clear benefit from pathway regulation, past research efforts in engineering cellular metabolism have focused on introducing new pathways and removing existing pathway regulation. Synthetic biology is a rapidly growing field that focuses on the development of new tools that support the design, construction, and optimization of biological systems. Recent advances have been made in the design of genetically-encoded biosensors and the application of this class of molecular tools for optimizing and regulating heterologous pathways. Biosensors to cellular metabolites can be taken directly from natural systems, engineered from natural sensors, or constructed entirely in vitro. When linked to reporters, such as antibiotic resistance markers, these metabolite sensors can be used to report on pathway productivity, allowing high-throughput screening for pathway optimization. Future directions will focus on the application of biosensors to introduce feedback control into metabolic pathways, providing dynamic control strategies to increase the efficient use of cellular resources and pathway reliability.     

Identification of a novel Ser/Thr protein phosphatase Ppq1 as a negative regulator of mating MAP kinase pathway in Saccharomyces cerevisiae

The specificity and efficiency of cell signaling is largely governed by the complex formation of signaling proteins. The precise spatio-temporal control of the complex assembly is crucial for proper signaling and cell survival. Protein phosphorylation is a key mechanism of signal processing in most of cell signaling networks. Phosphatases, along with kinases, control the phosphorylation state of many proteins and thus play a critical role in the precise regulation of signaling at each stage such as activation, propagation, and adaptation. Identification and functional analysis of pathway-specific phosphatase is, therefore, crucial for the understanding of cell signaling mechanisms. Here, we have developed a novel screening strategy to identify pathway-specific phosphatases, in which the entire repertoire of cell’s phosphatases was tethered to a signaling complex and the changes in signaling response were monitored. As a model target, we have chosen the mating MAP kinase pathway in the budding yeast, which is composed of three kinases and Ste5 scaffold protein. Using this strategy, a putative Ser/Thr phosphatase, Ppq1, was identified to be mating-specific. Results show that Ppq1 down-regulates mating signaling by targeting at or upstream of the terminal MAP kinase Fus3 in the cascade. The catalytic activity of Ppq1 as a phosphatase was confirmed in vitro and is necessary for its function in the regulation of mating signaling. Overall, the data suggest that Ppq1 functions as a negative regulator of mating MAPK pathway by dephosphorylating target pathway protein(s) and plays a key role in the control of the background signaling noise.     

Kinetically guided,ratiometric tuning of fatty acid biosynthesis

Cellular metabolism is a nonlinear reaction network in which dynamic shifts in enzyme concentration help regulate the flux of carbon to different products. Despite the apparent simplicity of these biochemical adjustments, their influence on metabolite biosynthesis tends to be context-dependent, difficult to predict, and challenging to exploit in metabolic engineering. This study combines a detailed kinetic model with a systematic set of in vitro and in vivo analyses to explore the use of enzyme concentration as a control parameter in fatty acid synthesis, an essential metabolic process with important applications in oleochemical production. Compositional analyses of a modeled and experimentally reconstituted fatty acid synthase (FAS) from Escherichia coli indicate that the concentration ratio of two native enzymes—a promiscuous thioesterase and a ketoacyl synthase—can tune the average length of fatty acids, an important design objective of engineered pathways. The influence of this ratio is sensitive to the concentrations of other FAS components, which can narrow or expand the range of accessible chain lengths. Inside the cell, simple changes in enzyme concentration can enhance product-specific titers by as much as 125-fold and elicit shifts in overall product profiles that rival those of thioesterase mutants. This work develops a kinetically guided approach for using ratiometric adjustments in enzyme concentration to control the product profiles of FAS systems and, broadly, provides a detailed framework for understanding how coordinated shifts in enzyme concentration can afford tight control over the outputs of nonlinear metabolic pathways.     

Experimental determination of group flux control coefficients in metabolic networks

Grouping of reactions around key metabolite branch points can facilitate the study of metabolic control of complex metabolic networks. This top-down Metabolic Control Analysis is exemplified through the introduction of group (flux, as well as concentration) control coefficients whose magnitudes provide a measure of the relative impact of each reaction group on the overall network flux, as well as on the overall network stability, following enzymatic amplification. In this article, we demonstrate the application of previously developed theory to the determination of group flux control coefficients. Experimental data for the changes in metabolic fluxes obtained in response to the introduction of six different environmental perturbations are used to determine the group flux control coefficients for three reaction groups formed around the phosphoenolpyruvate/pyruvate branch point. The consistency of the obtained group flux control coefficient estimates is systematically analyzed to ensure that all necessary conditions are satisfied. The magnitudes of the determined control coefficients suggest that the control of lysine production flux in Corynebacterium glutamicum cells at a growth base state resides within the lysine biosynthetic pathway that begins with the PEP/PYR carboxylation anaplorotic pathway. Copyright 1998 John Wiley & Sons, Inc.     

A ubiquitin-dependent mitophagy complex maintains mitochondrial function and insulin secretion in beta cells

Mitochondrial autophagy or mitophagy is a key component of mitochondrial quality control, which is necessary to maintain cellular bioenergetics. Pancreatic islet β-cells, which release insulin in response to circulating blood glucose levels, are particularly susceptible to mitochondrial dysfunction due to their high metabolic activity and energy requirements for insulin processing, maturation, and secretion. Therefore, dysregulated mitophagy has drawn interest in the etiology of β-cell failure in diabetes. We demonstrate that the pivotal β-cell mitophagy regulator, CLEC16A, is an E3 ligase that forms a ubiquitin-dependent tripartite complex with RNF41/NRDP1 and USP8. Maintenance of the CLEC16A-RNF41-USP8 mitophagy complex is necessary for maximal cellular respiration and insulin secretion. Further, we observe that diabetogenic metabolic stressors, including elevated glucose and fatty acids, destabilize the CLEC16A-RNF41-USP8 complex and lead to β-cell apoptosis. Thus, the β-cell mitophagy pathway requires ubiquitin signals to stabilize the CLEC16A-RNF41-USP8 complex and maintain mitochondrial quality control.     

The RpoE Stress Response Pathway Mediates Reduction of the Virulence of Enteropathogenic Escherichia coli by Zinc

Zinc supplements are an effective clinical treatment for infantile diarrheal disease caused by enteric pathogens. Previous studies demonstrated that zinc acts on enteropathogenic Escherichia coli (EPEC) bacteria directly to suppress several virulence-related genes at a concentration that can be achieved by oral delivery of dietary zinc supplements. Our in vitro studies showed that a micromolar concentration of zinc induced the envelope stress response and suppressed virulence in EPEC, providing a possible mechanistic explanation for zinc''s therapeutic action. In this report, we investigated the molecular and physiological changes in EPEC induced by zinc. We found that micromolar concentrations of zinc reduced the bacterial growth rate without affecting viability. We observed increased membrane permeability caused by zinc. Zinc upregulated the RpoE-dependent envelope stress response pathway and suppressed EPEC virulence gene expression. RpoE alone was sufficient to inhibit virulence factor expression and to attenuate attaching and effacing lesion formation on human host cells. By mutational analysis we demonstrate that the DNA-binding motif of RpoE is necessary for suppression of the LEE1, but not the LEE4, operon. Predictably, inhibition of the RpoE-mediated envelope stress response in combination with micromolar concentrations of zinc reduced EPEC viability. In conclusion, zinc induces the RpoE and stress response pathways in EPEC, and the alternate sigma factor RpoE downregulates EPEC LEE and non-LEE virulence genes by multiple mechanisms.     

Multi-enzyme systems and recombinant cells for synthesis of valuable saccharides: Advances and perspectives

Saccharides have recently attracted considerable attention because of their biological functions and potential applications in the pharmaceutical, cosmetic and food industries. Over the decades, a large amount of enzymes involved in saccharide synthesis have been discovered and characterised with the aid of available genome sequences. The advancement of metabolic engineering and synthetic biology strategies facilitated the artificial pathway design and construction for production of multiple sugars in vitro and in vivo based on those characterized enzymes. This review presented a panoramic view of enzymes related to saccharide synthesis and gave the detailed information. Furthermore, we provide an extensive overview of the recent advances in the construction of cell-free reaction systems and engineering of microbial cells for the production of natural or unnatural saccharides. In addition, the future trends in the synthesis of sugars with high structural diversity through the combination of multiple pathways are presented and evaluated.     

Ksp1 kinase regulates autophagy via the target of rapamycin complex 1 (TORC1) pathway

Macroautophagy (hereafter autophagy) is a bulk degradation system conserved in all eukaryotes, which engulfs cytoplasmic components within double-membrane vesicles to allow their delivery to, and subsequent degradation within, the vacuole/lysosome. Autophagy activity is tightly regulated in response to the nutritional state of the cell and also to maintain organelle homeostasis. In nutrient-rich conditions, Tor kinase complex 1 (TORC1) is activated to inhibit autophagy, whereas inactivation of this complex in response to stress leads to autophagy induction; however, it is unclear how the activity of TORC1 is controlled to allow precise adjustments in autophagy activity. In this study, we performed genetic analyses in Saccharomyces cerevisiae to identify factors that regulate TORC1 activity. We determined that the Ksp1 kinase functions in part as a negative regulator of autophagy; deletion of KSP1 facilitated dephosphorylation of Atg13, a TORC1 substrate, which correlates with enhanced autophagy. These results suggest that Ksp1 down-regulates autophagy activity via the TORC1 pathway. The suppressive function of Ksp1 is partially activated by the Ras/cAMP-dependent protein kinase A (PKA), which is another negative regulator of autophagy. Our study therefore identifies Ksp1 as a new component that functions as part of the PKA and TORC1 signaling network to control the magnitude of autophagy.     

Improving lycopene production in Escherichia coli by engineering metabolic control

Metabolic engineering has achieved encouraging success in producing foreign metabolites in a variety of hosts. However, common strategies for engineering metabolic pathways focus on amplifying the desired enzymes and deregulating cellular controls. As a result, uncontrolled or deregulated metabolic pathways lead to metabolic imbalance and suboptimal productivity. Here we have demonstrated the second stage of metabolic engineering effort by designing and engineering a regulatory circuit to control gene expression in response to intracellular metabolic states. Specifically, we recruited and altered one of the global regulatory systems in Escherichia coli, the Ntr regulon, to control the engineered lycopene biosynthesis pathway. The artificially engineered regulon, stimulated by excess glycolytic flux through sensing of an intracellular metabolite, acetyl phosphate, controls the expression of two key enzymes in lycopene synthesis in response to flux dynamics. This intracellular control loop significantly enhanced lycopene production while reducing the negative impact caused by metabolic imbalance. Although we demonstrated this strategy for metabolite production, it can be extended into other fields where gene expression must be closely controlled by intracellular physiology, such as gene therapy.     

13C metabolic flux analysis of three divergent extremely thermophilic bacteria: Geobacillus sp. LC300, Thermus thermophilus HB8, and Rhodothermus marinus DSM 4252

Thermophilic organisms are being increasingly investigated and applied in metabolic engineering and biotechnology. The distinct metabolic and physiological characteristics of thermophiles, including broad substrate range and high uptake rates, coupled with recent advances in genetic tool development, present unique opportunities for strain engineering. However, poor understanding of the cellular physiology and metabolism of thermophiles has limited the application of systems biology and metabolic engineering tools to these organisms. To address this concern, we applied high resolution ¹³C metabolic flux analysis to quantify fluxes for three divergent extremely thermophilic bacteria from separate phyla: Geobacillus sp. LC300, Thermus thermophilus HB8, and Rhodothermus marinus DSM 4252. We performed 18 parallel labeling experiments, using all singly labeled glucose tracers for each strain, reconstructed and validated metabolic network models, measured biomass composition, and quantified precise metabolic fluxes for each organism. In the process, we resolved many uncertainties regarding gaps in pathway reconstructions and elucidated how these organisms maintain redox balance and generate energy. Overall, we found that the metabolisms of the three thermophiles were highly distinct, suggesting that adaptation to growth at high temperatures did not favor any particular set of metabolic pathways. All three strains relied heavily on glycolysis and TCA cycle to generate key cellular precursors and cofactors. None of the investigated organisms utilized the Entner-Doudoroff pathway and only one strain had an active oxidative pentose phosphate pathway. Taken together, the results from this study provide a solid foundation for future model building and engineering efforts with these and related thermophiles.