Search NASA⌕ Search

SEARCH · Search NASA

Results for “BIOLOGICAL MODEL”

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 OSTI2017

PyDREAM: high-dimensional parameter inference for biological models in python

Biological models contain many parameters whose values are difficult to measure directly via experimentation and therefore require calibration against experimental data. Markov chain Monte Carlo (MCMC) methods are suitable to estimate multivariate posterior model parameter distributions, but these methods may exhibit slow or premature convergence in high-dimensional search spaces. Here, we present PyDREAM, a Python implementation of the (Multiple-Try) Differential Evolution Adaptive Metropolis [DREAM(ZS)] algorithm developed by Vrugt and ter Braak (2008) and Laloy and Vrugt (2012). PyDREAM achieves excellent performance for complex, parameter-rich models and takes full advantage of distributed computing resources, facilitating parameter inference and uncertainty estimation of CPU-intensive biological models.

59 BASIC BIOLOGICAL SCIENCES↗
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↗
NASA NTRS1980

Evaluation of biological models using Spacelab

Biological models of hypogravity effects are described, including the cardiovascular-fluid shift, musculoskeletal, embryological and space sickness models. These models predict such effects as loss of extracellular fluid and electrolytes, decrease in red blood cell mass, and the loss of muscle and bone mass in weight-bearing portions of the body. Experimentation in Spacelab by the use of implanted electromagnetic flow probes, by fertilizing frog eggs in hypogravity and fixing the eggs at various stages of early development and by assessing the role of the vestibulocular reflex arc in space sickness is suggested. It is concluded that the use of small animals eliminates the uncertainties caused by corrective or preventive measures employed with human subjects.

Tollinger, D.↗
DOE OSTI2018

Using both qualitative and quantitative data in parameter identification for systems biology models

Abstract In systems biology, qualitative data are often generated, but rarely used to parameterize models. We demonstrate an approach in which qualitative and quantitative data can be combined for parameter identification. In this approach, qualitative data are converted into inequality constraints imposed on the outputs of the model. These inequalities are used along with quantitative data points to construct a single scalar objective function that accounts for both datasets. To illustrate the approach, we estimate parameters for a simple model describing Raf activation. We then apply the technique to a more elaborate model characterizing cell cycle regulation in yeast. We incorporate both quantitative time courses (561 data points) and qualitative phenotypes of 119 mutant yeast strains (1647 inequalities) to perform automated identification of 153 model parameters. We quantify parameter uncertainty using a profile likelihood approach. Our results indicate the value of combining qualitative and quantitative data to parameterize systems biology models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

TH-A-19A-02: Expanding TOPAS Towards Biological Modeling

Purpose: To expand Monte Carlo simulations of TOPAS to directly calculate radiobiological effects using various biological models. Methods: We use the TOPAS toolkit version beta9 and expand the provided scoring system to directly include scorers for relative biological effect (RBE) modeling. A scorer is a class to record a given physical property in the Monte Carlo simulations. We developed a scorer for the dose averaged linear energy transfer (LETd). This scorer was used to compare RBE weighted dose distributions assuming a constant RBE of 1.1 with the RBE weighted dose using a variable RBE based on a model using LETd and dose as input. Additionally, the local effect model (LEM) was implemented in TOPAS to estimate RBE weighted dose distributions. Survival curves for the LEM were simulated for a monoenergetic proton beams impinging on phantoms and for full patient treatment plans. Results: We have implemented radiobiological effect models in TOPAS, providing users with the possibility to directly obtain parameters related to radiobiological dose and use those as additional information to assess the quality of treatment plans. We applied the scorers to clinical patient treatments and phantom simulations. We find that for the proton therapy plans investigated the regions with the highest LETd are generally outside of the target volume, coinciding with the distal edge of treatment fields. LETd-based variable RBE weighted dose models predict higher doses at the distal edge of over 10% of the prescribed dose than when assuming a constant RBE. Conclusion: By implementing radiobiological effect models into a Monte Carlo simulation package we combine accurate dose calculation with biological modeling. TOPAS thereby provides distributions of multiple physical quantities such as the dose and LETd, and estimates of biological parameters (RBE, survival) thus allowing direct comparison of biological effective dose predictions between different models.

60 APPLIED LIFE 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 OSTI2014

SU-E-T-146: Effects of Uncertainties of Radiation Sensitivity of Biological Modelling for Treatment Planning

Purpose: The aim of this study was to evaluate the distribution of uncertainty of cell survival by radiation, and assesses the usefulness of stochastic biological model applying for gaussian distribution. Methods: For single cell experiments, exponentially growing cells were harvested from the standard cell culture dishes by trypsinization, and suspended in test tubes containing 1 ml of MEM(2x10{sup 6} cells/ml). The hypoxic cultures were treated with 95% N{sub 2}−5% CO{sub 2} gas for 30 minutes. In vitro radiosensitization was also measured in EMT6/KU single cells to add radiosensitizer under hypoxic conditions. X-ray irradiation was carried out by using an Xray unit (Hitachi X-ray unit, model MBR-1505R3) with 0.5 mm Al/1.0 mm Cu filter, 150 kV, 4 Gy/min). In vitro assay, cells on the dish were irradiated with 1 Gy to 24 Gy, respectively. After irradiation, colony formation assays were performed. Variations of biological parameters were investigated at standard cell culture(n=16), hypoxic cell culture(n=45) and hypoxic cell culture(n=21) with radiosensitizers, respectively. The data were obtained by separate schedule to take account for the variation of radiation sensitivity of cell cycle. Results: At standard cell culture, hypoxic cell culture and hypoxic cell culture with radiosensitizers, median and standard deviation of alpha/beta ratio were 37.1±73.4 Gy, 9.8±23.7 Gy, 20.7±21.9 Gy, respectively. Average and standard deviation of D{sub 50} were 2.5±2.5 Gy, 6.1±2.2 Gy, 3.6±1.3 Gy, respectively. Conclusion: In this study, we have challenged to apply these uncertainties of parameters for the biological model. The variation of alpha values, beta values, D{sub 50} as well as cell culture might have highly affected by probability of cell death. Further research is in progress for precise prediction of the cell death as well as tumor control probability for treatment planning.

