Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “metabolic modeling”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources
PubMed2025

Reducing redundancy and enhancing accuracy through a phylogenetically-informed microbial community metabolic modeling approach.

MOTIVATION: Metabolic modeling has emerged as a powerful tool for predicting community functions. However, current modeling approaches face significant challenges in balancing the metabolic trade-offs between individual and community-level growth. In this study, we investigated the effect of metabolic relatedness among taxa on growth rate calculations by merging related taxa based on their metabolic similarity, introducing this approach as PhyloCOBRA. RESULTS: This approach enhanced the accuracy and efficiency of microbial community simulations by combining genome-scale metabolic models (GEMs) of closely related organisms, aligning with the concepts of niche differentiation and nestedness theory. To validate our approach, we implemented PhyloCOBRA within the MICOM and OptCom package (creating PhyloMICOM and PhyloOptCom, respectively), and applied it to metagenomic data from 186 individuals and four-species synthetic community (SynCom). Our results demonstrated significant improvement in the accuracy and reliability of growth rate predictions compared to the standard methods. Sensitivity analysis revealed that PhyloMICOM models were more robust to random noise, while Jaccard index calculations showed a reduction in redundancy, highlighting the enhanced specificity of the generated community models. Furthermore, PhyloMICOM reduced the computational complexity, addressing a key concern in microbial community simulations. This approach marks a significant advancement in community-scale metabolic modeling, offering a more stable, efficient, and ecologically relevant tool for simulating and understanding the intricate dynamics of microbial ecosystems. AVAILABILITY AND IMPLEMENTATION: PhyloCOBRA implementations are available as extensions to the MICOM packages and can be accessed at https://github.com/sepideh-mofidifar/PhyloCOBRA.

Phylogeny↗
PubMed2026

gMISpy: integration of complex regulatory networks and genome scale metabolic models.

MOTIVATION: Genome-scale metabolic models lack explicit regulatory mechanisms, limiting their predictive accuracy for genetic interventions. Current methods for computing genetic Minimal Cut Sets either ignore regulatory networks entirely or use simplified acyclic representations that cannot capture regulatory feedback loops, ubiquitous features critical in cellular modeling. RESULTS: We developed gMISpy, a Python package that that enables efficient computation of genetic Minimal Intervention Sets (gMISs) in integrated genome-scale metabolic and regulatory networks. gMISpy incorporates cyclic regulatory logic into our previous computational framework using layered Boolean networks and BoNesis framework, resulting in a more accurate modeling of how regulatory interactions affect metabolic genes. Benchmarking across four different regulatory networks with Human-GEM showed consistent improvements in prediction accuracy, with Matthews correlation coefficient gains ranging from 2.50% to 14.42%. Validation against cancer data from DepMap and Project Score confirmed that cyclic integration reduces false positives and better captures biological vulnerabilities compared to acyclic approaches. AVAILABILITY AND IMPLEMENTATION: https://github.com/PlanesLab/cyclic-gMISpy.

Software↗
PubMed1979

Metabolic model for cadmium in man.

A metabolic model for cadmium has been formulated in terms of a flow scheme for cadmium among eight body compartments. The mathematical description of the flow of cadmium between compartments consits of a number of differential equations, and the accumulation of cadmium may thus be calculated. Coeffcients for the flow of cadmium were estimated from empirical data both from animals and man. The modelling serves two main purposes: it provides a means of using present knowledge about metabolism in order to calculate expected accumulation in critical organs and other tissues and body fluids at certain intake levels and it makes it possible to define deficiencies in our present knowledge about metabolism. Data recently collected, partly as a result of considerations related to the model, have confirmed the assumptions of a very long biological half-time in other tissues and the small excretion of cadmium via bile.

Adult↗
PubMed2025

Evolution and applications of genome-scale metabolic models in yeast systems biology studies.

Genome-scale metabolic models (GEMs) can be used to simulate the metabolic network of an organism in a systematic and holistic way. Different yeast species, including Saccharomyces cerevisiae, have emerged as powerful cell factories for bioproduction. Recently, with the dedicated efforts from the scientific community, significant progress has been made in the development of yeast GEMs. Numerous versions of yeast GEMs and the derived multiscale models have been released, facilitating integrative omics analysis and rational strain design for different types of yeast cell factories. These advancements reflected the evolution and maturation of yeast GEMs together with a model ecosystem around them. This review will summarize the development and expansion of yeast GEMs and discuss their applications in yeast systems biology studies. It is anticipated that yeast GEMs will continue to play an increasingly important role in pioneering yeast physiological and metabolic studies in coming years.

