Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “genome-scale metabolic model”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

PET-FBA: A lightweight enzyme allocation and thermodynamics-constrained flux analysis approach to explore Escherichia coli metabolic adaptation to intracellular acidification

Escherichia coli employs diverse strategies to adapt to acidic environments that disrupt enzyme activity and the thermodynamic feasibility of essential reactions. To understand the impact of pH stress on cell metabolism, we present the PET-FBA (pH-, Enzyme protein allocation-, and Thermodynamics-constrained Flux Balance Analysis) framework. PET-FBA extends genome-scale modeling by integrating enzyme protein costs and reaction Gibbs free energy changes. Additionally, by incorporating pH-dependent enzyme kinetics in response to intracellular acidification, this framework enables the simulation of E. coli's metabolic adjustments across varying external pH levels. The model's accuracy is validated by comparing in silico growth simulations with experimental measurements under both anaerobic and aerobic conditions, as well as in silico gene knockouts of essential genes. By explicitly incorporating pH effects, our model accurately replicates the metabolic shift towards lactate production as the primary fermentation product at low pH in anaerobic conditions. This shift is only predicted when enzyme kinetics are dynamically adjusted as a function of pH. Further analysis revealed that this shift can be attributed to the reduced protein efficiency of the acetyl-CoA branch compared to lactate dehydrogenase under acidic stress, which then becomes crucial for maintaining NAD regeneration and cell growth at low pH. Furthermore, we identified strategies for enhancing cell growth under acidic anaerobic conditions by improving the enzyme activity of lactate dehydrogenase and pyruvate formate lyase, which increases NAD production efficiency and reduces enzyme protein allocation costs. Designed as a lightweight yet versatile framework, PET-FBA enables efficient genome-scale metabolic analysis. Using E. coli as a model system, our framework provides a systematic approach to understanding metabolic responses to environmental stress, pinpointing key metabolic bottlenecks, and identifying potential targets for strain optimization.

42 ENGINEERING↗
DOE OSTI2023

Genome-scale metabolic reconstruction of 7,302 human microorganisms for personalized medicine

The human microbiome influences the efficacy and safety of a wide variety of commonly prescribed drugs. Designing precision medicine approaches that incorporate microbial metabolism would require strain- and molecule-resolved, scalable computational modeling. Here, we extend our previous resource of genome-scale metabolic reconstructions of human gut microorganisms with a greatly expanded version. AGORA2 (assembly of gut organisms through reconstruction and analysis, version 2) accounts for 7,302 strains, includes strain-resolved drug degradation and biotransformation capabilities for 98 drugs, and was extensively curated based on comparative genomics and literature searches. The microbial reconstructions performed very well against three independently assembled experimental datasets with an accuracy of 0.72 to 0.84, surpassing other reconstruction resources and predicted known microbial drug transformations with an accuracy of 0.81. We demonstrate that AGORA2 enables personalized, strain-resolved modeling by predicting the drug conversion potential of the gut microbiomes from 616 patients with colorectal cancer and controls, which greatly varied between individuals and correlated with age, sex, body mass index and disease stages. AGORA2 serves as a knowledge base for the human microbiome and paves the way to personalized, predictive analysis of host–microbiome metabolic interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Evaluating proteome allocation of Saccharomyces cerevisiae phenotypes with resource balance analysis

Saccharomyces cerevisiae is an important model organism and a workhorse in bioproduction. Here, we reconstructed a compact and tractable genome-scale resource balance analysis (RBA) model (i.e., named scRBA) to analyze metabolic fluxes and proteome allocation in a computationally efficient manner. Resource capacity models such as scRBA provide the quantitative means to identify bottlenecks in biosynthetic pathways due to enzyme, compartment size, and/or ribosome availability limitations. ATP maintenance rate and in vivo apparent turnover numbers (k app ) were regressed from metabolic flux and protein concentration data to capture observed physiological growth yield and proteome efficiency and allocation, respectively. Estimated parameter values were found to vary with oxygen and nutrient availability. Altogether, this work (i) provides condition-specific model parameters to recapitulate phenotypes corresponding to different extracellular environments, (ii) alludes to the enhancing effect of substrate channeling and post-translational activation on in vivo enzyme efficiency in glycolysis and electron transport chain, and (iii) reveals that the Crabtree effect is underpinned by specific limitations in mitochondrial proteome capacity and secondarily ribosome availability rather than overall proteome capacity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A constraint-based framework for exploring the impact of multireaction dependencies on metabolic functions

