Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Models, Biological”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

SBML Level 3: an extensible format for the exchange and reuse of biological models

Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction-based models and packages that extend the core with features suited to other model types including constraint-based models, reaction-diffusion models, logical network models, and rule-based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single-cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Differential methods for assessing sensitivity in biological models

Differential sensitivity analysis is indispensable in fitting parameters, understanding uncertainty, and forecasting the results of both thought and lab experiments. Although there are many methods currently available for performing differential sensitivity analysis of biological models, it can be difficult to determine which method is best suited for a particular model. In this paper, we explain a variety of differential sensitivity methods and assess their value in some typical biological models. First, we explain the mathematical basis for three numerical methods: adjoint sensitivity analysis, complex perturbation sensitivity analysis, and forward mode sensitivity analysis. We then carry out four instructive case studies. (a) The CARRGO model for tumor-immune interaction highlights the additional information that differential sensitivity analysis provides beyond traditional naive sensitivity methods, (b) the deterministic SIR model demonstrates the value of using second-order sensitivity in refining model predictions, (c) the stochastic SIR model shows how differential sensitivity can be attacked in stochastic modeling, and (d) a discrete birth-death-migration model illustrates how the complex perturbation method of differential sensitivity can be generalized to a broader range of biological models. Finally, we compare the speed, accuracy, and ease of use of these methods. We find that forward mode automatic differentiation has the quickest computational time, while the complex perturbation method is the simplest to implement and the most generalizable.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Generative diffusion model surrogates for mechanistic agent-based biological models

Mechanistic, multicellular, agent-based models are commonly used to investigate tissue, organ, and organism-scale biology at single-cell resolution. The Cellular-Potts Model (CPM) is a powerful and popular framework for developing and interrogating these models. CPMs become computationally expensive at large space- and time- scales making application and investigation of developed models difficult. Surrogate models may allow for the accelerated evaluation of CPMs of complex biological systems. However, the stochastic nature of these models means each set of parameters may give rise to different model configurations, complicating surrogate model development. In this work, we leverage denoising diffusion probabilistic models (DDPMs) to train a generative AI surrogate of a CPM used to investigate in vitro vasculogenesis. We describe the use of an image classifier to learn the characteristics that define unique areas of a 2-dimensional parameter space. We then apply this classifier to aid in surrogate model selection and verification. Our CPM model surrogate generates model configurations 20,000 timesteps ahead of a reference configuration and demonstrates approximately a 22x reduction in computational time as compared to native code execution. Our work represents a step towards the implementation of DDPMs to develop digital twins of stochastic biological systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A validated thermal and biological model for predicting algal productivity in large scale outdoor cultivation systems

Proper assessment of the sustainability of algal products requires an understanding of algal growth rate with regional and temporal resolution. This study developed a bulk growth model that utilizes geospatial data, reactor geometry inputs, and a maximum of six strain-specific parameters to predict algal productivity with temporal and regional resolution. The model requires local weather data to calculate the time-resolved culture temperature which is combined with a biological model based on maximum theoretical carbon fixation to predict the time-resolved algal concentration. Validation of the model against experimental data illustrates an average ac- curacy of 0.27% ± 5.32% when comparing algal concentration at harvest over 45 days of cultivation data for Galdieria sulphuraria 5587.1, and Galdieria sulphuraria Soos grown in photobioreactors (PBRs) in Mesa, AZ (33.4152° N, 111.8315° W). For open raceway ponds (ORPs), an accuracy of 0.9% ± 2.35% was achieved when comparing algal concentration at harvest over 50 days of cultivation data for Chlorella vulgaris, Desmodesmus intermedius, Nannochloropsis oceanica grown in Atlanta, GA (33.7490° N, 84.3880° W) and Mesa, AZ (33.4152° N, 111.8315° W). The validation trials spanned three seasons (winter, spring, and summer) to ensure model ac- curacy throughout the year, and concentration between harvests was modeled with hourly resolution. The validated model is used to assess the productivity at several case-study locations in the United States with Chlorella vulgaris, grown in open raceway ponds and Galdieria sulphuraria grown in vertical flat panel photo- bioreactors. In conclusion, the model has been published as an open-source tool, in an effort to increase the fidelity of future studies that rely on outdoor algae cultivation for feedstock production and allow for time resolved results for future techno-economic and life cycle assessments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Using PyBioNetFit to leverage qualitative and quantitative data in biological model parameterization and uncertainty quantification

