OptStrain: a computational framework for redesign of microbial production systems

Machine learning for the advancement of genome-scale metabolic modeling

Constraint-based modeling (CBM) has evolved as the core systems biology tool to map the interrelations between genotype, phenotype, and external environment. The recent advancement of high-throughput experimental approaches and multi-omics strategies has generated a plethora of new and precise information from wide-ranging biological domains. On the other hand, the continuously growing field of machine learning (ML) and its specialized branch of deep learning (DL) provide essential computational architectures for decoding complex and heterogeneous biological data. In recent years, both multi-omics and ML have assisted in the escalation of CBM. Condition-specific omics data, such as transcriptomics and proteomics, helped contextualize the model prediction while analyzing a particular phenotypic signature. At the same time, the advanced ML tools have eased the model reconstruction and analysis to increase the accuracy and prediction power. However, the development of these multi-disciplinary methodological frameworks mainly occurs independently, which limits the concatenation of biological knowledge from different domains. Hence, we have reviewed the potential of integrating multi-disciplinary tools and strategies from various fields, such as synthetic biology, CBM, omics, and ML, to explore the biochemical phenomenon beyond the conventional biological dogma. How the integrative knowledge of these intersected domains has improved bioengineering and biomedical applications has also been highlighted. We categorically explained the conventional genome-scale metabolic model (GEM) reconstruction tools and their improvement strategies through ML paradigms. Further, the crucial role of ML and DL in omics data restructuring for GEM development has also been briefly discussed. Finally, the case-study-based assessment of the state-of-the-art method for improving biomedical and metabolic engineering strategies has been elaborated. Therefore, this review demonstrates how integrating experimental and in silico strategies can help map the ever-expanding knowledge of biological systems driven by condition-specific cellular information. This multiview approach will elevate the application of ML-based CBM in the biomedical and bioengineering fields for the betterment of society and the environment.

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.

Development and applications of genome-scale metabolic network models

The genome-scale metabolic network model is an effective tool for characterizing the gene-protein-response relationship in the entire metabolic pathway of an organism. By combining various algorithms, the genome-scale metabolic network model can effectively simulate the influence of a specific environment on the physiological state of cells, optimize the culture conditions of strains, and predict the targets of genetic modification to achieve targeted modification of strains. In this review, we summarize the whole process of model building, sort out the various tools that may be involved in the model building process, and explain the role of various algorithms in model analysis. In addition, we also summarized the application of GSMM in network characteristics, cell phenotypes, metabolic engineering, etc. Finally, we discuss the current challenges facing GSMM.

Flux Balance Analysis of Mammalian Cell Systems

Flux balance analysis (FBA) is a computational methodology to model and analyze the metabolic behavior of cells. In this chapter, we break down the key steps for formulating an FBA model and other FBA-derived methodologies in the context of mammalian cell biology, including strain design, developing cell line-specific models, and conducting flux sampling. We provide annotated COBRApy code for each step to show how it would work in practice.

Simulating compatible solute biosynthesis using a metabolic flux model of the biomining acidophile, Acidithiobacillus ferrooxidans ATCC 23270

Halotolerant, acidophilic, bioleaching microorganisms are crucial to biomining operations that utilize saline water. Compatible solutes play an important role in the adaptation of these microorganisms to saline environments. Acidithiobacillus ferrooxidans ATCC 23270, an iron- and sulfur-oxidizing acidophilic bacterium, synthesizes trehalose as its native compatible solute but is still sensitive to salinity. Recently, halotolerant bioleaching bacteria were found to use ectoine as their key compatible solute. Previously, bioleaching bacteria were recalcitrant to genetic manipulation; however, recent advancements in genetic tools and techniques allow successful genetic modification of A. ferrooxidans ATCC 23270. Therefore, this study aimed to test, in silico, the effect of native and synthetic compatible solute biosynthesis by A. ferrooxidans ATCC 23270 on its growth and metabolism. Metabolic network flux modelling was used to provide a computational framework for the prediction of metabolic fluxes during production of native and synthetic compatible solutes by A. ferrooxidans ATCC 23270, in silico. Complete pathways for trehalose biosynthesis by the bacterium are proposed and captured in the updated metabolic model including a newly discovered UDP-dependent trehalose synthesis pathway. Finally, the effect of nitrogen sources on compatible solute production was simulated and showed that using nitrogen gas as the sole nitrogen source enables the ectoine-producing 'engineered' microbe to oxidize up to 20% more ferrous iron in comparison to the native microbe that only produces trehalose. Therefore, the predictive outcomes of the model have the potential to guide the design and optimization of a halotolerant strain of A. ferrooxidans ATCC 23270 for saline bioleaching operations.

Regulatory-systemic approach in Aspergillus niger for bioleaching improvement by controlling precipitation

