SEARCH · PubMed Health

Results for “Genome scale models”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 37 records · Page 2Linked to original sources

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↗

PubMed2006

Identification of genome-scale metabolic network models using experimentally measured flux profiles.

Genome-scale metabolic network models can be reconstructed for well-characterized organisms using genomic annotation and literature information. However, there are many instances in which model predictions of metabolic fluxes are not entirely consistent with experimental data, indicating that the reactions in the model do not match the active reactions in the in vivo system. We introduce a method for determining the active reactions in a genome-scale metabolic network based on a limited number of experimentally measured fluxes. This method, called optimal metabolic network identification (OMNI), allows efficient identification of the set of reactions that results in the best agreement between in silico predicted and experimentally measured flux distributions. We applied the method to intracellular flux data for evolved Escherichia coli mutant strains with lower than predicted growth rates in order to identify reactions that act as flux bottlenecks in these strains. The expression of the genes corresponding to these bottleneck reactions was often found to be downregulated in the evolved strains relative to the wild-type strain. We also demonstrate the ability of the OMNI method to diagnose problems in E. coli strains engineered for metabolite overproduction that have not reached their predicted production potential. The OMNI method applied to flux data for evolved strains can be used to provide insights into mechanisms that limit the ability of microbial strains to evolve towards their predicted optimal growth phenotypes. When applied to industrial production strains, the OMNI method can also be used to suggest metabolic engineering strategies to improve byproduct secretion. In addition to these applications, the method should prove to be useful in general for reconstructing metabolic networks of ill-characterized microbial organisms based on limited amounts of experimental data.

Computational Biology↗

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↗

PubMed2025

A genome-scale metabolic reconstruction resource of 247,092 diverse human microbes spanning multiple continents, age groups, and body sites.

Genome-scale modeling of microbiome metabolism enables the simulation of diet-host-microbiome-disease interactions. However, current genome-scale reconstruction resources are limited in scope by computational challenges. We developed an optimized and highly parallelized reconstruction and analysis pipeline to build a resource of 247,092 microbial genome-scale metabolic reconstructions, deemed APOLLO. APOLLO spans 19 phyla, contains >60% of uncharacterized strains, and accounts for strains from 34 countries, all age groups, and multiple body sites. Using machine learning, we predicted with high accuracy the taxonomic assignment of strains based on the computed metabolic features. We then built 14,451 metagenomic sample-specific microbiome community models to systematically interrogate their community-level metabolic capabilities. We show that sample-specific metabolic pathways accurately stratify microbiomes by body site, age, and disease state. APOLLO is freely available, enables the systematic interrogation of the metabolic capabilities of largely still uncultured and unclassified species, and provides unprecedented opportunities for systems-level modeling of personalized host-microbiome co-metabolism.

Humans↗

PubMed2005

Large-scale 13C-flux analysis reveals mechanistic principles of metabolic network robustness to null mutations in yeast.

BACKGROUND: Quantification of intracellular metabolite fluxes by 13C-tracer experiments is maturing into a routine higher-throughput analysis. The question now arises as to which mutants should be analyzed. Here we identify key experiments in a systems biology approach with a genome-scale model of Saccharomyces cerevisiae metabolism, thereby reducing the workload for experimental network analyses and functional genomics. RESULTS: Genome-scale 13C flux analysis revealed that about half of the 745 biochemical reactions were active during growth on glucose, but that alternative pathways exist for only 51 gene-encoded reactions with significant flux. These flexible reactions identified in silico are key targets for experimental flux analysis, and we present the first large-scale metabolic flux data for yeast, covering half of these mutants during growth on glucose. The metabolic lesions were often counteracted by flux rerouting, but knockout of cofactor-dependent reactions, as in the adh1, ald6, cox5A, fum1, mdh1, pda1, and zwf1 mutations, caused flux responses in more distant parts of the network. By integrating computational analyses, flux data, and physiological phenotypes of all mutants in active reactions, we quantified the relative importance of 'genetic buffering' through alternative pathways and network redundancy through duplicate genes for genetic robustness of the network. CONCLUSIONS: The apparent dispensability of knockout mutants with metabolic function is explained by gene inactivity under a particular condition in about half of the cases. For the remaining 207 viable mutants of active reactions, network redundancy through duplicate genes was the major (75%) and alternative pathways the minor (25%) molecular mechanism of genetic network robustness in S. cerevisiae.

Carbon↗

PubMed2004

