Systems Biology of Metabolism.

Metabolism is highly complex and involves thousands of different connected reactions; it is therefore necessary to use mathematical models for holistic studies. The use of mathematical models in biology is referred to as systems biology. In this review, the principles of systems biology are described, and two different types of mathematical models used for studying metabolism are discussed: kinetic models and genome-scale metabolic models. The use of different omics technologies, including transcriptomics, proteomics, metabolomics, and fluxomics, for studying metabolism is presented. Finally, the application of systems biology for analyzing global regulatory structures, engineering the metabolism of cell factories, and analyzing human diseases is discussed.

Review and meta-analysis of the genetic Minimal Cut Set approach for gene essentiality prediction in cancer metabolism.

Cancer metabolism is a marvellously complex topic, in part, due to the reprogramming of its pathways to self-sustain the malignant phenotype in the disease, to the detriment of its healthy counterpart. Understanding these adjustments can provide novel targeted therapies that could disrupt and impair proliferation of cancerous cells. For this very purpose, genome-scale metabolic models (GEMs) have been developed, with Human1 being the most recent reconstruction of the human metabolism. Based on GEMs, we introduced the genetic Minimal Cut Set (gMCS) approach, an uncontextualized methodology that exploits the concepts of synthetic lethality to predict metabolic vulnerabilities in cancer. gMCSs define a set of genes whose knockout would render the cell unviable by disrupting an essential metabolic task in GEMs, thus, making cellular proliferation impossible. Here, we summarize the gMCS approach and review the current state of the methodology by performing a systematic meta-analysis based on two datasets of gene essentiality in cancer. First, we assess several thresholds and distinct methodologies for discerning highly and lowly expressed genes. Then, we address the premise that gMCSs of distinct length should have the same predictive power. Finally, we question the importance of a gene partaking in multiple gMCSs and analyze the importance of all the essential metabolic tasks defined in Human1. Our meta-analysis resulted in parameter evaluation to increase the predictive power for the gMCS approach, as well as a significant reduction of computation times by only selecting the crucial gMCS lengths, proposing the pertinency of particular parameters for the peak processing of gMCS.

Systems-level temporal immune-metabolic profile in Crimean-Congo hemorrhagic fever virus infection.

Crimean-Congo hemorrhagic fever (CCHF) caused by CCHF virus (CCHFV) is one of the epidemic-prone diseases prioritized by the World Health Organisation as public health emergency with an urgent need for accelerated research. The trajectory of host response against CCHFV is multifarious and remains unknown. Here, we reported the temporal spectrum of pathogenesis following the CCHFV infection using genome-wide blood transcriptomics analysis followed by advanced systems biology analysis, temporal immune-pathogenic alterations, and context-specific progressive and postinfection genome-scale metabolic models (GSMM) on samples collected during the acute (T0), early convalescent (T1), and convalescent-phase (T2). The interplay between the retinoic acid-inducible gene-I-like/nucleotide-binding oligomerization domain-like receptor and tumor necrosis factor signaling governed the trajectory of antiviral immune responses. The rearrangement of intracellular metabolic fluxes toward the amino acid metabolism and metabolic shift toward oxidative phosphorylation and fatty acid oxidation during acute CCHFV infection determine the pathogenicity. The upregulation of the tricarboxylic acid cycle during CCHFV infection, compared to the noninfected healthy control and between the severity groups, indicated an increased energy demand and cellular stress. The upregulation of glycolysis and pyruvate metabolism potentiated energy generation through alternative pathways associated with the severity of the infection. The downregulation of metabolic processes at the convalescent phase identified by blood cell transcriptomics and single-cell type proteomics of five immune cells (CD4+ and CD8+ T cells, CD14+ monocytes, B cells, and NK cells) potentially leads to metabolic rewiring through the recovery due to hyperactivity during the acute phase leading to post-viral fatigue syndrome.

Applications of genome-scale metabolic models to investigate microbial metabolic adaptations in response to genetic or environmental perturbations.

