Improved kinetic model of Escherichia coli central carbon metabolism in batch and continuous cultures

BayesianSSA: a Bayesian statistical model based on structural sensitivity analysis for predicting responses to enzyme perturbations in metabolic networks

BACKGROUND

Chemical bioproduction has attracted attention as a key technology in a decarbonized society. In computational design for chemical bioproduction, it is necessary to predict changes in metabolic fluxes when up-/down-regulating enzymatic reactions, that is, responses of the system to enzyme perturbations. Structural sensitivity analysis (SSA) was previously developed as a method to predict qualitative responses to enzyme perturbations on the basis of the structural information of the reaction network. However, the network structural information can sometimes be insufficient to predict qualitative responses unambiguously, which is a practical issue in bioproduction applications. To address this, in this study, we propose BayesianSSA, a Bayesian statistical model based on SSA. BayesianSSA extracts environmental information from perturbation datasets collected in environments of interest and integrates it into SSA predictions.

RESULTS

We applied BayesianSSA to synthetic and real datasets of the central metabolic pathway of Escherichia coli. Our result demonstrates that BayesianSSA can successfully integrate environmental information extracted from perturbation data into SSA predictions. In addition, the posterior distribution estimated by BayesianSSA can be associated with the known pathway reported to enhance succinate export flux in previous studies.

CONCLUSIONS

We believe that BayesianSSA will accelerate the chemical bioproduction process and contribute to advancements in the field.

Local flux coordination and global gene expression regulation in metabolic modeling

Genome-scale metabolic networks (GSMs) are fundamental systems biology representations of a cell's entire set of stoichiometrically balanced reactions. However, such static GSMs do not incorporate the functional organization of metabolic genes and their dynamic regulation (e.g., operons and regulons). Specifically, there are numerous topologically coupled local reactions through which fluxes are coordinated; the global growth state often dynamically regulates many gene expression of metabolic reactions via global transcription factor regulators. Here, we develop a GSM reconstruction method, Decrem, by integrating locally coupled reactions and global transcriptional regulation of metabolism by cell state. Decrem produces predictions of flux and growth rates, which are highly correlated with those experimentally measured in both wild-type and mutants of three model microorganisms Escherichia coli, Saccharomyces cerevisiae, and Bacillus subtilis under various conditions. More importantly, Decrem can also explain the observed growth rates by capturing the experimentally measured flux changes between wild-types and mutants. Overall, by identifying and incorporating locally organized and regulated functional modules into GSMs, Decrem achieves accurate predictions of phenotypes and has broad applications in bioengineering, synthetic biology, and microbial pathology.

A dynamic kinetic model captures cell-free metabolism for improved butanol production

Cell-free systems are useful tools for prototyping metabolic pathways and optimizing the production of various bioproducts. Mechanistically-based kinetic models are uniquely suited to analyze dynamic experimental data collected from cell-free systems and provide vital qualitative insight. However, to date, dynamic kinetic models have not been applied with rigorous biological constraints or trained on adequate experimental data to the degree that they would give high confidence in predictions and broadly demonstrate the potential for widespread use of such kinetic models. In this work, we construct a large-scale dynamic model of cell-free metabolism with the goal of understanding and optimizing butanol production in a cell-free system. Using a combination of parameterization methods, the resultant model captures experimental metabolite measurements across two experimental conditions for nine metabolites at timepoints between 0 and 24 h. We present analysis of the model predictions, provide recommendations for butanol optimization, and identify the aldehyde/alcohol dehydrogenase as the primary bottleneck in butanol production. Sensitivity analysis further reveals the extent to which various parameters are constrained, and our approach for probing valid parameter ranges can be applied to other modeling efforts.

Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework

Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, Linear Kinetics-Dynamic Flux Balance Analysis (LK-DFBA), that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this framework's success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that when genetic perturbations were implemented in the systems, the optimal constraint approach typically remained the same as for the wild-type regardless of the model topology or parameterization, indicating that just a single wild-type dataset could allow identification of the ideal constraint to enable model predictivity for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems.

