Itaconic acid production by co-feeding of Ustilago maydis: A combined approach of experimental data, design of experiments, and metabolic modeling

Itaconic acid is a platform chemical with a range of applications in polymer synthesis and is also discussed for biofuel production. While produced in industry from glucose or sucrose, co-feeding of glucose and acetate was recently discussed to increase itaconic acid production by the smut fungus Ustilago maydis. In this study, we investigate the optimal co-feeding conditions by interlocking experimental and computational methods. Flux balance analysis indicates that acetate improves the itaconic acid yield up to a share of 40% acetate on a carbon molar basis. A design of experiment results in the maximum yield of 0.14 itaconic acid per carbon source from 100 g L - 1 $\,\text{g L}{}^{-1}$ glucose and 12 g L - 1 $\,\text{g L}{}^{-1}$ acetate. The yield is improved by around 22% when compared to feeding of glucose as sole carbon source. To further improve the yield, gene deletion targets are discussed that were identified using the metabolic optimization tool OptKnock. The study contributes ideas to reduce land use for biotechnology by incorporating acetate as co-substrate, a C2-carbon source that is potentially derived from carbon dioxide.

Self-controlled in silico gene knockdown strategies to enhance the sustainable production of heterologous terpenoid by Saccharomyces cerevisiae

Microbial bioengineering is a growing field for producing plant natural products (PNPs) in recent decades, using heterologous metabolic pathways in host cells. Once heterologous metabolic pathways have been introduced into host cells, traditional metabolic engineering techniques are employed to enhance the productivity and yield of PNP biosynthetic routes, as well as to manage competing pathways. The advent of computational biology has marked the beginning of a novel epoch in strain design through in silico methods. These methods utilize genome-scale metabolic models (GEMs) and flux optimization algorithms to facilitate rational design across the entire cellular metabolic network. However, the implementation of in silico strategies can often result in an uneven distribution of metabolic fluxes due to the rigid knocking out of endogenous genes, which can impede cell growth and ultimately impact the accumulation of target products. In this study, we creatively utilized synthetic biology to refine in silico strain design for efficient PNPs production. OptKnock simulation was performed on the GEM of Saccharomyces cerevisiae OA07, an engineered strain for oleanolic acid (OA) bioproduction that has been reported previously. The simulation predicted that the single deletion of fol1, fol2, fol3, abz1, and abz2, or a combined knockout of hfd1, ald2 and ald3 could improve its OA production. Consequently, strains EK1∼EK7 were constructed and cultivated. EK3 (OA07△fol3), EK5 (OA07△abz1), and EK6 (OA07△abz2) had significantly higher OA titers in a batch cultivation compared to the original strain OA07. However, these increases were less pronounced in the fed-batch mode, indicating that gene deletion did not support sustainable OA production. To address this, we designed a negative feedback circuit regulated by malonyl-CoA, a growth-associated intermediate whose synthesis served as a bypass to OA synthesis, at fol3, abz1, abz2, and at acetyl-CoA carboxylase-encoding gene acc1, to dynamically and autonomously regulate the expression of these genes in OA07. The constructed strains R_3A, R_5A and R_6A had significantly higher OA titers than the initial strain and the responding gene-knockout mutants in either batch or fed-batch culture modes. Among them, strain R_3A stand out with the highest OA titer reported to date. Its OA titer doubled that of the initial strain in the flask-level fed-batch cultivation, and achieved at 1.23 ± 0.04 g L-1 in 96 h in the fermenter-level fed-batch mode. This indicated that the integration of optimization algorithm and synthetic biology approaches was efficiently rational for PNP-producing strain design.

FastKnock: an efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed by making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies (with a predefined maximum number of reaction deletions) for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using various Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple- and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more beneficial and important practical solutions. The availability of all the solutions provides the opportunity to further characterize, rank and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock .

Model validation and selection in metabolic flux analysis and flux balance analysis

13C-Metabolic Flux Analysis (13C-MFA) and Flux Balance Analysis (FBA) are widely used to investigate the operation of biochemical networks in both biological and biotechnological research. Both methods use metabolic reaction network models of metabolism operating at steady state so that reaction rates (fluxes) and the levels of metabolic intermediates are constrained to be invariant. They provide estimated (MFA) or predicted (FBA) values of the fluxes through the network in vivo, which cannot be measured directly. These fluxes can shed light on basic biology and have been successfully used to inform metabolic engineering strategies. Several approaches have been taken to test the reliability of estimates and predictions from constraint-based methods and to compare alternative model architectures. Despite advances in other areas of the statistical evaluation of metabolic models, such as the quantification of flux estimate uncertainty, validation and model selection methods have been underappreciated and underexplored. We review the history and state-of-the-art in constraint-based metabolic model validation and model selection. Applications and limitations of the χ2 -test of goodness-of-fit, the most widely used quantitative validation and selection approach in 13C-MFA, are discussed, and complementary and alternative forms of validation and selection are proposed. A combined model validation and selection framework for 13C-MFA incorporating metabolite pool size information that leverages new developments in the field is presented and advocated for. Finally, we discuss how adopting robust validation and selection procedures can enhance confidence in constraint-based modeling as a whole and ultimately facilitate more widespread use of FBA in biotechnology.

A Fermentation State Marker Rule Design Task in Metabolic Engineering

There are several ways in which mathematical modeling is used in fermentation control, but mechanistic mathematical genome-scale models of metabolism within the cell have not been applied or implemented so far. As part of the metabolic engineering task setting, we propose that metabolite fluxes and/or biomass growth rate be used to search for a fermentation steady state marker rule. During fermentation, the bioreactor control system can automatically detect the desired steady state using a logical marker rule. The marker rule identification can be also integrated with the production growth coupling approach, as presented in this study. A design of strain with marker rule is demonstrated on genome scale metabolic model iML1515 of Escherichia coli MG1655 proposing two gene deletions enabling a measurable marker rule for succinate production using glucose as a substrate. The marker rule example at glucose consumption 10.0 is: IF (specific growth rate μ is above 0.060 h-1, AND CO2 production under 1.0, AND ethanol production above 5.5), THEN succinate production is within the range 8.2-10, where all metabolic fluxes units are mmol ∗ gDW-1 ∗ h-1. An objective function for application in metabolic engineering, including productivity features and rule detecting sensor set characterizing parameters, is proposed. Two-phase approach to implementing marker rules in the cultivation control system is presented to avoid the need for a modeler during production.

OptEnvelope: A target point guided method for growth-coupled production using knockouts

Finding the best knockout strategy for coupling biomass growth and production of a target metabolite using a mathematic model of metabolism is a challenge in biotechnology. In this research, a three-step method named OptEnvelope is presented based on finding minimal set of active reactions for a target point in the feasible solution space (envelope) using a mixed-integer linear programming formula. The method initially finds the reduced desirable solution space envelope in the product versus biomass plot by removing all inactive reactions. Then, with reinsertion of the deleted reactions, OptEnvelope attempts to reduce the number of knockouts so that the desirable production envelope is preserved. Additionally, OptEnvelope searches for envelopes with higher minimum production rates or fewer knockouts by evaluating different target points within the desired solution space. It is possible to limit the maximal number of knockouts. The method was implemented on metabolic models of E. coli and S. cerevisiae to test the method benchmarking the capability of these industrial microbes for overproduction of acetate and glycerol under aerobic conditions and succinate and ethanol under anaerobic conditions. The results illustrate that OptEnvelope is capable to find multiple strong coupled envelopes located in the desired solution space because of its novel target point oriented strategy of envelope search. The results indicate that E. coli is more appropriate to produce acetate and succinate while S. cerevisiae is a better host for glycerol production. Gene deletions for some of the proposed reaction knockouts have been previously reported to increase the production of these metabolites in experiments. Both organisms are suitable for ethanol production, however, more knockouts for the adaptation of E. coli are required. OptEnvelope is available at https://github.com/lv-csbg/optEnvelope.