In the context of e-waste recycling by fungal bioleaching, nickel and cobalt precipitate as toxic metals by oxalic acid, whereas organic acids, such as citric, act as a high-performance chelating agent in dissolving these metals. Oxalic acid elimination requires an excess and uneconomical carbon source concentration in culture media. To resolve this issue, a novel and straightforward systems metabolic engineering method was devised to switch metabolic flux from oxalic acid to citric acid. In this technique, the genome-scale metabolic model of Aspergillus niger was applied to predicting flux variability and key reactions through the calculation of multiple optimal solutions for cellular regulation. Accordingly, BRENDA regulators and a novel molecular docking-oriented approach were defined a regulatory medium for this end. Then, ligands were evaluated in fungal culture to assess their impact on organic acid production for bioleaching of copper and nickel from waste telecommunication printed circuit boards. The protein structure of oxaloacetate hydrolase was modeled based on homology modeling for molecular docking. Metformin, glutathione, and sodium fluoride were found to be effective as inhibitors of oxalic acid production, enabling the production of 8100 ppm citric acid by controlling cellular metabolism. Indirect bioleaching demonstrated that nickel did not precipitate, and the bioleaching efficiency of copper and nickel increased from 40% and 24% to 61% and 100%, respectively. Bioleaching efficiency was evaluated qualitatively by FE-SEM, EDX, mapping, and XRD analysis. KEY POINTS: • A regulatory-systemic procedure for controlling cellular metabolism was introduced • Metformin inhibited oxalic acid, leading to 8100 ppm citric acid production • Bioleaching of copper and nickel in TPCBs improved by 21% and 76.

Genome-wide identification of overexpression and downregulation gene targets based on the sum of covariances of the outgoing reaction fluxes

In metabolic engineering, predicting gene overexpression targets remains challenging because both endogenous and heterologous genes in a large metabolic space can be candidates, in contrast to gene knockout targets that are confined to endogenous genes. We report the development of iBridge that identifies positive and negative metabolites exerting positive and negative impacts on product formation, respectively, based on the sum of covariances of their outgoing (consuming) reaction fluxes for a target chemical. Then, "bridge" reactions converting negative metabolites to positive metabolites are identified as overexpression targets, while the opposites as downregulation targets. Using iBridge, overexpression and downregulation targets are suggested for the production of 298 chemicals and validated for 36 chemicals experimentally demonstrated in previous studies. Finally, iBridge is employed to engineer Escherichia coli strains capable of producing 10.3 g/L of D-panthenol, a compound not previously produced, as well as putrescine and 4-hydroxyphenyllactate at enhanced titers, 63.7 and 8.3 g/L, respectively.

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.

Integrating Omics Data in Genome-Scale Metabolic Modeling: A Methodological Perspective for Precision Medicine

Recent advancements in omics technologies have generated a wealth of biological data. Integrating these data within mathematical models is essential to fully leverage their potential. Genome-scale metabolic models (GEMs) provide a robust framework for studying complex biological systems. GEMs have significantly contributed to our understanding of human metabolism, including the intrinsic relationship between the gut microbiome and the host metabolism. In this review, we highlight the contributions of GEMs and discuss the critical challenges that must be overcome to ensure their reproducibility and enhance their prediction accuracy, particularly in the context of precision medicine. We also explore the role of machine learning in addressing these challenges within GEMs. The integration of omics data with GEMs has the potential to lead to new insights, and to advance our understanding of molecular mechanisms in human health and disease.

Automating the design-build-test-learn cycle towards next-generation bacterial cell factories

Automation is playing an increasingly significant role in synthetic biology. Groundbreaking technologies, developed over the past 20 years, have enormously accelerated the construction of efficient microbial cell factories. Integrating state-of-the-art tools (e.g. for genome engineering and analytical techniques) into the design-build-test-learn cycle (DBTLc) will shift the metabolic engineering paradigm from an almost artisanal labor towards a fully automated workflow. Here, we provide a perspective on how a fully automated DBTLc could be harnessed to construct the next-generation bacterial cell factories in a fast, high-throughput fashion. Innovative toolsets and approaches that pushed the boundaries in each segment of the cycle are reviewed to this end. We also present the most recent efforts on automation of the DBTLc, which heralds a fully autonomous pipeline for synthetic biology in the near future.

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.

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.

Kinetic-model-based pathway optimization with application to reverse glycolysis in mammalian cells

Over the last two decades, model-based metabolic pathway optimization tools have been developed for the design of microorganisms to produce desired metabolites. However, few have considered more complex cellular systems such as mammalian cells, which requires the use of nonlinear kinetic models to capture the effects of concentration changes and cross-regulatory interactions. In this study, we develop a new two-stage pathway optimization framework based on kinetic models that incorporate detailed kinetics and regulation information. In Stage 1, a set of optimization problems are solved to identify and rank the enzymes that contribute the most to achieving the metabolic objective. Stage 2 then determines the optimal enzyme interventions for specified desired numbers of enzyme adjustments. It also incorporates multi-scenario optimization, which allows the simultaneous consideration of multiple physiological conditions. We apply the proposed framework to find enzyme adjustments that enable a reverse glucose flow in cultured mammalian cells, thereby eliminating the need for glucose feed in the late culture stage and enhancing process robustness. The computational results demonstrate the efficacy of the proposed approach; it not only captures the important regulations and key enzymes for reverse glycolysis but also identifies differences and commonalities in the metabolic requirements for different carbon sources.

Metabolomics and modelling approaches for systems metabolic engineering

