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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Discovery, characterization, and application of chromosomal integration sites for stable heterologous gene expression in Rhodotorula toruloides

Rhodotorula toruloides is a non-model, oleaginous yeast uniquely suited to produce acetyl-CoA-derived chemicals. However, the lack of well-characterized genomic integration sites has impeded the metabolic engineering of this organism. Here we report a set of computationally predicted and experimentally validated chromosomal integration sites in R. toruloides. We first implemented an in silico platform by integrating essential gene information and transcriptomic data to identify candidate sites that meet stringent criteria. We then conducted a full experimental characterization of these sites, assessing integration efficiency, gene expression levels, impact on cell growth, and long-term expression stability. Among the identified sites, 12 exhibited integration efficiencies of 50% or higher, making them sufficient for most metabolic engineering applications. Using selected high-efficiency sites, we achieved simultaneous double and triple integrations and efficiently integrated long functional pathways (up to 14.7 kb). Additionally, we developed a new inducible marker recycling system that allows multiple rounds of integration at our characterized sites. Here, we validated this system by performing five sequential rounds of GFP integration and three sequential rounds of MaFAR integration for fatty alcohol production, demonstrating, for the first time, precise gene copy number tuning in R. toruloides. These characterized integration sites should significantly advance metabolic engineering efforts and future genetic tool development in R. toruloides.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Hydrogen-Bonded Organic Frameworks: A Rising Class of Porous Molecular Materials

Hydrogen-bonded organic frameworks (HOFs) are a class of porous molecular materials that rely on the assembly of organic building blocks by means of hydrogen-bonding interactions to form two-dimensional (2D) and three-dimensional (3D) crystalline networks. The reversible nature of the hydrogen-bond formation endows HOFs with the attributes of solution processability and simple regeneration. High-quality single crystals of HOFs can be grown easily for unambiguous superstructure determination by single-crystal X-ray diffraction, which is crucial for the elucidation of superstructure–property relationships. During the past decade, considerable progress has been achieved in realizing stable HOFs with permanent porosities by focusing on the design of molecular building blocks in order to introduce rigidity, auxiliary [π···π] interactions, and interpenetration of their frameworks to sustain the extended networks. The applications of HOFs are far-reaching, spanning catalysis, energy, and biomedical products as well as the storage and separation of fine chemicals. In this paper, we, first of all, provide an overview of the chronological development of HOFs, starting from the seminal work by Marsh and Duchamp in 1969 on the crystal superstructure of the hydrogen-bonded networks of trimesic acid. We identify the development of novel hydrogen-bonding motifs such as diaminotriazine (DTA), the introduction of the concept of molecular tectonics, and the establishment of permanent porosity in HOFs as being some of the milestones, which incentivized the current burgeoning research endeavors on developing HOFs as multifunctional materials. This Account is focused primarily on surveying the strategies for constructing porous 3D HOFs based on organic building blocks with peripheral carboxyl groups. These strategies are presented in the following categories: (1) the polycatenation of 2D networks by trigonal building blocks to form global 3D frameworks, (2) the utilization of building blocks with 3D geometries—tetrahedral and trigonal prismatic—that are predisposed to form 3D networks, and (3) the docking by shape-fitting of geometrically labile building blocks. We emphasize how the molecular geometry of the building blocks plays an important role in modulating the superstructures of extended frameworks so as to address specific applications. Recognizing that the in silico design of HOFs is the ultimate goal of researchers in this field, we also discuss the recent advances in superstructure prediction that lead to the formation of porous supramolecular crystals and assess the complications in implementing computational methods for HOFs with complex superstructures. We hope this Account will inspire the development of new supramolecular designs and creative approaches to crystal engineering that aid and abet the assembly of multifunctional HOFs with customizable properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Steric and Electronic Influence of Excited-State Decay in Cu(I) MLCT Chromophores

For the past eleven years, a dedicated effort in our research group focused on fundamentally advancing the photophysical properties of cuprous bis-phenanthroline-based metal-to-ligand charge transfer (MLCT) excited states. We rationalized that by gaining control over the numerous factors limiting the more widespread use of Cu I MLCT photosensitizers, they would be readily adopted in numerous light-activated applications given the earth-abundance of copper and the extensive library of 1,10-phenanthrolines developed over the last century. Significant progress has been achieved by recognizing valuable structure-property concepts developed by other researchers in tandem with detailed ultrafast and conventional time-scale investigations, in-silico-inspired molecular designs to predict spectroscopic properties, and applying novel synthetic methodologies. Ultimately, we achieved a plateau in exerting cooperative steric influence to control Cu I MLCT excited state decay. This led to combining sterics with π-conjugation and/or inductive electronic effects to further exert control over molecular photophysical properties. The lessons gleaned from our studies of homoleptic complexes were recently extended to heteroleptic bis(phenanthrolines) featuring enhanced visible light absorption properties and long-lived room-temperature photoluminescence. Furthermore, this Account navigates the reader through our intellectual journey of decision-making, molecular and experimental design, and data interpretation in parallel with appropriate background information related to the quantitative characterization of molecular photophysics using Cu I MLCT chromophores as prototypical examples.

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

