K-FIT: An accelerated kinetic parameterization algorithm using steady-state fluxomic data
Not provided.
SEARCH · Search NASA
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.
Not provided.
Quantitative understanding of clostridial metabolism is of longstanding interest due to the importance of Clostridia as model anaerobes, biotechnology workhorses, and contributors to evolutionary history and ecosystem. Current computational methods such as flux balance analysis-based construction of clostridial metabolism in genome scale provide a fundamental framework for metabolic analysis. However, this method alone is inadequate to characterize cellular metabolic activity. Experiment-driven approaches including isotope tracer-based fluxomics in association with genetic and biochemical methods are needed to gain a more comprehensive understanding. Here we focus on typical examples where these integrated approaches have contributed to the identification of new metabolic pathways and quantification of metabolic fluxes in Clostridia. Here, we also highlight the opportunities and challenges of cutting-edge fluxomics approaches such as machine learning modeling, deuterium tracer approach, and high throughput flux phenotyping in exploring clostridial metabolism with respect to inorganic carbon utilization, redox cofactor interconversion, and other key metabolic features.
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.
Diazotrophic bacteria are an attractive biological alternative to synthetic nitrogen fertilizers due to their remarkable capacity to fix atmospheric nitrogen gas to ammonium via nitrogenase enzymes. However, how diazotrophic bacteria tailor central carbon catabolism to accommodate the energy requirement for nitrogenase activity is largely unknown. In this study, we used Azotobacter vinelandii DJ and an ammonium excreting mutant, AV3 (ΔNifL), to investigate central carbon metabolism fluxes and central cell bioenergetics in response to ammonium availability and nitrogenase activity. Enabled by the powerful and reliable methodology of 13 C-metabolic flux analysis, we show that the respiratory TCA cycle is upregulated in association with increased nitrogenase activity and causes a monotonic decrease in specific growth rate. Whereas the activity of the glycolytic Entner-Doudoroff pathway is positively correlated with the cell growth rate. These new observations are formulated into a 13 C-metabolic flux model which further improves the understanding and interpretation of intracellular bioenergetics. This analysis leads to the conclusion that, under aerobic conditions, respiratory TCA metabolism is responsible for the supply of additional ATP and reducing equivalents required for elevated nitrogenase activity. This study provides a quantitative relationship between central carbon and nitrogen metabolism in an aerobic diazotroph for the first time.
Metabolic flux analysis represents an essential perspective to understand cellular physiology and offers quantitative information to guide pathway engineering. A valuable approach for experimental elucidation of metabolic flux is dynamic flux analysis, which estimates the relative or absolute flow rates through a series of metabolic intermediates in a given pathway. It is based on kinetic isotope labeling experiments, liquid chromatography-mass spectrometry (LC-MS), and computational analysis that relate kinetic isotope trajectories of metabolites to pathway activity. Herein, we illustrate the mathematic principles underlying the dynamic flux analysis and mainly focus on describing the experimental procedures for data generation. This protocol is exemplified using cyanobacterial metabolism as an example, for which reliable labeling data for central carbon metabolites can be acquired quantitatively. This protocol is applicable to other microbial systems as well and can be readily adapted to address different metabolic processes.
Not Available
To render the resulting biomass more attractive and amenable for utilization as the basis for low-carbon intensity bioproducts, single-celled algae need to be biochemically and metabolically poised to assimilate and store the delivered carbon in the fastest and most efficient manner. Accelerating biochemical carbon storage, as primarily carbohydrates or lipids, is critical to achieve the high carbon capture potential that is assigned to algae. To guide strain optimization and engineering for maximizing carbon capture and storage, it is essential to elucidate the link between carbon metabolism and biomass composition. Most published metabolomics work in algae remains largely restricted to ideal and simplified environmental conditions in model organisms, thereby limiting their translation to outdoor implementation. In this work, we utilize 13 C isotopic labeling to characterize distinct intracellular metabolic fluxes before, during, and after nitrogen depletion-induced compositional shifts in Scenedesmus obliquus UTEX 393. The results indicate that a transition to carbohydrates is characterized by diverting flux to starch instead of replenishing the Calvin cycle for CO 2 fixation whereas the subsequent transition to lipids is fueled by NADPH produced by upregulating the phosphoenolpyruvate carboxylase (PEPC)–malic enzyme (ME) cycle flux. Our work highlights bottlenecks to carbohydrate- and lipid-rich biomass and can guide implementable strategies to control the fate of fixed carbon in S. obliquus.
Cyanobacteria with the capability of oxygen-evolving photosynthesis play a vital role in the global carbon, oxygen, and nitrogen cycles. From a biochemical perspective, all phenotypic traits and functions of cyanobacteria are derived from the potential of their metabolism and mirrored in metabolic fluxes through an integrated metabolic network. Metabolic flux, the rate of turnover of molecules through a metabolic pathway, is therefore of great importance for an in-depth understanding of the role of cyanobacteria in the biosphere. Recent interest in cyanobacteria metabolic fluxes is also motivated by the prospect to utilize their light-driven CO2 fixation for sustainable production of valuable chemicals and biofuels. Valorization of cyanobacterial metabolism will require engineering of the native metabolic network - a challenge that can benefit from a quantitative understanding of metabolic reactions taking place intracellularly. In this chapter, we will discuss some intriguing questions relevant to metabolic fluxes in cyanobacteria: How are metabolic fluxes predicted or measured? What do we know about the fluxome inside various cyanobacteria thus far? How are they different between strains as well as in different conditions? And how can we use fluxomics information to guide metabolic engineering and synthetic biology? This chapter will highlight updated knowledge on cyanobacterial metabolism. Specifically, we will focus on providing readers a quantitative vision to understand how this photosynthetic model works.
Metabolomics and fluxomics are core approaches to directly profile and interrogate cellular metabolism in response to various genetic or environmental perturbations. In order to accurately measure the abundance and isotope enrichment of intracellular metabolites, cell culture samples must be rapidly harvested and cold-quenched to preserve the in vivo metabolic state of the cells at the time of sample collection. When dealing with suspension cultures, this process is complicated by the need to separate the liquid culture media from cellular biomass prior to metabolite extraction. Here, we examine the efficacy of several commonly used metabolic quenching methods, using the model cyanobacterium Synechocystis sp. PCC 6803 as an example. Multiple 13 C-labeled compounds, including 13 C-bicarbonate, 13 C-glucose and 13 C-glutamine, were used as tracers during the sample collection and cold-quenching process to assess the extent of metabolic turnover after cells were harvested from culture flasks. We show that the combination of rapid filtration followed by 100% cold (–80°C) methanol quenching exhibits the highest quenching efficiency, while mixing cell samples with a partially frozen 30% methanol slurry (–24°C) followed by centrifugation is slightly less effective at quenching metabolism but enables less laborious sample processing. By contrast, rapidly mixing the cells with a saline ice slurry (~0°C) is less effective as indicated by high isotope-labeling rates after sample harvest, while mixing the cells with 60% cold methanol (–65°C) prior to centrifugation causes significant metabolite loss. Furthermore, this study demonstrates a rigorous, quantitative, and broadly applicable method for assessing the metabolic quenching efficacy of protocols used for sample collection in metabolomics and fluxomics studies.
Biomass crops engineered to accumulate energy-dense triacylglycerols (TAG or ‘vegetable oils’) in their vegetative tissues have emerged as potential feedstocks to meet the growing demand for renewable diesel and sustainable aviation fuel (SAF). Unlike oil palm and oilseed crops, the current commercial sources of TAG, vegetative tissues, such as leaves and stems, only transiently accumulate TAG. In this report, we used grain (Texas430 or TX430) and sugar-accumulating ‘sweet’ (Ramada) genotypes of sorghum, a high-yielding, environmentally resilient biomass crop, to accumulate TAG in leaves and stems. We initially tested several gene combinations for a ‘push-pull-protect’ strategy. The top TAG-yielding constructs contained five oil transgenes for a sorghum WRINKLED1 transcription factor (‘push’), a Cuphea viscosissima diacylglycerol acyltransferase (DGAT; ‘pull’), a modified sesame oleosin (‘protect’) and two combinations of specialized Cuphea lysophosphatidic acid acyltransferases and medium-chain acyl-acyl carrier protein thioesterases. Though intended to generate oils with medium-chain fatty acids, engineered lines accumulated oleic acid-rich oil to amounts of up to 2.5% DW in leaves and 2.0% DW in stems in the greenhouse, 36-fold and 49-fold increases relative to wild-type (WT) plants, respectively. Under field conditions, the top-performing event accumulated TAG to amount to 5.5% DW in leaves and 3.5% DW in stems, 78-fold and 58-fold increases, respectively, relative to WT TX430. Transcriptomic and fluxomic analyses revealed potential bottlenecks for increased TAG accumulation. Overall, our studies highlight the utility of a lab-to-field pipeline coupled with systems biology studies to deliver high vegetative oil sorghum for SAF and renewable diesel production.
Carbon capture and biochemical storage are some of the primary drivers of photosynthetic yield and productivity. To elucidate the mechanisms governing carbon allocation, we designed a photosynthetic light response test system for genetic and metabolic carbon assimilation tracking, using microalgae as simplified plant models. The systems biology mapping of high light-responsive photophysiology and carbon utilization dynamics between two variants of the same Picochlorum celeri species, TG1 and TG2 elucidated metabolic bottlenecks and transport rates of intermediates using instationary 13 C-fluxomics. Simultaneous global gene expression dynamics showed 73% of the annotated genes responding within one hour, elucidating a singular, diel-responsive transcription factor, closely related to the CCA1/LHY clock genes in plants, with significantly altered expression in TG2. Transgenic P. celeri TG1 cells expressing the TG2 CCA1/LHY gene, showed 15% increase in growth rates and 25% increase in storage carbohydrate content, supporting a coordinating regulatory function for a single transcription factor.
Cyanobacterial carboxysomes encapsulate carbonic anhydrase and ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO). Genetic deletion of the major structural proteins encoded within the ccm operon in Synechococcus sp. PCC 7002 ( ΔccmKLMN ) disrupts carboxysome formation and significantly affects cellular physiology. In this study, we employed both metabolite pool size analysis and isotopically nonstationary metabolic flux analysis (INST-MFA) to examine metabolic regulation in cells lacking carboxysomes. Under high CO 2 environments (1%), the ΔccmKLMN mutant could recover growth and had a similar central flux distribution as the control strain, with the exceptions of moderately decreased photosynthesis and elevated biomass protein content and photorespiration activity. Metabolite analyses indicated that the ΔccmKLMN strain had significantly larger pool sizes of pyruvate (>18 folds), UDPG (uridine diphosphate glucose), and aspartate as well as higher levels of secreted organic acids (e.g., malate and succinate). These results suggest that the ΔccmKLMN mutant is able to largely maintain a fluxome similar to the control strain by changing in intracellular metabolite concentrations and metabolite overflows under optimal growth conditions. When CO 2 was insufficient (0.2%), provision of acetate moderately promoted mutant growth. Interestingly, the removal of microcompartments may loosen the flux network and promote RuBisCO side-reactions, facilitating redirection of central metabolites to competing pathways (i.e., pyruvate to heterologous lactate production). This study provides important insights into metabolic regulation via enzyme compartmentation and cyanobacterial compensatory responses.
Soil Pseudomonas species, which thrive on lignin derivatives, are widely explored for biotechnology applications in lignin valorization. However, how the native metabolism coordinates phenolic carbon processing with required cofactor generation remains poorly understood. Here, we achieve quantitative understanding of this metabolic balance through a detailed multi-omics investigation of Pseudomonas putida KT2440 grown on four common phenolic acid substrates: ferulate, p-coumarate, vanillate, and 4-hydroxybenzoate. Relative to succinate, proteomics reveals > 140-fold increase in transport and catabolic proteins for aromatics, but metabolomics identifies bottlenecks in initial catabolism to maintain favorable cellular energy charge, which is compromised in mutants with resolved bottlenecks. Up to 30-fold increase in pyruvate carboxylase and glyoxylate shunt proteins implies a metabolic remodeling confirmed by kinetic 13 C-metabolomics. Quantitative analysis by 13 C-fluxomics demonstrates coupling of this remodeling with cofactor production. Specifically, anaplerotic carbon recycling through pyruvate carboxylase promotes tricarboxylic acid cycle fluxes to generate 50-60% NADPH yield and 60-80% NADH yield, resulting in up to 6-fold greater ATP surplus than with succinate metabolism; the glyoxylate shunt sustains cataplerotic flux through malic enzyme for the remaining NADPH yield. This quantitative blueprint affords cofactor imbalance predictions in proposed engineering of key metabolic nodes in lignin valorization pathways.
Summary Biomass crops engineered to accumulate energy‐dense triacylglycerols (TAG or ‘vegetable oils’) in their vegetative tissues have emerged as potential feedstocks to meet the growing demand for renewable diesel and sustainable aviation fuel (SAF). Unlike oil palm and oilseed crops, the current commercial sources of TAG, vegetative tissues, such as leaves and stems, only transiently accumulate TAG. In this report, we used grain (Texas430 or TX430) and sugar‐accumulating ‘sweet’ (Ramada) genotypes of sorghum, a high‐yielding, environmentally resilient biomass crop, to accumulate TAG in leaves and stems. We initially tested several gene combinations for a ‘push‐pull‐protect’ strategy. The top TAG‐yielding constructs contained five oil transgenes for a sorghum WRINKLED1 transcription factor (‘push’), a Cuphea viscosissima diacylglycerol acyltransferase (DGAT; ‘pull’), a modified sesame oleosin (‘protect’) and two combinations of specialized Cuphea lysophosphatidic acid acyltransferases and medium‐chain acyl‐acyl carrier protein thioesterases. Though intended to generate oils with medium‐chain fatty acids, engineered lines accumulated oleic acid‐rich oil to amounts of up to 2.5% DW in leaves and 2.0% DW in stems in the greenhouse, 36‐fold and 49‐fold increases relative to wild‐type (WT) plants, respectively. Under field conditions, the top‐performing event accumulated TAG to amount to 5.5% DW in leaves and 3.5% DW in stems, 78‐fold and 58‐fold increases, respectively, relative to WT TX430. Transcriptomic and fluxomic analyses revealed potential bottlenecks for increased TAG accumulation. Overall, our studies highlight the utility of a lab‐to‐field pipeline coupled with systems biology studies to deliver high vegetative oil sorghum for SAF and renewable diesel production.
Abstract Background Quantification of individual species in microbial co-cultures and consortia is critical to understanding and designing communities with prescribed functions. However, it is difficult to physically separate species or measure species-specific attributes in most multi-species systems. Anaerobic gut fungi (AGF) (Neocallimastigomycetes) are native to the rumen of large herbivores, where they exist as minority members among a wealth of prokaryotes. AGF have significant biotechnological potential owing to their diverse repertoire of potent lignocellulose-degrading carbohydrate-active enzymes (CAZymes), which indirectly bolsters activity of other rumen microbes through metabolic exchange. While decades of literature suggest that polysaccharide degradation and AGF growth are accelerated in co-culture with prokaryotes, particularly methanogens, methods have not been available to measure concentrations of individual species in co-culture. New methods to disentangle the contributions of AGF and rumen prokaryotes are sorely needed to calculate AGF growth rates and metabolic fluxes to prove this hypothesis and understand its causality for predictable co-culture design. Results We present a simple, microplate-based method to measure AGF and methanogen concentrations in co-culture based on fluorescence and absorbance spectroscopies. Using samples of < 2% of the co-culture volume, we demonstrate significant increases in AGF growth rate and xylan and glucose degradation rates in co-culture with methanogens relative to mono-culture. Further, we calculate significant differences in AGF metabolic fluxes in co-culture relative to mono-culture, namely increased flux through the energy-generating hydrogenosome organelle. While calculated fluxes highlight uncertainties in AGF primary metabolism that preclude definitive explanations for this shift, our method will enable steady-state fluxomic experiments to probe AGF metabolism in greater detail. Conclusions The method we present to measure AGF and methanogen concentrations enables direct growth measurements and calculation of metabolic fluxes in co-culture. These metrics are critical to develop a quantitative understanding of interwoven rumen metabolism, as well as the impact of co-culture on polysaccharide degradation and metabolite production. The framework presented here can inspire new methods to probe systems beyond AGF and methanogens. Simple modifications to the method will likely extend its utility to co-cultures with more than two organisms or those grown on solid substrates to facilitate the design and deployment of microbial communities for bioproduction and beyond.
Kinetic models of metabolic networks offer the promise of quantitative phenotype prediction. The mechanistic characterization of enzyme catalyzed reactions allows for tracing the effect of perturbations in metabolite concentrations and reaction fluxes in response to genetic and environmental perturbation that are beyond the scope of stoichiometric models. In this study, we develop a two-step computational pipeline for the rapid parameterization of kinetic models of metabolic networks using a curated metabolic model and available 13 C-labeling distributions under multiple genetic and environmental perturbations. The first step involves the elucidation of all intracellular fluxes in a core model of E. coli containing 74 reactions and 61 metabolites using 13 C-Metabolic Flux Analysis ( 13 C-MFA). Here, fluxes corresponding to the mid-exponential growth phase are elucidated for seven single gene deletion mutants from upper glycolysis, pentose phosphate pathway and the Entner-Doudoroff pathway. The computed flux ranges are then used to parameterize the same (i.e., k-ecoli74) core kinetic model for E. coli with 55 substrate-level regulations using the newly developed K-FIT parameterization algorithm. The K-FIT algorithm employs a combination of equation decomposition and iterative solution techniques to evaluate steady-state fluxes in response to genetic perturbations. k-ecoli74 predicted 86% of flux values for strains used during fitting within a single standard deviation of 13 C-MFA estimated values. By performing both tasks using the same network, errors associated with lack of congruity between the two networks are avoided, allowing for seamless integration of data with model building. Product yield predictions and comparison with previously developed kinetic models indicate shifts in flux ranges and the presence or absence of mutant strains delivering flux towards pathways of interest from training data significantly impact predictive capabilities. Using this workflow, the impact of completeness of fluxomic datasets and the importance of specific genetic perturbations on uncertainties in kinetic parameter estimation are evaluated.
NREL is developing robust microbial strains and scalable fermentation technologies for converting biomass sugars to fuels and biochemicals, and one-carbon (C1) gaseous substrates (e.g., CO2, CO, CH4) to high-value products. Core capabilities include: (1) Fermentation optimization of aerobic, microaerophilic, and anaerobic cultivation processes of gas and traditional liquid fermentation R&D; (2) Gas mass transfer R&D and CFD modeling for improved reactor design and gas mixing (3) Metabolic engineering in CRISPR editing tools to improve titer, rates, and yields; (4) Probing and optimizing carbon, redox, and energy flux using transcriptomics and fluxomic techniques; and (5) Consolidated bioprocessing (CBP) using cellulose and hemicellulose directly to minimize pretreatment.
Disclosed herein is Arrowland, a web-based software tool for inputting, managing and viewing multiomics data, such as transcriptomics, proteomics, metabolomics and fluxomics data in an interactive, intuitive and multiscale system.