Uses of Multi-Objective Flux Analysis for Optimization of Microbial Production of Secondary Metabolites

Secondary metabolites are not essential for the growth of microorganisms, but they play a critical role in how microbes interact with their surroundings. In addition to this important ecological role, secondary metabolites also have a variety of agricultural, medicinal, and industrial uses, and thus the examination of secondary metabolism of plants and microbes is a growing scientific field. While the chemical production of certain secondary metabolites is possible, industrial-scale microbial production is a green and economically attractive alternative. This is even more true, given the advances in bioengineering that allow us to alter the workings of microbes in order to increase their production of compounds of interest. This type of engineering requires detailed knowledge of the "chassis" organism's metabolism. Since the resources and the catalytic capacity of enzymes in microbes is finite, it is important to examine the tradeoffs between various bioprocesses in an engineered system and alter its working in a manner that minimally perturbs the robustness of the system while allowing for the maximum production of a product of interest. The in silico multi-objective analysis of metabolism using genome-scale models is an ideal method for such examinations.

FastKnock: An efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using three Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more useful and important practical solutions. The availability of all the solutions provides the opportunity to further characterize and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock.

Co-feeding enhances the yield of methyl ketones

The biotechnological production of methyl ketones is a sustainable alternative to fossil-derived chemical production. To date, the best host for microbial production of methyl ketones is a genetically engineered Pseudomonas taiwanensis VLB120 ∆6 pProd strain, achieving yields of 101 mgg-1 on glucose in batch cultivations. For competitiveness with the petrochemical production pathway, however, higher yields are necessary. Co-feeding can improve the yield by fitting the carbon-to-energy ratio to the organism and the target product. In this work, we developed co-feeding strategies for P. taiwanensis VLB120 ∆6 pProd by combined metabolic modeling and experimental work. In a first step, we conducted flux balance analysis with an expanded genome-scale metabolic model of iJN1463 and found ethanol as the most promising among five cosubstrates. Next, we performed cultivations with ethanol and found the highest reported yield in batch production of methyl ketones with P. taiwanensis VLB120 to date, namely, 154 mg g-1 methyl ketones. However, ethanol is toxic to the cell, which reflects in a lower substrate consumption and lower product concentrations when compared to production on glucose. Hence, we propose cofeeding ethanol with glucose and find that, indeed, higher concentrations than in ethanol-fed cultivation (0.84 g Laq-1 with glucose and ethanol as opposed to 0.48 g Laq-1 with only ethanol) were achieved, with a yield of 85 mg g-1. In a last step, comparing experimental with computational results suggested the potential for improving the methyl ketone yield by fed-batch cultivation, in which cell growth and methyl ketone production are separated into two phases employing optimal ethanol to glucose ratios.

ONE-SENTENCE SUMMARY

By combining computational and experimental work, we demonstrate that feeding ethanol in addition to glucose improves the yield of biotechnologically produced methyl ketones.

Gene Deletion Algorithms for Minimum Reaction Network Design by Mixed-Integer Linear Programming for Metabolite Production in Constraint-Based Models: gDel_minRN

Genome-scale constraint-based metabolic networks play an important role in the simulation of growth-coupled production, which means that cell growth and target metabolite production are simultaneously achieved. For growth-coupled production, a minimal reaction-network-based design is known to be effective. However, the obtained reaction networks often fail to be realized by gene deletions due to conflicts with gene-protein-reaction (GPR) relations. Here, we developed gDel_minRN that determines gene deletion strategies using mixed-integer linear programming to achieve growth-coupled production by repressing the maximum number of reactions via GPR relations. The results of computational experiments showed that gDel_minRN could determine the core parts, which include only 30% to 55% of whole genes, for stoichiometrically feasible growth-coupled production for many target metabolites, which include useful vitamins such as biotin (vitamin B7), riboflavin (vitamin B2), and pantothenate (vitamin B5). Since gDel_minRN calculates a constraint-based model of the minimum number of gene-associated reactions without conflict with GPR relations, it helps biological analysis of the core parts essential for growth-coupled production for each target metabolite. The source codes, implemented in MATLAB using CPLEX and COBRA Toolbox, are available on https://github.com/MetNetComp/gDel-minRN.

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

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

StrainDesign: a comprehensive Python package for computational design of metabolic networks

SUMMARY

Various constraint-based optimization approaches have been developed for the computational analysis and design of metabolic networks. Herein, we present StrainDesign, a comprehensive Python package that builds upon the COBRApy toolbox and integrates the most popular metabolic design algorithms, including nested strain optimization methods such as OptKnock, RobustKnock and OptCouple as well as the more general minimal cut sets approach. The optimization approaches are embedded in individual modules, which can also be combined for setting up more elaborate strain design problems. Advanced features, such as the efficient integration of GPR rules and the possibility to consider gene and reaction additions or regulatory interventions, have been generalized and are available for all modules. The package uses state-of-the-art preprocessing methods, supports multiple solvers and provides a number of enhanced tools for analyzing computed intervention strategies including 2D and 3D plots of user-selected metabolic fluxes or yields. Furthermore, a user-friendly interface for the StrainDesign package has been implemented in the GUI-based metabolic modeling software CNApy. StrainDesign provides thus a unique and rich framework for computational strain design in Python, uniting many algorithmic developments in the field and allowing modular extension in the future.

AVAILABILITY AND IMPLEMENTATION

