Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “biological system modeling”

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 271 records · Page 15

Multi-strain analysis of Pseudomonas putida reveals the metabolic and genetic diversity of the species

Pseudomonas putida is a gram-negative bacterial species increasingly utilized in biotechnology due to its robust growth, ability to degrade aromatic compounds, solvent tolerance, and genetic tractability. In this study, we report a comprehensive multi-strain analysis of 164 P. putida strains based on the reconstruction of a pan-putida metabolic network and the formulation of strain-specific genome-scale metabolic models (GEMs). We performed whole-genome sequencing and hybrid assembly for 40 strains, contributing a ~8% increase to the available genomic data for P. putida . Furthermore, high-throughput phenotypic profiling using the Biolog phenotype microarray system for 24 strains on 190 unique carbon sources, along with 15 aromatic compounds not present on Biolog plates, yielded 4,920 unique strain-phenotype measurements. These data were leveraged to curate GEMs for 24 representative strains, including a refined model for strain KT2440, which comprised 1,480 genes and 2,191 metabolites, achieving a prediction accuracy of 91.2% in carbon utilization. Systematic comparison of genomes and GEMs revealed both conserved core pathways and significant allelic and functional divergence across strains, highlighting strain-specific variation in aromatic degradation. While pathways for protocatechuate and phenylacetate degradation were widely conserved, metabolic capabilities for compounds such as ferulate, phenol, and cresols varied markedly, suggesting adaptation to distinct ecological niches. Alleleome analysis of enzymes, such as PcaI and PcaJ, revealed distinct, functionally similar clades, indicating possible convergent evolution or horizontal gene transfer. These results provide computable resources and informative models for selecting P. putida strains with desired traits for biomanufacturing and bioremediation and offer insights into the evolution and phylogeny of the P. putida species.

aromatics utilization↗

A machine learning Automated Recommendation Tool for synthetic biology

Abstract Synthetic biology allows us to bioengineer cells to synthesize novel valuable molecules such as renewable biofuels or anticancer drugs. However, traditional synthetic biology approaches involve ad-hoc engineering practices, which lead to long development times. Here, we present the Automated Recommendation Tool (ART), a tool that leverages machine learning and probabilistic modeling techniques to guide synthetic biology in a systematic fashion, without the need for a full mechanistic understanding of the biological system. Using sampling-based optimization, ART provides a set of recommended strains to be built in the next engineering cycle, alongside probabilistic predictions of their production levels. We demonstrate the capabilities of ART on simulated data sets, as well as experimental data from real metabolic engineering projects producing renewable biofuels, hoppy flavored beer without hops, fatty acids, and tryptophan. Finally, we discuss the limitations of this approach, and the practical consequences of the underlying assumptions failing.

59 BASIC BIOLOGICAL SCIENCES↗

Motility of acoustically powered micro-swimmers in a liquid crystalline environment

Suspensions of microswimmers in liquid crystals demonstrate remarkably complex dynamics and serve as a model system for studying active nematics. So far, experimental realization of microswimmers suspended in liquid crystalline media has relied on biological microorganisms that impose strict limitations on the compatible media and makes it difficult to regulate activity. Here, we demonstrate that acoustically powered bubble microswimmers can efficiently self-propel in a lyotropic liquid crystal. The velocity of the swimmers is controlled by the amplitude of the acoustic field. Histograms of swimming directions with respect to the local nematic field reveal a bimodal distribution: the swimmers tend to either fully align with or swim perpendicular to the director field of the liquid crystal, occasionally switching between these two states. The bubble-induced streaming from a swimmer locally melts the liquid crystal and produces topological defects at the tail of the swimmer. Further, we show that the defect proliferation rate increases with the angle between the swimmer's velocity and the local orientation of the director field.

36 MATERIALS SCIENCE↗

A multi-modal scanning system to digitize CBRNE emergency response scenes

