Search NASA⌕ Search

SEARCH · Search NASA

Results for “metabolic flux analysis”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mechanisms of metabolic adaptation in the duckweed Lemna gibba: an integrated metabolic, transcriptomic and flux analysis

Abstract Background Duckweeds are small, rapidly growing aquatic flowering plants. Due to their ability for biomass production at high rates they represent promising candidates for biofuel feedstocks. Duckweeds are also excellent model organisms because they can be maintained in well-defined liquid media, usually reproduce asexually, and because genomic resources are becoming increasingly available. To demonstrate the utility of duckweed for integrated metabolic studies, we examined the metabolic adaptation of growing Lemna gibba cultures to different nutritional conditions. Results To establish a framework for quantitative metabolic research in duckweeds we derived a central carbon metabolism network model of Lemna gibba based on its draft genome. Lemna gibba fronds were grown with nitrate or glutamine as nitrogen source. The two conditions were compared by quantification of growth kinetics, metabolite levels, transcript abundance, as well as by 13 C-metabolic flux analysis. While growing with glutamine, the fronds grew 1.4 times faster and accumulated more protein and less cell wall components compared to plants grown on nitrate. Characterization of photomixotrophic growth by 13 C-metabolic flux analysis showed that, under both metabolic growth conditions, the Calvin-Benson-Bassham cycle and the oxidative pentose-phosphate pathway are highly active, creating a futile cycle with net ATP consumption. Depending on the nitrogen source, substantial reorganization of fluxes around the tricarboxylic acid cycle took place, leading to differential formation of the biosynthetic precursors of the Asp and Gln families of proteinogenic amino acids. Despite the substantial reorganization of fluxes around the tricarboxylic acid cycle, flux changes could largely not be associated with changes in transcripts. Conclusions Through integrated analysis of growth rate, biomass composition, metabolite levels, and metabolic flux, we show that Lemna gibba is an excellent system for quantitative metabolic studies in plants. Our study showed that Lemna gibba adjusts to different nitrogen sources by reorganizing central metabolism. The observed disconnect between gene expression regulation and metabolism underscores the importance of metabolic flux analysis as a tool in such studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

The effects of photosynthetic rate on respiration in light, starch/sucrose partitioning, and other metabolic fluxes within photosynthesis

In the future, plants may encounter increased light and elevated CO 2 levels. How consequent alterations in photosynthetic rates will impact fluxes in photosynthetic carbon metabolism remains uncertain. Respiration in light ( R L ) is pivotal in plant carbon balance and a key parameter in photosynthesis models. Understanding the dynamics of photosynthetic metabolism and R L under varying environmental conditions is essential for optimizing plant growth and agricultural productivity. However, measuring R L under high light and high CO 2 (HLHC) conditions poses challenges using traditional gas exchange methods. In this study, we employed isotopically nonstationary metabolic flux analysis (INST-MFA) to estimate RL and investigate photosynthetic carbon flux, unveiling nuanced adjustments in Camelina sativa under HLHC. Despite numerous flux alterations in HLHC, RL remained stable. HLHC affects several factors influencing RL, such as starch and sucrose partitioning, v o /v c ratio, triose phosphate partitioning, and hexose kinase activity. Analysis of A/C i curve operational points reveals that HLHC’s major changes primarily stem from CO 2 suppressing photorespiration. Integration of these fluxes into a simplified model predicts changes in CBC labeling under HLHC. This study extends our prior discovery that incomplete CBC labeling is due to unlabeled carbon reimported during R L , offering insights into manipulating labeling through adjustments in photosynthetic rates.

Elevated CO2↗
DOE OSTI2026

Host analysis-guided selection and targeted engineering (HASTE) of Lipomyces tetrasporus for the conversion of CO2-derived feedstocks

