SEARCH · Search NASA

Results for “Systems Biology”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2021

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the “gold standard” coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller−Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Chemistry↗

DOE OSTI2021

The F-box protein gene exo-1 is a target for reverse engineering enzyme hypersecretion in filamentous fungi

Carbohydrate active enzymes (CAZymes) are vital for the lignocellulose-based biorefinery. The development of hypersecreting fungal protein production hosts is therefore a major aim for both academia and industry. However, despite advances in our understanding of their regulation, the number of promising candidate genes for targeted strain engineering remains limited. Here, we resequenced the genome of the classical hypersecreting Neurospora crassa mutant exo-1 and identified the causative point of mutation to reside in the F-box protein–encoding gene, NCU09899. The corresponding deletion strain displayed amylase and invertase activities exceeding those of the carbon catabolite derepressed strain ?cre-1, while glucose repression was still mostly functional in ?exo-1. Surprisingly, RNA sequencing revealed that while plant cell wall degradation genes are broadly misexpressed in ?exo-1, only a small fraction of CAZyme genes and sugar transporters are up-regulated, indicating that EXO-1 affects specific regulatory factors. Aiming to elucidate the underlying mechanism of enzyme hypersecretion, we found the high secretion of amylases and invertase in ?exo-1 to be completely dependent on the transcriptional regulator COL-26. Furthermore, misregulation of COL-26, CRE-1, and cellular carbon and nitrogen metabolism was confirmed by proteomics. Finally, we successfully transferred the hypersecretion trait of the exo-1 disruption by reverse engineering into the industrially deployed fungus Myceliophthora thermophila using CRISPR-Cas9. Our identification of an important F-box protein demonstrates the strength of classical mutants combined with next-generation sequencing to uncover unanticipated candidates for engineering. These data contribute to a more complete understanding of CAZyme regulation and will facilitate targeted engineering of hypersecretion in further organisms of interest.

Gabriel, Raphael↗

DOE OSTI2022

An uncharacteristically low-potential flavin governs the energy landscape of electron bifurcation

Significance Nature has long been an inspiration for materials design, as it exemplifies exquisite control of both matter and energy. Electron bifurcation, a mechanism employed in biological systems to drive thermodynamically unfavorable and energetically challenging chemical reactions, is one such example. A key feature of bifurcating enzymes is the ability of a single redox cofactor to distribute a pair of electrons across two spatially separated electron transfer pathways. Here, we report on the empirical determination of both the one-electron potential and two-electron potential of the bifurcating flavin cofactor in the NADH-dependent ferredoxin-NADP + oxidoreductase I (NfnSL) enzyme. Insights arising from the defined energy landscape of this bifurcation site may underlie the design of synthetic catalysts capable of generating high-energy intermediates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Learning to learn by using nonequilibrium training protocols for adaptable materials

Evolution in time-varying environments naturally leads to adaptable biological systems that can easily switch functionalities. Advances in the synthesis of environmentally responsive materials therefore open up the possibility of creating a wide range of synthetic materials which can also be trained for adaptability. We consider high-dimensional inverse problems for materials where any particular functionality can be realized by numerous equivalent choices of design parameters. By periodically switching targets in a given design algorithm, we can teach a material to perform incompatible functionalities with minimal changes in design parameters. We exhibit this learning strategy for adaptability in two simulated settings: elastic networks that are designed to switch deformation modes with minimal bond changes and heteropolymers whose folding pathway selections are controlled by a minimal set of monomer affinities. The resulting designs can reveal physical principles, such as nucleation-controlled folding, that enable such adaptability.

Science & Technology - Other Topics↗

DOE OSTI2023

Engineering activatable promoters for scalable and multi-input CRISPRa/i circuits