Data generated in studies of cellular regulatory systems are often qualitative. For example, measurements of signaling readouts in the presence and absence of mutations may reveal a rank ordering of responses across conditions but not the precise extents of mutation-induced differences. Qualitative data are often ignored by mathematical modelers or are considered in an ad hoc manner, as in the study of Kocieniewski and Lipniacki (2013) [Phys Biol 10: 035006], which was focused on the roles of MEK isoforms in ERK activation. In this earlier study, model parameter values were tuned manually to obtain consistency with a combination of qualitative and quantitative data. This approach is not reproducible, nor does it provide insights into parametric or prediction uncertainties. Here, starting from the same data and the same ordinary differential equation (ODE) model structure, we generate formalized statements of qualitative observations, making these observations more reusable, and we improve the model parameterization procedure by applying a systematic and automated approach enabled by the software package PyBioNetFit. We also demonstrate uncertainty quantification (UQ), which was absent in the original study. Our results show that PyBioNetFit enables qualitative data to be leveraged, together with quantitative data, in parameterization of systems biology models and facilitates UQ. These capabilities are important for reliable estimation of model parameters and model analyses in studies of cellular regulatory systems and reproducibility.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Parameter Estimation and Uncertainty Quantification for Systems Biology Models

Mathematical models can provide quantitative insights into immunoreceptor signaling, and other biological processes, but require parameterization and uncertainty quantification before reliable predictions become possible. We review currently available methods and software tools to address these problems. We consider gradient-based and gradient-free methods for point estimation of parameter values, and methods of profile likelihood, bootstrapping, and Bayesian inference for uncertainty quantification. We consider recent and potential future applications of these methods to systems-level modeling of immune-related phenomena.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Modeling biological nitrogen fixation in global natural terrestrial ecosystems

Biological nitrogen fixation plays an important role in the global nitrogen cycle. However, the fixation rate has been usually measured or estimated at a particular observational site. To quantify the fixation amount at the global scale, process-based models are needed. This study develops a biological nitrogen fixation model to quantitatively estimate the nitrogen fixation rate by plants in a natural environment. The revised nitrogen module better simulates the nitrogen cycle in comparison with our previous model that has not considered the fixation effects. The new model estimates that tropical forests have the highest fixation rate among all ecosystem types, which decreases from the Equator to the polar region. The estimated nitrogen fixation in global terrestrial ecosystems is 61.5 Tg N yr-1 with a range of 19.8–107.9 Tg N yr-1 in the 1990s. Our estimates are relatively low compared to some early estimates using empirical approaches but comparable to more recent estimates that involve more detailed processes in their modeling. Furthermore, the contribution of nitrogen made by biological nitrogen fixation depends on ecosystem type and climatic conditions. This study highlights that there are relatively large effects of biological nitrogen fixation on ecosystem nitrogen cycling. and the large uncertainty of the estimation calls for more comprehensive understanding of biological nitrogen fixation. More direct observational data for different ecosystems are in need to improve future quantification of fixation and its impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Predicting Performance of Macroalgae Farms with Hydrodynamic and Biological Modeling

As part of the ARPA-E MARINER program a set of simulation tools and findings were developed for the hydrodynamic, biological, and economic modeling of large scale offshore macroalgae farms. Results suggest the utility of the tools in understanding the complex interplay of design choices and environmental conditions on the structural loading and farm performance which drive the costs for macroalgae production.

09 BIOMASS FUELS↗
DOE OSTI2023

BMX: Biological modelling and interface exchange

Abstract High performance computing has a great potential to provide a range of significant benefits for investigating biological systems. These systems often present large modelling problems with many coupled subsystems, such as when studying colonies of bacteria cells. The aim to understand cell colonies has generated substantial interest as they can have strong economic and societal impacts through their roles in in industrial bioreactors and complex community structures, called biofilms, found in clinical settings. Investigating these communities through realistic models can rapidly exceed the capabilities of current serial software. Here, we introduce BMX, a software system developed for the high performance modelling of large cell communities by utilising GPU acceleration. BMX builds upon the AMRex adaptive mesh refinement package to efficiently model cell colony formation under realistic laboratory conditions. Using simple test scenarios with varying nutrient availability, we show that BMX is capable of correctly reproducing observed behavior of bacterial colonies on realistic time scales demonstrating a potential application of high performance computing to colony modelling. The open source software is available from the zenodo repository https://doi.org/10.5281/zenodo.8084270 under the BSD-2-Clause licence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Electromechanical memcapacitor model offering biologically plausible spiking.

In this article, we introduce a new nanoscale electromechanical device - a leaky memcapacitor - and show that it may be useful for the hardware implementation of spiking neurons. The leaky memcapacitor is a movableplate capacitor that becomes quite conductive when the plates come close to each other. The equivalent circuit of the leaky memcapacitor involves a memcapacitive and memristive system connected in parallel. In the leaky memcapacitor, resistance and capacitance depend on the same internal state variable, which is the displacement of the movable plate. We have performed a comprehensive analysis showing that several types of spiking observed in biological neurons can be implemented with the leaky memcapacitor. Significant attention is paid to the dynamic properties of the model. As in leaky memcapacitors the capacitive, leaking resistive, and reset functionalities are implemented naturally within the same device structure, their use will simplify the creation of spiking neural networks.

