An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR)

Global regulator IrrE on stress tolerance: a review

Stress tolerance is a vital attribute for all living beings to cope with environmental adversities. IrrE (also named PprI) from Deinococcus radiodurans enhances resistance to extreme radiation stress by functioning as a global regulator, mediating the transcription of genes involved in deoxyribonucleic acid (DNA) damage response (DDR). The expression of IrrE augmented the resilience of various species to heat, radiation, oxidation, osmotic stresses and inhibitors, encompassing bacterial, fungal, plant, and mammalian cells. Moreover, IrrE was employed in a global regulator engineering strategy to broaden its applications in stress tolerance. The regulatory impacts of heterologously expressed IrrE have been investigated at the molecular and systems level, including the regulation of genes, proteins, modules, or pathways involved in DNA repair, detoxification proteins, protective molecules, native regulators and other aspects. In this review, we discuss the regulatory role and mechanism of IrrE in the antiradiation response of D. radiodurans. Furthermore, the applications and regulatory effects of heterologous expression of IrrE to enhance abiotic stress tolerance are summarized in particular.

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.

Genome-scale metabolic network models for industrial microorganisms metabolic engineering: Current advances and future prospects

The construction of high-performance microbial cell factories (MCFs) is the centerpiece of biomanufacturing. However, the complex metabolic regulatory network of microorganisms poses great challenges for the efficient design and construction of MCFs. The genome-scale metabolic network models (GSMs) can systematically simulate the metabolic regulation process of microorganisms in silico, providing effective guidance for the rapid design and construction of MCFs. In this review, we summarized the development status of 16 important industrial microbial GSMs, and further outline the technologies or methods that continuously promote high-quality GSMs construction from five aspects: I) Databases and modeling tools facilitate GSMs reconstruction; II) evolving gap-filling technologies; III) constraint-based model reconstruction; IV) advances in algorithms; and V) developed visualization tools. In addition, we also summarized the applications of GSMs in guiding metabolic engineering from four aspects: I) exploring and explaining metabolic features; II) predicting the effects of genetic perturbations on metabolism; III) predicting the optimal phenotype; IV) guiding cell factories construction in practical experiment. Finally, we discussed the development of GSMs, aiming to provide a reference for efficiently reconstructing GSMs and guiding metabolic engineering.

Dynamic flux balance analysis of high cell density fed-batch culture of Escherichia coli BL21 (DE3) with mass spectrometry-based spent media analysis

Dynamic flux balance analysis (FBA) allows estimation of intracellular reaction rates using organism-specific genome-scale metabolic models (GSMM) and by assuming instantaneous pseudo-steady states for processes that are inherently dynamic. This technique is well-suited for industrial bioprocesses employing complex media characterized by a hierarchy of substrate uptake and product secretion. However, knowledge of exchange rates of many components of the media would be required to obtain meaningful results. Here, we performed spent media analysis using mass spectrometry coupled with liquid and gas chromatography for a fed-batch, high-cell density cultivation of Escherichia coli BL21(DE3) expressing a recombinant protein. Time course measurements thus obtained for 246 metabolites were converted to instantaneous exchange rates. These were then used as constraints for dynamic FBA using a previously reported GSMM, thus providing insights into how the flux map evolves through the process. Changes in tri-carboxylic acid cycle fluxes correlated with the increased demand for energy during recombinant protein production. The results show how amino acids act as hubs for the synthesis of other cellular metabolites. Our results provide a deeper understanding of an industrial bioprocess and will have implications in further optimizing the process.

Cross-regulation between proteome reallocation and metabolic flux redistribution governs bacterial growth transition kinetics

Bacteria need to adjust their metabolism and protein synthesis simultaneously to adapt to changing nutrient conditions. It's still a grand challenge to predict how cells coordinate such adaptation due to the cross-regulation between the metabolic fluxes and the protein synthesis. Here we developed a dynamic Constrained Allocation Flux Balance Analysis method (dCAFBA), which integrates flux-controlled proteome allocation and protein limited flux balance analysis. This framework can predict the redistribution dynamics of metabolic fluxes without requiring detailed enzyme parameters. We reveal that during nutrient up-shifts, the calculated metabolic fluxes change in agreement with experimental measurements of enzyme protein dynamics. During nutrient down-shifts, we uncover a switch of metabolic bottleneck from carbon uptake proteins to metabolic enzymes, which disrupts the coordination between metabolic flux and their enzyme abundance. Our method provides a quantitative framework to investigate cellular metabolism under varying environments and reveals insights into bacterial adaptation strategies.

Microbial Pathway Thermodynamics: Stoichiometric Models Unveil Anabolic and Catabolic Processes

The biotechnological exploitation of microorganisms enables the use of metabolism for the production of economically valuable substances, such as drugs or food. It is, thus, unsurprising that the investigation of microbial metabolism and its regulation has been an active research field for many decades. As a result, several theories and techniques were developed that allow for the prediction of metabolic fluxes and yields as biotechnologically relevant output parameters. One important approach is to derive macrochemical equations that describe the overall metabolic conversion of an organism and basically treat microbial metabolism as a black box. The opposite approach is to include all known metabolic reactions of an organism to assemble a genome-scale metabolic model. Interestingly, both approaches are rather successful at characterizing and predicting the expected product yield. Over the years, macrochemical equations especially have been extensively characterized in terms of their thermodynamic properties. However, a common challenge when characterizing microbial metabolism by a single equation is to split this equation into two, describing the two modes of metabolism, anabolism and catabolism. Here, we present strategies to systematically identify separate equations for anabolism and catabolism. Based on metabolic models, we systematically identify all theoretically possible catabolic routes and determine their thermodynamic efficiency. We then show how anabolic routes can be derived, and we use these to approximate biomass yield. Finally, we challenge the view of metabolism as a linear energy converter, in which the free energy gradient of catabolism drives the anabolic reactions.

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

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

Genome-Scale Model of Rhizopus microsporus: Metabolic integration of a fungal holobiont with its bacterial and viral endosymbionts

