Application of Genome-Scale Metabolic Models in Metabolic Engineering

Genome-Scale Metabolic Models in Fungal Pathogens: Past, Present, and Future

Fungi are diverse organisms with various characteristics and functions. Some play a role in recycling essential elements, such as nitrogen and carbon, while others are utilized in the food and drink production industry. Some others are known to cause diseases in various organisms, including humans. Fungal pathogens cause superficial, subcutaneous, and systemic infections. Consequently, many scientists have focused on studying the factors contributing to the development of human diseases. Therefore, multiple approaches have been assessed to examine the biology of these intriguing organisms. The genome-scale metabolic models (GEMs) have demonstrated many advantages to microbial metabolism studies and the ability to propose novel therapeutic alternatives. Despite significant advancements, much remains to be elucidated regarding the use of this tool for investigating fungal metabolism. This review aims to compile the data provided by the published GEMs of human fungal pathogens. It gives specific examples of the most significant contributions made by these models, examines the advantages and difficulties associated with using such models, and explores the novel approaches suggested to enhance and refine their development.

Systematic analysis of microorganisms' metabolism for selective targeting

Selective drugs with a relatively narrow spectrum can reduce the side effects of treatments compared to broad-spectrum antibiotics by specifically targeting the pathogens responsible for infection. Furthermore, combating an infectious pathogen, especially a drug-resistant microorganism, is more efficient by attacking multiple targets. Here, we combined synthetic lethality with selective drug targeting to identify multi-target and organism-specific potential drug candidates by systematically analyzing the genome-scale metabolic models of six different microorganisms. By considering microorganisms as targeted or conserved in groups ranging from one to six members, we designed 665 individual case studies. For each case, we identified single essential reactions as well as double, triple, and quadruple synthetic lethal reaction sets that are lethal for targeted microorganisms and neutral for conserved ones. As expected, the number of obtained solutions for each case depends on the genomic similarity between the studied microorganisms. Mapping the identified potential drug targets to their corresponding pathways highlighted the importance of key subsystems such as cell envelope biosynthesis, glycerophospholipid metabolism, membrane lipid metabolism, and the nucleotide salvage pathway. To assist in the validation and further investigation of our proposed potential drug targets, we introduced two sets of targets that can theoretically address a substantial portion of the 665 cases. We expect that the obtained solutions provide valuable insights into designing narrow-spectrum drugs that selectively cause system-wide damage only to the target microorganisms.

A genome-scale metabolic model of the effect of dissolved oxygen on 1,3-propanediol fermentation by Klebsiella pneumoniae

Although 1,3-propanediol (1,3-PD) is usually considered an anaerobic fermentation product from glycerol by Klebsiella pneumoniae, microaerobic conditions proved to be more conducive to 1,3-PD production. In this study, a genome-scale metabolic model (GSMM) specific to K. pneumoniae KG2, a high 1.3-PD producer, was constructed. The iZY1242 model contains 2090 reactions, 1242 genes and 1433 metabolites. The model was not only able to accurately characterise cell growth, but also accurately simulate the fed-batch 1,3-PD fermentation process. Flux balance analyses by iZY1242 was performed to dissect the mechanism of stimulated 1,3-PD production under microaerobic conditions, and the maximum yield of 1,3-PD on glycerol was 0.83 mol/mol under optimal microaerobic conditions. Combined with experimental data, the iZY1242 model is a useful tool for establishing the best conditions for microaeration fermentation to produce 1,3-PD from glycerol in K. pneumoniae.

Next-Generation Genome-Scale Metabolic Modeling through Integration of Regulatory Mechanisms