Metabolic engineering involves the manipulation of microbes to produce desirable compounds through genetic engineering or synthetic biology approaches. Metabolomics involves the quantitation of intracellular and extracellular metabolites, where mass spectrometry and nuclear magnetic resonance based analytical instrumentation are often used. Here, the experimental designs, sample preparations, metabolite quenching and extraction are essential to the quantitative metabolomics workflow. The resultant metabolomics data can then be used with computational modelling approaches, such as kinetic and constraint-based modelling, to better understand underlying mechanisms and bottlenecks in the synthesis of desired compounds, thereby accelerating research through systems metabolic engineering. Constraint-based models, such as genome scale models, have been used successfully to enhance the yield of desired compounds from engineered microbes, however, unlike kinetic or dynamic models, constraint-based models do not incorporate regulatory effects. Nevertheless, the lack of time-series metabolomic data generation has hindered the usefulness of dynamic models till today. In this review, we show that improvements in automation, dynamic real-time analysis and high throughput workflows can drive the generation of more quality data for dynamic models through time-series metabolomics data generation. Spatial metabolomics also has the potential to be used as a complementary approach to conventional metabolomics, as it provides information on the localization of metabolites. However, more effort must be undertaken to identify metabolites from spatial metabolomics data derived through imaging mass spectrometry, where machine learning approaches could prove useful. On the other hand, single-cell metabolomics has also seen rapid growth, where understanding cell-cell heterogeneity can provide more insights into efficient metabolic engineering of microbes. Moving forward, with potential improvements in automation, dynamic real-time analysis, high throughput workflows, and spatial metabolomics, more data can be produced and studied using machine learning algorithms, in conjunction with dynamic models, to generate qualitative and quantitative predictions to advance metabolic engineering efforts.

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.

Underground metabolism as a rich reservoir for pathway engineering

Motivation

Bioproduction of value-added compounds is frequently achieved by utilizing enzymes from other species. However, expression of such heterologous enzymes can be detrimental due to unexpected interactions within the host cell. Recently, an alternative strategy emerged, which relies on recruiting side activities of host enzymes to establish new biosynthetic pathways. Although such low-level ‘underground’ enzyme activities are prevalent, it remains poorly explored whether they may serve as an important reservoir for pathway engineering.

Results

Here, we use genome-scale modeling to estimate the theoretical potential of underground reactions for engineering novel biosynthetic pathways in Escherichia coli. We found that biochemical reactions contributed by underground enzyme activities often enhance the in silico production of compounds with industrial importance, including several cases where underground activities are indispensable for production. Most of these new capabilities can be achieved by the addition of one or two underground reactions to the native network, suggesting that only a few side activities need to be enhanced during implementation. Remarkably, we find that the contribution of underground reactions to the production of value-added compounds is comparable to that of heterologous reactions, underscoring their biotechnological potential. Taken together, our genome-wide study demonstrates that exploiting underground enzyme activities could be a promising addition to the toolbox of industrial strain development.

Availability and implementation

The data and scripts underlying this article are available on GitHub at https://github.com/pappb/Kovacs-et-al-Underground-metabolism.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Application of the Metabolic Modeling Pipeline in KBase to Categorize Reactions, Predict Essential Genes, and Predict Pathways in an Isolate Genome

The DOE Systems Biology Knowledgebase (KBase) platform offers a range of powerful tools for the reconstruction, refinement, and analysis of genome-scale metabolic models built from microbial isolate genomes. In this chapter, we describe and demonstrate these tools in action with an analysis of isoprene production in the Bacillus subtilis DSM genome. Two different methods are applied to build initial metabolic models for the DSM genome, then the models are gapfilled in three different growth conditions. Next, flux balance analysis (FBA) and flux variability analysis (FVA) techniques are applied to both study the growth of these models in minimal media and classify reactions within each model based on essentiality and functionality. The models are applied with the FBA method to predict essential genes, which are then compared to an updated list of essential genes obtained for B. subtilis 168, a very similar strain to the DSM isolate. The models are also applied to simulate Biolog growth conditions, and these results are compared with Biolog data collected for B. subtilis 168. Finally, the DSM metabolic models are applied to explore the pathways and genes responsible for producing isoprene in this strain. These studies demonstrate the accuracy and utility of models generated from the KBase pipelines, as well as exploring the tools available for analyzing these models.

Escherichia coli Data-Driven Strain Design Using Aggregated Adaptive Laboratory Evolution Mutational Data

Microbes are being engineered for an increasingly large and diverse set of applications. However, the designing of microbial genomes remains challenging due to the general complexity of biological systems. Adaptive Laboratory Evolution (ALE) leverages nature’s problem-solving processes to generate optimized genotypes currently inaccessible to rational methods. The large amount of public ALE data now represents a new opportunity for data-driven strain design. This study describes how novel strain designs, or genome sequences not yet observed in ALE experiments or published designs, can be extracted from aggregated ALE data and demonstrates this by designing, building, and testing three novel Escherichia coli strains with fitnesses comparable to ALE mutants. These designs were achieved through a meta-analysis of aggregated ALE mutations data (63 Escherichia coli K-12 MG1655 based ALE experiments, described by 93 unique environmental conditions, 357 independent evolutions, and 13 957 observed mutations), which additionally revealed global ALE mutation trends that inform on ALE-derived strain design principles. Such informative trends anticipate ALE-derived strain designs as largely gene-centric, as opposed to noncoding, and composed of a relatively small number of beneficial variants (approximately 6). These results demonstrate how strain design efforts can be enhanced by the meta-analysis of aggregated ALE data.

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

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

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.

