Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “METABOLISM”

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 109 records · Page 6

Metabolic Cost of Experimental Exercises

Although the type and duration of activity during decompression was well documented, the metabolic cost of 1665 subject-exposures with 8 activity profiles from 17 altitude decompression sickness (DCS) protocols at Brooks City-Base, TX from 1983-2005 was not determined. Female and male human volunteers (30 planned, 4 completed) performed activity profiles matching those 8 activity profiles at ground level with continuous monitoring of metabolic cost. A Cosmed K4b2 Cardio Pulmonary Exercise Testing device was used to measure oxygen uptake (VO2) during the profiles. The results show levels of metabolic cost to the females for the profiles tested varied from 4.3 to 25.5 ml/kg/min and from 3.0 to 12.0 ml/kg/min to the males. The increase in VO2 from seated rest to the most strenuous of the 8 activity profiles was 3.6-fold for the females and 2.8-fold for the males. These preliminary data on 4 subjects indicate close agreement of oxygen uptake for activity performed during many subject-exposures as published earlier. The relatively low average oxygen uptake required to perform the most strenuous activity may imply the need for adjustment of modeling efforts using metabolic cost as a risk factor. Better definition of metabolic cost during exposure to altitude, a critical factor in DCS risk, may allow refinement of DCS prediction models.

Webb, James T.↗

Radiation Exposure Alters Expression of Metabolic Enzyme Genes in Mice

Most administered pharmaceuticals are metabolized by the liver. The health of the liver, especially the rate of its metabolic enzymes, determines the concentration of circulating drugs as well as the duration of their efficacy. Most pharmaceuticals are metabolized by the liver, and clinically-used medication doses are given with normal liver function in mind. A drug overdose can result in the case of a liver that is damaged and removing pharmaceuticals from the circulation at a rate slower than normal. Alternatively, if liver function is elevated and removing drugs from the system more quickly than usual, it would be as if too little drug had been given for effective treatment. Because of the importance of the liver in drug metabolism, we want to understand the effects of spaceflight on the enzymes of the liver and exposure to cosmic radiation is one aspect of spaceflight that can be modeled in ground experiments. Additionally, it has been previous noted that pre-exposure to small radiation doses seems to confer protection against later and larger radiation doses. This protective power of pre-exposure has been called a priming effect or radioadaptation. This study is an effort to examine the drug metabolizing effects of radioadaptation mechanisms that may be triggered by early exposure to low radiation doses.

Wotring, V. E.↗

Advantage of Animal Models with Metabolic Flexibility for Space Research Beyond Low Earth Orbit

As the world's space agencies and commercial entities continue to expand beyond Low Earth Orbit (LEO), novel approaches to carry out biomedical experiments with animals are required to address the challenge of adaptation to space flight and new planetary environments. The extended time and distance of space travel along with reduced involvement of Earth-based mission support increases the cumulative impact of the risks encountered in space. To respond to these challenges, it becomes increasingly important to develop the capability to manage an organism's self-regulatory control system, which would enable survival in extraterrestrial environments. To significantly reduce the risk to animals on future long duration space missions, we propose the use of metabolically flexible animal models as "pathfinders," which are capable of tolerating the environmental extremes exhibited in spaceflight, including altered gravity, exposure to space radiation, chemically reactive planetary environments and temperature extremes. In this report we survey several of the pivotal metabolic flexibility studies and discuss the importance of utilizing animal models with metabolic flexibility with particular attention given to the ability to suppress the organism's metabolism in spaceflight experiments beyond LEO. The presented analysis demonstrates the adjuvant benefits of these factors to minimize damage caused by exposure to spaceflight and extreme planetary environments. Examples of microorganisms and animal models with dormancy capabilities suitable for space research are considered in the context of their survivability under hostile or deadly environments outside of Earth. Potential steps toward implementation of metabolic control technology in spaceflight architecture and its benefits for animal experiments and manned space exploration missions are discussed.

metabilic flexibility↗

Modeling Electrolytic O2 Recovery from Metabolic CO2 for Advanced Closed Loop Life Support Systems in Extraterrestrial Human Missions