Exploring the impacts of conformer selection methods on ion mobility collision cross section predictions

The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy threshold reductions, and similarity reduction. Generating 50,000 conformers each for 18 molecules, we used the In Silico Chemical Library Engine (ISiCLE) to calculate the collision cross sections for the entire dataset. First, we employed Monte Carlo simulations to understand the variability between conformer structures as generated using simulated annealing. Then we employed Monte Carlo simulations to the aforementioned conformer selection techniques applied on the simulated molecular property—the ion mobility collision cross section. Based on our analyses, we found Boltzmann weighting to be a good tradeoff between precision and theoretical accuracy. Combining multiple techniques revealed that energy thresholds and root-mean-squared deviation-based similarity reductions can save considerable computational expense while maintaining property prediction accuracy. Molecular dynamic conformer generation tools like AMBER can continue to generate new lowest energy conformers even after tens of thousands of generations, decreasing precision between runs. This reduced precision can be ameliorated and theoretical accuracy increased by running density functional theory geometry optimization on carefully selected conformers.

Nielson, Felicity F.↗
DOE OSTI2021

CFM-ID 4.0: More Accurate ESI MS/MS Spectral Prediction and Compound Identification

In the field of metabolomics, mass spectrometry (MS) is the method most commonly used for identifying and annotating metabolites. As this typically involves matching a given MS spectrum against an experimentally acquired reference spectral library, this approach is limited by the coverage and size of such libraries (which typically number in the thousands). These experimental libraries can be greatly extended by predicting the MS spectra of known chemical structures (which number in the millions) to create computational reference spectral libraries. To facilitate the generation of predicted spectral reference libraries we developed CFM-ID, a computer program that can accurately predict ESI-MS/MS spectrum for a given compound structure. CFM-ID is one of the best-performing methods for compound-to-mass-spectrum prediction, and also one of the top tools for in silico mass-spectrum-to-compound identification. This work improves CFM-ID’s ability to predict ESI-MS/MS spectra from compounds by: (1) learning parameters from features based on the molecular topology, (2) adding a new approach to ring cleavage that models such cleavage as a sequence of simple chemical bond dissociations and (3) expanding its hand-written rule-based predictor to cover more chemical classes, including acylcarnitines, acylcholines, flavonols, flavones, flavanones, and flavonoid glycosides. We demonstrate that this new version of CFM-ID (version 4.0) is significantly more accurate than previous CFM-ID versions, in terms of both EI-MS/MS spectral prediction and compound identification. CFM-ID 4.0 is available at http://cfmid4.wishartlab.com/ as a webservice and docker images can be downloaded at https://hub.docker.com/r/wishartlab/cfmid

Wang, Fei↗
DOE OSTI2024

Exploring the Effects of Intersubunit Interface Mutations on Virus-Like Particle Structure and Stability

Virus-like particles (VLPs) from bacteriophage MS2 provide a platform to study protein self-assembly and create engineered systems for drug delivery. Here, we aim to understand the impact of intersubunit interface mutations on the local and global structure and function of MS2-based VLPs. In previous work, our lab identified locally supercharged double mutants [T71K/G73R] that concentrate positive charge at capsid pores, enhancing uptake into mammalian cells. To study the effects of particle size on cellular internalization, we combined these double mutants with a single point mutation [S37P] that was previously reported to switch particle geometry from T = 3 to T = 1 icosahedral symmetry. These new variants retained their enhanced cellular uptake activity and could deliver small-molecule drugs with efficacy levels similar to our first-generation capsids. Surprisingly, these engineered triple mutants exhibit increased thermostability and unexpected geometry, producing T = 3 particles instead of the anticipated T = 1 assemblies. Transmission electron microscopy revealed various capsid assembly states, including wild-type (T = 3), T = 1, and rod-like particles, that could be accessed using different combinations of these point mutations. Molecular dynamics experiments recapitulated the structural rationale in silico for the single point mutation [S37P] forming a T = 1 virus-like particle and showed that this assembly state was not favored when combined with mutations that favor rod-like architectures. Through this work, we investigated how interdimer interface dynamics influence VLP size and morphology and how these properties affect particle function in applications such as drug delivery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Sequence Design of Random Heteropolymers as Protein Mimics

