Computing complex metabolic intervention strategies using constrained minimal cut sets

Adding metabolic tasks to human GEM models to improve the study of gene targets and their associated toxicities

Genetic minimal cut sets (gMCS) are genes that must be deactivated simultaneously to avoid unwanted states in a metabolic model. The concept of gMCS can be applied to two different scenarios. First, it can be used to identify potential gene toxicities in generic or healthy cell models. Second, it can be used to develop genetic strategies to target cancer cells and prevent their proliferation. Up to now, gMCS have been evaluated using the traditional procedure of preventing biomass production. This paper proposes an additional way: using essential metabolic tasks, which any human cell should perform, to enlarge the set of unwanted states. Including this addition can significantly improve the study of toxicities and reveal targets that can be used to treat unhealthy cells. Excluding metabolic tasks can cause important information to be overlooked, which could impact the study's success. Regarding toxicities, using the generic Human model, the number of detected generic toxicities with metabolic tasks increases from 106 to 281 (136 gMCSs of length 1 and 49 of length 2). We have used the following context-specific models to evaluate specific toxicities in different healthy tissues: blood, pancreas, liver, heart, and kidney. Again, considering metabolic tasks, we have found new toxicities (lengths 1 and 2) whose inactivation could damage these healthy tissues.Our research strategy has been applied to identify new cancer drug targets in two myeloma cell lines. We obtained new therapeutic targets of lengths 1 and 2 for each cell line. After analyzing the data, we conclude that incorporating metabolic tasks into cancer models can reveal important therapeutic targets previously disregarded by the conventional method of inhibiting biomass production. This approach also improves the evaluation of potential drug toxicities.

Logic programming-based Minimal Cut Sets reveal consortium-level therapeutic targets for chronic wound infections

Minimal Cut Sets (MCSs) identify sets of reactions which, when removed from a metabolic network, disable certain cellular functions. The traditional search for MCSs within genome-scale metabolic models (GSMMs) targets cellular growth, identifies reaction sets resulting in a lethal phenotype if disrupted, and retrieves a list of corresponding gene, mRNA, or enzyme targets. Using the dual link between MCSs and Elementary Flux Modes (EFMs), our logic programming-based tool aspefm was able to compute MCSs of any size from GSMMs in acceptable run times. The tool demonstrated better performance when computing large-sized MCSs than the mixed-integer linear programming methods. We applied the new MCSs methodology to a medically-relevant consortium model of two cross-feeding bacteria, Staphylococcus aureus and Pseudomonas aeruginosa. aspefm constraints were used to bias the computation of MCSs toward exchanged metabolites that could complement lethal phenotypes in individual species. We found that interspecies metabolite exchanges could play an essential role in rescuing single-species growth, for instance inosine could complement lethal reaction knock-outs in the purine synthesis, glycolysis, and pentose phosphate pathways of both bacteria. Finally, MCSs were used to derive a list of promising enzyme targets for consortium-level therapeutic applications that cannot be circumvented via interspecies metabolite exchange.

Computational analysis into the potential of azo dyes as a feedstock for actinorhodin biosynthesis in Pseudomonas putida

Fermentation-based biosynthesis in synthetic biology relies heavily on sugar-derived feedstocks, a limited and carbon-intensive commodity. Unconventional feedstocks from less-noble sources such as waste are being utilized to produce high-value chemical products. Azo dyes, a major pollutant commonly discharged by food, textile, and pharmaceutical industries, present significant health and environmental risks. We explore the potential of engineering Pseudomonas putida KT2440 to utilize azo dyes as a substrate to produce a polyketide, actinorhodin (ACT). Using the constrained minimal cut sets (cMCS) approach, we identified metabolic interventions that optimize ACT biosynthesis and compare the growth-coupling solutions attainable on an azo dye compared to glucose. Our results predicted that azo dyes could perform better as a feedstock for ACT biosynthesis than glucose as it allowed growth-coupling regimes that are unfeasible with glucose and generated an 18.28% higher maximum ACT flux. By examining the flux distributions enabled in different carbon sources, we observed that carbon fluxes from aromatic compounds like azo dyes have a unique capability to leverage gluconeogenesis to support both growth and production of secondary metabolites that produce excess NADH. Carbon sources are commonly chosen based on the host organism, availability, cost, and environmental implications. We demonstrated that careful selection of carbon sources is also crucial to ensure that the resulting flux distribution is suitable for further metabolic engineering of microbial cell factories.

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 .

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.

Assessing and reducing phenotypic instability in cyanobacteria

Cyanobacteria have promising potential as sustainable cell factories. However, one challenge that is still largely unreported in scaling-up cyanobacteria bioproduction is phenotypic instability, where the emergence and selection of nonproducing cells leading to loss in production has longer evolutionary timescales to take place in industrial-scale bioreactors. Quantifying phenotypic instability early on in strain development allows researchers to make informed decisions on whether to proceed with scalable designs, or if present, devise countermeasures to reduce instability. One particularly effective strategy to mitigate instability is the use of genome-scale metabolic models to design growth-coupled production strains. In silico studies have predicted that creating certain cofactor imbalances or removing recycling reactions in cyanobacteria can be exploited to stably produce a wide variety of metabolites.

Trimming Gene Deletion Strategies for Growth-Coupled Production in Constraint-Based Metabolic Networks: TrimGdel