Efficient and cost-competitive bioproduction calls for utilizing CO2-derived feedstocks, such as products from electro-reduction of CO2 and hydrolysate from lignocellulosic biomass. However, efficiently using all their carbon components, including acetate, glucose, and xylose, remains a challenge. Here, we characterize Lipomyces tetrasporus, a novel, robust yeast strain capable of effectively assimilating these carbon sources. We used an integrated systems biology approach combining ¹³C metabolic flux analysis, dynamic labeling experiments, and RNA sequencing. We conducted the first metabolic flux analysis for glucose, xylose, and acetate catabolism in this species. Dynamic labeling revealed a highly active TCA cycle during acetate metabolism, evidenced by rapid citrate and malate accumulation. The strain demonstrated strong NADH/NADPH production and acetyl-CoA synthase activity. Using insights and gene targets from this analysis, we engineered L. tetrasporus for malate production. The engineered strain produced 7.5 g/L malic acid (0.25 g/g yield) in shake flasks with glucose-acetate media and 28.8 g/L malic acid at a yield of 0.20 g/g in fed-batch mode with corn-stover hydrolysate. Together, these insights and rational strain engineering establish L. tetrasporus as a versatile, Crabtree-negative platform that is an energy-CO2-bioproduction nexus for channeling CO2 carbon into value-added bioproducts.

Xiao, Zhengyang↗
DOE OSTI2022

Integrated flux and pool size analysis in plant central metabolism reveals unique roles of glycine and serine during photorespiration

Photorespiration is an essential process juxtaposed between plant carbon and nitrogen metabolism that responds to dynamic environments. Photorespiration recycles inhibitory intermediates arising from oxygenation reactions catalysed by Rubisco back into the C 3 cycle, but it is unclear what proportions of its nitrogen-containing intermediates (glycine and serine) are exported into other metabolisms in vivo and how these pool sizes affect net CO 2 gas exchange during photorespiratory transients. Here, to address this uncertainty, we measured rates of amino acid export from photorespiration using isotopically non-stationary metabolic flux analysis. This analysis revealed that ~23–41% of the photorespiratory carbon was exported from the pathway as serine under various photorespiratory conditions. Furthermore, we determined that the build-up and relaxation of glycine pools constrained a large portion of photosynthetic acclimation during photorespiratory transients. Furthermore, these results reveal the unique and important roles of glycine and serine in successfully maintaining various photorespiratory fluxes that occur under environmental fluctuations in nature and providing carbon and nitrogen for metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Metabolic flux, metabolite, and transcript analysis uncover reprogramming of metabolism toward higher seed oil

Overexpression of WRINKLED1 (WRI1), a master regulator of glycolysis and fatty acid biosynthesis, together with DIACYLGLYCEROL ACYLTRANSFERASE1 (DGAT1), which catalyzes the final step of triacylglycerol assembly, is a promising strategy for enhancing seed oil content. However, how these regulators coordinate system-wide metabolic reprogramming at the levels of gene expression, metabolite pools, and fluxes remains poorly understood. To address this, we performed 13 C-metabolic flux analysis, metabolomics, and transcriptomics on in vitro cultured pennycress (Thlaspi arvense L.) embryos overexpressing the native WRI1 and DGAT1 homologs. Here, in cultured embryos, WRI1/DGAT1 overexpression increased triacylglycerol accumulation by 28% while reducing protein content by 34%, relative to the wild type. Embryos showed ∼20-fold and 50-fold upregulation of WRI1 and DGAT1 along with induction of WRI1 target genes in glycolysis and fatty acid biosynthesis. Genes associated with photosynthesis and Calvin cycle functions were also upregulated, whereas genes encoding ribosomal proteins and seed storage proteins were strongly repressed, consistent with the observed lipid–protein tradeoff. Flux analysis revealed that enhanced triacylglycerol biosynthesis is supported by increased flux through the Rubisco shunt and cytosolic pyruvate kinase, while the oxidative pentose phosphate pathway and malic enzyme contributed little to NADPH or pyruvate supply. Metabolomic profiling revealed extensive perturbations in glycolytic intermediates, tricarboxylic acid cycle metabolites, and amino acids. In plant grown seeds, WRI1/DGAT1 lines also showed a modest but significant increase in total lipid content. Collectively, these findings reveal how WRI1 and DGAT1 reprogram central metabolism to enhance oil accumulation, with relevance to mature seeds.