Rhizopus microsporus often lives in association with bacterial and viral symbionts that alter its biology. This fungal model represents an example of the complex interactions established among diverse organisms in functional holobionts. We constructed a Genome-Scale Model (GSM) of the fungal-bacterial-viral holobiont (iHol). We employed a constraint-based method to calculate the metabolic fluxes to decipher the metabolic interactions of the symbionts with their host. Our computational analyses of iHol simulate the holobiont's growth and the production of the toxin rhizoxin. Analyses of the calculated fluxes between R. microsporus in symbiotic (iHol) versus asymbiotic conditions suggest that changes in the lipid and nucleotide metabolism of the host are necessary for the functionality of the holobiont. Glycerol plays a pivotal role in the fungal-bacterial metabolic interaction, as its production does not compromise fungal growth, and Mycetohabitans bacteria can efficiently consume it. Narnavirus RmNV-20S and RmNV-23S affected the nucleotide metabolism without impacting the fungal-bacterial symbiosis. Our analyses highlighted the metabolic stability of Mycetohabitans throughout its co-evolution with the fungal host. We also predicted changes in reactions of the bacterial metabolism required for the active production of rhizoxin. This iHol is the first GSM of a fungal holobiont.

The EcoCyc Database (2023)

EcoCyc is a bioinformatics database available online at EcoCyc.org that describes the genome and the biochemical machinery of Escherichia coli K-12 MG1655. The long-term goal of the project is to describe the complete molecular catalog of the E. coli cell, as well as the functions of each of its molecular parts, to facilitate a system-level understanding of E. coli. EcoCyc is an electronic reference source for E. coli biologists and for biologists who work with related microorganisms. The database includes information pages on each E. coli gene product, metabolite, reaction, operon, and metabolic pathway. The database also includes information on the regulation of gene expression, E. coli gene essentiality, and nutrient conditions that do or do not support the growth of E. coli. The website and downloadable software contain tools for the analysis of high-throughput data sets. In addition, a steady-state metabolic flux model is generated from each new version of EcoCyc and can be executed online. The model can predict metabolic flux rates, nutrient uptake rates, and growth rates for different gene knockouts and nutrient conditions. Data generated from a whole-cell model that is parameterized from the latest data on EcoCyc are also available. This review outlines the data content of EcoCyc and of the procedures by which this content is generated.

Evaluating E. coli genome-scale metabolic model accuracy with high-throughput mutant fitness data

The Escherichia coli genome-scale metabolic model (GEM) is an exemplar systems biology model for the simulation of cellular metabolism. Experimental validation of model predictions is essential to pinpoint uncertainty and ensure continued development of accurate models. Here, we quantified the accuracy of four subsequent E. coli GEMs using published mutant fitness data across thousands of genes and 25 different carbon sources. This evaluation demonstrated the utility of the area under a precision-recall curve relative to alternative accuracy metrics. An analysis of errors in the latest (iML1515) model identified several vitamins/cofactors that are likely available to mutants despite being absent from the experimental growth medium and highlighted isoenzyme gene-protein-reaction mapping as a key source of inaccurate predictions. A machine learning approach further identified metabolic fluxes through hydrogen ion exchange and specific central metabolism branch points as important determinants of model accuracy. This work outlines improved practices for the assessment of GEM accuracy with high-throughput mutant fitness data and highlights promising areas for future model refinement in E. coli and beyond.

Local flux coordination and global gene expression regulation in metabolic modeling

Genome-scale metabolic networks (GSMs) are fundamental systems biology representations of a cell's entire set of stoichiometrically balanced reactions. However, such static GSMs do not incorporate the functional organization of metabolic genes and their dynamic regulation (e.g., operons and regulons). Specifically, there are numerous topologically coupled local reactions through which fluxes are coordinated; the global growth state often dynamically regulates many gene expression of metabolic reactions via global transcription factor regulators. Here, we develop a GSM reconstruction method, Decrem, by integrating locally coupled reactions and global transcriptional regulation of metabolism by cell state. Decrem produces predictions of flux and growth rates, which are highly correlated with those experimentally measured in both wild-type and mutants of three model microorganisms Escherichia coli, Saccharomyces cerevisiae, and Bacillus subtilis under various conditions. More importantly, Decrem can also explain the observed growth rates by capturing the experimentally measured flux changes between wild-types and mutants. Overall, by identifying and incorporating locally organized and regulated functional modules into GSMs, Decrem achieves accurate predictions of phenotypes and has broad applications in bioengineering, synthetic biology, and microbial pathology.

Linear Response Theory of Evolved Metabolic Systems

Predicting cellular metabolic states is a central problem in biophysics. Conventional approaches, however, sensitively depend on the microscopic details of individual metabolic systems. In this Letter, we derived a universal linear relationship between the metabolic responses against nutrient conditions and metabolic inhibition, with the aid of a microeconomic theory. The relationship holds in arbitrary metabolic systems as long as the law of mass conservation stands, as supported by extensive numerical calculations. It offers quantitative predictions without prior knowledge of systems.

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

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

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.

A Protocol for the Automatic Construction of Highly Curated Genome-Scale Models of Human Metabolism

Genome-scale metabolic models (GEMs) have emerged as a tool to understand human metabolism from a holistic perspective with high relevance in the study of many diseases and in the metabolic engineering of human cell lines. GEM building relies on either automated processes that lack manual refinement and result in inaccurate models or manual curation, which is a time-consuming process that limits the continuous update of reliable GEMs. Here, we present a novel algorithm-aided protocol that overcomes these limitations and facilitates the continuous updating of highly curated GEMs. The algorithm enables the automatic curation and/or expansion of existing GEMs or generates a highly curated metabolic network based on current information retrieved from multiple databases in real time. This tool was applied to the latest reconstruction of human metabolism (Human1), generating a series of the human GEMs that improve and expand the reference model and generating the most extensive and comprehensive general reconstruction of human metabolism to date. The tool presented here goes beyond the current state of the art and paves the way for the automatic reconstruction of a highly curated, up-to-date GEM with high potential in computational biology as well as in multiple fields of biological science where metabolism is relevant.

Computer-Based Design of a Cell Factory for High-Yield Cytidine Production

