Results for “physiological computing”

DOE OSTI2025

Structural and Mechanistic Advances in the Chemistry of Methyl-Coenzyme M Reductase (MCR)

Methane represents 34% of U.S. energy consumption and is a major greenhouse gas related to the global carbon cycle and energy production. However, current industrial practices significantly increase atmospheric methane levels, necessitating a deeper understanding of its biosynthesis and oxidation. Methyl-coenzyme M reductase (MCR) is central to biological methane metabolism, catalyzing the final step of methanogenesis and the first step in anaerobic methane oxidation. It is also a key target for strategies to capture and transform methane into value-added chemicals.The active site of MCR is a buried Ni-based cofactor only accessible by the substrates via a 50 Å long tunnel. Although the Ni(I) state is required to initiate catalysis, capturing this state remains a challenge for the current structural techniques. Recent advances in structural biology using X-ray Free-Electron Laser serial crystallography have provided insights into MCR's inactive Ni(II) state at room temperature and show promise for capturing its active Ni(I) form.Our team has established several critical aspects of the MCR mechanism using a combination of experimental and computational studies. MCR uses CH 3 -SCoM and CoBSH as substrates, producing methane and a disulfide product CoMSSCoB. Kinetic analysis showed that productive substrate binding requires CH 3 -SCoM to bind first, inducing conformational changes that optimize the active site for subsequent CoBSH binding. Following substrate binding, four proposed methane production/oxidation mechanisms were examined, establishing whether the reaction proceeds through an organometallic methyl-nickel(III), methyl anion ion, or methyl radical intermediate. Experimental measurements using CoBSH analogs successfully slowed the reaction, allowing for mechanistic insight that demonstrated the methyl radical pathway, where the initial interactions involve homolytic cleavage of the methyl-sulfur bond, generating a methyl radical that quickly abstracts the thiol hydrogen atom of CoBSH to form methane. Computational studies further confirmed that, compared to other mechanisms, the methyl radical mechanism is thermodynamically more favorable and accessible under physiological conditions.Spectroscopic and computational studies challenged the conventional understanding of substrate binding in MCR by proposing an alternative positioning of CH 3 -SCoM and CoMSSCoB in the active site pocket. The research suggested that CH 3 -SCoM (substrate) and CoMSSCoB (product) bind via their sulfonate groups to the Ni(I) center of cofactor F 430 . This binding allows for the reaction without substrate reorganization in the pocket but would require a long-range electron transfer. Altogether, the work summarized in this review reflects our current understanding of the enzyme's catalytic mechanism and structural dynamics. This is essential for developing efficient methane conversion technologies that could mitigate its environmental impact while harnessing energy-storage properties.

03 NATURAL GAS↗

DOE OSTI2025

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗

DOE OSTI2023

A computational modeling framework for pre-clinical evaluation of cardiac mapping systems

There are a variety of difficulties in evaluating clinical cardiac mapping systems, most notably the inability to record the transmembrane potential throughout the entire heart during patient procedures which prevents the comparison to a relevant “gold standard”. Cardiac mapping systems are comprised of hardware and software elements including sophisticated mathematical algorithms, both of which continue to undergo rapid innovation. The purpose of this study is to develop a computational modeling framework to evaluate the performance of cardiac mapping systems. The framework enables rigorous evaluation of a mapping system’s ability to localize and characterize (i.e., focal or reentrant) arrhythmogenic sources in the heart. The main component of our tool is a library of computer simulations of various dynamic patterns throughout the entire heart in which the type and location of the arrhythmogenic sources are known. Our framework allows for performance evaluation for various electrode configurations, heart geometries, arrhythmias, and electrogram noise levels and involves blind comparison of mapping systems against a “silver standard” comprised of computer simulations in which the precise transmembrane potential patterns throughout the heart are known. A feasibility study was performed using simulations of patterns in the human left atria and three hypothetical virtual catheter electrode arrays. Activation times (AcT) and patterns (AcP) were computed for three virtual electrode arrays: two basket arrays with good and poor contact and one high-resolution grid with uniform spacing. The average root mean squared difference of AcTs of electrograms and those of the nearest endocardial action potential was less than 1 ms and therefore appears to be a poor performance metric. In an effort to standardize performance evaluation of mapping systems a novel performance metric is introduced based on the number of AcPs identified correctly and those considered spurious as well as misclassifications of arrhythmia type; spatial and temporal localization accuracy of correctly identified patterns was also quantified. This approach provides a rigorous quantitative analysis of cardiac mapping system performance. Proof of concept of this computational evaluation framework suggests that it could help safeguard that mapping systems perform as expected as well as provide estimates of system accuracy.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Coarse-Grained Modeling Elucidates Differential Metabolism of Saccharomyces cerevisiae under Varied Nutrient Limitations