Abstract Metabolism operates under physico-chemical constraints that result in multireaction dependencies. Understanding how multireaction dependencies affect metabolic phenotypes remains challenging, hindering their biotechnological applications. Here, we propose the concept of a forcedly balanced complex that allows to efficiently determine the effects of specific multireaction dependencies on metabolic network functions in constrained-based models. Using this concept, we found that the fraction of multireaction dependencies induced by forcedly balanced complexes in genome-scale metabolic networks followed power law with exponential cut-off. We identified forcedly balanced complexes that are lethal in cancer but have little effect on growth in healthy tissue models. In addition, these forcedly balanced complexes are largely specific to models of particular cancer types. Therefore, multireaction dependencies resulting from forced balancing of complexes represent an innovative means to control cancers that, we argue, can be implemented via transporter engineering. The presented constraint-based approaches pave the way for using multireaction dependencies in metabolic engineering for diverse biotechnological applications.

Küken, Anika↗
DOE OSTI2023

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗
DOE OSTI2025

Genome-scale Design and Engineering of Non-model Yeast Organisms for Production of Biofuels and Bioproducts

The overall goal of this project was to develop genome-scale design and engineering tools for two non-model yeast organisms including Rhodotorula toruloides and Issatchenkia orientalis to produce high-levels of fatty acids-derived products and organic acids, respectively. The project was performed between 9/15/2017 and 9/14/2024 (the last two-years were no-cost extensions). The team consisted of Huimin Zhao (Lead PI) and Christopher Rao (Co-PI) from the University of Illinois at Urbana-Champaign (UIUC), Costas Maranas (Co-PI) from the Pennsylvania State University, Joshua Rabinowitz (Co-PI) and Martin Wuhr (Co-PI) from Princeton University, and Yasuo Yoshikuni (Co-PI) from the DOE Joint Genome Institute. The team has made great progress in both tool development and fundamental understanding of these two non-model yeasts. In total, there were 40 research publications (one of them is still under review) and one patent application as well as numerous oral presentations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Metabolic flux and resource balance in the oleaginous yeast Rhodotorula toruloides

The yeast Rhodotorula toruloides is a promising bioproduction organism due to its high lipid yields and ability to grow on cheap and abundant substrates. Quantitative, systems-level assessment of its metabolic activity is accordingly merited. Resource-balance analysis (RBA) models capture not only reaction stoichiometry but also enzyme requirements for catalysis, providing valuable tools for understanding metabolic trade-offs and optimizing metabolic engineering strategies. Here, in this work, we present systems-level measurements of R. toruloides metabolic flux based on isotope tracing and metabolic flux analysis. In combination with new proteomic measurements, these flux data are used to parameterize a genome-scale resource balance model rtRBA. We find that S. cerevisiae and R. toruloides grow at nearly indistinguishable rates using similar biosynthetic but dramatically different central metabolic programs. R. toruloides consumes one-fifth as much glucose, which it metabolizes primarily via the pentose phosphate pathway and TCA cycle unlike primarily glycolysis in S. cerevisiae . Overall, across these two divergent yeasts, protein abundances aligned more closely than metabolic flux. Resource balance modeling of these metabolic programs predicts superior theoretical yields but lower productivities in R. toruloides than S. cerevisiae for industrial chemicals, highlighting the value of rapid glucose uptake for productivity but respiratory metabolism for yields.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

An automated workflow that generates atom mappings for large‐scale metabolic models and its application to Arabidopsis thaliana

