Exploring the gap between dynamic and constraint-based models of metabolism

Systems biology provides new approaches for metabolic engineering through the development of models and methods for simulation and optimization of microbial metabolism. Here we explore the relationship between two modeling frameworks in common use namely, dynamic models with kinetic rate laws and constraint-based flux models. We compare and analyze dynamic and constraint-based formulations of the same model of the central carbon metabolism of Escherichia coli. Our results show that, if unconstrained, the space of steady states described by both formulations is the same. However, the imposition of parameter-range constraints can be mapped into kinetically feasible regions of the solution space for the dynamic formulation that is not readily transferable to the constraint-based formulation. Therefore, with partial kinetic parameter knowledge, dynamic models can be used to generate constraints that reduce the solution space below that identified by constraint-based models, eliminating infeasible solutions and increasing the accuracy of simulation and optimization methods.     

Improving the accuracy of flux balance analysis through the implementation of carbon availability constraints for intracellular reactions

Constraint-based modeling methods, such as Flux Balance Analysis (FBA), have been extensively used to decipher complex, information rich -omics datasets to elicit system-wide behavioral patterns of cellular metabolism. FBA has been successfully used to gain insight in a wide range of applications, such as range of substrate utilization, product yields and to design metabolic engineering strategies to improve bioprocess performance. A well-known challenge associated with large genome-scale metabolic networks is that they result in underdetermined problem formulations. Consequently, rather than unique solutions, FBA and related methods examine ranges of reaction flux values that are consistent with the studied physiological conditions. The wider the reported flux ranges, the higher the uncertainty in the determination of basic reaction properties, limiting interpretability of and confidence in the results. Herein, we propose a new, computationally efficient approach that refines flux range predictions by constraining reaction fluxes on the basis of the elemental balance of carbon. We compared carbon constraint FBA (ccFBA) against experimentally-measured intracellular fluxes using the latest CHO GEM (iCHO1766) and were able to substantially improve the accuracy of predicted flux values compared with FBA. ccFBA can be used as a stand-alone method but is also compatible with and complimentary to other constraint-based approaches.     

生理代谢参数RQ在指导发酵过程优化中的应用

过程在线参数检测是进行发酵过程工程优化控制的基础。呼吸熵RQ是微生物胞内微观代谢流在宏观代谢参数上的响应,反映了微生物培养过程中底物的利用情况、产物和副产物的合成情况、及微观代谢途径通量的变化。结合发酵葡萄糖酸、2,3-丁二醇谷氨酸、柠檬酸、头孢菌素C、毕赤酵母α-干扰素产品的工业生产过程,分析了以微生物胞内微观代谢与宏观的生理参数RQ为指导的发酵工艺优化策略。RQ值在发酵过程中可以根据尾气数据进行在线采集,对指导通过宏观代谢参数的调控来最优化微生物胞内的代谢途经通量,提高目的产物的产率具有非常重要的意义。     

Flux balance analysis of the ammonia-oxidizing bacterium Nitrosomonas europaea ATCC19718 unravels specific metabolic activities while degrading toxic compounds

The ammonia-oxidizing bacterium Nitrosomonas europaea has been widely recognized as an important player in the nitrogen cycle as well as one of the most abundant members in microbial communities for the treatment of industrial or sewage wastewater. Its natural metabolic versatility and extraordinary ability to degrade environmental pollutants (e.g., aromatic hydrocarbons such as benzene and toluene) enable it to thrive under various harsh environmental conditions. Constraint-based metabolic models constructed from genome sequences enable quantitative insight into the central and specialized metabolism within a target organism. These genome-scale models have been utilized to understand, optimize, and design new strategies for improved bioprocesses. Reduced modeling approaches have been used to elucidate Nitrosomonas europaea metabolism at a pathway level. However, genome-scale knowledge about the simultaneous oxidation of ammonia and pollutant metabolism of N. europaea remains limited. Here, we describe the reconstruction, manual curation, and validation of the genome-scale metabolic model for N. europaea, iGC535. This reconstruction is the most accurate metabolic model for a nitrifying organism to date, reaching an average prediction accuracy of over 90% under several growth conditions. The manually curated model can predict phenotypes under chemolithotrophic and chemolithoorganotrophic conditions while oxidating methane and wastewater pollutants. Calculated flux distributions under different trophic conditions show that several key pathways are affected by the type of carbon source available, including central carbon metabolism and energy production.     