Microorganisms such as Saccharomyces cerevisiae have a native ability to adapt their metabolism to varying nutrient conditions. Understanding their responses to nutrient limitations is critical for decoding cellular physiology and designing strategies for metabolic engineering. While the influence of carbon availability on yeast metabolism has been extensively studied, the role of nitrogen availability remains relatively underexplored. In this study, we utilized a coarse-grained kinetic model to systematically analyze and compare the effects of carbon and nitrogen limitations on yeast metabolism. Our model successfully revealed the differential metabolic characteristics of S. cerevisiae under carbon- and nitrogen-limited chemostat conditions. It also highlighted the significance of protein activity regulation at varying carbon-to-nitrogen ratios, and elucidated distinct strategies employed to maintain ATP homeostasis. Furthermore, this study provides a computational tool for investigating yeast physiology under nutrient limitations and offers quantitative and mechanistic insights into yeast metabolism.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Better together: Elements of successful scientific software development in a distributed collaborative community

Many scientific disciplines rely on computational methods for data analysis, model generation, and prediction. Implementing these methods is often accomplished by researchers with domain expertise but without formal training in software engineering or computer science. This arrangement has led to underappreciation of sustainability and maintainability of scientific software tools developed in academic environments. Some software tools have avoided this fate, including the scientific library Rosetta. We use this software and its community as a case study to show how modern software development can be accomplished successfully, irrespective of subject area. Rosetta is one of the largest software suites for macromolecular modeling, with 3.1 million lines of code and many state-of-the-art applications. Since the mid 1990s, the software has been developed collaboratively by the RosettaCommons, a community of academics from over 60 institutions worldwide with diverse backgrounds including chemistry, biology, physiology, physics, engineering, mathematics, and computer science. Developing this software suite has provided us with more than two decades of experience in how to effectively develop advanced scientific software in a global community with hundreds of contributors. Here we illustrate the functioning of this development community by addressing technical aspects (like version control, testing, and maintenance), community-building strategies, diversity efforts, software dissemination, and user support. We demonstrate how modern computational research can thrive in a distributed collaborative community. The practices described here are independent of subject area and can be readily adopted by other software development communities

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Biokinetic and Dosimetric Models [Slides]

Presentation describes internal dosimetry as an intersection of disciplines involving Physiology, Anatomy, Physics, Mathematics and Computer Science, which some people call "The art of Internal Dosimetry". The basic components for internal dose calculations are: Biokinetic (metabolic) models describing the intake, distribution, retention and excretion of radionuclides in the body, Dosimetric models describing the interaction of radiation within the several body organs and tissues, and System of dose limitation.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2020

Accelerating massively parallel hemodynamic models of coarctation of the aorta using neural networks

Comorbidities such as anemia or hypertension and physiological factors related to exertion can influence a patient’s hemodynamics and increase the severity of many cardiovascular diseases. Observing and quantifying associations between these factors and hemodynamics can be difficult due to the multitude of co-existing conditions and blood flow parameters in real patient data. Machine learning-driven, physics-based simulations provide a means to understand how potentially correlated conditions may affect a particular patient. Here, we use a combination of machine learning and massively parallel computing to predict the effects of physiological factors on hemodynamics in patients with coarctation of the aorta. We first validated blood flow simulations against in vitro measurements in 3D-printed phantoms representing the patient’s vasculature. We then investigated the effects of varying the degree of stenosis, blood flow rate, and viscosity on two diagnostic metrics – pressure gradient across the stenosis (ΔP) and wall shear stress (WSS) - by performing the largest simulation study to date of coarctation of the aorta (over 70 million compute hours). Using machine learning models trained on data from the simulations and validated on two independent datasets, we developed a framework to identify the minimal training set required to build a predictive model on a per-patient basis. We then used this model to accurately predict ΔP (mean absolute error within 1.18 mmHg) and WSS (mean absolute error within 0.99 Pa) for patients with this disease.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Formation of a constructed microbial community in a nutrient-rich environment indicates bacterial interspecific competition

