Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematical biology”

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 469 records · Page 26

Biosynthesis of Haloterpenoids in Red Algae via Microbial-like Type I Terpene Synthases

Red algae or seaweeds produce highly distinctive halogenated terpenoid compounds, including the pentabromochlorinated monoterpene halomon that was once heralded as a promising anticancer agent. The first dedicated step in the biosynthesis of these natural product molecules is expected to be catalyzed by terpene synthase (TS) enzymes. Recent work has demonstrated an emerging class of type I TSs in red algal terpene biosynthesis. However, only one such enzyme from a notoriously haloterpenoid-producing red alga (Laurencia pacifica) has been functionally characterized and the product structure is not related to halogenated terpenoids. Herein, we report 10 new type I TSs from the red algae Portieria hornemannii, Plocamium pacificum, L. pacifica, and Laurencia subopposita that produce a diversity of halogenated mono- and sesquiterpenes. Here we used a combination of genome sequencing, terpenoid metabolomics, in vitro biochemistry, and bioinformatics to establish red algal TSs in all four species, including those associated with the selective production of key halogenated terpene precursors myrcene, trans-β-ocimene, and germacrene D-4-ol. These results expand on a small but growing number of characterized red algal TSs and offer insight into the biosynthesis of iconic halogenated algal compounds that are not without precedence elsewhere in biology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deep Learning on Multimodal Chemical and Whole Slide Imaging Data for Predicting Prostate Cancer Directly from Tissue Images

Prostate cancer is one of the most common cancers globally and is the second most common cancer in the male population in the US. Here we develop a study based on correlating the hematoxylin and eosin (H&E)-stained biopsy data with MALDI mass-spectrometric imaging data of the corresponding tissue to determine the cancerous regions and their unique chemical signatures and variations of the predicted regions with original pathological annotations. We obtain features from high-resolution optical micrographs of whole slide H&E stained data through deep learning and spatially register them with mass spectrometry imaging (MSI) data to correlate the chemical signature with the tissue anatomy of the data. We then use the learned correlation to predict prostate cancer from observed H&E images using trained coregistered MSI data. This multimodal approach can predict cancerous regions with ~80% accuracy, which indicates a correlation between optical H&E features and chemical information found in MSI. Further, we show that such paired multimodal data can be used for training feature extraction networks on H&E data which bypasses the need to acquire expensive MSI data and eliminates the need for manual annotation saving valuable time. Two chemical biomarkers were also found to be predicting the ground truth cancerous regions. This study shows promise in generating improved patient treatment trajectories by predicting prostate cancer directly from readily available H&E-stained biopsy images aided by coregistered MSI data.

60 APPLIED LIFE SCIENCES↗

EcoPLOT: dynamic analysis of biogeochemical data

Motivation: We have created EcoPLOT (parameterized linkage of omics-driven technologies), a web-app for the dynamic, interactive analysis of biogeochemical datasets that combines state-of-the-art analysis tools to statistically and graphically explore environmental, geochemical and microbiome datasets. Using the iterative random forest, a machine learning algorithm, EcoPLOT allows for the de novo discovery of drivers which exhibit significant impact on plant, microbial or soil dynamics. Availability and implementation: EcoPLOT is built entirely within the R language. It can be accessed through any system where R is installed, including Windows, Mac and most Linux systems. EcoPLOT is free to use and can be accessed at https://github.com/cdsanchez18/EcoPLOT.

59 BASIC BIOLOGICAL SCIENCES↗

Identifying Spectral Descriptors for Protonation in BaZr 0.8 Y 0.2 O 3-x with Electron Energy Loss Spectroscopy