59 BASIC BIOLOGICAL 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

A High-Quality Genome-Scale Model for Rhodococcus opacus Metabolism

Rhodococcus opacus is a bacterium that has a high tolerance to aromatic compounds and can produce significant amounts of triacylglycerol (TAG). Here, we present iGR1773, the first genome-scale model (GSM) of R. opacus PD630 metabolism based on its genomic sequence and associated data. The model includes 1773 genes, 3025 reactions, and 1956 metabolites, was developed in a reproducible manner using CarveMe, and was evaluated through Metabolic Model tests (MEMOTE). We combine the model with two Constraint-Based Reconstruction and Analysis (COBRA) methods that use transcriptomics data to predict growth rates and fluxes: E-Flux2 and SPOT (Simplified Pearson Correlation with Transcriptomic data). Growth rates are best predicted by E-Flux2. Flux profiles are more accurately predicted by E-Flux2 than flux balance analysis (FBA) and parsimonious FBA (pFBA), when compared to 44 central carbon fluxes measured by 13C-Metabolic Flux Analysis (13C-MFA). Under glucose-fed conditions, E-Flux2 presents an R2 value of 0.54, while predictions based on pFBA had an inferior R2 of 0.28. We attribute this improved performance to the extra activity information provided by the transcriptomics data. For phenol-fed metabolism, in which the substrate first enters the TCA cycle, E-Flux2’s flux predictions display a high R2 of 0.96 while pFBA showed an R2 of 0.93. We also show that glucose metabolism and phenol metabolism function with similar relative ATP maintenance costs. These findings demonstrate that iGR1773 can help the metabolic engineering community predict aromatic substrate utilization patterns and perform computational strain design.

Roell, Garrett W.↗
DOE OSTI2026

Metabolic rewiring and biomass redistribution enable optimized mixotrophic growth in Chlamydomonas

Aquatic photosynthetic systems account for approximately one-half of all global carbon assimilation and could be a significant source of renewable fuels and feedstocks. However, rapid growth and biomass production in algae have not always translated into high product yields, partly because central metabolism is context specific, with metabolic fluxes being influenced by nutrient conditions and other environmental factors. In the green microalga Chlamydomonas reinhardtii (Chlamydomonas), mixotrophic cultures (acetate + light) grow far faster than phototrophic (light only) or heterotrophic (acetate + dark) cultures, even though acetate partially suppresses photosynthesis. Here, an isotopic dilution strategy with unlabeled acetate was combined with 13 CO 2 transient labeling to perform isotopically nonstationary metabolic flux analysis (INST-MFA) and to directly compare autotrophic and mixotrophic metabolism in Chlamydomonas supported by data from transcriptomics, proteomics, and metabolomics. INST-MFA indicated that acetate induces a synergistic rewiring of metabolism, conserving carbon by using the glyoxylate cycle and suppressing gluconeogenesis, the latter of which was discordant with omics results and prior models. Additionally, our data provide a plausible rationale for the well-known suppression of photosynthesis by acetate. We propose that reduced total protein content in mixotrophic versus phototrophic cells, much of which is attributed to reduced levels of photosynthetic proteins, decreases the costly metabolic burden of protein synthesis and represents a growth rate optimization strategy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

In Vivo Thermodynamic Analysis of Glycolysis in Clostridium thermocellum and Thermoanaerobacterium saccharolyticum Using 13 C and 2 H Tracers

