Simulating the biophysical world, filling the gap between the micro and macro reality [Slides]
Abstract not provided.
59 BASIC BIOLOGICAL SCIENCES↗
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.
Columnar Localization and Laminar Origin of Cortical Surface Electrical Potentials
Electrocorticography (ECoG) methodologically bridges basic neuroscience and understanding of human brains in health and disease. However, the localization of ECoG signals across the surface of the brain and the spatial distribution of their generating neuronal sources are poorly understood. To address this gap, we recorded from rat auditory cortex using customized μECoG, and simulated cortical surface electrical potentials with a full-scale, biophysically detailed cortical column model. Experimentally, μECoG-derived auditory representations were tonotopically organized and signals were anisotropically localized to less than or equal to ±200 μm, that is, a single cortical column. Biophysical simulations reproduce experimental findings and indicate that neurons in cortical layers V and VI contribute ∼85% of evoked high-gamma signal recorded at the surface. Cell number and synchrony were the primary biophysical properties determining laminar contributions to evoked μECoG signals, whereas distance was only a minimal factor. Thus, evoked μECoG signals primarily originate from neurons in the infragranular layers of a single cortical column. SIGNIFICANCE STATEMENT ECoG methodologically bridges basic neuroscience and understanding of human brains in health and disease. However, the localization of ECoG signals across the surface of the brain and the spatial distribution of their generating neuronal sources are poorly understood. We investigated the localization and origins of sensory-evoked ECoG responses. We experimentally found that ECoG responses were anisotropically localized to a cortical column. Biophysically detailed simulations revealed that neurons in layers V and VI were the primary sources of evoked ECoG responses. These results indicate that evoked ECoG high-gamma responses are primarily generated by the population spike rate of pyramidal neurons in layers V and VI of single cortical columns and highlight the possibility of understanding how microscopic sources produce mesoscale signals.
Baratham, Vyassa L.↗
Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in situ and remote sensing data sets
The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).
Cheng, Yanyan↗
Ceramide-1-phosphate transfer protein enhances lipid transport by disrupting hydrophobic lipid–membrane contacts
Cellular distributions of the sphingolipid ceramide-1-phosphate (C1P) impact essential biological processes. C1P levels are spatiotemporally regulated by ceramide-1-phosphate transfer protein (CPTP), which efficiently shuttles C1P between organelle membranes. Yet, how CPTP rapidly extracts and inserts C1P into a membrane remains unknown. Here, we devise a multiscale simulation approach to elucidate biophysical details of CPTP-mediated C1P transport. We find that CPTP binds a membrane poised to extract and insert C1P and that membrane binding promotes conformational changes in CPTP that facilitate C1P uptake and release. By significantly disrupting a lipid’s local hydrophobic environment in the membrane, CPTP lowers the activation free energy barrier for passive C1P desorption and enhances C1P extraction from the membrane. Upon uptake of C1P, further conformational changes may aid membrane unbinding in a manner reminiscent of the electrostatic switching mechanism used by other lipid transfer proteins. Insertion of C1P into an acceptor membrane, eased by a decrease in membrane order by CPTP, restarts the transfer cycle. Most notably, we provide molecular evidence for CPTP’s ability to catalyze C1P extraction by breaking hydrophobic C1P–membrane contacts with compensatory hydrophobic lipid–protein contacts. Our work, thus, provides biophysical insights into how CPTP efficiently traffics C1P between membranes to maintain sphingolipid homeostasis and, additionally, presents a simulation method aptly suited for uncovering the catalytic mechanisms of other lipid transfer proteins.
59 BASIC BIOLOGICAL SCIENCES↗
Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in-situ and remote sensing datasets
The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).
54 ENVIRONMENTAL SCIENCES↗
Identification of preferred multimodal ligand‐binding regions on IgG1 F C using nuclear magnetic resonance and molecular dynamics simulations
Abstract In this study, the binding of multimodal chromatographic ligands to the IgG1 F C domain were studied using nuclear magnetic resonance and molecular dynamics simulations. Nuclear magnetic resonance experiments carried out with chromatographic ligands and a perdeuterated 15 N‐labeled F C domain indicated that while single‐mode ion exchange ligands interacted very weakly throughout the F C surface, multimodal ligands containing negatively charged and aromatic moieties interacted with specific clusters of residues with relatively high affinity, forming distinct binding regions on the F C . The multimodal ligand‐binding sites on the F C were concentrated in the hinge region and near the interface of the C H 2 and C H 3 domains. Furthermore, the multimodal binding sites were primarily composed of positively charged, polar, and aliphatic residues in these regions, with histidine residues exhibiting some of the strongest binding affinities with the multimodal ligand. Interestingly, comparison of protein surface property data with ligand interaction sites indicated that the patch analysis on F C corroborated molecular‐level binding information obtained from the nuclear magnetic resonance experiments. Finally, molecular dynamics simulation results were shown to be qualitatively consistent with the nuclear magnetic resonance results and to provide further insights into the binding mechanisms. An important contribution to multimodal ligand‐F C binding in these preferred regions was shown to be electrostatic interactions and π–π stacking of surface‐exposed histidines with the ligands. This combined biophysical and simulation approach has provided a deeper molecular‐level understanding of multimodal ligand–F C interactions and sets the stage for future analyses of even more complex biotherapeutics.
Gudhka, Ronak B.↗
Catalyst: Fast and flexible modeling of reaction networks
We introduce Catalyst.jl, a flexible and feature-filled Julia library for modeling and high-performance simulation of chemical reaction networks (CRNs). Catalyst supports simulating stochastic chemical kinetics (jump process), chemical Langevin equation (stochastic differential equation), and reaction rate equation (ordinary differential equation) representations for CRNs. Through comprehensive benchmarks, we demonstrate that Catalyst simulation runtimes are often one to two orders of magnitude faster than other popular tools. More broadly, Catalyst acts as both a domain-specific language and an intermediate representation for symbolically encoding CRN models as Julia-native objects. This enables a pipeline of symbolically specifying, analyzing, and modifying CRNs; converting Catalyst models to symbolic representations of concrete mathematical models; and generating compiled code for numerical solvers. Leveraging ModelingToolkit.jl and Symbolics.jl, Catalyst models can be analyzed, simplified, and compiled into optimized representations for use in numerical solvers. Finally, we demonstrate Catalyst’s broad extensibility and composability by highlighting how it can compose with a variety of Julia libraries, and how existing open-source biological modeling projects have extended its intermediate representation.
59 BASIC BIOLOGICAL SCIENCES↗
Uncertainties in estimating global potential yields and their impacts for long-term modeling
Estimating realistic potential yields by crop type and region is challenging; such yields depend on both biophysical characteristics (e.g., soil characteristics, climate, etc.), and the crop management practices available in any site or region (e.g., mechanization, irrigation, crop cultivars). A broad body of literature has assessed potential yields for selected crops and regions, using several strategies. In this study we first analyze future potential yields of major crop types globally by two different estimation methods, one of which is based on historical observed yields (“Empirical”), while the other is based on biophysical conditions (“Simulated”). Potential yields by major crop and region are quite different between the two methods; in particular, Simulated potential yields are typically 200% higher than Empirical potential yields in tropical regions for major crops. Applying both of these potential yields in yield gap closure scenarios in a global agro-economic model, GCAM, the two estimates of future potential yields lead to very different outcomes for the agricultural sector globally. In the Simulated potential yield closure scenario, Africa, Asia, and South America see comparatively favorable outcomes for agricultural sustainability over time: low land use change emissions, low crop prices, and high levels of self-sufficiency. In contrast, the Empirical potential yield scenario is characterized by a heavy reliance on production and exports in temperate regions that currently practice industrial agriculture. At the global level, this scenario has comparatively high crop commodity prices, and more land allocated to crop production (and associated land use change emissions) than either the baseline or Simulated potential yield scenarios. This study highlights the importance of the choice of methods of estimating potential yields for agro-economic modeling.
60 APPLIED LIFE SCIENCES↗
Global Life Cycle and Techno-Economic Assessment of Algal-Based Biofuels
Techno-economic analyses (TEAs) and life cycle assessments (LCAs) of algal biofuels often focus on locations in suboptimal latitudes for algal cultivation, which can under-represent the sustainability potential of the technology. This study identifies the optimal global productivity potential, environmental impacts, and economic viability of algal biofuels by using validated biophysical and sustainability modeling. The biophysical model simulates growth rates of Scenedesmus obliquusbased on temperature, photoinhibition, and respiration effects at 6685 global locations. Region-specific labor costs, construction factors, and tax rates allow for spatially resolved TEA, while the LCA includes regional impacts of electricity, hydrogen, and nutrient markets across ten environmental categories. The analysis identifies optimal locations for algal biofuel production in terms of environmental impacts and economic viability which are shown to follow biomass yields. Modeling results highlight the global variability of productivity with maximum yields ranging between 24.8 and 27.5 g m -2 d -1 in equatorial regions. Environmental impact results show favorable locations tracked with low-carbon electricity grids, with the well-to-wheels global warming potential (GWP) ranging from 31 to 45 g CO 2eq MJ -1 in South America and Central Africa. When including direct land use change impacts, the GWP ranged between 44 and 55 g CO 2eq MJ -1 in these high-productivity regions. Low-carbon electricity also favors air quality and eutrophication impacts. The TEA shows that minimum algal fuel prices of $\$1.89$-$\$2.15$ per liter of gasoline-equivalent are possible in southeast Asia and Venezuela. Furthermore, this discussion focuses on the challenges and opportunities to reduce fuel prices and the environmental impacts of algal biofuels in various global regions.
09 BIOMASS FUELS↗
Biophysical principles of choanoflagellate self-organization
Inspired by the patterns of multicellularity in choanoflagellates, the closest living relatives of animals, we quantify the biophysical processes underlying the morphogenesis of rosette colonies in the choanoflagellate Salpingoeca rosetta . We find that rosettes reproducibly transition from an early stage of 2-dimensional (2D) growth to a later stage of 3D growth, despite the underlying variability of the cell lineages. Our perturbative experiments demonstrate the fundamental importance of a basally secreted extracellular matrix (ECM) for rosette morphogenesis and show that the interaction of the ECM with cells in the colony physically constrains the packing of proliferating cells and, thus, controls colony shape. Simulations of a biophysically inspired model that accounts for the size and shape of the individual cells, the fraction of ECM, and its stiffness relative to that of the cells suffices to explain our observations and yields a morphospace consistent with observations across a range of multicellular choanoflagellate colonies. Overall, our biophysical perspective on rosette development complements previous genetic perspectives and, thus, helps illuminate the interplay between cell biology and physics in regulating morphogenesis.
59 BASIC BIOLOGICAL SCIENCES↗
OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations
We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.
59 BASIC BIOLOGICAL SCIENCES↗
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↗
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↗
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↗
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↗