Search NASA⌕ Search

SEARCH · Search NASA

Results for “biophysical simulations”

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 37 records · Page 2

Integrated Biophysical Modeling and Image Analysis: Application to Neuro-Oncology

Central nervous system (CNS) tumors come with vastly heterogeneous histologic, molecular, and radiographic landscapes, rendering their precise characterization challenging. The rapidly growing fields of biophysical modeling and radiomics have shown promise in better characterizing the molecular, spatial, and temporal heterogeneity of tumors. Integrative analysis of CNS tumors, including clinically acquired multi-parametric magnetic resonance imaging (mpMRI) and the inverse problem of calibrating biophysical models to mpMRI data, assists in identifying macroscopic quantifiable tumor patterns of invasion and proliferation, potentially leading to improved (a) detection/segmentation of tumor subregions and (b) computer-aided diagnostic/prognostic/predictive modeling. This article presents a summary of (a) biophysical growth modeling and simulation,(b) inverse problems for model calibration, (c) these models' integration with imaging workflows, and (d) their application to clinically relevant studies. We anticipate that such quantitative integrative analysis may even be beneficial in a future revision of the World Health Organization (WHO) classification for CNS tumors, ultimately improving patient survival prospects.

60 APPLIED LIFE SCIENCES↗

Machine Learning‐Augmented Molecular Dynamics Simulations (MD) Reveal Insights Into the Disconnect Between Affinity and Activation of ZTP Riboswitch Ligands

Abstract The challenge of targeting RNA with small molecules necessitates a better understanding of RNA–ligand interaction mechanisms. However, the dynamic nature of nucleic acids, their ligand‐induced stabilization, and how conformational changes influence gene expression pose significant difficulties for experimental investigation. This work employs a combination of computational and experimental methods to address these challenges. By integrating structure‐informed design, crystallography, and machine learning‐augmented all‐atom molecular dynamics simulations (MD), we synthesized, biophysically and biochemically characterized, and studied the dissociation of a library of small molecule activators of the 5‐aminoimidazole–4–carboxamide ribonucleotide triphosphate (ZTP) riboswitch, a ligand‐binding RNA motif that regulates bacterial gene expression. We uncovered key interaction mechanisms, revealing valuable insights into the role of ligand binding kinetics on riboswitch activation. Further, we established that ligand on‐rates determine activation potency as opposed to binding affinity and elucidated RNA structural differences, which provide mechanistic insights into the interplay of RNA structure on riboswitch activation.

Chemistry↗

Quantitative cross-species translators of cardiac myocyte electrophysiology: Model training, experimental validation, and applications

Animal experimentation is key in the evaluation of cardiac efficacy and safety of novel therapeutic compounds. However, interspecies differences in the mechanisms regulating excitation-contraction coupling can limit the translation of experimental findings from animal models to human physiology and undermine the assessment of drugs’ efficacy and safety. Here, we built a suite of translators for quantitatively mapping electrophysiological responses in ventricular myocytes across species. We trained these statistical operators using a broad dataset obtained by simulating populations of our biophysically detailed computational models of action potential and Ca 2+ transient in mouse, rabbit, and human. We then tested our translators against experimental data describing the response to stimuli, such as ion channel block, change in beating rate, and β-adrenergic challenge. We demonstrate that this approach is well suited to predicting the effects of perturbations across different species or experimental conditions and suggest its integration into mechanistic studies and drug development pipelines.

60 APPLIED LIFE SCIENCES↗

Multi-Factor Impact Analysis of Agricultural Production in Bangladesh with Climate Change

Diverse vulnerabilities of Bangladesh's agricultural sector in 16 sub-regions are assessed using experiments designed to investigate climate impact factors in isolation and in combination. Climate information from a suite of global climate models (GCMs) is used to drive models assessing the agricultural impact of changes in temperature, precipitation, carbon dioxide concentrations, river floods, and sea level rise for the 2040-2069 period in comparison to a historical baseline. Using the multi-factor impacts analysis framework developed in Yu et al. (2010), this study provides new sub-regional vulnerability analyses and quantifies key uncertainties in climate and production. Rice (aman, boro, and aus seasons) and wheat production are simulated in each sub-region using the biophysical Crop Environment REsource Synthesis (CERES) models. These simulations are then combined with the MIKE BASIN hydrologic model for river floods in the Ganges-Brahmaputra-Meghna (GBM) Basins, and the MIKE21Two-Dimensional Estuary Model to determine coastal inundation under conditions of higher mean sea level. The impacts of each factor depend on GCM configurations, emissions pathways, sub-regions, and particular seasons and crops. Temperature increases generally reduce production across all scenarios. Precipitation changes can have either a positive or a negative impact, with a high degree of uncertainty across GCMs. Carbon dioxide impacts on crop production are positive and depend on the emissions pathway. Increasing river flood areas reduce production in affected sub-regions. Precipitation uncertainties from different GCMs and emissions scenarios are reduced when integrated across the large GBM Basins' hydrology. Agriculture in Southern Bangladesh is severely affected by sea level rise even when cyclonic surges are not fully considered, with impacts increasing under the higher emissions scenario.