Environmental perturbations are encountered by microorganisms regularly and will require metabolic adaptations to ensure an organism can survive in the newly presenting conditions. In order to study the mechanisms of metabolic adaptation in such conditions, various experimental and computational approaches have been used. Genome-scale metabolic models (GEMs) are one of the most powerful approaches to study metabolism, providing a platform to study the systems level adaptations of an organism to different environments which could otherwise be infeasible experimentally. In this review, we are describing the application of GEMs in understanding how microbes reprogram their metabolic system as a result of environmental variation. In particular, we provide the details of metabolic model reconstruction approaches, various algorithms and tools for model simulation, consequences of genetic perturbations, integration of '-omics' datasets for creating context-specific models and their application in studying metabolic adaptation due to the change in environmental conditions.

Yeast9: a consensus genome-scale metabolic model for S. cerevisiae curated by the community.

Genome-scale metabolic models (GEMs) can facilitate metabolism-focused multi-omics integrative analysis. Since Yeast8, the yeast-GEM of Saccharomyces cerevisiae, published in 2019, has been continuously updated by the community. This has increased the quality and scope of the model, culminating now in Yeast9. To evaluate its predictive performance, we generated 163 condition-specific GEMs constrained by single-cell transcriptomics from osmotic pressure or reference conditions. Comparative flux analysis showed that yeast adapting to high osmotic pressure benefits from upregulating fluxes through central carbon metabolism. Furthermore, combining Yeast9 with proteomics revealed metabolic rewiring underlying its preference for nitrogen sources. Lastly, we created strain-specific GEMs (ssGEMs) constrained by transcriptomics for 1229 mutant strains. Well able to predict the strains' growth rates, fluxomics from those large-scale ssGEMs outperformed transcriptomics in predicting functional categories for all studied genes in machine learning models. Based on those findings we anticipate that Yeast9 will continue to empower systems biology studies of yeast metabolism.

Metastatic triple negative breast cancer adapts its metabolism to destination tissues while retaining key metabolic signatures.

Triple negative breast cancer (TNBC) metastases are assumed to exhibit similar functions in different organs as in the original primary tumor. However, studies of metastasis are often limited to a comparison of metastatic tumors with primary tumors of their origin, and little is known about the adaptation to the local environment of the metastatic sites. We therefore used transcriptomic data and metabolic network analyses to investigate whether metastatic tumors adapt their metabolism to the metastatic site and found that metastatic tumors adopt a metabolic signature with some similarity to primary tumors of their destinations. The extent of adaptation, however, varies across different organs, and metastatic tumors retain metabolic signatures associated with TNBC. Our findings suggest that a combination of anti-metastatic approaches and metabolic inhibitors selected specifically for different metastatic sites, rather than solely targeting TNBC primary tumors, may constitute a more effective treatment approach.

LooplessFluxSampler: an efficient toolbox for sampling the loopless flux solution space of metabolic models.

BACKGROUND: Uniform random sampling of mass-balanced flux solutions offers an unbiased appraisal of the capabilities of metabolic networks. Unfortunately, it is impossible to avoid thermodynamically infeasible loops in flux samples when using convex samplers on large metabolic models. Current strategies for randomly sampling the non-convex loopless flux space display limited efficiency and lack theoretical guarantees. RESULTS: Here, we present LooplessFluxSampler, an efficient algorithm for exploring the loopless mass-balanced flux solution space of metabolic models, based on an Adaptive Directions Sampling on a Box (ADSB) algorithm. ADSB is rooted in the general Adaptive Direction Sampling (ADS) framework, specifically the Parallel ADS, for which theoretical convergence and irreducibility results are available for sampling from arbitrary distributions. By sampling directions that adapt to the target distribution, ADSB traverses more efficiently the sample space achieving faster mixing than other methods. Importantly, the presented algorithm is guaranteed to target the uniform distribution over convex regions, and it provably converges on the latter distribution over more general (non-convex) regions provided the sample can have full support. CONCLUSIONS: LooplessFluxSampler enables scalable statistical inference of the loopless mass-balanced solution space of large metabolic models. Grounded in a theoretically sound framework, this toolbox provides not only efficient but also reliable results for exploring the properties of the almost surely non-convex loopless flux space. Finally, LooplessFluxSampler includes a Markov Chain diagnostics suite for assessing the quality of the final sample and the performance of the algorithm.