The International Space Station (ISS) is currently equipped with a complex, heavy, and power consuming system that recovers approximately 50% of O2 from metabolic CO2. Future long duration human missions to the Moon and Mars will necessitate a sustainable and highly efficient metabolic oxygen recovery system capable of yielding a minimum of 75% O2 recovery. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to significantly increase current metabolic O2 recovery efficiency, expand mission sustainability, and reduce complexity of the system. The novel design combines CO2 conversion to O2 along with C2H4 as byproduct and water electrolysis (currently conducted in two separate units) into a single compact unit that runs at standard conditions and is theoretically capable of generating O2 with a theoretical maximum metabolic CO2 conversion of 73% while consuming less than metabolic water. This paper presents a comprehensive multi-physic 3D model developed at MSFC on CO2 conversion to O2 and C2H4 at standard conditions via MFECR. The 3D spatial domain of the model is a replica of the actual MFECR’s 3D drawing generated for the MFECR fabrication and operated to recover O2 from CO2 yielding C2H4 as byproduct. Electrochemical (EC) physics that includes EC multicomponent reaction mechanisms, mass transport, and electrical current density distributions is coupled in the model with all the other physics phenomena involved in the MFECR’s process, such as two-phase flow, free and porous fluid regimes, multicomponent mass transfer, heat transfer, and DC electrical current generation along with Joule heating effect. The EC reaction sections of the MFECR consists of two porous gas diffusion electrodes (GDE) and an electrolyte serpentine channel sandwiched in the middle. The CO2 feeds the cathode serpentine channel and part of the O2 product is fed back to the anode serpentine chamber. An alkaline solution feeds the electrolyte serpentine chamber wetting the GDEs of both, the anode and cathode allowing the OH- ionic transport between them. The EC reactions in the cathode’s GDE yield C2H4 from CO2 and H2 from water while the EC reaction in the anode’s GDE yields O2. The authors will present in this paper the validation of the model using experimental data and the utilization of the validated model in building a reliable simulator that will not only assist the authors on the MFECR design but also the optimization of its operation in the ISS and future spatial human missions.

Jesus A Dominguez↗

Probing interspecies metabolic interactions within a synthetic binary microbiome using genome-scale modeling

Metabolic interactions within a microbial community play a key role in determining the structure, function, and composition of the community. However, due to the complexity and intractability of natural microbiomes, limited knowledge is available on interspecies interactions within a community. In this work, using a binary synthetic microbiome, a methanotroph-photoautotroph (M-P) coculture, as the model system, we examined different genome-scale metabolic modeling (GEM) approaches to gain a better understanding of the metabolic interactions within the coculture, how they contribute to the enhanced growth observed in the coculture, and how they evolve over time. Using batch growth data of the model M-P coculture, we compared three GEM approaches for microbial communities. Two of the methods are existing approaches: SteadyCom, a steady state GEM, and dynamic flux balance analysis (DFBA) Lab, a dynamic GEM. We also proposed an improved dynamic GEM approach, DynamiCom, for the M-P coculture. SteadyCom can predict the metabolic interactions within the coculture but not their dynamic evolutions; DFBA Lab can predict the dynamics of the coculture but cannot identify interspecies interactions. DynamiCom was able to identify the cross-fed metabolite within the coculture, as well as predict the evolution of the interspecies interactions over time. A new dynamic GEM approach, DynamiCom, was developed for a model M-P coculture. Constrained by the predictions from a validated kinetic model, DynamiCom consistently predicted the top metabolites being exchanged in the M-P coculture, as well as the establishment of the mutualistic N-exchange between the methanotroph and cyanobacteria. The interspecies interactions and their dynamic evolution predicted by DynamiCom are supported by ample evidence in the literature on methanotroph, cyanobacteria, and other cyanobacteria-heterotroph cocultures.

59 BASIC BIOLOGICAL SCIENCES↗

Integrated lipidomic and proteomic profiling reveals metabolic network disruption by SARS-CoV-2 variants

