Search NASA⌕ Search

SEARCH · Search NASA

Results for “in silico”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

HypoRiPPAtlas as an Atlas of hypothetical natural products for mass spectrometry database search

Abstract Recent analyses of public microbial genomes have found over a million biosynthetic gene clusters, the natural products of the majority of which remain unknown. Additionally, GNPS harbors billions of mass spectra of natural products without known structures and biosynthetic genes. We bridge the gap between large-scale genome mining and mass spectral datasets for natural product discovery by developing HypoRiPPAtlas, an Atlas of hypothetical natural product structures, which is ready-to-use for in silico database search of tandem mass spectra. HypoRiPPAtlas is constructed by mining genomes using seq2ripp, a machine-learning tool for the prediction of ribosomally synthesized and post-translationally modified peptides (RiPPs). In HypoRiPPAtlas, we identify RiPPs in microbes and plants. HypoRiPPAtlas could be extended to other natural product classes in the future by implementing corresponding biosynthetic logic. This study paves the way for large-scale explorations of biosynthetic pathways and chemical structures of microbial and plant RiPP classes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Deciphering metabolic differentiation during Bacillus subtilis sporulation

Abstract The bacteriumBacillus subtilisundergoes asymmetric cell division during sporulation, producing a mother cell and a smaller forespore connected by the SpoIIQ-SpoIIIA (or Q-A) channel. The two cells differentiate metabolically, and the forespore becomes dependent on the mother cell for essential building blocks. Here, we investigate the metabolic interactions between mother cell and forespore using genome-scale metabolic and expression models as well as experiments. Our results indicate that nucleotides are synthesized in the mother cell and transported in the form of nucleoside di- or tri-phosphates to the forespore via the Q-A channel. However, if the Q-A channel is inactivated later in sporulation, then glycolytic enzymes can form an ATP and NADH shuttle, providing the forespore with energy and reducing power. Our integrated in silico and in vivo approach sheds light into the intricate metabolic interactions underlying cell differentiation inB. subtilis, and provides a foundation for future studies of metabolic differentiation.

Science & Technology - Other Topics↗
DOE OSTI2025

Functional protein mining with conformal guarantees

Molecular structure prediction and homology detection offer promising paths to discovering protein function and evolutionary relationships. However, current approaches lack statistical reliability assurances, limiting their practical utility for selecting proteins for further experimental and in-silico characterization. To address this challenge, we introduce a statistically principled approach to protein search leveraging principles from conformal prediction, offering a framework that ensures statistical guarantees with user-specified risk and provides calibrated probabilities (rather than raw ML scores) for any protein search model. Our method (1) lets users select many biologically-relevant loss metrics (i.e. false discovery rate) and assigns reliable functional probabilities for annotating genes of unknown function; (2) achieves state-of-the-art performance in enzyme classification without training new models; and (3) robustly and rapidly pre-filters proteins for computationally intensive structural alignment algorithms. Our framework enhances the reliability of protein homology detection and enables the discovery of uncharacterized proteins with likely desirable functional properties.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Design of diverse, functional mitochondrial targeting sequences across eukaryotic organisms using variational autoencoder

Mitochondria play a key role in energy production and metabolism, making them a promising target for metabolic engineering and disease treatment. However, despite the known influence of passenger proteins on localization efficiency, only a few protein-localization tags have been characterized for mitochondrial targeting. To address this limitation, we leverage a Variational Autoencoder to design novel mitochondrial targeting sequences. In silico analysis reveals that a high fraction of the generated peptides (90.14%) are functional and possess features important for mitochondrial targeting. We characterize artificial peptides in four eukaryotic organisms and, as a proof-of-concept, demonstrate their utility in increasing 3-hydroxypropionic acid titers through pathway compartmentalization and improving 5-aminolevulinate synthase delivery by 1.62-fold and 4.76-fold, respectively. Moreover, we employ latent space interpolation to shed light on the evolutionary origins of dual-targeting sequences. Overall, our work demonstrates the potential of generative artificial intelligence for both fundamental research and practical applications in mitochondrial biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Enhancing tandem mass spectrometry-based metabolite annotation with online chemical labeling