When simulating genome-scale metabolite production using constraint-based metabolic networks, it is often necessary to find gene deletion strategies which lead to growth-coupled production, which means that target metabolites are produced when cell growth is maximized. Existing methods are effective when the number of gene deletions is relatively small, but when the number of required gene deletions exceeds approximately 1% of whole genes, the time required for the calculation is often unfeasible. Therefore, a complementing algorithm that is effective even when the required number of gene deletions is approximately 1% to 5% of whole genes would be helpful because the number of deletable genes in a strain is increasing with advances in genetic engineering technology. In this study, the author developed an algorithm, TrimGdel, which first computes a strategy with many gene deletions that results in growth-coupled production and then gradually reduces the number of gene deletions while ensuring the original production rate and growth rate. The results of the computer experiments showed that TrimGdel can calculate stoichiometrically feasible gene deletion strategies, especially those whose sizes are 1 to 5% of whole genes, which lead to growth-coupled production of many target metabolites, which include useful vitamins such as biotin and pantothenate, for which existing methods could not.

Staphylococcus epidermidis RP62A’s Metabolic Network: Validation and Intervention Strategies

Increasingly, systems biology is gaining relevance in basic and applied research. The combination of computational biology with wet laboratory methods produces synergy that results in an exponential increase in knowledge of biological systems. The study of microorganisms such as Staphylococcus epidermidis RP62A enables the researcher to understand better their metabolic networks, which allows the design of effective strategies to treat infections caused by this species or others. S. epidermidis is the second most commoncause of infection in patients with joint implants, so treating its proliferation seems vital for public health. There are different approaches to the analysis of metabolic networks. Flux balance analysis (FBA) is one of the most widespread streams of research. It allows the study of large metabolic networks, the study their structural properties, the optimization of metabolic flux, and the search for intervention strategies to modify the state of the metabolic network. This work presents the validation of the Staphylococcus epidermidis RP62A metabolic network model elaborated by Díaz Calvo et al. Then, we elaborate further on the network analysis’s essential reactions. The full set of essential reactions (including a previously unobserved one) was computed, and we classified them into equivalence classes. Some proposals to intervene in the network and design knock-outs by studying minimal cut sets of small length are also introduced. In particular, minimal cut sets related to the medium (including exchange reactions associated with medium metabolites) have been computed. In this sense, the unique external MCS (composed of cysteine and sulfate ion) has been found, and all hybrid MCS (based on knocking out both internal and exchange reactions) of length two have also been computed. The paper also points out the possible importance of these new intervention strategies.

The landscape of receptor-mediated precision cancer combination therapy via a single-cell perspective

Mining a large cohort of single-cell transcriptomics data, here we employ combinatorial optimization techniques to chart the landscape of optimal combination therapies in cancer. We assume that each individual therapy can target any one of 1269 genes encoding cell surface receptors, which may be targets of CAR-T, conjugated antibodies or coated nanoparticle therapies. We find that in most cancer types, personalized combinations composed of at most four targets are then sufficient for killing at least 80% of tumor cells while sparing at least 90% of nontumor cells in the tumor microenvironment. However, as more stringent and selective killing is required, the number of targets needed rises rapidly. Emerging individual targets include PTPRZ1 for brain and head and neck cancers and EGFR in multiple tumor types. In sum, this study provides a computational estimate of the identity and number of targets needed in combination to target cancers selectively and precisely.

Microbial production of advanced biofuels

Concerns over climate change have necessitated a rethinking of our transportation infrastructure. One possible alternative to carbon-polluting fossil fuels is biofuels produced by engineered microorganisms that use a renewable carbon source. Two biofuels, ethanol and biodiesel, have made inroads in displacing petroleum-based fuels, but their uptake has been limited by the amounts that can be used in conventional engines and by their cost. Advanced biofuels that mimic petroleum-based fuels are not limited by the amounts that can be used in existing transportation infrastructure but have had limited uptake due to costs. In this Review, we discuss engineering metabolic pathways to produce advanced biofuels, challenges with substrate and product toxicity with regard to host microorganisms and methods to engineer tolerance, and the use of functional genomics and machine learning approaches to produce advanced biofuels and prospects for reducing their costs.

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.

Modelling microbial communities: Harnessing consortia for biotechnological applications

Microbes propagate and thrive in complex communities, and there are many benefits to studying and engineering microbial communities instead of single strains. Microbial communities are being increasingly leveraged in biotechnological applications, as they present significant advantages such as the division of labour and improved substrate utilisation. Nevertheless, they also present some interesting challenges to surmount for the design of efficient biotechnological processes. In this review, we discuss key principles of microbial interactions, followed by a deep dive into genome-scale metabolic models, focussing on a vast repertoire of constraint-based modelling methods that enable us to characterise and understand the metabolic capabilities of microbial communities. Complementary approaches to model microbial communities, such as those based on graph theory, are also briefly discussed. Taken together, these methods provide rich insights into the interactions between microbes and how they influence microbial community productivity. We finally overview approaches that allow us to generate and test numerous synthetic community compositions, followed by tools and methodologies that can predict effective genetic interventions to further improve the productivity of communities. With impending advancements in high-throughput omics of microbial communities, the stage is set for the rapid expansion of microbial community engineering, with a significant impact on biotechnological processes.

