Search NASA⌕ Search

SEARCH · Search NASA

Results for “Models, Molecular”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molecular Modeling and Molecular Dynamics Simulation of a Packed and Intact Bacterial Microcompartment

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins and are packed full of enzyme cargoes and metabolites as they fulfill their function. Coupling together recent SAXS and proteomics work, it is possible to develop molecular models for these microcompartments and interrogate enzyme and metabolite dynamics within. Our primary goal of this study is to quantify the permeability of metabolite glyceraldehyde-3-phosphate (G3P) and dihydroxyacetone phosphate (DHAP) across the BMC shell through classical molecular dynamics simulation. The Haliangium ochraceum model of BMC shell (PDB: 6MZX) was used to model an intact BMC of approximately 10 million atoms. Working at this scale presented its own challenges in managing large data sets, with multiple challenges and hardware advances discussed that facilitated this work. Over approximately 750 ns of aggregate simulation, we see multiple permeation events for these metabolites that were added at high concentration through the pores present within BMC shell tiles. When compared to independent permeability estimates for the same metabolites determined through replica exchange umbrella sampling simulations, the permeabilities varied by approximately 3 orders of magnitude. Regardless, the permeability coefficients for both G3P and DHAP are highly similar and very high, such that only very small concentration gradients can be maintained across the BMC shell between the cytosol and BMC interior. The large simulation systems also facilitated comparisons for molecular diffusivity in the crowded environment within the BMC shell. By our estimates, the viscosity within a packed BMC shell is at least 10-fold higher than it would be in neat solution and is the real driver for varying permeability estimates we obtained through simulation. These findings will be used as design inputs for future bioengineering efforts to make products from BMCs, highlighting how permeable BMC shells can be.

Diffusion↗
DOE OSTI2018

C2 adsorption in zeolites: in silico screening and sensitivity to molecular models

Efficient separation of ethane and ethylene has been a long-standing challenge for the chemical industry. In this work, we use molecular modeling to identify zeolite and zeotype frameworks that have the potential to be the next-generation solution for the separation of these C2 compounds. Using two different united-atom versions of the transferable potentials for phase equilibria (TraPPE) force field, the zeolitic structures in the database of the International Zeolite Association are screened for the separation of ethane and ethylene. A detailed analysis, with regards to accessibility of favorable sites and sensitivity to molecular models (also considering the explicit-hydrogen TraPPE model for ethane), is carried out on the top-performing structures. This study provides insights on the performance and limitations of molecular models for predicting mixture adsorption in zeolites.

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

Building molecular model series from heterogeneous CryoEM structures using Gaussian mixture models and deep neural networks

Cryogenic electron microscopy (CryoEM) produces structures of macromolecules at near-atomic resolution. However, building molecular models with good stereochemical geometry from those structures can be challenging and time-consuming, especially when many structures are obtained from datasets with conformational heterogeneity. Here we present a model refinement protocol that automatically generates series of molecular models from CryoEM datasets, which describe the dynamics of the macromolecular system and have near-perfect geometry scores. This method makes it easier to interpret the movement of the protein complex from heterogeneity analysis and to compare the structural dynamics observed from CryoEM data with results from other experimental and simulation techniques.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

2023 Roadmap on molecular modelling of electrochemical energy materials

New materials for electrochemical energy storage and conversion are the key to the electrification and sustainable development of our modern societies. Molecular modelling based on the principles of quantum mechanics and statistical mechanics as well as empowered by machine learning techniques can help us to understand, control and design electrochemical energy materials at atomistic precision. Therefore, this roadmap, which is a collection of authoritative opinions, serves as a gateway for both the experts and the beginners to have a quick overview of the current status and corresponding challenges in molecular modelling of electrochemical energy materials for batteries, supercapacitors, CO 2 reduction reaction, and fuel cell applications.

25 ENERGY STORAGE↗
DOE OSTI2020

Biological evaluation of Keggin‐type polyoxometalates on tyrosinase: Kinetics and molecular modeling

Abstract Abnormal overexpression of tyrosinase activity can lead to the production of hyperpigmentation in human skin and enzymatic browning in fruits and vegetables. Herein, the inhibition and mechanism of the H 3 PMo 12 O 40 and two transition metal‐substituted Keggin‐type polyoxometalates (Na 7 PMo 11 CoO 40 and Na 7 PMo 11 ZnO 40 ) on tyrosinase were studied by kinetics and molecular modeling. Kinetic studies indicated that all compounds had more potent inhibitory activities than standard arbutin, and H 3 PMo 12 O 40 (IC 50 = 0.443 ± 0.006 m m ) is ~15‐fold stronger inhibition than arbutin. Additionally, all compounds inhibited tyrosinase in a reversible competitive manner. Intriguingly, molecular modeling elucidated that three compounds competitively bind to tyrosinase mainly through more interactions with Cu 2+ ions and the amino acid residue capable of forming cation‐π and hydrogen bonding, forming a reversible non‐covalent complex. Molecular simulation study correlated well with the biological activity of three compounds in vitro. This work provided new insights into design and synthesis of polyoxometalates as tyrosinase inhibitors in the field of medicine, cosmetic, and food.