Insights into yeast response to chemotherapeutic agent through time series genome-scale metabolic models.

Organism-specific genome-scale metabolic models (GSMMs) can unveil molecular mechanisms within cells and are commonly used in diverse applications, from synthetic biology, biotechnology, and systems biology to metabolic engineering. There are limited studies incorporating time-series transcriptomics in GSMM simulations. Yeast is an easy-to-manipulate model organism for tumor research. Here, a novel approach (TS-GSMM) was proposed to integrate time-series transcriptomics with GSMMs to narrow down the feasible solution space of all possible flux distributions and attain time-series flux samples. The flux samples were clustered using machine learning techniques, and the clusters' functional analysis was performed using reaction set enrichment analysis. A time series transcriptomics response of Yeast cells to a chemotherapeutic reagent-doxorubicin-was mapped onto a Yeast GSMM. Eleven flux clusters were obtained with our approach, and pathway dynamics were displayed. Induction of fluxes related to bicarbonate formation and transport, ergosterol and spermidine transport, and ATP production were captured. Integrating time-series transcriptomics data with GSMMs is a promising approach to reveal pathway dynamics without any kinetic modeling and detects pathways that cannot be identified through transcriptomics-only analysis. The codes are available at https://github.com/karabekmez/TS-GSMM.

The role of bacterial metabolism in human gut colonization.

Can we anticipate the emergence of the next pandemic antibiotic-resistant bacterial clone? Addressing such an ambitious question relies on our ability to comprehensively understand the ecological and epidemiological factors fostering the evolution of high-risk clones. Among these factors, the ability to persistently colonize and thrive in the human gut is crucial for most high-risk clones. Nonetheless, the causes and mechanisms facilitating successful gut colonization remain obscure. Here, we review recent evidence that suggests that bacterial metabolism plays a pivotal role in determining the ability of high-risk clones to colonize the human gut. Subsequently, we outline novel approaches that enable the exploration of microbial metabolism at an unprecedented scale and level of detail. A thorough understanding of the constraints and opportunities of bacterial metabolism in gut colonization will foster our ability to predict the emergence of high-risk clones and take appropriate containment strategies.

Integrating Genome-Scale Metabolic Models with Patient Plasma Metabolome to Study Endothelial Metabolism In Situ.

Patient blood samples are invaluable in clinical omics databases, yet current methodologies often fail to fully uncover the molecular mechanisms driving patient pathology. While genome-scale metabolic models (GEMs) show promise in systems medicine by integrating various omics data, having only exometabolomic data remains a limiting factor. To address this gap, we introduce a comprehensive pipeline integrating GEMs with patient plasma metabolome. This pipeline constructs case-specific GEMs using literature-based and patient-specific metabolomic data. Novel computational methods, including adaptive sampling and an in-house developed algorithm for the rational exploration of the sampled space of solutions, enhance integration accuracy while improving computational performance. Model characterization involves task analysis in combination with clustering methods to identify critical cellular functions. The new pipeline was applied to a cohort of trauma patients to investigate shock-induced endotheliopathy using patient plasma metabolome data. By analyzing endothelial cell metabolism comprehensively, the pipeline identified critical therapeutic targets and biomarkers that can potentially contribute to the development of therapeutic strategies. Our study demonstrates the efficacy of integrating patient plasma metabolome data into computational models to analyze endothelial cell metabolism in disease contexts. This approach offers a deeper understanding of metabolic dysregulations and provides insights into diseases with metabolic components and potential treatments.

A reactive species reactions module for integration into genome-scale metabolic models for improved insights: Application to cancer.

