XomicsToModel: omics data integration and generation of thermodynamically consistent metabolic models.

Constraint-based models (CBMs) are used to study metabolic network structure and function in organisms ranging from microbes to multicellular eukaryotes. Published CBMs are usually generic rather than context-specific, meaning that they do not capture differences in reaction activities, which, in turn, determine metabolic capabilities, between cell types, tissues, environments, or other conditions. Only a subset of a CBM’s metabolic reactions and capabilities are likely to be active in any given context, and several methods have therefore been developed to extract context-specific models from generic CBMs through integration of omics data. We tested the ability of six model extraction methods (MEMs) to create functionally accurate context-specific models of Atlantic salmon using a generic CBM (SALARECON) and liver transcriptomics data from contexts differing in water salinity (life stage) and dietary lipids. Three MEMs (iMAT, INIT, and GIMME) outperformed the others in terms of functional accuracy, which we defined as the extracted models’ ability to perform context-specific metabolic tasks inferred directly from the data, and one MEM (GIMME) was faster than the others. Context-specific versions of SALARECON consistently outperformed the generic version, showing that context-specific modeling better captures salmon metabolism. Thus, we demonstrate that results from human studies also hold for a non-mammalian animal and major livestock species.

Genome-scale metabolic network model (GEM) is a fundamental framework in systems metabolic engineering. GEM is built upon extensive experimental data and literature information on gene annotation and function, metabolites and enzymes so that it contains all known metabolic reactions within an organism. Constraint-based analysis of GEM enables the identification of phenotypic properties of an organism and hypothesis-driven engineering of cellular functions to achieve objectives. Along with the advances in omics, high-throughput technology and computational algorithms, the scope and applications of GEM have substantially expanded. In particular, various computational algorithms have been developed to predict beneficial gene deletion and amplification targets and used to guide the strain development process for the efficient production of industrially important chemicals. Furthermore, an Escherichia coli GEM was integrated with a pathway prediction algorithm and used to evaluate all possible routes for the production of a list of commodity chemicals in E. coli. Combined with the wealth of experimental data produced by high-throughput techniques, much effort has been exerted to add more biological contexts into GEM through the integration of omics data and regulatory network information for the mechanistic understanding and improved prediction capabilities. In this paper, we review the recent developments and applications of GEM focusing on the GEM-based computational algorithms available for microbial metabolic engineering.

Genome-scale metabolic network reconstructions are considered a key step in quantifying the genotype-phenotype relationship. We present a novel gap-filling approach, MetabolIc Reconstruction via functionAl GEnomics (MIRAGE), which identifies missing network reactions by integrating metabolic flux analysis and functional genomics data. MIRAGE's performance is demonstrated on the reconstruction of metabolic network models of E. coli and Synechocystis sp. and validated via existing networks for these species. Then, it is applied to reconstruct genome-scale metabolic network models for 36 sequenced cyanobacteria amenable for constraint-based modeling analysis and specifically for metabolic engineering. The reconstructed network models are supplied via standard SBML files.

Cheese taste and flavour properties result from complex metabolic processes occurring in microbial communities. A deeper understanding of such mechanisms makes it possible to improve both industrial production processes and end-product quality through the design of microbial consortia. In this work, we caracterise the metabolism of a three-species community consisting of Lactococcus lactis, Lactobacillus plantarum and Propionibacterium freudenreichii during a seven-week cheese production process. Using genome-scale metabolic models and omics data integration, we modeled and calibrated individual dynamics using monoculture experiments, and coupled these models to capture the metabolism of the community. This model accurately predicts the dynamics of the community, enlightening the contribution of each microbial species to organoleptic compound production. Further metabolic exploration revealed additional possible interactions between the bacterial species. This work provides a methodological framework for the prediction of community-wide metabolism and highlights the added value of dynamic metabolic modeling for the comprehension of fermented food processes.