ABSTRACT Understanding the organizational principles of microbial communities is essential for interpreting ecosystem stability. Previous studies have investigated the formation of bacterial communities under nutrient-poor conditions or obligate relationships to observe cooperative interactions among different species. How microorganisms form stabilized communities in nutrient-rich environments, without obligate metabolic interdependency for growth, is still not fully disclosed. In this study, three bacterial strains isolated from the Populus deltoides rhizosphere were co-cultured in complex medium, and their growth behavior was tracked. These strains co-exist in mixed culture over serial transfer for multiple growth-dilution cycles. Competition is proposed as an emergent interaction relationship among the three bacteria based on their significantly decreased growth levels. The effects of different initial inoculum ratios, up to three orders of magnitude, on community structure were investigated, and the final compositions of the mixed communities with various starting composition indicate that community structure is not dependent on the initial inoculum ratio. Furthermore, the competitive relationships within the community were not altered by different initial inoculum ratios. The community structure was simulated by generalized Lotka-Volterra and dynamic flux balance analysis to provide mechanistic predictions into emergence of community structure under a nutrient-rich environment. Metaproteomic analyses provide support for the metabolite exchanges predicted by computational modeling and for highly altered physiologies when microbes are grown in co-culture. These findings broaden our understanding of bacterial community dynamics and metabolic diversity in higher-order interactions and could be significant in the management of rhizospheric bacterial communities. IMPORTANCE Bacteria naturally co-exist in multispecies consortia, and the ability to engineer such systems can be useful in biotechnology. Despite this, few studies have been performed to understand how bacteria form a stable community and interact with each other under nutrient-rich conditions. In this study, we investigated the effects of initial inoculum ratios on bacterial community structure using a complex medium and found that the initial inoculum ratio has no significant impact on resultant community structure or on interaction patterns between community members. The microbial population profiles were simulated using computational tools in order to understand intermicrobial relationships and to identify potential metabolic exchanges that occur during stabilization of the bacterial community. Studying microbial community assembly processes is essential for understanding fundamental ecological principles in microbial ecosystems and can be critical in predicting microbial community structure and function.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

National User Resource for Biological Accelerator Mass Spectrometry

The National User Resource for Biological Accelerator Mass Spectrometry (User Resource) will provide isotopic analysis (primarily radiocarbon or 14C) by accelerator mass spectrometry (AMS) for NIH- funded researchers across the United States and will be the only User Resource of its type in the United States. The User Resource will provide measurement capability and expertise to a research community that requires highly sensitive, quantitative isotope analyses. Since commissioning a new accelerator mass spectrometer in June 2014, we have measured over 4000 samples a year for collaborators and service users. The User Resource will enable us to continue to meet these research needs, as well as provide for new users whose research programs would benefit from AMS as a measurement tool. The User Resource’s forte will be ultra-high sensitivity quantitation of radiocarbon and selected other radioisotopes for research studies where isotopes are required. Radioisotope labeling studies have been and will continue to be an important tool for addressing many complex biomedical science problems. AMS is a specialized and unique type of mass spectrometry that provides absolute quantitation of radiocarbon and other relevant radioisotopes with extreme sensitivity, having limits of detection in real samples on the order of a few attomol/mg of sample at measurement precisions of ~3%. It is the only instrumental method capable of quantifying radioisotope-labeled agents routinely in real-world samples with such precision and sensitivity. The sensitivity of AMS allows for the quantification of radiolabeled metabolites in extremely complex matrices of cells and organisms at very low concentrations and in small samples. AMS allows studies to be conducted without perturbing metabolism leading to more relevant quantification of metabolic rates and pathways. In addition, it enables quantification of pharmacokinetic and metabolic properties of toxicants at environmentally relevant concentrations in model systems as well as the ability to quantify pharmacokinetics and other molecular endpoints directly in humans. Such quantitative assessments can 1) improve risk assessment for toxicants, 2) address safety and efficacy considerations for therapeutic entities, 3) deepen understanding of xenobiotic and intermediary metabolism, 4) help understand the interactions between critical molecular pathways, and 5) improve efforts to model and predict various metabolic and biological states. These capabilities have been applied in a number of areas including research in carcinogenesis, toxicology, nutrition, pharmacology/drug development and basic biological science. As a NIGMS National Resource the National User Resource for Biological Accelerator Mass Spectrometry will help NIH funded scientists achieve a deeper understanding of the etiology of human health concerns by (1) enabling the quantification of pharmacokinetics and other molecular endpoints directly in humans; (2) offering the ability to conduct quantitative studies using biologics such as proteins or lipids, and thereby reducing the amount of radioisotope usage in biomedical labs; and (3) enabling more relevant studies of metabolic pathways in health and disease through the use of much lower, more biologically-relevant, concentrations of metabolic substrates in cells and intact organisms. Such studies support NIGMS’s basic biomedical research areas that contribute to the understanding of fundamental cellular and physiological principles and enable research supported by the Biophysics, Biomedical Technology, and Computational Biosciences (BBCB); Genetics and Molecular, Cellular, and Developmental Biology (GMCDB); Pharmacology, Physiology, Biological Chemistry (PPBC) and Training, Workforce Development, and Diversity (TWD) Divisions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