GeneReg: a constraint-based approach for design of feasible metabolic engineering strategies at the gene level

Motivation

Large-scale metabolic models are widely used to design metabolic engineering strategies for diverse biotechnological applications. However, the existing computational approaches focus on alteration of reaction fluxes and often neglect the manipulations of gene expression to implement these strategies.

Results

Here, we find that the association of genes with multiple reactions leads to infeasibility of engineering strategies at the flux level, since they require contradicting manipulations of gene expression. Moreover, we identify that all of the existing approaches to design gene knockout strategies do not ensure that the resulting design may also require other gene alterations, such as up- or downregulations, to match the desired flux distribution. To address these issues, we propose a constraint-based approach, termed GeneReg, that facilitates the design of feasible metabolic engineering strategies at the gene level and that is readily applicable to large-scale metabolic networks. We show that GeneReg can identify feasible strategies to overproduce ethanol in Escherichia coli and lactate in Saccharomyces cerevisiae, but overproduction of the TCA cycle intermediates is not feasible in five organisms used as cell factories under default growth conditions. Therefore, GeneReg points at the need to couple gene regulation and metabolism to design rational metabolic engineering strategies.

Availability and implementation

https://github.com/MonaRazaghi/GeneReg

Supplementary information

Supplementary data are available at Bioinformatics online.

Stoichiometric Modeling of Artificial String Chemistries Reveals Constraints on Metabolic Network Structure

Uncovering the general principles that govern the structure of metabolic networks is key to understanding the emergence and evolution of living systems. Artificial chemistries can help illuminate this problem by enabling the exploration of chemical reaction universes that are constrained by general mathematical rules. Here, we focus on artificial chemistries in which strings of characters represent simplified molecules, and string concatenation and splitting represent possible chemical reactions. We developed a novel Python package, ARtificial CHemistry NEtwork Toolbox (ARCHNET), to study string chemistries using tools from the field of stoichiometric constraint-based modeling. In addition to exploring the topological characteristics of different string chemistry networks, we developed a network-pruning algorithm that can generate minimal metabolic networks capable of producing a specified set of biomass precursors from a given assortment of environmental nutrients. We found that the composition of these minimal metabolic networks was influenced more strongly by the metabolites in the biomass reaction than the identities of the environmental nutrients. This finding has important implications for the reconstruction of organismal metabolic networks and could help us better understand the rise and evolution of biochemical organization. More generally, our work provides a bridge between artificial chemistries and stoichiometric modeling, which can help address a broad range of open questions, from the spontaneous emergence of an organized metabolism to the structure of microbial communities.

Integrated knowledge mining, genome-scale modeling, and machine learning for predicting Yarrowia lipolytica bioproduction

Predicting bioproduction titers from microbial hosts has been challenging due to complex interactions between microbial regulatory networks, stress responses, and suboptimal cultivation conditions. This study integrated knowledge mining, feature extraction, genome-scale modeling (GSM), and machine learning (ML) to develop a model for predicting Yarrowia lipolytica chemical titers (i.e., organic acids, terpenoids, etc.). First, Y. lipolytica production data, including cultivation conditions, genetic engineering strategies, and product information, was manually collected from literature (~100 papers) and stored as either numerical (e.g., substrate concentrations) or categorical (e.g., bioreactor modes) variables. For each case recorded, central pathway fluxes were estimated using GSMs and flux balance analysis (FBA) to provide metabolic features. Second, a ML ensemble learner was trained to predict strain production titers. Accurate predictions on the test data were obtained for instances with production titers >1 g/L (R² = 0.87). However, the model had reduced predictability for low performance strains (0.01-1 g/L, R² = 0.29) potentially due to biosynthesis bottlenecks not captured in the features. Feature ranking indicated that the FBA fluxes, the number of enzyme steps, the substrate inputs, and thermodynamic barriers (i.e., Gibbs free energy of reaction) were the most influential factors. Third, the model was evaluated on other oleaginous yeasts and indicated there were conserved features for some hosts that can be potentially exploited by transfer learning. The platform was also designed to assist computational strain design tools (such as OptKnock) to screen genetic targets for improved microbial production in light of experimental conditions.

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.

Towards a Systems Biology Approach to Understanding the Lichen Symbiosis: Opportunities and Challenges of Implementing Network Modelling

