Search NASA⌕ Search

SEARCH · Search NASA

Results for “in silico design”

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 19 records

In Silico Design of Methyl-Driven Overhauser Dynamic Nuclear Polarization Agents

Overhauser effect (OE) dynamic nuclear polarization (DNP) has drawn attention owing to its enhanced performance at ultrahigh magnetic fields. The lack of design principles has, nevertheless, limited the development of OE polarizing agents when compared to those used for cross-effect DNP. Here, we measured the 19 F OE DNP performance of a series of CF 3 -functionalized Blatter-type radicals. Using density functional theory calculations, we accurately predict the methyl-driven OE performance, paving the way for computer-aided design of optimized OE DNP polarizing agents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

“Freedom of design” in chemical compound space: towards rational in silico design of molecules with targeted quantum-mechanical properties

The rational design of molecules with targeted quantum-mechanical (QM) properties requires an advanced understanding of the structure–property/property–property relationships (SPR/PPR) that exist across chemical compound space (CCS). In this work, we analyze these fundamental relationships in the sector of CCS spanned by small (primarily organic) molecules using the recently developed QM7-X dataset, a systematic, extensive, and tightly converged collection of 42 QM properties corresponding to ≈4.2M equilibrium and non-equilibrium molecular structures containing up to seven heavy/non-hydrogen atoms (including C, N, O, S, and Cl). By characterizing and enumerating progressively more complex manifolds of molecular property space—the corresponding high-dimensional space defined by the properties of each molecule in this sector of CCS—our analysis reveals that one has a substantial degree of flexibility or “freedom of design” when searching for a single molecule with a desired pair of properties or a set of distinct molecules sharing an array of properties. To explore how this intrinsic flexibility manifests in the molecular design process, we used multi-objective optimization to search for molecules with simultaneously large polarizabilities and HOMO–LUMO gaps; analysis of the resulting Pareto fronts identified non-trivial paths through CCS consisting of sequential structural and/or compositional changes that yield molecules with optimal combinations of these properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In silico design of microporous polymers for chemical separations and storage

Polymers of intrinsic microporosity (PIMs) are a family of materials with potential to be effective and scalable solutions for challenging adsorbent and membrane applications. The broad range of repeat unit chemistry, microporous structural features, and polymer processing makes exploration of the expansive PIM design space inefficient via chemical and materials intuition alone. Computational techniques such as molecular simulations and machine learning can provide a leap in capabilities to address this polymer design challenge and will be central to the future development of PIMs. In this work, we highlight recent microporous material studies that arrived at key results by employing computational techniques and provide our perspective on the prospects for in silico design and development of PIMs.

adsorption↗

In-Silico Design of Next Generation Cellulose-Derived Packaging Materials (CRADA Final Report)

Developing sustainable solutions for single-use packaging is an important objective to combat the environmental crisis of plastics pollution. Most embodiments of cellulose-based packaging materials, including CellophaneTM, are completely biodegradable in both terrestrial and marine environments. However, petroleum-derived alternatives offer some performance advantages for metrics such as moisture barriers and mechanical properties. This project leverages molecular dynamics simulation to investigate how molecular modifications to cellulose-based polymer assemblies impact their material properties. An important performance criterion for the modified materials was to retain biodegradability; thus, modifications by naturally occurring, biodegradable additives were the focus of this study. Specifically, we developed models with xylan and lignin of varying monomeric compositions into the cellulose matrix. The mechanical properties were investigated by performing stress-strain simulations, and the water barrier and hydrophobicity were investigated by simulating the water contact angle. Our findings indicate that the incorporation of xylan into the cellulose matrix tends to increase the mechanical properties with an optimal loading of ~27 wt%. We also predict that orienting the nanoscale directionality of the xylan chains such that they are perpendicular to the cellulose fibrils will dramatically increase mechanical strength. In contrast, the incorporation of lignin tends to weaken the composite at all loadings investigated. Simulations of water contact angle predicted that coating polymers on the surface of the cellulose assembly creates a more hydrophobic surface than incorporating them throughout the matrix. Of the coatings investigated, lignin resulted in the most hydrophobic surface, followed by pectin and keratin, which both imparted modest increases in hydrophobicity. Future experimental work done by Futamura will focus on designing material prototypes to capitalize on the predictions of performance enhancement obtained from molecular modeling. While substantial progress was made by the simulations performed in this project, there still exists a vast parameter space that we were unable to investigate, including branching, functional group decoration, and degree of polymerization of polymer additives. However, the methods developed in this initial investigation will facilitate more rapid evaluation of the impact of molecular characteristics on the performance of biopolymer composite materials and thereby accelerate future materials discovery efforts in this area.

36 MATERIALS SCIENCE↗

High-throughput native mass spectrometry as experimental validation for in silico drug design