Systems Biology↗
PubMed2026

WILDkCAT: extract, retrieve, and predict enzyme turnover numbers of constraint-based metabolic models.

SUMMARY: Accurate enzyme turnover numbers are essential for building enzyme-constrained genome-scale metabolic models. However, collecting and curating these parameters remains a major bottleneck. Indeed, kcat values are scattered across multiple databases, reported under varying experimental conditions, and often missing for many enzymes. To address this challenge, we present WILDkCAT, a Python-based pipeline that enables the retrieval of kcat values from wild-type enzyme measured under user-specified pH and temperature ranges for a given metabolic model. The application to Escherichia coli (iML1515) and Homo sapiens (Human-GEM) models demonstrated the ability of WILDkCAT to retrieve substantial kcat coverage and its applicability across diverse genome-scale models. AVAILABILITY AND IMPLEMENTATION: WILDkCAT is available at https://github.com/sysbiolux/WILDkCAT and from PyPI. WILDkCAT works on all major operating systems and computer architectures. The documentation is available at https://sysbiolux.github.io/WILDkCAT.

Software↗
PubMed2025

A genome-scale metabolic model for the denitrifying bacterium Thauera sp. MZ1T accurately predicts degradation of pollutants and production of polymers.

The denitrifying bacterium Thauera sp. MZ1T, a common member of microbial communities in wastewater treatment facilities, can produce different compounds from a range of carbon (C) and nitrogen (N) sources under aerobic and anaerobic conditions. In these different conditions, Thauera modifies its metabolism to produce different compounds that influence the microbial community. In particular, Thauera sp. MZ1T produces different exopolysaccharides with floc-forming properties, impacting the physical disposition of wastewater consortia and the efficiency of nutrient assimilation by the microbial community. Under N-limiting conditions, Thauera sp. MZ1T decreases its growth rate and accelerates the accumulation of polyhydroxyalkanoate-related (PHA) compounds including polyhydroxybutyrate (PHB), which plays a fundamental role as C and energy storage in this β-proteobacterium. However, the metabolic mechanisms employed by Thauera sp. MZ1T to assimilate and catabolize many of the different C and N sources under aerobic and anaerobic conditions remain unknown. Systems biology approaches such as genome-scale metabolic modeling have been successfully used to unveil complex metabolic mechanisms for various microorganisms. Here, we developed a comprehensive metabolic model (M-model) for Thauera sp. MZ1T (iThauera861), consisting of 1,744 metabolites, 2,384 reactions, and 861 genes. We validated the model experimentally using over 70 different C and N sources under both aerobic and anaerobic conditions. iThauera861 achieved a prediction accuracy of 95% for growth on various C and N sources and close to 85% for assimilation of aromatic compounds under denitrifying conditions. The M-model was subsequently deployed to determine the effects of substrates, oxygen presence, and the C:N ratio on the production of PHB and exopolysaccharides (EPS), showing the highest polymer yields are achieved with nucleotides and amino acids under aerobic conditions. This comprehensive M-model will help reveal the metabolic processes by which this ubiquitous species influences communities in wastewater treatment systems and natural environments.

Thauera↗
PubMed2025

In silico encounters: harnessing metabolic modelling to understand plant-microbe interactions.

Understanding plant-microbe interactions is vital for developing sustainable agricultural practices and mitigating the consequences of climate change on food security. Plant-microbe interactions can improve nutrient acquisition, reduce dependency on chemical fertilizers, affect plant health, growth, and yield, and impact plants' resistance to biotic and abiotic stresses. These interactions are largely driven by metabolic exchanges and can thus be understood through metabolic network modelling. Recent developments in genomics, metagenomics, phenotyping, and synthetic biology now enable researchers to harness the potential of metabolic modelling at the genome scale. Here, we review studies that utilize genome-scale metabolic modelling to study plant-microbe interactions in symbiotic, pathogenic, and microbial community systems. This review catalogues how metabolic modelling has advanced our understanding of the plant host and its associated microorganisms as a holobiont. We showcase how these models can contextualize heterogeneous datasets and serve as valuable tools to dissect and quantify underlying mechanisms. Finally, we consider studies that employ metabolic models as a testbed for in silico design of synthetic microbial communities with predefined traits. We conclude by discussing broader implications of the presented studies, future perspectives, and outstanding challenges.