Reactive species (RS) play significant roles in many disease contexts. Despite their crucial roles in diseases including cancer, the RS are not adequately modeled in the genome-scale metabolic (GSM) models, which are used to understand cell metabolism in disease contexts. We have developed a scalable RS reactions module that can be integrated with any Recon 3D-derived human metabolic model, or after fine-tuning, with any metabolic model. With RS-integration, the GSM models of three cancers (basal-like triple negative breast cancer (TNBC), high grade serous ovarian carcinoma (HGSOC) and colorectal cancer (CRC)) built from Recon 3D, precisely highlighted the increases/decreases in fluxes (dysregulation) occurring in important pathways of these cancers. These dysregulations were not prominent in the standard cancer models without the RS module. Further, the results from these RS-integrated cancer GSM models suggest the following decreasing order in the ease of ferroptosis-targeting to treat the cancers: TNBC > HGSOC > CRC.

High-quality genome-scale metabolic network reconstruction of probiotic bacterium Escherichia coli Nissle 1917.

BACKGROUND: Escherichia coli Nissle 1917 (EcN) is a probiotic bacterium used to treat various gastrointestinal diseases. EcN is increasingly being used as a chassis for the engineering of advanced microbiome therapeutics. To aid in future engineering efforts, our aim was to construct an updated metabolic model of EcN with extended secondary metabolite representation. RESULTS: An updated high-quality genome-scale metabolic model of EcN, iHM1533, was developed based on comparison with 55 E. coli/Shigella reference GEMs and manual curation, including expanded secondary metabolite pathways (enterobactin, salmochelins, aerobactin, yersiniabactin, and colibactin). The model was validated and improved using phenotype microarray data, resulting in an 82.3% accuracy in predicting growth phenotypes on various nutrition sources. Flux variability analysis with previously published 13C fluxomics data validated prediction of the internal central carbon fluxes. A standardised test suite called Memote assessed the quality of iHM1533 to have an overall score of 89%. The model was applied by using constraint-based flux analysis to predict targets for optimisation of secondary metabolite production. Modelling predicted design targets from across amino acid metabolism, carbon metabolism, and other subsystems that are common or unique for influencing the production of various secondary metabolites. CONCLUSION: iHM1533 represents a well-annotated metabolic model of EcN with extended secondary metabolite representation. Phenotype characterisation and the iHM1533 model provide a better understanding of the metabolic capabilities of EcN and will help future metabolic engineering efforts.

Yeast metabolic innovations emerged via expanded metabolic network and gene positive selection.

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome-scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence-level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.

Pareto optimal metabolic engineering for the growth-coupled overproduction of sustainable chemicals.

Our research aims to help industrial biotechnology develop a sustainable economy using green technology based on microorganisms and synthetic biology through two case studies that improve metabolic capacity in yeast models Yarrowia lipolytica (Y. lipolytica) and Saccharomyces cerevisiae (S. cerevisiae). We aim to increase the production capacity of beta-carotene (ß-carotene) and succinic acid, which are among the highest market demands due to their versatile use in numerous consumer products. We performed simulations to identify in silico ranking of strains based on multiple objectives: the growth rate of yeast microorganisms, the number of used chromosomes, and the production capability of ß-carotene (for Y. lipolytica) and succinate (for S. cerevisiae). Our multiobjective optimization methodology identified notable gene deletions by searching a vast solution space to highlight near-optimal strains on Pareto Fronts, balancing the above-cited three objectives. Moreover, preserving the metabolic constraints and the essential genes, this study produced robust results: seven significant strains of Y. lipolytica and seven strains of S. cerevisiae. We examined gene knockout to study the function of genes and pathways. In fact, by studying the frequently silenced genes, we found that when the GPH1 gene is knocked out in S. cerevisiae, the isocitrate lyase enzyme is activated, which converts the isocitrate into succinate. Our goals are to simplify and facilitate the in vitro processes. Hence, we present strains with the least possible number of knockout genes and solutions in which the genes are turned off on the same chromosome. Therefore, we present results where the constraints mentioned above are met, like the strains where only two genes are switched off and other strains where half of the knockout genes are on the same chromosome. This study offers solutions for developing an efficient in vitro mutagenesis for microorganisms and demonstrates the efficiency of multiobjective optimization in automatizing metabolic engineering processes.

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

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