A handheld system developed to digitize a contextual understanding of the scene at a chemical, biological, radiological, nuclear and/or explosives (CBRNE) events is described. The system uses LiDAR and cameras to create a colorized 3D model of the environment, which helps domain experts that are supporting responders in the field. To generate the digitized model, a responder scans any suspicious objects and the surroundings by carrying the system through the scene. The scanning system provides a real-time user interface to inform the user about scanning progress and to indicate any areas that may have been missed either by the LiDAR sensors or the cameras. Currently, the collected data are post-processed on a different device, building a colorized triangular mesh of the encountered scene, with the intention of moving this pipeline to the scanner at a later point. The mesh is sufficiently compressed to be sent over a reduced bandwidth connection to a remote analyst. Furthermore, the system tracks fiducial markers attached to diagnostic equipment that is placed around the suspicious object. The resulting tracking information can be transmitted to remote analysts to further facilitate their supporting efforts. The paper will discuss the system's design, software components, the user interface used for scanning a scene, the necessary procedures for calibration of the sensors, and the processing steps of the resulting data. The discussion will close by evaluating the system's performance on 11 scenes.

Salathe, Marco↗

Effects of small marine energy deployments on oceanographic systems

The placement and operation of marine energy deployments in the ocean have the potential to change flow patterns, decrease wave heights, and/or remove energy from the oceanographic system. Changes in oceanographic systems resulting from harvesting marine energy, particularly tidal and wave energy, may be of concern. These changes include alterations in nearfield and farfield physical processes, as well as potential secondary environmental effects such as changes in sediment transport patterns, biological processes, or coastal erosion. Knowledge of changes in oceanographic systems associated with marine energy is primarily available from numerical modeling studies, informed by some laboratory tests and very few field measurements. A literature review was conducted using the Tethys knowledge base and other online sources, building on conclusions from the Ocean Energy Systems-Environmental State of the Science report. Potential changes in oceanographic systems that may be caused by marine energy differ between tidal and wave devices because of different extraction mechanisms and siting locations. Numerical models show that tidal extraction on the order of hundreds of megawatts or with significant channel blockage is required to create changes in oceanographic processes that exceed natural variability. Effects from wave energy extraction in arrays are localized and dependent on array spacing and proximity to the shore. Here the available evidence supports the conclusion that the risk of significant environmental effects from such changes could be retired (i.e., less investigation required for every project) for small deployments—those representative of the state of the industry in 2021. Determining changes in oceanographic systems to be low risk for small deployments can thereby streamline environmental consenting by reducing monitoring needs at this early stage in the industry.

54 ENVIRONMENTAL SCIENCES↗

A multi-modal scanning system to digitize CBRNE emergency response scenes

A handheld system developed to digitize a contextual understanding of the scene at a chemical, biological, radiological, nuclear and/or explosives (CBRNE) events is described. The system uses LiDAR and cameras to create a colorized 3D model of the environment, which helps domain experts that are supporting responders in the field. To generate the digitized model, a responder scans any suspicious objects and the surroundings by carrying the system through the scene. The scanning system provides a real-time user interface to inform the user about scanning progress and to indicate any areas that may have been missed either by the LiDAR sensors or the cameras. Currently, the collected data are post-processed on a different device, building a colorized triangular mesh of the encountered scene, with the intention of moving this pipeline to the scanner at a later point. The mesh is sufficiently compressed to be sent over a reduced bandwidth connection to a remote analyst. Furthermore, the system tracks fiducial markers attached to diagnostic equipment that is placed around the suspicious object. The resulting tracking information can be transmitted to remote analysts to further facilitate their supporting efforts. The paper will discuss the system's design, software components, the user interface used for scanning a scene, the necessary procedures for calibration of the sensors, and the processing steps of the resulting data. The discussion will close by evaluating the system's performance on 11 scenes.

Prins, Nicholas↗

Substantial and overlooked greenhouse gas emissions from deep Arctic lake sediment - supporting data and code