Clostridium thermocellum and Thermoanaerobacterium saccharolyticum are thermophilic anaerobic bacteria with complementary metabolic capabilities that utilize distinct glycolytic pathways for the conversion of cellulosic sugars to biofuels. We integrated quantitative metabolomics with 2 H and 13 C metabolic flux analysis to investigate the in vivo reversibility and thermodynamics of the central metabolic networks of these two microbes. We found that the glycolytic pathway in C. thermocellum operates remarkably close to thermodynamic equilibrium, with an overall drop in Gibbs free energy 5-fold lower than that of T. saccharolyticum or anaerobically grown Escherichia coli . The limited thermodynamic driving force of glycolysis in C. thermocellum could be attributed in large part to the small free energy of the phosphofructokinase reaction producing fructose bisphosphate. The ethanol fermentation pathway was also substantially more reversible in C. thermocellummthan in T. saccharolyticum. These observations help explain the comparatively low ethanol titers of C. thermocellum and suggest engineering interventions that can be used to increase its ethanol productivity and glycolytic rate. In addition to thermodynamic analysis, we used our isotope tracer data to reconstruct the T. saccharolyticum central metabolic network, revealing exclusive use of the Embden-Meyerhof-Parnas (EMP) pathway for glycolysis, a bifurcated tricarboxylic acid (TCA) cycle, and a sedoheptulose bisphosphate bypass active within the pentose phosphate pathway. Thermodynamics constitutes a key determinant of flux and enzyme efficiency in metabolic networks. Here, we provide new insights into the divergent thermodynamics of the glycolytic pathways of C. thermocellum and T. saccharolyticum, two industrially relevant thermophilic bacteria whose metabolism still is not well understood. We report that while the glycolytic pathway in T. saccharolyticum is as thermodynamically favorable as that found in model organisms, such as E. coli or Saccharomyces cerevisiae , the glycolytic pathway of C. thermocellum operates near equilibrium. The use of a near-equilibrium glycolytic pathway, with potentially increased ATP yield, by this cellulolytic microbe may represent an evolutionary adaptation to growth on cellulose, but it has the drawback of being highly susceptible to product feedback inhibition. The results of this study will facilitate future engineering of high-performance strains capable of transforming cellulosic biomass to biofuels at high yields and titers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Re-Programing Glucose Catabolism in the Microalga Chlorella sorokiniana under Light Condition

The microalga Chlorella sorokiniana has attracted much attention for lipid production and wastewater treatment. It can perform photosynthesis and organic carbon utilization concurrently. To understand its phototrophic metabolism, a biomass compositional analysis, a 13 C metabolic flux analysis, and metabolite pool size analyses were performed. Under dark condition, the oxidative pentose phosphate pathway (OPP) was the major route for glucose catabolism (88% carbon flux) and a cyclic OPP–glycolytic route for glucose catabolism was formed. Under light condition, fluxes in the glucose catabolism, tricarboxylic acid (TCA) cycle, and anaplerotic reaction (CO 2 fixation via phosphoenolpyruvate carboxylase) were all suppressed. Meanwhile, the RuBisCO reaction became active and the ratio of its carbon fixation to glucose carbon utilization was determined as 7:100. Moreover, light condition significantly reduced the pool sizes of sugar phosphate metabolites (such as E4P, F6P, and S7P) and promoted biomass synthesis (which reached 0.155 h –1 ). In addition, light condition increased glucose consumption rates, leading to higher ATP and NADPH production and a higher protein content (43% vs. 30%) in the biomass during the exponential growth phase.

13C-MFA↗
DOE OSTI2020

Autotrophic and mixotrophic metabolism of an anammox bacterium revealed by in vivo 13 C and 2 H metabolic network mapping