Genome-scale metabolic models (GEMs) are powerful tools for understanding metabolism from a systems-level perspective. However, GEMs in their most basic form fail to account for cellular regulation. A diverse set of mechanisms regulate cellular metabolism, enabling organisms to respond to a wide range of conditions. This limitation of GEMs has prompted the development of new methods to integrate regulatory mechanisms, thereby enhancing the predictive capabilities and broadening the scope of GEMs. Here, we cover integrative models encompassing six types of regulatory mechanisms: transcriptional regulatory networks (TRNs), post-translational modifications (PTMs), epigenetics, protein-protein interactions and protein stability (PPIs/PS), allostery, and signaling networks. We discuss 22 integrative GEM modeling methods and how these have been used to simulate metabolic regulation during normal and pathological conditions. While these advances have been remarkable, there remains a need for comprehensive and widespread integration of regulatory constraints into GEMs. We conclude by discussing challenges in constructing GEMs with regulation and highlight areas that need to be addressed for the successful modeling of metabolic regulation. Next-generation integrative GEMs that incorporate multiple regulatory mechanisms and their crosstalk will be invaluable for discovering cell-type and disease-specific metabolic control mechanisms.

The Plasma Membrane at the Cornerstone Between Flexibility and Adaptability: Implications for Saccharomyces cerevisiae as a Cell Factory

In the last decade, microbial-based biotechnological processes are paving the way toward sustainability as they implemented the use of renewable feedstocks. Nonetheless, the viability and competitiveness of these processes are often limited due to harsh conditions such as: the presence of feedstock-derived inhibitors including weak acids, non-uniform nature of the substrates, osmotic pressure, high temperature, extreme pH. These factors are detrimental for microbial cell factories as a whole, but more specifically the impact on the cell’s membrane is often overlooked. The plasma membrane is a complex system involved in major biological processes, including establishing and maintaining transmembrane gradients, controlling uptake and secretion, intercellular and intracellular communication, cell to cell recognition and cell’s physical protection. Therefore, when designing strategies for the development of versatile, robust and efficient cell factories ready to tackle the harshness of industrial processes while delivering high values of yield, titer and productivity, the plasma membrane has to be considered. Plasma membrane composition comprises diverse macromolecules and it is not constant, as cells adapt it according to the surrounding environment. Remarkably, membrane-specific traits are emerging properties of the system and therefore it is not trivial to predict which membrane composition is advantageous under certain conditions. This review includes an overview of membrane engineering strategies applied to Saccharomyces cerevisiae to enhance its fitness under industrially relevant conditions as well as strategies to increase microbial production of the metabolites of interest.

Machine learning applied for metabolic flux-based control of micro-aerated fermentations in bioreactors

Various bio-based processes depend on controlled micro-aerobic conditions to achieve a satisfactory product yield. However, the limiting oxygen concentration varies according to the micro-organism employed, while for industrial applications, there is no cost-effective way of measuring it at low levels. This study proposes a machine learning procedure within a metabolic flux-based control strategy (SUPERSYS_MCU) to address this issue. The control strategy used simulations of a genome-scale metabolic model to generate a surrogate model in the form of an artificial neural network, to be used in a micro-aerobic fermentation strategy (MF-ANN). The meta-model provided setpoints to the controller, allowing adjustment of the inlet air flow to control the oxygen uptake rate. The strategy was evaluated in micro-aerobic batch cultures employing industrial Saccharomyces cerevisiae yeast, with defined medium and glucose as the carbon source, as a case study. The performance of the proposed control scheme was compared with a conventional fermentation and with three previously reported micro-aeration strategies, including respiratory quotient-based control and constant air flow rate. Due to maintenance of the oxidative balance at the anaerobiosis threshold, the MF-ANN provided volumetric ethanol productivity of 4.16 g·L^-1 ·h^-1 and a yield of 0.48 g_ethanol .g_substrate ^-1 , which were higher than the values achieved for the other conditions studied (maximum of 3.4 g·L^-1 ·h^-1 and 0.35-0.40 g_ethanol ·g_substrate ^-1 , respectively). Due to its modular character, the MF-ANN strategy could be adapted to other micro-aerated bioprocesses.

