Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
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 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↗
DOE OSTI2025

Optical Absorption Spectroscopy Probes Water Wire and Its Ordering in a Hydrogen-Bond Network

Water wires, quasi-one-dimensional chains composed of hydrogen-bonded (H-bonded) water molecules, play a fundamental role in numerous chemical, physical, and physiological processes. Yet direct experimental detection of water wires has been elusive so far. Based on advanced many-body theory that includes electron-hole interactions, we report that optical absorption spectroscopy can serve as a sensitive probe of water wires and their ordering. In both liquid and solid water, the main peak of the spectrum is discovered to be a charge-transfer exciton. In water, the charge-transfer exciton is strongly coupled to the H-bonding environment where the exciton is excited between H-bonded water molecules with a large spectral intensity. In regular ice, the spectral weight of the charge-transfer exciton is enhanced by a collective excitation occurring on proton-ordered water wires, whose spectral intensity scales with the ordering length of water wire. The spectral intensity and excitonic interaction strength reaches its maximum in ice XI, where the long-range ordering length yields the most pronounced spectral signal. Our findings suggest that water wires, which widely exist in important physiological and biological systems and other phases of ice, can be directly probed by this approach. Published by the American Physical Society 2025

Tang, Fujie (ORCID:0000000197615359)↗
DOE OSTI2025

Efficient analysis of small-angle scattering curves for large biomolecular assemblies using Monte Carlo methods

Structure elucidation from small-angle scattering curves of large biomolecular assemblies is notoriously challenging. This is because the simulation of high-resolution features in the structure of large macromolecular assemblies, such as de novo protein assemblies, is computationally demanding when it needs to cover a broad range of length scales. Conventional methods, such as the numerical approximation to the Debye equation or the use of spherical harmonics, do not scale well as the size of the assembly increases, which limits their application to small structures (e.g. individual proteins). This work explores the effectiveness of a Monte Carlo method to simulate and fit scattering curves for large biomolecular assemblies spanning over ranges covering atomic and molecular detail (e.g. spacing and orientation of proteins in an assembly) as well as large-scale (hundreds of nanometres) features. Owing to its speed and scalability, it can be combined with a fitting algorithm to extract structural features from experimental small-angle scattering curves in biomolecular assemblies that are otherwise intractable for interpretation. This work first demonstrates the effectiveness of the tool using experimental small-angle X-ray scattering (SAXS) data from tile-like proteins that assemble into 1D tube-like macromolecular structures. Here, the diameter distribution of tubes is extracted from SAXS fits, and this is quantitatively compared with distributions from electron microscopy. SAXS data are also obtained from 2D sheet-like protein assemblies, and the proposed method is used to quantify structural features such as the separation distance between protein building blocks and the flexing of the sheet. An open-source implementation of the methodology is provided for use in a broad range of biological systems involving multi-scale scattering analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