This data package contains data, descriptions, and code-based analyses that were used to support conclusions drawn in “Substantial and overlooked greenhouse gas emissions from deep Arctic lake sediment”, by Freitas et al. (2025) (https://doi.org/10.1038/s41561-024-01614-y). The study evaluated greenhouse gas production along a deep sediment core (20 m) taken in 2018 from below Goldstream Lake, a field site approximately 15 km north of Fairbanks, Alaska.The file “ESSDive_NFreitas_2024_flmd.csv” includes an overview of all other csv files in this data package, namely: sediment descriptions (depth and type of sediment), sediment characterizations (bulk density, gravimetric water content, total carbon, etc.), calculated respiration and temperature sensitivity values associated with year-long incubations of the sediment core, and the R code used to process the dataset. Additional details regarding the content of these files and how the calculations were performed are described in the Methods section of this archive. The “data_dictionary_ESSDive_NFreitas_2024_dd.csv” is a data dictionary for all files included in the data package. Each row in the data dictionary represents a column name in a given file. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).Data package updates- August 2024: Title was updated, R code now includes a process for calculating potential production across the whole sediment column, the "sediment_column_production_potentials_ESSDive_NFreitas_2024.csv" is the associated data output file, and the file level metadata and data dictionaries were updated to reflect the contents of the additional csv file.- October 2024: Title was updated, and a statement was added to Step 2 of the Methods that describes that samples were collected (and exported) in a responsible manner and in accordance with relevant permits and local laws.- January 2025: Updated the associated manuscript details in the Abstract (manuscript title, publication date, DOI link) and in the Related References section (full citation).

54 ENVIRONMENTAL SCIENCES↗

Mechanistic basis of post-treatment control of SIV after anti-α4β7 antibody therapy

Treating macaques with an anti-α4β7 antibody under the umbrella of combination antiretroviral therapy (cART) during early SIV infection can lead to viral remission, with viral loads maintained at < 50 SIV RNA copies/ml after removal of all treatment in a subset of animals. Depletion of CD8+ lymphocytes in controllers resulted in transient recrudescence of viremia, suggesting that the combination of cART and anti-α4β7 antibody treatment led to a state where ongoing immune responses kept the virus undetectable in the absence of treatment. A previous mathematical model of HIV infection and cART incorporates immune effector cell responses and exhibits the property of two different viral load set-points. While the lower set-point could correspond to the attainment of long-term viral remission, attaining the higher set-point may be the result of viral rebound. Here we expand that model to include possible mechanisms of action of an anti-α4β7 antibody operating in these treated animals. We show that the model can fit the longitudinal viral load data from both IgG control and anti-α4β7 antibody treated macaques, suggesting explanations for the viral control associated with cART and an anti-α4β7 antibody treatment. This effective perturbation to the virus-host interaction can also explain observations in other nonhuman primate experiments in which cART and immunotherapy have led to post-treatment control or resetting of the viral load set-point. Interestingly, because the viral kinetics in the various treated animals differed—some animals exhibited large fluctuations in viral load after cART cessation—the model suggests that anti-α4β7 treatment could act by different primary mechanisms in different animals and still lead to post-treatment viral control. This outcome is nonetheless in accordance with a model with two stable viral load set-points, in which therapy can perturb the system from one set-point to a lower one through different biological mechanisms.

60 APPLIED LIFE SCIENCES↗

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the “gold standard” coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller−Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Chemistry↗

Plant Characteristics, Porewater, Gas Flux, and Soil Biogeochemistry at Council Road Site Mile Marker 71, Seward Peninsula, Alaska, 2023

Data collected at Council, AK (64°51’35.0”N 163°41’59.1”W) during a summer campaign in 2023. Water data consists of soil porewater collected by centrifuging soil cores and also by field collection with porewater samplers (rhizons). Gas data consists of CO2, CH4 and N2O surface soil fluxes measured with a portable FTIR analyzer. Plant and root data consists of biomass, root length, diameter and mass. Soil data consists of total C and N. Air, water and soil samples span two main locations: a thermokarst wetland and an upland tussock. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).This dataset was generated to broadly address the following research question: how will climate change (i.e., thawing permafrost, landscape change) alter the ecosystem flux (sink versus source) of important greenhouse gases such as CO2, CH4 and N2O?Description of the contents of this data package: This dataset contains 5 different individual .csv files containing plant, soil, water and gas data. No software is needed to utilize them. PFTCover: Plant functional type ground cover in 1x1 meter plots. SoilCores: Solidphase and porewater phase soil biogeochemical variablesPlantData: Above and belowground plant traits. GasFlux: Surface plant-soil gas measurementsFieldPorewater: Field collected porewater biogeochemical variables