Modelling Cell Metabolism: A Review on Constraint-Based Steady-State and Kinetic Approaches

Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic therapy. Remarkable successes in these fields emerged from heuristics approaches, for instance, with the introduction of effective strategies for genetic modifications, drug developments and optimization of bioprocess management. However, heuristics approaches have showed significant shortcomings, such as to describe regulation of metabolic pathways and to extrapolate experimental conditions. In the specific case of bioprocess management, such shortcomings limit their capacity to increase product quality, while maintaining desirable productivity and reproducibility levels. For instance, since heuristics approaches are not capable of prediction of the cellular functions under varying experimental conditions, they may lead to sub-optimal processes. Also, such approaches used for bioprocess control often fail in regulating a process under unexpected variations of external conditions. Therefore, methodologies inspired by the systematic mathematical formulation of cell metabolism have been used to address such drawbacks and achieve robust reproducible results. Mathematical modelling approaches are effective for both the characterization of the cell physiology, and the estimation of metabolic pathways utilization, thus allowing to characterize a cell population metabolic behavior. In this article, we present a review on methodology used and promising mathematical modelling approaches, focusing primarily to investigate metabolic events and regulation. Proceeding from a topological representation of the metabolic networks, we first present the metabolic modelling approaches that investigate cell metabolism at steady state, complying to the constraints imposed by mass conservation law and thermodynamics of reactions reversibility. Constraint-based models (CBMs) are reviewed highlighting the set of assumed optimality functions for reaction pathways. We explore models simulating cell growth dynamics, by expanding flux balance models developed at steady state. Then, discussing a change of metabolic modelling paradigm, we describe dynamic kinetic models that are based on the mathematical representation of the mechanistic description of nonlinear enzyme activities. In such approaches metabolic pathway regulations are considered explicitly as a function of the activity of other components of metabolic networks and possibly far from the metabolic steady state. We have also assessed the significance of metabolic model parameterization in kinetic models, summarizing a standard parameter estimation procedure frequently employed in kinetic metabolic modelling literature. Finally, some optimization practices used for the parameter estimation are reviewed.

Opportunities and Challenges for Microbial Synthesis of Fatty Acid-Derived Chemicals (FACs)

Global warming and uneven distribution of fossil fuels worldwide concerns have spurred the development of alternative, renewable, sustainable, and environmentally friendly resources. From an engineering perspective, biosynthesis of fatty acid-derived chemicals (FACs) is an attractive and promising solution to produce chemicals from abundant renewable feedstocks and carbon dioxide in microbial chassis. However, several factors limit the viability of this process. This review first summarizes the types of FACs and their widely applications. Next, we take a deep look into the microbial platform to produce FACs, give an outlook for the platform development. Then we discuss the bottlenecks in metabolic pathways and supply possible solutions correspondingly. Finally, we highlight the most recent advances in the fast-growing model-based strain design for FACs biosynthesis.

Answer Set Programming for Computing Constraints-Based Elementary Flux Modes: Application to Escherichia coli Core Metabolism

Elementary Flux Modes (EFMs) provide a rigorous basis to systematically characterize the steady state, cellular phenotypes, as well as metabolic network robustness and fragility. However, the number of EFMs typically grows exponentially with the size of the metabolic network, leading to excessive computational demands, and unfortunately, a large fraction of these EFMs are not biologically feasible due to system constraints. This combinatorial explosion often prevents the complete analysis of genome-scale metabolic models. Traditionally, EFMs are computed by the double description method, an efficient algorithm based on matrix calculation; however, only a few constraints can be integrated into this computation. They must be monotonic with regard to the set inclusion of the supports; otherwise, they must be treated in post-processing and thus do not save computational time. We present aspefm, a hybrid computational tool based on Answer Set Programming (ASP) and Linear Programming (LP) that permits the computation of EFMs while implementing many different types of constraints. We apply our methodology to the Escherichia coli core model, which contains 226×106 EFMs. In considering transcriptional and environmental regulation, thermodynamic constraints, and resource usage considerations, the solution space is reduced to 1118 EFMs that can be computed directly with aspefm. The solution set, for E. coli growth on O2 gradients spanning fully aerobic to anaerobic, can be further reduced to four optimal EFMs using post-processing and Pareto front analysis.

Speeding up the core algorithm for the dual calculation of minimal cut sets in large metabolic networks

Background

The concept of minimal cut sets (MCS) has become an important mathematical framework for analyzing and (re)designing metabolic networks. However, the calculation of MCS in genome-scale metabolic models is a complex computational problem. The development of duality-based algorithms in the last years allowed the enumeration of thousands of MCS in genome-scale networks by solving mixed-integer linear problems (MILP). A recent advancement in this field was the introduction of the MCS² approach. In contrast to the Farkas-lemma-based dual system used in earlier studies, the MCS² approach employs a more condensed representation of the dual system based on the nullspace of the stoichiometric matrix, which, due to its reduced dimension, holds promise to further enhance MCS computations.

Results

In this work, we introduce several new variants and modifications of duality-based MCS algorithms and benchmark their effects on the overall performance. As one major result, we generalize the original MCS² approach (which was limited to blocking the operation of certain target reactions) to the most general case of MCS computations with arbitrary target and desired regions. Building upon these developments, we introduce a new MILP variant which allows maximal flexibility in the formulation of MCS problems and fully leverages the reduced size of the nullspace-based dual system. With a comprehensive set of benchmarks, we show that the MILP with the nullspace-based dual system outperforms the MILP with the Farkas-lemma-based dual system speeding up MCS computation with an averaged factor of approximately 2.5. We furthermore present several simplifications in the formulation of constraints, mainly related to binary variables, which further enhance the performance of MCS-related MILP. However, the benchmarks also reveal that some highly condensed formulations of constraints, especially on reversible reactions, may lead to worse behavior when compared to variants with a larger number of (more explicit) constraints and involved variables.