Lichen associations, a classic model for successful and sustainable interactions between micro-organisms, have been studied for many years. However, there are significant gaps in our understanding about how the lichen symbiosis operates at the molecular level. This review addresses opportunities for expanding current knowledge on signalling and metabolic interplays in the lichen symbiosis using the tools and approaches of systems biology, particularly network modelling. The largely unexplored nature of symbiont recognition and metabolic interdependency in lichens could benefit from applying a holistic approach to understand underlying molecular mechanisms and processes. Together with ‘omics’ approaches, the application of signalling and metabolic network modelling could provide predictive means to gain insights into lichen signalling and metabolic pathways. First, we review the major signalling and recognition modalities in the lichen symbioses studied to date, and then describe how modelling signalling networks could enhance our understanding of symbiont recognition, particularly leveraging omics techniques. Next, we highlight the current state of knowledge on lichen metabolism. We also discuss metabolic network modelling as a tool to simulate flux distribution in lichen metabolic pathways and to analyse the co-dependence between symbionts. This is especially important given the growing number of lichen genomes now available and improved computational tools for reconstructing such models. We highlight the benefits and possible bottlenecks for implementing different types of network models as applied to the study of lichens.

Reaction lumping in metabolic networks for application with thermodynamic metabolic flux analysis

Thermodynamic metabolic flux analysis (TMFA) can narrow down the space of steady-state flux distributions, but requires knowledge of the standard Gibbs free energy for the modelled reactions. The latter are often not available due to unknown Gibbs free energy change of formation [Formula: see text], of metabolites. To optimize the usage of data on thermodynamics in constraining a model, reaction lumping has been proposed to eliminate metabolites with unknown [Formula: see text]. However, the lumping procedure has not been formalized nor implemented for systematic identification of lumped reactions. Here, we propose, implement, and test a combined procedure for reaction lumping, applicable to genome-scale metabolic models. It is based on identification of groups of metabolites with unknown [Formula: see text] whose elimination can be conducted independently of the others via: (1) group implementation, aiming to eliminate an entire such group, and, if this is infeasible, (2) a sequential implementation to ensure that a maximal number of metabolites with unknown [Formula: see text] are eliminated. Our comparative analysis with genome-scale metabolic models of Escherichia coli, Bacillus subtilis, and Homo sapiens shows that the combined procedure provides an efficient means for systematic identification of lumped reactions. We also demonstrate that TMFA applied to models with reactions lumped according to the proposed procedure lead to more precise predictions in comparison to the original models. The provided implementation thus ensures the reproducibility of the findings and their application with standard TMFA.

In silico Design for Systems-Based Metabolic Engineering for the Bioconversion of Valuable Compounds From Industrial By-Products

Selecting appropriate metabolic engineering targets to build efficient cell factories maximizing the bioconversion of industrial by-products to valuable compounds taking into account time restrictions is a significant challenge in industrial biotechnology. Microbial metabolism engineering following a rational design has been widely studied. However, it is a cost-, time-, and laborious-intensive process because of the cell network complexity; thus, it is important to use tools that allow predicting gene deletions. An in silico experiment was performed to model and understand the metabolic engineering effects on the cell factory considering a second complexity level by transcriptomics data integration. In this study, a systems-based metabolic engineering target prediction was used to increase glycerol bioconversion to succinic acid based on Escherichia coli. Transcriptomics analysis suggests insights on how to increase cell glycerol utilization to further design efficient cell factories. Three E. coli models were used: a core model, a second model based on the integration of transcriptomics data obtained from growth in an optimized culture media, and a third one obtained after integration of transcriptomics data from adaptive laboratory evolution (ALE) experiments. A total of 2,402 strains were obtained with fumarase and pyruvate dehydrogenase being frequently predicted for all the models, suggesting these reactions as essential to increase succinic acid production. Finally, based on using flux balance analysis (FBA) results for all the mutants predicted, a machine learning method was developed to predict new mutants as well as to propose optimal metabolic engineering targets and mutants based on the measurement of the importance of each knockout's (feature's) contribution. Glycerol has become an interesting carbon source for industrial processes due to biodiesel business growth since it has shown promising results in terms of biomass/substrate yields. The combination of transcriptome, systems metabolic modeling, and machine learning analyses revealed the versatility of computational models to predict key metabolic engineering targets in a less cost-, time-, and laborious-intensive process. These data provide a platform to improve the prediction of metabolic engineering targets to design efficient cell factories. Our results may also work as a guide and platform for the selection/engineering of microorganisms for the production of interesting chemical compounds.

Genome-scale metabolic modeling of P. thermoglucosidasius NCIMB 11955 reveals metabolic bottlenecks in anaerobic metabolism

Parageobacillus thermoglucosidasius represents a thermophilic, facultative anaerobic bacterial chassis, with several desirable traits for metabolic engineering and industrial production. To further optimize strain productivity, a systems level understanding of its metabolism is needed, which can be facilitated by a genome-scale metabolic model. Here, we present p-thermo, the most complete, curated and validated genome-scale model (to date) of Parageobacillus thermoglucosidasius NCIMB 11955. It spans a total of 890 metabolites, 1175 reactions and 917 metabolic genes, forming an extensive knowledge base for P. thermoglucosidasius NCIMB 11955 metabolism. The model accurately predicts aerobic utilization of 22 carbon sources, and the predictive quality of internal fluxes was validated with previously published ¹³C-fluxomics data. In an application case, p-thermo was used to facilitate more in-depth analysis of reported metabolic engineering efforts, giving additional insight into fermentative metabolism. Finally, p-thermo was used to resolve a previously uncharacterised bottleneck in anaerobic metabolism, by identifying the minimal required supplemented nutrients (thiamin, biotin and iron(III)) needed to sustain anaerobic growth. This highlights the usefulness of p-thermo for guiding the generation of experimental hypotheses and for facilitating data-driven metabolic engineering, expanding the use of P. thermoglucosidasius as a high yield production platform.