Plants↗
PubMed1978

Cadmium concentrations in human liver, blood, and bile: comparison with a metabolic model.

Cadmium concentrations in liver biopsies, blood, and bile were measured by atomic absorption spectrophotometry in 23 patients in connection with routine gallstone operations. On a group basis cadmium in blood was a good indicator of cadmium in liver, and the estimated linear relationship agreed well with calculations from a formerly proposed metabolic model. Cadmium in bile was also analyzed, and an average of about 2.5 ng of Cd/g wet weight was found. This is about 10 times more than would have been expected form the metabolic model and suggests that bile might be an important excretion route for cadmium. Definite conclusions cannot be drawn, however, since the results could not be cross-checked with neutron activation analysis, due to insufficient sensitivity of the latter method.

Adult↗
PubMed2025

HUMESS: integrating quantitative transcriptomic analysis and metabolic modeling to unveil condition-specific gene signatures.

SUMMARY: Transcriptomic analysis is a key tool for exploring gene expression, but the complexity of biological systems often limits its insights. In particular, the lack of intermodal or multi-layered analysis hinders the ability to fully capture key cellular functions such as metabolism from transcriptomic data alone. Here, we introduce a novel approach that informs transcriptomic data analysis with metabolic network modeling to address this. Unlike traditional methods, HUman MEtabolism Specific Signature (HUMESS) uses genome-scale metabolic modeling and flux analysis to highlight reactions and involved genes based on their metabolic significance, offering a deeper understanding of transcriptomic data. Our computational pipeline, supported by a user-friendly Rshiny application, enhances gene expression analysis by uncovering metabolic phenotypic signatures. AVAILABILITY AND IMPLEMENTATION: HUMESS is open source and available under GitLab https://gitlab.univ-nantes.fr/bird_pipeline_registry/humess with the complete documentation available at https://gitlab.univ-nantes.fr/bird_pipeline_registry/humess/-/wikis/Home. A zenodo archive is also available at the following DOI: https://doi.org/10.5281/zenodo.15487717. An RShiny application has been developed to facilitate the exploration and analysis of HUMESS's results. The app is available online at the following address: https://shiny-bird.univ-nantes.fr/app/shinymess but can also be installed locally, available under GitLab https://gitlab.univ-nantes.fr/pare-l/shinymess.

Humans↗
PubMed2025

RBC-GEM: A genome-scale metabolic model for systems biology of the human red blood cell.

Advancements with cost-effective, high-throughput omics technologies have had a transformative effect on both fundamental and translational research in the medical sciences. These advancements have facilitated a departure from the traditional view of human red blood cells (RBCs) as mere carriers of hemoglobin, devoid of significant biological complexity. Over the past decade, proteomic analyses have identified a growing number of different proteins present within RBCs, enabling systems biology analysis of their physiological functions. Here, we introduce RBC-GEM, one of the most comprehensive, curated genome-scale metabolic reconstructions of a specific human cell type to-date. It was developed through meta-analysis of proteomic data from 29 studies published over the past two decades resulting in an RBC proteome composed of more than 4,600 distinct proteins. Through workflow-guided manual curation, we have compiled the metabolic reactions carried out by this proteome to form a genome-scale metabolic model (GEM) of the RBC. RBC-GEM is hosted on a version-controlled GitHub repository, ensuring adherence to the standardized protocols for metabolic reconstruction quality control and data stewardship principles. RBC-GEM represents a metabolic network is a consisting of 820 genes encoding proteins acting on 1,685 unique metabolites through 2,723 biochemical reactions: a 740% size expansion over its predecessor. We demonstrated the utility of RBC-GEM by creating context-specific proteome-constrained models derived from proteomic data of stored RBCs for 616 blood donors, and classified reactions based on their simulated abundance dependence. This reconstruction as an up-to-date curated GEM can be used for contextualization of data and for the construction of a computational whole-cell models of the human RBC.

Humans↗
PubMed2026

High-Quality Genome Assembly, Metabolome, Pangenome, and Metabolic Models of Megasphaera hexanoica KCCM 43214T.