Anaerobic ammonium-oxidizing (anammox) bacteria mediate a key step in the biogeochemical nitrogen cycle and have been applied worldwide for the energy-efficient removal of nitrogen from wastewater. However, outside their core energy metabolism, little is known about the metabolic networks driving anammox bacterial anabolism and use of different carbon and energy substrates beyond genome-based predictions. Here, we experimentally resolved the central carbon metabolism of the anammox bacterium Candidatus ‘Kuenenia stuttgartiensis’ using time-series 13 C and 2 H isotope tracing, metabolomics, and isotopically nonstationary metabolic flux analysis. Our findings confirm predicted metabolic pathways used for CO 2 fixation, central metabolism, and amino acid biosynthesis in K. stuttgartiensis , and reveal several instances where genomic predictions are not supported by in vivo metabolic fluxes. This includes the use of the oxidative branch of an incomplete tricarboxylic acid cycle for alpha-ketoglutarate biosynthesis, despite the genome not having an annotated citrate synthase. We also demonstrate that K. stuttgartiensis is able to directly assimilate extracellular formate via the Wood–Ljungdahl pathway instead of oxidizing it completely to CO 2 followed by reassimilation. In contrast, our data suggest that K. stuttgartiensis is not capable of using acetate as a carbon or energy source in situ and that acetate oxidation occurred via the metabolic activity of a low-abundance microorganism in the bioreactor’s side population. Together, these findings provide a foundation for understanding the carbon metabolism of anammox bacteria at a systems-level and will inform future studies aimed at elucidating factors governing their function and niche differentiation in natural and engineered ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Systematic engineering for production of anti-aging sunscreen compound in Pseudomonas putida

Sunscreen has been used for thousands of years to protect skin from ultraviolet radiation. However, the use of modern commercial sunscreen containing oxybenzone, ZnO, and TiO 2 has raised concerns due to their negative effects on human health and the environment. In this study, we aim to establish an efficient microbial platform for production of shinorine, a UV light absorbing compound with anti-aging properties. First, we methodically selected an appropriate host for shinorine production by analyzing central carbon flux distribution data from prior studies alongside predictions from genome-scale metabolic models (GEMs). We enhanced shinorine productivity through CRISPRi-mediated downregulation and utilized shotgun proteomics to pinpoint potential competing pathways. Simultaneously, we improved the shinorine biosynthetic pathway by refining its design, optimizing promoter usage, and altering the strength of ribosome binding sites. Finally, we conducted amino acid feeding experiments under various conditions to identify the key limiting factors in shinorine production. The study combines meta-analysis of 13 C-metabolic flux analysis, GEMs, synthetic biology, CRISPRi-mediated gene downregulation, and omics analysis to improve shinorine production, demonstrating the potential of Pseudomonas putida KT2440 as platform for shinorine production.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Accurate flux predictions using tissue-specific gene expression in plant metabolic modeling

The accurate prediction of complex phenotypes such as metabolic fluxes in living systems is a grand challenge for systems biology and central to efficiently identifying biotechnological interventions that can address pressing industrial needs. The application of gene expression data to improve the accuracy of metabolic flux predictions using mechanistic modeling methods such as flux balance analysis (FBA) has not been previously demonstrated in multi-tissue systems, despite their biotechnological importance. We hypothesized that a method for generating metabolic flux predictions informed by relative expression levels between tissues would improve prediction accuracy. Relative gene expression levels derived from multiple transcriptomic and proteomic datasets were integrated into FBA predictions of a multi-tissue, diel model of Arabidopsis thaliana’s central metabolism. This integration dramatically improved the agreement of flux predictions with experimentally based flux maps from 13 C metabolic flux analysis compared with a standard parsimonious FBA approach. Disagreement between FBA predictions and MFA flux maps was measured using weighted averaged percent error values, and for parsimonious FBA this was 169%–180% for high light conditions and 94%–103% for low light conditions, depending on the gene expression dataset used. This fell to 10%-13% and 9%-11% upon incorporating expression data into the modeling process, which also substantially altered the predicted carbon and energy economy of the plant.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

BayFlux v1.0

This is an open source software library for performing Bayesian genome scale 13C Metabolic Flux Analysis (Bayesian 13C MFA). This is a new type of metabolic modeling that quantifies flux in a probabilistic manner at the genome scale.