Hydrogen economy is of paramount importance in the global transition to a sustainable, clean energy source that contributes to decarbonization efforts. In particular, the proton conducting proton ceramic fuel cells (PCFCs) play a crucial role in promoting hydrogen energy technology [1-5]. In a PCFC, the electrolyte is typically a solid oxide material that enables proton transport and can operate at temperatures lower than those of traditional oxygen ion conducting fuel cells [6]. The reduced operating temperature of PCFCs contributes to their durability, scalability, and efficiency [7,8]. BaZr 0.8 Y 0.2 O 3-x (BZY) is a promising proton conducting solid oxide electrolyte [9,10]. The emphasis on proton conductivity entails the importance of understanding proton content in the system. However, previous studies have largely relied on bulk techniques such as electrochemical impedance spectroscopy [7,11], Karl-Fischer titration [12,13], and thermogravimetric analysis [14] to obtain proton concentration. This is due to the small atomic size and light mass of hydrogen species making direct detection challenging. Bulk methods may be useful in estimating the proton content; however, it overlooks possible proton concentration gradient or segregation that may occur across or within a nanoparticle, especially when proton incorporation occurs nonuniformly through steam exposure on powder samples. In this paper, BZY is protonated at an estimated 0.15 mol of protons via steam exposure. Electron energy loss spectroscopy (EELS) with nanoscale spatial resolution is explored to identify proxy signals for proton detection using a JEOL ARM300 microscope operated at 300 kV with a Gatan GIF Continuum detector.

08 HYDROGEN↗

Structure analysis of the telomere resolvase from the Lyme disease spirochete Borrelia garinii reveals functional divergence of its C-terminal domain

Borrelia spirochetes are the causative agents of Lyme disease and relapsing fever, two of the most common tick-borne illnesses. A characteristic feature of these spirochetes is their highly segmented genomes which consists of a linear chromosome and a mixture of up to approximately 24 linear and circular extrachromosomal plasmids. The complexity of this genomic arrangement requires multiple strategies for efficient replication and partitioning during cell division, including the generation of hairpin ends found on linear replicons mediated by the essential enzyme ResT, a telomere resolvase. Using an integrative structural biology approach employing advanced modelling, circular dichroism, X-ray crystallography and small-angle X-ray scattering, we have generated high resolution structural data on ResT from B. garinii. Our data provides the first high-resolution structures of ResT from Borrelia spirochetes and revealed active site positioning in the catalytic domain. We also demonstrate that the C-terminal domain of ResT is required for both transesterification steps of telomere resolution, and is a requirement for DNA binding, distinguishing ResT from other telomere resolvases from phage and bacteria. These results advance our understanding of the molecular function of this essential enzyme involved in genome maintenance in Borrelia pathogens.

59 BASIC BIOLOGICAL SCIENCES↗

Construction and Modeling of a Coculture Microplate for Real-Time Measurement of Microbial Interactions

The dynamic structures of microbial communities emerge from the complex network of interactions between their constituent microorganisms. Quantitative measurements of these interactions are important for understanding and engineering ecosystem structure. Here, we present the development and application of the BioMe plate, a redesigned microplate device in which pairs of wells are separated by porous membranes. BioMe facilitates the measurement of dynamic microbial interactions and integrates easily with standard laboratory equipment. We first applied BioMe to recapitulate recently characterized, natural symbiotic interactions between bacteria isolated from the Drosophila melanogaster gut microbiome. Specifically, the BioMe plate allowed us to observe the benefit provided by two Lactobacillus strains to an Acetobacter strain. We next explored the use of BioMe to gain quantitative insight into the engineered obligate syntrophic interaction between a pair of Escherichia coli amino acid auxotrophs. We integrated experimental observations with a mechanistic computational model to quantify key parameters associated with this syntrophic interaction, including metabolite secretion and diffusion rates. This model also allowed us to explain the slow growth observed for auxotrophs growing in adjacent wells by demonstrating that, under the relevant range of parameters, local exchange between auxotrophs is essential for efficient growth. The BioMe plate provides a scalable and flexible approach for the study of dynamic microbial interactions.

3D printed device↗

Addressing delayed case reporting in infectious disease forecast modeling