Pyrimidine ribonucleotide de novo biosynthesis pathway (PRdnBP) is an important pathway to produce pyrimidine nucleosides. We attempted to systematically investigate PRdnBP in Escherichia coli with genome-scale metabolic models and utilized the models to guide strain design. The balance of central carbon metabolism and PRdnBP affected the production of cytidine from glucose. Using Bayesian metabolic flux analysis, the effect of modified PRdnBP on the metabolic network was analyzed. The acetate overflow became coupled with PRdnBP flux, while they were originally independent under oxygen-sufficient conditions. The coupling between cytidine production and acetate secretion in the modified strain was weakened by arcA deletion, which resulted in further improving the efficient accumulation of cytidine. In total, 1.28 g/L of cytidine with a yield of 0.26 g/g glucose was produced. The yield of cytidine produced by E. coli is higher than previous reports. Our strategy provides an effective attempt to find metabolic bottlenecks in genetically engineered bacteria by using flux coupling analysis.

Past, Present, and Future of Genome Modification in Escherichia coli

Escherichia coli K-12 is one of the most well-studied species of bacteria. This species, however, is much more difficult to modify by homologous recombination (HR) than other model microorganisms. Research on HR in E. coli has led to a better understanding of the molecular mechanisms of HR, resulting in technical improvements and rapid progress in genome research, and allowing whole-genome mutagenesis and large-scale genome modifications. Developments using λ Red (exo, bet, and gam) and CRISPR-Cas have made E. coli as amenable to genome modification as other model microorganisms, such as Saccharomyces cerevisiae and Bacillus subtilis. This review describes the history of recombination research in E. coli, as well as improvements in techniques for genome modification by HR. This review also describes the results of large-scale genome modification of E. coli using these technologies, including DNA synthesis and assembly. In addition, this article reviews recent advances in genome modification, considers future directions, and describes problems associated with the creation of cells by design.

Adaptation on xylose improves glucose-xylose co-utilization and ethanol production in a carbon catabolite repression (CCR) compromised ethanologenic strain

BACKGROUND

Sugar hydrolysates from lignocellulosic biomass are majorly composed of glucose and xylose that can be fermented to biofuels. Bacteria, despite having the natural ability to consume xylose are unable to consume it in presence of glucose due to a carbon catabolite repression (CCR) mechanism. This leads to overall reduced productivity as well as incomplete xylose utilization due to ethanol build-up from glucose utilization. In our effort to develop a strain for simultaneous fermentation of glucose and xylose into ethanol, we deleted ptsG in ethanologenic E. coli SSK42 to make it deficient in CCR and performed adaptive laboratory evolution to achieve accelerated growth rate, sugar consumption and ethanol production. Finally, we performed proteomics study to identify changes that might have been responsible for the observed improved phenotype of the evolved strain.

RESULTS

The parental strain of SSK42, i.e., wild-type E. coli B, did not co-utilize glucose and xylose as expected. After deleting the ptsG gene encoding the EIIBC^Glc subunit of PTS system, glucose consumption is severely affected in wild-type E. coli B. However, the ethanologenic, SSK42 strain, which was evolved in our earlier study on both glucose and xylose, didn't show such a drastic effect of EIIBC^Glc deletion, instead consumed glucose first, followed by xylose without delay for switching from one sugar to another. To improve growth on xylose and co-utilization capabilities, the ptsG deleted SSK42 was evolved on xylose. The strain evolved for 78 generations, strain SCD78, displayed significant co-utilization of glucose and xylose sugars. At the bioreactor level, the strain SCD78 produced 3-times the ethanol titer of the parent strain with significant glucose-xylose co-utilization. The rate of glucose and xylose consumption also increased 3.4-fold and 3-fold, respectively. Proteome data indicates significant upregulation of TCA cycle proteins, respiration-related proteins, and some transporters, which may have a role in increasing the total sugar consumption and co-utilization of sugars.

CONCLUSION

Through adaptive evolution, we have obtained a strain that has a significant glucose-xylose co-utilization phenotype with 3-fold higher total sugar consumption rate and ethanol production rate compared to the unevolved strain. This study also points out that adaptation on xylose is enough to impart glucose-xylose co-utilization property in CCR compromised ethanologenic strain SSK42.

The transcription factor network of E. coli steers global responses to shifts in RNAP concentration

The robustness and sensitivity of gene networks to environmental changes is critical for cell survival. How gene networks produce specific, chronologically ordered responses to genome-wide perturbations, while robustly maintaining homeostasis, remains an open question. We analysed if short- and mid-term genome-wide responses to shifts in RNA polymerase (RNAP) concentration are influenced by the known topology and logic of the transcription factor network (TFN) of Escherichia coli. We found that, at the gene cohort level, the magnitude of the single-gene, mid-term transcriptional responses to changes in RNAP concentration can be explained by the absolute difference between the gene's numbers of activating and repressing input transcription factors (TFs). Interestingly, this difference is strongly positively correlated with the number of input TFs of the gene. Meanwhile, short-term responses showed only weak influence from the TFN. Our results suggest that the global topological traits of the TFN of E. coli shape which gene cohorts respond to genome-wide stresses.

Genome-Wide Association Study Reveals Host Factors Affecting Conjugation in Escherichia coli

The emergence and dissemination of antibiotic resistance threaten the treatment of common bacterial infections. Resistance genes are often encoded on conjugative elements, which can be horizontally transferred to diverse bacteria. In order to delay conjugative transfer of resistance genes, more information is needed on the genetic determinants promoting conjugation. Here, we focus on which bacterial host factors in the donor assist transfer of conjugative plasmids. We introduced the broad-host-range plasmid pKJK10 into a diverse collection of 113 Escherichia coli strains and measured by flow cytometry how effectively each strain transfers its plasmid to a fixed E. coli recipient. Differences in conjugation efficiency of up to 2.7 and 3.8 orders of magnitude were observed after mating for 24 h and 48 h, respectively. These differences were linked to the underlying donor strain genetic variants in genome-wide association studies, thereby identifying candidate genes involved in conjugation. We confirmed the role of fliF, fliK, kefB and ucpA in the donor ability of conjugative elements by validating defects in the conjugation efficiency of the corresponding lab strain single-gene deletion mutants. Based on the known cellular functions of these genes, we suggest that the motility and the energy supply, the intracellular pH or salinity of the donor affect the efficiency of plasmid transfer. Overall, this work advances the search for targets for the development of conjugation inhibitors, which can be administered alongside antibiotics to more effectively treat bacterial infections.

Synthesis of C2-C4 diols from bioresources: Pathways and metabolic intervention strategies