Garcia Martin, Hector↗
DOE OSTI2021

Computational Framework for Machine-Learning-Enabled 13 C Fluxomics

13 C metabolic flux analysis (MFA) has emerged as a powerful tool for synthetic biology. This optimization-based approach suffers long computation time and unstable solutions depending on the initial guess. Here, we develop a machine-learning-based framework for 13 C fluxomics. Specifically, training and test data sets are generated by metabolic network decomposition and flux sampling, in which flux ratios at metabolic nodes and simulated labeling patterns of metabolites are used as training targets and features, respectively. To improve prediction accuracy and simplify the model, automated processes are developed for flux ratio selection based on solvability and feature screening based on importance. We found that predictive performance can be significantly improved using both amino acids and central carbon metabolites in comparison with amino acids alone. Together with measured external fluxes, the predicted flux ratios determine the mass balance system, yielding global flux distributions. This approach is validated by flux estimation using both simulated and experimental data in comparison with canonical 13 C MFA. The approach represents a reliable fluxomics method readily applicable to high-throughput metabolic phenotyping, which highlights the advances of intelligent learning algorithms in synthetic biology, specifically in the Test and Learn stage of the Design-Build-Test-Learn cycle.

13C metabolic flux analysis↗
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↗
DOE OSTI2023

Progress in understanding and improving oil content and quality in seeds

The world’s population is projected to increase by two billion by 2050, resulting in food and energy insecurity. Oilseed crops have been identified as key to address these challenges: they produce and store lipids in the seeds as triacylglycerols that can serve as a source of food/feed, renewable fuels, and other industrially-relevant chemicals. Therefore, improving seed oil content and composition has generated immense interest. Research efforts aiming to unravel the regulatory pathways involved in fatty acid synthesis and to identify targets for metabolic engineering have made tremendous progress. This review provides a summary of the current knowledge of oil metabolism and discusses how photochemical activity and unconventional pathways can contribute to high carbon conversion efficiency in seeds. It also highlights the importance of 13 C-metabolic flux analysis as a tool to gain insights on the pathways that regulate oil biosynthesis in seeds. Finally, a list of key genes and regulators that have been recently targeted to enhance seed oil production are reviewed and additional possible targets in the metabolic pathways are proposed to achieve desirable oil content and quality.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

BayFlux: A Bayesian method to quantify metabolic Fluxes and their uncertainty at the genome scale

Metabolic fluxes, the number of metabolites traversing each biochemical reaction in a cell per unit time, are crucial for assessing and understanding cell function. 13 C Metabolic Flux Analysis ( 13 C MFA) is considered to be the gold standard for measuring metabolic fluxes. 13 C MFA typically works by leveraging extracellular exchange fluxes as well as data from 13 C labeling experiments to calculate the flux profile which best fit the data for a small, central carbon, metabolic model. However, the nonlinear nature of the 13 C MFA fitting procedure means that several flux profiles fit the experimental data within the experimental error, and traditional optimization methods offer only a partial or skewed picture, especially in “non-gaussian” situations where multiple very distinct flux regions fit the data equally well. Here, we present a method for flux space sampling through Bayesian inference (BayFlux), that identifies the full distribution of fluxes compatible with experimental data for a comprehensive genome-scale model. This Bayesian approach allows us to accurately quantify uncertainty in calculated fluxes. We also find that, surprisingly, the genome-scale model of metabolism produces narrower flux distributions (reduced uncertainty) than the small core metabolic models traditionally used in 13 C MFA. The different results for some reactions when using genome-scale models vs core metabolic models advise caution in assuming strong inferences from 13 C MFA since the results may depend significantly on the completeness of the model used. Based on BayFlux, we developed and evaluated novel methods (P- 13 C MOMA and P- 13 C ROOM) to predict the biological results of a gene knockout, that improve on the traditional MOMA and ROOM methods by quantifying prediction uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