Infectious disease forecasting is of great interest to the public health community and policymakers, since forecasts can provide insight into disease dynamics in the near future and inform interventions. Due to delays in case reporting, however, forecasting models may often underestimate the current and future disease burden. In this paper, we propose a general framework for addressing reporting delay in disease forecasting efforts with the goal of improving forecasts. We propose strategies for leveraging either historical data on case reporting or external internet-based data to estimate the amount of reporting error. We then describe several approaches for adapting general forecasting pipelines to account for under- or over-reporting of cases. We apply these methods to address reporting delay in data on dengue fever cases in Puerto Rico from 1990 to 2009 and to reports of influenza-like illness (ILI) in the United States between 2010 and 2019. Through a simulation study, we compare method performance and evaluate robustness to assumption violations. Our results show that forecasting accuracy and prediction coverage almost always increase when correction methods are implemented to address reporting delay. Some of these methods required knowledge about the reporting error or high quality external data, which may not always be available. Provided alternatives include excluding recently-reported data and performing sensitivity analysis. This work provides intuition and guidance for handling delay in disease case reporting and may serve as a useful resource to inform practical infectious disease forecasting efforts.

59 BASIC BIOLOGICAL SCIENCES↗

Quantification of early nonpharmaceutical interventions aimed at slowing transmission of Coronavirus Disease 2019 in the Navajo Nation and surrounding states (Arizona, Colorado, New Mexico, and Utah)

During an early period of the Coronavirus Disease 2019 (COVID-19) pandemic, the Navajo Nation, much like New York City, experienced a relatively high rate of disease transmission. Yet, between January and October 2020, it experienced only a single period of growth in new COVID-19 cases, which ended when cases peaked in May 2020. The daily number of new cases slowly decayed in the summer of 2020 until late September 2020. In contrast, the surrounding states of Arizona, Colorado, New Mexico, and Utah all experienced at least two periods of growth in the same time frame, with second surges beginning in late May to early June. Here, we investigated these differences in disease transmission dynamics with the objective of quantifying the contributions of non-pharmaceutical interventions (NPIs) (e.g., behaviors that limit disease transmission). We considered a compartmental model accounting for distinct periods of NPIs to analyze the epidemic in each of the five regions. We used Bayesian inference to estimate region-specific model parameters from regional surveillance data (daily reports of new COVID-19 cases) and to quantify uncertainty in parameter estimates and model predictions. Our results suggest that NPIs in the Navajo Nation were sustained over the period of interest, whereas in the surrounding states, NPIs were relaxed, which allowed for subsequent surges in cases. Our region-specific model parameterizations allow us to quantify the impacts of NPIs on disease incidence in the regions of interest.

60 APPLIED LIFE SCIENCES↗

ADEPT: A Pedagogical Framework for Integrating Agentic AI with Deterministic Scientific Workflows

The integration of Large Language Models (LLMs) into scientific research promises to accelerate discovery, yet a significant gap remains between the dynamic reasoning of Artificial Intelligence (AI) agents and the static, deterministic nature of canonical scientific workflows. This paper introduces ADEPT (Agentic Discovery and Exploration Platform for Tools), a reference architecture and pedagogical framework explicitly designed to bridge this gap. ADEPT's primary mission is to provide a transparent, "glass-box" environment where researchers and engineers can learn to effectively wrap established scientific software (e.g., BLAST, Nextflow pipelines) and compose it into reliable, agent-driven workflows. We describe its modular, multi-server architecture, which leverages the Model Context Protocol (MCP) for tool serving, LangGraph for robust agentic orchestration, and a secure nsjail-based sandbox for safe code execution. By prioritizing architectural clarity, safety, and modularity, ADEPT serves as an extensible blueprint for building trustworthy AI-augmented systems and fosters the collaborative development necessary to responsibly employ agentic AI for science. We provide practical examples of how to adapt and extend this framework, highlighting its utility in workforce development and AI-readiness capabilities across research and development projects.

97 MATHEMATICS AND COMPUTING↗

