[Graph-based and constraint-based heterologous metabolic pathway design methods and application].

It is among the goals in metabolic engineering to construct microbial cell factories producing high-yield and high value-added target products, and an important solution is to design efficient synthetic pathway for the target products. However, due to the difference in metabolic capacity among microbial chassises, the available substrate and the yielded products are limited. Therefore, it is urgent to design related metabolic pathways to improve the production capacity. Existing metabolic engineering approaches to designing heterologous pathways are mainly based on biological experience, which are inefficient. Moreover, the yielded results are in no way comprehensive. However, systems biology provides new methods for heterologous pathway design, particularly the graph-based and constraint-based methods. Based on the databases containing rich metabolism information, they search for and uncover possible metabolic pathways with designated strategy (graph-based method) or algorithm (constraint-based method) and then screen out the optimal pathway to guide the modification of strains. In this paper, we reviewed the databases and algorithms for pathway design, and the applications in metabolic engineering and discussed the strengths and weaknesses of existing algorithms in practical application, hoping to provide a reference for the selection of optimal methods for the design of product synthesis pathway.

Graph-based and constraint-based heterologous metabolic pathway design methods and application / 生物工程学报

It is among the goals in metabolic engineering to construct microbial cell factories producing high-yield and high value-added target products, and an important solution is to design efficient synthetic pathway for the target products. However, due to the difference in metabolic capacity among microbial chassises, the available substrate and the yielded products are limited. Therefore, it is urgent to design related metabolic pathways to improve the production capacity. Existing metabolic engineering approaches to designing heterologous pathways are mainly based on biological experience, which are inefficient. Moreover, the yielded results are in no way comprehensive. However, systems biology provides new methods for heterologous pathway design, particularly the graph-based and constraint-based methods. Based on the databases containing rich metabolism information, they search for and uncover possible metabolic pathways with designated strategy (graph-based method) or algorithm (constraint-based method) and then screen out the optimal pathway to guide the modification of strains. In this paper, we reviewed the databases and algorithms for pathway design, and the applications in metabolic engineering and discussed the strengths and weaknesses of existing algorithms in practical application, hoping to provide a reference for the selection of optimal methods for the design of product synthesis pathway.

An energetic profile of Corynebacterium glutamicum underpinned by measured biomass yield on ATP.

The supply and usage of energetic cofactors in metabolism is a central concern for systems metabolic engineering, particularly in case of energy intensive products. One of the most important parameters for systems wide balancing of energetic cofactors is the ATP requirement for biomass formation YATP/Biomass. Despite its fundamental importance, YATP/Biomass values for non-fermentative organisms are still rough estimates deduced from theoretical considerations. For the first time, we present an approach for the experimental determination of YATP/Biomass using comparative 13C metabolic flux analysis (13C MFA) of a wild type strain and an ATP synthase knockout mutant. We show that the energetic profile of a cell can then be deduced from a genome wide stoichiometric model and experimental maintenance data. Particularly, the contributions of substrate level phosphorylation (SLP) and electron transport phosphorylation (ETP) to ATP generation become available which enables the overall energetic efficiency of a cell to be characterized. As a model organism, the industrial platform organism Corynebacterium glutamicum is used. C. glutamicum uses a respiratory type of energy metabolism, implying that ATP can be synthesized either by SLP or by ETP with the membrane-bound F1FO-ATP synthase using the proton motive force (pmf) as driving force. The presence of two terminal oxidases, which differ in their proton translocation efficiency by a factor of three, further complicates energy balancing for this organism. By integration of experimental data and network models, we show that in the wild type SLP and ETP contribute equally to ATP generation. Thus, the role of ETP in respiring bacteria may have been overrated in the past. Remarkably, in the genome wide setting 65% of the pmf is actually not used for ATP synthesis. However, it turns out that, compared to other organisms C. glutamicum still uses its energy budget rather efficiently.