54 ENVIRONMENTAL SCIENCES↗

Geometric frustration in the myosin superlattice of vertebrate muscle

Geometric frustration results from an incompatibility between minimum energy arrangements and the geometry of a system, and gives rise to interesting and novel phenomena. Here, we report geometric frustration in a native biological macromolecular system---vertebrate muscle. We analyse the disorder in the myosin filament rotations in the myofibrils of vertebrate striated (skeletal and cardiac) muscle, as seen in thin-section electron micrographs, and show that the distribution of rotations corresponds to an archetypical geometrically frustrated system---the triangular Ising antiferromagnet. Spatial correlations are evident out to at least six lattice spacings. The results demonstrate that geometric frustration can drive the development of structure in complex biological systems, and may have implications for the nature of the actin--myosin interactions involved in muscle contraction. Identification of the distribution of myosin filament rotations with an Ising model allows the extensive results on the latter to be applied to this system. It shows how local interactions (between adjacent myosin filaments) can determine long-range order and, conversely, how observations of long-range order (such as patterns seen in electron micrographs) can be used to estimate the energetics of these local interactions. Furthermore, since diffraction by a disordered system is a function of the second-order statistics, the derived correlations allow more accurate diffraction calculations, which can aid in interpretation of X-ray diffraction data from muscle specimens for structural analysis.

59 BASIC BIOLOGICAL SCIENCES↗

Ex vivo to in vivo model of malignant peripheral nerve sheath tumors for precision oncology

Malignant peripheral nerve sheath tumors (MPNST) are aggressive soft tissue sarcomas that often develop in patients with neurofibromatosis type 1 (NF1). To address the critical need for novel therapeutics in MPNST, we aimed to establish an ex vivo 3D platform that accurately captured the genomic diversity of MPNST and could be utilized in a medium-throughput manner for drug screening studies to be validated in vivo using patient-derived xenografts (PDX). Genomic analysis was performed on all PDX-tumor pairs. Selected PDX were harvested for assembly into 3D microtissues. Based on prior work in our labs, we evaluated drugs (trabectedin, olaparib, and mirdametinib) ex vivo and in vivo. For 3D microtissue studies, cell viability was the endpoint as assessed by Zeiss Axio Observer. For PDX drug studies, tumor volume was measured twice weekly. Bulk RNA sequencing was performed to identify pathways enriched in cells. We developed 13 NF1-associated MPNST-PDX and identified mutations or structural abnormalities in NF1 (100%), SUZ12 (85%), EED (15%), TP53 (15%), CDKN2A (85%), and chromosome 8 gain (77%). We successfully assembled PDX into 3D microtissues, categorized as robust (>90% viability at 48 h), good (>50%), or unusable (<50%). We evaluated drug response to “robust” or “good” microtissues, namely MN-2, JH-2-002, JH-2-079-c, and WU-225. Drug response ex vivo predicted drug response in vivo, and enhanced drug effects were observed in select models. These data support the successful establishment of a novel 3D platform for drug discovery and MPNST biology exploration in a system representative of the human condition.

60 APPLIED LIFE SCIENCES↗

Privacy Preservation from High-Performance Computing to Autonomous Science [Industrial and Governmental Activities]

High-Performance Computing (HPC) and Leadership-Class Supercomputing are driving forces behind scientific advancements, enabling researchers to tackle complex challenges in physics, chemistry, biology, and engineering. These systems power vast simulations and data analyses, fueling discoveries in fields ranging from materials science to climate modeling. However, their use often involves processing sensitive data—such as proprietary industry simulations, biomedical records, and national security computations—posing significant privacy concerns. In conclusion, this issue is amplified in collaborative environments like Department of Energy (DOE) user facilities, where HPC resources are shared across institutions to foster innovation.