Flux analysis and metabolomics for systematic metabolic engineering of microorganisms

Rational engineering of metabolism is important for bio-production using microorganisms. Metabolic design based on in silico simulations and experimental validation of the metabolic state in the engineered strain helps in accomplishing systematic metabolic engineering. Flux balance analysis (FBA) is a method for the prediction of metabolic phenotype, and many applications have been developed using FBA to design metabolic networks. Elementary mode analysis (EMA) and ensemble modeling techniques are also useful tools for in silico strain design. The metabolome and flux distribution of the metabolic pathways enable us to evaluate the metabolic state and provide useful clues to improve target productivity. Here, we reviewed several computational applications for metabolic engineering by using genome-scale metabolic models of microorganisms. We also discussed the recent progress made in the field of metabolomics and ¹³C-metabolic flux analysis techniques, and reviewed these applications pertaining to bio-production development. Because these in silico or experimental approaches have their respective advantages and disadvantages, the combined usage of these methods is complementary and effective for metabolic engineering.     

Shrinking the metabolic solution space using experimental datasets

Constraint-based models of metabolism have been used in a variety of studies on drug discovery, metabolic engineering, evolution, and multi-species interactions. These genome-scale models can be generated for any sequenced organism since their main parameters (i.e., reaction stoichiometry) are highly conserved. Their relatively low parameter requirement makes these models easy to develop; however, these models often result in a solution space with multiple possible flux distributions, making it difficult to determine the precise flux state in the cell. Recent research efforts in this modeling field have investigated how additional experimental data, including gene expression, protein expression, metabolite concentrations, and kinetic parameters, can be used to reduce the solution space. This mini-review provides a summary of the data-driven computational approaches that are available for reducing the solution space and thereby improve predictions of intracellular fluxes by constraint-based models.     

Inferring Metabolic States in Uncharacterized Environments Using Gene-Expression Measurements

The large size of metabolic networks entails an overwhelming multiplicity in the possible steady-state flux distributions that are compatible with stoichiometric constraints. This space of possibilities is largest in the frequent situation where the nutrients available to the cells are unknown. These two factors: network size and lack of knowledge of nutrient availability, challenge the identification of the actual metabolic state of living cells among the myriad possibilities. Here we address this challenge by developing a method that integrates gene-expression measurements with genome-scale models of metabolism as a means of inferring metabolic states. Our method explores the space of alternative flux distributions that maximize the agreement between gene expression and metabolic fluxes, and thereby identifies reactions that are likely to be active in the culture from which the gene-expression measurements were taken. These active reactions are used to build environment-specific metabolic models and to predict actual metabolic states. We applied our method to model the metabolic states of Saccharomyces cerevisiae growing in rich media supplemented with either glucose or ethanol as the main energy source. The resulting models comprise about 50% of the reactions in the original model, and predict environment-specific essential genes with high sensitivity. By minimizing the sum of fluxes while forcing our predicted active reactions to carry flux, we predicted the metabolic states of these yeast cultures that are in large agreement with what is known about yeast physiology. Most notably, our method predicts the Crabtree effect in yeast cells growing in excess glucose, a long-known phenomenon that could not have been predicted by traditional constraint-based modeling approaches. Our method is of immediate practical relevance for medical and industrial applications, such as the identification of novel drug targets, and the development of biotechnological processes that use complex, largely uncharacterized media, such as biofuel production.     

Systematic engineering of TCA cycle for optimal production of a four-carbon platform chemical 4-hydroxybutyric acid in Escherichia coli