In this project, we developed automated workflows for both experimental validation and computational prediction of protein-ligand interactions. The ultimate goal is to establish an integrated pipeline for high throughput design of inhibitors to enzymes relevant to all areas of biological research. Our experimental approach is based on native mass spectrometry (native MS), which measures accurate masses and quantify the relative abundance of protein-ligand complexes to define binding affinity. We set up an in-house built autosampler with highly flexible configurations to minimize the manual steps for high throughput native MS. In parallel, we also performed manual native MS to characterize the binding of substrates and inhibitors of SARS-Cov-2 nonstructural protein nsp10/16 in order to optimize the experimental parameters for future automation. On the computational side, we streamlined the pipeline to achieve minimal manual intervention for predicting enzyme inhibitors via simulation, using the same nsp10/16 system as an example. Using the native MS method we examined 8 top-ranked designed compounds, 2 of which showed weak binding of ~50 µM. The information from native MS experiment provided critical insights and the foundation for a fully integrated workflow for enzyme inhibitor design.

59 BASIC BIOLOGICAL SCIENCES↗

gRNA-SeqRET: a universal tool for targeted and genome-scale gRNA design and sequence extraction for prokaryotes and eukaryotes

High-throughput genetic screening is frequently employed to rapidly associate gene with phenotype and establish sequence-function relationships. With the advent of CRISPR technology, and the ability to functionally interrogate previously genetically recalcitrant organisms, non-model organisms can be investigated using pooled guide RNA (gRNA) libraries and sequencing-based assays to quantitatively assess fitness of every targeted locus in parallel. To aid the construction of pooled gRNA assemblies, we have developed an in silico design workflow for gRNA selection using the gRNA Sequence Region Extraction Tool (gRNA-SeqRET). Built upon the previously developed CCTop, gRNA-SeqRET enables automated, scalable design of gRNA libraries that target user-specified regions or whole genomes of any prokaryote or eukaryote. Additionally, gRNA-SeqRET automates the bulk extraction of any regions of sequence relative to genes or other features, aiding in the design of homology arms for insertion or deletion constructs. We also assess in silico the application of a designed gRNA library to other closely related genomes and demonstrate that for very closely related organisms Average Nucleotide Identity (ANI) > 95% a large fraction of the library may be of relevance. The gRNA-SeqRET web application pipeline can be accessed at https://grna.jgi.doe.gov. The source code is comprised of freely available software tools and customized Python scripts, and is available at https://bitbucket.org/berkeleylab/grnadesigner/src/master/ under a modified BSD open-source license (https://bitbucket.org/berkeleylab/grnadesigner).

59 BASIC BIOLOGICAL SCIENCES↗

De novo designed protein inhibitors of amyloid aggregation and seeding

Neurodegenerative diseases are characterized by the pathologic accumulation of aggregated proteins. Known as amyloid, these fibrillar aggregates include proteins such as tau and amyloid-β (Aβ) in Alzheimer’s disease (AD) and alpha-synuclein (αSyn) in Parkinson’s disease (PD). The development and spread of amyloid fibrils within the brain correlates with disease onset and progression, and inhibiting amyloid formation is a possible route toward therapeutic development. Recent advances have enabled the determination of amyloid fibril structures to atomic-level resolution, improving the possibility of structure-based inhibitor design. In this work, we use these amyloid structures to design inhibitors that bind to the ends of fibrils, “capping” them so as to prevent further growth. Using de novo protein design, we develop a library of miniprotein inhibitors of 35 to 48 residues that target the amyloid structures of tau, Aβ, and αSyn. Biophysical characterization of top in silico designed inhibitors shows they form stable folds, have no sequence similarity to naturally occurring proteins, and specifically prevent the aggregation of their targeted amyloid-prone proteins in vitro. The inhibitors also prevent the seeded aggregation and toxicity of fibrils in cells. In vivo evaluation reveals their ability to reduce aggregation and rescue motor deficits in Caenorhabditis elegans models of PD and AD.

59 BASIC BIOLOGICAL SCIENCES↗

Chemistry-mediated Ostwald ripening in carbon-rich C/O systems at extreme conditions

Abstract There is significant interest in establishing a capability for tailored synthesis of next-generation carbon-based nanomaterials due to their broad range of applications and high degree of tunability. High pressure (e.g., shockwave-driven) synthesis holds promise as an effective discovery method, but experimental challenges preclude elucidating the processes governing nanocarbon production from carbon-rich precursors that could otherwise guide efforts through the prohibitively expansive design space. Here we report findings from large scale atomistically-resolved simulations of carbon condensation from C/O mixtures subjected to extreme pressures and temperatures, made possible by machine-learned reactive interatomic potentials. We find that liquid nanocarbon formation follows classical growth kinetics driven by Ostwald ripening (i.e., growth of large clusters at the expense of shrinking small ones) and obeys dynamical scaling in a process mediated by carbon chemistry in the surrounding reactive fluid. The results provide direct insight into carbon condensation in a representative system and pave the way for its exploration in higher complexity organic materials. They also suggest that simulations using machine-learned interatomic potentials could eventually be employed as in-silico design tools for new nanomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Architector for high-throughput cross-periodic table 3D complex building