SUMMARY Quantification of reaction fluxes of metabolic networks can help us understand how the integration of different metabolic pathways determines cellular functions. Yet, intracellular fluxes cannot be measured directly but are estimated with metabolic flux analysis (MFA), which relies on the patterns of isotope labeling of metabolites in the network. The application of MFA also requires a stoichiometric model with atom mappings that are currently not available for the majority of large‐scale metabolic network models, particularly of plants. While automated approaches such as the Reaction Decoder Toolkit (RDT) can produce atom mappings for individual reactions, tracing the flow of individual atoms of the entire reactions across a metabolic model remains challenging. Here we establish an automated workflow to obtain reliable atom mappings for large‐scale metabolic models by refining the outcome of RDT, and apply the workflow to metabolic models of Arabidopsis thaliana . We demonstrate the accuracy of RDT through a comparative analysis with atom mappings from a large database of biochemical reactions, MetaCyc. We further show the utility of our automated workflow by simulating 15 N isotope enrichment and identifying nitrogen (N)‐containing metabolites which show enrichment patterns that are informative for flux estimation in future 15 N‐MFA studies of A. thaliana . The automated workflow established in this study can be readily expanded to other species for which metabolic models have been established and the resulting atom mappings will facilitate MFA and graph‐theoretic structural analyses with large‐scale metabolic networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

maranasgroup/scRBA

Saccharomyces cerevisiae is an important model organism and a workhorse in bioproduction. Here, we reconstructed a compact and tractable genome-scale resource balance analysis (RBA) model (i.e., named scRBA) to analyze metabolic fluxes and proteome allocation in a computationally efficient manner. Resource capacity models such as scRBA provide the quantitative means to identify bottlenecks in biosynthetic pathways due to enzyme, compartment size, and/or ribosome availability limitations. ATP maintenance rate and in vivo apparent turnover numbers (kapp) were regressed from metabolic flux and protein concentration data to capture observed physiological growth yield and proteome efficiency and allocation, respectively. Estimated parameter values were found to vary with oxygen and nutrient availability. Overall, this work (i) provides condition-specific model parameters to recapitulate phenotypes corresponding to different extracellular environments, (ii) alludes to the enhancing effect of substrate channeling and post-translational activation on in vivo enzyme efficiency in glycolysis and electron transport chain, and (iii) reveals that the Crabtree effect is underpinned by specific limitations in mitochondrial proteome capacity and secondarily ribosome availability rather than overall proteome capacity.

Dinh, Hoang↗
DOE OSTI2022

Advances in engineering algae for biofuel production

While algae demonstrate potential as a sustainable fuel source, low productivities limit the economic realization of algal biofuels. High-throughput strain engineering, omics-informed genome-scale modeling, and microbiome engineering are key technologies for enabling algal biofuels. High-throughput strain engineering efforts generate improved traits, including high biomass productivity and lipid content, in diverse algal species. Genome-scale models, constructed with the aid of omics data, provide insight into metabolic limitations and guide rational algal strain engineering efforts. As outdoor cultivation systems introduce exogenous organisms, microbiome engineering seeks to eliminate harmful organisms and introduce beneficial species. In conclusion, optimizing algal biomass production and lipid content using these technologies may overcome the productivity barrier for the commercialization of algal biofuels.

09 BIOMASS FUELS↗
DOE OSTI2022

Syntrophic Co-Cultures of Clostridium Organisms to Produce Higher Alcohols & Other C6-C8 Metabolites (Final Report)

The goal of this project was to advance the systems biology understanding and predictive modeling of synthetic & syntrophic Clostridium microbial consortia, focusing on elucidation of metabolic networks and environmental signals in the consortia. The project has direct applicability to lignocellulosic-biomass based production of higher alcohols as advanced biofuels and C6-C8 metabolites, that can be used as chemicals or serve as biofuel precursors. This project studied and optimized three synthetic syntrophic systems. It examined the population dynamics using flow-cytometry, time-lapse microscopy and PCR analysis. 13C-based tracer analyses was used to examine the metabolite exchange between the syntrophic cell populations and the impact of those interactions on the transcriptome of the individual populations. To enhance our analytical and predictive capabilities, genome-scale models (GSMs) for these syntrophies was developed. RNAseq data for these syntrophic coculture systems were acquired to enable a molecular level understanding of the syntrophies aiming to identify the genetic networks of each organism in the co-culture and compare those against the networks of pure cultures.

09 BIOMASS FUELS↗
DOE OSTI2025

Competition and cooperation: The plasticity of bacterial interactions across environments