Abstract Metabolite identification in non-targeted mass spectrometry-based metabolomics remains a major challenge due to limited spectral library coverage and difficulties in predicting metabolite fragmentation patterns. Here, we introduce Multiplexed Chemical Metabolomics (MCheM), which employs orthogonal post-column derivatization reactions integrated into a unified mass spectrometry data framework. MCheM generates orthogonal structural information that substantially improves metabolite annotation through in silico spectrum matching and open-modification searches, offering a powerful new toolbox for the structure elucidation of unknown metabolites at scale.

Science & Technology - Other Topics↗
DOE OSTI2025

Data-driven organic solubility prediction at the limit of aleatoric uncertainty

Abstract Small molecule solubility is a critically important property which affects the efficiency, environmental impact, and phase behavior of synthetic processes. Experimental determination of solubility is a time- and resource-intensive process and existing methods for in silico estimation of solubility are limited by their generality, speed, and accuracy. This work presents two models derived from the FASTPROP and CHEMPROP architectures and trained on BigSolDB which are capable of predicting solubility at arbitrary temperatures for a wide range of small molecules in organic solvent. Both extrapolate to unseen solutes 2–3 times more accurately than the current state-of-the-art model and we demonstrate that they are approaching the aleatoric limit (0.5–1$$\log S$$ log S ) of available test data, suggesting that further improvements in prediction accuracy require more accurate datasets. The FASTPROP-derived model (called FASTSOLV) and the CHEMPROP-based model are open source, freely accessible via a Python package and web interface, highly reproducible, and up to 2 orders of magnitude faster than current alternatives.

Science & Technology - Other Topics↗
DOE OSTI2024

Metabolic interactions shape emergent biofilm structures in a conceptual model of gut mucosal bacterial communities

Abstract The gut microbiome plays a major role in human health; however, little is known about the structural arrangement of microbes and factors governing their distribution. In this work, we present an in silico agent-based model (ABM) to conceptually simulate the dynamics of gut mucosal bacterial communities. We explored how various types of metabolic interactions, including competition, neutralism, commensalism, and mutualism, affect community structure, through nutrient consumption and metabolite exchange. Results showed that, across scenarios with different initial species abundances, cross-feeding promotes species coexistence. Morphologically, competition and neutralism resulted in segregation, while mutualism and commensalism fostered high intermixing. In addition, cooperative relations resulted in community properties with little sensitivity to the selective uptake of metabolites produced by the host. Moreover, metabolic interactions strongly influenced colonization success following the invasion of newcomer species. These results provide important insights into the utility of ABM in deciphering complex microbiome patterns.

Biotechnology & Applied Microbiology↗
DOE OSTI2025

Deep-learning atomistic semi-empirical pseudopotential model for nanomaterials

The semi-empirical pseudopotential method (SEPM) has been widely applied to provide computational insights into the electronic structure, photophysics, and charge carrier dynamics of nanoscale materials. We present “DeepPseudopot”, a machine-learned atomistic pseudopotential model that extends the SEPM framework by combining a flexible neural network representation of the local pseudopotential with parameterized non-local and spin-orbit coupling terms. Trained on bulk quasiparticle band structures and deformation potentials from GW calculations, the model captures many-body and relativistic effects with very high accuracy across diverse semiconducting materials, as illustrated for silicon and group III-V semiconductors. DeepPseudopot’s accuracy, efficiency, and transferability make it well-suited for data-driven in silico design and discovery of novel optoelectronic nanomaterials.