60 APPLIED LIFE 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↗
NASA NTRS1990

Coupling of Coastal Zone Color Scanner data to a physical-biological model of the southeastern U.S. continental shelf ecosystem. I - CZCS data description and Lagrangian particle tracing experiments. II - An Eulerian model. III - Nutrient and phytoplankton fluxes and CZCS data assimilation

Surface phytoplankton biomass of the southeastern U.S. continental shelf area is discussed based on coastal zone color scanner (CZCS) images obtained in April 1980. Data of chlorophyll distributions are analyzed in conjunction with concurrent flow and temperature fields. Lagrangian particle tracing experiments show that the particles move consistently with the evolution of the chlorophyll patterns. A four-component physical-biological model for a horizontal plane at a nominal depth of 17 m is presented. Model simulations using various physical-biological dynamics and boundary conditions show that the variability of chlorophyll distributions is controlled by horizontal advection. Phytoplankton and nutrient fluxes, calculated using the model, show considerable variability with time. The chlorophyll distributions obtained from the CZCS images are assimilated into the model to improve the phytoplankton flux estimates.

Ishizaka, Joji↗
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 OSTI2015

Fast Biological Modeling for Voxel-based Heavy Ion Treatment Planning Using the Mechanistic Repair-Misrepair-Fixation Model and Nuclear Fragment Spectra

Purpose: The physical and biological differences between heavy ions and photons have not been fully exploited and could improve treatment outcomes. In carbon ion therapy, treatment planning must account for physical properties, such as the absorbed dose and nuclear fragmentation, and for differences in the relative biological effectiveness (RBE) of ions compared with photons. We combined the mechanistic repair-misrepair-fixation (RMF) model with Monte Carlo-generated fragmentation spectra for biological optimization of carbon ion treatment plans. Methods and Materials: Relative changes in double-strand break yields and radiosensitivity parameters with particle type and energy were determined using the independently benchmarked Monte Carlo damage simulation and the RMF model to estimate the RBE values for primary carbon ions and secondary fragments. Depth-dependent energy spectra were generated with the Monte Carlo code FLUKA for clinically relevant initial carbon ion energies. The predicted trends in RBE were compared with the published experimental data. Biological optimization for carbon ions was implemented in a 3-dimensional research treatment planning tool. Results: We compared the RBE and RBE-weighted dose (RWD) distributions of different carbon ion treatment scenarios with and without nuclear fragments. The inclusion of fragments in the simulations led to smaller RBE predictions. A validation of RMF against measured cell survival data reported in published studies showed reasonable agreement. We calculated and optimized the RWD distributions on patient data and compared the RMF predictions with those from other biological models. The RBE values in an astrocytoma tumor ranged from 2.2 to 4.9 (mean 2.8) for a RWD of 3 Gy(RBE) assuming (α/β){sub X} = 2 Gy. Conclusions: These studies provide new information to quantify and assess uncertainties in the clinically relevant RBE values for carbon ion therapy based on biophysical mechanisms. We present results from the first biological optimization of carbon ion radiation therapy beams on patient data using a combined RMF and Monte Carlo damage simulation modeling approach. The presented method is advantageous for fast biological optimization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
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 OSTI2014

WE-G-BRE-02: Biological Modeling of Gold Nanoparticle Radiosensitization for Proton Therapy

Purpose: The aim of this work is to investigate the radiosensitization effect of gold nanoparticles (GNP) in a proton beam. A computational model was built using the Local Effect Model (LEM) to predict the biological outcome of gold nanoparticle (GNP) sensitization. We present the results using a clinical proton beam, 6MV photon beam and two kilovoltage photon beams. Methods: First, Monte Carlo simulations were carried out using TOPAS (TOol for PArticle Simulation) to obtain the spatial dose distribution in the vicinity of GNPs. The dose distribution was then used as an input for LEM, which predicts dose-response curves for high linear energy transfer radiation using the track structure. The cell survival curves were evaluated for three particle sources (proton beam, MV photon beam and kV photon beam), various treatment depths for each particle source, various GNP uptakes and two different GNP sizes. Results: For proton therapy, the GNP sensitization effect is highly dependent on the treatment depth due to the energy-dependent interaction probability. We predict that if GNPs can be taken up by the cell nucleus, proton therapy can be significantly enhanced. If GNPs are only internalized into the cytoplasm, proton therapy can still be enhanced by GNPs and if GNPs are not internalized into cells, there will be no direct damage to the nucleus. For the same GNP uptake and concentration, the cell survival at 2Gy is reduced by 80% using kilovoltage photons, 50% using protons and only 2% using clinical MV photons. Finally, for the same weight of GNPs taken up by the cells, 10 nm GNPs causes 3 times more damage than 50 nm GNPs. Conclusion: We showed that GNPs have potential to be used to enhance radiation therapy for clinical proton beams.

07 ISOTOPE AND RADIATION SOURCES↗
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↗