Finding branched pathways in metabolic network via atom group tracking

Finding non-standard or new metabolic pathways has important applications in metabolic engineering, synthetic biology and the analysis and reconstruction of metabolic networks. Branched metabolic pathways dominate in metabolic networks and depict a more comprehensive picture of metabolism compared to linear pathways. Although progress has been developed to find branched metabolic pathways, few efforts have been made in identifying branched metabolic pathways via atom group tracking. In this paper, we present a pathfinding method called BPFinder for finding branched metabolic pathways by atom group tracking, which aims to guide the synthetic design of metabolic pathways. BPFinder enumerates linear metabolic pathways by tracking the movements of atom groups in metabolic network and merges the linear atom group conserving pathways into branched pathways. Two merging rules based on the structure of conserved atom groups are proposed to accurately merge the branched compounds of linear pathways to identify branched pathways. Furthermore, the integrated information of compound similarity, thermodynamic feasibility and conserved atom groups is also used to rank the pathfinding results for feasible branched pathways. Experimental results show that BPFinder is more capable of recovering known branched metabolic pathways as compared to other existing methods, and is able to return biologically relevant branched pathways and discover alternative branched pathways of biochemical interest. The online server of BPFinder is available at http://114.215.129.245:8080/atomic/. The program, source code and data can be downloaded from https://github.com/hyr0771/BPFinder.

Computational search of branched metabolic pathways is a fundamental problem in metabolic engineering and metabolic network analysis, which provides a systematic way of understanding the metabolism and discovering alternative pathways for synthesis of useful biomolecules. We propose BPFinder, a novel computational approach to identify branched metabolic pathways via atom group tracking. Different from other pathfinding methods using atom tracking, BPFinder tracks the movement of atom groups in metabolic network to find linear atom group conserving pathways, and merge the found linear pathways by the selected branched compounds to generate branched pathways. Based on the structure of conserved atom groups in branched compounds, we design two merging rules for branched compounds: overlapping rule and non-overlapping rule. The user can flexibly adopt these rules to accurately find the branched pathways that contain overlapping/non-overlapping conserved atom groups. BPFinder also enables the user to combine the information of compound similarity, Gibbs free energy of reactions, and conserved atom groups to sort resulting pathways. Compared with other existing methods, BPFinder can more accurately recover the known branched pathways. The alternative branched pathways returned by BPFinder reveal that the user can flexibly utilize our proposed merging rules to discover biochemically meaningful pathways of interest.

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.

NetFlow: A tool for isolating carbon flows in genome-scale metabolic networks

Genome-scale stoichiometric models (GSMs) have been widely utilized to predict and understand cellular metabolism. GSMs and the flux predictions resulting from them have proven indispensable to fields ranging from metabolic engineering to human disease. Nonetheless, it is challenging to parse these flux predictions due to the inherent size and complexity of the GSMs. Several previous approaches have reduced this complexity by identifying key pathways contained within the genome-scale flux predictions. However, a reduction method that overlays carbon atom transitions on stoichiometry and flux predictions is lacking. To fill this gap, we developed NetFlow, an algorithm that leverages genome-scale carbon mapping to extract and quantitatively distinguish biologically relevant metabolic pathways from a given genome-scale flux prediction. NetFlow extends prior approaches by utilizing both full carbon mapping and context-specific flux predictions. Thus, NetFlow is uniquely able to quantitatively distinguish between biologically relevant pathways of carbon flow within the given flux map. NetFlow simulates ¹³C isotope labeling experiments to calculate the extent of carbon exchange, or carbon yield, between every metabolite in the given GSM. Based on the carbon yield, the carbon flow to or from any metabolite or between any pair of metabolites of interest can be isolated and readily visualized. The resulting pathways are much easier to interpret, which enables an in-depth mechanistic understanding of the metabolic phenotype of interest. Here, we first demonstrate NetFlow with a simple network. We then depict the utility of NetFlow on a model of central carbon metabolism in E. coli. Specifically, we isolated the production pathway for succinate synthesis in this model and the metabolic mechanism driving the predicted increase in succinate yield in a double knockout of E. coli. Finally, we describe the application of NetFlow to a GSM of lycopene-producing E. coli, which enabled the rapid identification of the mechanisms behind the measured increases in lycopene production following single, double, and triple knockouts.

PATHcre8: A Tool That Facilitates the Searching for Heterologous Biosynthetic Routes