National User Resource for Biological Accelerator Mass Spectrometry (Final Report)

The National User Resource for Biological Accelerator Mass Spectrometry (User Resource) will provide isotopic analysis (primarily radiocarbon or 14C) by accelerator mass spectrometry (AMS) for NIH- funded researchers across the United States and will be the only User Resource of its type in the United States. The User Resource will provide measurement capability and expertise to a research community that requires highly sensitive, quantitative isotope analyses. Since commissioning a new accelerator mass spectrometer in June 2014, we have measured over 4000 samples a year for collaborators and service users. The User Resource will enable us to continue to meet these research needs, as well as provide for new users whose research programs would benefit from AMS as a measurement tool. The User Resource’s forte will be ultra-high sensitivity quantitation of radiocarbon and selected other radioisotopes for research studies where isotopes are required. Radioisotope labeling studies have been and will continue to be an important tool for addressing many complex biomedical science problems. AMS is a specialized and unique type of mass spectrometry that provides absolute quantitation of radiocarbon and other relevant radioisotopes with extreme sensitivity, having limits of detection in real samples on the order of a few attomol/mg of sample at measurement precisions of ~3%. It is the only instrumental method capable of quantifying radioisotope-labeled agents routinely in real-world samples with such precision and sensitivity. The sensitivity of AMS allows for the quantification of radiolabeled metabolites in extremely complex matrices of cells and organisms at very low concentrations and in small samples. AMS allows studies to be conducted without perturbing metabolism leading to more relevant quantification of metabolic rates and pathways. In addition, it enables quantification of pharmacokinetic and metabolic properties of toxicants at environmentally relevant concentrations in model systems as well as the ability to quantify pharmacokinetics and other molecular endpoints directly in humans. Such quantitative assessments can 1) improve risk assessment for toxicants, 2) address safety and efficacy considerations for therapeutic entities, 3) deepen understanding of xenobiotic and intermediary metabolism, 4) help understand the interactions between critical molecular pathways, and 5) improve efforts to model and predict various metabolic and biological states. These capabilities have been applied in a number of areas including research in carcinogenesis, toxicology, nutrition, pharmacology/drug development and basic biological science. As a NIGMS National Resource the National User Resource for Biological Accelerator Mass Spectrometry will help NIH funded scientists achieve a deeper understanding of the etiology of human health concerns by (1) enabling the quantification of pharmacokinetics and other molecular endpoints directly in humans; (2) offering the ability to conduct quantitative studies using biologics such as proteins or lipids, and thereby reducing the amount of radioisotope usage in biomedical labs; and (3) enabling more relevant studies of metabolic pathways in health and disease through the use of much lower, more biologically-relevant, concentrations of metabolic substrates in cells and intact organisms. Such studies support NIGMS’s basic biomedical research areas that contribute to the understanding of fundamental cellular and physiological principles and enable research supported by the Biophysics, Biomedical Technology, and Computational Biosciences (BBCB); Genetics and Molecular, Cellular, and Developmental Biology (GMCDB); Pharmacology, Physiology, Biological Chemistry (PPBC) and Training, Workforce Development, and Diversity (TWD) Divisions. Over the next five years, our goals are to: 1. Improve the efficiency of operation for AMS measurements through installation of new interfaces to our AMS systems, technical modifications to improve gas accepting ion source efficiency and upgrading our data analysis software for improved ease of use and data reporting. 2. Increase the accessibility and visibility of ultra-sensitive 14C measurements for the biomedical research community by training of new investigators and expanding our national user base. 3. Provide high throughput, ultra-sensitive 14C analysis for the NIGMS and NIH user community.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Multigene engineering in plants: Technologies, applications, and future prospects