Genome-scale in silico models of E. coli have multiple equivalent phenotypic states: assessment of correlated reaction subsets that comprise network states.

The constraint-based analysis of genome-scale metabolic and regulatory networks has been successful in predicting phenotypes and useful for analyzing high-throughput data sets. Within this modeling framework, linear optimization has been used to study genome-scale metabolic models, resulting in the enumeration of single optimal solutions describing the best use of the network to support growth. Here mixed-integer linear programming was used to calculate and study a subset of the alternate optimal solutions for a genome-scale metabolic model of Escherichia coli (iJR904) under a wide variety of environmental conditions. Analysis of the calculated sets of optimal solutions found that: (1) only a small subset of reactions in the network have variable fluxes across optima; (2) sets of reactions that are always used together in optimal solutions, correlated reaction sets, showed moderate agreement with the currently known transcriptional regulatory structure in E. coli and available expression data, and (3) reactions that are used under certain environmental conditions can provide clues about network regulatory needs. In addition, calculation of suboptimal flux distributions, using flux variability analysis, identified reactions which are used under significantly more environmental conditions suboptimally than optimally. Together these results demonstrate the utilization of reactions in genome-scale models under a variety of different growth conditions.

Aerobiosis↗

PubMed2025

Simple scaling laws control the genetic architectures of human complex traits.

Genome-wide association studies have revealed that the genetic architectures of complex traits vary widely, including in terms of the numbers, effect sizes, and allele frequencies of significant hits. However, at present we lack a principled way of understanding the similarities and differences among traits. Here, we describe a probabilistic model that combines the effects of mutation, drift, and stabilizing selection at individual sites with a genome-scale model of phenotypic variation. In this model, the architecture of a trait arises from the distribution of selection coefficients of mutations and from two scaling parameters. We fit this model for 95 highly polygenic quantitative traits of different kinds from the UK Biobank. Notably, we infer that all these traits have fairly similar, though not identical, distributions of selection coefficients. This similarity suggests that differences in architectures of highly polygenic traits arise mainly from the two scaling parameters: the mutational target size and heritability per site, which vary by orders of magnitude among traits. When these two scale factors are accounted for, we find that the architectures of all 95 traits are very similar.

Humans↗

PubMed2006

Experimental and computational assessment of conditionally essential genes in Escherichia coli.

Genome-wide gene essentiality data sets are becoming available for Escherichia coli, but these data sets have yet to be analyzed in the context of a genome scale model. Here, we present an integrative model-driven analysis of the Keio E. coli mutant collection screened in this study on glycerol-supplemented minimal medium. Out of 3,888 single-deletion mutants tested, 119 mutants were unable to grow on glycerol minimal medium. These conditionally essential genes were then evaluated using a genome scale metabolic and transcriptional-regulatory model of E. coli, and it was found that the model made the correct prediction in approximately 91% of the cases. The discrepancies between model predictions and experimental results were analyzed in detail to indicate where model improvements could be made or where the current literature lacks an explanation for the observed phenotypes. The identified set of essential genes and their model-based analysis indicates that our current understanding of the roles these essential genes play is relatively clear and complete. Furthermore, by analyzing the data set in terms of metabolic subsystems across multiple genomes, we can project which metabolic pathways are likely to play equally important roles in other organisms. Overall, this work establishes a paradigm that will drive model enhancement while simultaneously generating hypotheses that will ultimately lead to a better understanding of the organism.

Bacterial Proteins↗

PubMed2006

ScrumPy: metabolic modelling with Python.

ScrumPy is a software package used for the definition and analysis of metabolic models. It is written using the Python programming language that is also used as a user interface. ScrumPy has features for both kinetic and structural modelling, but the emphasis is on structural modelling and those features of most relevance to analysis of large (genome-scale) models. The aim is at describing ScrumPy's functionality to readers with some knowledge of metabolic modelling, but implementation, programming and other computational details are omitted. ScrumPy is released under the Gnu Public Licence, and available for download from http://mudshark.brookes.ac.uk/ ScrumPy.

Adaptation, Physiological↗

PubMed2006

GEM System: automatic prototyping of cell-wide metabolic pathway models from genomes.

