Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular Modeling”

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 397 records · Page 22

3D-Scaffold: A Deep Learning Framework to Generate 3D Coordinates of Drug-like Molecules with Desired Scaffolds

The prerequisite of therapeutic drug design is to identify novel molecules with desired biophysical and biochemical properties. Deep generative models have demonstrated their ability to find such molecules by exploring a huge chemical space efficiently. An effective way to obtain molecules with desired target properties is the preservation of critical scaffolds in the generation process. To this end, we propose a domain aware generative framework called 3D-Scaffold that takes 3D coordinates of a desired scaffold as an input and generates 3D coordinates of novel therapeutic candidates as an output while always preserving the desired scaffolds in generated structures. We show that our framework generates predominantly valid, unique, novel, and experimentally synthesizable molecules that have drug-like properties similar to the molecules in the training set. Using domain specific datasets, we generate covalent and non-covalent antiviral inhibitors. Therefore, to measure the success of our framework in generating therapeutic candidates, generated structures were subjected to high throughput virtual screening via docking simulations, which shows favorable interaction against SARS-CoV-2 main protease and non-structural protein endoribonuclease (NSP15) targets. Most importantly, our model performs well with relatively small volumes of training data and generalizes to new scaffolds, making it applicable to other domain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elucidating molecular level interfacial interactions between a de novo protein and nucleated calcite with solid-state NMR

Biomineralization is the process by which organisms use biomolecules to produce hierarchically structured organic–inorganic composites. Using biology as inspiration, a protein construct (FD31) was previously designed to accelerate formation of nano-calcite with an unconventional {110} face. Here, to understand the molecular interactions essential for protein aided calcite nucleation, solid-state nuclear magnetic resonance (ssNMR) spectroscopy was used in this work to characterize the FD31–calcite interface at the atomic level. Glutamic acid side chains designed to interact directly with calcium ions on the surface were found to have dynamics on the sub-millisecond timescale, indicating possible interactions between the protein and surface waters that were not included in the original model. Dipolar ssNMR recoupling techniques also showed that the protein backbone is ∼2 Å closer to the surface than in the original docking model. Refined molecular simulations were done in the presence of explicit waters, which resulted in the protein backbone closer to the surface than in the original docking structure, providing better agreement with experiment and highlighting the important role played by water in FD31–calcite interactions. While this work provides the first experimental confirmation that FD31 interactions with calcite are localized to the surface of the protein designed to serve as a template, these studies do indicate a more dynamic binding and closer binding mode between FD31 and the nucleated surface than originally proposed. In all, this enhanced molecular insight into the FD31–calcite interface has advanced our fundamental understanding of the atomic interactions at the organic–inorganic interface and will aid in the design of biological templates for the nucleation of inorganic crystals.