Queueing theory model of Krebs cycle

MOTIVATION

Queueing theory can be effective in simulating biochemical reactions taking place in living cells, and the article paves a step toward development of a comprehensive model of cell metabolism. Such a model could help to accelerate and reduce costs for developing and testing investigational drugs reducing number of laboratory animals needed to evaluate drugs.

RESULTS

The article presents a Krebs cycle model based on queueing theory. The model allows for tracking of metabolites concentration changes in real time. To validate the model, a drug-induced inhibition affecting activity of enzymes involved in Krebs cycle was simulated and compared with available experimental data.

AVAILABILITYAND IMPLEMENTATION

The source code is freely available for download at https://github.com/UTP-WTIiE/KrebsCycleUsingQueueingTheory, implemented in C# supported in Linux or MS Windows.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Recent advances in metabolic engineering-integration of in silico design and experimental analysis of metabolic pathways

Microorganisms are widely used to produce valuable compounds. Because thousands of metabolic reactions occur simultaneously and many metabolic reactions are related to target production and cell growth, the development of a rational design method for metabolic pathway modification to optimize target production is needed. In this paper, recent advances in metabolic engineering are reviewed, specifically considering computational pathway modification design and experimental evaluation of metabolic fluxes by ¹³C-metabolic flux analysis. Computational tools for seeking effective gene deletion targets and recruiting heterologous genes are described in flux balance analysis approaches. A kinetic model and adaptive laboratory evolution were applied to identify and eliminate the rate-limiting step in metabolic pathways. Data science-based approaches for process monitoring and control are described to maximize the performance of engineered cells in bioreactors.

Building kinetic models for metabolic engineering

Kinetic formalisms of metabolism link metabolic fluxes to enzyme levels, metabolite concentrations and their allosteric regulatory interactions. Though they require the identification of physiologically relevant values for numerous parameters, kinetic formalisms uniquely establish a mechanistic link across heterogeneous omics datasets and provide an overarching vantage point to effectively inform metabolic engineering strategies. Advances in computational power, gene annotation coverage, and formalism standardization have led to significant progress over the past few years. However, careful interpretation of model predictions, limited metabolic flux datasets, and assessment of parameter sensitivity remain as challenges. In this review we highlight fundamental considerations which influence model quality and prediction, advances in methodologies, and success stories of deploying kinetic models to guide metabolic engineering.

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms

Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

Identification of a rate-limiting step in a metabolic pathway using the kinetic model and in vitro experiment

Identification of the rate-limiting step in a metabolic pathway is an important challenge in metabolic engineering for enhancing pathway flow. Although specific enzyme activities (V_max) provide valuable clues for the identification, it is time-consuming and difficult to measure multiple enzymes in the pathway because different assay protocols are required for each enzyme. In the present study, we propose a method to simultaneously determine the V_max values of multiple enzymes using a kinetic model with a time course of the intermediate concentrations through an in vitro experiment. To demonstrate this method, nine glycolysis reactions for converting glucose-6-phosphate (G6P) to pyruvate in Escherichia coli were considered. In a reaction mixture containing G6P and cofactors, glycolysis was initiated by adding a crude cell extract obtained from stationary phase cells. The V_max values were optimized to minimize the difference between the measured and simulated time-courses using a kinetic model. Metabolic control analysis using the kinetic model with the estimated V_max values revealed that fructose bisphosphate aldolase (FBA) was the rate-limiting step in the upper part of glycolysis. The addition of FBA in the reaction mixture successfully increased the glycolytic flux in vitro. Furthermore, in vivo, the specific glucose consumption rate of an FBA overexpression strain was 1.4 times higher than that of the control strain during the stationary phase. These results confirmed that FBA was the rate-limiting step in glycolysis under the stationary phase. This approach provides V_max values of multiple enzymes in a pathway for metabolic control analysis with a kinetic model.