Carbon dioxide evolution rate as a method to monitor and control an aerobic biological waste treatment system

An experimental system was developed to study the microbial growth kinetic of an undefined mixed culture in an erobic biological waste treatment process. The experimental results were used to develop a mathematical model that can predict the performance of a bioreactor. The bioreactor will be used to regeneratively treat waste material which is expected to be generated during a long term manned space mission. Since the presence of insoluble particles in the chemically undefined complex media made estimating biomass very difficult in the real system, a clean system was devised to study the microbial growth from the soluble substrate.

Lee, S. S.↗

Nested Neural Networks

Report presents analysis of nested neural networks, consisting of interconnected subnetworks. Analysis based on simplified mathematical models more appropriate for artificial electronic neural networks, partly applicable to biological neural networks. Nested structure allows for retrieval of individual subpatterns. Requires fewer wires and connection devices than fully connected networks, and allows for local reconstruction of damaged subnetworks without rewiring entire network.

Baram, Yoram↗

Investigation of evolution-related aspects of bacterial rhodopsins

We have investigated evolution-related aspects of bacterial rhodopsins, the unique retinal-based energy transducing systems of halophilic archae. The approach was to describe both structural and functional aspects: the structure by sequencing genes to explore which regions are conserved, and the function by comparing proton and chloride transport in the closely related systems, bacteriorhodopsin and halorhodopsin, respectively. In the latter, we have made a good start toward the ultimate goal of separating the attributes of the general principles of retinal-based ionic pumps from those of the specific ion specificities, by determining the thermodynamics of the internal steps of the protein-mediated active transport process, as well as some of the intraprotein ion-transfer steps. Our present emphasis is on continuing to acquire the tools for studying what distinguishes proton transport from chloride transport. We consider it important, therefore, that we have been able to provide firm mathematical grounds for the kinetics analyses which underlies these studies. Our molecular biological studies have received a great boost from the expression vector for the bop gene based on a halobacterial plasmid, that we recently developed.

Source record↗

Effects of target fragmentation on evaluation of LET spectra from space radiations: implications for space radiation protection studies

We present calculations of linear energy transfer (LET) spectra in low earth orbit from galactic cosmic rays and trapped protons using the HZETRN/BRYNTRN computer code. The emphasis of our calculations is on the analysis of the effects of secondary nuclei produced through target fragmentation in the spacecraft shield or detectors. Recent improvements in the HZETRN/BRYNTRN radiation transport computer code are described. Calculations show that at large values of LET (> 100 keV/micrometer) the LET spectra seen in free space and low earth orbit (LEO) are dominated by target fragments and not the primary nuclei. Although the evaluation of microdosimetric spectra is not considered here, calculations of LET spectra support that the large lineal energy (y) events are dominated by the target fragments. Finally, we discuss the situation for interplanetary exposures to galactic cosmic rays and show that current radiation transport codes predict that in the region of high LET values the LET spectra at significant shield depths (> 10 g/cm2 of Al) is greatly modified by target fragments. These results suggest that studies of track structure and biological response of space radiation should place emphasis on short tracks of medium charge fragments produced in the human body by high energy protons and neutrons.

NASA Discipline Radiation Health↗

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

Abstract In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled 60 Co (Type M) and 131 I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method ( expm ). The differences generally observed between the implicit solvers and expm are attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

61 RADIATION PROTECTION AND DOSIMETRY↗

Environmental parameters of the Tennessee River in Alabama. 2: Physical, chemical, and biological parameters