Ruane, Alex C.↗

Investigating the Diurnal Radiative, Turbulent, and Biophysical Processes in the Amazonian Canopy‐Atmosphere Interface by Combining LES Simulations and Observations

Abstract We investigate the atmospheric diurnal variability inside and above the Amazonian rainforest for a representative day during the dry season. To this end, we combine high‐resolution large‐eddy simulations that are constrained and evaluated against a comprehensive observation set, including CO 2 concentrations, gathered during GoAmazon2014/15. We design systematic numerical experiments to quantify whether a multilayer approach in solving the explicit canopy improves our canopy‐atmosphere representation. We particularly focus on the relationship between photosynthesis and plant transpiration, and their distribution at leaf and canopy scales. We found the variability of photosynthesis drivers like vapor pressure deficit and leaf temperature to be about 3 times larger for sunlit leaves compared to shaded leaves. This leads to a large spread on leaf stomatal conductance values with minimum and maximum values varying more than 100%. Regarding the turbulent structure, we find wind‐driven stripe‐like shapes at the canopy top and structures resembling convective cells at the canopy. Wind‐related variables provide the best spatiotemporal agreement between model and observations. The potential temperature and heat flux profiles agree with an observed decoupling near the canopy top interface, although with less variability and cold biases of up to 3 K. The increasing complexity on the biophysical processes leads to the largest disagreements for evaporation, CO 2 plant assimilation and soil efflux. The model is able to capture the correct dependences and trends with the magnitudes still differing. We finally discuss the need to revise leaf and soil models and to complete the observations at leaf and canopy levels.

54 ENVIRONMENTAL SCIENCES↗

Modeling Bird Migration under Climate Change: A Mechanistic Approach

How will migrating birds respond to changes in the environment under climate change? What are the implications for migratory success under the various accelerated climate change scenarios as forecast by the Intergovernmental Panel on Climate Change? How will reductions or increased variability in the number or quality of wetland stop-over sites affect migratory bird species? The answers to these questions have important ramifications for conservation biology and wildlife management. Here, we describe the use of continental scale simulation modeling to explore how spatio-temporal changes along migratory flyways affect en-route migration success. We use an individually based, biophysical, mechanistic, bird migration model to simulate the movement of shorebirds in North America as a tool to study how such factors as drought and wetland loss may impact migratory success and modify migration patterns. Our model is driven by remote sensing and climate data and incorporates important landscape variables. The energy budget components of the model include resting, foraging, and flight, but presently predation is ignored. Results/Conclusions We illustrate our model by studying the spring migration of sandpipers through the Great Plains to their Arctic breeding grounds. Why many species of shorebirds have shown significant declines remains a puzzle. Shorebirds are sensitive to stop-over quality and spacing because of their need for frequent refueling stops and their opportunistic feeding patterns. We predict bird "hydrographs that is, stop-over frequency with latitude, that are in agreement with the literature. Mean stop-over durations predicted from our model for nominal cases also are consistent with the limited, but available data. For the shorebird species simulated, our model predicts that shorebirds exhibit significant plasticity and are able to shift their migration patterns in response to changing drought conditions. However, the question remains as to whether this behavior can be maintained over increasing and sustained environmental change. Also, the problem is much more complex than described by the current processes captured in our model. We have taken some important and interesting steps, and our model does demonstrate how local scale information about individual stop-over sites can be linked into the migratory flyway as a whole. We are incorporating additional, species specific, mechanistic processes to better reflect different climate change scenarios

Smith, James A.↗

Insights into the dual cleavage activity of the GH16 laminarinase enzyme class on β-1,3 and β-1,4 glycosidic bonds