Conclusions

Our results further enhance the algorithmic toolbox for MCS calculations and are of general importance for theoretical developments as well as for practical applications of the MCS framework.

Software and Methods for Computational Flux Balance Analysis

As genetic engineering of organisms has grown easier and more precise, computational modeling of metabolic systems has played an increasingly important role in both guiding experimental interventions and in understanding the results of metabolic perturbations.

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.

Computing irreversible minimal cut sets in genome-scale metabolic networks via flux cone projection

MOTIVATION

Minimal cut sets (MCSs) for metabolic networks are sets of reactions which, if they are removed from the network, prevent a target reaction from carrying flux. To compute MCSs different methods exist, which may fail to find sufficiently many MCSs for larger genome-scale networks.

RESULTS

Here we introduce irreversible minimal cut sets (iMCSs). These are MCSs that consist of irreversible reactions only. The advantage of iMCSs is that they can be computed by projecting the flux cone of the metabolic network on the set of irreversible reactions, which usually leads to a smaller cone. Using oriented matroid theory, we show how the projected cone can be computed efficiently and how this can be applied to find iMCSs even in large genome-scale networks.

AVAILABILITY AND IMPLEMENTATION

Software is freely available at https://sourceforge.net/projects/irreversibleminimalcutsets/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Automatic construction of metabolic models with enzyme constraints

Background

In order to improve the accuracy of constraint-based metabolic models, several approaches have been developed which intend to integrate additional biological information. Two of these methods, MOMENT and GECKO, incorporate enzymatic (kcat) parameters and enzyme mass constraints to further constrain the space of feasible metabolic flux distributions. While both methods have been proven to deliver useful extensions of metabolic models, they may considerably increase size and complexity of the models and there is currently no tool available to fully automate generation and calibration of such enzyme-constrained models from given stoichiometric models.

Results

In this work we present three major developments. We first conceived short MOMENT (sMOMENT), a simplified version of the MOMENT approach, which yields the same predictions as MOMENT but requires significantly fewer variables and enables direct inclusion of the relevant enzyme constraints in the standard representation of a constraint-based model. When measurements of enzyme concentrations are available, these can be included as well leading in the extreme case, where all enzyme concentrations are known, to a model representation that is analogous to the GECKO approach. Second, we developed the AutoPACMEN toolbox which allows an almost fully automated creation of sMOMENT-enhanced stoichiometric metabolic models. In particular, this includes the automatic read-out and processing of relevant enzymatic data from different databases and the reconfiguration of the stoichiometric model with embedded enzymatic constraints. Additionally, tools have been developed to adjust (kcat and enzyme pool) parameters of sMOMENT models based on given flux data. We finally applied the new sMOMENT approach and the AutoPACMEN toolbox to generate an enzyme-constrained version of the E. coli genome-scale model iJO1366 and analyze its key properties and differences with the standard model. In particular, we show that the enzyme constraints improve flux predictions (e.g., explaining overflow metabolism and other metabolic switches) and demonstrate, for the first time, that these constraints can markedly change the spectrum of metabolic engineering strategies for different target products.

Conclusions

The methodological and tool developments presented herein pave the way for a simplified and routine construction and analysis of enzyme-constrained metabolic models.

Systems and synthetic biology tools for advanced bioproduction hosts

The genomic revolution ushered in an era of discovery and characterization of enzymes from novel organisms that fueled engineering of microbes to produce commodity and high-value compounds. Over the past decade advances in synthetic biology tools in recent years contributed to significant progress in metabolic engineering efforts to produce both biofuels and bioproducts resulting in several such related items being brought to market. These successes represent a burgeoning bio-economy; however, significant resources and time are still necessary to progress a system from proof-of-concept to market. In order to fully realize this potential, methods that examine biological systems in a comprehensive, systematic and high-throughput manner are essential. Recent success in synthetic biology has coincided with the development of systems biology and analytical approaches that kept pace and scaled with technology development. Here, we review a selection of systems biology methods and their use in synthetic biology approaches for microbial biotechnology platforms.

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.

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.

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.

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 ".

A mathematical framework for yield (vs. rate) optimization in constraint-based modeling and applications in metabolic engineering

BACKGROUND

The optimization of metabolic rates (as linear objective functions) represents the methodical core of flux-balance analysis techniques which have become a standard tool for the study of genome-scale metabolic models. Besides (growth and synthesis) rates, metabolic yields are key parameters for the characterization of biochemical transformation processes, especially in the context of biotechnological applications. However, yields are ratios of rates, and hence the optimization of yields (as nonlinear objective functions) under arbitrary linear constraints is not possible with current flux-balance analysis techniques. Despite the fundamental importance of yields in constraint-based modeling, a comprehensive mathematical framework for yield optimization is still missing.

RESULTS

We present a mathematical theory that allows one to systematically compute and analyze yield-optimal solutions of metabolic models under arbitrary linear constraints. In particular, we formulate yield optimization as a linear-fractional program. For practical computations, we transform the linear-fractional yield optimization problem to a (higher-dimensional) linear problem. Its solutions determine the solutions of the original problem and can be used to predict yield-optimal flux distributions in genome-scale metabolic models. For the theoretical analysis, we consider the linear-fractional problem directly. Most importantly, we show that the yield-optimal solution set (like the rate-optimal solution set) is determined by (yield-optimal) elementary flux vectors of the underlying metabolic model. However, yield- and rate-optimal solutions may differ from each other, and hence optimal (biomass or product) yields are not necessarily obtained at solutions with optimal (growth or synthesis) rates. Moreover, we discuss phase planes/production envelopes and yield spaces, in particular, we prove that yield spaces are convex and provide algorithms for their computation. We illustrate our findings by a small example and demonstrate their relevance for metabolic engineering with realistic models of E. coli.