Dynamic, multi-input gene regulatory networks (GRNs) are ubiquitous in nature. Multilayer CRISPR-based genetic circuits hold great promise for building GRNs akin to those found in naturally occurring biological systems. We develop an approach for creating high-performing activatable promoters that can be assembled into deep, wide, and multi-input CRISPR-activation and -interference (CRISPRa/i) GRNs. By integrating sequence-based design and in vivo screening, we engineer activatable promoters that achieve up to 1,000-fold dynamic range in an Escherichia coli-based cell-free system. These components enable CRISPRa GRNs that are six layers deep and four branches wide. We show the generalizability of the promoter engineering workflow by improving the dynamic range of the light-dependent EL222 optogenetic system from 6-fold to 34-fold. Additionally, high dynamic range promoters enable CRISPRa systems mediated by small molecules and protein–protein interactions. We apply these tools to build input-responsive CRISPRa/i GRNs, including feedback loops, logic gates, multilayer cascades, and dynamic pulse modulators. Our work provides a generalizable approach for the design of high dynamic range activatable promoters and enables classes of gene regulatory functions in cell-free systems.

42 ENGINEERING↗

DOE OSTI2023

Gas-phase microactuation using kinetically controlled surface states of ultrathin catalytic sheets

Biological systems convert chemical energy into mechanical work by using protein catalysts that assume kinetically controlled conformational states. Synthetic chemomechanical systems using chemical catalysis have been reported, but they are slow, require high temperatures to operate, or indirectly perform work by harnessing reaction products in liquids (e.g., heat or protons). Here, we introduce a bioinspired chemical strategy for gas-phase chemomechanical transduction that sequences the elementary steps of catalytic reactions on ultrathin (<10 nm) platinum sheets to generate surface stresses that directly drive microactuation (bending radii of 700 nm) at ambient conditions (T = 20 °C; P total = 1 atm). When fueled by hydrogen gas and either oxygen or ozone gas, we show how kinetically controlled surface states of the catalyst can be exploited to achieve fast actuation (600 ms/cycle) at 20 °C. Finally, we also show that the approach can integrate photochemically controlled reactions and can be used to drive the reconfiguration of microhinges and complex origami- and kirigami-based microstructures.

42 ENGINEERING↗

DOE OSTI2025

Engineering spin coherence in core-shell diamond nanocrystals

Fluorescent diamond nanocrystals can host spin qubit sensors capable of probing the physical properties of biological systems with nanoscale spatial resolution. Sub-100 nm diamond nanosensors can readily be delivered into intact cells and even living organisms. However, applications beyond current proof-of-principle experiments require a substantial increase in sensitivity, which is limited by surface induced charge instability and electron-spin dephasing. In this work, we utilize engineered core-shell structures to achieve a drastic increase in qubit coherence times (T 2 ) from 1.1 to 35 μs in bare nanodiamonds to upward of 52 to 87 μs. We use electron-paramagnetic-resonance results to present a band bending model and connect silica encapsulation to the removal of deleterious mid-gap surface states that are negatively affecting the qubit’s spin properties. Combined with a 1.9-fold increase in particle luminescence these advances correspond to up to two-order-of-magnitude reduction in integration time. Probing qubit dynamics at a single particle level further reveals that the noise characteristics fundamentally change from a bath with spins that rearrange their spatial configuration during the course of an experiment to a more dilute static bath. The observed results shed light on the underlying mechanisms governing fluorescence and spin properties in diamond nanocrystals and offer an effective noise mitigation strategy based on engineered core-shell structures.

core-shell↗

DOE OSTI2026

Three-dimensional high-content imaging of unstained soft tissue with subcellular resolution using a laboratory-based X-ray microscope