Towards dynamic genome-scale models.

The analysis of the dynamic behaviour of genome-scale models of metabolism (GEMs) currently presents considerable challenges because of the difficulties of simulating such large and complex networks. Bacterial GEMs can comprise about 5000 reactions and metabolites, and encode a huge variety of growth conditions; such models cannot be used without sophisticated tool support. This article is intended to aid modellers, both specialist and non-specialist in computerized methods, to identify and apply a suitable combination of tools for the dynamic behaviour analysis of large-scale metabolic designs. We describe a methodology and related workflow based on publicly available tools to profile and analyse whole-genome-scale biochemical models. We use an efficient approximative stochastic simulation method to overcome problems associated with the dynamic simulation of GEMs. In addition, we apply simulative model checking using temporal logic property libraries, clustering and data analysis, over time series of reaction rates and metabolite concentrations. We extend this to consider the evolution of reaction-oriented properties of subnets over time, including dead subnets and functional subsystems. This enables the generation of abstract views of the behaviour of these models, which can be large-up to whole genome in size-and therefore impractical to analyse informally by eye. We demonstrate our methodology by applying it to a reduced model of the whole-genome metabolism of Escherichia coli K-12 under different growth conditions. The overall context of our work is in the area of model-based design methods for metabolic engineering and synthetic biology.

The view of microbes as energy converters illustrates the trade-off between growth rate and yield.

The application of thermodynamics to microbial growth has a long tradition that originated in the middle of the 20th century. This approach reflects the view that self-replication is a thermodynamic process that is not fundamentally different from mechanical thermodynamics. The key distinction is that a free energy gradient is not converted into mechanical (or any other form of) energy but rather into new biomass. As such, microbes can be viewed as energy converters that convert a part of the energy contained in environmental nutrients into chemical energy that drives self-replication. Before the advent of high-throughput sequencing technologies, only the most central metabolic pathways were known. However, precise measurement techniques allowed for the quantification of exchanged extracellular nutrients and heat of growing microbes with their environment. These data, together with the absence of knowledge of metabolic details, drove the development of so-called black-box models, which only consider the observable interactions of a cell with its environment and neglect all details of how exactly inputs are converted into outputs. Now, genome sequencing and genome-scale metabolic models (GEMs) provide us with unprecedented detail about metabolic processes inside the cell. However, mostly due to computational complexity issues, the derived modelling approaches make surprisingly little use of thermodynamic concepts. Here, we review classical black-box models and modern approaches that integrate thermodynamics into GEMs. We also illustrate how the description of microbial growth as an energy converter can help to understand and quantify the trade-off between microbial growth rate and yield.

Evaluating accessibility, usability and interoperability of genome-scale metabolic models for diverse yeasts species.

Metabolic network reconstructions have become an important tool for probing cellular metabolism in the field of systems biology. They are used as tools for quantitative prediction but also as scaffolds for further knowledge contextualization. The yeast Saccharomyces cerevisiae was one of the first organisms for which a genome-scale metabolic model (GEM) was reconstructed, in 2003, and since then 45 metabolic models have been developed for a wide variety of relevant yeasts species. A systematic evaluation of these models revealed that-despite this long modeling history-the sequential process of tracing model files, setting them up for basic simulation purposes and comparing them across species and even different versions, is still not a generalizable task. These findings call the yeast modeling community to comply to standard practices on model development and sharing in order to make GEMs accessible and useful for a wider public.

Metabolic flux sampling predicts strain-dependent differences related to aroma production among commercial wine yeasts.