The StrainDesign package can be retrieved from PyPi, Anaconda and GitHub (https://github.com/klamt-lab/straindesign) and is also part of the latest CNApy package.

Genome-scale metabolic network models: from first-generation to next-generation

Over the last two decades, thousands of genome-scale metabolic network models (GSMMs) have been constructed. These GSMMs have been widely applied in various fields, ranging from network interaction analysis, to cell phenotype prediction. However, due to the lack of constraints, the prediction accuracy of first-generation GSMMs was limited. To overcome these limitations, the next-generation GSMMs were developed by integrating omics data, adding constrain condition, integrating different biological models, and constructing whole-cell models. Here, we review recent advances of GSMMs from the first generation to the next generation. Then, we discuss the major application of GSMMs in industrial biotechnology, such as predicting phenotypes and guiding metabolic engineering. In addition, human health applications, including understanding biological mechanisms, discovering biomarkers and drug targets, are also summarized. Finally, we address the challenges and propose new trend of GSMMs. KEY POINTS: •This mini-review updates the literature on almost all published GSMMs since 1999. •Detailed insights into the development of the first- and next-generation GSMMs. •The application of GSMMs is summarized, and the prospects of integrating machine learning are emphasized.

gcFront: a tool for determining a Pareto front of growth-coupled cell factory designs

Motivation

A widely applicable strategy to create cell factories is to knockout (KO) genes or reactions to redirect cell metabolism so that chemical synthesis is made obligatory when the cell grows at its maximum rate. Synthesis is thus growth-coupled, and the stronger the coupling the more deleterious any impediments in synthesis are to cell growth, making high producer phenotypes evolutionarily robust. Additionally, we desire that these strains grow and synthesize at high rates. Genome-scale metabolic models can be used to explore and identify KOs that growth-couple synthesis, but these are rare in an immense design space, making the search difficult and slow.

Results

To address this multi-objective optimization problem, we developed a software tool named gcFront—using a genetic algorithm it explores KOs that maximize cell growth, product synthesis and coupling strength. Moreover, our measure of coupling strength facilitates the search so that gcFront not only finds a growth-coupled design in minutes but also outputs many alternative Pareto optimal designs from a single run—granting users flexibility in selecting designs to take to the lab.

Availability and implementation

gcFront, with documentation and a workable tutorial, is freely available at GitHub: https://github.com/lLegon/gcFront and archived at Zenodo, DOI: 10.5281/zenodo.5557755.

Supplementary information

Supplementary data are available at Bioinformatics online.

In Silico Design Strategies for the Production of Target Chemical Compounds Using Iterative Single-Level Linear Programming Problems

The optimization of metabolic reaction modifications for the production of target compounds is a complex computational problem whose execution time increases exponentially with the number of metabolic reactions. Therefore, practical technologies are needed to identify reaction deletion combinations to minimize computing times and promote the production of target compounds by modifying intracellular metabolism. In this paper, a practical metabolic design technology named AERITH is proposed for high-throughput target compound production. This method can optimize the production of compounds of interest while maximizing cell growth. With this approach, an appropriate combination of metabolic reaction deletions can be identified by solving a simple linear programming problem. Using a standard CPU, the computation time could be as low as 1 min per compound, and the system can even handle large metabolic models. AERITH was implemented in MATLAB and is freely available for non-profit use.

OptDesign: Identifying Optimum Design Strategies in Strain Engineering for Biochemical Production

Computational tools have been widely adopted for strain optimization in metabolic engineering, contributing to numerous success stories of producing industrially relevant biochemicals. However, most of these tools focus on single metabolic intervention strategies (either gene/reaction knockout or amplification alone) and rely on hypothetical optimality principles (e.g., maximization of growth) and precise gene expression (e.g., fold changes) for phenotype prediction. This paper introduces OptDesign, a new two-step strain design strategy. In the first step, OptDesign selects regulation candidates that have a noticeable flux difference between the wild type and production strains. In the second step, it computes optimal design strategies with limited manipulations (combining regulation and knockout), leading to high biochemical production. The usefulness and capabilities of OptDesign are demonstrated for the production of three biochemicals in Escherichia coli using the latest genome-scale metabolic model iML1515, showing highly consistent results with previous studies while suggesting new manipulations to boost strain performance. The source code is available at https://github.com/chang88ye/OptDesign.

Examining organic acid production potential and growth-coupled strategies in Issatchenkia orientalis using constraint-based modeling

Growth-coupling product formation can facilitate strain stability by aligning industrial objectives with biological fitness. Organic acids make up many building block chemicals that can be produced from sugars obtainable from renewable biomass. Issatchenkia orientalis is a yeast strain tolerant to acidic conditions and is thus a promising host for industrial production of organic acids. Here, we use constraint-based methods to assess the potential of computationally designing growth-coupled production strains for I. orientalis that produce 22 different organic acids under aerobic or microaerobic conditions. We explore native and engineered pathways using glucose or xylose as the carbon substrates as proxy constituents of hydrolyzed biomass. We identified growth-coupled production strategies for 37 of the substrate-product pairs, with 15 pairs achieving production for any growth rate. We systematically assess the strain design solutions and categorize the underlying principles involved.

Enhancement of Ethanol Production Using a Hybrid of Firefly Algorithm and Dynamic Flux Balance Analysis

Many high-demand industrial products are generated by microorganisms, including fuels, food, vitamins, and other chemicals. Metabolic engineering is the method of circumventing cellular control to manufacture a desirable product or to create a new product that the host cells do not normally need to produce. One of the objectives of microorganism metabolic engineering is to maximise the production of a desired product. However, owing to the structure of the regulatory cellular and metabolic network, identifying specific genes to be knocked out is difficult. The development of optimization algorithms often confronts issues such as easily trapping in local maxima and handling multivariate and multimodal functions inefficiently. To predict the gene knockout list that can generate high yields of desired product, a hybrid of Firefly Algorithm and Dynamic Flux Balance Analysis (FADFBA) is proposed. This paper focuses on the ethanol production of Escherichia coli (E. coli). The findings of the experiments include gene lists, ethanol production, growth rate, and the performance of FADFBA.

Metabolic systems biology and multi-omics of cyanobacteria: Perspectives and future directions

Cyanobacteria are oxygenic photoautotrophs whose metabolism contains key biochemical pathways to fix atmospheric CO₂ and synthesize various metabolites. The development of bioengineering tools has enabled the manipulation of cyanobacterial chassis to produce various valuable bioproducts photosynthetically. However, effective utilization of cyanobacteria as photosynthetic cell factories needs a detailed understanding of their metabolism and its interaction with other cellular processes. Implementing systems and synthetic biology tools has generated a wealth of information on various metabolic pathways. However, to design effective engineering strategies for further improvement in growth, photosynthetic efficiency, and enhanced production of target biochemicals, in-depth knowledge of their carbon/nitrogen metabolism, pathway fluxe distribution, genetic regulation and integrative analyses are necessary. In this review, we discuss the recent advances in the development of genome-scale metabolic models (GSMMs), omics analyses (metabolomics, transcriptomics, proteomics, fluxomics), and integrative modeling approaches to showcase the current understanding of cyanobacterial metabolism.

Engineering Bacteroides thetaiotaomicron to produce non-native butyrate based on a genome-scale metabolic model-guided design

Bacteroides thetaiotaomicron represents a major symbiont of the human gut microbiome that is increasingly viewed as a promising candidate strain for microbial therapeutics. Here, we engineer B. thetaiotaomicron for heterologous production of non-native butyrate as a proof-of-concept biochemical at therapeutically relevant concentrations. Since B. thetaiotaomicron is not a natural producer of butyrate, we heterologously expressed a butyrate biosynthetic pathway in the strain, which led to the production of butyrate at the final concentration of 12 mg/L in a rich medium. Further optimization of butyrate production was achieved by a round of metabolic engineering guided by an expanded genome-scale metabolic model (GEM) of B. thetaiotaomicron. The in silico knock-out simulation of the expanded model showed that pta and ldhD were the potent knock-out targets to enhance butyrate production. The maximum titer and specific productivity of butyrate in the pta-ldhD double knockout mutant increased by nearly 3.4 and 4.8 folds, respectively. To our knowledge, this is the first engineering attempt that enabled butyrate production from a non-butyrate producing commensal B. thetaiotaomicron. The study also highlights that B. thetaiotaomicron can serve as an effective strain for live microbial therapeutics in human.

Systematizing the different notions of growth-coupled product synthesis and a single framework for computing corresponding strain designs

A widely used design principle for metabolic engineering of microorganisms aims to introduce interventions that enforce growth-coupled product synthesis such that the product of interest becomes a (mandatory) by-product of growth. However, different variants and partially contradicting notions of growth-coupled production (GCP) exist. Herein, we propose an ontology for the different degrees of GCP and clarify their relationships. Ordered by coupling degree, we distinguish four major classes: potentially, weakly, and directionally growth-coupled production (pGCP, wGCP, dGCP) as well as substrate-uptake coupled production (SUCP). We then extend the framework of Minimal Cut Sets (MCS), previously used to compute dGCP and SUCP strain designs, to allow inclusion of implicit optimality constraints, a feature required to compute pGCP and wGCP designs. This extension closes the gap between MCS-based and bilevel-based strain design approaches and enables computation (and comparison) of designs for all GCP classes within a single framework. By computing GCP strain designs for a range of products, we illustrate the hierarchical relationships between the different coupling degrees. We find that feasibility of coupling is not affected by the chosen GCP degree and that strongest coupling (SUCP) requires often only one or two more interventions than wGCP and dGCP. Finally, we show that the principle of coupling can be generalized to couple product synthesis with other cellular functions than growth, for example, with net ATP formation. This work provides important theoretical results and algorithmic developments and a unified terminology for computational strain design based on GCP.

Plug-in repressor library for precise regulation of metabolic flux in Escherichia coli

In metabolic engineering, enhanced production of value-added chemicals requires precise flux control between growth-essential competing and production pathways. Although advances in synthetic biology have facilitated the exploitation of a number of genetic elements for precise flux control, their use requires expensive inducers, or more importantly, needs complex and time-consuming processes to design and optimize appropriate regulator components, case-by-case. To overcome this issue, we devised the plug-in repressor libraries for target-specific flux control, in which expression levels of the repressors were diversified using degenerate 5' untranslated region (5' UTR) sequences employing the UTR Library Designer. After we validated a wide expression range of the repressor libraries, they were applied to improve the production of lycopene from glucose and 3-hydroxypropionic acid (3-HP) from acetate in Escherichia coli via precise flux rebalancing to enlarge precursor pools. Consequently, we successfully achieved optimal carbon fluxes around the precursor nodes for efficient production. The most optimized strains were observed to produce 2.59 g/L of 3-HP and 11.66 mg/L of lycopene, which were improved 16.5-fold and 2.82-fold, respectively, compared to those produced by the parental strains. These results indicate that carbon flux rebalancing using the plug-in library is a powerful strategy for efficient production of value-added chemicals in E. coli.