Developing computational tools that can facilitate the rational design of cell factories producing desired products at increased yields is challenging, as the tool needs to take into account that the preferred host organism usually has compounds that are consumed by competing reactions that reduce the yield of the desired product. On the other hand, the preferred host organisms may not have the native metabolic reactions needed to produce the compound of interest; thus, the computational tool needs to identify the metabolic reactions that will most efficiently produce the desired product. In this regard, we developed the generic tool PATH^cre8 to facilitate an optimized search for heterologous biosynthetic pathway routes. PATH^cre8 finds and ranks biosynthesis routes in a large number of organisms, including Cyanobacteria. The tool ranks the pathways based on feature scores that reflect reaction thermodynamics, the potentially toxic products in the pathway (compound toxicity), intermediate products in the pathway consumed by competing reactions (product consumption), and host-specific information such as enzyme copy number. A comparison with several other similar tools shows that PATH^cre8 is more efficient in ranking functional pathways. To illustrate the effectiveness of PATH^cre8, we further provide case studies focused on isoprene production and the biodegradation of cocaine. PATH^cre8 is free for academic and nonprofit users and can be accessed at https://www.cbrc.kaust.edu.sa/pathcre8/.

Metabolic modelling approaches for describing and engineering microbial communities

Microbes do not live in isolation but in microbial communities. The relevance of microbial communities is increasing due to growing awareness of their influence on a huge number of environmental, health and industrial processes. Hence, being able to control and engineer the output of both natural and synthetic communities would be of great interest. However, most of the available methods and biotechnological applications involving microorganisms, both in vivo and in silico, have been developed in the context of isolated microbes. In vivo microbial consortia development is extremely difficult and costly because it implies replicating suitable environments in the wet-lab. Computational approaches are thus a good, cost-effective alternative to study microbial communities, mainly via descriptive modelling, but also via engineering modelling. In this review we provide a detailed compilation of examples of engineered microbial communities and a comprehensive, historical revision of available computational metabolic modelling methods to better understand, and rationally engineer wild and synthetic microbial communities.

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.

An insight to flux-balance analysis for biochemical networks

Systems biology is one of the integrated ways to study biological systems and is more favourable than the earlier used approaches. It includes metabolic pathway analysis, modelling, and regulatory as well as signal transduction for getting insights into cellular behaviour. Among the various techniques of modelling, simulation, analysis of networks and pathways, flux-based analysis (FBA) has been recognised because of its extensibility as well as simplicity. It is widely accepted because it is not like a mechanistic simulation which depends on accurate kinetic data. The study of fluxes through the network is informative and can give insights even in the absence of kinetic data. FBA is one of the widely used tools to study biochemical networks and needs information of reaction stoichiometry, growth requirements, specific measurement parameters of the biological system, in particular the reconstruction of the metabolic network for the genome-scale, many of which have already been built previously. It defines the boundaries of flux distributions which are possible and achievable with a defined set of genes. This review article gives an insight into FBA, from the extension of flux balancing to mathematical representation followed by a discussion about the formulation of flux-balance analysis problems, defining constraints for the stoichiometry of the pathways and the tools that can be used in FBA such as FASIMA, COBRA toolbox, and OptFlux. It also includes broader areas in terms of applications which can be covered by FBA as well as the queries which can be addressed through FBA.

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.

Genome-scale metabolic reconstruction of the non-model yeast Issatchenkia orientalis SD108 and its application to organic acids production

Many platform chemicals can be produced from renewable biomass by microorganisms, with organic acids making up a large fraction. Intolerance to the resulting low pH growth conditions, however, remains a challenge for the industrial production of organic acids by microorganisms. Issatchenkia orientalis SD108 is a promising host for industrial production because it is tolerant to acidic conditions as low as pH 2.0. With the goal to systematically assess the metabolic capabilities of this non-model yeast, we developed a genome-scale metabolic model for I. orientalis SD108 spanning 850 genes, 1826 reactions, and 1702 metabolites. In order to improve the model’s quantitative predictions, organism-specific macromolecular composition and ATP maintenance requirements were determined experimentally and implemented. We examined its network topology, including essential genes and flux coupling analysis and drew comparisons with the Yeast 8.3 model for Saccharomyces cerevisiae. We explored the carbon substrate utilization and examined the organism’s production potential for the industrially-relevant succinic acid, making use of the OptKnock framework to identify gene knockouts which couple production of the targeted chemical to biomass production. The genome-scale metabolic model iIsor850 is a data-supported curated model which can inform genetic interventions for overproduction.

•

Genome-scale metabolic model iIsor850 describes metabolism of I. orientalis SD108.

•

Customized biomass reaction highlights differences with S. cerevisiae.

•

Chemostat data elucidate growth-associated ATP maintenance.

•

Substrate utilization and CRISPR/Cas9 gene knockout phenotypes validate model.

•

Model pinpoints candidate gene deletions coupling succinic acid production to growth.

Metabolic pathway engineering: Perspectives and applications

BACKGROUND

Metabolic engineering aims at contriving microbes as biocatalysts for enhanced and cost-effective production of countless secondary metabolites. These secondary metabolites can be treated as the resources of industrial chemicals, pharmaceuticals and fuels. Plants are also crucial targets for metabolic engineers to produce necessary secondary metabolites. Metabolic engineering of both microorganism and plants also contributes towards drug discovery. In order to implement advanced metabolic engineering techniques efficiently, metabolic engineers should have detailed knowledge about cell physiology and metabolism. Principle behind methodologies: Genome-scale mathematical models of integrated metabolic, signal transduction, gene regulatory and protein-protein interaction networks along with experimental validation can provide such knowledge in this context. Incorporation of omics data into these models is crucial in the case of drug discovery. Inverse metabolic engineering and metabolic control analysis (MCA) can help in developing such models. Artificial intelligence methodology can also be applied for efficient and accurate metabolic engineering.