With increasing interest in studying biological systems across spatial scales—from centimeters down to nanometers—histology continues to be the gold standard for tissue imaging at cellular resolution, providing an essential bridge between macroscopic and nanoscopic analysis. However, its inherently destructive and two-dimensional nature limits its ability to capture the full three-dimensional complexity of tissue architecture. Here, we show that phase-contrast X-ray microscopy can enable three-dimensional virtual histology with subcellular resolution. This technique provides direct quantification of electron density without restrictive assumptions, allowing for direct characterization of cellular nuclei in a standard laboratory setting. By combining high spatial resolution and soft tissue contrast, with automated segmentation of cell nuclei, we demonstrated virtual Hematoxylin and Eosin (H&E) staining using machine learning-based style transfer, yielding volumetric datasets compatible with existing histopathological analysis tools. Furthermore, by integrating electron density and the sensitivity to nanometric features of the dark field contrast channel, we achieve stain-free, high-content imaging capable of distinguishing nuclei and extracellular matrix.

3D virtual histology↗

DOE OSTI2019

A biokinetic model for trivalent or hexavalent chromium in adult humans

Chromium exists in several oxidation states, with the trivalent state (Cr(III)) being the dominant naturally occurring form. Chromium in other oxidation states tends to be converted to the trivalent oxide in the natural environment and in biological systems. Chromium(III) has been shown to be an essential nutrient for humans and several non-human species. Chromium(VI), the second most stable form of chromium, is an important environmental contaminant that is mostly of industrial origin and is associated with lung cancer and nose tumours in chromium workers. This paper proposes a biokinetic model for chromium that addresses the distinctive behaviours of Cr(III) and Cr(VI) following uptake to blood of an adult human. Additionally, the model is based on biokinetic data derived from relatively short-term studies involving administration of chromium tracers to adult human subjects or laboratory animals, supplemented with data on the long-term distribution of chromium in adult humans as estimated from autopsy measurements. The model is part of a comprehensive update of biokinetic models of the International Commission on Radiological Protection, used to project or evaluate radiation doses from occupational intake of radionuclides.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2025

Poplar: a phylogenomics pipeline

Motivation Generating phylogenomic trees from the genomic data is essential in understanding biological systems. Each step of this complex process has received extensive attention and has been significantly streamlined over the years. Given the public availability of data, obtaining genomes for a wide selection of species is straightforward. However, analyzing that data to generate a phylogenomic tree is a multistep process with legitimate scientific and technical challenges, often requiring a significant input from a domain-area scientist. Results We present Poplar, a new, streamlined computational pipeline, to address the computational logistical issues that arise when constructing the phylogenomic trees. It provides a framework that runs state-of-the-art software for essential steps in the phylogenomic pipeline, beginning from a genome with or without an annotation, and resulting in a species tree. Running Poplar requires no external databases. In the execution, it enables parallelism for execution for clusters and cloud computing. The trees generated by Poplar match closely with state-of-the-art published trees. The usage and performance of Poplar is far simpler and quicker than manually running a phylogenomic pipeline. Availability and implementation Freely available on GitHub at https://github.com/sandialabs/poplar. Implemented using Python and supported on Linux.