Escherichia coli as a platform microbial host for systems metabolic engineering

Bio-based production of industrially important chemicals and materials from non-edible and renewable biomass has become increasingly important to resolve the urgent worldwide issues including climate change. Also, bio-based production, instead of chemical synthesis, of food ingredients and natural products has gained ever increasing interest for health benefits. Systems metabolic engineering allows more efficient development of microbial cell factories capable of sustainable, green, and human-friendly production of diverse chemicals and materials. Escherichia coli is unarguably the most widely employed host strain for the bio-based production of chemicals and materials. In the present paper, we review the tools and strategies employed for systems metabolic engineering of E. coli. Next, representative examples and strategies for the production of chemicals including biofuels, bulk and specialty chemicals, and natural products are discussed, followed by discussion on materials including polyhydroxyalkanoates (PHAs), proteins, and nanomaterials. Lastly, future perspectives and challenges remaining for systems metabolic engineering of E. coli are discussed.

Genetic and metabolic engineering challenges of C1-gas fermenting acetogenic chassis organisms

Unabated mining and utilisation of petroleum and petroleum resources and their conversion to essential fuels and chemicals have drastic environmental consequences, contributing to global warming and climate change. In addition, fossil fuels are finite resources, with a fast-approaching shortage. Accordingly, research efforts are increasingly focusing on developing sustainable alternatives for chemicals and fuels production. In this context, bioprocesses, relying on microorganisms, have gained particular interest. For example, acetogens use the Wood-Ljungdahl pathway to grow on single carbon C1-gases (CO2 and CO) as their sole carbon source and produce valuable products such as acetate or ethanol. These autotrophs can, therefore, be exploited for large-scale fermentation processes to produce industrially relevant chemicals from abundant greenhouse gases. In addition, genetic tools have recently been developed to improve these chassis organisms through synthetic biology approaches. This review will focus on the challenges of genetically and metabolically modifying acetogens. It will first discuss the physical and biochemical obstacles complicating successful DNA transfer in these organisms. Current genetic tools developed for several acetogens, crucial for strain engineering to consolidate and expand their catalogue of products, will then be described. Recent tool applications for metabolic engineering purposes to allow redirection of metabolic fluxes or production of non-native compounds will lastly be covered.

Automated engineering of synthetic metabolic pathways for efficient biomanufacturing

Metabolic engineering involves the engineering and optimization of processes from single-cell to fermentation in order to increase production of valuable chemicals for health, food, energy, materials and others. A systems approach to metabolic engineering has gained traction in recent years thanks to advances in strain engineering, leading to an accelerated scaling from rapid prototyping to industrial production. Metabolic engineering is nowadays on track towards a truly manufacturing technology, with reduced times from conception to production enabled by automated protocols for DNA assembly of metabolic pathways in engineered producer strains. In this review, we discuss how the success of the metabolic engineering pipeline often relies on retrobiosynthetic protocols able to identify promising production routes and dynamic regulation strategies through automated biodesign algorithms, which are subsequently assembled as embedded integrated genetic circuits in the host strain. Those approaches are orchestrated by an experimental design strategy that provides optimal scheduling planning of the DNA assembly, rapid prototyping and, ultimately, brings forward an accelerated Design-Build-Test-Learn cycle and the overall optimization of the biomanufacturing process. Achieving such a vision will address the increasingly compelling demand in our society for delivering valuable biomolecules in an affordable, inclusive and sustainable bioeconomy.

NIHBA: a network interdiction approach for metabolic engineering design

MOTIVATION

Flux balance analysis (FBA) based bilevel optimization has been a great success in redesigning metabolic networks for biochemical overproduction. To date, many computational approaches have been developed to solve the resulting bilevel optimization problems. However, most of them are of limited use due to biased optimality principle, poor scalability with the size of metabolic networks, potential numeric issues or low quantity of design solutions in a single run.

RESULTS

Here, we have employed a network interdiction model free of growth optimality assumptions, a special case of bilevel optimization, for computational strain design and have developed a hybrid Benders algorithm (HBA) that deals with complicating binary variables in the model, thereby achieving high efficiency without numeric issues in search of best design strategies. More importantly, HBA can list solutions that meet users' production requirements during the search, making it possible to obtain numerous design strategies at a small runtime overhead (typically ∼1 h, e.g. studied in this article).

AVAILABILITY AND IMPLEMENTATION

Source code implemented in the MATALAB Cobratoolbox is freely available at https://github.com/chang88ye/NIHBA.

CONTACT

math4neu@gmail.com or natalio.krasnogor@ncl.ac.uk.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Modeling a co-culture of Clostridium autoethanogenum and Clostridium kluyveri to increase syngas conversion to medium-chain fatty-acids

•

We model a co-culture of C. autoethanogenum and C. kluyveri for syngas fermentation.

•

Biomass species ratio affects ethanol and acetate profiles in the co-culture.

•

The model predicts that addition of succinate increases caproate production.

•

Genetic interventions  in C. autoethanogenum could increase caproate production.