Glycoside hydrolases (GH) are involved in the degradation of a wide diversity of carbohydrates and present several biotechnological applications. Many GH families are composed of enzymes with a single well-defined specificity. In contrast, enzymes from the GH16 family can act on a range of different polysaccharides, including β-glucans and galactans. SCLam, a GH16 member derived from a soil metagenome, an endo-β-1,3(4)-glucanase (EC 3.2.1.6), can cleave both β-1,3 and β-1,4 glycosidic bonds in glucans, such as laminarin, barley β-glucan, and cello-oligosaccharides. A similar cleavage pattern was previously reported for other GH16 family members. However, the molecular mechanisms for this dual cleavage activity on (1,3)- and (1,4)-β-D-glycosidic bonds by laminarinases have not been elucidated. In this sense, we determined the X-ray structure of a presumably inactive form of SCLam co-crystallized with different oligosaccharides. The solved structures revealed general bound products that are formed due to residual activities of hydrolysis and transglycosylation. Biochemical and biophysical analyses and molecular dynamics simulations help to rationalize differences in activity towards different substrates. Our results depicted a bulky aromatic residue near the catalytic site critical to select the preferable configuration of glycosidic bonds in the binding cleft. Altogether, these data contribute to understanding the structural basis of recognition and hydrolysis of β-1,3 and β-1,4 glycosidic linkages of the laminarinase enzyme class, which is valuable for future studies on the GH16 family members and applications related to biomass conversion into feedstocks and bioproducts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Contrasting the Biophysical and Radiative Effects of Rising CO 2 Concentrations on Ozone Dry Deposition Fluxes

Abstract The dry deposition of ozone from the atmosphere to ecosystems is an important coupling mechanism between atmospheric chemistry and terrestrial biogeochemical processes. In most Earth system models, dry deposition is simulated using a resistor‐in‐series approach that aims to parameterize the governing biological, chemical, and physical processes through a series of functional approximations. Here, we evaluate the influence of carbon cycle‐climate responses on this parameterization using the results of the Energy Exascale Earth System Model v1.1 Biogeochemistry simulation campaign. This simulation campaign was designed in part to explore the biophysical and radiative effects of rising historical CO 2 concentrations on the Earth system. We find that while the global annual ozone dry deposition is relatively insensitive to these effects, regionally the influence can be up to 10%. The strongest regional sensitivities in ozone dry deposition are predominantly in higher latitudes over land in the northern hemisphere and are dominated by the radiative effect of CO 2 , with little net influence of biophysical responses. Of all the impacts of the radiative effect of CO 2 , we point to the potential importance of accurately representing ozone deposition to snow in Earth System Models and provide recommendations for future simulation campaigns.

54 ENVIRONMENTAL SCIENCES↗

A Stochastic Model of Space Radiation Transport as a Tool in the Development of Time-Dependent Risk Assessment

A new computer model, the GCR Event-based Risk Model code (GERMcode), was developed to describe biophysical events from high-energy protons and heavy ions that have been studied at the NASA Space Radiation Laboratory (NSRL) [1] for the purpose of simulating space radiation biological effects. In the GERMcode, the biophysical description of the passage of heavy ions in tissue and shielding materials is made with a stochastic approach that includes both ion track structure and nuclear interactions. The GERMcode accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model [2]. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections

Kim, Myung-Hee Y.↗

GERMcode: A Stochastic Model for Space Radiation Risk Assessment

A new computer model, the GCR Event-based Risk Model code (GERMcode), was developed to describe biophysical events from high-energy protons and high charge and energy (HZE) particles that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the GERMcode, the biophysical description of the passage of HZE particles in tissue and shielding materials is made with a stochastic approach that includes both particle track structure and nuclear interactions. The GERMcode accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections. For NSRL applications, the GERMcode evaluates a set of biophysical properties, such as the Poisson distribution of particles or delta-ray hits for a given cellular area and particle dose, the radial dose on tissue, and the frequency distribution of energy deposition in a DNA volume. By utilizing the ProE/Fishbowl ray-tracing analysis, the GERMcode will be used as a bi-directional radiation transport model for future spacecraft shielding analysis in support of Mars mission risk assessments. Recent radiobiological experiments suggest the need for new approaches to risk assessment that include time-dependent biological events due to the signaling times for activation and relaxation of biological processes in cells and tissue. Thus, the tracking of the temporal and spatial distribution of events in tissue is a major goal of the GERMcode in support of the simulation of biological processes important in GCR risk assessments. In order to validate our approach, basic radiobiological responses such as cell survival curves, mutation, chromosomal aberrations, and representative mouse tumor induction curves are implemented into the GERMcode. Extension of these descriptions to other endpoints related to non-targeted effects and biochemical pathway responses will be discussed.

Kim, Myung-Hee Y.↗

3D Visualization of Monte-Carlo Simulation's of HZE Track Structure and Initial Chemical Species