Chi, Guoxiang↗
DOE OSTI2024

Early Career Researcher Support for Foundations of Molecular Modeling and Simulation Conference

This project enabled 14 early career researchers, including graduate students and postdocs, to attend the 8th triennial conference on Foundations of Molecular Modeling and Simulation (FOMMS 2022), which took place in Delavan, WI (Lake Geneva area), from July 17 to 21, 2022. Molecular-level modeling and simulation, which serve as critical tools for advancing science and engineering in a myriad of applications, were at the forefront of discussions. These ranged from the design of new materials for various energy applications to insights into the interaction between drug molecules and protein receptors within the body.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Molecular model of J-aggregated pseudoisocyanine fibers

Aggregated cyanines form ordered supramolecular structures with the potential to transport energy efficiently over long distances, a hallmark of photosynthetic light-harvesting complexes. In concentrated aqueous solution, pseudoisocyanine (PIC) spontaneously forms fibers with a chiral J-band red-shifted 1600 cm -1 from the monomeric 0–0 transition. A cryogenic transmission electron microscopy analysis of these fibers show an average fiber width of 2.89 nm, although the molecular-level structure of the aggregate is currently unknown. To determine a molecular model for these PIC fibers, the calculated spectra and dynamics using a Frenkel exciton model are compared to experiment. A chiral aggregate model in which the PIC monomers are neither parallel nor orthogonal to the long axis of the fiber is shown to replicate the experimental spectra most closely. This model can be physically realized by the sequential binding of PIC dimers and monomers to the ends of the fiber. Furthermore, these insights into the molecular aggregation model for aqueous PIC can also be applied to other similar cyanine-based supramolecular complexes with the potential for long-range energy transport, a key building block for the rational design of novel excitonic systems

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

The comparison of semiempirical and ab initio molecular modeling methods in activity and property evaluation of selected antimicrobial sulfonamides

A group of sulfonamides exhibiting antimicrobial activity was analyzed with the application of quantitative structure–activity/property relationships method. The purpose of this study was to show the common and differentiating characteristics of the analyzed chemical structures alike physicochemically as well as pharmacologically based on the quantum-chemical calculations both in vacuo and in the aquatic environment together with their microbiological activity and chromatographic retention data. The semiempirical and ab initio level of in silico molecular modeling was performed for calculations of statistically significant molecular descriptors to compare obtained results. The relationship between the structure and biological activity and physicochemical parameters data was able to class and describe analyzed molecules and the applied chemometric approaches (principal component analysis, factor analysis, and multiple regression analysis) revealed the influential features of the tested structures responsible for the antimicrobial activity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Modeling Molecular Hydrogen in Low-metallicity Galaxies

We use a suite of hydrodynamics simulations of the interstellar medium (ISM) within a galactic disk, which includes radiative transfer, a nonequilibrium model of molecular hydrogen, and a realistic model for star formation and feedback, to study the structure of the ISM and H 2 abundance as a function of local ISM properties. We show that the star formation rate and structure of the ISM are sensitive to the metallicity of the gas with a progressively smoother density distribution with decreasing metallicity. In addition to the well-known trend of the H I–H 2 transition shifting to higher densities with decreasing metallicity, the maximum achieved molecular fraction in the ISM drops drastically at Z ≲ 0.2 Z ⊙ as the formation time of H 2 becomes much longer than a typical lifetime of dense regions of the ISM. We present accurate fitting formulae for both volumetric and projected ƒ H 2 measured on different scales as a function of gas metallicity, UV radiation field, and gas density. We show that when the formulae are applied to the patches in the simulated galaxy, the overall molecular gas mass is reproduced to better than a factor of ≲1.5 across the entire range of metallicities and scales. We also show that the presented fit is considerably more accurate than any of the previous ƒ H 2 models and fitting formulae in the low-metallicity regime. The fit can thus be used for modeling molecular gas in low-resolution simulations and semi-analytic models of galaxy formation in the dwarf and high-redshift regimes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Critical behavior in a chiral molecular model

Understanding the condensed-phase behavior of chiral molecules is important in biology as well as in a range of technological applications, such as the manufacture of pharmaceuticals. Here, we use molecular dynamics simulations to study a chiral four-site molecular model that exhibits a second-order symmetry-breaking phase transition from a supercritical racemic liquid into subcritical D-rich and L-rich liquids. We determine the infinite-size critical temperature using the fourth-order Binder cumulant, and we show that the finite-size scaling behavior of the order parameter is compatible with the 3D Ising universality class. We also study the spontaneous D-rich to L-rich transition at a slightly subcritical temperature of T = 0.985Tc, and our findings indicate that the free energy barrier for this transformation increases with system size as N2/3, where N is the number of molecules, consistent with a surface-dominated phenomenon. The critical behavior observed herein suggests a mechanism for chirality selection in which a liquid of chiral molecules spontaneously forms a phase enriched in one of the two enantiomers as the temperature is lowered below the critical point. Furthermore, the increasing free energy barrier with system size indicates that fluctuations between the L-rich and D-rich phases are suppressed as the size of the system increases, trapping it in one of the two enantiomerically enriched phases. Such a process could provide the basis for an alternative explanation for the origin of biological homochirality. We also conjecture the possibility of observing nucleation at subcritical temperatures under the action of a suitable chiral external field.