To address climate change and environmental problems, it is becoming increasingly important to establish biorefineries for the production of chemicals from renewable non-food biomass. Here we report the development of Escherichia coli strains capable of overproducing a four-carbon platform chemical 4-hybroxybutyric acid (4-HB). Because 4-HB production is significantly affected by aeration level, genome-scale metabolic model-based engineering strategies were designed under aerobic and microaerobic conditions with emphasis on oxidative/reductive TCA branches and glyoxylate shunt. Several different metabolic engineering strategies were employed to develop strains suitable for fermentation both under aerobic and microaerobic conditions. It was found that microaerobic condition was more efficient than aerobic condition in achieving higher titer and productivity of 4-HB. The final engineered strain produced 103.4 g/L of 4-HB by microaerobic fed-batch fermentation using glycerol. The aeration-dependent optimization strategy of TCA cycle will be useful for developing microbial strains producing other reduced derivative chemicals of TCA cycle intermediates.     

Leveraging knowledge engineering and machine learning for microbial bio-manufacturing

Genome scale modeling (GSM) predicts the performance of microbial workhorses and helps identify beneficial gene targets. GSM integrated with intracellular flux dynamics, omics, and thermodynamics have shown remarkable progress in both elucidating complex cellular phenomena and computational strain design (CSD). Nonetheless, these models still show high uncertainty due to a poor understanding of innate pathway regulations, metabolic burdens, and other factors (such as stress tolerance and metabolite channeling). Besides, the engineered hosts may have genetic mutations or non-genetic variations in bioreactor conditions and thus CSD rarely foresees fermentation rate and titer. Metabolic models play important role in design-build-test-learn cycles for strain improvement, and machine learning (ML) may provide a viable complementary approach for driving strain design and deciphering cellular processes. In order to develop quality ML models, knowledge engineering leverages and standardizes the wealth of information in literature (e.g., genomic/phenomic data, synthetic biology strategies, and bioprocess variables). Data driven frameworks can offer new constraints for mechanistic models to describe cellular regulations, to design pathways, to search gene targets, and to estimate fermentation titer/rate/yield under specified growth conditions (e.g., mixing, nutrients, and O₂). This review highlights the scope of information collections, database constructions, and machine learning techniques (such as deep learning and transfer learning), which may facilitate “Learn and Design” for strain development.     

基因组规模代谢网络模型的约束算法及其应用

工业微生物是微生物制造的核心,而微生物制造过程中微生物的生产性能是提高微生物制造过程效率的关键。结合基于约束条件的优化算法的基因组规模代谢网络模型是全局规模上认识、调控和优化工业微生物生理功能的重要平台。本文在简述基因组规模代谢网络模型构建的基础上,详细介绍了基于约束条件的优化算法的原理、分类和应用实例,为提高微生物制造过程效率奠定了坚实的基础。     

Computational Analysis of Reciprocal Association of Metabolism and Epigenetics in the Budding Yeast: A Genome-Scale Metabolic Model (GSMM) Approach

Metaboloepigenetics is a newly coined term in biological sciences that investigates the crosstalk between epigenetic modifications and metabolism. The reciprocal relation between biochemical transformations and gene expression regulation has been experimentally demonstrated in cancers and metabolic syndromes. In this study, we explored the metabolism-histone modifications crosstalk by topological analysis and constraint-based modeling approaches in the budding yeast. We constructed nine models through the integration of gene expression data of four mutated histone tails into a genome-scale metabolic model of yeast. Accordingly, we defined the centrality indices of the lowly expressed enzymes in the undirected enzyme-centric network of yeast by CytoHubba plug-in in Cytoscape. To determine the global effects of histone modifications on the yeast metabolism, the growth rate and the range of possible flux values of reactions, we used constraint-based modeling approach. Centrality analysis shows that the lowly expressed enzymes could affect and control the yeast metabolic network. Besides, constraint-based modeling results are in a good agreement with the experimental findings, confirming that the mutations in histone tails lead to non-lethal alterations in the yeast, but have diverse effects on the growth rate and reveal the functional redundancy.     

Genome-scale metabolic models as platforms for strain design and biological discovery