Koning, Elizabeth [Sandia National Laboratories (S↗

DOE OSTI2024

Environment-specific virocell metabolic reprogramming

Abstract Viruses impact microbial systems through killing hosts, horizontal gene transfer, and altering cellular metabolism, consequently impacting nutrient cycles. A virus-infected cell, a “virocell,” is distinct from its uninfected sister cell as the virus commandeers cellular machinery to produce viruses rather than replicate cells. Problematically, virocell responses to the nutrient-limited conditions that abound in nature are poorly understood. Here we used a systems biology approach to investigate virocell metabolic reprogramming under nutrient limitation. Using transcriptomics, proteomics, lipidomics, and endo- and exo-metabolomics, we assessed how low phosphate (low-P) conditions impacted virocells of a marine Pseudoalteromonas host when independently infected by two unrelated phages (HP1 and HS2). With the combined stresses of infection and nutrient limitation, a set of nested responses were observed. First, low-P imposed common cellular responses on all cells (virocells and uninfected cells), including activating the canonical P-stress response, and decreasing transcription, translation, and extracellular organic matter consumption. Second, low-P imposed infection-specific responses (for both virocells), including enhancing nitrogen assimilation and fatty acid degradation, and decreasing extracellular lipid relative abundance. Third, low-P suggested virocell-specific strategies. Specifically, HS2-virocells regulated gene expression by increasing transcription and ribosomal protein production, whereas HP1-virocells accumulated host proteins, decreased extracellular peptide relative abundance, and invested in broader energy and resource acquisition. These results suggest that although environmental conditions shape metabolism in common ways regardless of infection, virocell-specific strategies exist to support viral replication during nutrient limitation, and a framework now exists for identifying metabolic strategies of nutrient-limited virocells in nature.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

eQuilibrator 3.0: a database solution for thermodynamic constant estimation

Abstract eQuilibrator (equilibrator.weizmann.ac.il) is a database of biochemical equilibrium constants and Gibbs free energies, originally designed as a web-based interface. While the website now counts around 1,000 distinct monthly users, its design could not accommodate larger compound databases and it lacked a scalable Application Programming Interface (API) for integration into other tools developed by the systems biology community. Here, we report on the recent updates to the database as well as the addition of a new Python-based interface to eQuilibrator that adds many new features such as a 100-fold larger compound database, the ability to add novel compounds, improvements in speed and memory use, and correction for Mg2+ ion concentrations. Moreover, the new interface can compute the covariance matrix of the uncertainty between estimates, for which we show the advantages and describe the application in metabolic modelling. We foresee that these improvements will make thermodynamic modelling more accessible and facilitate the integration of eQuilibrator into other software platforms.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Stereoselective recognition of morphine enantiomers by μ -opioid receptor

Stereospecific recognition of chiral molecules plays a crucial role in biological systems. The μ-opioid receptor (MOR) exhibits binding affinity towards (-)-morphine, a well-established gold standard in pain management, while it shows minimal binding affinity for the (+)-morphine enantiomer, resulting in a lack of analgesic activity. Understanding how MOR stereoselectively recognizes morphine enantiomers has remained a puzzle in neuroscience and pharmacology for over half-a-century due to the lack of direct observation techniques. To unravel this mystery, we constructed the binding and unbinding processes of morphine enantiomers with MOR via molecular dynamics simulations to investigate the thermodynamics and kinetics governing MOR's stereoselective recognition of morphine enantiomers. Our findings reveal that the binding of (-)-morphine stabilizes MOR in its activated state, exhibiting a deep energy well and a prolonged residence time. In contrast, (+)-morphine fails to sustain the activation state of MOR. Furthermore, the results suggest that specific residues, namely D114 2.50 and D147 3.32 , are deprotonated in the active state of MOR bound to (-)-morphine. This work highlights that the selectivity in molecular recognition goes beyond binding affinities, extending into the realm of residence time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Proteomic and metabolic disturbances in lignin-modified Brachypodium distachyon

Abstract Lignin biosynthesis begins with the deamination of phenylalanine and tyrosine (Tyr) as a key branch point between primary and secondary metabolism in land plants. Here, we used a systems biology approach to investigate the global metabolic responses to lignin pathway perturbations in the model grass Brachypodium distachyon. We identified the lignin biosynthetic protein families and found that ammonia-lyases (ALs) are among the most abundant proteins in lignifying tissues in grasses. Integrated metabolomic and proteomic data support a link between lignin biosynthesis and primary metabolism mediated by the ammonia released from ALs that is recycled for the synthesis of amino acids via glutamine. RNA interference knockdown of lignin genes confirmed that the route of the canonical pathway using shikimate ester intermediates is not essential for lignin formation in Brachypodium, and there is an alternative pathway from Tyr via sinapic acid for the synthesis of syringyl lignin involving yet uncharacterized enzymatic steps. Our findings support a model in which plant ALs play a central role in coordinating the allocation of carbon for lignin synthesis and the nitrogen available for plant growth. Collectively, these data also emphasize the value of integrative multiomic analyses to advance our understanding of plant metabolism.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Aryl hydrocarbon receptor-dependent toxicity by retene requires metabolic competence

Polycyclic aromatic hydrocarbons (PAHs) are a class of organic compounds frequently detected in the environment with widely varying toxicities. Many PAHs activate the aryl hydrocarbon receptor (AHR), inducing the expression of a battery of genes, including xenobiotic metabolizing enzymes like cytochrome P450s (CYPs); however, not all PAHs act via this mechanism. We screened several parent and substituted PAHs in in vitro AHR activation assays to classify their unique activity. Retene (1-methyl-7-isopropylphenanthrene) displays Ahr2-dependent teratogenicity in zebrafish, but did not activate human AHR or zebrafish Ahr2, suggesting a retene metabolite activates Ahr2 in zebrafish to induce developmental toxicity. To investigate the role of metabolism in retene toxicity, studies were performed to determine the functional role of cyp1a, cyp1b1, and the microbiome in retene toxicity, identify the zebrafish window of susceptibility, and measure retene uptake, loss, and metabolite formation in vivo. Cyp1a-null fish were generated using CRISPR-Cas9. Cyp1a-null fish showed increased sensitivity to retene toxicity, whereas Cyp1b1-null fish were less susceptible, and microbiome elimination had no significant effect. Zebrafish required exposure to retene between 24 and 48 hours post fertilization (hpf) to exhibit toxicity. After static exposure, retene concentrations in zebrafish embryos increased until 24 hpf, peaked between 24 and 36 hpf, and decreased rapidly thereafter. We detected retene metabolites at 36 and 48 hpf, indicating metabolic onset preceding toxicity. This study highlights the value of combining molecular and systems biology approaches with mechanistic and predictive toxicology to interrogate the role of biotransformation in AHR-dependent toxicity.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Bayesian differential programming for robust systems identification under uncertainty

This paper presents a machine learning framework for Bayesian systems identification from noisy, sparse and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of recent developments in differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling. This allows an efficient inference of the posterior distributions over plausible models with quantified uncertainty, while the use of sparsity-promoting priors enables the discovery of interpretable and parsimonious representations for the underlying latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed methods, including nonlinear oscillators, predator–prey systems and examples from systems biology. Taken together, our findings put forth a flexible and robust workflow for data-driven model discovery under uncertainty. All codes and data accompanying this article are available at https://bit.ly/34FOJMj .

Science & Technology - Other Topics↗

DOE OSTI2025

Capacitive response of biological membranes

We present a minimal model to analyze the capacitive response of a biological membrane subjected to a step voltage via blocking electrodes. Through a perturbative analysis of the underlying electrolyte transport equations, we show that the leading-order relaxation of the transmembrane potential is governed by a capacitive timescale, τ_{C}=λ_{D}L/D(2+Γδ^{M}/L/4+Γδ^{M}/λ_{D}), where λ_{D} is the Debye screening length, L is the electrolyte width, Γ is the ratio of the permittivity of the electrolyte to the membrane, δ^{M} is the membrane thickness, and D is the ionic diffusivity. This timescale is considerably shorter than the traditional RC timescale λ_{D}L/D for a bare electrolyte due to the membrane's low permittivity and finite thickness. Beyond the linear regime, however, salt diffusion in the bulk electrolyte drives a secondary, nonlinear relaxation process of the transmembrane potential over a longer timescale τ_{L}=L^{2}/4π^{2}D. A simple equivalent-circuit model accurately captures the linear behavior, and the perturbation expansion remains applicable across the entire range of observed physiological transmembrane potentials. Together, these findings underscore the importance of the faster capacitive timescale and nonlinear effects on the bulk diffusion timescale in determining transmembrane potential dynamics for a range of biological systems.

Farhadi, Jafar↗