Expanding Metabolic Capabilities Using Novel Pathway Designs: Computational Tools and Case Studies.

Design and selection of efficient metabolic pathways is critical for the success of metabolic engineering endeavors. Convenient pathways should not only produce the target metabolite in high yields but also are required to be thermodynamically feasible under production conditions, and to prefer efficient enzymes. To support the design and selection of such pathways, different computational approaches have been proposed for exploring the feasible pathway space under many of the above constraints. In this review, an overview of recent constraint-based optimization frameworks for metabolic pathway prediction, as well as relevant pathway engineering case studies that highlight the importance of rational metabolic designs is presented. Despite the availability and suitability of in silico design tools for metabolic pathway engineering, scarce-although increasing-application of computational outcomes is found. Finally, challenges and limitations hindering the broad adoption and successful application of these tools in metabolic engineering projects are discussed.

Approaches to Computational Strain Design in the Multiomics Era.

Modern omics analyses are able to effectively characterize the genetic, regulatory, and metabolic phenotypes of engineered microbes, yet designing genetic interventions to achieve a desired phenotype remains challenging. With recent developments in genetic engineering techniques, timelines associated with building and testing strain designs have been greatly reduced, allowing for the first time an efficient closed loop iteration between experiment and analysis. However, the scale and complexity associated with multi-omics datasets complicates manual biological reasoning about the mechanisms driving phenotypic changes. Computational techniques therefore form a critical part of the Design-Build-Test-Learn (DBTL) cycle in metabolic engineering. Traditional statistical approaches can reduce the dimensionality of these datasets and identify common motifs among high-performing strains. While successful in many studies, these methods do not take full advantage of known connections between genes, proteins, and metabolic networks. There is therefore a growing interest in model-aided design, in which modeling frameworks from systems biology are used to integrate experimental data and generate effective and non-intuitive design predictions. In this mini-review, we discuss recent progress and challenges in this field. In particular, we compare methods augmenting flux balance analysis with additional constraints from fluxomic, genomic, and metabolomic datasets and methods employing kinetic representations of individual metabolic reactions, and machine learning. We conclude with a discussion of potential future directions for improving strain design predictions in the omics era and remaining experimental and computational hurdles.

Integration of enzymatic data in Bacillus subtilis genome-scale metabolic model improves phenotype predictions and enables in silico design of poly-γ-glutamic acid production strains.

BACKGROUND: Genome-scale metabolic models (GEMs) allow predicting metabolic phenotypes from limited data on uptake and secretion fluxes by defining the space of all the feasible solutions and excluding physio-chemically and biologically unfeasible behaviors. The integration of additional biological information in genome-scale models, e.g., transcriptomic or proteomic profiles, has the potential to improve phenotype prediction accuracy. This is particularly important for metabolic engineering applications where more accurate model predictions can translate to more reliable model-based strain design. RESULTS: Here we present a GEM with Enzymatic Constraints using Kinetic and Omics data (GECKO) model of Bacillus subtilis, which uses publicly available proteomic data and enzyme kinetic parameters for central carbon (CC) metabolic reactions to constrain the flux solution space. This model allows more accurate prediction of the flux distribution and growth rate of wild-type and single-gene/operon deletion strains compared to a standard genome-scale metabolic model. The flux prediction error decreased by 43% and 36% for wild-type and mutants respectively. The model additionally increased the number of correctly predicted essential genes in CC pathways by 2.5-fold and significantly decreased flux variability in more than 80% of the reactions with variable flux. Finally, the model was used to find new gene deletion targets to optimize the flux toward the biosynthesis of poly-γ-glutamic acid (γ-PGA) polymer in engineered B. subtilis. We implemented the single-reaction deletion targets identified by the model experimentally and showed that the new strains have a twofold higher γ-PGA concentration and production rate compared to the ancestral strain. CONCLUSIONS: This work confirms that integration of enzyme constraints is a powerful tool to improve existing genome-scale models, and demonstrates the successful use of enzyme-constrained models in B. subtilis metabolic engineering. We expect that the new model can be used to guide future metabolic engineering efforts in the important industrial production host B. subtilis.