Heavy ions biophysics is important for space radiation risk assessment [1] and hadron-therapy [2]. The characteristic of heavy ions tracks include a very high energy deposition region close to the track (<20 nm) denoted as the track core, and an outer penumbra region consisting of individual secondary electrons (6-rays). A still open question is the radiobiological effects of 6- rays relative to the track core. Of importance is the induction of double-strand breaks (DSB) [3] and oxidative damage to the biomolecules and the tissue matrix, considered the most important lesions for acute and long term effects of radiation. In this work, we have simulated a 56Fe26+ ion track of 1 GeV/amu with our Monte-Carlo code RITRACKS [4]. The simulation results have been used to calculate the energy depiction and initial chemical species in a "voxelized" space, which is then visualized in 3D. Several voxels with dose >1000 Gy are found in the penumbra, some located ~0.1 mm from the track core. In computational models, the DSB induction probability is calculated with radial dose [6], which may not take into account the higher RBE of electron track ends for DSB induction. Therefore, these simulations should help improve models of DSB induction and our understanding of heavy ions biophysics.

Plante, Ianik↗

Extended ecosystem signatures with application to Eos synergism requirements

The primary objective is to define the advantages of synergistically combining optical and microwave remote sensing measurements for the determination of biophysical properties important in ecosystem modeling. This objective was approached in a stepwise fashion starting with ground-based observations of controlled agricultural and orchard canopies and progressing to airborne observations of more natural forest ecosystems. This observational program is complemented by a parallel effort to model the visible reflectance and microwave scattering properties of composite vegetation canopies. The goals of the modeling studies are to verify our basic understanding of the sensor-scene interaction physics and to provide the basis for development of inverse models optimized for retrieval of key biophysical properties. These retrieval algorithms can then be used to simulate the expected performance of various aspects of Eos including the need for simultaneous SAR and HIRIS observations or justification for other (non-synchronous) relative timing constraints and the frequency, polarization, and angle of incidence requirements for accurate biophysical parameter extractions. This program completed a very successful series of truck-mounted experiments, made remarkable progress in development and validation of optical reflectance and microwave scattering models for vegetation, extended the scattering models to accommodate discontinuous and periodic canopies, developed inversion approaches for surface and canopy properties, and disseminated these results widely through symposia and journal publications. In addition, the third generation of the computer code for the microwave scattering models was provided to a number of other US, Canadian, Australian, and European investigators who are currently presenting and publishing results using the MIMICS research code.

Ulaby, Fawwaz T.↗

Biomembrane Structure and Material Properties Studied With Neutron Scattering

Cell membranes and their associated structures are dynamical supramolecular structures where different physiological processes take place. Detailed knowledge of their static and dynamic structures is therefore needed, to better understand membrane biology. The structure–function relationship is a basic tenet in biology and has been pursued using a range of different experimental approaches. In this review, we will discuss one approach, namely the use of neutron scattering techniques as applied, primarily, to model membrane systems composed of lipid bilayers. An advantage of neutron scattering, compared to other scattering techniques, is the differential sensitivity of neutrons to isotopes of hydrogen and, as a result, the relative ease of altering sample contrast by substituting protium for deuterium. This property makes neutrons an ideal probe for the study of hydrogen-rich materials, such as biomembranes. In this review article, we describe isotopic labeling studies of model and viable membranes, and discuss novel applications of neutron contrast variation in order to gain unique insights into the structure, dynamics, and molecular interactions of biological membranes. We specifically focus on how small-angle neutron scattering data is modeled using different contrast data and molecular dynamics simulations. We also briefly discuss neutron reflectometry and present a few recent advances that have taken place in neutron spin echo spectroscopy studies and the unique membrane mechanical data that can be derived from them, primarily due to new models used to fit the data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insight into Molecular Basis and Dynamics of Full-length CRaf Kinase in Cellular Signaling Mechanism

This study presents the first large-scale simulation using an initial structure predicted by AI/ML algorithms for the 648-amino-acid CRaf kinase, which plays a key role in cellular signaling. Simulation results show the evolution of the predicted structure into much more compact structures with inter-domain interactions that shed insights into auto-inhibition mechanism, paradoxical effect, activation, and recruitment pathways in the CRaf kinase. Newly identified epitopes in the CRaf may suggest additional drug targets. The results were published in Biophysical Journal, DOI:10.1016/j.bpj.2024.06.028.

59 BASIC BIOLOGICAL SCIENCES↗

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗

Herbicide‐resistant weed management with robots: A weed ecological–economic model