CONCLUSION

In this review, we discuss, at the beginning, the perspectives of metabolic engineering and its application on microorganism and plant leading to drug discovery. At the end, we elaborate why inverse metabolic engineering and MCA are closely related to modern metabolic engineering. In addition, some crucial steps ensuring efficient and optimal metabolic engineering strategies have been discussed. Moreover, we explore the use of genomics data for the activation of silent metabolic clusters and how it can be integrated with metabolic engineering. Finally, we exhibit a few applications of artificial intelligence to metabolic engineering.

Harnessing Natural Modularity of Metabolism with Goal Attainment Optimization to Design a Modular Chassis Cell for Production of Diverse Chemicals

Modular design is key to achieve efficient and robust systems across engineering disciplines. Modular design potentially offers advantages to engineer microbial systems for biocatalysis, bioremediation, and biosensing, overcoming the slow and costly design-build-test-learn cycles in the conventional cell engineering approach. These systems consist of a modular (chassis) cell compatible with exchangeable modules that enable programmed functions such as overproduction of a desirable chemical. We previously proposed a multiobjective optimization framework coupled with metabolic flux models to design modular cells and solved it using multiobjective evolutionary algorithms. However, such approach might not achieve solution optimality and hence limits design options for experimental implementation. In this study, we developed the goal attainment formulation compatible with optimization algorithms that guarantee solution optimality. We applied goal attainment to design an Escherichia coli modular cell capable of synthesizing all molecules in a biochemically diverse library at high yields and rates with only a few genetic manipulations. To elucidate modular organization of the designed cells, we developed a flux variance clustering (FVC) method by identifying reactions with high flux variance and clustering them to identify metabolic modules. Using FVC, we identified reaction usage patterns for different modules in the modular cell, revealing that its broad pathway compatibility is enabled by the natural modularity and flexible flux capacity of endogenous core metabolism. Overall, this study not only sheds light on modularity in metabolic networks from their topology and metabolic functions but also presents a useful synthetic biology toolbox to design modular cells with broad applications.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

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.

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.

Biochemical Characteristics and a Genome-Scale Metabolic Model of an Indian Euryhaline Cyanobacterium with High Polyglucan Content

Marine cyanobacteria are promising microbes to capture and convert atmospheric CO₂ and light into biomass and valuable industrial bio-products. Yet, reports on metabolic characteristics of non-model cyanobacteria are scarce. In this report, we show that an Indian euryhaline Synechococcus sp. BDU 130192 has biomass accumulation comparable to a model marine cyanobacterium and contains approximately double the amount of total carbohydrates, but significantly lower protein levels compared to Synechococcus sp. PCC 7002 cells. Based on its annotated chromosomal genome sequence, we present a genome scale metabolic model (GSMM) of this cyanobacterium, which we have named as iSyn706. The model includes 706 genes, 908 reactions, and 900 metabolites. The difference in the flux balance analysis (FBA) predicted flux distributions between Synechococcus sp. PCC 7002 and Synechococcus sp. BDU130192 strains mimicked the differences in their biomass compositions. Model-predicted oxygen evolution rate for Synechococcus sp. BDU130192 was found to be close to the experimentally-measured value. The model was analyzed to determine the potential of the strain for the production of various industrially-useful products without affecting growth significantly. This model will be helpful to researchers interested in understanding the metabolism as well as to design metabolic engineering strategies for the production of industrially-relevant compounds.

Genome-Scale Metabolic Modeling with Protein Expressions of Normal and Cancerous Colorectal Tissues for Oncogene Inference

Although cancer has historically been regarded as a cell proliferation disorder, it has recently been considered a metabolic disease. The first discovery of metabolic alterations in cancer cells refers to Otto Warburg’s observations. Cancer metabolism results in alterations in metabolic fluxes that are evident in cancer cells compared with most normal tissue cells. This study applied protein expressions of normal and cancer cells to reconstruct two tissue-specific genome-scale metabolic models. Both models were employed in a tri-level optimization framework to infer oncogenes. Moreover, this study also introduced enzyme pseudo-coding numbers in the gene association expression to avoid performing posterior decision-making that is necessary for the reaction-based method. Colorectal cancer (CRC) was the topic of this case study, and 20 top-ranked oncogenes were determined. Notably, these dysregulated genes were involved in various metabolic subsystems and compartments. We found that the average similarity ratio for each dysregulation is higher than 98%, and the extent of similarity for flux changes is higher than 93%. On the basis of surveys of PubMed and GeneCards, these oncogenes were also investigated in various carcinomas and diseases. Most dysregulated genes connect to catalase that acts as a hub and connects protein signaling pathways, such as those involving TP53, mTOR, AKT1, MAPK1, EGFR, MYC, CDK8, and RAS family.