Microbial fermentation of synthesis gas (syngas) is becoming more attractive for sustainable production of commodity chemicals. To date, syngas fermentation focuses mainly on the use of Clostridium species for the production of small organic molecules such as ethanol and acetate. The co-cultivation of syngas-fermenting microorganisms with chain-elongating bacteria can expand the range of possible products, allowing, for instance, the production of medium-chain fatty acids (MCFA) and alcohols from syngas. To explore these possibilities, we report herein a genome-scale, constraint-based metabolic model to describe growth of a co-culture of Clostridium autoethanogenum and Clostridium kluyveri on syngas for the production of valuable compounds. Community flux balance analysis was used to gain insight into the metabolism of the two strains and their interactions, and to reveal potential strategies enabling production of butyrate and hexanoate.

The model suggests that one strategy to optimize the production of medium-chain fatty-acids from syngas would be the addition of succinate. According to the prediction, addition of succinate would increase the pool of crotonyl-CoA and the ethanol/acetate uptake ratio in C. kluyveri, resulting in a flux of up to 60% of electrons into hexanoate. Another potential way to further optimize butyrate and hexanoate production would be an increase of C. autoethanogenum ethanol production. Blocking either acetaldehyde dehydrogenase or formate dehydrogenase (ferredoxin) activity or formate transport, in the C. autoethanogenum metabolic model could potentially lead to an up to 150% increase in ethanol production.

An extended and generalized framework for the calculation of metabolic intervention strategies based on minimal cut sets

The concept of minimal cut sets (MCS) provides a flexible framework for analyzing properties of metabolic networks and for computing metabolic intervention strategies. In particular, it has been used to support the targeted design of microbial strains for bio-based production processes. Herein we present a number of major extensions that generalize the existing MCS approach and broaden its scope for applications in metabolic engineering. We first introduce a modified approach to integrate gene-protein-reaction associations (GPR) in the metabolic network structure for the computation of gene-based intervention strategies. In particular, we present a set of novel compression rules for GPR associations, which effectively speedup the computation of gene-based MCS by a factor of up to one order of magnitude. These rules are not specific for MCS and as well applicable to other computational strain design methods. Second, we enhance the MCS framework by allowing the definition of multiple target (undesired) and multiple protected (desired) regions. This enables precise tailoring of the metabolic solution space of the designed strain with unlimited flexibility. Together with further generalizations such as individual cost factors for each intervention, direct combinations of reaction/gene deletions and additions as well as the possibility to search for substrate co-feeding strategies, the scope of the MCS framework could be broadly extended. We demonstrate the applicability and performance benefits of the described developments by computing (gene-based) Escherichia coli strain designs for the bio-based production of 2,3-butanediol, a chemical, that has recently received much attention in the field of metabolic engineering. With our extended framework, we could identify promising strain designs that were formerly unpredictable, including those based on substrate co-feeding.

Characterizing and ranking computed metabolic engineering strategies

MOTIVATION

The computer-aided design of metabolic intervention strategies has become a key component of an integrated metabolic engineering approach and a broad range of methods and algorithms has been developed for this task. Many of these algorithms enforce coupling of growth with product synthesis and may return thousands of possible intervention strategies from which the most suitable strategy must then be selected.

RESULTS

This work focuses on how to evaluate and rank, in a meaningful way, a given pool of computed metabolic engineering strategies for growth-coupled product synthesis. Apart from straightforward criteria, such as a preferably small number of necessary interventions, a reasonable growth rate and a high product yield, we present several new criteria useful to pick the most suitable intervention strategy. Among others, we investigate the robustness of the intervention strategies by searching for metabolites that may disrupt growth coupling when accumulated or secreted and by checking whether the interventions interrupt pathways at their origin (preferable) or at downstream steps. We also assess thermodynamic properties of the pathway(s) favored by the intervention strategy. Furthermore, strategies that have a significant overlap with alternative solutions are ranked higher because they provide flexibility in implementation. We also introduce the notion of equivalence classes for grouping intervention strategies with identical solution spaces. Our ranking procedure involves in total ten criteria and we demonstrate its applicability by assessing knockout-based intervention strategies computed in a genome-scale model of E.coli for the growth-coupled synthesis of l-methionine and of the heterologous product 1,4-butanediol.

AVAILABILITY AND IMPLEMENTATION

The MATLAB scripts that were used to characterize and rank the example intervention strategies are available at http://www2.mpi-magdeburg.mpg.de/projects/cna/etcdownloads.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Key Challenges in Designing CHO Chassis Platforms

Following the success of and the high demand for recombinant protein-based therapeutics during the last 25 years, the pharmaceutical industry has invested significantly in the development of novel treatments based on biologics. Mammalian cells are the major production systems for these complex biopharmaceuticals, with Chinese hamster ovary (CHO) cell lines as the most important players. Over the years, various engineering strategies and modeling approaches have been used to improve microbial production platforms, such as bacteria and yeasts, as well as to create pre-optimized chassis host strains. However, the complexity of mammalian cells curtailed the optimization of these host cells by metabolic engineering. Most of the improvements of titer and productivity were achieved by media optimization and large-scale screening of producer clones. The advances made in recent years now open the door to again consider the potential application of systems biology approaches and metabolic engineering also to CHO. The availability of a reference genome sequence, genome-scale metabolic models and the growing number of various “omics” datasets can help overcome the complexity of CHO cells and support design strategies to boost their production performance. Modular design approaches applied to engineer industrially relevant cell lines have evolved to reduce the time and effort needed for the generation of new producer cells and to allow the achievement of desired product titers and quality. Nevertheless, important steps to enable the design of a chassis platform similar to those in use in the microbial world are still missing. In this review, we highlight the importance of mammalian cellular platforms for the production of biopharmaceuticals and compare them to microbial platforms, with an emphasis on describing novel approaches and discussing still open questions that need to be resolved to reach the objective of designing enhanced modular chassis CHO cell lines.

Enhanced succinic acid production by Mannheimia employing optimal malate dehydrogenase

Succinic acid (SA), a dicarboxylic acid of industrial importance, can be efficiently produced by metabolically engineered Mannheimia succiniciproducens. Malate dehydrogenase (MDH) is one of the key enzymes for SA production, but has not been well characterized. Here we report biochemical and structural analyses of various MDHs and development of hyper-SA producing M. succiniciproducens by introducing the best MDH. Corynebacterium glutamicum MDH (CgMDH) shows the highest specific activity and least substrate inhibition, whereas M. succiniciproducens MDH (MsMDH) shows low specific activity at physiological pH and strong uncompetitive inhibition toward oxaloacetate (ki of 67.4 and 588.9 μM for MsMDH and CgMDH, respectively). Structural comparison of the two MDHs reveals a key residue influencing the specific activity and susceptibility to substrate inhibition. A high-inoculum fed-batch fermentation of the final strain expressing cgmdh produces 134.25 g L^-1 of SA with the maximum productivity of 21.3 g L^-1 h^-1, demonstrating the importance of enzyme optimization in strain development.

MOMO - multi-objective metabolic mixed integer optimization: application to yeast strain engineering

BACKGROUND

In this paper, we explore the concept of multi-objective optimization in the field of metabolic engineering when both continuous and integer decision variables are involved in the model. In particular, we propose a multi-objective model that may be used to suggest reaction deletions that maximize and/or minimize several functions simultaneously. The applications may include, among others, the concurrent maximization of a bioproduct and of biomass, or maximization of a bioproduct while minimizing the formation of a given by-product, two common requirements in microbial metabolic engineering.