Physical, chemical and biological water quality data from five sites in the Tennessee River, two in Guntersville Reservoir and three in Wheeler Reservoir were correlated with climatological data for three annual cycles. Two of the annual cycles are for the years prior to the Browns Ferry Nuclear Power Plant operations and one is for the first 14 months of Plant operations. A comparison of the results of the annual cycles indicates that two distinct physical conditions in the reservoirs occur, one during the warm months when the reservoirs are at capacity and one during the colder winter months when the reservoirs have been drawn-down for water storage during the rainy months and for weed control. The wide variations of physical and chemical parameters to which the biological organisms are subjected on an annual basis control the biological organisms and their population levels. A comparison of the parameters of the site below the Power plant indicates that the heated effluent from the plant operating with two of the three reactors has not had any effect on the organisms at this site. Recommendations given include the development of prediction mathematical models (statistical analysis) for the physical and chemical parameters under specific climatological conditions which affect biological organisms. Tabulated data of chemical analysis of water and organism populations studied is given.

Rosing, L. M.↗

Network Security via Biometric Recognition of Patterns of Gene Expression

Molecular biology provides the ability to implement forms of information and network security completely outside the bounds of legacy security protocols and algorithms. This paper addresses an approach which instantiates the power of gene expression for security. Molecular biology provides a rich source of gene expression and regulation mechanisms, which can be adopted to use in the information and electronic communication domains. Conventional security protocols are becoming increasingly vulnerable due to more intensive, highly capable attacks on the underlying mathematics of cryptography. Security protocols are being undermined by social engineering and substandard implementations by IT organizations. Molecular biology can provide countermeasures to these weak points with the current security approaches. Future advances in instruments for analyzing assays will also enable this protocol to advance from one of cryptographic algorithms to an integrated system of cryptographic algorithms and real-time expression and assay of gene expression products.

cybersecurity↗

Solving the sample size problem for resource selection functions

Abstract Sample size sufficiency is a critical consideration for estimating resource selection functions (RSFs) from GPS‐based animal telemetry. Cited thresholds for sufficiency include a number of captured animals and as many relocations per animal N as possible. These thresholds render many RSF‐based studies misleading if large sample sizes were truly insufficient, or unpublishable if small sample sizes were sufficient but failed to meet reviewer expectations. We provide the first comprehensive solution for RSF sample size by deriving closed‐form mathematical expressions for the number of animals M and the number of relocations per animal N required for model outputs to a given degree of precision. The sample sizes needed depend on just 3 biologically meaningful quantities: habitat selection strength, variation in individual selection and a novel measure of landscape complexity, which we define rigorously. The mathematical expressions are calculable for any environmental dataset at any spatial scale and are applicable to any study involving resource selection (including sessile organisms). We validate our analytical solutions using globally relevant empirical data including 5,678,623 GPS locations from 511 animals from 10 species (omnivores, carnivores and herbivores living in boreal, temperate and tropical forests, montane woodlands, swamps and Arctic tundra). Our analytic expressions show that the required M and N must decline with increasing selection strength and increasing landscape complexity, and this decline is insensitive to the definition of availability used in the analysis. Our results demonstrate that the most biologically relevant effects on the utilization distribution (i.e. those landscape conditions with the greatest absolute magnitude of resource selection) can often be estimated with much fewer than animals. We identify several critical steps in implementing these equations, including (a) a priori selection of expected model coefficients and (b) regular sampling of background (pseudoabsence) data within a given definition of availability. We discuss possible methods to identify a priori expectations for habitat selection coefficients, effects of scale on RSF estimation and caveats for rare species applications. We argue that these equations should be a mandatory component for all future RSF studies.

Street, Garrett M.↗

A python package for analyzing Resilience of Complex Systems (pyRoCS) v.0.0

SAND2024-01040O PyRoCS software synthesizes mathematical equations from several domains—including information theory, ecology, and engineering sciences—to support resilience analysis for complex systems. Resilience is the ability of the complex system being analyzed to withstand, operate through, and recover from a disruption. The complex system can be a physical system such as an electric grid, an organization such as a company, or even a subfunction of an organization. Existing mathematical equations for resilience analysis are found within multiple domains including information theory, biological sciences, and complex systems. This package synthesizes and refactors equations from these various domains to make them more generalizable for application across different types of complex systems relevant for resilience analysis. Users will be able to apply these equations to characterize different components of complex systems based on available data. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Verzi, Stephen↗