Zhang, Zixi↗
DOE OSTI2021

Following O and OH in He/O 2 and He/H 2 O gas mixtures—from the gas phase through the liquid phase to modifications on a biological sample

Abstract Applied cold atmospheric plasma allows for the controlled delivery of reactive oxygen and nitrogen species tailored for specific applications. Through the manipulation of the plasma parameters, feed gases, and careful consideration of the environment surrounding the treatment target, selective chemistries that preferentially influence the target can be produced and delivered. To demonstrate this, the COST reference microscale atmospheric pressure plasma jet is used to study the generation and transport of O and ⋅ OH from the gas phase through the liquid to the biological model target cysteine. Relative and absolute species densities of ⋅ OH and O are measured in the gas phase through laser induced fluorescence (LIF) and two-photon absorption LIF respectively. The transport of these species is followed into the liquid phase by hydrogen peroxide quantification and visualized by a fluorescence assay. Modifications to the model biological sample cysteine exposed to ⋅ OH and H 2 O 2 dominated chemistry (He/H 2 O (0.25%)) and O dominated chemistry (He/O 2 (0.6%)) is measured by FTIR spectroscopy. The origin of these species that modify cysteine is considered through the use of heavy water (H 2 18 O) and mass spectrometry. It is found that the reaction pathways differ significantly for He/O 2 and He/H 2 O. Hydrogen peroxide is formed mainly in the liquid phase in the presence of a substrate for He/O 2 whereas for He/H 2 O it forms in the gas phase. The liquid chemistry resulting from the He/O 2 admixture mainly targets the sulfur moiety of cysteine for oxidation up to irreversible oxidation states, while He/H 2 O treatment leads preferentially to reversible oxidation products. The more O or OH/H 2 O 2 dominated chemistry produced by the two gas admixtures studied offers the possibility to select species for target modification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Hypergraph Models of Biological Networks to Identify Genes Critical to Pathogenic Viral Response

Motivation: Representing biological networks as graphs is a powerful approach to reveal underlying patterns, signatures, and critical components from high-throughput biomolecular data. However, graphs do not natively capture the multi-way relationships present among genes and proteins in biological systems such as protein complexes, metabolic reactions, and signal transduction pathways. Hypergraphs are generalizations of graphs that naturally model multi-way interactions in data, and we therefore seek to understand how they can more faithfully identify, and potentially predict, complex relationships in genomic expression data sets. Results: We compiled a novel data set of transcriptional host response to pathogenic viral infections and formulated relationships between genes as a hypergraph where hyperedges are differentially expressed genes and vertices represent conditions. We find that hypergraph betweenness centrality is a superior method for identification of genes important to viral response when compared with graph centrality. Our results demonstrate the utility of using hypergraphs to represent complex biological systems, and highlight potentially interesting biological results about host response to highly pathogenic viruses.

systems biology, hypergraph, viral infection, biol↗
DOE OSTI2021

Biological Response Models: Predicting Injury and Mortality of Fish During Downstream Passage through Hydropower Facilities

This report describes the U.S. Department of Energy’s HydroPASSAGE project efforts to develop and collect biological response models for integration into the Biological Performance Assessment (BioPA) toolset and the Hydropower Biological Evaluation Toolset (HBET). These models help understand how fish are likely to respond during dam passage when exposed to hydraulic and physical stressors associated with turbines and other hydropower structures. When fish pass through hydropower facilities, they may encounter several stressors, of which the three most common are collision, rapid decompression, and fluid shear. Specialized equipment has been used to develop 99 biological response models for exposure to blade strike, rapid decompression, or fluid shear. The models were gathered from the literature or developed as part of this effort and include models for 31 different species of fish that have various predicted endpoints (i.e., injury or mortality). Among these models, considerable variation in susceptibility to the stressors has been observed from one species to another, and a species’ susceptibility to one stressor does not necessarily indicate similar susceptibility to another. Although several species have been examined, it is still unclear how many other species, which may have different morphological traits, may respond to these stressors, so further examination of the different species is needed. These models can and have been applied in several different cases, often using the BioPA toolset and HBET, to better understand the potential for injury and mortality that may occur during fish passage at hydropower facilities, including specific applications, such as turbine replacement, the installation of new turbines, or changing operations of currently installed turbines. As hydropower is continually developed to meet the electricity needs of society, tools such as HBET and BioPA, used with the integrated biological response models, will aid in the development of technologies and strategies that avoid, minimize, mitigate, or manage environmental effects.

13 HYDRO ENERGY↗