RESULTS

Production of ethanol by the widely used cell factory Saccharomyces cerevisiae was adopted as a case study to demonstrate the usefulness of the proposed approach in identifying genetic manipulations that improve productivity and yield of this economically highly relevant bioproduct. We did an in vivo validation and we could show that some of the predicted deletions exhibit increased ethanol levels in comparison with the wild-type strain.

CONCLUSIONS

The multi-objective programming framework we developed, called MOMO, is open-source and uses POLYSCIP (Available at http://polyscip.zib.de/). as underlying multi-objective solver. MOMO is available at http://momo-sysbio.gforge.inria.fr.

Distributed flux balance analysis simulations of serial biomass fermentation by two organisms

Intelligent biorefinery design that addresses both the composition of the biomass feedstock as well as fermentation microorganisms could benefit from dedicated tools for computational simulation and computer-assisted optimization. Here we present the BioLego Vn2.0 framework, based on Microsoft Azure Cloud, which supports large-scale simulations of biomass serial fermentation processes by two different organisms. BioLego enables the simultaneous analysis of multiple fermentation scenarios and the comparison of fermentation potential of multiple feedstock compositions. Thanks to the effective use of cloud computing it further allows resource intensive analysis and exploration of media and organism modifications. We use BioLego to obtain biological and validation results, including (1) exploratory search for the optimal utilization of corn biomasses-corn cobs, corn fiber and corn stover-in fermentation biorefineries; (2) analysis of the possible effects of changes in the composition of K. alvarezi biomass on the ethanol production yield in an anaerobic two-step process (S. cerevisiae followed by E. coli); (3) analysis of the impact, on the estimated ethanol production yield, of knocking out single organism reactions either in one or in both organisms in an anaerobic two-step fermentation process of Ulva sp. into ethanol (S. cerevisiae followed by E. coli); and (4) comparison of several experimentally measured ethanol fermentation rates with the predictions of BioLego.

Bioinoculant capability enhancement through metabolomics and systems biology approaches

Bioinoculants are eco-friendly microorganisms, and their products are utilized for improving the potential of soil and fulfill the nutrients requirement for the host plant. The agricultural yield has increased due to the use of bioinoculants over chemical-based fertilizers, and thus it generates interest in understanding the innovation process by various methods. By gene-editing tool, the desired gene product can be changed for engineered microbial inoculants. We have also described various modern biotechnological tools like constraint-based modeling, OptKnock, flux balance analysis and modeling of the biological network for enhancing the bioinoculant capability. These fluxes give the fascinating perception of the metabolic network in the absence of comprehensive kinetic information. These tools also help in the stimulation of the metabolic networks by incorporation of enzyme-encoding genes. The present review explains the use of systems biology and gene-editing tools for improving the capability of bioinoculants. Moreover, this review also emphasizes on the challenges and future perspective of systems biology and its multidisciplinary facets.

Metabolic model guided strain design of cyanobacteria

Cyanobacteria are oxygenic photoautotrophs that serve as potential platforms for the production of biochemicals from cheap and renewable raw materials - sunlight, water, and carbon dioxide. Systems level analysis of the metabolic network of these organisms could enable the successful engineering of these organisms for the enhanced production of target chemicals. Metabolic modeling techniques including both stoichiometric and kinetic modeling with a genome-wide coverage enable a global assessment of metabolic capabilities. Recent studies guided by such modeling techniques have engineered strains for the enhanced production of valuable chemicals such as ethanol, n-butanol, 1,3-propanediol, glycerol, limonene, and isoprene from CO₂.

Determination of growth-coupling strategies and their underlying principles

BACKGROUND

Metabolic coupling of product synthesis and microbial growth is a prominent approach for maximizing production performance. Growth-coupling (GC) also helps stabilizing target production and allows the selection of superior production strains by adaptive laboratory evolution. To support the implementation of growth-coupling strain designs, we seek to identify biologically relevant, metabolic principles that enforce strong growth-coupling on the basis of reaction knockouts.

RESULTS

We adapted an established bilevel programming framework to maximize the minimally guaranteed production rate at a fixed, medium growth rate. Using this revised formulation, we identified various GC intervention strategies for metabolites of the central carbon metabolism, which were examined for GC generating principles under diverse conditions. Curtailing the metabolism to render product formation an essential carbon drain was identified as one major strategy generating strong coupling of metabolic activity and target synthesis. Impeding the balancing of cofactors and protons in the absence of target production was the underlying principle of all other strategies and further increased the GC strength of the aforementioned strategies.

CONCLUSION

Maximizing the minimally guaranteed production rate at a medium growth rate is an attractive principle for the identification of strain designs that couple growth to target metabolite production. Moreover, it allows for controlling the inevitable compromise between growth coupling strength and the retaining of microbial viability. With regard to the corresponding metabolic principles, generating a dependency between the supply of global metabolic cofactors and product synthesis appears to be advantageous in enforcing strong GC for any metabolite. Deriving such strategies manually, is a hard task, due to which we suggest incorporating computational metabolic network analyses in metabolic engineering projects seeking to determine GC strain designs.

Seeing the wood for the trees: a forest of methods for optimization and omic-network integration in metabolic modelling

Metabolic modelling has entered a mature phase with dozens of methods and software implementations available to the practitioner and the theoretician. It is not easy for a modeller to be able to see the wood (or the forest) for the trees. Driven by this analogy, we here present a 'forest' of principal methods used for constraint-based modelling in systems biology. This provides a tree-based view of methods available to prospective modellers, also available in interactive version at http://modellingmetabolism.net, where it will be kept updated with new methods after the publication of the present manuscript. Our updated classification of existing methods and tools highlights the most promising in the different branches, with the aim to develop a vision of how existing methods could hybridize and become more complex. We then provide the first hands-on tutorial for multi-objective optimization of metabolic models in R. We finally discuss the implementation of multi-view machine learning approaches in poly-omic integration. Throughout this work, we demonstrate the optimization of trade-offs between multiple metabolic objectives, with a focus on omic data integration through machine learning. We anticipate that the combination of a survey, a perspective on multi-view machine learning and a step-by-step R tutorial should be of interest for both the beginner and the advanced user.

Genome-driven cell engineering review: in vivo and in silico metabolic and genome engineering

Producing 'designer cells' with specific functions is potentially feasible in the near future. Recent developments, including whole-cell models, genome design algorithms and gene editing tools, have advanced the possibility of combining biological research and mathematical modelling to further understand and better design cellular processes. In this review, we will explore computational and experimental approaches used for metabolic and genome design. We will highlight the relevance of modelling in this process, and challenges associated with the generation of quantitative predictions about cell behaviour as a whole: although many cellular processes are well understood at the subsystem level, it has proved a hugely complex task to integrate separate components together to model and study an entire cell. We explore these developments, highlighting where computational design algorithms compensate for missing cellular information and underlining where computational models can complement and reduce lab experimentation. We will examine issues and illuminate the next steps for genome engineering.

Comparison of pathway analysis and constraint-based methods for cell factory design

BACKGROUND

Computational strain optimisation methods (CSOMs) have been successfully used to exploit genome-scale metabolic models, yielding strategies useful for allowing compound overproduction in metabolic cell factories. Minimal cut sets are particularly interesting since their definition allows searching for intervention strategies that impose strong growth-coupling phenotypes, and are not subject to optimality bias when compared with simulation-based CSOMs. However, since both types of methods have different underlying principles, they also imply different ways to formulate metabolic engineering problems, posing an obstacle when comparing their outputs.

RESULTS

In this work, we perform an in-depth analysis of potential strategies that can be obtained with both methods, providing a critical comparison of performance, robustness, predicted phenotypes as well as strategy structure and size. To this end, we devised a pipeline including enumeration of strategies from evolutionary algorithms (EA) and minimal cut sets (MCS), filtering and flux analysis of predicted mutants to optimize the production of succinic acid in Saccharomyces cerevisiae. We additionally attempt to generalize problem formulations for MCS enumeration within the context of growth-coupled product synthesis. Strategies from evolutionary algorithms show the best compromise between acceptable growth rates and compound overproduction. However, constrained MCSs lead to a larger variety of phenotypes with several degrees of growth-coupling with production flux. The latter have proven useful in revealing the importance, in silico, of the gamma-aminobutyric acid shunt and manipulation of cofactor pools in growth-coupled designs for succinate production, mechanisms which have also been touted as potentially useful for metabolic engineering.

CONCLUSIONS

The two main groups of CSOMs are valuable for finding growth-coupled mutants. Despite the limitations in maximum growth rates and large strategy sizes, MCSs help uncover novel mechanisms for compound overproduction and thus, analyzing outputs from both methods provides a richer overview on strategies that can be potentially carried over in vivo.

OptCouple: Joint simulation of gene knockouts, insertions and medium modifications for prediction of growth-coupled strain designs

Biological production of chemicals is an attractive alternative to petrochemical-based production, due to advantages in environmental impact and the spectrum of feasible targets. However, engineering microbial strains to overproduce a compound of interest can be a long, costly and painstaking process. If production can be coupled to cell growth it is possible to use adaptive laboratory evolution to increase the production rate. Strategies for coupling production to growth, however, are often not trivial to find. Here we present OptCouple, a constraint-based modeling algorithm to simultaneously identify combinations of gene knockouts, insertions and medium supplements that lead to growth-coupled production of a target compound. We validated the algorithm by showing that it can find novel strategies that are growth-coupled in silico for a compound that has not been coupled to growth previously, as well as reproduce known growth-coupled strain designs for two different target compounds. Furthermore, we used OptCouple to construct an alternative design with potential for higher production. We provide an efficient and easy-to-use implementation of the OptCouple algorithm in the cameo Python package for computational strain design.

OptRAM: In-silico strain design via integrative regulatory-metabolic network modeling

The ultimate goal of metabolic engineering is to produce desired compounds on an industrial scale in a cost effective manner. To address challenges in metabolic engineering, computational strain optimization algorithms based on genome-scale metabolic models have increasingly been used to aid in overproducing products of interest. However, most of these strain optimization algorithms utilize a metabolic network alone, with few approaches providing strategies that also include transcriptional regulation. Moreover previous integrated approaches generally require a pre-existing regulatory network. In this study, we developed a novel strain design algorithm, named OptRAM (Optimization of Regulatory And Metabolic Networks), which can identify combinatorial optimization strategies including overexpression, knockdown or knockout of both metabolic genes and transcription factors. OptRAM is based on our previous IDREAM integrated network framework, which makes it able to deduce a regulatory network from data. OptRAM uses simulated annealing with a novel objective function, which can ensure a favorable coupling between desired chemical and cell growth. The other advance we propose is a systematic evaluation metric of multiple solutions, by considering the essential genes, flux variation, and engineering manipulation cost. We applied OptRAM to generate strain designs for succinate, 2,3-butanediol, and ethanol overproduction in yeast, which predicted high minimum predicted target production rate compared with other methods and previous literature values. Moreover, most of the genes and TFs proposed to be altered by OptRAM in these scenarios have been validated by modification of the exact genes or the target genes regulated by the TFs, for overproduction of these desired compounds by in vivo experiments cataloged in the LASER database. Particularly, we successfully validated the predicted strain optimization strategy for ethanol production by fermentation experiment. In conclusion, OptRAM can provide a useful approach that leverages an integrated transcriptional regulatory network and metabolic network to guide metabolic engineering applications.

The genetic basis for adaptation of model-designed syntrophic co-cultures

Understanding the fundamental characteristics of microbial communities could have far reaching implications for human health and applied biotechnology. Despite this, much is still unknown regarding the genetic basis and evolutionary strategies underlying the formation of viable synthetic communities. By pairing auxotrophic mutants in co-culture, it has been demonstrated that viable nascent E. coli communities can be established where the mutant strains are metabolically coupled. A novel algorithm, OptAux, was constructed to design 61 unique multi-knockout E. coli auxotrophic strains that require significant metabolite uptake to grow. These predicted knockouts included a diverse set of novel non-specific auxotrophs that result from inhibition of major biosynthetic subsystems. Three OptAux predicted non-specific auxotrophic strains—with diverse metabolic deficiencies—were co-cultured with an L-histidine auxotroph and optimized via adaptive laboratory evolution (ALE). Time-course sequencing revealed the genetic changes employed by each strain to achieve higher community growth rates and provided insight into mechanisms for adapting to the syntrophic niche. A community model of metabolism and gene expression was utilized to predict the relative community composition and fundamental characteristics of the evolved communities. This work presents new insight into the genetic strategies underlying viable nascent community formation and a cutting-edge computational method to elucidate metabolic changes that empower the creation of cooperative communities.

Many basic characteristics underlying the establishment of cooperative growth in bacterial communities have not been studied in detail. The presented work sought to understand the adaptation of syntrophic communities by first employing a new computational method to generate a comprehensive catalog of E. coli auxotrophic mutants. Many of the knockouts in the catalog had the predicted effect of disabling a major biosynthetic process. As a result, these strains were predicted to be capable of growing when supplemented with many different individual metabolites (i.e., a non-specific auxotroph), but the strains would require a high amount of metabolic cooperation to grow in community. Three such non-specific auxotroph mutants from this catalog were co-cultured with a proven auxotrophic partner in vivo and evolved via adaptive laboratory evolution. In order to successfully grow, each strain in co-culture had to evolve under a pressure to grow cooperatively in its new niche. The non-specific auxotrophs further had to adapt to significant homeostatic changes in cell’s metabolic state caused by knockouts in metabolic genes. The genomes of the successfully growing communities were sequenced, thus providing unique insights into the genetic changes accompanying the formation and optimization of the viable communities. A computational model was further developed to predict how finite protein availability, a fundamental constraint on cell metabolism, could impact the composition of the community (i.e., the relative abundances of each community member).

Use of genome-scale models to get new insights into the marine actinomycete genus Salinispora

BACKGROUND

There is little published regarding metabolism of Salinispora species. In continuation with efforts performed towards this goal, this study is focused on new insights into the metabolism of the three-identified species of Salinispora using constraints-based modeling. At present, only one manually curated genome-scale metabolic model (GSM) for Salinispora tropica strain CNB-440T has been built despite the role of Salinispora strains in drug discovery.

RESULTS

Here, we updated, and expanded the scope of the model of Salinispora tropica CNB-440T, and GSMs were constructed for two sequenced type strains covering the three-identified species. We also constructed a Salinispora core model that contains the genes shared by 93 sequenced strains and a few non-conserved genes associated with essential reactions. The models predicted no auxotrophies for essential amino acids, which was corroborated experimentally using a defined minimal medium (DMM). Experimental observations suggest possible sulfur accumulation. The Core metabolic content shows that the biosynthesis of specialised metabolites is the less conserved subsystem. Sets of reactions were analyzed to explore the differences between the reconstructions. Unique reactions associated to each GSM were mainly due to genome sequence data except for the ST-CNB440 reconstruction. In this case, additional reactions were added from experimental evidence. This reveals that by reaction content the ST-CNB440 model is different from the other species models. The differences identified in reaction content between models gave rise to different functional predictions of essential nutrient usage by each species in DMM. Furthermore, models were used to evaluate in silico single gene knockouts under DMM and complex medium. Cluster analysis of these results shows that ST-CNB440, and SP-CNR114 models are more similar when considering predicted essential genes.

CONCLUSIONS

Models were built for each of the three currently identified Salinispora species, and a core model representing the conserved metabolic capabilities of Salinispora was constructed. Models will allow in silico metabolism studies of Salinispora strains, and help researchers to guide and increase the production of specialised metabolites. Also, models can be used as templates to build GSMs models of closely related organisms with high biotechnology potential.

Yeast Genome-Scale Metabolic Models for Simulating Genotype-Phenotype Relations

Understanding genotype-phenotype dependency is a universal aim for all life sciences. While the complete genotype-phenotype relations remain challenging to resolve, metabolic phenotypes are moving within the reach through genome-scale metabolic model simulations. Genome-scale metabolic models are available for commonly investigated yeasts, such as model eukaryote and domesticated fermentation species Saccharomyces cerevisiae, and automatic reconstruction methods facilitate obtaining models for any sequenced species. The models allow for investigating genotype-phenotype relations through simulations simultaneously considering the effects of nutrient availability, and redox and energy homeostasis in cells. Genome-scale models also offer frameworks for omics data integration to help to uncover how the translation of genotypes to the apparent phenotypes is regulated at different levels. In this chapter, we provide an overview of the yeast genome-scale metabolic models and the simulation approaches for using these models to interrogate genotype-phenotype relations. We review the methodological approaches according to the underlying biological reasoning in order to inspire formulating novel questions and applications that the genome-scale metabolic models could contribute to. Finally, we discuss current challenges and opportunities in the genome-scale metabolic model simulations.

Data Management in Computational Systems Biology: Exploring Standards, Tools, Databases, and Packaging Best Practices

Computational systems biology involves integrating heterogeneous datasets in order to generate models. These models can assist with understanding and prediction of biological phenomena. Generating datasets and integrating them into models involves a wide range of scientific expertise. As a result these datasets are often collected by one set of researchers, and exchanged with others researchers for constructing the models. For this process to run smoothly the data and models must be FAIR-findable, accessible, interoperable, and reusable. In order for data and models to be FAIR they must be structured in consistent and predictable ways, and described sufficiently for other researchers to understand them. Furthermore, these data and models must be shared with other researchers, with appropriately controlled sharing permissions, before and after publication. In this chapter we explore the different data and model standards that assist with structuring, describing, and sharing. We also highlight the popular standards and sharing databases within computational systems biology.

Identification of growth-coupled production strains considering protein costs and kinetic variability

Conversion of renewable biomass to useful molecules in microbial cell factories can be approached in a rational and systematic manner using constraint-based reconstruction and analysis. Filtering for high confidence in silico designs is critical because in vivo construction and testing of strains is expensive and time consuming. As such, a workflow was devised to analyze the robustness of growth-coupled production when considering the biosynthetic costs of the proteome and variability in enzyme kinetic parameters using a genome-scale model of metabolism and gene expression (ME-model). A collection of 2632 unfiltered knockout designs in Escherichia coli was evaluated by the workflow. A ME-model was used in the workflow to test the designs' growth-coupled production in addition to a less complex genome-scale metabolic model (M-model). The workflow identified 634 M-model growth-coupled designs which met the filtering criteria and 42 robust designs, which met growth-coupled production criteria using both M and ME-models. Knockouts were found to follow a pattern of controlling intermediate metabolite consumption such as pyruvate consumption and high flux subsystems such as glycolysis. Kinetic parameter sampling using the ME-model revealed how enzyme efficiency and pathway tradeoffs can affect growth-coupled production phenotypes.

Computing Minimum Reaction Modifications in a Boolean Metabolic Network

In metabolic network modification, we newly add enzymes or/and knock-out genes to maximize the biomass production with minimum side-effect. Although this problem has been studied for various problem settings via mathematical models including flux balance analysis, elementary mode, and Boolean models, some important problem settings still remain to be studied. In this paper, we consider the Boolean Reaction Modification (BRM) problem, where a host metabolic network and a reference metabolic network are given in the Boolean model. The host network initially produces some toxic compounds and cannot produce some necessary compounds, but the reference network can produce the necessary compounds, and we should minimize the total number of removed reactions from the host network and added reactions from the reference network so that the toxic compounds are not producible, but the necessary compounds are producible in the resulting host network. We developed integer linear programming (ILP)-based methods for BRM, and compared them with OptStrain and SimOptStrain. The results show that our method performed better for reducing the total number of added and removed reactions, while OptStrain and SimOptStrain performed better for optimizing the production of the target compound. Our developed software is freely available at "http://sunflower.kuicr.kyoto-u.ac.jp/~rogi/solBRM/solBRM.html ".

Genome-wide analysis of fitness data and its application to improve metabolic models

Background

Synthetic biology and related techniques enable genome scale high-throughput investigation of the effect on organism fitness of different gene knock-downs/outs and of other modifications of genomic sequence.

Results

We develop statistical and computational pipelines and frameworks for analyzing high throughput fitness data over a genome scale set of sequence variants. Analyzing data from a high-throughput knock-down/knock-out bacterial study, we investigate differences and determinants of the effect on fitness in different conditions. Comparing fitness vectors of genes, across tens of conditions, we observe that fitness consequences strongly depend on genomic location and more weakly depend on gene sequence similarity and on functional relationships. In analyzing promoter sequences, we identified motifs associated with conditions studied in bacterial media such as Casaminos, D-glucose, Sucrose, and other sugars and amino-acid sources.

We also use fitness data to infer genes associated with orphan metabolic reactions in the iJO1366 E. coli metabolic model. To do this, we developed a new computational method that integrates gene fitness and gene expression profiles within a given reaction network neighborhood to associate this reaction with a set of genes that potentially encode the catalyzing proteins. We then apply this approach to predict candidate genes for 107 orphan reactions in iJO1366. Furthermore - we validate our methodology with known reactions using a leave-one-out approach. Specifically, using top-20 candidates selected based on combined fitness and expression datasets, we correctly reconstruct 39.7% of the reactions, as compared to 33% based on fitness and to 26% based on expression separately, and to 4.02% as a random baseline. Our model improvement results include a novel association of a gene to an orphan cytosine nucleosidation reaction.

Conclusion

Our pipeline for metabolic modeling shows a clear benefit of using fitness data for predicting genes of orphan reactions. Along with the analysis pipelines we developed, it can be used to analyze similar high-throughput data.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2341-9) contains supplementary material, which is available to authorized users.