Close, Emily G. S. [Pacific Northwest National Lab↗

DataSet for Elucidating molecular level interfacial interactions between a de novo protein and nucleated calcite with solid-state NMR

Biomineralization is the process by which organisms use biomolecules to produce hierarchically structured organic-inorganic composites. Using biology as inspiration, a protein construct (FD31) was previously designed to accelerate formation of nano-calcite with an unconventional {110} face. To understand the molecular interactions essential for protein aided calcite nucleation, solid-state nuclear magnetic resonance (ssNMR) spectroscopy was used in this work to characterize the FD31-calcite interface at the atomic level. Glutamic acid side chains designed to interact directly with calcium ions on the surface were found to have dynamics on the sub-millisecond timescale, indicating possible interactions between the protein and surface waters that were not included in the original model. Dipolar ssNMR recoupling techniques also showed that the protein backbone is ~2 Å closer to the surface than in the original docking model. Refined molecular simulations were done in the presence of explicit waters, which resulted in the protein backbone closer to the surface than in the original docking structure, providing better agreement with experiment and highlighting the important role played by water in FD31-calcite interactions. These studies provide the first experimental evidence to confirm that FD31 interactions with calcite are localized to the surface of the protein designed to serve as a template. However, these studies do indicate a more dynamic binding and closer binding mode between FD31 and the nucleated surface than originally proposed. In all, this enhanced molecular insight into the FD31-calcite interface has advanced our fundamental understanding of the atomic interactions at the organic-inorganic interface and will aid in the design of biological templates for the nucleation of inorganic crystals.

Saccuzzo Close, Emily Grace [Pacific Northwest Nat↗

Thermochemical Data Fusion Using Graph Representation Learning

Large databases are required for “Big Data” applications in catalysis and materials science. Thermochemical databases can be created by combining data from various sources and by correcting low-fidelity datasets to higher accuracy with minimal computation. To achieve this “data fusion”, thermochemical quantities of interest, calculated at various levels of density functional theory (DFT), need to be mapped to the same, high levels of theory. In this work, a graph theoretical, statistical framework is proposed for such tasks. Subgraph frequencies are shown to provide a natural representation for learning these fusion maps. The maps are linear and are learnt with automated descriptor selection. Using a dataset of as few as ~1% from the QM9 database of 133,885 molecules, these models can predict multiple thermochemical quantities at a higher level of theory with an accuracy of 1 kcal/mol. Here, the method is explainable, generalizable, and provides a diagnostic tool for outlier identification

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strategies for using membrane-based separations to extract critical metals from waste streams

Critical metals are currently extracted by mining followed by their purification. These processes are costly and not environmentally very desirable. In this perspective paper we discuss the potential of extracting these critical metals from a range waste-streams available in abundance globally. These waste streams include brine from desalination plants, effluents from oil drilling and hydraulic fracturing, as well as discharges from various industrial processes such as metal finishing, electroplating, mining, and chemical manufacturing. We show that with a range of new separation processes being developed their separation is showing potential of being both technologically and economically feasible. We also show how high performance computing can be combined with computational models to screen and accelerate the development of new technologies for extracting critical metals from waste streams.

36 MATERIALS SCIENCE↗

A Molecular Description of Hydrogel Forming Polymers for Cement-Based Printing Paste Applications

This research endeavors to link the physical and chemical characteristics of select polymer hydrogels to differences in printability when used as printing aids in cement-based printing pastes. A variety of experimental probes including differential scanning calorimetry (DSC), NMR-diffusion ordered spectroscopy (DOSY), quasi-elastic neutron scattering (QENS) using neutron backscattering spectroscopy, and X-ray powder diffraction (XRD), along with molecular dynamic simulations, were used. Conjectures based on objective measures of printability and physical and chemical-molecular characteristics of the polymer gels are emerging that should help target printing aid selection and design, and mix formulation. Molecular simulations were shown to link higher hydrogen bond probability and larger radius of gyration to higher viscosity gels. Furthermore, the higher viscosity gels also produced higher elastic properties, as measured by neutron backscattering spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PDB ‐101: Educational resources supporting molecular explorations through biology and medicine

Abstract The Protein Data Bank (PDB) archive is a rich source of information in the form of atomic‐level three‐dimensional (3D) structures of biomolecules experimentally determined using macromolecular crystallography, nuclear magnetic resonance (NMR) spectroscopy, and electron microscopy (3DEM). Originally established in 1971 as a resource for protein crystallographers to freely exchange data, today PDB data drive research and education across scientific disciplines. In 2011, the online portal PDB‐101 was launched to support teachers, students, and the general public in PDB archive exploration ( pdb101.rcsb.org ). Maintained by the Research Collaboratory for Structural Bioinformatics PDB, PDB‐101 aims to help train the next generation of PDB users and to promote the overall importance of structural biology and protein science to nonexperts. Regularly published features include the highly popular Molecule of the Month series, 3D model activities, molecular animation videos, and educational curricula. Materials are organized into various categories (Health and Disease, Molecules of Life, Biotech and Nanotech, and Structures and Structure Determination) and searchable by keyword. A biennial health focus frames new resource creation and provides topics for annual video challenges for high school students. Web analytics document that PDB‐101 materials relating to fundamental topics (e.g., hemoglobin, catalase) are highly accessed year‐on‐year. In addition, PDB‐101 materials created in response to topical health matters (e.g., Zika, measles, coronavirus) are well received. PDB‐101 shows how learning about the diverse shapes and functions of PDB structures promotes understanding of all aspects of biology, from the central dogma of biology to health and disease to biological energy.

Zardecki, Christine↗

Decoding the Effect of Anion Identity on the Solubility of N -(2-(2-Methoxyethoxy)ethyl)phenothiazine (MEEPT)

Variations in the solubility of redox-active organic molecules (ROM) of interest for nonaqueous redox flow batteries (RFB), especially as the ROM state-of-charge changes during charge-discharge cycling, present significant molecular design challenges. The situation is further complicated as ROM solubility can be regulated by the choice of electrolyte salt and solvent that together with the ROM comprise the catholyte or anolyte (redox electrolyte) formulation, presenting materials design challenges. The ROM N-(2-(2-methoxyethoxy)ethyl)phenothiazine (MEEPT) is a viscous liquid at room temperature and is miscible in several organic solvents, including acetonitrile and propylene carbonate. The MEEPT radical cation (MEEPT +· ) paired with tetrafluoroborate (BF 4 - ) in acetonitrile presents a 0.5 M solubility, a dramatic decrease when compared to the viscous liquid of neutral MEEPT. Here, in this study, we present a joint experimental, regression modeling, and molecular dynamics (MD) simulations investigation to explore MEEPT-X (where X represents the counteranion) salt solubility variability as a function of concentration and counteranion chemistry in acetonitrile. We find a strong dependence of the salt solubility on the counteranion and relate these findings to explicit intermolecular interactions between MEEPT +· and the counteranion in the electrolyte solution.

Perera, Anton S. [Univ. of Kentucky, Lexington, KY↗

A lightweight method for evaluating in situ workflow efficiency

Performance evaluation is crucial to understanding the behavior of scientific workflows. In this study, we target an emerging type of workflow, called in situ workflows. These workflows tightly couple components such as simulation and analysis to improve overall workflow performance. To understand the tradeoffs of various configurable parameters for coupling these heterogeneous tasks, namely simulation stride, and component placement, separately monitoring each component is insufficient to gain insights into the entire workflow behavior. Through an analysis of the state-of-the-art research, we propose a lightweight metric, derived from a defined in situ step, for assessing resource usage efficiency of an in situ workflow execution. By applying this metric to a synthetic workflow, which is parameterized to emulate behaviors of a molecular dynamics simulation, we explore two possible scenarios (Idle Simulation and Idle Analyzer) for the characterization of in situ workflow execution. In addition to preliminary results from a recently published study [11], we further exploit the proposed metric to evaluate a practical in situ workflow with a real molecular dynamics application, i.e., GROMACS. Here, experimental results show that the in transit placement (analytics on dedicated nodes) sustains a higher frequency for performing in situ analysis compared to the helper-core configuration (analytics co-allocated with simulation).

97 MATHEMATICS AND COMPUTING↗

AlphaFold Protein Structure Database for Sequence-Independent Molecular Replacement

Crystallographic phasing recovers the phase information that is lost during a diffraction experiment. Molecular replacement is a commonly used phasing method for crystal structures in the protein data bank. In one form it uses a protein sequence to search a structure database to find suitable templates for phasing. However, sequence information is not always available, such as when proteins are crystallized with unknown binding partner proteins or when the crystal is of a contaminant. The recent development of AlphaFold published the predicted protein structures for every protein from twenty distinct species. In this work, we tested whether AlphaFold-predicted E. coli protein structures were accurate enough to enable sequence-independent phasing of diffraction data from two crystallization contaminants of unknown sequence. Using each of more than 4000 predicted structures as a search model, robust molecular replacement solutions were obtained, which allowed the identification and structure determination of YncE and YadF. Our results demonstrate the general utility of the AlphaFold-predicted structure database with respect to sequence-independent crystallographic phasing.

59 BASIC BIOLOGICAL SCIENCES↗

Polymers for Extreme Conditions Designed Using Syntax-Directed Variational Autoencoders

We report the design/discovery of new materials is highly nontrivial owing to the near-infinite possibilities of material candidates and multiple required property/performance objectives. Thus, machine learning tools are now commonly employed to virtually screen material candidates with desired properties by learning a theoretical mapping from material-to-property space, referred to as the forward problem. However, this approach is inefficient and severely constrained by the candidates that the human imagination can conceive. Thus, in this work on polymers, we tackle the materials discovery challenge by solving the inverse problem: directly generating candidates that satisfy desired property/performance objectives. We utilize syntax-directed variational autoencoders (VAE) in tandem with Gaussian process regression (GPR) models to discover polymers expected to be robust under three extreme conditions: (1) high temperatures, (2) high electric field, and (3) high temperature and high electric field, useful for critical structural, electrical, and energy storage applications. This approach to learn from and augment) human ingenuity is general and can be extended to discover polymers with other targeted properties and performance measures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Comparison of Exact and Model Exchange–Correlation Potentials for Molecules

Accurate exchange–correlation (XC) potentials for three-dimensional systems–via solution of the inverse density functional theory (DFT) problem–are now available to test the quality of DFT approximations. Herein, the exact XC potential for seven molecules–dihydrogen at four different bond-lengths, lithium hydride, water, and ortho-benzyne–are computed from full configuration interaction reference densities. These are compared to model XC potentials from nonlocal (B3LYP, HSE06, SCAN0, and M08-HX) and semilocal/local (SCAN, PBE, and PW92) XC functionals. Whereas for most systems, relative errors in the ground-state densities are O(10 –3 –10 –2 ) , the model XC potentials have much higher errors of O(10 –1 –10 0 ) . Among the model XC functionals, SCAN0 offers the best agreement with the exact XC potential, underlining the significance of satisfying exact conditions as well as including nonlocal effects in XC functionals. This work indicates that tests against the exact XC potential will provide a promising new direction for building more accurate XC functionals for DFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of software methods for estimating protein-protein relative binding affinities

A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83–98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.

59 BASIC BIOLOGICAL SCIENCES↗

PDB‐101: Molecular Explorations through Biology and Medicine

PDB‐101 is an online portal for teachers, students, and the general public to promote exploration of the structural biology of proteins and nucleic acids ( pdb101.rcsb.org ). Learning about the diverse shapes and functions of these biological macromolecules helps to understand all aspects of biomedicine and agriculture, from protein synthesis to health and disease to biological energy. Why PDB‐101? Researchers around the world are studying these molecules at the atomic level. These 3D structures are freely available at the Protein Data Bank (PDB), the central storehouse of biomolecular structures. This website builds introductory materials to help beginners get started in the basics of biomolecular structure and function (“101”, as in an entry level course) as well as resources for extended learning. Since 2011, PDB‐101 has been developed by the RCSB PDB , a global resource for the advancement of research and education in biology and medicine. Along with our Worldwide PDB collaborators, RCSB PDB curates, annotates, and makes publicly available the PDB data deposited by scientists around the globe. The RCSB PDB then provides a window to these data through a rich online resource with powerful searching, reporting, and visualization tools for researchers. This information is then streamlined for students and teachers at PDB‐101. Features include the ongoing Molecule of the Month series, educational materials such as paper models, posters, molecular animations, educational curricula and more. The section “Guide to Understanding PDB Data” is a primer for detailed PDB‐specific information: PDB Data, Visualizing Structures, Reading Coordinate Files, scientific methods for structure determination, and more. PDB‐101 also runs annual Video Challenges for high school students. Participants create short videos that tell molecular stories that connect structural biology and medicine. Previous topics have included HIV/AIDS, diabetes, and antimicrobial resistance. The 2022 challenge will focus on Molecular Mechanisms of Cancer. PDB‐101 activities are evaluated using user surveys, feedback from in‐person activities, and website analytics. In 2020, PDB‐101 hosted >850,000 users and >2.6 million page views.

Zardecki, Christine↗

Accurate Machine Learning for Predicting the Viscosities of Deep Eutectic Solvents

Deep eutectic solvents (DESs) are emerging as environmentally friendly designer solvents for mass transport and heat transfer processes in industrial applications; however, the lack of accurate tools to predict and thus control their viscosities under both a range of environmental factors and formulations hinders their general application. While DESs may serve as designer solvents, with nearly unlimited combinations, this unfortunately makes it experimentally infeasible to comprehensively measure the viscosities of all DESs of potential industrial interest. To assist in the design of DESs, we have developed several new machine learning (ML) models that accurately and rapidly predict the viscosities of a diverse group of DESs at different temperatures and molar ratios using, to date, one of the most comprehensive data sets containing the properties of over 670 DESs over a wide range of temperatures (278.15–385.25 K). Three ML models, including support vector regression (SVR), feed forward neural networks (FFNNs), and categorical boosting (CatBoost), were developed to predict DES viscosity as a function of temperature and molar ratio and contrasted with multilinear and two-factor polynomial regression baselines. Further, quantum chemistry-based, COSMO-RS-derived sigma profile (σ-profile) features were used as inputs for the ML models. The CatBoost model is excellent at externally predicting DES viscosity, as indicated by high R 2 (0.99) and low root-mean-square-error (RMSE) and average absolute relative deviations (AARD) (5.22%) values for the testing data sets, and 98% of the data points lie within the 15% of AARD deviations. Furthermore, SHapley additive explanation (SHAP) analysis was employed to interpret the ML results and rationalize the viscosity predictions. The result is an ML approach that accurately predicts viscosity and will aid in accelerating the design of appropriate DESs for industrial applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intramolecular hydrogen bonding analysis

The quasi-atomic orbital (QUAO) bonding analysis is used to study intramolecular hydrogen bonding (IMHB) in salicylic acid and an intermediate that is crucial to the synthesis of aspirin. The bonding analysis rigorously explores IMHB through directly accessing information that is intrinsic to the molecular wave function, thereby bypassing the need for intrinsically biased methods. The variables that affect the strength of IMHB are determined using kinetic bond orders, QUAO populations, and QUAO hybridizations. Important properties include both the interatomic distance between hydrogen and oxygen participating in the IMHB and the hybridization on the oxygen. The bonding analysis further shows that each intramolecular hydrogen bond is a four-electron three-center bond. Here, the bonding analysis is used to understand how aromatic reactivity changes due to the effect of functional groups on the aromatic ring.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stoichiometry dependent properties of cerium hydride: An active learning developed interatomic potential study

Cerium hydride has a variety of interesting properties, including a known lattice contraction and densification with increasing hydrogen content. However, precise stoichiometric control is not experimentally straightforward and ab initio approaches are not computationally feasible for many properties such as melting and low temperature diffusion. Therefore, we develop a machine-learned interatomic potential for cerium hydride that is valid for H to Ce ratios from 2.0 to 3.0. A query-by-committee active learning approach is used to develop the training set. Leveraging classical molecular dynamics simulations, we assess a range of properties and provide fundamental mechanisms for the trends with stoichiometry. Finally, a majority of the properties follow the trend of lattice contraction, being governed by the stronger lattice binding induced by adding octahedral atoms.

36 MATERIALS SCIENCE↗