A proteome-integrated, carbon source dependent genetic regulatory network in Saccharomyces cerevisiae

Integrated regulatory networks can be powerful tools to examine and test properties of cellular systems, such as modelling environmental effects on the molecular bioeconomy, where protein levels are altered in response to changes in growth conditions. Although extensive regulatory pathways and protein interaction data sets exist which represent such networks, few have formally considered quantitative proteomics data to validate and extend them. We generate and consider such data here using a label-free proteomics strategy to quantify alterations in protein abundance for S. cerevisiae when grown on minimal media using glucose, galactose, maltose and trehalose as sole carbon sources. Using a high quality-controlled subset of proteins observed to be differentially abundant, we constructed a proteome-informed network, comprising 1850 transcription factor interactions and 37 chaperone interactions, which defines the major changes in the cellular proteome when growing under different carbon sources. Analysis of the differentially abundant proteins involved in the regulatory network pointed to their significant roles in specific metabolic pathways and function, including glucose homeostasis, amino acid biosynthesis, and carbohydrate metabolic process. We noted strong statistical enrichment in the differentially abundant proteome of targets of known transcription factors associated with stress responses and altered carbon metabolism. This shows how such integrated analysis can lend further experimental support to annotated regulatory interactions, since the proteomic changes capture both magnitude and direction of gene expression change at the level of the affected proteins. Overall this study highlights the power of quantitative proteomics to help define regulatory systems pertinent to environmental conditions.

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

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

Discovery and Evaluation of Biosynthetic Pathways for the Production of Five Methyl Ethyl Ketone Precursors

The limited supply of fossil fuels and the establishment of new environmental policies shifted research in industry and academia toward sustainable production of the second generation of biofuels, with methyl ethyl ketone (MEK) being one promising fuel candidate. MEK is a commercially valuable petrochemical with an extensive application as a solvent. However, as of today, a sustainable and economically viable production of MEK has not yet been achieved despite several attempts of introducing biosynthetic pathways in industrial microorganisms. We used BNICE.ch as a retrobiosynthesis tool to discover all novel pathways around MEK. Out of 1325 identified compounds connecting to MEK with one reaction step, we selected 3-oxopentanoate, but-3-en-2-one, but-1-en-2-olate, butylamine, and 2-hydroxy-2-methylbutanenitrile for further study. We reconstructed 3 679 610 novel biosynthetic pathways toward these 5 compounds. We then embedded these pathways into the genome-scale model of E. coli, and a set of 18 622 were found to be the most biologically feasible ones on the basis of thermodynamics and their yields. For each novel reaction in the viable pathways, we proposed the most similar KEGG reactions, with their gene and protein sequences, as candidates for either a direct experimental implementation or as a basis for enzyme engineering. Through pathway similarity analysis we classified the pathways and identified the enzymes and precursors that were indispensable for the production of the target molecules. These retrobiosynthesis studies demonstrate the potential of BNICE.ch for discovery, systematic evaluation, and analysis of novel pathways in synthetic biology and metabolic engineering studies.

Methods for automated genome-scale metabolic model reconstruction

In the era of next-generation sequencing and ubiquitous assembly and binning of metagenomes, new putative genome sequences are being produced from isolate and microbiome samples at ever-increasing rates. Genome-scale metabolic models have enormous utility for supporting the analysis and predictive characterization of these genomes based on sequence data. As a result, tools for rapid automated reconstruction of metabolic models are becoming critically important for supporting the analysis of new genome sequences. Many tools and algorithms have now emerged to support rapid model reconstruction and analysis. Here, we are comparing and contrasting the capabilities and output of a variety of these tools, including ModelSEED, Raven Toolbox, PathwayTools, SuBliMinal Toolbox and merlin.

Systems Metabolic Engineering of Escherichia coli