Kinetic modeling of microbial growth, enzyme activity, and gene deletions: An integrated model of β-glucosidase function in Cellvibrio japonicus

Understanding the complex growth and metabolic dynamics in microorganisms requires advanced kinetic models containing both metabolic reactions and enzymatic regulation to predict phenotypic behaviors under different conditions and perturbations. Most current kinetic models lack gene expression dynamics and are separately calibrated to distinct media, which consequently makes them unable to account for genetic perturbations or multiple substrates. This challenge limits our ability to gain a comprehensive understanding of microbial processes towards advanced metabolic optimizations that are desired for many biotechnology applications. Here, we present an integrated computational and experimental approach for the development and optimization of mechanistic kinetic models for microbial growth and metabolic and enzymatic dynamics. Our approach integrates growth dynamics, gene expression, protein secretion, and gene-deletion phenotypes. We applied this methodology to build a dynamic model of the growth kinetics in batch culture of the bacterium Cellvibrio japonicus grown using either cellobiose or glucose media. The model parameters were inferred from an experimental data set using an evolutionary computation method. The resulting model was able to explain the growth dynamics of C. japonicus using either cellobiose or glucose media and was also able to accurately predict the metabolite concentrations in the wild-type strain as well as in β-glucosidase gene deletion mutant strains. We validated the model by correctly predicting the non-diauxic growth and metabolite consumptions of the wild-type strain in a mixed medium containing both cellobiose and glucose, made further predictions of mutant strains growth phenotypes when using cellobiose and glucose media, and demonstrated the utility of the model for designing industrially-useful strains. Importantly, the model is able to explain the role of the different β-glucosidases and their behavior under genetic perturbations. This integrated approach can be extended to other metabolic pathways to produce mechanistic models for the comprehensive understanding of enzymatic functions in multiple substrates.

Prioritization of Mur family drug targets against A. baumannii and identification of their homologous proteins through molecular phylogeny, primary sequence, and structural analysis

Background

The World Health Organization (WHO) report stated that Acinetobacter baumannii had been classified as one of the most important pathogenic bacteria causing nosocomial infection in hospital patients due to multi-drug resistance (MDR). It is vital to find out new bacterial drug targets and annotated their structure and function for the exploration of new anti-bacterial agents. The present study utilized a systematic route to prioritize the potential drug targets that belong to Mur family of Acinetobacter baumannii and identify their homologous proteins using a computational approach such as sequence similarity search, multiple sequence alignment, phylogenetic analysis, protein sequence, and protein structure analysis.

Results

From the results of protein sequence analysis of eight Mur family proteins, they divided into three main enzymatic classes namely transferases (MurG, MurA and MraY), ligases (MurC, MurD, MurE, and MurF), and oxidoreductase (MurB). Based on the results of intra-comparative protein sequence analysis and enzymatic classification, we have chosen MurB, MurE, and MurG as the prioritized drug targets from A. baumannii and subjected them for further detailed studies of inter-species comparison. This inter-species comparison help us to explore the sequential and structural properties of homologous proteins in other species and hence, opens a gateway for new target identification and using common inhibitor for different bacterial species caused by various diseases. The pairwise sequence alignment results between A. baumannii’s MurB with A. calcoaceticus’s MurB, A. baumannii’s MurE with A. seifertii’s MurE, and A. baumannii’s MurG with A. pittii’s MurG showed that every group of the proteins are highly similar with each other and they showed sequence identity of 95.7% and sequence similarity of 97.2%.

Conclusion

Together with the results of secondary and three-dimensional structure predictions explained that three selected proteins (MurB, MurE, and MurG) from A. baumannii and their related proteins (AcMurB, AsMurE, and ApMurG) belong to mixed αβ class and they are very similar.

Simultaneous cross-evaluation of heterogeneous E. coli datasets via mechanistic simulation