BackgroundClostridium difficile is a bacterium which can infect various animal species, including humans. Infection with this bacterium is a leading healthcare-associated illness. A better understanding of this organism and the relationship between its genotype and phenotype is essential to the search for an effective treatment. Genome-scale metabolic models contain all known biochemical reactions of a microorganism and can be used to investigate this relationship.ResultsWe present icdf834, an updated metabolic network of C. difficile that builds on iMLTC806cdf and features 1227 reactions, 834 genes, and 807 metabolites. We used this metabolic network to reconstruct the metabolic landscape of this bacterium. The standard metabolic model cannot account for changes in the bacterial metabolism in response to different environmental conditions. To account for this limitation, we also integrated transcriptomic data, which details the gene expression of the bacterium in a wide array of environments. Importantly, to bridge the gap between gene expression levels and protein abundance, we accounted for the synonymous codon usage bias of the bacterium in the model. To our knowledge, this is the first time codon usage has been quantified and integrated into a metabolic model. The metabolic fluxes were defined as a function of protein abundance. To determine potential therapeutic targets using the model, we conducted gene essentiality and metabolic pathway sensitivity analyses and calculated flux control coefficients. We obtained 92.3% accuracy in predicting gene essentiality when compared to experimental data for C. difficile R20291 (ribotype 027) homologs. We validated our context-specific metabolic models using sensitivity and robustness analyses and compared model predictions with literature on C. difficile. The model predicts interesting facets of the bacterium’s metabolism, such as changes in the bacterium’s growth in response to different environmental conditions.ConclusionsAfter an extensive validation process, we used icdf834 to obtain state-of-the-art predictions of therapeutic targets for C. difficile. We show how context-specific metabolic models augmented with codon usage information can be a beneficial resource for better understanding C. difficile and for identifying novel therapeutic targets. We remark that our approach can be applied to investigate and treat against other pathogens.

The advent of model-enabled workflows in systems biology allows for the integration of experimental data types with genome-scale models to discover new features of biology. This work demonstrates such a workflow, aimed at establishing a metabolomics platform applied to study the differences in metabolomes between anaerobic and aerobic growth of Escherichia coli. Constraint-based modeling was utilized to deduce a target list of compounds for downstream method development. An analytical and experimental methodology was developed and tailored to the compound chemistry and growth conditions of interest. This included the construction of a rapid sampling apparatus for use with anaerobic cultures. The resulting genome-scale data sets for anaerobic and aerobic growth were validated by comparison to previous small-scale studies comparing growth of E. coli under the same conditions. The metabolomics data were then integrated with the E. coli genome-scale metabolic model (GEM) via a sensitivity analysis that utilized reaction thermodynamics to reconcile simulated growth rates and reaction directionalities. This analysis highlighted several optimal network usage inconsistencies, including the incorrect use of the beta-oxidation pathway for synthesis of fatty acids. This analysis also identified enzyme promiscuity for the pykA gene, that is critical for anaerobic growth, and which has not been previously incorporated into metabolic models of E coli.

Due to tumour heterogeneity, cellular metabolic activities and tumour proliferation rates can vary across patients, limiting the predictive and prognostic power of population-based metabolic biomarkers. Genome-scale metabolic models have been developed to address this challenge. These models can simulate patient-specific metabolic reaction rates and mechanistically elucidate cellular metabolic alterations. Simultaneously, phenotypes and other omics data can be integrated to provide valuable insights into specific tumour behaviour and molecular mechanisms underlying cancer biology. However, integrative approaches that combine multiple data modalities with metabolic modelling remain largely undeveloped. We propose a unified framework based on interpretable multi-modal machine learning, integrating different data modalities, including transcriptomics, clinical data, and genome-scale metabolic modelling, to stratify breast cancer patients. By integrating validation analysis across data scales, from bulk to single-cell and spatial data, our framework combines mechanistic knowledge from metabolic modelling with machine learning to characterise molecular and metabolic dysregulations across breast cancer risk groups. This unique mechanistic interpretation approach, combining data-driven and biologically-driven knowledge, provides a comprehensive understanding of breast cancer biology to improve personalised cancer therapy and reveal key biomarkers across multiple scales, including patient-specific, cell-specific, and spot-specific prognostic signatures.

A genome-scale metabolic model (GEM) represents metabolic pathways of an organism in a mathematical form and can be built using biochemistry and genome annotation data. GEMs are invaluable for understanding organisms since they analyze the metabolic capabilities and behaviors quantitatively and can predict phenotypes. The development of high-throughput data collection techniques led to an immense increase in omics data such as metagenomics, which expand our knowledge on the human microbiome, but this also created a need for systematic analysis of these data. In recent years, GEMs have also been reconstructed for microbial species, including human gut microbiota, and methods for the analysis of microbial communities have been developed to examine the interaction between the organisms or the host. The purpose of this review is to provide a comprehensive guide for the applications of GEMs in microbial community analysis. Starting with GEM repositories, automatic GEM reconstruction tools, and quality control of models, this review will give insights into microbe-microbe and microbe-host interaction predictions and optimization of microbial community models. Recent studies that utilize microbial GEMs and personalized models to infer the influence of microbiota on human diseases such as inflammatory bowel diseases (IBD) or Parkinson's disease are exemplified. Being powerful system biology tools for both species-level and community-level analysis of microbes, GEMs integrated with omics data and machine learning techniques will be indispensable for studying the microbiome and their effects on human physiology as well as for deciphering the mechanisms behind human diseases.