BACKGROUND: Successful realization of a "systems biology" approach to analyzing cells is a grand challenge for our understanding of life. However, current modeling approaches to cell simulation are labor-intensive, manual affairs, and therefore constitute a major bottleneck in the evolution of computational cell biology. RESULTS: We developed the Genome-based Modeling (GEM) System for the purpose of automatically prototyping simulation models of cell-wide metabolic pathways from genome sequences and other public biological information. Models generated by the GEM System include an entire Escherichia coli metabolism model comprising 968 reactions of 1195 metabolites, achieving 100% coverage when compared with the KEGG database, 92.38% with the EcoCyc database, and 95.06% with iJR904 genome-scale model. CONCLUSION: The GEM System prototypes qualitative models to reduce the labor-intensive tasks required for systems biology research. Models of over 90 bacterial genomes are available at our web site.

Chromosome Mapping↗

PubMed2004

Integrated analysis of metabolic phenotypes in Saccharomyces cerevisiae.

BACKGROUND: The yeast Saccharomyces cerevisiae is an important microorganism for both industrial processes and scientific research. Consequently, there have been extensive efforts to characterize its cellular processes. In order to fully understand the relationship between yeast's genome and its physiology, the stockpiles of diverse biological data sets that describe its cellular components and phenotypic behavior must be integrated at the genome-scale. Genome-scale metabolic networks have been reconstructed for several microorganisms, including S. cerevisiae, and the properties of these networks have been successfully analyzed using a variety of constraint-based methods. Phenotypic phase plane analysis is a constraint-based method which provides a global view of how optimal growth rates are affected by changes in two environmental variables such as a carbon and an oxygen uptake rate. Some applications of phenotypic phase plane analysis include the study of optimal growth rates and of network capacity and function. RESULTS: In this study, the Saccharomyces cerevisiae genome-scale metabolic network was used to formulate a phenotypic phase plane that displays the maximum allowable growth rate and distinct patterns of metabolic pathway utilization for all combinations of glucose and oxygen uptake rates. In silico predictions of growth rate and secretion rates and in vivo data for three separate growth conditions (aerobic glucose-limited, oxidative-fermentative, and microaerobic) were concordant. CONCLUSIONS: Taken together, this study examines the function and capacity of yeast's metabolic machinery and shows that the phenotypic phase plane can be used to accurately predict metabolic phenotypes and to interpret experimental data in the context of a genome-scale model.

Fermentation↗

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↗

PubMed2026

GENKI: A generative framework for scalable and robust metabolic kinetic modeling.

GENKI (Generative ENsemble KPI-Informed) is a variational autoencoder-based framework for large-scale kinetic modeling of metabolism. Developed for metabolic engineering applications, GENKI is designed to improve the recovery of kinetically feasible models that reproduce experimentally observed phenotypes under genetic and environmental perturbations. The framework is trained on feasible kinetic model ensembles and uses phenotype-based key performance indicators (KPIs), derived from multi-omics and bioprocess data, to label and enrich models according to their agreement with mutant and condition-specific observations. This enables targeted generation of biologically relevant parameter sets with improved predictive performance. Crucially, GENKI recovers kinetic parameter sets that jointly reproduce wild-type and multiple perturbed physiologies within a single model. We apply GENKI to large-scale kinetic models of Escherichia coli and Saccharomyces cerevisiae under enzyme perturbations and oxygen shifts. In both systems, GENKI enriches kinetic ensembles with models that more accurately reproduce experimentally observed physiologies across multiple perturbations and conditions. GENKI therefore provides a practical framework for perturbation-aware kinetic model refinement within iterative Design-Build-Test-Learn workflows.

DBTL↗

PubMed2006

Modeling methanogenesis with a genome-scale metabolic reconstruction of Methanosarcina barkeri.

We present a genome-scale metabolic model for the archaeal methanogen Methanosarcina barkeri. We characterize the metabolic network and compare it to reconstructions from the prokaryotic, eukaryotic and archaeal domains. Using the model in conjunction with constraint-based methods, we simulate the metabolic fluxes and resulting phenotypes induced by different environmental and genetic conditions. This represents the first large-scale simulation of either a methanogen or an archaeal species. Model predictions are validated by comparison to experimental growth measurements and phenotypes of M. barkeri on different substrates. The predicted growth phenotypes for wild type and mutants of the methanogenic pathway have a high level of agreement with experimental findings. We further examine the efficiency of the energy-conserving reactions in the methanogenic pathway, specifically the Ech hydrogenase reaction, and determine a stoichiometry for the nitrogenase reaction. This work demonstrates that a reconstructed metabolic network can serve as an analysis platform to predict cellular phenotypes, characterize methanogenic growth, improve the genome annotation and further uncover the metabolic characteristics of methanogenesis.

Genome, Bacterial↗

PubMed2026