Lin, Kailai [University of California, Berkeley, C↗
DOE OSTI2024

Exploring the frontiers of condensed-phase chemistry with a general reactive machine learning potential

Abstract Atomistic simulation has a broad range of applications from drug design to materials discovery. Machine learning interatomic potentials (MLIPs) have become an efficient alternative to computationally expensive ab initio simulations. For this reason, chemistry and materials science would greatly benefit from a general reactive MLIP, that is, an MLIP that is applicable to a broad range of reactive chemistry without the need for refitting. Here we develop a general reactive MLIP (ANI-1xnr) through automated sampling of condensed-phase reactions. ANI-1xnr is then applied to study five distinct systems: carbon solid-phase nucleation, graphene ring formation from acetylene, biofuel additives, combustion of methane and the spontaneous formation of glycine from early earth small molecules. In all studies, ANI-1xnr closely matches experiment (when available) and/or previous studies using traditional model chemistry methods. As such, ANI-1xnr proves to be a highly general reactive MLIP for C, H, N and O elements in the condensed phase, enabling high-throughput in silico reactive chemistry experimentation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The persistence of memory in ionic conduction probed by nonlinear optics

Predicting practical rates of transport in condensed phases enables the rational design of materials, devices and processes. This is especially critical to developing low-carbon energy technologies such as rechargeable batteries. For ionic conduction, the collective mechanisms, variation of conductivity with timescales and confinement, and ambiguity in the phononic origin of translation, call for a direct probe of the fundamental steps of ionic diffusion: ion hops. However, such hops are rare-event large-amplitude translations, and are challenging to excite and detect. Here we use single-cycle terahertz pumps to impulsively trigger ionic hopping in battery solid electrolytes. This is visualized by an induced transient birefringence, enabling direct probing of anisotropy in ionic hopping on the picosecond timescale. The relaxation of the transient signal measures the decay of orientational memory, and the production of entropy in diffusion. We extend experimental results using in silico transient birefringence to identify vibrational attempt frequencies for ion hopping. Using nonlinear optical methods, we probe ion transport at its fastest limit, distinguish correlated conduction mechanisms from a true random walk at the atomic scale, and demonstrate the connection between activated transport and the thermodynamics of information.

25 ENERGY STORAGE↗
DOE OSTI2026

Trafficking of a nitrogenase FeMo-cofactor assembly intermediate

The maturation of the unique FeMo-cofactor of molybdenum nitrogenase is a multistep process requiring the sequential action of a series of maturase complexes. As a final step, the NifEN complex forms FeMo-cofactor from the precursor NifB-co, also called L-cluster, through replacement of an apical iron ion by molybdenum and the attachment of an organic homocitrate ligand. NifB-co is delivered by a small cofactor chaperone, NifX, and initially bound near the surface of the maturase NifEN. Here, we report high-resolution cryo-electron microscopy structures of NifEN in complex with NifX, showing NifB-co binding to NifEN in full detail, capturing both interacting partners in the act of cluster transfer. In a dynamic transfer complex, the large metal cluster is coordinated by single residues from both NifEN and NifX. In silico studies concur with these structures but suggest a third, internal conversion site where cluster maturation likely takes place.

metalloproteins↗
DOE OSTI2020

Database of pharmacokinetic time-series data and parameters for 144 environmental chemicals

Time courses of compound concentrations in plasma are used in chemical safety analysis to evaluate the relationship between external administered doses and internal tissue exposures. This type of experimental data is rarely available for the thousands of non-pharmaceutical chemicals to which people may potentially be unknowingly exposed but is necessary to properly assess the risk of such exposures. In vitro assays and in silico models are often used to craft an understanding of a chemical’s pharmacokinetics; however, the certainty of the quantitative application of these estimates for chemical safety evaluations cannot be determined without in vivo data for external validation. To address this need, we present a public database of chemical time-series concentration data from 567 studies in humans or test animals for 144 environmentally-relevant chemicals and their metabolites (187 analytes total). All major administration routes are incorporated, with concentrations measured in blood/plasma, tissues, and excreta. We also include calculated pharmacokinetic parameters for some studies, and a bibliography of additional source documents to support future extraction of time-series. In addition to pharmacokinetic model calibration and validation, these data may be used for analyses of diferential chemical distribution across chemicals, species, doses, or routes, and for meta-analyses on pharmacokinetic studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules

We introduce QM7-X, a comprehensive dataset of 42 physicochemical properties for ≈4.2 million equilibrium and non-equilibrium structures of small organic molecules with up to seven non-hydrogen (C, N, O, S, Cl) atoms. To span this fundamentally important region of chemical compound space (CCS), QM7-X includes an exhaustive sampling of (meta-)stable equilibrium structures—comprised of constitutional/structural isomers and stereoisomers, e.g., enantiomers and diastereomers (including cis-/trans- and conformational isomers)—as well as 100 non-equilibrium structural variations thereof to reach a total of ≈4.2 million molecular structures. Computed at the tightly converged quantum-mechanical PBE0+MBD level of theory, QM7-X contains global (molecular) and local (atom-in-a-molecule) properties ranging from ground state quantities (such as atomization energies and dipole moments) to response quantities (such as polarizability tensors and dispersion coefficients). By providing a systematic, extensive, and tightly-converged dataset of quantum-mechanically computed physicochemical properties, we expect that QM7-X will play a critical role in the development of next-generation machine-learning based models for exploring greater swaths of CCS and performing in silico design of molecules with targeted properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Structural characterization of a GNAT family acetyltransferase from Elizabethkingia anophelis bound to acetyl-CoA reveals a new dimeric interface

General control non-repressible 5 (GCN5)-related N-acetyltransferases (GNATs) catalyse the acetylation of a diverse range of substrates, thereby orchestrating a variety of biological processes within prokaryotes and eukaryotes. GNAT enzymes can catalyze the transfer of an acetyl group from acetyl coenzyme A to substrates such as aminoglycoside antibiotics, amino acids, polyamines, peptides, vitamins, catecholamines, and large macromolecules including proteins. Although GNATs generally exhibit low to moderate sequence identity, they share a conserved catalytic fold and conserved structural motifs. In this current study we characterize the high-resolution X-ray crystallographic structure of a GNAT enzyme bound with acetyl-CoA from Elizabethkingia anophelis , an important multi-drug resistant bacterium. The tertiary structure is comprised of six α-helices and nine β-strands, and is similar with other GNATs. We identify a new and uncharacterized GNAT dimer interface, which is conserved in at least two other unpublished GNAT structures. This suggests that GNAT enzymes can form at least five different types of dimers, in addition to a range of other oligomers including trimer, tetramer, hexamer, and dodecamer assemblies. The high-resolution structure presented in this study is suitable for future in-silico docking and structure–activity relationship studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A data-driven approach to modeling cancer cell mechanics during microcirculatory transport

In order to understand the effect of cellular level features on the transport of circulating cancer cells in the microcirculation, there has been an increasing reliance on high-resolution in silico models. Accurate simulation of cancer cells flowing with blood cells requires resolving cellular-scale interactions in 3D, which is a significant computational undertaking warranting a cancer cell model that is both computationally efficient yet sufficiently complex to capture relevant behavior. Given that the characteristics of metastatic spread are known to depend on cancer type, it is crucial to account for mechanistic behavior representative of a specific cancer’s cells. To address this gap, in the present work we develop and validate a means by which an efficient and popular membrane model-based approach can be used to simulate deformable cancer cells and reproduce experimental data from specific cell lines. Here, cells are modeled using the immersed boundary method (IBM) within a lattice Boltzmann method (LBM) fluid solver, and the finite element method (FEM) is used to model cell membrane resistance to deformation. Through detailed comparisons with experiments, we (i) validate this model to represent cancer cells undergoing large deformation, (ii) outline a systematic approach to parameterize different cell lines to optimally fit experimental data over a range of deformations, and (iii) provide new insight into nucleated vs. non-nucleated cell models and their ability to match experiments. While many works have used the membrane-model based method employed here to model generic cancer cells, no quantitative comparisons with experiments exist in the literature for specific cell lines undergoing large deformation. Here, we describe a phenomenological, data-driven approach that can not only yield good agreement for large deformations, but explicitly detail how it can be used to represent different cancer cell lines. This model is readily incorporated into cell-resolved hemodynamic transport simulations, and thus offers significant potential to complement experiments towards providing new insights into various aspects of cancer progression.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Development of a CNS-permeable reactivator for nerve agent exposure: an iterative, multi-disciplinary approach

Nerve agents have experienced a resurgence in recent times with their use against civilian targets during the attacks in Syria (2012), the poisoning of Sergei and Yulia Skripal in the United Kingdom (2018) and Alexei Navalny in Russia (2020), strongly renewing the importance of antidote development against these lethal substances. The current standard treatment against their effects relies on the use of small molecule-based oximes that can efficiently restore acetylcholinesterase (AChE) activity. Despite their efficacy in reactivating AChE, the action of drugs like 2-pralidoxime (2-PAM) is primarily limited to the peripheral nervous system (PNS) and, thus, provides no significant protection to the central nervous system (CNS). This lack of action in the CNS stems from their ionic nature that, on one end makes them very powerful reactivators and on the other renders them ineffective at crossing the Blood Brain Barrier (BBB) to reach the CNS. In this report, we describe the use of an iterative approach composed of parallel chemical and in silico syntheses, computational modeling, and a battery of detailed in vitro and in vivo assays that resulted in the identification of a promising, novel CNS-permeable oxime reactivator. Additional experiments to determine acute and chronic toxicity are ongoing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Crystal structure of ultra-humanized anti-pTau Fab reveals how germline substitutions humanize CDRs without loss of binding’

Administration of therapeutic antibodies can elicit adverse immune responses in patients through the generation of anti-drug antibodies that, in turn, reduce the efficacy of the therapeutic. Removal of foreign amino acid content by humanization can lower the immunogenic risk of the therapeutic mAb. We previously developed the ultra-humanization technology “Augmented Binary Substitution” (ABS) which enables single-step CDR germlining of antibodies. The application of ABS to a chicken anti-pTau antibody generated an ultra-humanized variant, anti-pTau C21-ABS, with increased human amino acid content in the CDRs and reduced in-silico predicted immunogenicity risk. Here, we report the high-resolution crystal structure of anti-pTau C21-ABS Fab in complex with the pTau peptide (7KQK). This study examines how ultra-humanization, via CDR germlining, is facilitated while maintaining near-identical antigen affinity (within 1.6-fold). The co-complex structure reveals that the ABS molecule targets the same antigenic epitope, accommodated by structurally-similar changes in the paratope. These findings confirm that ABS enables the germlining of amino acids within CDRs by exploiting CDR plasticity, to reduce non-human amino acid CDR content, with few alterations to the overall mechanism of binding.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Development and validation of a mathematical model of heart rate response to fluid perturbation

Physiological closed-loop controlled (PCLC) medical devices monitor and automatically adjust the patient’s condition by using physiological variables as feedback, ideally with minimal human intervention to achieve the target levels set by a clinician. PCLC devices present a challenge when it comes to evaluating their performance, where conducting large clinical trials can be expensive. Virtual physiological patients simulated by validated mathematical models can be utilized to obtain preclinical evidence of safety and assess the performance of the PCLC medical device during normal and worst-case conditions that are unlikely to happen in a limited clinical trial. A physiological variable that plays a major role during fluid resuscitation is heart rate (HR). For in silico assessment of PCLC medical devices regarding fluid perturbation, there is currently no mathematical model of HR validated in terms of its predictive capability performance. This paper develops and validates a mathematical model of HR response using data collected from sheep subjects undergoing hemorrhage and fluid infusion. The model proved to be accurate in estimating the HR response to fluid perturbation, where averaged between 21 calibration datasets, the fitting performance showed a normalized root mean square error (NRMSE) of 7.41 ± 2.8%. The model was also evaluated in terms of model predictive capability performance via a leave-one-out procedure (21 subjects) and an independent validation dataset (6 subjects). Two different virtual cohort generation tools were used in each validation analysis. The generated envelope of virtual subjects robustly met the defined acceptance criteria, in which 95% of the testing datasets presented simulated HR patterns that were within a deviation of 50% from the observed data. In addition, out of 16000 and 18522 simulated subjects for the leave one-out and independent datasets, the model was able to generate at least one virtual subject that was close to the real subject within an error margin of 9.56 ± 3.15% and 11.1 ± 1.22% NRMSE, respectively. In conclusion, the model can generate valid virtual HR physiological responses to fluid perturbation and be incorporated into future non-clinical simulated testing setups for assessing PCLC devices intended for fluid resuscitation.

59 BASIC BIOLOGICAL SCIENCES↗