Systems metabolic engineering, which recently emerged as metabolic engineering integrated with systems biology, synthetic biology, and evolutionary engineering, allows engineering of microorganisms on a systemic level for the production of valuable chemicals far beyond its native capabilities. Here, we review the strategies for systems metabolic engineering and particularly its applications in Escherichia coli. First, we cover the various tools developed for genetic manipulation in E. coli to increase the production titers of desired chemicals. Next, we detail the strategies for systems metabolic engineering in E. coli, covering the engineering of the native metabolism, the expansion of metabolism with synthetic pathways, and the process engineering aspects undertaken to achieve higher production titers of desired chemicals. Finally, we examine a couple of notable products as case studies produced in E. coli strains developed by systems metabolic engineering. The large portfolio of chemical products successfully produced by engineered E. coli listed here demonstrates the sheer capacity of what can be envisioned and achieved with respect to microbial production of chemicals. Systems metabolic engineering is no longer in its infancy; it is now widely employed and is also positioned to further embrace next-generation interdisciplinary principles and innovation for its upgrade. Systems metabolic engineering will play increasingly important roles in developing industrial strains including E. coli that are capable of efficiently producing natural and nonnatural chemicals and materials from renewable nonfood biomass.

Integrative analysis of human omics data using biomolecular networks

High-throughput '-omics' technologies have given rise to an increasing abundance of genome-scale data detailing human biology at the molecular level. Although these datasets have already made substantial contributions to a more comprehensive understanding of human physiology and diseases, their interpretation becomes increasingly cryptic and nontrivial as they continue to expand in size and complexity. Systems biology networks offer a scaffold upon which omics data can be integrated, facilitating the extraction of new and physiologically relevant information from the data. Two of the most prevalent networks that have been used for such integrative analyses of omics data are genome-scale metabolic models (GEMs) and protein-protein interaction (PPI) networks, both of which have demonstrated success among many different omics and sample types. This integrative approach seeks to unite 'top-down' omics data with 'bottom-up' biological networks in a synergistic fashion that draws on the strengths of both strategies. As the volume and resolution of high-throughput omics data continue to grow, integrative network-based analyses are expected to play an increasingly important role in their interpretation.

Improving the flux distributions simulated with genome-scale metabolic models of Saccharomyces cerevisiae

Genome-scale metabolic models (GEMs) can be used to evaluate genotype-phenotype relationships and their application to microbial strain engineering is increasing in popularity. Some of the algorithms used to simulate the phenotypes of mutant strains require the determination of a wild-type flux distribution. However, the accuracy of this reference, when calculated with flux balance analysis, has not been studied in detail before.

Here, the wild-type simulations of selected GEMs for Saccharomyces cerevisiae have been analysed and most of the models tested predicted erroneous fluxes in central pathways, especially in the pentose phosphate pathway. Since the problematic fluxes were mostly related to areas of the metabolism consuming or producing NADPH/NADH, we have manually curated all reactions including these cofactors by forcing the use of NADPH/NADP⁺ in anabolic reactions and NADH/NAD⁺ for catabolic reactions. The curated models predicted more accurate flux distributions and performed better in the simulation of mutant phenotypes.

•

The flux distributions of the genome-scale models of Saccharomyces cerevisiae were evaluated

•

Most of the tested models showed fluxes inconsistent with experimental data

•

A manual curation process was performed on all reactions including NADH or NADPH

•

The curated models showed flux distributions more consistent with experimental data

•

Phenotype simulations improved when the curated flux distributions were used

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.

BioMet Toolbox 2.0: genome-wide analysis of metabolism and omics data

Analysis of large data sets using computational and mathematical tools have become a central part of biological sciences. Large amounts of data are being generated each year from different biological research fields leading to a constant development of software and algorithms aimed to deal with the increasing creation of information. The BioMet Toolbox 2.0 integrates a number of functionalities in a user-friendly environment enabling the user to work with biological data in a web interface. The unique and distinguishing feature of the BioMet Toolbox 2.0 is to provide a web user interface to tools for metabolic pathways and omics analysis developed under different platform-dependent environments enabling easy access to these computational tools.