Abstract Rare-earth and actinide complexes are critical for a wealth of clean-energy applications. Three-dimensional (3D) structural generation and prediction for these organometallic systems remains a challenge, limiting opportunities for computational chemical discovery. Here, we introduce Architector , a high-throughput in-silico synthesis code for s-, p-, d-, and f-block mononuclear organometallic complexes capable of capturing nearly the full diversity of the known experimental chemical space. Beyond known chemical space, Architector performs in-silico design of new complexes including any chemically accessible metal-ligand combinations. Architector leverages metal-center symmetry, interatomic force fields, and tight binding methods to build many possible 3D conformers from minimal 2D inputs including metal oxidation and spin state. Over a set of more than 6,000 x-ray diffraction (XRD)-determined complexes spanning the periodic table, we demonstrate quantitative agreement between Architector-predicted and experimentally observed structures. Further, we demonstrate out-of-the box conformer generation and energetic rankings of non-minimum energy conformers produced from Architector , which are critical for exploring potential energy surfaces and training force fields. Overall, Architector represents a transformative step towards cross-periodic table computational design of metal complex chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

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↗

HIV Env assembling is tightly modulated by lipid composition and cytoplasmic tail integrity

HIV envelope glycoprotein (Env) is a homotrimeric transmembrane protein which binds to receptors CD4, CCR5- CXCR4, leading to a cascade of conformational changes ultimately resulting in the entrance of the virus into the host cell. Its whole structure comprises an extracellular region, a transmembrane domain and a cytoplasmic region. While many efforts are carried in order to develop an effective vaccine specifically targeting the extracellular domain, less efforts are focus into regions close to the transmembrane region; mostly limited by the lack of structural knowledge within this domain. Thus, it is desirable to have a molecular description of the trimeric association of the transmembrane domain, as well as the membrane properties affecting their structure. We have carried molecular dynamic (MD) simulations of HIV transmembrane region in different lipid environments, mimicking the complex viral membrane. Our results indicate that the lipid environment heavily affects the dynamics and association of the transmembrane helices, modulating its structure. In addition, we have investigated the effect of the cytoplasmic tail (CT). Unexpectedly, this region has a tremendous effect on the lateral organization of the transmembrane. These results clearly highlights the importance of this region in the overall structure of the Env protein and point towards its key role for the in-silico design of effective vaccines towards this domain.

59 BASIC BIOLOGICAL SCIENCES↗

Final Technical Report: In-Silico Heterogeneous Catalyst Design for GHG Reduction via Bulk Chemicals

This project has advanced the understanding of ammonia synthesis by developing novel thermochemical catalysts, setting a new benchmark for energy and carbon efficiency in large-scale bulk chemicals production. This research leverages a proprietary computational discovery platform and high-throughput experimental validation, significantly accelerating catalyst optimization and enabling ammonia synthesis at lower temperature and pressure conditions than conventional processes. As a result, green ammonia can be produced at a cost that meets or exceeds established targets, supporting a viable pathway toward decarbonized ammonia for fuel, fertilizers, and hydrogen supply applications.

Grose, Jacob E [Copernic Catalysts, Inc.]↗

Interleukin-2 superkines by computational design

Significance While computational engineering of therapeutic proteins is a desirable goal, in practice the optimization of protein–protein interactions requires substantial experimental intervention. We present here a computational approach that focuses on stabilizing core protein structures rather than engineering the protein–protein interface. Using this approach, we designed thermostabilized interleukin-2 (IL-2) variants that bind tightly to their receptor without experimental optimization, mimicking the properties of the yeast-display engineered IL-2 variant “super-2.” Our results suggest that structure-guided stabilization may be a general method for in silico affinity maturation of protein–protein interactions.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

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↗

A novel approach to build algal consortia for sustainable biomass production

In the last decade, microalgae have reemerged as a feedstock for biofuels and a diverse suite of bioproducts. Yet, considerable challenges must be overcome before algal biofuels and bioproducts become technoeconomically viable. At present, single algal strain selected for particular phenotypic traits, such as maximum specific growth rate or lipid content, are commonly scaled for cultivation in open, outdoor raceway ponds due to the low capital costs of these systems. Although this monoculture approach may maximize the production of end products, monocultures are particularly susceptible to crashes associated with environmental and biological variability. An approach that has been proposed to generate more productive and stable microalgal crops is the use of eco-engineered communities, or consortia. Yet, attempts to construct productive consortia have not been consistently successful. We argue that failures stem from the lack of an eco-engineering approach to design species combinations. Here, we used an in silico method to build consortia before testing their performance against monocultures. Focusing on consortia of Nannochloropsis and Microchloropsis, we measured growth of strains along gradients of light, temperature, and salinity and used a functional dispersion approach to generate over 8000 functionally-diverse consortia combinations. We tested the 50 most functionally diverse consortia in a laboratory experiment and found that consortia overwhelmingly outperformed monocultures. Indeed, overyielding (OY) and a positive net biodiversity effect (NBE) was found respetively in 8%-86% and 88-92% of consortia combinations over the different experimental phases. To our knowledge, this is the first application of an in silico approach to design functionally diverse consortia before laboratory and field testing. Furthermore, our results highlight the importance of employing a functional diversity approach for consortia design.

09 BIOMASS FUELS↗

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. In this study, we aim to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

60 APPLIED LIFE SCIENCES↗