Chemistry↗
DOE OSTI2021

CASTELO: clustered atom subtypes aided lead optimization—a combined machine learning and molecular modeling method

Background: Drug discovery is a multi-stage process that comprises two costly major steps: pre-clinical research and clinical trials. Among its stages, lead optimization easily consumes more than half of the pre-clinical budget. We propose a combined machine learning and molecular modeling approach that partially automates lead optimization workflow in silico, providing suggestions for modification hot spots. Results: The initial data collection is achieved with physics-based molecular dynamics simulation. Contact matrices are calculated as the preliminary features extracted from the simulations. To take advantage of the temporal information from the simulations, we enhanced contact matrices data with temporal dynamism representation, which are then modeled with unsupervised convolutional variational autoencoder (CVAE). Finally, conventional and CVAE-based clustering methods are compared with metrics to rank the submolecular structures and propose potential candidates for lead optimization. Conclusion: With no need for extensive structure-activity data, our method provides new hints for drug modification hotspots which can be used to improve drug potency and reduce the lead optimization time. It can potentially become a valuable tool for medicinal chemists.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Molecular Models of Atomically Dispersed Uranium at MoS 2 Surfaces Reveal Cooperative Mechanism of Water Reduction

Single atoms of uranium supported on molybdenum sulfide surfaces (U@MoS 2 ) have been recently demonstrated to facilitate the hydrogen evolution reaction (HER) through electrocatalysis. Theoretical calculations have predicted uranium hydroxide moieties bound to edge-sulfur atoms of MoS 2 as a proposed transition state involved in the HER process. However, the isolation of relevant intermediates involved in this process remains a challenge, rendering mechanistic hypotheses unverified. The present work describes the isolation and characterization of a uranium-hydroxide intermediate on molybdenum sulfide surfaces using [(Cp* 3 Mo 3 S 4 )UCp*], a molecular model of a reduced uranium center supported at MoS 2 . Mechanistic investigations highlight the metalloligand cooperativity with uranium involved in the water activation pathway. The corresponding uranium-oxo analogue, [(Cp* 3 Mo 3 S 4 )Cp*U(═O)], was also accessed from the hydroxide cluster via hydrogen atom transfer and from [(Cp* 3 Mo 3 S 4 )UCp*] through an alternative direct oxygen atom transfer. These results provide an atomistic perspective on the reactivity of low-valent uranium at molybdenum sulfide surfaces toward water, modeling key intermediates associated with the HER of U@MoS 2 catalysts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2017

Azomethines, isoxazole, N-substituted pyrazoles and pyrimidine containing curcumin derivatives: Urease inhibition and molecular modeling studies

Curcumin has shown large number of pharmacological properties against different phenotypes of various disease models. Different synthetic routes have been employed to develop its various derivatives for diverse biological functions. In this study, curcumin derived azomethine, isoxazole, pyrimidines and N-substituted pyrazoles were synthesized to investigate their urease enzyme inhibition. The structures of newly synthesized compounds were described by IR, MS, {sup 1}H NMR and {sup 13}C NMR spectral data. Urease enzyme inhibition was evaluated through in vitro assays in which compound 8b was found to be the most potent (IC{sub 50} = 2.44 ± 0.07 μM) among the tested compounds. The compounds with diazine ring system except the 4d showed better urease inhibition (IC{sub 50} = 11.43 ± 0.21–19.63 ± 0.28 μM) than the standard urease inhibitor thiourea (IC{sub 50} = 22.61 ± 0.23 μM). Similarly enzyme kinetics data revealed that compounds 3c-3e and 8b were competitive inhibitors with Ki values of 20.0, 19.87, 20.23 and 19.11 μM respectively while the compounds 4b, 4c and 4e were mixed type of inhibitors with Ki values 6.72, 19.69 and 6.72 μM respectively. Molecular docking studies were also performed to identify the plausible binding modes of the most active compounds. - Highlights: • Azomethines, Isoxazole, N-substituted pyrazoles and pyrimidine with curcumin scaffold were synthesized. • Anti-urease studies of target compounds were performed. • Molecular modeling of most active compounds at active site of urease. • Thiourea conjugated with curcumin showed excellent urease inhibition with IC{sub 50} = 2.44 ± 0.07 μM.

60 APPLIED LIFE SCIENCES↗