Random heteropolymers (RHPs) have been computationally designed and experimentally shown to recapitulate protein-like phase behavior and function. However, unlike proteins, RHP sequences are only statistically defined and cannot be sequenced. Recent developments in reversible-deactivation radical polymerization allowed simulated polymer sequences based on the well-established Mayo–Lewis equation to more accurately reflect ground-truth sequences that are experimentally synthesized. This led to opportunities to perform bioinformatics-inspired analysis on simulated sequences to guide the design, synthesis, and interpretation of RHPs. We compared batches on the order of 10000 simulated RHP sequences that vary by synthetically controllable and measurable RHP characteristics such as chemical heterogeneity and average degree of polymerization. Our analysis spans across 3 levels: segments along a single chain, sequences within a batch, and batch-averaged statistics. We discuss simulator fidelity and highlight the importance of robust segment definition. Examples are presented that demonstrate the use of simulated sequence analysis for in-silico iterative design to mimic protein hydrophobic/hydrophilic segment distributions in RHPs and compare RHP and protein sequence segments to explain experimental results of RHPs that mimic protein function. To facilitate the community use of this workflow, the simulator and analysis modules have been made available through an open source toolkit, the RHPapp.

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

Computational Prediction of Coiled–Coil Protein Gelation Dynamics and Structure

Protein hydrogels represent an important and growing biomaterial for a multitude of applications, including diagnostics and drug delivery. We have previously explored the ability to engineer the thermoresponsive supramolecular assembly of coiled–coil proteins into hydrogels with varying gelation properties, where we have defined important parameters in the coiled–coil hydrogel design. Using Rosetta energy scores and Poisson–Boltzmann electrostatic energies, we iterate a computational design strategy to predict the gelation of coiled–coil proteins while simultaneously exploring five new coiled–coil protein hydrogel sequences. Provided this library, we explore the impact of in silico energies on structure and gelation kinetics, where we also reveal a range of blue autofluorescence that enables hydrogel disassembly and recovery. As a result of this library, we identify the new coiled–coil hydrogel sequence, Q5, capable of gelation within 24 h at 4 °C, a more than 2-fold increase over that of our previous iteration Q2. The fast gelation time of Q5 enables the assessment of structural transition in real time using small-angle X-ray scattering (SAXS) that is correlated to coarse-grained and atomistic molecular dynamics simulations revealing the supramolecular assembling behavior of coiled–coils toward nanofiber assembly and gelation. This work represents the first system of hydrogels with predictable self-assembly, autofluorescent capability, and a molecular model of coiled–coil fiber formation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Cellular Automaton Simulation for Predicting Phase Evolution in Solid-State Reactions

New computational tools for solid-state synthesis recipe design are needed in order to accelerate the experimental realization of novel functional materials proposed by high-throughput materials discovery workflows. This work contributes a cellular automaton simulation framework for predicting the time-dependent evolution of intermediate and product phases during solid-state reactions as a function of precursor choice and amount, reaction atmosphere, and heating profile. The simulation captures the effects of reactant particle spatial distribution, particle melting, and reaction atmosphere. Reaction rates based on rudimentary kinetics are estimated using density functional theory data from the Materials Project and machine learning estimators for the melting point and the vibrational entropy component of the Gibbs free energy. The resulting simulation framework allows for the prediction of the likely outcome of a reaction recipe before any experiments are performed. We analyze five experimental solid-state recipes for BaTiO 3 , CaZrN 2 , and YMnO 3 found in the literature to illustrate the performance of the model in capturing reaction selectivity and reaction pathways as a function of temperature and precursor choice. This simulation framework offers an easier way to optimize existing recipes, aid in the identification of intermediates, and design effective recipes for yet unrealized inorganic solids in silico .

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

Machine Learning Models to Predict Inhibition of the Bile Salt Export Pump

Cholestatic liver injury is frequently associated with drug inhibition of bile salt transporters, such as the bile salt export pump (BSEP). Reliable in silico models to predict BSEP inhibition directly from chemical structures would significantly reduce costs during drug discovery and could help avoid injury to patients. We report our development of classification and regression models for BSEP inhibition with substantially improved performance over previously published models. We assessed the performance effects of different methods of chemical featurization, data set partitioning, and class labeling and identified the methods producing models that generalized best to novel chemical entities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Ligand- and Structure-Based Analysis of Deep Learning-Generated Potential α2a Adrenoceptor Agonists