The emerging bioeconomy presents a promising solution to both economic and environmental challenges. Within the bioeconomy, plants serve as a renewable, sustainable, and cost-effective source of foods, fuels, chemicals, and materials. However, traditional breeding and single-gene engineering approaches fall short in addressing complex traits (e.g., drought tolerance, disease resistance, yield, nutrient use efficiency) which are controlled by multiple genes. The complexity of plant biology often necessitates the use of multigene engineering (MGE), which involves simultaneous ectopic expression, up/down-regulation, or editing of multiple genes, to enhance plant traits relevant to the bioeconomy. These genes may be associated with distinct traits or function as components of specific metabolic and regulatory pathways. This review summarizes current technologies for MGE within the synthetic biology-driven Design-Build-Test-Learn (DBTL) framework, detailing its four key stages: Design – gene construct development; Build – DNA assembly and plant transformation; Test – the molecular, biochemical, and physiological characterization of engineered plants; and Learn – computational modeling to refine, multiplex and iterate the process. Despite good progress in the applications of MGE in biofortification, metabolic engineering, and stress resilience, challenges remain in construct stability, coordinated gene expression, and regulatory predictability. We identified optimization paths and future directions to accelerate MGE deployment in sustainable agriculture, with possible societal benefits including reduced production costs, increased yield, and improved food and nutritional security.

AI-aided plant engineering↗

DOE OSTI2020

Phosphoproteome Response to Dithiothreitol Reveals Unique Versus Shared Features of Saccharomyces cerevisiae Stress Responses

To cope with sudden changes in the external environment, the budding yeast Saccharomyces cerevisiae orchestrates a multifaceted response that spans many levels of physiology. Several studies have interrogated the transcriptome response to endoplasmic reticulum (ER) stress and the role of regulators such as the Ire1 kinase and Hac1 transcription factors. However, less is known about responses to ER stress at other levels of physiology. Here, we used quantitative phosphoproteomics and computational network inference to uncover the yeast phosphoproteome response to the reducing agent dithiothreitol (DTT) and the upstream signaling network that controls it. We profiled wild-type cells and mutants lacking IRE1 or MAPK kinases MKK1 and MKK2, before and at various times after DTT treatment. In addition to revealing downstream targets of these kinases, our inference approach predicted new regulators in the DTT response, including cell-cycle regulator Cdc28 and osmotic-response kinase Rck2, which we validated computationally. Our results also revealed similarities and surprising differences in responses to different stress conditions, especially in the response of protein kinase A targets. Furthermore, these results have implications for the breadth of signaling programs that can give rise to common stress response signatures.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Cabin Thermal Management Analysis for SuperTruck II Next-Generation Hybrid Electric Truck Design

In this article, we present a multistage, coupled thermal management simulation approach, informed by physical testing where available, to aid design decisions for PACCAR's SuperTruck II hybrid truck cabin concept. Focus areas include cabin insulation, battery sizing, and sleeper curtain position, as well as heating, ventilating, and air-conditioning (HVAC) component and accessory configurations, to maintain or improve thermal comfort while saving energy. The authors analyzed weather data and determined the national vehicle miles traveled weighted temperature and solar conditions for long-haul trucks. Example weather day profiles were selected to approximate the 5th and 95th percentile weighted conditions. A daylong drive cycle was developed to impose appropriate external wind conditions during rest and driving periods. Using the National Renewable Energy Laboratory's vehicle HVAC modeling and simulation tool VTCab, HVAC load design trade-off studies for the new truck geometry concept were completed. Parameters analyzed included effects of paint color, insulation, glass transmissivity, and curtain location. Simulation results helped with early design material selections for efficient cabin climate control. A detailed three-dimensional computer-aided engineering (CAE), computational fluid dynamics (CFD), radiation, and human physiology co-simulation, referred to in this article as CAE Thermal-CFD, was used to evaluate thermal comfort and energy impacts of diffuser configurations and air supply settings in driving and hoteling modes. Analysis revealed that it is more difficult to heat the cabin in hoteling mode during the winter than to cool the space in the summer. This seasonal load profile drives the requirement of additional energy storage for heating comfort. To determine the battery capacity requirement, multiday HVAC operation drive cycle simulations were then completed, showing that a 15-kWh battery would be required for HVAC operation during hoteling. Results helped reduce cabin thermal loads, determine component sizing requirements, and improve occupant comfort to save fuel and contribute to the economic viability of the hybrid system.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Structural organization of the spongy mesophyll