The rapid evolution of SARS-CoV-2 has produced myriad viral strains with increasing transmissibility and capacity for immune evasion. While effective vaccination campaigns have reduced the fatalities associated with SARS-CoV-2, infections continue, and a detailed understanding of how this virus manipulates host biochemical pathways remains elusive. We asked both whether the patterns of host lipid rewiring remained consistent across variants and whether the changes in the abundance of lipid classes are related to changes in the expression of the enzymes involved in their biosynthesis. We compared global nontargeted lipidomics on A549-ACE2 cells infected with the delta variant (B.1.617.2), or the omicron (B.1.1.529) variant to our previous results of global nontargeted lipidomics on A549-ACE2 cells infected with the original WA1 strain and further performed quantitative proteomics to assess changes in the host proteome. We found that metabolic rewiring, both on the lipid and the enzymatic level, is remarkably consistent across all three variants. We further mapped changes in the expression of host metabolic enzymes, linking enzyme expression to alterations in the abundance of specific lipids during infection. This analysis identified key proteins related to virus-mediated changes in lipid abundance, including fatty acid synthase (FASN), lysosomal acid lipase (LIPA), and ORMDL, a regulator of sphingolipid biosynthesis. These integrated lipidomic and proteomic experiments shed light on the importance of the complex network of host metabolism networks that support SARS-CoV-2 infection and suggest that lipid metabolism may be a promising avenue for uncovering conserved therapeutic targets.

SARS-CoV-2↗

Adaptive laboratory evolution and metabolic engineering of Cupriavidus necator for improved catabolism of volatile fatty acids

Bioconversion of high-volume waste streams into value-added products will be an integral component of the growing bioeconomy. Volatile fatty acids (VFAs) (e.g., butyrate, valerate, and hexanoate) are an emerging and promising waste-derived feedstock for microbial carbon upcycling. Cupriavidus necator H16 is a favorable host for conversion of VFAs into various bioproducts due to its diverse carbon metabolism, ease of metabolic engineering, and use at industrial scales. Here, in this study, we report that a common strategy to improve product titers in C. necator, deletion of the polyhydroxybutyrate (PHB) biosynthetic operon, results in a significant growth defect on VFA substrates. Using adaptive laboratory evolution, we identify mutations to the regulator gene phaR, the two-component response regulator-histidine kinase pair encoded by H16_A1372/H16_A1373, and the tripartite transporter assembly encoded by H16_A2296-A2298 as causative for improved growth on VFA substrates. Deletion of phaR and H16_A1373 led to significantly reduced NADH abundance accompanied by large changes to expression of genes involved in carbon metabolism, balance of electron carriers, and oxidative stress tolerance that may be responsible for improved growth of these engineered strains. These results provide insight into the role of PHB biosynthesis in carbon and energy metabolism and highlight a key role for the regulator PhaR in global regulatory networks. By combining mutations, we generated platform strains with significant growth improvements on VFAs, which can enable improved conversion of waste-derived VFA substrates to target bioproducts.

09 BIOMASS FUELS↗

Cell-Free Systems Biology: Characterizing Central Metabolism of Clostridium thermocellum with a Three-Enzyme Cascade Reaction

Genetic approaches have been traditionally used to understand microbial metabolism, but this process can be slow in nonmodel organisms due to limited genetic tools. An alternative approach is to study metabolism directly in the cell lysate. This avoids the need for genetic tools and is routinely used to study individual enzymatic reactions but is not generally used to study systems-level properties of metabolism. Here we demonstrate a new approach that we call “cell-free systems biology”, where we use well-characterized enzymes and multienzyme cascades to serve as sources or sinks of intermediate metabolites. This allows us to isolate subnetworks within metabolism and study their systems-level properties. To demonstrate this, we worked with a threeenzyme cascade reaction that converts pyruvate to 2,3-butanediol. Although it has been previously used in cell-free systems, its pH dependence was not well characterized, limiting its utility as a sink for pyruvate. We showed that improved proton accounting allowed better prediction of pH changes and that active pH control allowed 2,3-butanediol titers of up to 2.1 M (189 g/L) from acetoin and 1.6 M (144 g/L) from pyruvate. The improved proton accounting provided a crucial insight that preventing the escape of CO 2 from the system largely eliminated the need for active pH control, dramatically simplifying our experimental setup. We then used this cascade reaction to understand limits to product formation in Clostridium thermocellum, an organism with potential applications for cellulosic biofuel production. We showed that the fate of pyruvate is largely controlled by electron availability and that reactions upstream of pyruvate limit overall product formation.

09 BIOMASS FUELS↗

Developmentally-specific physiological and metabolic responses support drought resilience in switchgrass and constrains biofuel yield