For many years, antibiotics reliably protected mankind against bacterial infections, including the respiratory tract colonizer Corynebacterium tuberculostearicum. However, the spread of antimicrobial resistance necessitates the search for new treatment options, where the microbiota may play a crucial role. One way to investigate the complex nature of bacteria and their interactions with human hosts or microbiota is through genome-scale metabolic models (GEMs). Constructing GEMs is labor-intensive and time-consuming. We introduce the Python package Mass and Charge Curation (MCC), which implements a new automated algorithm to facilitate mass and charge balancing-one of the most time-consuming reconstruction steps. This package manipulates reconstructions by consolidating data from multiple resources and updating the notes field with relevant changes. It also visually compares draft and curated models, ensuring high-quality metabolic reconstructions. Using MCC, we developed a metabolic reconstruction of C. tuberculostearicum strain DSM 44922. The model was improved based on standardization policies, resulting in a functional, well-annotated, high-quality product. We also simulated the organism's growth in synthetic nasal medium 3 (SNM3). The high-quality model iCTUB2024RM consistently resembles growth behavior under realistic conditions in an artificial human nasal environment, enhancing the understanding of C. tuberculostearicum and its potential impact on health and disease. The curation process of this model led to the development of the MCC package, which facilitates the mass and charge balancing of arbitrary flux balance constraints models in the SBML format.IMPORTANCEThe rise of antibiotic resistance has made it essential to explore alternative treatments for bacterial infections, particularly those caused by respiratory tract colonizers like Corynebacterium tuberculostearicum. Understanding the metabolic behavior of these bacteria and their interactions with the human host or microbiota is crucial. Genome-scale metabolic models (GEMs) are powerful tools for investigating these interactions, but they are time-consuming to build. Our new Python package, Mass and Charge Curation, automates a crucial step in the GEM reconstruction process-mass and charge balancing-making it more efficient and reliable. By applying this tool, we developed a high-quality, functional metabolic model for C. tuberculostearicum (iCTUB2024RM), which provides deeper insights into the organism's growth in a simulated human nasal environment. This work offers a foundation for future research into microbial communities and their role in human health.

Vitis vinifera, also known as grapevine, is widely cultivated and commercialized, particularly to produce wine. As wine quality is directly linked to fruit quality, studying grapevine metabolism is important to understand the processes underlying grape composition. Genome-scale metabolic models (GSMMs) have been used for the study of plant metabolism and advances have been made, allowing the integration of omics datasets with GSMMs. On the other hand, Machine learning (ML) has been used to analyze and integrate omics data, and while the combination of ML with GSMMs has shown promising results, it is still scarcely used to study plants. Here, the first GSSM of V. vinifera was reconstructed and validated, comprising 7199 genes, 5399 reactions, and 5141 metabolites across 8 compartments. Tissue-specific models for the stem, leaf, and berry of the Cabernet Sauvignon cultivar were generated from the original model, through the integration of RNA-Seq data. These models have been merged into diel multi-tissue models to study the interactions between tissues at light and dark phases. The potential of combining ML with GSMMs was explored by using ML to analyze the fluxomics data generated by green and mature grape GSMMs and provide insights regarding the metabolism of grapes at different developmental stages. Therefore, the models developed in this work are useful tools to explore different aspects of grapevine metabolism and understand the factors influencing grape quality.

The liver is one of the largest organs of the adult body and most of its tissue is formed by hepatocyte cells, the main site of the metabolic conversions underlying its diverse physiological functions. Hepatocellular carcinoma is one of the most important human cancers. Genome-Scale Metabolic Models (GSMMs), mathematical representations of the cell metabolism in different organisms including humans, are useful tools to simulate metabolic phenotypes and understand metabolic diseases. In the last years, a few algorithms have been developed to generate tissue-specific metabolic models that allow the simulation of phenotypes for distinct cell types/tissues. This work based on general template GSMMs, which are integrated with available omics data. In this work, we propose to develop a pipeline for the systematic evaluation of these algorithms in the creation of models for regular hepatocytes and cancer cell lines, addressing the comparison of the final models obtained.