The extensive heterogeneity of biological data poses challenges to analysis and interpretation. Construction of a large-scale mechanistic model of Escherichia coli enabled us to integrate and cross-evaluate a massive, heterogeneous dataset based on measurements reported by various groups over decades. We identified inconsistencies with functional consequences across the data, including that the total output of the ribosomes and RNA polymerases described by data are not sufficient for a cell to reproduce measured doubling times, that measured metabolic parameters are neither fully compatible with each other nor with overall growth, and that essential proteins are absent during the cell cycle-and the cell is robust to this absence. Finally, considering these data as a whole leads to successful predictions of new experimental outcomes, in this case protein half-lives.

Computer-Aided Rational Design of Efficient NADPH Production System by Escherichia coli pgi Mutant Using a Mixture of Glucose and Xylose

Lignocellulosic biomass can be hydrolyzed into two major sugars of glucose and xylose, and thus the strategy for the efficient consumption of both sugars is highly desirable. NADPH is the essential molecule for the production of industrially important value-added chemicals, and thus its availability is quite important. Escherichia coli mutant lacking the pgi gene encoding phosphoglucose isomerase (Pgi) has been preferentially used to overproduce the NADPH. However, there exists a disadvantage that the cell growth rate becomes low for the mutant grown on glucose. This limits the efficient NADPH production, and therefore, it is quite important to investigate how addition of different carbon source such as xylose (other than glucose) effectively improves the NADPH production. In this study, we have developed a kinetic model to propose an efficient NADPH production system using E. coli pgi-knockout mutant with a mixture of glucose and xylose. The proposed system adds xylose to glucose medium to recover the suppressed growth of the pgi mutant, and determines the xylose content to maximize the NADPH productivity. Finally, we have designed a mevalonate (MVA) production system by implementing ArcA overexpression into the pgi-knockout mutant using a mixture of glucose and xylose. In addition to NADPH overproduction, the accumulation of acetyl-CoA (AcCoA) is necessary for the efficient MVA production. In the present study, therefore, we considered to overexpress ArcA, where ArcA overexpression suppresses the TCA cycle, causing the overflow of AcCoA, a precursor of MVA. We predicted the xylose content that maximizes the MVA production. This approach demonstrates the possibility of a great progress in the computer-aided rational design of the microbial cell factories for useful metabolite production.

Autocatalytic networks in biology: structural theory and algorithms

Self-sustaining autocatalytic networks play a central role in living systems, from metabolism at the origin of life, simple RNA networks and the modern cell, to ecology and cognition. A collectively autocatalytic network that can be sustained from an ambient food set is also referred to more formally as a 'reflexively autocatalytic food-generated' (RAF) set. In this paper, we first investigate a simplified setting for studying RAFs, which is nevertheless relevant to real biochemistry and which allows an exact mathematical analysis based on graph-theoretic concepts. This, in turn, allows for the development of efficient (polynomial-time) algorithms for questions that are computationally intractable (NP-hard) in the general RAF setting. We then show how this simplified setting for RAF systems leads naturally to a more general notion of RAFs that are 'generative' (they can be built up from simpler RAFs) and for which efficient algorithms carry over to this more general setting. Finally, we show how classical RAF theory can be extended to deal with ensembles of catalysts as well as the assignment of rates to reactions according to which catalysts (or combinations of catalysts) are available.

The Rocky Road From Fed-Batch to Continuous Processing With E. coli