Switchgrass (Panicum virgatum) is a promising bioenergy crop due in part to its resilience to drought stress. However, the significance of drought timing remains poorly understood, both from a plant biology perspective and its impact on downstream biofuel production. This study determines the developmental stage-specific physiological and metabolic responses of switchgrass to drought stress and its implications for biofuel production using a custom-built programmable irrigation system. Vegetative, flowering, and senescence-stage drought significantly reduced carbon dioxide assimilation, and stomatal conductance without affecting biomass yield. Metabolic profiling revealed significant accumulation of glucose, fructose, quinic acid, shikimate and GABA during vegetative-stage drought, while flowering and senescence stages exhibited limited metabolic changes. Similarly, specialized metabolites also displayed distinct developmental patterns, with vegetative-stage drought driving the most pronounced metabolic alterations. Thermochemically-treated and hydrolyzed switchgrass biomass from vegetative-stage drought showed elevated lignocellulose-derived compounds and saponins with the latter most positively correlating with fermentation lag times. Conversely, senescence-stage drought enhanced ethanol yields while lowering saponin levels in the hydrolysates. While vegetative-stage drought enhanced physiological resilience, it compromises downstream biofuel production by introducing fermentation inhibitors, particularly saponins.

biofuel↗

Understanding the dynamic nature of plant lipid anabolic and catabolic metabolism is key to sustainable oilseed engineering

Plant-derived oils are essential sources of reduced carbon and various fatty acid (FA) structures for food, biofuels, and the oleochemical industry. Despite extensive efforts, engineering mainstream oilseed crops to produce high levels of industrially valuable unusual FAs (UFAs) remains challenging. This review synthesizes recent advances in the understanding of lipid metabolic networks, emphasizing how species-specific regulation of FA synthesis, activation, and delivery influences triacylglycerol (TAG) assembly to govern the efficiency of UFA accumulation. Key insights reveal that acyl flux through anabolic and catabolic branches of lipid metabolism is tightly controlled by enzyme substrate selectivities, diacylglycerol (DAG) pool compartmentalization, and metabolic context, including lipid remodeling and degradation pathways. Engineering success is often constrained by incompatibilities between UFA biosynthetic enzymes and endogenous host metabolism, leading to flux imbalances, futile cycles, and undesired phenotypes. We highlight emerging strategies to overcome these barriers, such as the use of UFA-selective acyltransferases, coordinated manipulation of DAG source pools, suppression of competing endogenous enzymes, and exploitation of TAG remodeling mechanisms. This integrated synthesis provides a conceptual framework for logic-based engineering of oilseeds with enhanced UFA content by offering new avenues for sustainable biomanufacturing of valuable lipids.

acyltransferase specificity↗

Exercise alters molecular profiles of inflammation and substrate metabolism in human white adipose tissue

White adipose tissue (WAT) plays a significant role in whole body energy homeostasis, and its excess typifies obesity. In addition to WAT quantity, perturbations in the basic cellular processes of WAT (i.e., quality) are also associated with obesity and metabolic disease. Exercise training alleviates metabolic perturbations associated with obesity; however, the underlying molecular mechanisms that drive these metabolic adaptations in WAT are not well described. For this work, abdominal subcutaneous WAT biopsies were collected after an acute bout of exercise (1 day after) at baseline and following 3 wk of supervised aerobic training in sedentary overweight women (n = 6) without alterations in body weight and fat mass. RNA-seq, global proteomics, and phosphoproteomics in WAT revealed training-induced changes in 1,527 transcripts, 154 proteins, and 144 phosphosites, respectively. Training decreased abundance of transcripts and proteins involved in inflammation and components of the extracellular matrix and increased abundance of transcripts and proteins related to fatty acid esterification and lipolysis. In summary, short-term aerobic training significantly reduces local inflammation and increases lipid metabolism in WAT of sedentary overweight women—independent of alterations in body and fat mass. As such, some of the health benefits of aerobic training may occur through molecular alterations in WAT (i.e., enhanced quality) rather than a sheer reduction in WAT quantity.

60 APPLIED LIFE SCIENCES↗

Walker Branch Watershed: Daily Stream Metabolism and Organic Carbon Spiraling Metrics in the West Fork of Walker Branch, Tennessee, USA, 2004-2010