Summary Many plant leaves have two layers of photosynthetic tissue: the palisade and spongy mesophyll. Whereas palisade mesophyll consists of tightly packed columnar cells, the structure of spongy mesophyll is not well characterized and often treated as a random assemblage of irregularly shaped cells. Using micro‐computed tomography imaging, topological analysis, and a comparative physiological framework, we examined the structure of the spongy mesophyll in 40 species from 30 genera with laminar leaves and reticulate venation. A spectrum of spongy mesophyll diversity encompassed two dominant phenotypes: first, an ordered, honeycomblike tissue structure that emerged from the spatial coordination of multilobed cells, conforming to the physical principles of Euler’s law; and second, a less‐ordered, isotropic network of cells. Phenotypic variation was associated with transitions in cell size, cell packing density, mesophyll surface‐area‐to‐volume ratio, vein density, and maximum photosynthetic rate. These results show that simple principles may govern the organization and scaling of the spongy mesophyll in many plants and demonstrate the presence of structural patterns associated with leaf function. This improved understanding of mesophyll anatomy provides new opportunities for spatially explicit analyses of leaf development, physiology, and biomechanics.

3D↗

DOE OSTI2021

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. Overall, this protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Examining metastatic behavior within 3D bioprinted vasculature for the validation of a 3D computational flow model

Understanding the dynamics of circulating tumor cell (CTC) behavior within the vasculature has remained an elusive goal in cancer biology. To elucidate the contribution of hydrodynamics in determining sites of CTC vascular colonization, the physical forces affecting these cells must be evaluated in a highly controlled manner. To this end, we have bioprinted endothelialized vascular beds and perfused these constructs with metastatic mammary gland cells under physiological flow rates. By pairing these in vitro devices with an advanced computational flow model, we found that the bioprinted analog was readily capable of evaluating the accuracy and integrated complexity of a computational flow model, while also highlighting the discrete contribution of hydrodynamics in vascular colonization. This intersection of these two technologies, bioprinting and computational simulation, is a key demonstration in the establishment of an experimentation pipeline for the understanding of complex biophysical events.

42 ENGINEERING↗

DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗

DOE OSTI2025

EvoNet: A phylogenomic and systems biology approach to identify genes underlying plant survival in marginal, low‐N soils

The DOE‐BER “EvoNet” project investigates the genetic and molecular basis of plant resilience in extreme environments. We do this by identifying key genes that enable “extreme survivor” species to thrive in the nitrogen-poor soils of Chile’s hyper-arid Atacama Desert. Our collections focus on 32 Atacama extremophile species, including seven grass species with potential biofuel applications. To identify genes-of-importance to survival we compared genomic and transcriptomic profiles of extremophile species that thrive in the Atacama to those of closely related “sister” species from nitrogen-rich arid and mesic regions of California. Deep RNA sequencing and de novo transcriptome assembly across these triplet species sets supported a phylogenomic framework for identifying positively selected genes associated with adaptive divergence. Our integrative analysis combined ecological and environmental data, metagenomics, evolutionary and systems biology, and metabolomics. This enabled us to create an unprecedented framework for systematically understanding how non-model plants have adapted to survive in extreme conditions. Our resulting database of positively selected ortholog groups in the extremophile plants offers promising targets for engineering crop and biofuel species with enhanced resilience to drought and extreme weather. Additionally, our newest dataset explores and exploits a complementary metabolomic approach. This new aspect provides innovative strategies to manipulate plant cell metabolism, further supporting efforts to improve agricultural productivity in the face of extreme climates. Importantly, our combined evolutionary- and metabolomic-based strategies focused on convergent patterns of adaptation, providing a genetic and metabolomic toolkit for improving crop and biofuel resilience across diverse plant species. Finally, our novel exploration of ecological and evolutionary dynamics delivered to the community a phylogenomic computational pipeline called “PhyloGeneious.” Our continued adaptations of this pipeline are publicly available to expedite evolutionary genomic research for future scientific discoveries. In total, our DOE-BER has provided genomic, metabolomic, and computational strategies to understand how extremophile plants provide evolutionary and physiological targets for improving agricultural and biofuel production.

59 BASIC BIOLOGICAL SCIENCES↗