Diols are important platform chemicals with extensive industrial applications in biopolymer synthesis, cosmetics, and fuels. The increased dependence on non-renewable sources to meet the energy requirement of the population raised issues regarding fossil fuel depletion and environmental impacts. The utilization of biological methods for the synthesis of diols by utilizing renewable resources such as glycerol and agro-residual wastes gained attention worldwide because of its advantages. Among these, biotransformation of 1,3-propanediol (1,3-PDO) and 2,3-butanediol (2,3-BDO) were extensively studied and at present, these diols are produced commercially in large scale with high yield. Many important isomers of C2-C4 diols lack natural synthetic pathways and development of chassis strains for the synthesis can be accomplished by adopting synthetic biology approaches. This current review depicts an overall idea about the pathways involved in C2-C4 diol production, metabolic intervention strategies and technologies in recent years.

DCcov: Repositioning of drugs and drug combinations for SARS-CoV-2 infected lung through constraint-based modeling

The 2019 coronavirus disease (COVID-19) became a worldwide pandemic with currently no approved effective antiviral drug. Flux balance analysis (FBA) is an efficient method to analyze metabolic networks. Here, FBA was applied on human lung cells infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to reposition metabolic drugs and drug combinations against the virus replication within the host tissue. Making use of expression datasets of infected lung tissue, genome-scale COVID-19-specific metabolic models were reconstructed. Then, host-specific essential genes and gene pairs were determined through in silico knockouts that permit reducing the viral biomass production without affecting the host biomass. Key pathways that are associated with COVID-19 severity in lung tissue are related to oxidative stress, ferroptosis, and pyrimidine metabolism. By in silico screening of Food and Drug Administration (FDA)-approved drugs on the putative disease-specific essential genes and gene pairs, 85 drugs and 52 drug combinations were predicted as promising candidates for COVID-19 (https://github.com/sysbiolux/DCcov).

Insights on the Advancements of In Silico Metabolic Studies of Succinic Acid Producing Microorganisms: A Review with Emphasis on Actinobacillus succinogenes

Succinic acid (SA) is one of the top candidate value-added chemicals that can be produced from biomass via microbial fermentation. A considerable number of cell factories have been proposed in the past two decades as native as well as non-native SA producers. Actinobacillus succinogenes is among the best and earliest known natural SA producers. However, its industrial application has not yet been realized due to various underlying challenges. Previous studies revealed that the optimization of environmental conditions alone could not entirely resolve these critical problems. On the other hand, microbial in silico metabolic modeling approaches have lately been the center of attention and have been applied for the efficient production of valuable commodities including SA. Then again, literature survey results indicated the absence of up-to-date reviews assessing this issue, specifically concerning SA production. Hence, this review was designed to discuss accomplishments and future perspectives of in silico studies on the metabolic capabilities of SA producers. Herein, research progress on SA and A. succinogenes, pathways involved in SA production, metabolic models of SA-producing microorganisms, and status, limitations and prospects on in silico studies of A. succinogenes were elaborated. All in all, this review is believed to provide insights to understand the current scenario and to develop efficient mathematical models for designing robust SA-producing microbial strains.

Path to improving the life cycle and quality of genome-scale models of metabolism

Genome-scale models of metabolism (GEMs) are key computational tools for the systems-level study of metabolic networks. Here, we describe the "GEM life cycle," which we subdivide into four stages: inception, maturation, specialization, and amalgamation. We show how different types of GEM reconstruction workflows fit in each stage and proceed to highlight two fundamental bottlenecks for GEM quality improvement: GEM maturation and content removal. We identify common characteristics contributing to increasing quality of maturing GEMs drawing from past independent GEM maturation efforts. We then shed some much-needed light on the latent and unrecognized but pervasive issue of content removal, demonstrating the substantial effects of model pruning on its solution space. Finally, we propose a novel framework for content removal and associated confidence-level assignment which will help guide future GEM development efforts, reduce duplication of effort across groups, potentially aid automated reconstruction platforms, and boost the reproducibility of model development.

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.

Active and machine learning-based approaches to rapidly enhance microbial chemical production

In order to make renewable fuels and chemicals from microbes, new methods are required to engineer microbes more intelligently. Computational approaches, to engineer strains for enhanced chemical production typically rely on detailed mechanistic models (e.g., kinetic/stoichiometric models of metabolism)-requiring many experimental datasets for their parameterization-while experimental methods may require screening large mutant libraries to explore the design space for the few mutants with desired behaviors. To address these limitations, we developed an active and machine learning approach (ActiveOpt) to intelligently guide experiments to arrive at an optimal phenotype with minimal measured datasets. ActiveOpt was applied to two separate case studies to evaluate its potential to increase valine yields and neurosporene productivity in Escherichia coli. In both the cases, ActiveOpt identified the best performing strain in fewer experiments than the case studies used. This work demonstrates that machine and active learning approaches have the potential to greatly facilitate metabolic engineering efforts to rapidly achieve its objectives.

Tolerance to NADH/NAD+ imbalance anticipates aging and anti-aging interventions

Redox couples coordinate cellular function, but the consequences of their imbalances are unclear. This is somewhat associated with the limitations of their experimental quantification. Here we circumvent these difficulties by presenting an approach that characterizes fitness-based tolerance profiles to redox couple imbalances using an in silico representation of metabolism. Focusing on the NADH/NAD⁺ redox couple in yeast, we demonstrate that reductive disequilibria generate metabolic syndromes comparable to those observed in cancer cells. The tolerance of yeast mutants to redox disequilibrium can also explain 30% of the variability in their experimentally measured chronological lifespan. Moreover, by predicting the significance of some metabolites to help stand imbalances, we correctly identify nutrients underlying mechanisms of pathology, lifespan-protecting molecules, or caloric restriction mimetics. Tolerance to redox imbalances becomes, in this way, a sound framework to recognize properties of the aging phenotype while providing a consistent biological rationale to assess anti-aging interventions.

•

We simulate how imbalances in NADH/NAD+ ratio modify cellular metabolic behavior

•

This reveals a mechanism to understand metabolic alterations at low growth rates

•

Tolerance to imbalance explains experimentally measured lifespan in yeast

•

We predict lifespan-protecting metabolites in yeast, animal, and human models

Hydrogen-driven microbial biogas upgrading: Advances, challenges and solutions

As a clean and renewable energy, biogas is an important alternative to fossil fuels. However, the high carbon dioxide (CO₂) content in biogas limits its value as a fuel. 'Biogas upgrading' is an advanced process which removes CO₂ from biogas, thereby converting biogas to biomethane, which has a higher commercial value. Microbial technologies offer a sustainable and cost-effective way to upgrade biogas, removing CO₂ using hydrogen (H₂) as electron donor, generated by surplus electricity from renewable wind or solar energy. Hydrogenotrophic methanogens can be applied to convert CO₂ with H₂ to methane (CH₄), or alternatively, homoacetogens can convert both CO₂ and H₂ into value-added chemicals. Here, we comprehensively review the current state of biogas generation and utilization, and describe the advances in biological, H₂-dependent biogas upgrading technologies, with particular attention to key challenges associated with the processes, e.g., metabolic limitations, low H₂ transfer rate, and finite CO₂ conversion rate. We also highlight several new strategies for overcoming technical barriers to achieve efficient CO₂ conversion, including process optimization to eliminate metabolic limitation, novel reactor designs to improve H₂ transfer rate and utilization efficiency, and employing advanced genetic engineering tools to generate more efficient microorganisms. The insights offered in this review will promote further exploration into microbial, H₂-driven biogas upgrading, towards addressing the global energy crisis and climate change associated with use of fossil fuels.

Genome-scale metabolic modelling when changes in environmental conditions affect biomass composition

Genome-scale metabolic modeling is an important tool in the study of metabolism by enhancing the collation of knowledge, interpretation of data, and prediction of metabolic capabilities. A frequent assumption in the use of genome-scale models is that the in vivo organism is evolved for optimal growth, where growth is represented by flux through a biomass objective function (BOF). While the specific composition of the BOF is crucial, its formulation is often inherited from similar organisms due to the experimental challenges associated with its proper determination.

A cell’s macro-molecular composition is not fixed and it responds to changes in environmental conditions. As a consequence, initiatives for the high-fidelity determination of cellular biomass composition have been launched. Thus, there is a need for a mathematical and computational framework capable of using multiple measurements of cellular biomass composition in different environments. Here, we propose two different computational approaches for directly addressing this challenge: Biomass Trade-off Weighting (BTW) and Higher-dimensional-plane InterPolation (HIP).

In lieu of experimental data on biomass composition-variation in response to changing nutrient environment, we assess the properties of BTW and HIP using three hypothetical, yet biologically plausible, BOFs for the Escherichia coli genome-scale metabolic model iML1515. We find that the BTW and HIP formulations have a significant impact on model performance and phenotypes. Furthermore, the BTW method generates larger growth rates in all environments when compared to HIP. Using acetate secretion and the respiratory quotient as proxies for phenotypic changes, we find marked differences between the methods as HIP generates BOFs more similar to a reference BOF than BTW. We conclude that the presented methods constitute a conceptual step in developing genome-scale metabolic modelling approaches capable of addressing the inherent dependence of cellular biomass composition on nutrient environments.

Changes in the environment promote changes in an organism’s metabolism. To achieve balanced growth states for near-optimal function, cells respond through metabolic rearrangements, which may influence the biosynthesis of metabolic precursors for building a cell’s molecular constituents. Therefore, it is necessary to take the dependence of biomass composition on environmental conditions into consideration. While measuring the biomass composition for some environments is possible, and should be done, it cannot be completed for all possible environments.

In this work, we propose two main approaches, BTW and HIP, for addressing the challenge of estimating biomass composition in response to environmental changes. We evaluate the phenotypic consequences of BTW and HIP by characterizing their effect on growth, secretion potential, respiratory efficiency, and gene essentiality of a cell.

Our work constitutes a first conceptual step in accounting for the influence of growth conditions on biomass composition, and in turn the biomass composition’s effect on metabolic phenotypic traits, within constraint-based modelling. As such, we believe it will improve the relevance of constraint-based methods in metabolic engineering and drug discovery, since the biosynthetic potential of microbes for generating industrially relevant products or drugs often is closely linked to their biomass composition.

Improving ethanol production by studying the effect of pH using a modified metabolic model and a systemic approach

pH is an important factor affecting the growth and production of microorganisms; especially, its effect on ethanologenic microorganisms. It can change the ionization state of metabolites via the change in the charge of their functional groups that may lead to metabolic alteration. Here, we estimated the ionization state of metabolites and balanced the charge of reactions in genome-scale metabolic models of Saccharomyces cerevisiae, Escherichia coli, and Zymomonas mobilis at pH levels 5, 6, and 7. The robustness analysis was first implemented to anticipate the effect of proton exchange flux on growth rates for the constructed metabolic models at various pH. In accordance with previous experimental reports, the models predict that Z. mobilis is more sensitive to pH rather than S. cerevisiae and the yeast is more regulated by pH rather than E. coli. Then, a systemic approach was proposed to predict the pH effect on metabolic change and to find effective reactions on ethanol production in S. cerevisiae. The correlated reactions with ethanol production at predicted optimal pH in a range of proton exchange rates determined by robustness analysis were identified using the Pearson correlation coefficient. Then, fluxes of these reactions were applied to cluster the various pHs by principal component analysis and to identify the role of these reactions on metabolic differentiation because of pH change. Finally, 12 reactions were selected for up and downregulation to improve ethanol production. Enzyme regulators of the selected reactions were identified using the BRENDA database and 11 selected regulators were screened and optimized via Plackett-Burman and two-level full factorial designs, respectively. The proposed approach has enhanced yields of ethanol from 0.18 to 0.36 mol/mol carbon. Hence, not only a comprehensive approach for understanding the effect of pH on metabolism was proposed in this study, but also it successfully introduced key manipulations for ethanol overproduction.

Understanding FBA Solutions under Multiple Nutrient Limitations

Genome-scale stoichiometric modeling methods, in particular Flux Balance Analysis (FBA) and variations thereof, are widely used to investigate cell metabolism and to optimize biotechnological processes. Given (1) a metabolic network, which can be reconstructed from an organism’s genome sequence, and (2) constraints on reaction rates, which may be based on measured nutrient uptake rates, FBA predicts which reactions maximize an objective flux, usually the production of cell components. Although FBA solutions may accurately predict the metabolic behavior of a cell, the actual flux predictions are often hard to interpret. This is especially the case for conditions with many constraints, such as for organisms growing in rich nutrient environments: it remains unclear why a certain solution was optimal. Here, we rationalize FBA solutions by explaining for which properties the optimal combination of metabolic strategies is selected. We provide a graphical formalism in which the selection of solutions can be visualized; we illustrate how this perspective provides a glimpse of the logic that underlies genome-scale modeling by applying our formalism to models of various sizes.

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.

New insights into the degradation of synthetic pollutants in contaminated environments

The environment is contaminated by synthetic contaminants owing to their extensive applications globally. Hence, the removal of synthetic pollutants (SPs) from the environment has received widespread attention. Different remediation technologies have been investigated for their abilities to eliminate SPs from the ecosystem; these include photocatalysis, sonochemical techniques, nanoremediation, and bioremediation. SPs, which can be organic or inorganic, can be degraded by microbial metabolism at contaminated sites. Owing to their diverse metabolisms, microbes can adapt to a wide variety of environments. Several microbial strains have been reported for their bioremediation potential concerning synthetic chemical compounds. The selection of potential strains for large-scale removal of organic pollutants is an important research priority. Additionally, novel microbial consortia have been found to be capable of efficient degradation owing to their combined and co-metabolic activities. Microbial engineering is one of the most prominent and promising techniques for providing new opportunities to develop proficient microorganisms for various biological processes; here, we have targeted the SP-degrading mechanisms of microorganisms. This review provides an in-depth discussion of microbial engineering techniques that are used to enhance the removal of both organic and inorganic pollutants from different contaminated environments and under different conditions. The degradation of these pollutants is investigated using abiotic and biotic approaches; interestingly, biotic approaches based on microbial methods are preferable owing to their high potential for pollutant removal and cost-effectiveness.

Controlling network ensembles

The field of optimal control typically requires the assumption of perfect knowledge of the system one desires to control, which is an unrealistic assumption for biological systems, or networks, typically affected by high levels of uncertainty. Here, we investigate the minimum energy control of network ensembles, which may take one of a number of possible realizations. We ensure the controller derived can perform the desired control with a tunable amount of accuracy and we study how the control energy and the overall control cost scale with the number of possible realizations. Our focus is in characterizing the solution of the optimal control problem in the limit in which the systems are drawn from a continuous distribution, and in particular, how to properly pose the weighting terms in the objective function. We verify the theory in three examples of interest: a unidirectional chain network with uncertain edge weights and self-loop weights, a network where each edge weight is drawn from a given distribution, and the Jacobian of the dynamics corresponding to the cell signaling network of autophagy in the presence of uncertain parameters.

Application of the control usually requires complete knowledge of the system, which is rare for biological networks characterized by uncertainty. Klickstein et al. propose an optimal control for uncertain systems represented by network ensembles where only weight distributions for edges are known.

L-Alanine Prototrophic Suppressors Emerge from L-Alanine Auxotroph through Stress-Induced Mutagenesis in Escherichia coli

In Escherichia coli, L-alanine is synthesized by three isozymes: YfbQ, YfdZ, and AvtA. When an E. coli L-alanine auxotrophic isogenic mutant lacking the three isozymes was grown on L-alanine-deficient minimal agar medium, L-alanine prototrophic mutants emerged considerably more frequently than by spontaneous mutation; the emergence frequency increased over time, and, in an L-alanine-supplemented minimal medium, correlated inversely with L-alanine concentration, indicating that the mutants were derived through stress-induced mutagenesis. Whole-genome analysis of 40 independent L-alanine prototrophic mutants identified 16 and 18 clones harboring point mutation(s) in pyruvate dehydrogenase complex and phosphotransacetylase-acetate kinase pathway, which respectively produce acetyl coenzyme A and acetate from pyruvate. When two point mutations identified in L-alanine prototrophic mutants, in pta (D656A) and aceE (G147D), were individually introduced into the original L-alanine auxotroph, the isogenic mutants exhibited almost identical growth recovery as the respective cognate mutants. Each original- and isogenic-clone pair carrying the pta or aceE mutation showed extremely low phosphotransacetylase or pyruvate dehydrogenase activity, respectively. Lastly, extracellularly-added pyruvate, which dose-dependently supported L-alanine auxotroph growth, relieved the L-alanine starvation stress, preventing the emergence of L-alanine prototrophic mutants. Thus, L-alanine starvation-provoked stress-induced mutagenesis in the L-alanine auxotroph could lead to intracellular pyruvate increase, which eventually induces L-alanine prototrophy.

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

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

Balancing glucose and oxygen uptake rates to enable high amorpha-4,11-diene production in Escherichia coli via the methylerythritol phosphate pathway

Amorpha-4,11-diene (AMD4,11) is a precursor to artemisinin, a potent antimalarial drug that is traditionally extracted from the shrubs of Artemisia annua. Despite significant prior efforts to produce artemisinin and its precursors through biotechnology, there remains a dire need for more efficient biosynthetic routes for its production. Here, we describe the optimization of key process conditions for an Escherichia coli strain producing AMD4,11 via the native methylerythritol phosphate (MEP) pathway. By studying the interplay between glucose uptake rates and oxygen demand, we were able to identify optimal conditions for increasing carbon flux through the MEP pathway by manipulating the availability of NADPH required for terpenoid production. Installation of an optimal q_O2 /q_glucose led to a 6.7-fold increase in product titers and a 6.5-fold increase in carbon yield.

The ModelSEED Biochemistry Database for the integration of metabolic annotations and the reconstruction, comparison and analysis of metabolic models for plants, fungi and microbes

For over 10 years, ModelSEED has been a primary resource for the construction of draft genome-scale metabolic models based on annotated microbial or plant genomes. Now being released, the biochemistry database serves as the foundation of biochemical data underlying ModelSEED and KBase. The biochemistry database embodies several properties that, taken together, distinguish it from other published biochemistry resources by: (i) including compartmentalization, transport reactions, charged molecules and proton balancing on reactions; (ii) being extensible by the user community, with all data stored in GitHub; and (iii) design as a biochemical 'Rosetta Stone' to facilitate comparison and integration of annotations from many different tools and databases. The database was constructed by combining chemical data from many resources, applying standard transformations, identifying redundancies and computing thermodynamic properties. The ModelSEED biochemistry is continually tested using flux balance analysis to ensure the biochemical network is modeling-ready and capable of simulating diverse phenotypes. Ontologies can be designed to aid in comparing and reconciling metabolic reconstructions that differ in how they represent various metabolic pathways. ModelSEED now includes 33,978 compounds and 36,645 reactions, available as a set of extensible files on GitHub, and available to search at https://modelseed.org/biochem and KBase.

The ModelSEED Biochemistry Database for the integration of metabolic annotations and the reconstruction, comparison and analysis of metabolic models for plants, fungi and microbes

For over ten years, ModelSEED has been a primary resource for the construction of draft genome-scale metabolic models based on annotated microbial or plant genomes. Now being released, the biochemistry database serves as the foundation of biochemical data underlying ModelSEED and KBase. The biochemistry database embodies several properties that, taken together, distinguish it from other published biochemistry resources by: (i) including compartmentalization, transport reactions, charged molecules and proton balancing on reactions;; (ii) being extensible by the user community, with all data stored in GitHub; and (iii) design as a biochemical “Rosetta Stone” to facilitate comparison and integration of annotations from many different tools and databases. The database was constructed by combining chemical data from many resources, applying standard transformations, identifying redundancies, and computing thermodynamic properties. The ModelSEED biochemistry is continually tested using flux balance analysis to ensure the biochemical network is modeling-ready and capable of simulating diverse phenotypes. Ontologies can be designed to aid in comparing and reconciling metabolic reconstructions that differ in how they represent various metabolic pathways. ModelSEED now includes 33,978 compounds and 36,645 reactions, available as a set of extensible files on GitHub, and available to search at https://modelseed.org and KBase.

Genome-scale metabolic modeling underscores the potential of Cutaneotrichosporon oleaginosus ATCC 20509 as a cell factory for biofuel production

BACKGROUND

Cutaneotrichosporon oleaginosus ATCC 20509 is a fast-growing oleaginous basidiomycete yeast that is able to grow in a wide range of low-cost carbon sources including crude glycerol, a byproduct of biodiesel production. When glycerol is used as a carbon source, this yeast can accumulate more than 50% lipids (w/w) with high concentrations of mono-unsaturated fatty acids.

RESULTS

To increase our understanding of this yeast and to provide a knowledge base for further industrial use, a FAIR re-annotated genome was used to build a genome-scale, constraint-based metabolic model containing 1553 reactions involving 1373 metabolites in 11 compartments. A new description of the biomass synthesis reaction was introduced to account for massive lipid accumulation in conditions with high carbon-to-nitrogen (C/N) ratio in the media. This condition-specific biomass objective function is shown to better predict conditions with high lipid accumulation using glucose, fructose, sucrose, xylose, and glycerol as sole carbon source.

CONCLUSION

Contributing to the economic viability of biodiesel as renewable fuel, C. oleaginosus ATCC 20509 can effectively convert crude glycerol waste streams in lipids as a potential bioenergy source. Performance simulations are essential to identify optimal production conditions and to develop and fine tune a cost-effective production process. Our model suggests ATP-citrate lyase as a possible target to further improve lipid production.

Unlocking Elementary Conversion Modes: ecmtool Unveils All Capabilities of Metabolic Networks

The metabolic capabilities of cells determine their biotechnological potential, fitness in ecosystems, pathogenic threat levels, and function in multicellular organisms. Their comprehensive experimental characterization is generally not feasible, particularly for unculturable organisms. In principle, the full range of metabolic capabilities can be computed from an organism's annotated genome using metabolic network reconstruction. However, current computational methods cannot deal with genome-scale metabolic networks. Part of the problem is that these methods aim to enumerate all metabolic pathways, while computation of all (elementally balanced) conversions between nutrients and products would suffice. Indeed, the elementary conversion modes (ECMs, defined by Urbanczik and Wagner) capture the full metabolic capabilities of a network, but the use of ECMs has not been accessible until now. We explain and extend the theory of ECMs, implement their enumeration in ecmtool, and illustrate their applicability. This work contributes to the elucidation of the full metabolic footprint of any cell.

•

Elementary conversion modes (ECMs) specify all metabolic capabilities of any organism

•

Ecmtool computes all ECMs from a reconstructed metabolic network

•

ECM enumeration enables metabolic characterization of larger networks than ever

•

Focusing on ECMs between relevant metabolites even enables genome-scale enumeration

Understanding the metabolic capabilities of cells is of profound importance. Microbial metabolism shapes global cycles of elements and cleans polluted soils. Human and pathogen metabolism affects our health. Recent advances allow for automatic reconstruction of reaction networks for any organism, which is already used in synthetic biology, (food) microbiology, and agriculture to compute optimal yields from resources to products. However, computational tools are limited to optimal states or subnetworks, leaving many capabilities of organisms hidden. Our program, ecmtool, creates a blueprint of any organism's metabolic functionalities, drastically improving insights obtained from genome sequences. Ecmtool may become essential in exploratory research, especially for studying cells that are not culturable in laboratory conditions. Ideally, elementary conversion mode enumeration will someday be a standard step after metabolic network reconstruction, achieving the metabolic characterization of all known organisms.

Metabolic flux configuration determination using information entropy

Constraint-based models use steady-state mass balances to define a solution space of flux configurations, which can be narrowed down by measuring as many fluxes as possible. Due to loops and redundant pathways, this process typically yields multiple alternative solutions. To address this ambiguity, flux sampling can estimate the probability distribution of each flux, or a flux configuration can be singled out by further minimizing the sum of fluxes according to the assumption that cellular metabolism favors states where enzyme-related costs are economized. However, flux sampling is susceptible to artifacts introduced by thermodynamically infeasible cycles and is it not clear if the economy of fluxes assumption (EFA) is universally valid. Here, we formulated a constraint-based approach, MaxEnt, based on the principle of maximum entropy, which in this context states that if more than one flux configuration is consistent with a set of experimentally measured fluxes, then the one with the minimum amount of unwarranted assumptions corresponds to the best estimation of the non-observed fluxes. We compared MaxEnt predictions to Escherichia coli and Saccharomyces cerevisiae publicly available flux data. We found that the mean square error (MSE) between experimental and predicted fluxes by MaxEnt and EFA-based methods are three orders of magnitude lower than the median of 1,350,000 MSE values obtained using flux sampling. However, only MaxEnt and flux sampling correctly predicted flux through E. coli’s glyoxylate cycle, whereas EFA-based methods, in general, predict no flux cycles. We also tested MaxEnt predictions at increasing levels of overflow metabolism. We found that MaxEnt accuracy is not affected by overflow metabolism levels, whereas the EFA-based methods show a decreasing performance. These results suggest that MaxEnt is less sensitive than flux sampling to artifacts introduced by thermodynamically infeasible cycles and that its predictions are less susceptible to overfitting than EFA-based methods.

Genome-Scale Metabolic Modeling for Unraveling Molecular Mechanisms of High Threat Pathogens

Pathogens give rise to a wide range of diseases threatening global health and hence drawing public health agencies' attention to establish preventative and curative solutions. Genome-scale metabolic modeling is ever increasingly used tool for biomedical applications including the elucidation of antibiotic resistance, virulence, single pathogen mechanisms and pathogen-host interaction systems. With this approach, the sophisticated cellular system of metabolic reactions inside the pathogens as well as between pathogen and host cells are represented in conjunction with their corresponding genes and enzymes. Along with essential metabolic reactions, alternate pathways and fluxes are predicted by performing computational flux analyses for the growth of pathogens in a very short time. The genes or enzymes responsible for the essential metabolic reactions in pathogen growth are regarded as potential drug targets, as a priori guide to researchers in the pharmaceutical field. Pathogens alter the key metabolic processes in infected host, ultimately the objective of these integrative constraint-based context-specific metabolic models is to provide novel insights toward understanding the metabolic basis of the acute and chronic processes of infection, revealing cellular mechanisms of pathogenesis, identifying strain-specific biomarkers and developing new therapeutic approaches including the combination drugs. The reaction rates predicted during different time points of pathogen development enable us to predict active pathways and those that only occur during certain stages of infection, and thus point out the putative drug targets. Among others, fatty acid and lipid syntheses reactions are recent targets of new antimicrobial drugs. Genome-scale metabolic models provide an improved understanding of how intracellular pathogens utilize the existing microenvironment of the host. Here, we reviewed the current knowledge of genome-scale metabolic modeling in pathogen cells as well as pathogen host interaction systems and the promising applications in the extension of curative strategies against pathogens for global preventative healthcare.

Genome-scale metabolic models highlight stage-specific differences in essential metabolic pathways in Trypanosoma cruzi

Chagas disease is a neglected tropical disease and a leading cause of heart failure in Latin America caused by a protozoan called Trypanosoma cruzi. This parasite presents a complex multi-stage life cycle. Anti-Chagas drugs currently available are limited to benznidazole and nifurtimox, both with severe side effects. Thus, there is a need for alternative and more efficient drugs. Genome-scale metabolic models (GEMs) can accurately predict metabolic capabilities and aid in drug discovery in metabolic genes. This work developed an extended GEM, hereafter referred to as iIS312, of the published and validated T. cruzi core metabolism model. From iIS312, we then built three stage-specific models through transcriptomics data integration, and showed that epimastigotes present the most active metabolism among the stages (see S1-S4 GEMs). Stage-specific models predicted significant metabolic differences among stages, including variations in flux distribution in core metabolism. Moreover, the gene essentiality predictions suggest potential drug targets, among which some have been previously proven lethal, including glutamate dehydrogenase, glucokinase and hexokinase. To validate the models, we measured the activity of enzymes in the core metabolism of the parasite at different stages, and showed the results were consistent with model predictions. Our results represent a potential step forward towards the improvement of Chagas disease treatment. To our knowledge, these stage-specific models are the first GEMs built for the stages Amastigote and Trypomastigote. This work is also the first to present an in silico GEM comparison among different stages in the T. cruzi life cycle.

Genome-scale reconstruction of Paenarthrobacter aurescens TC1 metabolic model towards the study of atrazine bioremediation

Atrazine is an herbicide and a pollutant of great environmental concern that is naturally biodegraded by microbial communities. Paenarthrobacter aurescens TC1 is one of the most studied degraders of this herbicide. Here, we developed a genome scale metabolic model for P. aurescens TC1, iRZ1179, to study the atrazine degradation process at organism level. Constraint based flux balance analysis and time dependent simulations were used to explore the organism’s phenotypic landscape. Simulations aimed at designing media optimized for supporting growth and enhancing degradation, by passing the need in strain design via genetic modifications. Growth and degradation simulations were carried with more than 100 compounds consumed by P. aurescens TC1. In vitro validation confirmed the predicted classification of different compounds as efficient, moderate or poor stimulators of growth. Simulations successfully captured previous reports on the use of glucose and phosphate as bio-stimulators of atrazine degradation, supported by in vitro validation. Model predictions can go beyond supplementing the medium with a single compound and can predict the growth outcomes for higher complexity combinations. Hence, the analysis demonstrates that the exhaustive power of the genome scale metabolic reconstruction allows capturing complexities that are beyond common biochemical expertise and knowledge and further support the importance of computational platforms for the educated design of complex media. The model presented here can potentially serve as a predictive tool towards achieving optimal biodegradation efficiencies and for the development of ecologically friendly solutions for pollutant degradation.

Biosystems Design by Machine Learning

Biosystems such as enzymes, pathways, and whole cells have been increasingly explored for biotechnological applications. However, the intricate connectivity and resulting complexity of biosystems poses a major hurdle in designing biosystems with desirable features. As -omics and other high throughput technologies have been rapidly developed, the promise of applying machine learning (ML) techniques in biosystems design has started to become a reality. ML models enable the identification of patterns within complicated biological data across multiple scales of analysis and can augment biosystems design applications by predicting new candidates for optimized performance. ML is being used at every stage of biosystems design to help find nonobvious engineering solutions with fewer design iterations. In this review, we first describe commonly used models and modeling paradigms within ML. We then discuss some applications of these models that have already shown success in biotechnological applications. Moreover, we discuss successful applications at all scales of biosystems design, including nucleic acids, genetic circuits, proteins, pathways, genomes, and bioprocesses. Finally, we discuss some limitations of these methods and potential solutions as well as prospects of the combination of ML and biosystems design.