Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Metabolic modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A standardized workflow for kinetic metabolic model curation and dissemination

Kinetic metabolic models provide invaluable insights into cellular metabolism, supporting applications in synthetic biology, metabolic engineering, and systems biology. However, reproducibility and utility of these models hinge on clear and rigorous documentation, standardized annotation, and accessible visualization. This paper presents a workflow for building, annotating, visualizing, and sharing kinetic metabolic models. Our method integrates community standards and open-source tools to ensure reproducibility, interoperability, and user accessibility. This procedure enables researchers to produce reusable and well-documented kinetic models, advancing their role as powerful tools in metabolic research.

Cook, Margaret [Univ. of Washington, Seattle, WA (↗
DOE OSTI2025

TUTORIAL: A new custom metabolic model for iron-oxidizing bacteria

In this tutorial narrative, we introduce a novel template developed to enable the creation of stoichiometric genome-scale metabolic models for iron-oxidizing bacteria. We demonstrate the development of this template by applying it to Sideroxydans lithotrophicus ES-1, and validate our model using transcriptomic data (Published in Zhou et al., 2022 AEM). Below, we further show that our template facilitates the modeling of mixotrophic iron-oxidizing bacteria and metagenome-assembled genomes (MAGs), by applying our template to the MAG of the mixotrophic iron oxidizer Leptothrix ochracea (Published in Tothero et al, 2024). This work represents the first instance of a generalized and adaptable template for modeling diverse iron-oxidizing microbial systems, expanding the accessibility and applicability of metabolic modeling in this field.

genome-scale model↗
DOE OSTI2025

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 ( i Thauera861), 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. i Thauera861 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.

Tec-Campos, Diego (ORCID:0000000188194150)↗
DOE OSTI2024

Quantifying microbial roles in environmental iron oxidation via an integrated kinetics, `omics and metabolic modeling study (Final Report)

Iron oxyhydroxides are extremely reactive components of environmental systems, and therefore exert a strong influence on biogeochemical cycles. These oxyhydroxides strongly adsorb many biologically-relevant elements, including organic carbon and phosphate, as well as a wide range of metals including uranium and actinide species. Thus, the formation mechanism of iron oxyhydroxides are key to understanding both nutrient and contaminant cycling. Microorganisms can catalyze iron oxidation and promote the formation of Fe biominerals and thus are increasingly recognized as important players in biogeochemical cycling. However, it is completely unknown how much of environmental iron oxidation is biologically mediated versus abiotic, and various challenges in studying microbial iron oxidation have hindered accurate incorporation into hydrobiogeochemical models. The overarching goal of our work was to quantify and constrain microbial iron oxidation rates and use ‘omics to gain insight into the controls on this process, while developing tools to enable integration of biotic iron oxidation into hydrobiogeochemical models. Our work focused on the Savannah River Site (SRS) in South Carolina, where extensive microbial iron oxidation has been observed. At Tims Branch, part of the Argonne National Laboratory Wetland Hydrobiogeochemistry Science Focus Area (Argonne SFA), where groundwater discharges into a stream, iron-oxidizing microbial mats form and appear to be a major sink of uranium. In the wetlands that surround Tims Branch, there are wide swaths of iron microbial mats and flocs (mobilized mat). We measured biotic and abiotic iron oxidation rates using mats and water sampled from these sites and found that iron oxidation is primarily carried out by chemolithotrophic microorganisms. The resulting rate constants can be incorporated into models. These mats were characterized by metagenomics and metatranscriptomics, which showed that aerobic chemolithotrophs were the dominant iron-oxidizing bacteria (FeOB), and these included Gallionellaceae and Leptothrix, and possibly Rhodoferax, which is known as an Fe-reducer but may also oxidize Fe(II). This demonstrated that diverse FeOB can coexist and suggests that there are a range of niches and therefore drivers of chemolithotrophic iron oxidation. Analysis of reconstructed genomes strongly suggests that a major factor in diversity is carbon source, as genomes contained varied pathways for autotrophy and heterotrophy. We performed an in-depth analysis of Leptothrix ochracea genomes, since this sheath-former is one of the primary mat builders, yet its physiology remained unresolved. A combination of genomics, transcriptomics, and metabolic modeling suggest that L. ochracea grows mixotrophically using a combination of Fe(II) and organics for energy and both inorganic and organic carbon to create biomass. This contrasts with the largely autotrophic Gallionellaceae (Gallionella, Sideroxydans, and Ferriphaselus) also present in the mats and flocs. Remarkably, multiple FeOB, both Leptothrix and Gallionellaceae, showed activity in response to Fe(II) in live mat incubations. We tracked the gene expression of individual MAGs to Fe(II) and found that various autotrophic and heterotrophic FeOB responded to Fe(II), increasing expression of both carbon fixation and organic utilization genes. The results of the integrated field, kinetics, and omics studies give detailed insight into 1) the taxa that oxidize Fe, and 2) how they connect Fe, C, and N cycles. Towards the goal of connecting omics data to hydrobiogeochemical models, we worked with the KBase team to create a template metabolic model for chemolithotrophic iron oxidation. We initially modeled the well-characterized isolate Gallionellaceae Sideroxydans lithotrophicus, and also applied the model to the mixotroph L. ochracea. In all, we have characterized diverse FeOB in a representative wetland system and solved key problems that enable better incorporation of iron-oxidizing microbes into hydrobiogeochemical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modification and analysis of context-specific genome-scale metabolic models: methane-utilizing microbial chassis as a case study

ABSTRACT Context-specific genome-scale model (CS-GSM) reconstruction is becoming an efficient strategy for integrating and cross-comparing experimental multi-scale data to explore the relationship between cellular genotypes, facilitating fundamental or applied research discoveries. However, the application of CS modeling for non-conventional microbes is still challenging. Here, we present a graphical user interface that integrates COBRApy, EscherPy, and RIPTiDe, Python-based tools within the BioUML platform, and streamlines the reconstruction and interrogation of the CS genome-scale metabolic frameworks via Jupyter Notebook. The approach was tested using -omics data collected for Methylotuvimicrobium alcaliphilum 20Z R , a prominent microbial chassis for methane capturing and valorization. We optimized the previously reconstructed whole genome-scale metabolic network by adjusting the flux distribution using gene expression data. The outputs of the automatically reconstructed CS metabolic network were comparable to manually optimized i IA409 models for Ca-growth conditions. However, the CS model questions the reversibility of the phosphoketolase pathway and suggests higher flux via primary oxidation pathways. The model also highlighted unresolved carbon partitioning between assimilatory and catabolic pathways at the formaldehyde-formate node. Only a very few genes and only one enzyme with a predicted function in C1 metabolism, a homolog of the formaldehyde oxidation enzyme ( fae1-2 ), showed a significant change in expression in La-growth conditions. The CS-GSM predictions agreed with the experimental measurements under the assumption that the Fae1-2 is a part of the tetrahydrofolate-linked pathway. The cellular roles of the tungsten (W)-dependent formate dehydrogenase ( fdhAB ) and fae homologs ( fae1-2 and fae3 ) were investigated via mutagenesis. The phenotype of the f dhAB mutant followed the model prediction. Furthermore, a more significant reduction of the biomass yield was observed during growth in La-supplemented media, confirming a higher flux through formate. M. alcaliphilum 20Z R mutants lacking fae1-2 did not display any significant defects in methane or methanol-dependent growth. However, contrary to fae1, the fae1-2 homolog failed to restore the formaldehyde-activating enzyme function in complementation tests. Overall, the presented data suggest that the developed computational workflow supports the reconstruction and validation of CS-GSM networks of non-model microbes. IMPORTANCE The interrogation of various types of data is a routine strategy to explore the relationship between genotype and phenotype. An efficient approach for integrating and cross-comparing experimental multi-scale data in the context of whole-genome-based metabolic network reconstruction becomes a powerful tool that facilitates fundamental and applied research discoveries. The present study describes the reconstruction of a context-specific (CS) model for the methane-utilizing bacterium, Methylotuvimicrobium alcaliphilum 20Z R . M. alcaliphilum 20Z R is becoming an attractive microbial platform for the production of biofuels, chemicals, pharmaceuticals, and bio-sorbents for capturing atmospheric methane. We demonstrate that this pipeline can help reconstruct metabolic models that are similar to manually curated networks. Furthermore, the model is able to highlight previously overlooked pathways, thus advancing fundamental knowledge of non-model microbial systems or promoting their development toward biotechnological or environmental implementations.

Kulyashov, M. A.↗
DOE OSTI2026

Synergy and antagonism in a genome-scale model of metabolic hijacking by bacteriophages

Bacteriophage auxiliary metabolic genes (AMGs) alter host metabolism by hijacking reactions, but previous studies mostly inferred their roles from annotations, ignoring system-wide impacts and phage production. Here we integrate AMGs and phage assembly into a genome-scale metabolic model of Prochloroccocus marinus MED4 infected by P-HM2. We show that 17 directly hijacked reactions substantially affect more than 30% of the reactions in MED4 metabolism, including carbon fixation, photosynthesis, and nucleotide synthesis, distinguishing these AMGs as either phage aligned—shifting feasible reaction velocities in accordance with maximal phage production—or phage antialigned. Pareto optimization reveals that phage-aligned reactions alter phage-host growth trade-offs, while phage-antialigned reactions do not. We experimentally validate our predictions of system-level AMG impacts by measuring the N-dependent effect of P-HM2 cp12 expression on growth in a model relative of the genetically intractable MED4, Synechococcus elongatus. We also show that AMGs’ indirect impacts are synergistically and antagonistically coupled, providing systems-level insight into AMG perturbations and highlighting how nontrivial cascading effects shape host metabolism.

Rozum, Jordan C. [Pacific Northwest National Labor↗
DOE OSTI2026

Integrative Modeling and Analysis of Fungal Central Carbon Metabolism

Over a thousand fungal genomes have been sequenced, yet manually curated genome-scale metabolic models (GEMs) are available for only a limited number of species. Moreover, these models have often been developed independently, leading to inconsistencies in namespaces, compartment definitions, and pathway representations that hinder comparative analysis, the systematic reuse of prior curation efforts, and the integration of consolidated metabolic knowledge. Here, we present the Consolidated Fungal Core Metabolism Model (CFCMM), constructed by integrating thirteen published fungal models spanning Ascomycota, Mucoromycota, and both Crabtree-positive and Crabtree-negative yeasts. We harmonized metabolites and reactions into a non-redundant shared ModelSEED ontological space, standardized compartmentalization, and refined gene–protein–reaction (GPR) rules. Using pathway-level visualization and systematic gap detection, we further improved the integrated network through literature-guided curation to correct stoichiometry, stereospecificity, and pathway architecture. Orthologous protein family reconstruction and functional annotation workflows were used to validate and inform GPR associations, with particular emphasis on ambiguous enzyme superfamilies and membrane-associated components. Using the resulting CFCMM, we built high-quality central carbon core models for each fungus and performed flux balance analysis to quantify ATP-yield variation under aerobic and anaerobic conditions, explicitly evaluating scenarios driven by differences in electron transport chain (ETC) composition. Simulations reproduced the expected fermentative yield of approximately 2 mmol ATP per mmol glucose under anaerobic conditions and separated the thirteen fungi into two bioenergetic groups under aerobic respiration based on Complex I status, with predicted yields of approximately 30 versus 22 mmol ATP per mmol glucose. Forcing flux through the alternative oxidase bypass further reduced ATP yields to approximately 12 and 4 mmol ATP per mmol glucose in Complex I-containing and Complex I-lacking fungi, respectively. Collectively, this work provides a manually curated, ModelSEED-consistent, and extensible fungal core metabolic template, deployed in DOE KBase as a resource for automated reconstruction of central carbon core models from any sequenced fungal genome. In addition, the CFCMM provides modular components for developing GEMs with more accurate energy predictions and enables robust comparative analyses of fungal bioenergetics and core metabolic diversity

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Model of metabolism and gene expression predicts proteome allocation in Pseudomonas putida

Abstract The genome-scale model of metabolism and gene expression (ME-model) forPseudomonas putidaKT2440,iPpu1676-ME, provides a comprehensive representation of biosynthetic costs and proteome allocation. Compared to a metabolic-only model,iPpu1676-ME significantly expands on gene expression, macromolecular assembly, and cofactor utilization, enabling accurate growth predictions without additional constraints. Multi-omics analysis using RNA sequencing and ribosomal profiling data revealed translational prioritization inP. putida, with core pathways, such as nicotinamide biosynthesis and queuosine metabolism, exhibiting higher translational efficiency, while secondary pathways displayed lower priority. Notably, the ME-model significantly outperformed the M-model in alignment with multi-omics data, thereby validating its predictive capacity. Thus,iPpu1676-ME offers valuable insights intoP. putida’s proteome allocation and presents a powerful tool for understanding resource allocation in this industrially relevant microorganism.

Mathematical & Computational Biology↗
DOE OSTI2026

MONet/1000 Soils metagenome pathway modelling narrative w/ auto batch import

This narrative performs metabolic modeling and flux balance analysis (FBA) using metagenome-assembled genomes (MAGs) from the 1000 Soils samples, as described by Song et al. (2026, accepted). The set of MAGs (in FASTA format) is converted into a set of assembly objects compatible with functional annotation via RASTtk, yielding a set of genome objects that undergo metabolic modeling via OMEGGA. These genome objects are then used to conduct FBA, generating tables of metabolite uptake rates across the MAGs under investigation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

pnnl-predictive-phenomics/csc052-gem

Genome-Scale Metabolic Model Continuous Validation with Memote for CarbStor Community Member Bacillus These repositories contain the continuous validation environment for an organism-specific genome-scale metabolic model (GEM) using Memote. Memote is a software tool that provides a suite of tests to ensure the quality and consistency of metabolic models. By integrating Memote into a continuous integration (CI) workflow, we can automatically validate updates to the GEM, ensuring that model modifications improve or maintain the model's integrity.

Torres, Victor E.↗
DOE OSTI2024

pnnl-predictive-phenomics/csc040-gem

Genome-Scale Metabolic Model Continuous Validation with Memote for CarbStor Community Member Rhodococcus These repositories contain the continuous validation environment for an organism-specific genome-scale metabolic model (GEM) using Memote. Memote is a software tool that provides a suite of tests to ensure the quality and consistency of metabolic models. By integrating Memote into a continuous integration (CI) workflow, we can automatically validate updates to the GEM, ensuring that model modifications improve or maintain the model's integrity.

McNaughton, Andrew [@PNNL]↗
DOE OSTI2024

pnnl-predictive-phenomics/csc043-gem

Genome-Scale Metabolic Model Continuous Validation with Memote for CarbStor Community Member Paenibacillus These repositories contain the continuous validation environment for an organism-specific genome-scale metabolic model (GEM) using Memote. Memote is a software tool that provides a suite of tests to ensure the quality and consistency of metabolic models. By integrating Memote into a continuous integration (CI) workflow, we can automatically validate updates to the GEM, ensuring that model modifications improve or maintain the model's integrity.

Zucker, Jeremy [Pacific Northwest National Laborat↗
DOE OSTI2024

pnnl-predictive-phenomics/csc031-gem

Genome-Scale Metabolic Model of CarbStor Community member Microbacterium (csc031) Continuous Validation with Memote These repositories contain the continuous validation environment for an organism-specific genome-scale metabolic model (GEM) using Memote. Memote is a software tool that provides a suite of tests to ensure the quality and consistency of metabolic models. By integrating Memote into a continuous integration (CI) workflow, we can automatically validate updates to the GEM, ensuring that model modifications improve or maintain the model's integrity

McNaughton, Andrew [@PNNL]↗
DOE OSTI2024

pnnl-predictive-phenomics/csc009-gem

Genome-Scale Metabolic Model of CarbStore Community member Curtobacterium (csc009) Continuous Validation with Memote These repositories contain the continuous validation environment for an organism-specific genome-scale metabolic model (GEM) using Memote. Memote is a software tool that provides a suite of tests to ensure the quality and consistency of metabolic models. By integrating Memote into a continuous integration (CI) workflow, we can automatically validate updates to the GEM, ensuring that model modifications improve or maintain the model's integrity

Lin, Tesia↗
DOE OSTI2025

Engineering the hyperthermophilic archaeon Pyrococcus furiosus for 1-propanol production

ABSTRACT Society relies heavily on chemicals traditionally produced through the refinement of fossil fuels. The conversion of renewable biomass to value-added chemicals by microbes, particularly hyperthermophiles (T opt ≥80°C), offers a renewable alternative to this traditional approach. Herein, we describe the engineering of the hyperthermophilic archaeon Pyrococcus furiosus , which grows optimally (T opt ) at 100°C, for the conversion of sugar to 1-propanol. This was accomplished by constructing a hybrid metabolic pathway consisting of two native and seven heterologously produced enzymes to convert acetyl-CoA from carbohydrate metabolism to 1-propanol. A total of eleven foreign genes from two other organisms were utilized, one from the thermophilic bacterium Thermoanaerobacter sp. strain X514 and 10 from the thermoacidophilic archaeon Metallosphaera sedula, both of which grow optimally near 70°C. The recombinant P. furiosus strain produced 1-propanol at similar concentrations (up to ~1 mM) when incubated at 75°C to activate the gene products of Thermoanaerobacter sp. strain X514 and M. sedula and by initially incubating at 95°C for P. furiosus growth and then subsequently returning to 75°C to promote 1-propanol formation. Note that 1-propanol was not produced if the culture was grown only at 95°C. This work has the potential for future optimization through harnessing the genome-scale metabolic model of P. furiosus that was used herein to identify engineering targets to increase 1-propanol titer. IMPORTANCE As petroleum reserves become increasingly strained, the development of renewable alternatives to traditional chemical synthesis becomes more important. In this work, a high-temperature biological system for sugar to 1-propanol conversion was demonstrated by metabolic engineering of the hyperthermophilic archaeon Pyrococcus furiosus (T opt 100°C). The engineered strain produced 1-propanol by temperature shifting from 75°C to 95°C and then back to 75°C to accommodate the temperature ranges for native and foreign proteins associated with 1-propanol biosynthesis. Genome-scale metabolic modeling informed the carbon and reductant flux in the system, identified potential factors limiting 1-propanol production, and revealed potential optimization targets.

O'Quinn, Hailey C. (ORCID:0009000274865755)↗
DOE OSTI2026

Conjugation-based genome engineering enables rapid prototyping and bioproduction in non-model bacteria

Abstract Non-model bacteria offer unique metabolic capabilities for sustainable bioproduction, yet their limited genetic accessibility hinders systematic strain development. Here we present conjugation-based serine recombinase-assisted genome engineering (cSAGE), a broad-host-range platform that enables predictable, iterative genomic integration in transformation-resistant bacteria. cSAGE combines conjugative DNA delivery, standardized low-copy vectors, orthogonal recombinases, and modular genetic parts to support rapid pathway assembly and cross-host benchmarking. Using purple nonsulfur bacteria as a testbed, we integrate promoter engineering, multi-payload genome modification, and genome-scale metabolic modeling to empirically evaluate host-dependent pathway performance. Applying this workflow, we identify strain-specific differences in photosynthetic conversion of lignin-derived p -coumarate to the thermoplastic precursor p -vinylphenol. By enabling genome engineering and functional comparison across diverse bacteria using a single plasmid system, cSAGE provides a general framework for non-model strain prototyping and biotransformation discovery.

Guzman, Michael S. [Department of Chemical Enginee↗
DOE OSTI2025

Developing a media formulation to sustain ex vivo chloroplast function

Chloroplasts are critical organelles in plants and algae responsible for accumulating biomass through photosynthetic carbon fixation and cellular maintenance through metabolism in the cell. Chloroplasts are increasingly appreciated for their role in biomanufacturing, as they can produce many useful molecules, and a deeper understanding of chloroplast regulation and function would provide more insight for the biotechnological applications of these organelles. However, traditional genetic approaches to manipulate chloroplasts are slow, and generation of transgenic organisms to study their function can take weeks to months, significantly delaying the pace of research. To develop chloroplasts themselves as a quicker and more defined platform, we isolated chloroplasts from the green algae, Chlamydomonas reinhardtii, and examined their photosynthetic function after extraction. Combined with a metabolic modeling approach using flux-balance analysis, we identified key metabolic reactions essential to chloroplast function and leveraged this information into reagents that can be used in a “chloroplast media” capable of maintaining chloroplast photosynthetic function over time ex vivo compared to buffer alone. We envision this could serve as a model platform to enable more rapid design-build-test-learn cycles to study and improve chloroplast function in combination with genetic modifications and potentially as a starting point for the bottom-up design of a synthetic organelle-containing cell.

Chlamydomonas reinhardtii↗