CONCLUSIONS

We develop a comprehensive mathematical framework for yield optimization in metabolic models. Our theory is particularly useful for the study and rational modification of cell factories designed under given yield and/or rate requirements.

Sequential computation of elementary modes and minimal cut sets in genome-scale metabolic networks using alternate integer linear programming

Motivation

Elementary (flux) modes (EMs) have served as a valuable tool for investigating structural and functional properties of metabolic networks. Identification of the full set of EMs in genome-scale networks remains challenging due to combinatorial explosion of EMs in complex networks. It is often, however, that only a small subset of relevant EMs needs to be known, for which optimization-based sequential computation is a useful alternative. Most of the currently available methods along this line are based on the iterative use of mixed integer linear programming (MILP), the effectiveness of which significantly deteriorates as the number of iterations builds up. To alleviate the computational burden associated with the MILP implementation, we here present a novel optimization algorithm termed alternate integer linear programming (AILP).

Results

Our algorithm was designed to iteratively solve a pair of integer programming (IP) and linear programming (LP) to compute EMs in a sequential manner. In each step, the IP identifies a minimal subset of reactions, the deletion of which disables all previously identified EMs. Thus, a subsequent LP solution subject to this reaction deletion constraint becomes a distinct EM. In cases where no feasible LP solution is available, IP-derived reaction deletion sets represent minimal cut sets (MCSs). Despite the additional computation of MCSs, AILP achieved significant time reduction in computing EMs by orders of magnitude. The proposed AILP algorithm not only offers a computational advantage in the EM analysis of genome-scale networks, but also improves the understanding of the linkage between EMs and MCSs.

Availability and Implementation

The software is implemented in Matlab, and is provided as supplementary information .

Contact

hyunseob.song@pnnl.gov.

Supplementary information

Supplementary data are available at Bioinformatics online.

Designing Optimized Production Hosts by Metabolic Modeling

Many of the complex and expensive production steps in the chemical industry are readily available in living cells. In order to overcome the metabolic limits of these cells, the optimal genetic intervention strategies can be computed by the use of metabolic modeling. Elementary flux mode analysis (EFMA) is an ideal tool for this task, as it does not require defining a cellular objective function. We present two EFMA-based methods to optimize production hosts: (1) the standard approach that can only be used for small and medium scale metabolic networks and (2) the advanced dual system approach that can be utilized to directly compute intervention strategies in a genome-scale metabolic model.

From elementary flux modes to elementary flux vectors: Metabolic pathway analysis with arbitrary linear flux constraints

Elementary flux modes (EFMs) emerged as a formal concept to describe metabolic pathways and have become an established tool for constraint-based modeling and metabolic network analysis. EFMs are characteristic (support-minimal) vectors of the flux cone that contains all feasible steady-state flux vectors of a given metabolic network. EFMs account for (homogeneous) linear constraints arising from reaction irreversibilities and the assumption of steady state; however, other (inhomogeneous) linear constraints, such as minimal and maximal reaction rates frequently used by other constraint-based techniques (such as flux balance analysis [FBA]), cannot be directly integrated. These additional constraints further restrict the space of feasible flux vectors and turn the flux cone into a general flux polyhedron in which the concept of EFMs is not directly applicable anymore. For this reason, there has been a conceptual gap between EFM-based (pathway) analysis methods and linear optimization (FBA) techniques, as they operate on different geometric objects. One approach to overcome these limitations was proposed ten years ago and is based on the concept of elementary flux vectors (EFVs). Only recently has the community started to recognize the potential of EFVs for metabolic network analysis. In fact, EFVs exactly represent the conceptual development required to generalize the idea of EFMs from flux cones to flux polyhedra. This work aims to present a concise theoretical and practical introduction to EFVs that is accessible to a broad audience. We highlight the close relationship between EFMs and EFVs and demonstrate that almost all applications of EFMs (in flux cones) are possible for EFVs (in flux polyhedra) as well. In fact, certain properties can only be studied with EFVs. Thus, we conclude that EFVs provide a powerful and unifying framework for constraint-based modeling of metabolic networks.

Development of a Framework for Metabolic Pathway Analysis-Driven Strain Optimization Methods

Genome-scale metabolic models (GSMMs) have become important assets for rational design of compound overproduction using microbial cell factories. Most computational strain optimization methods (CSOM) using GSMMs, while useful in metabolic engineering, rely on the definition of questionable cell objectives, leading to some bias. Metabolic pathway analysis approaches do not require an objective function. Though their use brings immediate advantages, it has mostly been restricted to small scale models due to computational demands. Additionally, their complex parameterization and lack of intuitive tools pose an important challenge towards making these widely available to the community. Recently, MCSEnumerator has extended the scale of these methods, namely regarding enumeration of minimal cut sets, now able to handle GSMMs. This work proposes a tool implementing this method as a Java library and a plugin within the OptFlux metabolic engineering platform providing a friendly user interface. A standard enumeration problem and pipeline applicable to GSMMs is proposed, making use by the community simpler. To highlight the potential of these approaches, we devised a case study for overproduction of succinate, providing a phenotype analysis of a selected strategy and comparing robustness with a selected solution from a bi-level CSOM.