Genome-scale metabolic models (GEMs) have been developed and used in guiding systems’ metabolic engineering strategies for strain design and development. This strategy has been used in fermentative production of bio-based industrial chemicals and fuels from alternative carbon sources. However, computer-aided hypotheses building using established algorithms and software platforms for biological discovery can be integrated into the pipeline for strain design strategy to create superior strains of microorganisms for targeted biosynthetic goals. Here, I described an integrated workflow strategy using GEMs for strain design and biological discovery. Specific case studies of strain design and biological discovery using Escherichia coli genome-scale model are presented and discussed. The integrated workflow presented herein, when applied carefully would help guide future design strategies for high-performance microbial strains that have existing and forthcoming genome-scale metabolic models.     

Genome-scale analysis of Mannheimia succiniciproducens metabolism

Mannheimia succiniciproducens MBEL55E isolated from bovine rumen is a capnophilic gram-negative bacterium that efficiently produces succinic acid, an industrially important four carbon dicarboxylic acid. In order to design a metabolically engineered strain which is capable of producing succinic acid with high yield and productivity, it is essential to optimize the whole metabolism at the systems level. Consequently, in silico modeling and simulation of the genome-scale metabolic network was employed for genome-scale analysis and efficient design of metabolic engineering experiments. The genome-scale metabolic network of M. succiniciproducens consisting of 686 reactions and 519 metabolites was constructed based on reannotation and validation experiments. With the reconstructed model, the network structure and key metabolic characteristics allowing highly efficient production of succinic acid were deciphered; these include strong PEP carboxylation, branched TCA cycle, relative weak pyruvate formation, the lack of glyoxylate shunt, and non-PTS for glucose uptake. Constraints-based flux analyses were then carried out under various environmental and genetic conditions to validate the genome-scale metabolic model and to decipher the altered metabolic characteristics. Predictions based on constraints-based flux analysis were mostly in excellent agreement with the experimental data. In silico knockout studies allowed prediction of new metabolic engineering strategies for the enhanced production of succinic acid. This genome-scale in silico model can serve as a platform for the systematic prediction of physiological responses of M. succiniciproducens to various environmental and genetic perturbations and consequently for designing rational strategies for strain improvement.     

A guide to metabolic flux analysis in metabolic engineering: Methods,tools and applications

The field of metabolic engineering is primarily concerned with improving the biological production of value-added chemicals, fuels and pharmaceuticals through the design, construction and optimization of metabolic pathways, redirection of intracellular fluxes, and refinement of cellular properties relevant for industrial bioprocess implementation. Metabolic network models and metabolic fluxes are central concepts in metabolic engineering, as was emphasized in the first paper published in this journal, “Metabolic fluxes and metabolic engineering” (Metabolic Engineering, 1: 1–11, 1999). In the past two decades, a wide range of computational, analytical and experimental approaches have been developed to interrogate the capabilities of biological systems through analysis of metabolic network models using techniques such as flux balance analysis (FBA), and quantify metabolic fluxes using constrained-based modeling approaches such as metabolic flux analysis (MFA) and more advanced experimental techniques based on the use of stable-isotope tracers, i.e. ¹³C-metabolic flux analysis (¹³C-MFA). In this review, we describe the basic principles of metabolic flux analysis, discuss current best practices in flux quantification, highlight potential pitfalls and alternative approaches in the application of these tools, and give a broad overview of pragmatic applications of flux analysis in metabolic engineering practice.     

Dynamic flux balance analysis of batch fermentation: effect of genetic manipulations on ethanol production

In silico optimization of bioethanol production from lignocellulosic biomasses is investigated by combining process systems engineering approach and systems biology approach. Lignocellulosic biomass is an attractive sustainable carbon source for fermentative production of bioethanol. For enhanced ethanol production, metabolic engineering of wild-type strains—that can metabolize both hexose and pentose sugars or microbial consortia consisting of substrate-selective microbes—may be advantageous. This study presents a detailed in silico analysis of bioethanol production from glucose-xylose mixtures of various compositions by batch mono-culture and co-culture fermentation of specialized microbes. Dynamic flux balance models based on available genome-scale reconstructions of the microorganisms have been used to analyze bioethanol production, and the maximization of ethanol productivity is addressed by computing optimal aerobic–anaerobic switching times. Effects of ten metabolic engineering strategies that have been suggested in the literature for ethanol overproduction, have been evaluated for their efficiency in enhancing the ethanol productivity in the context of batch mono-culture and co-culture processes.     