Developing a new production host from a blueprint: Bacillus pumilus as an industrial enzyme producer

BACKGROUND

Since volatile and rising cost factors such as energy, raw materials and market competitiveness have a significant impact on the economic efficiency of biotechnological bulk productions, industrial processes need to be steadily improved and optimized. Thereby the current production hosts can undergo various limitations. To overcome those limitations and in addition increase the diversity of available production hosts for future applications, we suggest a Production Strain Blueprinting (PSB) strategy to develop new production systems in a reduced time lapse in contrast to a development from scratch.To demonstrate this approach, Bacillus pumilus has been developed as an alternative expression platform for the production of alkaline enzymes in reference to the established industrial production host Bacillus licheniformis.

RESULTS

To develop the selected B. pumilus as an alternative production host the suggested PSB strategy was applied proceeding in the following steps (dedicated product titers are scaled to the protease titer of Henkel's industrial production strain B. licheniformis at lab scale): Introduction of a protease production plasmid, adaptation of a protease production process (44%), process optimization (92%) and expression optimization (114%). To further evaluate the production capability of the developed B. pumilus platform, the target protease was substituted by an α-amylase. The expression performance was tested under the previously optimized protease process conditions and under subsequently adapted process conditions resulting in a maximum product titer of 65% in reference to B. licheniformis protease titer.

CONCLUSIONS

In this contribution the applied PSB strategy performed very well for the development of B. pumilus as an alternative production strain. Thereby the engineered B. pumilus expression platform even exceeded the protease titer of the industrial production host B. licheniformis by 14%. This result exhibits a remarkable potential of B. pumilus to be the basis for a next generation production host, since the strain has still a large potential for further genetic engineering. The final amylase titer of 65% in reference to B. licheniformis protease titer suggests that the developed B. pumilus expression platform is also suitable for an efficient production of non-proteolytic enzymes reaching a final titer of several grams per liter without complex process modifications.

Generalized framework for context-specific metabolic model extraction methods

Genome-scale metabolic models (GEMs) are increasingly applied to investigate the physiology not only of simple prokaryotes, but also eukaryotes, such as plants, characterized with compartmentalized cells of multiple types. While genome-scale models aim at including the entirety of known metabolic reactions, mounting evidence has indicated that only a subset of these reactions is active in a given context, including: developmental stage, cell type, or environment. As a result, several methods have been proposed to reconstruct context-specific models from existing genome-scale models by integrating various types of high-throughput data. Here we present a mathematical framework that puts all existing methods under one umbrella and provides the means to better understand their functioning, highlight similarities and differences, and to help users in selecting a most suitable method for an application.

Software platforms to facilitate reconstructing genome-scale metabolic networks

System-level analyses of microbial metabolism are facilitated by genome-scale reconstructions of microbial biochemical networks. A reconstruction provides a structured representation of the biochemical transformations occurring within an organism, as well as the genes necessary to carry out these transformations, as determined by the annotated genome sequence and experimental data. Network reconstructions also serve as platforms for constraint-based computational techniques, which facilitate biological studies in a variety of applications, including evaluation of network properties, metabolic engineering and drug discovery. Bottom-up metabolic network reconstructions have been developed for dozens of organisms, but until recently, the pace of reconstruction has failed to keep up with advances in genome sequencing. To address this problem, a number of software platforms have been developed to automate parts of the reconstruction process, thereby alleviating much of the manual effort previously required. Here, we review four such platforms in the context of established guidelines for network reconstruction. While many steps of the reconstruction process have been successfully automated, some manual evaluation of the results is still required to ensure a high-quality reconstruction. Widespread adoption of these platforms by the scientific community is underway and will be further enabled by exchangeable formats across platforms.