This dataset contains daily metabolism estimates of gross primary production (GPP), ecosystem respiration (ER), and net ecosystem production (NEP), in addition to organic carbon spiraling length (SOC) and mineralization velocity (VfOC) estimates at the West Fork of Walker Branch, a small headwater stream, in the Walker Branch Watershed, Tennessee, USA. Observations were made from 2004-2010 (2004-01-01 to 2010-12-31). These data were generated to assess seasonal and interannual variability in metabolism and organic carbon spiraling and to explore potential driver variables, as analyses of intra- and interannual variability in metabolism and organic carbon spiraling are currently limited, leaving knowledge gaps in the driving mechanisms of and future changes to stream metabolism and carbon processing under climate change. Additionally, measurements of discharge (Q), stream width, stream- and canopy-level photosynthetically active radiation (PAR), water temperature, and precipitation from this time frame are included. This dataset contains one data file in comma separated (*.csv) format.

54 ENVIRONMENTAL SCIENCES↗

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↗

Metabolic skinflint or spendthrift? Insights into ground sloth integument and thermophysiology revealed by biophysical modeling and clumped isotope paleothermometry

Abstract Remains of megatheres have been known since the 18th -century and were among the first megafaunal vertebrates to be studied. While several examples of preserved integument show a thick coverage of fur for smaller ground sloths living in cold climates such as Mylodon and Nothrotheriops , comparatively very little is known about megathere skin. Assuming a typical placental mammal metabolism, it was previously hypothesized that megatheres would have had little-to-no fur as they achieved giant body sizes. Here the “hairless model of integument” is tested using geochemical analyses to estimate body temperature to generate novel models of ground sloth metabolism, fur coverage, and paleoclimate with Niche Mapper software. The simulations assuming metabolic activity akin to those of modern xenarthrans suggest that sparse fur coverage would have resulted in cold stress across most latitudinal ranges inhabited by extinct ground sloths. Specifically, Eremotherium predominantly required dense 10 mm fur with implications for seasonal changes of coat depth in northernmost latitudes and sparse fur in the tropics; Megatherium required dense 30 mm fur year-round in its exclusive range of cooler, drier climates; Mylodon and Nothrotheriops required dense 10–50 mm fur to avoid thermal stress, matching the integument remains of both genera, and further implying the use of behavioral thermoregulation. Moreover, clumped isotope paleothermometry data from the preserved teeth of four genera of ground sloth yielded reconstructed body temperatures lower than those previously reported for large terrestrial mammals (29 ± 2°–32 ± 3° C). This combination of low metabolisms and thick fur allowed ground sloths to inhabit various environments.

Deak, Michael D.↗

Alkaloids are associated with increased microbial diversity and metabolic function in poison frogs

Shifts in host-associated microbiomes can impact both host and microbes. It is of interest to understand how perturbations, like the introduction of exogenous chemicals, impact microbiomes. In poison frogs (family Dendrobatidae), the skin microbiome is exposed to alkaloids that the frogs sequester for defense. These alkaloids are antimicrobial; however, their effect on the frogs’ skin microbiome is unknown. To test this, we characterized microbial communities from field-collected dendrobatid frogs. Then, we conducted a laboratory experiment to monitor the effect of the alkaloid decahydroquinoline (DHQ) on the microbiome of two frog species with contrasting alkaloid loads in nature. In both datasets, we found that alkaloid-exposed microbiomes were more phylogenetically diverse, with an increase in diversity among rare taxa. Further, to better understand the isolate-specific response to alkaloids, we cultured microbial isolates from poison frog skin and found that many isolates exhibited enhanced growth or were not impacted by the addition of DHQ. To further explore the microbial response to alkaloids, we sequenced the metagenomes from high- and low-alkaloid frogs and observed a greater diversity of genes associated with nitrogen and carbon metabolism in high-alkaloid frogs. From these data, we hypothesized that some strains may metabolize the alkaloids. We used stable isotope tracing coupled to nanoSIMS (nanoscale secondary ion mass spectrometry), which supported the idea that some of these isolates are able to metabolize DHQ. Together, these data suggest that poison frog alkaloids open new niches for skin-associated microbes with specific adaptations, such as alkaloid metabolism, that enable survival in this environment.

59 BASIC BIOLOGICAL SCIENCES↗

A physiologically based pharmacokinetic (PBPK) model to align dosimetry of the isobutyl metabolic series in rats and humans