Bacteria live in diverse communities, forming complex networks of interacting species. A central question in bacterial ecology is whether species engage in cooperative or competitive interactions. But this question often neglects the role of the environment. Here, we use genome-scale metabolic networks from two different open-access collections (AGORA and CarveMe) to assess pairwise interactions of different microbes in varying environmental conditions (provision of different environmental compounds). By computationally simulating thousands of environments for 10,000 pairs of bacteria from each collection, we found that most pairs were able to both compete and cooperate depending on the availability of environmental resources. This modeling approach allowed us to determine commonalities between environments that could facilitate the potential for cooperation or competition between a pair of species. Namely, cooperative interactions, especially obligate, were most common in less diverse environments. Further, as compounds were removed from the environment, we found interactions tended to degrade towards obligacy. However, we also found that on average at least one compound could be removed from an environment to switch the interaction from competition to facultative cooperation or vice versa. Together our approach indicates a high degree of plasticity in microbial interactions in response to the availability of environmental resources.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Revealing systematic changes in the transcriptome during the transition from exponential growth to stationary phase

ABSTRACT The composition of bacterial transcriptomes is determined by the transcriptional regulatory network (TRN). The TRN regulates the transition from one physiological state to another. Here, we use independent component analysis to monitor the composition of the transcriptome during the transition from the exponential growth phase to the stationary phase. With Escherichia coli K-12 MG1655 as a model strain, we trigger the transition using carbon, nitrogen, and sulfur starvation. We find that (i) the transition to the stationary phase accompanies common transcriptome changes, including increased stringent responses and reduced production of cellular building blocks and energy regardless of the limiting element; (ii) condition-specific changes are strongly associated with transcriptional regulators ( e.g. , Crp, NtrC, CysB, Cbl) responsible for metabolizing the limiting element; and (iii) the shortage of each limiting element differentially affects the production of amino acids and extracellular polymers. This study demonstrates how the combination of genome-scale datasets and new data analytics reveals the fundamental characteristics of a key transition in the life cycle of bacteria. IMPORTANCE Nutrient limitations are critical environmental perturbations in bacterial physiology. Despite its importance, a detailed understanding of how bacterial transcriptomes are adjusted has been limited. By utilizing independent component analysis (ICA) to decompose transcriptome data, this study reveals key regulatory events that enable bacteria to adapt to nutrient limitations. The findings not only highlight common responses, such as the stringent response, but also condition-specific regulatory shifts associated with carbon, nitrogen, and sulfur starvation. The insights gained from this work advance our knowledge of bacterial physiology, gene regulation, and metabolic adaptation.

Lim, Hyun Gyu (ORCID:0000000204692388)↗
DOE OSTI2024

Transporter annotations are holding up progress in metabolic modeling

Mechanistic, constraint-based models of microbial isolates or communities are a staple in the metabolic analysis toolbox, but predictions about microbe-microbe and microbe-environment interactions are only as good as the accuracy of transporter annotations. A number of hurdles stand in the way of comprehensive functional assignments for membrane transporters. These include general or non-specific substrate assignments, ambiguity in the localization, directionality and reversibility of a transporter, and the many-to-many mapping of substrates, transporters and genes. In this perspective, we summarize progress in both experimental and computational approaches used to determine the function of transporters and consider paths forward that integrate both. Investment in accurate, high-throughput functional characterization is needed to train the next-generation of predictive tools toward genome-scale metabolic network reconstructions that better predict phenotypes and interactions. More reliable predictions in this domain will benefit fields ranging from personalized medicine to metabolic engineering to microbial ecology.

Casey, John↗
DOE OSTI2026

Data for High Yield Production of 3-Hydroxypropionic Acid Using Issatchenkia orientalis

Biomanufacturing provides a more sustainable alternative to fossil-based chemical manufacturing. 3-Hydroxypropionic acid (3HP) is a top Department of Energy value-added chemical and precursor to bioplastics, yet cost-effective microbial production remains elusive. Here, we establish the acid-tolerant yeast Issatchenkia orientalis as a robust host for low-pH 3HP biosynthesis. Genome-scale modeling identifies the β-alanine pathway as optimal, offering the highest theoretical yield and lowest oxygen requirement. Thermodynamic analysis confirms its favorability under acidic conditions. Using sequence similarity network analysis, we discover highly active aspartate 1-decarboxylase (PAND), β-alanine-pyruvate aminotransferase (BAPAT), and 3HP dehydrogenase (YDFG), which significantly improve the pathway efficiency. Next, to further elevate the production, pathway optimization through multi-copy PAND integration, byproduct elimination (knockouts of pyruvate decarboxylase and glycerol-3-phosphate dehydrogenase), and reinforcement of aspartate flux by overexpression of pyruvate carboxylase and aspartate amino transferase improves the titer to 29 g/L in shake flasks. Fed-batch fermentation at pH 4 with low-cost corn steep liquor medium further increases the production to 92 g/L with 0.7 g/g yield and 0.55 g/L/h productivity. Techno-economic analysis indicates that such performance could potentially enable a financially viable process for sustainable acrylic acid production. This work establishes I. orientalis as a next-generation platform for cost-effective 3HP production and paves the way toward industrial commercialization.