Megasphaera hexanoica KCCM 43214T, isolated from cow rumen, is capable of producing medium-chain carboxylic acids such as hexanoate and octanoate. In this study, we present a high-quality genome assembly, along with intracellular metabolomic profiling and pangenomic analysis. Illumina sequencing generated 2.3 Gbp from 15,293,634 reads with a GC content of 49.5%, while PacBio HiFi sequencing produced 331.5 Mbp across 45,266 reads, with an average read length of 7,323 bp and a HiFi read N50 of 8,214 bp. Hybrid assembly of short and long reads resulted in a single 2.88 Mbp contig, containing 2,835 protein-coding genes. Genome-scale metabolic models were constructed to evaluate its metabolic capabilities under specific growth conditions. Intracellular metabolomic analysis of cells grown in medium containing fructose and lactate revealed key metabolic activities associated with chain elongation. Pangenomic analysis across nine annotated genomes identified 6,721 orthologous genes using OrthoMCL, emphasizing the genetic and functional diversity within the Megasphaera genus. This dataset offers valuable insights into the metabolism and biotechnological potential of M. hexanoica KCCM 43214T.

Metabolome↗
PubMed2026

hypeR-GEM: connecting metabolite signatures to enzyme-coding genes via genome-scale metabolic models.

MOTIVATION: Enrichment analysis is a cornerstone of "omics" data interpretation, enabling researchers to connect analysis results to biological processes and generate testable hypotheses. Enrichment analysis in metabolomics poses distinct challenges for interpretation and multi-omics integration due to the lack of well-defined and consistent connections to well-curated gene-centered biological knowledge repositories. To address these challenges, we developed hypeR-GEM, a methodology and associated R package that adapts gene set enrichment analysis to metabolomics. hypeR-GEM leverages genome-scale metabolic models (GEMs) to infer reaction-based links between metabolites and enzyme-coding genes, enabling the mapping of metabolite signatures to gene signatures and their subsequent annotation via gene set enrichment analysis. RESULTS: We validated hypeR-GEM using paired metabolomics-proteomics and metabolomics-transcriptomics datasets by assessing whether genes mapped from metabolites significantly overlapped with differentially expressed proteins or transcripts. We further evaluated whether pathways enriched via hypeR-GEM-mapped genes corresponded to those derived from paired proteomic or transcriptomic data. In most datasets analyzed, both the predicted enzyme-coding genes and the associated enriched pathways showed significant concordance with independently derived omics signatures, supporting the utility and robustness of hypeR-GEM. Finally, we applied hypeR-GEM to the analysis of age-associated metabolic signatures from the New England Centenarian Study. The results revealed consistent enrichment of lipid-related pathways, aligning with the well-established role of lipid metabolism in aging, and highlighted additional pathways not captured in the metabolites' annotation, demonstrating hypeR-GEM's practical utility in a real-world use case. AVAILABILITY AND IMPLEMENTATION: The hypeR-GEM R package, documentation, and workflow examples are freely available at https://github.com/montilab/hypeR-GEM and archived at https://doi.org/10.5281/zenodo.20586748.

Metabolomics↗
PubMed2025

Biophysical metabolic modeling of complex bacterial colony morphology.

Microbial colony growth is shaped by the physics of biomass propagation and nutrient diffusion and by the metabolic reactions that organisms activate as a function of the surrounding environment. While microbial colonies have been explored using minimal models of growth and motility, full integration of biomass propagation and metabolism is still lacking. Here, building upon our framework for computation of microbial ecosystems in time and space (COMETS), we combine dynamic flux balance modeling of metabolism with collective biomass propagation and demographic fluctuations to provide nuanced simulations of E. coli colonies. Simulations produced realistic colony morphology, consistent with our experiments. They characterize the transition between smooth and furcated colonies and the decay of genetic diversity. Furthermore, we demonstrate that under certain conditions, biomass can accumulate along "metabolic rings" that are reminiscent of coffee-stain rings but have a completely different origin. Our approach is a key step toward predictive microbial ecosystems modeling. A record of this paper's transparent peer review process is included in the supplemental information.

Models, Biological↗
PubMed2026

Dynamic metabolic modelling of ATP allocation during viral infection.