The α2a adrenoceptor is a particularly medically relevant subtype of the G protein-couple receptor family. Unfortunately, high-throughput techniques aimed at producing novel drug leads targeted to this receptor have been largely unsuccessful due to the complex pharmacology of adrenergic receptors. As such, cutting-edge in silico ligand- and structure-based assessment and de novo deep learning methods are well-positioned to provide new insights into protein-ligand interactions and potential active compounds. In this work, we (i) collect a dataset of α2a adrenoceptor agonists and provide it as a resource for the drug design community; (ii) use the dataset as a basis to generate candidate active structures via deep learning; and (iii) apply computational ligand- and structure-based analysis techniques to gain new insight into α2a adrenoceptor agonists and assess the quality of the computer-generated compounds. Additionally, we further describe how such assessment techniques can be applied to putative chemical probes with a case study involved proposed medetomidine-based probes.

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

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

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

Protein–Ligand Binding Free-Energy Calculations with ARROW–A Purely First-Principles Parameterized Polarizable Force Field

Protein–ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)–a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein–ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ~0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

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

Molecular Modeling of Surfactant Interaction on Phospholipid Bilayers Mimicking Corneal Epithelium

Surfactants found in consumer products can compromise eye corneal membrane integrity upon accidental exposure. Traditional in vitro and in vivo approaches to evaluate membrane–surfactant interaction pose experimental limitations such as species variability, reproducibility, and most often do not provide the overall picture. These limitations motivate the use of in silico models to study phenomena like cellular disruption assays caused by surfactants at the molecular scale. In this work, coarse-grained molecular dynamics simulations have been employed to investigate how nonionic alcohol ethoxylate (AE) and anionic surfactant alcohol ethoxy sulfate (AES) interact with lipid bilayer liposomes that mimic corneal epithelial cell membranes. The spherical liposome is composed of 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-di(9Z-octadecenoyl)-sn-glycero-3-phosphoethanolamine (DOPE), 1,2-di(9Z-octadecenoyl)-sn-glycero-3-phospho-l-serine (DOPS), and cholesterol, resembling the composition of the corneal epithelial cells’ membrane bilayer. The simulation consisted of varying degrees of representative surfactant compositions and two initial types of surfactant configurations within or outside the liposome. Our results reveal that both surfactants induce outer leaflet bulging, agreeing with membrane solubilization models. The more highly ethoxylated surfactant, AE, caused more consistent inner leaflet disruption than AES, resulting in significantly more water permeation and membrane thinning. In addition, both surfactants increase the lateral diffusion of lipids within the membrane layers, with higher ethoxylated AE showing a stronger effect than AES. This study demonstrates how surfactant structure and localization influence bilayer membrane integrity, offering mechanistic insights into the irritation potential, thus guiding the rational design of effective surfactant-based formulations.

Lipids↗
DOE OSTI2021

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

In the explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

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

Univariate Prediction of Hammett Parameters and Select Relative Reaction Rates Using Loewdin Atomic Charges

Loewdin charges from density functional theory calculations were used here to obtain general, univariate linear correlations for the prediction of experimental Hammett parameters and relative reaction rates. While previous studies have established that Hirshfeld and CM5 charges perform strongly as univariate predictors, the near-ubiquitous Loewdin charges have not yet been evaluated. To this end, we assess the predictive capability of Loewdin charges for three chemical systems. First, we show that Loewdin charges outperform Hirshfeld and CM5 charges for Hammett parameter prediction. Second, we see that Loewdin charges generally perform comparably to Hirshfeld charges for predicting the relative rates of olefin cleavage by photoexcited nitroarenes. The single case of poor correlation, between relative rates and the Loewdin charges on nitrogen sites, is ameliorated when considering the net charge on the NO 2 group. Third, we show that Loewdin, Hirshfeld, and CM5 charges all perform very well for generating correlations for relative reaction rates for C–H activation of 9-(4-X-phenyl)-9H-fluorene substrates by a transition metal catalyst. The equations generated throughout the study enable the prediction of Hammett parameters and relative reaction rates. Finally, these tools can accelerate synthetic and experimental studies by enabling the in silico prediction of uncharacterized chemical properties.

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

Kinetic Separation of Siloxanes in Metal–Organic Frameworks

We present an in silico assessment of metal–organic frameworks (MOFs) for the kinetic separation of linear and cyclic siloxanes. We employed molecular dynamics simulations investigating both rigid and flexible 1D MOF frameworks to identify a specific range of pore parameters that enables the diffusion of linear siloxanes but leads to slow diffusion of cyclic siloxanes. We then extended our analysis to flexible 3D MOFs to select adsorbents for the kinetic separation of cyclic and linear siloxanes. Based on synthesizability metrics we identified four 3D MOFs capable of discriminating between cyclic and linear siloxanes. One of the MOFs with structure code WIYFAM stood out with the ability to distinguish between cyclic and linear siloxanes and facilitate the diffusion of all linear siloxanes investigated in this study. One of the other MOFs, IRMOF-6, is found to be capable of not only discriminating between cyclic and linear siloxanes, but even between shorter and longer linear siloxanes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