Model-driven analysis reveals oxidative stress adaptation enabling efficient energy utilization in a Crabtree-negative Saccharomyces cerevisiae.

Although abolishing the Crabtree effect in Saccharomyces cerevisiae through a pyruvate dehydrogenase bypass eliminates carbon loss through ethanol overflow metabolism, it compromises growth rates. While the Crabtree effect has been a valuable natural adaptation, it is energetically inferior to respiration and is generally undesirable in cell factories engineered to produce assimilatory compounds. Restoring growth efficiency in Crabtree-negative strains remains a central challenge. Through adaptive laboratory evolution of the engineered strain (sZJD23) and subsequent reverse engineering, a variant (sZJD28) with markedly improved growth was identified. This improvement is driven primarily by a mutation in MED2 (encoding a Mediator complex subunit) and, to a lesser extent, a mutation in GPD1 (encoding glycerol-3-phosphate dehydrogenase). By integrating quantitative proteomics with enzyme-constrained genome-scale modelling, we demonstrate that these mutations jointly enable a more efficient mode of oxidative stress adaptation and energy utilization. The GPD1 mutation suppresses a protein-costly, suboptimal NAD⁺-recycling strategy reliant on glycerol synthesis, while the MED2 mutation reshapes the oxidative stress response towards peroxisomal detoxification. Collectively, these adjustments optimize metabolic flux distribution and reduce protein costs in energy metabolism, thereby increasing ATP availability. Our findings reveal how coordinated mutations in regulatory and metabolic genes restore growth fitness in engineered Crabtree-negative yeast.

Saccharomyces cerevisiae↗

PubMed2006

Multiple mapping method: a novel approach to the sequence-to-structure alignment problem in comparative protein structure modeling.

A major bottleneck in comparative protein structure modeling is the quality of input alignment between the target sequence and the template structure. A number of alignment methods are available, but none of these techniques produce consistently good solutions for all cases. Alignments produced by alternative methods may be superior in certain segments but inferior in others when compared to each other; therefore, an accurate solution often requires an optimal combination of them. To address this problem, we have developed a new approach, Multiple Mapping Method (MMM). The algorithm first identifies the alternatively aligned regions from a set of input alignments. These alternatively aligned segments are scored using a composite scoring function, which determines their fitness within the structural environment of the template. The best scoring regions from a set of alternative segments are combined with the core part of the alignments to produce the final MMM alignment. The algorithm was tested on a dataset of 1400 protein pairs using 11 combinations of two to four alignment methods. In all cases MMM showed statistically significant improvement by reducing alignment errors in the range of 3 to 17%. MMM also compared favorably over two alignment meta-servers. The algorithm is computationally efficient; therefore, it is a suitable tool for genome scale modeling studies.

Algorithms↗

PubMed2006

Systems approach to refining genome annotation.

Genome-scale models of Escherichia coli K-12 MG1655 metabolism have been able to predict growth phenotypes in most, but not all, defined growth environments. Here we introduce the use of an optimization-based algorithm that predicts the missing reactions that are required to reconcile computation and experiment when they disagree. The computer-generated hypotheses for missing reactions were verified experimentally in five cases, leading to the functional assignment of eight ORFs (yjjLMN, yeaTU, dctA, idnT, and putP) with two new enzymatic activities and four transport functions. This study thus demonstrates the use of systems analysis to discover metabolic and transport functions and their genetic basis by a combination of experimental and computational approaches.

Biological Transport↗

PubMed2003

Optknock: a bilevel programming framework for identifying gene knockout strategies for microbial strain optimization.

The advent of genome-scale models of metabolism has laid the foundation for the development of computational procedures for suggesting genetic manipulations that lead to overproduction. In this work, the computational OptKnock framework is introduced for suggesting gene deletion strategies leading to the overproduction of chemicals or biochemicals in E. coli. This is accomplished by ensuring that a drain towards growth resources (i.e., carbon, redox potential, and energy) must be accompanied, due to stoichiometry, by the production of a desired product. Computational results for gene deletions for succinate, lactate, and 1,3-propanediol (PDO) production are in good agreement with mutant strains published in the literature. While some of the suggested deletion strategies are straightforward and involve eliminating competing reaction pathways, many others suggest complex and nonintuitive mechanisms of compensating for the removed functionalities. Finally, the OptKnock procedure, by coupling biomass formation with chemical production, hints at a growth selection/adaptation system for indirectly evolving overproducing mutants.

Algorithms↗