Viral pathogens, like SARS-CoV-2, hijack the host's macromolecular production machinery, imposing an energetic burden that is distributed across cellular metabolism. To explore the dynamic metabolic tension between the host's survival and viral replication, we developed a computational framework that uses genome-scale models to perform dynamic flux balance analysis of human cell metabolism during virus infections. Relative to previous models, our framework addresses the physiology of viral infections of non-proliferating host cells through two new features. First, by incorporating the lipid content of SARS-CoV-2 biomass, we discovered activation of previously overlooked pathways giving rise to new predictions of possible drug targets. Furthermore, we introduce a dynamic model that simulates the partitioning of resources between the virus and the host cell, capturing the extent to which the competition depletes the human cells from essential ATP. By incorporating viral dynamics into our COMETS framework for spatio-temporal modelling of metabolism, we provide a mechanistic, dynamic and generalizable starting point for bridging systems biology modelling with viral pathogenesis. This framework could be extended to broadly incorporate phage dynamics in microbial systems and ecosystems.

Humans↗
PubMed2022

Genome-scale metabolic modeling reveals increased reliance on valine catabolism in clinical isolates of Klebsiella pneumoniae.

Infections due to carbapenem-resistant Enterobacteriaceae have recently emerged as one of the most urgent threats to hospitalized patients within the United States and Europe. By far the most common etiological agent of these infections is Klebsiella pneumoniae, frequently manifesting in hospital-acquired pneumonia with a mortality rate of ~50% even with antimicrobial intervention. We performed transcriptomic analysis of data collected previously from in vitro characterization of both laboratory and clinical isolates which revealed shifts in expression of multiple master metabolic regulators across isolate types. Metabolism has been previously shown to be an effective target for antibacterial therapy, and genome-scale metabolic network reconstructions (GENREs) have provided a powerful means to accelerate identification of potential targets in silico. Combining these techniques with the transcriptome meta-analysis, we generated context-specific models of metabolism utilizing a well-curated GENRE of K. pneumoniae (iYL1228) to identify novel therapeutic targets. Functional metabolic analyses revealed that both composition and metabolic activity of clinical isolate-associated context-specific models significantly differs from laboratory isolate-associated models of the bacterium. Additionally, we identified increased catabolism of L-valine in clinical isolate-specific growth simulations. These findings warrant future studies for potential efficacy of valine transaminase inhibition as a target against K. pneumoniae infection.

Humans↗
PubMed1977

[Experimental studies on clinical application of 13N-ammonia; a metabolic model and functional image based on animal experiment].

Ammonia has been known to be an inducing agent of hepatic coma associated with liver cirrhosis. In order to establish a non-invasive method of investigating the portal circulation and the metabolism of ammonia at liver the following studies were performed, by using cyclotron produced 13N-labeled ammonia and delayed line-camera interphased with on-line computer system. 1) Animal experiment. Dynamic scintigraphy of thoraco-abdominal region of a rabbit after intrasigmoidal as well as intravenous administration of 13N-ammonia were performed. 2) Theoretical consideration on the functional image of liver: According to the results of the animal experiments, a model of ammonia metabolism having the minimal complexity necessary to represent data in the case of intrasigmoidal administration was proposed. Starting from the solution to the differential equation describing the model, relevant parameters characterizing the dynamic curve of liver were duded.

Ammonia↗
PubMed2026

Patient-specific modeling identifies metabolic interventions for reversing glucose use reprogramming in alcohol-associated hepatitis.

Alcoholic hepatitis (AH) is an acute form of alcohol-associated liver disease with very few treatment options. Recent studies highlighted liver metabolic reprogramming in AH as an indicator of severity. We aim at identifying new intervention points to reverse liver metabolic dysregulation across varying degrees of AH. We develop 89 personalized genome-scale metabolic models by integrating a generic human cellular metabolic model with liver transcriptomics data from AH patients with varying disease severity and healthy controls. We grade the AH patients based on the model-predicted level of glycolysis reprogramming and validate the results using published metabolomics data. We test in silico gene knockdown interventions to reverse the aberrant metabolic reprogramming in AH. Knockdown of two glycolytic genes, Hkdc1 and Pkm, significantly rebalance the metabolic fluxes toward a healthy liver metabolic phenotype. We use machine learning on the glycolysis fluxes to develop a quantitative glucose use reprogramming score, which correlates with AH severity and patient-specific responses to in silico gene knockdown interventions. The score was independently validated using a published AH liver transcriptomics dataset. We propose a cellular metabolism-based therapy targeting Hkdc1 and Pkm in the glycolysis pathway as a potential treatment for reversing the aberrant glucose metabolism in AH.

Humans↗