Here, we developed a physiologically based pharmacokinetic (PBPK) model in rats and humans for the isobutyl metabolic series including isobutyl acetate, isobutanol, isobutyraldehyde, and isobutyric acid. Chemical manufactures routinely use these compounds as solvents, for chemical synthesis, as potential biofuels, de-icing fluids, and additives for food and/or fragrance in consumer products. Human exposure to isobutyl compounds can occur through inhalation or oral routes. We previously developed a PBPK model for the propyl metabolic series and utilized it as a framework to create the isobutyl PBPK model due to the chemical similarities between the two series. To support model development, we measured in vitro metabolism of isobutyl acetate in rat and human blood and liver S9 fractions. Compared to rats, humans demonstrated faster isobutyl acetate hydrolysis in liver S9 fractions, while the hydrolysis rates in blood were similar between the two species. We used concentrations of isobutyl compounds measured in air and blood from rats exposed to isobutyl acetate and isobutanol as well as other published data to further parameterize the model. Following exposure to either isobutyl acetate or isobutanol, we observed isobutanol concentrations highest among the isobutyl compounds in the blood of rats. Overall, the model accurately predicts measured time course concentrations of isobutyl acetate, isobutanol, and isobutyric acid in available data in rats and humans. Sensitivity analyses identified alveolar ventilation rates, isobutyl metabolism rates, and cardiac output as the most sensitive parameters affecting concentrations of isobutyl compounds in blood. The isobutyl PBPK model enables comparisons of internal dose metrics across various isobutyl compound exposures and species and allows for calculation of equivalent external exposures that result in the same dose metric. Regulators can employ this PBPK model to predict and align internal dose metrics of isobutyl compounds for risk assessment purposes.

2-methyl-1-propanol↗

Metabolic complexity drives divergence in microbial communities

Microbial communities are shaped by environmental metabolites, but the principles that govern whether different communities will converge or diverge in any given condition remain unknown, posing fundamental questions about the feasibility of microbiome engineering. Here, in this work, we studied the longitudinal assembly dynamics of a set of natural microbial communities grown in laboratory conditions of increasing metabolic complexity. We found that different microbial communities tend to become similar to each other when grown in metabolically simple conditions, but they diverge in composition as the metabolic complexity of the environment increases, a phenomenon we refer to as the divergence-complexity effect. A comparative analysis of these communities revealed that this divergence is driven by community diversity and by the assortment of specialist taxa capable of degrading complex metabolites. An ecological model of community dynamics indicates that the hierarchical structure of metabolism itself, where complex molecules are enzymatically degraded into progressively simpler ones that then participate in cross-feeding between community members, is necessary and sufficient to recapitulate our experimental observations. In addition to helping understand the role of the environment in community assembly, the divergence-complexity effect can provide insight into which environments support multiple community states, enabling the search for desired ecosystem functions towards microbiome engineering applications.

59 BASIC BIOLOGICAL SCIENCES↗

Coupling flux balance analysis with reactive transport modeling through machine learning for rapid and stable simulation of microbial metabolic switching

Integrating genome-scale metabolic networks with reactive transport models (RTMs) provides a detailed description of the dynamic changes in microbial growth and metabolism. Despite promising demonstrations in the past, computational inefficiency has been pointed out as a critical issue to overcome because it requires repeated application of linear programming (LP) to obtain flux balance analysis (FBA) solutions in every time step and spatial grid. To address this challenge, we propose a new simulation method where we train and validate artificial neural networks (ANNs) using randomly sampled FBA solutions and incorporate the resulting surrogate FBA model (represented as algebraic equations) into RTMs as source/sink terms. We demonstrate the efficiency of our method via a case study of Shewanella oneidensis MR-1. During aerobic growth on lactate, S. oneidensis produces metabolic byproducts (such as pyruvate and acetate), which are subsequently consumed as alternative carbon sources when the preferred nutrients are depleted. To effectively simulate these complex dynamics, we used a cybernetic approach that models metabolic switches as the outcome of dynamic competition among multiple growth options. In both zero-dimensional batch and one-dimensional column configurations, the ANN-based surrogate models achieved substantial reduction of computational time by several orders of magnitude compared to the original LP-based FBA models. Moreover, the ANN models produced robust solutions without any special measures to prevent numerical instability. These developments significantly promote our ability to utilize genome-scale networks in complex, multi-physics, and multi-dimensional ecosystem modeling.

59 BASIC BIOLOGICAL SCIENCES↗