BACKGROUND: Metabolomics coupled with genome-scale metabolic modeling approaches have been employed recently to quantitatively analyze the physiological states of various organisms, including Saccharomyces cerevisiae. Although yeast physiology in laboratory strains is well-studied, the metabolic states under industrially relevant scenarios such as winemaking are still not sufficiently understood, especially as there is considerable variation in metabolism between commercial strains. To study the potential causes of strain-dependent variation in the production of volatile compounds during enological conditions, random flux sampling and statistical methods were used, along with experimental extracellular metabolite flux data to characterize the differences in predicted intracellular metabolic states between strains. RESULTS: It was observed that four selected commercial wine yeast strains (Elixir, Opale, R2, and Uvaferm) produced variable amounts of key volatile organic compounds (VOCs). Principal component analysis was performed on extracellular metabolite data from the strains at three time points of cell cultivation (24, 58, and 144 h). Separation of the strains was observed at all three time points. Furthermore, Uvaferm at 24 h, for instance, was most associated with propanol and ethyl hexanoate. R2 was found to be associated with ethyl acetate and Opale could be associated with isobutanol while Elixir was most associated with phenylethanol and phenylethyl acetate. Constraint-based modeling (CBM) was employed using the latest genome-scale metabolic model of yeast (Yeast8) and random flux sampling was performed with experimentally derived fluxes at various stages of growth as constraints for the model. The flux sampling simulations allowed us to characterize intracellular metabolic flux states and illustrate the key parts of metabolism that likely determine the observed strain differences. Flux sampling determined that Uvaferm and Elixir are similar while R2 and Opale exhibited the highest degree of differences in the Ehrlich pathway and carbon metabolism, thereby causing strain-specific variation in VOC production. The model predictions also established the top 20 fluxes that relate to phenotypic strain variation (e.g. at 24 h). These fluxes indicated that Opale had a higher median flux for pyruvate decarboxylase reactions compared with the other strains. Conversely, R2 which was lower in all VOCs, had higher median fluxes going toward central metabolism. For Elixir and Uvaferm, the differences in metabolism were most evident in fluxes pertaining to transaminase and hexokinase associated reactions. The applied analysis of metabolic divergence unveiled strain-specific differences in yeast metabolism linked to fusel alcohol and ester production. CONCLUSIONS: Overall, this approach proved useful in elucidating key reactions in amino acid, carbon, and glycerophospholipid metabolism which suggest genetic divergence in activity in metabolic subsystems among these wine strains related to the observed differences in VOC formation. The findings in this study could steer more focused research endeavors in developing or selecting optimal aroma-producing yeast stains for winemaking and other types of alcoholic fermentations.

Metabolic engineering design to enhance (R,R)-2,3-butanediol production from glycerol in Bacillus subtilis based on flux balance analysis.

BACKGROUND: Glycerol is a desirable alternative substrate for 2,3-butanediol (2,3-BD) production for sustainable development in biotechnological industries and non-food competitive feedstock. B. subtilis, a "generally recognized as safe" organism that is highly tolerant to fermentation products, is an ideal platform microorganism to engineer the pathways for the production of valuable bio-based chemicals, but it has never been engineered to improve 2,3-BD production from glycerol. In this study, we aimed to enhance 2,3-BD production from glycerol in B. subtilis through in silico analysis. Genome-scale metabolic model (GSM) simulations was used to design and develop the metabolic pathways of B. subtilis. Flux balance analysis (FBA) simulation was used to evaluate the effects of step-by-step gene knockouts to improve 2,3-BD production from glycerol in B. subtilis. RESULTS: B. subtilis was bioengineered to enhance 2,3-BD production from glycerol using FBA in a published GSM model of B. subtilis, iYO844. Four genes, ackA, pta, lctE, and mmgA, were knocked out step by step, and the effects thereof on 2,3-BD production were evaluated. While knockout of ackA and pta had no effect on 2,3-BD production, lctE knockout led to a substantial increase in 2,3-BD production. Moreover, 2,3-BD production was improved by mmgA knockout, which had never been investigated. In addition, comparisons between in silico simulations and fermentation profiles of all B. subtilis strains are presented in this study. CONCLUSIONS: The strategy developed in this study, using in silico FBA combined with experimental validation, can be used to optimize metabolic pathways for enhanced 2,3-BD production from glycerol. It is expected to provide a novel platform for the bioengineering of strains to enhance the bioconversion of glycerol into other highly valuable chemical products.