Optimal knockout strategies in genome-scale metabolic networks using particle swarm optimization

BACKGROUND

Knockout strategies, particularly the concept of constrained minimal cut sets (cMCSs), are an important part of the arsenal of tools used in manipulating metabolic networks. Given a specific design, cMCSs can be calculated even in genome-scale networks. We would however like to find not only the optimal intervention strategy for a given design but the best possible design too. Our solution (PSOMCS) is to use particle swarm optimization (PSO) along with the direct calculation of cMCSs from the stoichiometric matrix to obtain optimal designs satisfying multiple objectives.

RESULTS

To illustrate the working of PSOMCS, we apply it to a toy network. Next we show its superiority by comparing its performance against other comparable methods on a medium sized E. coli core metabolic network. PSOMCS not only finds solutions comparable to previously published results but also it is orders of magnitude faster. Finally, we use PSOMCS to predict knockouts satisfying multiple objectives in a genome-scale metabolic model of E. coli and compare it with OptKnock and RobustKnock.

CONCLUSIONS

PSOMCS finds competitive knockout strategies and designs compared to other current methods and is in some cases significantly faster. It can be used in identifying knockouts which will force optimal desired behaviors in large and genome scale metabolic networks. It will be even more useful as larger metabolic models of industrially relevant organisms become available.

EColiCore2: a reference network model of the central metabolism of Escherichia coli and relationships to its genome-scale parent model

Genome-scale metabolic modeling has become an invaluable tool to analyze properties and capabilities of metabolic networks and has been particularly successful for the model organism Escherichia coli. However, for several applications, smaller metabolic (core) models are needed. Using a recently introduced reduction algorithm and the latest E. coli genome-scale reconstruction iJO1366, we derived EColiCore2, a model of the central metabolism of E. coli. EColiCore2 is a subnetwork of iJO1366 and preserves predefined phenotypes including optimal growth on different substrates. The network comprises 486 metabolites and 499 reactions, is accessible for elementary-modes analysis and can, if required, be further compressed to a network with 82 reactions and 54 metabolites having an identical solution space as EColiCore2. A systematic comparison of EColiCore2 with its genome-scale parent model iJO1366 reveals that several key properties (flux ranges, reaction essentialities, production envelopes) of the central metabolism are preserved in EColiCore2 while it neglects redundancies along biosynthetic routes. We also compare calculated metabolic engineering strategies in both models and demonstrate, as a general result, how intervention strategies found in a core model allow the identification of valid strategies in a genome-scale model. Overall, EColiCore2 holds promise to become a reference model of E. coli's central metabolism.

Exploring complex cellular phenotypes and model-guided strain design with a novel genome-scale metabolic model of Clostridium thermocellum DSM 1313 implementing an adjustable cellulosome

BACKGROUND

Clostridium thermocellum is a gram-positive thermophile that can directly convert lignocellulosic material into biofuels. The metabolism of C. thermocellum contains many branches and redundancies which limit biofuel production, and typical genetic techniques are time-consuming. Further, the genome sequence of a genetically tractable strain C. thermocellum DSM 1313 has been recently sequenced and annotated. Therefore, developing a comprehensive, predictive, genome-scale metabolic model of DSM 1313 is desired for elucidating its complex phenotypes and facilitating model-guided metabolic engineering.

RESULTS

We constructed a genome-scale metabolic model iAT601 for DSM 1313 using the KEGG database as a scaffold and an extensive literature review and bioinformatic analysis for model refinement. Next, we used several sets of experimental data to train the model, e.g., estimation of the ATP requirement for growth-associated maintenance (13.5 mmol ATP/g DCW/h) and cellulosome synthesis (57 mmol ATP/g cellulosome/h). Using our tuned model, we investigated the effect of cellodextrin lengths on cell yields, and could predict in silico experimentally observed differences in cell yield based on which cellodextrin species is assimilated. We further employed our tuned model to analyze the experimentally observed differences in fermentation profiles (i.e., the ethanol to acetate ratio) between cellobiose- and cellulose-grown cultures and infer regulatory mechanisms to explain the phenotypic differences. Finally, we used the model to design over 250 genetic modification strategies with the potential to optimize ethanol production, 6155 for hydrogen production, and 28 for isobutanol production.

CONCLUSIONS

Our developed genome-scale model iAT601 is capable of accurately predicting complex cellular phenotypes under a variety of conditions and serves as a high-quality platform for model-guided strain design and metabolic engineering to produce industrial biofuels and chemicals of interest.

The steady-state assumption in oscillating and growing systems

The steady-state assumption, which states that the production and consumption of metabolites inside the cell are balanced, is one of the key aspects that makes an efficient analysis of genome-scale metabolic networks possible. It can be motivated from two different perspectives. In the time-scales perspective, we use the fact that metabolism is much faster than other cellular processes such as gene expression. Hence, the steady-state assumption is derived as a quasi-steady-state approximation of the metabolism that adapts to the changing cellular conditions. In this article we focus on the second perspective, stating that on the long run no metabolite can accumulate or deplete. In contrast to the first perspective it is not immediately clear how this perspective can be captured mathematically and what assumptions are required to obtain the steady-state condition. By presenting a mathematical framework based on the second perspective we demonstrate that the assumption of steady-state also applies to oscillating and growing systems without requiring quasi-steady-state at any time point. However, we also show that the average concentrations may not be compatible with the average fluxes. In summary, we establish a mathematical foundation for the steady-state assumption for long time periods that justifies its successful use in many applications. Furthermore, this mathematical foundation also pinpoints unintuitive effects in the integration of metabolite concentrations using nonlinear constraints into steady-state models for long time periods.