Optimization of bioprocess productivity based on metabolic-genetic network models with bilevel dynamic programming.

One of the main goals of metabolic engineering is to obtain high levels of a microbial product through genetic modifications. To improve the productivity of such a process, the dynamic implementation of metabolic engineering strategies has been proven to be more beneficial compared to static genetic manipulations in which the gene expression is not controlled over time, by resolving the trade-off between growth and production. In this work, a bilevel optimization framework based on constraint-based models is applied to identify an optimal strategy for dynamic genetic and process level manipulations to increase productivity. The dynamic enzyme-cost flux balance analysis (deFBA) as underlying metabolic network model captures the network dynamics and enables the analysis of temporal regulation in the metabolic-genetic network. We apply our computational framework to maximize ethanol productivity in a batch process with Escherichia coli. The results highlight the importance of integrating the genetic level and enzyme production and degradation processes for obtaining optimal dynamic gene and process manipulations.

Dynamic optimization of metabolic networks coupled with gene expression.

The regulation of metabolic activity by tuning enzyme expression levels is crucial to sustain cellular growth in changing environments. Metabolic networks are often studied at steady state using constraint-based models and optimization techniques. However, metabolic adaptations driven by changes in gene expression cannot be analyzed by steady state models, as these do not account for temporal changes in biomass composition. Here we present a dynamic optimization framework that integrates the metabolic network with the dynamics of biomass production and composition. An approximation by a timescale separation leads to a coupled model of quasi-steady state constraints on the metabolic reactions, and differential equations for the substrate concentrations and biomass composition. We propose a dynamic optimization approach to determine reaction fluxes for this model, explicitly taking into account enzyme production costs and enzymatic capacity. In contrast to the established dynamic flux balance analysis, our approach allows predicting dynamic changes in both the metabolic fluxes and the biomass composition during metabolic adaptations. Discretization of the optimization problems leads to a linear program that can be efficiently solved. We applied our algorithm in two case studies: a minimal nutrient uptake network, and an abstraction of core metabolic processes in bacteria. In the minimal model, we show that the optimized uptake rates reproduce the empirical Monod growth for bacterial cultures. For the network of core metabolic processes, the dynamic optimization algorithm predicted commonly observed metabolic adaptations, such as a diauxic switch with a preference ranking for different nutrients, re-utilization of waste products after depletion of the original substrate, and metabolic adaptation to an impending nutrient depletion. These examples illustrate how dynamic adaptations of enzyme expression can be predicted solely from an optimization principle.

Biomarker signature identification in "omics" data with multi-class outcome.

Biomarker signature identification in "omics" data is a complex challenge that requires specialized feature selection algorithms. The objective of these algorithms is to select the smallest set(s) of molecular quantities that are able to predict a given outcome (target) with maximal predictive performance. This task is even more challenging when the outcome comprises of multiple classes; for example, one may be interested in identifying the genes whose expressions allow discrimination among different types of cancer (nominal outcome) or among different stages of the same cancer, e.g. Stage 1, 2, 3 and 4 of Lung Adenocarcinoma (ordinal outcome). In this work, we consider a particular type of successful feature selection methods, named constraint-based, local causal discovery algorithms. These algorithms depend on performing a series of conditional independence tests. We extend these algorithms for the analysis of problems with continuous predictors and multi-class outcomes, by developing and equipping them with an appropriate conditional independence test procedure for both nominal and ordinal multi-class targets. The test is based on multinomial logistic regression and employs the log-likelihood ratio test for model selection. We present a comparative, experimental evaluation on seven real-world, high-dimensional, gene-expression datasets. Within the scope of our analysis the results indicate that the new conditional independence test allows the identification of smaller and better performing signatures for multi-class outcome datasets, with respect to the current alternatives for performing the independence tests.