Bioproducts↗
DOE OSTI2024

Harnessing photoautotroph-methanotroph interactions for biogas conversion to fuels and chemicals using binary consortia (Project Final Technical Report)

Industrial, municipal, and agricultural waste streams containing stranded organic carbon represent a significant and underutilized feedstock to produce fuels and chemicals. With anaerobic digestion deployed at large scales to capture organic waste streams, over 6 million tons of biogas are available today. However, the utilization of biogas represents a significant challenge due to its low pressure and presence of contaminants such as H 2 S, ammonia, and volatile organic carbon compounds. To tap into this immense potential, effective biotechnologies that co-utilize both CO 2 and CH 4 are needed. Recent studies demonstrated that, in nature, microbial communities have developed a highly efficient way to recover energy and capture carbon from both CH 4 and CO 2 through metabolic coupling of methane oxidation to oxygenic photosynthesis. Using two synthetic methanotroph – photoautotroph (M-P) co-cultures that exhibit stable growth under a broad range of cultivation conditions, in this project we proposed to harness the interspecies interactions within these cocultures for biogas conversion to fuels and chemicals. To facilitate this overarching objective, we aim to develop experimental and computational tools to gain qualitative and quantitative understandings on the interactions and dynamics of the coculture at both systems and molecular levels, and to validate our findings through experiments and mutant development. The fundamental understanding on the interactions and dynamics of the photoautotroph-methanotroph will lay the foundation for the design and optimization of synthetic binary consortia for production of fuels and chemicals from biogas. We expect the knowledge gained from this project may be generally applicable to other cross-feeding binary consortium, and the tools developed can be adapted to study the interactions and dynamics of other multi-organism platforms.

09 BIOMASS FUELS↗
DOE OSTI2026

High yield production of 3-hydroxypropionic acid using Issatchenkia orientalis

Biomanufacturing provides a more sustainable alternative to fossil-based chemical manufacturing. 3-Hydroxypropionic acid (3HP) is a top Department of Energy value-added chemical and precursor to bioplastics, yet cost-effective microbial production remains elusive. Here, we establish the acid-tolerant yeast Issatchenkia orientalis as a robust host for low-pH 3HP biosynthesis. Genome-scale modeling identifies the β-alanine pathway as optimal, offering the highest theoretical yield and lowest oxygen requirement. Thermodynamic analysis confirms its favorability under acidic conditions. Using sequence similarity network analysis, we discover highly active aspartate 1-decarboxylase (PAND), β-alanine-pyruvate aminotransferase (BAPAT), and 3HP dehydrogenase (YDFG), which significantly improve the pathway efficiency. Next, to further elevate the production, pathway optimization through multi-copy PAND integration, byproduct elimination (knockouts of pyruvate decarboxylase and glycerol-3-phosphate dehydrogenase), and reinforcement of aspartate flux by overexpression of pyruvate carboxylase and aspartate amino transferase improves the titer to 29 g/L in shake flasks. Fed-batch fermentation at pH 4 with low-cost corn steep liquor medium further increases the production to 92 g/L with 0.7 g/g yield and 0.55 g/L/h productivity. Techno-economic analysis indicates that such performance could potentially enable a financially viable process for sustainable acrylic acid production. This work establishes I. orientalis as a next-generation platform for cost-effective 3HP production and paves the way toward industrial commercialization.

Biotechnology↗
DOE OSTI2023

Comparative study of two Saccharomyces cerevisiae strains with kinetic models at genome-scale

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by 13 C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

59 BASIC BIOLOGICAL SCIENCES↗