Escherichia coli still serves as a beloved workhorse for the production of many biopharmaceuticals as it fulfills essential criteria, such as having fast doubling times, exhibiting a low risk of contamination, and being easy to upscale. Most industrial processes in E. coli are carried out in fed-batch mode. However, recent trends show that the biotech industry is moving toward time-independent processing, trying to improve the space-time yield, and especially targeting constant quality attributes. In the 1950s, the term "chemostat" was introduced for the first time by Novick and Szilard, who followed up on the previous work performed by Monod. Chemostat processing resulted in a major hype 10 years after its official introduction. However, enthusiasm decreased as experiments suffered from genetic instabilities and physiology issues. Major improvements in strain engineering and the usage of tunable promotor systems facilitated chemostat processes. In addition, critical process parameters have been identified, and the effects they have on diverse quality attributes are understood in much more depth, thereby easing process control. By pooling the knowledge gained throughout the recent years, new applications, such as parallelization, cascade processing, and population controls, are applied nowadays. However, to control the highly heterogeneous cultivation broth to achieve stable productivity throughout long-term cultivations is still tricky. Within this review, we discuss the current state of E. coli fed-batch process understanding and its tech transfer potential within continuous processing. Furthermore, the achievements in the continuous upstream applications of E. coli and the continuous downstream processing of intracellular proteins will be discussed.

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.

From Escherichia coli mutant 13C labeling data to a core kinetic model: A kinetic model parameterization pipeline

Kinetic models of metabolic networks offer the promise of quantitative phenotype prediction. The mechanistic characterization of enzyme catalyzed reactions allows for tracing the effect of perturbations in metabolite concentrations and reaction fluxes in response to genetic and environmental perturbation that are beyond the scope of stoichiometric models. In this study, we develop a two-step computational pipeline for the rapid parameterization of kinetic models of metabolic networks using a curated metabolic model and available 13C-labeling distributions under multiple genetic and environmental perturbations. The first step involves the elucidation of all intracellular fluxes in a core model of E. coli containing 74 reactions and 61 metabolites using 13C-Metabolic Flux Analysis (13C-MFA). Here, fluxes corresponding to the mid-exponential growth phase are elucidated for seven single gene deletion mutants from upper glycolysis, pentose phosphate pathway and the Entner-Doudoroff pathway. The computed flux ranges are then used to parameterize the same (i.e., k-ecoli74) core kinetic model for E. coli with 55 substrate-level regulations using the newly developed K-FIT parameterization algorithm. The K-FIT algorithm employs a combination of equation decomposition and iterative solution techniques to evaluate steady-state fluxes in response to genetic perturbations. k-ecoli74 predicted 86% of flux values for strains used during fitting within a single standard deviation of 13C-MFA estimated values. By performing both tasks using the same network, errors associated with lack of congruity between the two networks are avoided, allowing for seamless integration of data with model building. Product yield predictions and comparison with previously developed kinetic models indicate shifts in flux ranges and the presence or absence of mutant strains delivering flux towards pathways of interest from training data significantly impact predictive capabilities. Using this workflow, the impact of completeness of fluxomic datasets and the importance of specific genetic perturbations on uncertainties in kinetic parameter estimation are evaluated.

Dynamic metabolic control: towards precision engineering of metabolism

Advances in metabolic engineering have led to the synthesis of a wide variety of valuable chemicals in microorganisms. The key to commercializing these processes is the improvement of titer, productivity, yield, and robustness. Traditional approaches to enhancing production use the “push–pull-block” strategy that modulates enzyme expression under static control. However, strains are often optimized for specific laboratory set-up and are sensitive to environmental fluctuations. Exposure to sub-optimal growth conditions during large-scale fermentation often reduces their production capacity. Moreover, static control of engineered pathways may imbalance cofactors or cause the accumulation of toxic intermediates, which imposes burden on the host and results in decreased production. To overcome these problems, the last decade has witnessed the emergence of a new technology that uses synthetic regulation to control heterologous pathways dynamically, in ways akin to regulatory networks found in nature. Here, we review natural metabolic control strategies and recent developments in how they inspire the engineering of dynamically regulated pathways. We further discuss the challenges of designing and engineering dynamic control and highlight how model-based design can provide a powerful formalism to engineer dynamic control circuits, which together with the tools of synthetic biology, can work to enhance microbial production.