Modeling <Emphasis Type="Italic">Lactococcus lactis</Emphasis> using a genome-scale flux model

Background  

Genome-scale flux models are useful tools to represent and analyze microbial metabolism. In this work we reconstructed the metabolic network of the lactic acid bacteria Lactococcus lactis and developed a genome-scale flux model able to simulate and analyze network capabilities and whole-cell function under aerobic and anaerobic continuous cultures. Flux balance analysis (FBA) and minimization of metabolic adjustment (MOMA) were used as modeling frameworks.     

Genome-scale model for Clostridium acetobutylicum: Part II. Development of specific proton flux states and numerically determined sub-systems

A regulated genome-scale model for Clostridium acetobutylicum ATCC 824 was developed based on its metabolic network reconstruction. To aid model convergence and limit the number of flux-vector possible solutions (the size of the phenotypic solution space), modeling strategies were developed to impose a new type of constraint at the endo-exo-metabolome interface. This constraint is termed the specific proton flux state, and its use enabled accurate prediction of the extracellular medium pH during vegetative growth of batch cultures. The specific proton flux refers to the influx or efflux of free protons (per unit biomass) across the cell membrane. A specific proton flux state encompasses a defined range of specific proton fluxes and includes all metabolic flux distributions resulting in a specific proton flux within this range. Effective simulation of time-course batch fermentation required the use of independent flux balance solutions from an optimum set of specific proton flux states. Using a real-coded genetic algorithm to optimize temporal bounds of specific proton flux states, we show that six separate specific proton flux states are required to model vegetative-growth metabolism and accurately predict the extracellular medium pH. Further, we define the apparent proton flux stoichiometry per weak acids efflux and show that this value decreases from approximately 3.5 mol of protons secreted per mole of weak acids at the start of the culture to approximately 0 at the end of vegetative growth. Calculations revealed that when specific weak acids production is maximized in vegetative growth, the net proton exchange between the cell and environment occurs primarily through weak acids efflux (apparent proton flux stoichiometry is 1). However, proton efflux through cation channels during the early stages of acidogenesis was found to be significant. We have also developed the concept of numerically determined sub-systems of genome-scale metabolic networks here as a sub-network with a one-dimensional null space basis set. A numerically determined sub-system was constructed in the genome-scale metabolic network to study the flux magnitudes and directions of acetylornithine transaminase, alanine racemase, and D-alanine transaminase. These results were then used to establish additional constraints for the genome-scale model.     

In silico analysis of bioethanol production from glucose/xylose mixtures during fed-batch fermentation of co-culture and mono-culture systems

This study presents a detailed in silico analysis of bioethanol production from glucose/xylose mixtures of various compositions by fed-batch co-culture and mono-culture fermentation of specialized microbes. The mono-culture consists of recombinant Saccharomyces cerevisise that can metabolize both hexose and pentose sugars while the co-culture system consists of substrate-selective microbes. Dynamic flux balance models based on available genome-scale reconstructions of the microorganisms have been used to analyze bioethanol production in fed-batch culture with constant feed rates and the maximization of ethanol productivity is addressed by computing optimal aerobic-anaerobic switching times. The simulation results clearly point to the superior performance of fed-batch fermentation of microbial co-culture against fed-batch fermentation of mono-culture for bioethanol production from glucose/xylose mixtures. A set of potential genetic engineering strategies for enhancement of S. cerevisiae and Escherichia coli strains performance have been identified. Such in silico predictions using genome-scale models provide valuable guidance for conducting in vivo metabolic engineering experiments.     

Metabolic flux analysis and metabolic engineering of microorganisms

Recent advances in metabolic flux analysis including genome-scale constraints-based flux analysis and its applications in metabolic engineering are reviewed. Various computational aspects of constraints-based flux analysis including genome-scale stoichiometric models, additional constraints used for the improved accuracy, and several algorithms for identifying the target genes to be manipulated are described. Also, some of the successful applications of metabolic flux analysis in metabolic engineering are reviewed. Finally, we discuss the limitations that need to be overcome to make the results of genome-scale flux analysis more realistically represent the real cell metabolism.