Kotevska, Olivera [Oak Ridge National Laboratory (↗

Resilience of the slow component in timescale-separated synchronized oscillators

Physiological networks are usually made of a large number of biological oscillators evolving on a multitude of different timescales. Phase oscillators are particularly useful in the modelling of the synchronization dynamics of such systems. If the coupling is strong enough compared to the heterogeneity of the internal parameters, synchronized states might emerge where phase oscillators start to behave coherently. Here, we focus on the case where synchronized oscillators are divided into a fast and a slow component so that the two subsets evolve on separated timescales. We assess the resilience of the slow component by, first, reducing the dynamics of the fast one using Mori-Zwanzig formalism. Second, we evaluate the variance of the phase deviations when the oscillators in the two components are subject to noise with possibly distinct correlation times. From the general expression for the variance, we consider specific network structures and show how the noise transmission between the fast and slow components is affected. Interestingly, we find that oscillators that are among the most robust when there is only a single timescale, might become the most vulnerable when the system undergoes a timescale separation. We also find that layered networks seem to be insensitive to such timescale separations.

59 BASIC BIOLOGICAL SCIENCES↗

Modular cell-free expression plasmids to accelerate biological design in cells

Industrial biotechnology aims to produce high-value products from renewable resources. This can be challenging because model microorganisms—organisms that are easy to use like Escherichia coli—often lack the machinery required to utilize desired feedstocks like lignocellulosic biomass or syngas. Non-model organisms, such as Clostridium, are industrially proven and have desirable metabolic features but have several hurdles to mainstream use. Namely, these species grow more slowly than conventional laboratory microbes, and genetic tools for engineering them are far less prevalent. To address these hurdles for accelerating cellular design, cell-free synthetic biology has matured as an approach for characterizing non-model organisms and rapidly testing metabolic pathways in vitro. Unfortunately, cell-free systems can require specialized DNA architectures with minimal regulation that are not compatible with cellular expression. In this work, we develop a modular vector system that allows for T7 expression of desired enzymes for cell-free expression and direct Golden Gate assembly into Clostridium expression vectors. Utilizing the Joint Genome Institute’s DNA Synthesis Community Science Program, we designed and synthesized these plasmids and genes required for our projects allowing us to shuttle DNA easily between our in vitro and in vivo experiments. We next validated that these vectors were sufficient for cell-free expression of functional enzymes, performing on par with the previous state-of-the-art. Lastly, we demonstrated automated six-part DNA assemblies for Clostridium autoethanogenum expression with efficiencies ranging from 68% to 90%. We anticipate this system of plasmids will enable a framework for facile testing of biosynthetic pathways in vitro and in vivo by shortening development cycles.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling COVID-19 Impacts on Prison Population

The goal of this project was to develop a unique, high-level, concept model that will (1) help examine the consequences COVID-19 has had on prison population and flow dynamics in Illinois and (2) explore (to the extent possible) how that will unfold in the short- and mid-term future. Argonne has unique expertise in dynamic systems modeling and highly specialized simulation capabilities, which will be helpful in developing such a model. The ultimate objective of this CRADA is to develop with the Illinois Public Sentencing Policy Advisory Council (SPAC) the capability to identify and understand factors and interactions associated with prison population and flow dynamics in Illinois. This understanding will then potentially lead to improvements in the allocation of various resources and services. This approach would also be useful to other states.

59 BASIC BIOLOGICAL SCIENCES↗

Designing, synthesizing, and modeling active fluids

We review recent advances in the design, synthesis, and modeling of active fluids. Active fluids have been at the center of many technological innovations and theoretical advances over the past two decades. Research on this new class of fluids has been inspired by the fascinating and remarkably efficient strategies that biological systems employ, leading to the development of biomimetic nano- and micro-machines and swimmers. The review encompasses active fluids on both the nano- and micro-scale. We start with examples of biological active systems before we discuss how experimentalists leverage novel propulsion mechanisms to power nano- and micro-machines. Furthermore, we then examine how the study of these far-from-equilibrium systems has prompted the development of new simulation methods and theoretical models in nonequilibrium physics to account for their mechanical, thermodynamic, and emergent properties. Recent advances in the field have paved the way for the design, synthesis, and modeling of autonomous systems at the nano- and micro-scale and opened the door to the development of soft matter robotics.

36 MATERIALS SCIENCE↗