Which sets of elementary flux modes form thermodynamically feasible flux distributions?

Elementary flux modes (EFMs) are non-decomposable steady-state fluxes through metabolic networks. Every possible flux through a network can be described as a superposition of EFMs. The definition of EFMs is based on the stoichiometry of the network, and it has been shown previously that not all EFMs are thermodynamically feasible. These infeasible EFMs cannot contribute to a biologically meaningful flux distribution. In this work, we show that a set of thermodynamically feasible EFMs need not be thermodynamically consistent. We use first principles of thermodynamics to define the feasibility of a flux distribution and present a method to compute the largest thermodynamically consistent sets (LTCSs) of EFMs. An LTCS contains the maximum number of EFMs that can be combined to form a thermodynamically feasible flux distribution. As a case study we analyze all LTCSs found in Escherichia coli when grown on glucose and show that only one LTCS shows the required phenotypical properties. Using our method, we find that in our E. coli model < 10% of all EFMs are thermodynamically relevant.

Cyanobacterial Biofuels: Strategies and Developments on Network and Modeling

Cyanobacteria, the phototrophic microorganisms, have attracted much attention recently as a promising source for environmentally sustainable biofuels production. However, barriers for commercial markets of cyanobacteria-based biofuels concern the economic feasibility. Miscellaneous strategies for improving the production performance of cyanobacteria have thus been developed. Among these, the simple ad hoc strategies resulting in failure to optimize fully cell growth coupled with desired product yield are explored. With the advancement of genomics and systems biology, a new paradigm toward systems metabolic engineering has been recognized. In particular, a genome-scale metabolic network reconstruction and modeling is a crucial systems-based tool for whole-cell-wide investigation and prediction. In this review, the cyanobacterial genome-scale metabolic models, which offer a system-level understanding of cyanobacterial metabolism, are described. The main process of metabolic network reconstruction and modeling of cyanobacteria are summarized. Strategies and developments on genome-scale network and modeling through the systems metabolic engineering approach are advanced and employed for efficient cyanobacterial-based biofuels production.

In Vivo Validation of In Silico Predicted Metabolic Engineering Strategies in Yeast: Disruption of α-Ketoglutarate Dehydrogenase and Expression of ATP-Citrate Lyase for Terpenoid Production

Background

Engineering of the central carbon metabolism of Saccharomyces cerevisiae to redirect metabolic flux towards cytosolic acetyl-CoA has become a central topic in yeast biotechnology. A cell factory with increased flux into acetyl-CoA can be used for heterologous production of terpenoids for pharmaceuticals, biofuels, fragrances, or other acetyl-CoA derived compounds. In a previous study, we identified promising metabolic engineering targets in S. cerevisiae using an in silico stoichiometric metabolic network analysis. Here, we validate selected in silico strategies in vivo.

Results

Patchoulol was produced by yeast via a heterologous patchoulol synthase of Pogostemon cablin. To increase the metabolic flux from acetyl-CoA towards patchoulol, a truncated HMG-CoA reductase was overexpressed and farnesyl diphosphate synthase was fused with patchoulol synthase. The highest increase in production could be achieved by modifying the carbon source; sesquiterpenoid titer increased from glucose to ethanol by a factor of 8.4. Two strategies predicted in silico were chosen for validation in this work. Disruption of α-ketoglutarate dehydrogenase gene (KGD1) was predicted to redirect the metabolic flux via the pyruvate dehydrogenase bypass towards acetyl-CoA. The metabolic flux was redirected as predicted, however, the effect was dependent on cultivation conditions and the flux was interrupted at the level of acetate. High amounts of acetate were produced. As an alternative pathway to synthesize cytosolic acetyl-CoA, ATP-citrate lyase was expressed as a polycistronic construct, however, in vivo performance of the enzyme needs to be optimized to increase terpenoid production.

Conclusions

Stoichiometric metabolic network analysis can be used successfully as a metabolic prediction tool. However, this study highlights that kinetics, regulation and cultivation conditions may interfere, resulting in poor in vivo performance. Main sites of regulation need to be released and improved enzymes are essential to meet the required activities for an increased product formation in vivo.

Designing minimal microbial strains of desired functionality using a genetic algorithm

Background

The rational, in silico prediction of gene-knockouts to turn organisms into efficient cell factories is an essential and computationally challenging task in metabolic engineering. Elementary flux mode analysis in combination with constraint minimal cut sets is a particularly powerful method to identify optimal engineering targets, which will force an organism into the desired metabolic state. Given an engineering objective, it is theoretically possible, although computationally impractical, to find the best minimal intervention strategies.

Results

We developed a genetic algorithm (GA-MCS) to quickly find many (near) optimal intervention strategies while overcoming the above mentioned computational burden. We tested our algorithm on Escherichia coli metabolic networks of three different sizes to find intervention strategies satisfying three different engineering objectives.

Conclusions

We show that GA-MCS finds all practically relevant targets for any (non)-linear engineering objective. Our algorithm also found solutions comparable to previously published results. We show that for large networks optimal solutions are found within a fraction of the time used for a complete enumeration.

In Silico Constraint-Based Strain Optimization Methods: the Quest for Optimal Cell Factories