The heavy reliance on herbicides for weed control has led to an increase in resistant weeds in the United States. Robotic weed control is emerging as an alternative technology for removing weeds mechanically using artificial intelligence. We develop an integrated weed ecological and economic dynamic (I‐WEED) model to examine the biophysical and economic drivers of adopting robotic weed management and simulate the optimal timing and intensity of robotic adoption within and across growing seasons. We specify a cohort‐based weed growth model that relates yield damages to effective weed density and treats the susceptibility of weeds to herbicides as a renewable resource that can be regenerated by using mechanical weeding robots, due to a fitness cost that makes resistant weeds less prolific. Compared to myopic weed management which ignores resistance development, forward‐looking management leads to earlier adoption of robots and treating robots as complements instead of substitutes to herbicides. This weed management results in adopting fewer robots, deploying robots on a smaller portion of the land, higher profitability, and lower yield loss in the long run, relative to myopic management. Counterintuitively, myopic management leads to a lower resistance level through its higher robot adoption intensity. We also find that a lower level of initial weed seed resistance and/or a higher fitness cost result in a higher level of resistance because they create incentives for farmers to delay the adoption of robotic weed control. Our analysis shows the importance of jointly considering the interactions between weed ecology and economics in analyzing the incentives and effects of robotic weed management on weed resistance.

agricultural robotics↗

Assessing impacts of selective logging on water, energy, and carbon budgets and ecosystem dynamics in Amazon forests using the Functionally Assembled Terrestrial Ecosystem Simulator

Abstract. Tropical forest degradation from logging, fire, and fragmentation not only alters carbon stocks and carbon fluxes, but also impacts physical land surface properties such as albedo and roughness length. Such impacts are poorly quantified to date due to difficulties in accessing and maintaining observational infrastructures, as well as the lack of proper modeling tools for capturing the interactions among biophysical properties, ecosystem demography, canopy structure, and biogeochemical cycling in tropical forests. As a first step to address these limitations, we implemented a selective logging module into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) by mimicking the ecological, biophysical, and biogeochemical processes following a logging event. The model can specify the timing and aerial extent of logging events, splitting the logged forest patch into disturbed and intact patches; determine the survivorship of cohorts in the disturbed patch; and modifying the biomass and necromass (total mass of coarse woody debris and litter) pools following logging. We parameterized the logging module to reproduce a selective logging experiment at the Tapajós National Forest in Brazil and benchmarked model outputs against available field measurements. Our results suggest that the model permits the coexistence of early and late successional functional types and realistically characterizes the seasonality of water and carbon fluxes and stocks, the forest structure and composition, and the ecosystem succession following disturbance. However, the current version of FATES overestimates water stress in the dry season and therefore fails to capture seasonal variation in latent and sensible heat fluxes. Moreover, we observed a bias towards low stem density and leaf area when compared to observations, suggesting that improvements are needed in both carbon allocation and establishment of trees. The effects of logging were assessed by different logging scenarios to represent reduced impact and conventional logging practices, both with high and low logging intensities. The model simulations suggest that in comparison to old-growth forests the logged forests rapidly recover water and energy fluxes in 1 to 3 years. In contrast, the recovery times for carbon stocks, forest structure, and composition are more than 30 years depending on logging practices and intensity. This study lays the foundation to simulate land use change and forest degradation in FATES, which will be an effective tool to directly represent forest management practices and regeneration in the context of Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Modeling bacterial microcompartment architectures for enhanced cyanobacterial carbon fixation

The carboxysome is a bacterial microcompartment (BMC) which plays a central role in the cyanobacterial CO 2 -concentrating mechanism. These proteinaceous structures consist of an outer protein shell that partitions Rubisco and carbonic anhydrase from the rest of the cytosol, thereby providing a favorable microenvironment that enhances carbon fixation. The modular nature of carboxysomal architectures makes them attractive for a variety of biotechnological applications such as carbon capture and utilization. In silico approaches, such as molecular dynamics (MD) simulations, can support future carboxysome redesign efforts by providing new spatio-temporal insights on their structure and function beyond in vivo experimental limitations. However, specific computational studies on carboxysomes are limited. Fortunately, all BMC (including the carboxysome) are highly structurally conserved which allows for practical inferences to be made between classes. Here, we review simulations on BMC architectures which shed light on (1) permeation events through the shell and (2) assembly pathways. These models predict the biophysical properties surrounding the central pore in BMC-H shell subunits, which in turn dictate the efficiency of substrate diffusion. Meanwhile, simulations on BMC assembly demonstrate that assembly pathway is largely dictated kinetically by cargo interactions while final morphology is dependent on shell factors. Overall, these findings are contextualized within the wider experimental BMC literature and framed within the opportunities for carboxysome redesign for biomanufacturing and enhanced carbon fixation.

59 BASIC BIOLOGICAL SCIENCES↗