Metabolic modeling provides the means to define metabolic processes at a systems level; however, genome-scale metabolic models often remain incomplete in their description of metabolic networks and may include reactions that are experimentally unverified. This shortcoming is exacerbated in reconstructed models of newly isolated algal species, as there may be little to no biochemical evidence available for the metabolism of such isolates. The phenotype microarray (PM) technology (Biolog, Hayward, CA, USA) provides an efficient, high-throughput method to functionally define cellular metabolic activities in response to a large array of entry metabolites. The platform can experimentally verify many of the unverified reactions in a network model as well as identify missing or new reactions in the reconstructed metabolic model. The PM technology has been used for metabolic phenotyping of non-photosynthetic bacteria and fungi, but it has not been reported for the phenotyping of microalgae. Here, we introduce the use of PM assays in a systematic way to the study of microalgae, applying it specifically to the green microalgal model species Chlamydomonas reinhardtii. The results obtained in this study validate a number of existing annotated metabolic reactions and identify a number of novel and unexpected metabolites. The obtained information was used to expand and refine the existing COBRA-based C. reinhardtii metabolic network model iRC1080. Over 254 reactions were added to the network, and the effects of these additions on flux distribution within the network are described. The novel reactions include the support of metabolism by a number of d-amino acids, l-dipeptides, and l-tripeptides as nitrogen sources, as well as support of cellular respiration by cysteamine-S-phosphate as a phosphorus source. The protocol developed here can be used as a foundation to functionally profile other microalgae such as known microalgae mutants and novel isolates.

Genome-scale metabolic models (GEMs) are widely used in systems biology to investigate metabolism and predict perturbation responses. Automatic GEM reconstruction tools generate GEMs with different properties and predictive capacities for the same organism. Since different models can excel at different tasks, combining them can increase metabolic network certainty and enhance model performance. Here, we introduce GEMsembler, a Python package designed to compare cross-tool GEMs, track the origin of model features, and build consensus models containing any subset of the input models. GEMsembler provides comprehensive analysis functionality, including identification and visualization of biosynthesis pathways, growth assessment, and an agreement-based curation workflow. GEMsembler-curated consensus models built from four Lactiplantibacillus plantarum and Escherichia coli automatically reconstructed models outperform the gold-standard models in auxotrophy and gene essentiality predictions. Optimizing gene-protein-reaction (GPR) combinations from consensus models improves gene essentiality predictions, even in the manually curated gold-standard models. GEMsembler explains model performance by highlighting relevant metabolic pathways and GPR alternatives, informing experiments to resolve model uncertainty. Thus, GEMsembler facilitates building more accurate and biologically informed metabolic models for systems biology applications.IMPORTANCEGenome-scale metabolic models (GEMs) capture our knowledge of cellular metabolism as encoded in the genome, enabling us to describe and predict how cells function under different conditions. While several automated tools can generate these models directly from genome data, the resulting models often contain gaps and uncertainties, highlighting areas where our metabolic knowledge is incomplete. Here, we introduce a new tool called GEMsembler, which integrates GEMs constructed by different methods, evaluate model uncertainty, and build consensus models, harnessing the unique features of each approach. These consensus models more accurately reflect experimentally observed metabolic traits, such as nutrient requirements and condition-specific gene essentiality. GEMsembler facilitates comprehensive analysis of model structure and function, helping to pinpoint knowledge gaps and prioritize experiments to address them. By synthesizing information from diverse sources, GEMsembler accelerates the development of more reliable and biologically meaningful models, advancing research in metabolic engineering, pathogen biology, and microbial community studies.

Metabolically re-modeling the drug pipeline

Bacillus subtilis is a well-characterized microorganism and a model for the study of Gram-positive bacteria. The bacterium can produce proteins at high densities and yields, which has made it valuable for industrial bioproduction. Like other cell factories, metabolic modeling of B. subtilis has discovered ways to optimize its metabolism toward various applications. The first genome-scale metabolic model (M-model) of B. subtilis was published more than a decade ago and has been applied extensively to understand metabolism, to predict growth phenotypes, and served as a template to reconstruct models for other Gram-positive bacteria. However, M-models are ill-suited to simulate the production and secretion of proteins as well as their proteomic response to stress. Thus, a new generation of metabolic models, known as metabolism and gene expression models (ME-models), has been initiated. Here, we describe the reconstruction and validation of a ME model of B. subtilis, iJT964-ME. This model achieved higher performance scores on the prediction of gene essentiality as compared to the M-model. We successfully validated the model by integrating physiological and omics data associated with gene expression responses to ethanol and salt stress. The model further identified the mechanism by which tryptophan synthesis is upregulated under ethanol stress. Further, we employed iJT964-ME to predict amylase production rates under two different growth conditions. We analyzed these flux distributions and identified key metabolic pathways that permitted the increase in amylase production. Models like iJT964-ME enable the study of proteomic response to stress and the illustrate the potential for optimizing protein production in bacteria.