Shifting from chemical to biotechnological processes is one of the cornerstones of 21st century industry. The production of a great range of chemicals via biotechnological means is a key challenge on the way toward a bio-based economy. However, this shift is occurring at a pace slower than initially expected. The development of efficient cell factories that allow for competitive production yields is of paramount importance for this leap to happen. Constraint-based models of metabolism, together with in silico strain design algorithms, promise to reveal insights into the best genetic design strategies, a step further toward achieving that goal. In this work, a thorough analysis of the main in silico constraint-based strain design strategies and algorithms is presented, their application in real-world case studies is analyzed, and a path for the future is discussed.

Computing smallest intervention strategies for multiple metabolic networks in a boolean model

This article considers the problem whereby, given two metabolic networks N1 and N2, a set of source compounds, and a set of target compounds, we must find the minimum set of reactions whose removal (knockout) ensures that the target compounds are not producible in N1 but are producible in N2. Similar studies exist for the problem of finding the minimum knockout with the smallest side effect for a single network. However, if technologies of external perturbations are advanced in the near future, it may be important to develop methods of computing the minimum knockout for multiple networks (MKMN). Flux balance analysis (FBA) is efficient if a well-polished model is available. However, that is not always the case. Therefore, in this article, we study MKMN in Boolean models and an elementary mode (EM)-based model. Integer linear programming (ILP)-based methods are developed for these models, since MKMN is NP-complete for both the Boolean model and the EM-based model. Computer experiments are conducted with metabolic networks of clostridium perfringens SM101 and bifidobacterium longum DJO10A, respectively known as bad bacteria and good bacteria for the human intestine. The results show that larger networks are more likely to have MKMN solutions. However, solving for these larger networks takes a very long time, and often the computation cannot be completed. This is reasonable, because small networks do not have many alternative pathways, making it difficult to satisfy the MKMN condition, whereas in large networks the number of candidate solutions explodes. Our developed software minFvskO is available online.

Elucidating central metabolic redox obstacles hindering ethanol production in Clostridium thermocellum

Clostridium thermocellum is an anaerobic, Gram-positive, thermophilic bacterium that has generated great interest due to its ability to ferment lignocellulosic biomass to ethanol. However, ethanol production is low due to the complex and poorly understood branched metabolism of C. thermocellum, and in some cases overflow metabolism as well. In this work, we developed a predictive stoichiometric metabolic model for C. thermocellum which incorporates the current state of understanding, with particular attention to cofactor specificity in the atypical glycolytic enzymes and the complex energy, redox, and fermentative pathways with the goal of aiding metabolic engineering efforts. We validated the model's capability to encompass experimentally observed phenotypes for the parent strain and derived mutants designed for significant perturbation of redox and energy pathways. Metabolic flux distributions revealed significant alterations in key metabolic branch points (e.g., phosphoenol pyruvate, pyruvate, acetyl-CoA, and cofactor nodes) in engineered strains for channeling electron and carbon fluxes for enhanced ethanol synthesis, with the best performing strain doubling ethanol yield and titer compared to the parent strain. In silico predictions of a redox-imbalanced genotype incapable of growth were confirmed in vivo, and a mutant strain was used as a platform to probe redox bottlenecks in the central metabolism that hinder efficient ethanol production. The results highlight the robustness of the redox metabolism of C. thermocellum and the necessity of streamlined electron flux from reduced ferredoxin to NAD(P)H for high ethanol production. The model was further used to design a metabolic engineering strategy to phenotypically constrain C. thermocellum to achieve high ethanol yields while requiring minimal genetic manipulations. The model can be applied to design C. thermocellum as a platform microbe for consolidated bioprocessing to produce ethanol and other reduced metabolites.

Rational design of efficient modular cells

The modular cell design principle is formulated to devise modular (chassis) cells. These cells can be assembled with exchangeable production modules in a plug-and-play fashion to build microbial cell factories for efficient combinatorial biosynthesis of novel molecules, requiring minimal iterative strain optimization steps. A modular cell is designed to be auxotrophic, containing core metabolic pathways that are necessary but insufficient to support cell growth and maintenance. To be functional, it must tightly couple with an exchangeable production module containing auxiliary metabolic pathways that not only complement cell growth but also enhance production of targeted molecules. We developed a MODCELL (modular cell) framework based on metabolic pathway analysis to implement the modular cell design principle. MODCELL identifies genetic modifications and requirements to construct modular cell candidates and their associated exchangeable production modules. By defining the degree of similarity and coupling metrics, MODCELL can evaluate which exchangeable production module(s) can be tightly coupled with a modular cell candidate. We first demonstrated how MODCELL works in a step-by-step manner for example metabolic networks, and then applied it to design modular Escherichia coli cells for efficient combinatorial biosynthesis of five alcohols (ethanol, propanol, isopropanol, butanol and isobutanol) and five butyrate esters (ethyl butyrate, propyl butyrate, isopropyl butyrate, butyl butyrate and isobutyl butyrate) from pentose sugars (arabinose and xylose) and hexose sugars (glucose, mannose, and galactose) under anaerobic conditions. We identified three modular cells, MODCELL1, MODCELL2 and MODCELL3, that can couple well with Group 1 of modules (ethanol, isobutanol, butanol, ethyl butyrate, isobutyl butyrate, butyl butyrate), Group 2 (isopropanol, isopropyl butyrate), and Group 3 (propanol, isopropanol), respectively. We validated the design of MODCELL1 for anaerobic production of ethanol, butanol, and ethyl butyrate using experimental data available in literature.