Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular dynamics computer simulation”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Structural and Dynamic Heterogeneity of Deep Eutectic Solvents Composed of Choline Chloride and Ortho-Phenol Derivatives

Structural, thermal, and dynamic properties of four deep eutectic solvents comprising choline chloride paired with ortho-phenolic derivative hydrogen-bond donors were probed using experiments and molecular simulations. The hydrogen-bond donors include phenol, catechol, o-chlorophenol, and o-cresol, in a 3:1 mixture with the hydrogen-bond acceptor choline chloride. Density, viscosity, and pulsed-field gradient NMR diffusivity measurements were conducted over a range of temperatures. Classical and ab initio molecular dynamics simulation results match experimental data reasonably well. Furthermore, the simulation results were then used to perform a more detailed analysis of the local structure and dynamics of these systems.

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

Structural Characterization of TRAF6 N-Terminal for Therapeutic Uses and Computational Studies on New Derivatives

Tumor necrosis factor receptor-associated factors (TRAFs) are a protein family with a wide variety of roles and binding partners. Among them, TRAF6, a ubiquitin ligase, possesses unique receptor binding specificity and shows diverse functions in immune system regulation, cellular signaling, central nervous system, and tumor formation. TRAF6 consists of an N-terminal Really Interesting New Gene (RING) domain, multiple zinc fingers, and a C-terminal TRAF domain. TRAF6 is an important therapeutic target for various disorders and structural studies of this protein are crucial for the development of next-generation therapeutics. Here, we presented a TRAF6 N-terminal structure determined at the Turkish light source “Turkish DeLight” to be 3.2 Å resolution at cryogenic temperature (PDB ID: 8HZ2). This structure offers insight into the domain organization and zinc-binding, which are critical for protein function. Since the RING domain and the zinc fingers are key targets for TRAF6 therapeutics, structural insights are crucial for future research. Separately, we rationally designed numerous new compounds and performed molecular docking studies using this template (PDB ID:8HZ2). According to the results, 10 new compounds formed key interactions with essential residues and zinc ion in the N-terminal region of TRAF6. Molecular dynamic (MD) simulations were performed for 300 ns to evaluate the stability of three docked complexes (compounds 256, 322, and 489). Compounds 256 and 489 was found to possess favorable bindings with TRAF6. These new compounds also showed moderate to good pharmacokinetic profiles, making them potential future drug candidates as TRAF6 inhibitors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Environmental Contributions to Proton Sharing in Protein Low-Barrier Hydrogen Bonds

Hydrogen bonds (H-bonds) are central to biomolecular structure and dynamics. Although H-bonds are typically characterized by well-defined proton positions, proton delocalization has been proposed to play a role in facilitating enzyme catalysis and allostery in some systems. Experimentally locating protons is difficult, hampering the study of proton mobility in H-bonds. We used neutron crystallography, atomic resolution X-ray bond length analysis, and large quantum mechanics/molecular mechanics-Born–Oppenheimer molecular dynamics (QM/MM-BOMD) simulations to comprehensively characterize the shared proton/deuteron in a Glu–Asp low-barrier hydrogen bond (LBHB) in the bacterial protein YajL that is a conventional H-bond in the homologous disease-associated human protein DJ-1. X-ray bond length analysis of protiated and perdeuterated DJ-1 and YajL shows no significant effect of deuteron substitution on these carboxylic acid-carboxylate H-bonds but does reveal an effect at the active site glutamic acid near a cysteine thiolate. Residues in an H-bonded network that might favor LBHB formation in YajL were interrogated by the mutation of homologous residues in DJ-1. A distal DJ-1 substitution increases proton delocalization in the Glu–Asp H-bond, demonstrating that mutations within extended H-bond networks can modulate proton transfer barriers in carboxylic acid-carboxylate H-bonds. In addition, proton mobility in the H-bond is correlated with dimer-spanning motions in the QM/MM-BOMD simulations of YajL and DJ-1. Our results show that proton delocalization can be tuned using combined bioinformatic, structural, and computational information, opening the possibility of using engineered proton delocalization as a probe of H-bonding environments and as a tool to test hypotheses about LBHB function.

Lin, Jiusheng [University of Nebraska, Lincoln, NE↗
DOE OSTI2021

InteraChem: Virtual Reality Visualizer for Reactive Interactive Molecular Dynamics

Interactive molecular dynamics in virtual reality (IMD-VR) simulations provide a digital molecular playground for students as an alternative or complement to traditional molecular modeling kits or 2D illustrations. Previous IMD-VR studies have used molecular mechanics to enable simulations of macromolecules such as proteins and nanostructures for the classroom setting with considerable success. Here, we present the InteraChem molecular visualizer, intended for reactive IMD-VR simulation using semiempirical and ab initio methods. InteraChem visualizes not only the molecular geometry but also (1) isosurfaces such as molecular orbitals and electrostatic potentials and (2) two-dimensional graphs of time-varying simulation quantities such as kinetic/potential energy, internal coordinates, and user-applied force. Additionally, InteraChem employs speech recognition to facilitate user interaction and introduces a novel “atom happiness” visualization using emojis to indicate the energetic feasibility of a particular bonding arrangement. Furthermore, we include a set of accompanying exercises that we have used to teach chemical reactivity in small molecular systems.

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

Molecular Viscosity and Diffusivity Effects in Transitional and Shock-Driven Mixing Flows

Here, this paper investigates the importance of molecular viscosity and diffusivity for the prediction of transitional and shock-driven mixing flows featuring high and low Reynolds and Mach number regions. Two representative problems are computed with implicit large-eddy simulations using the inviscid Euler equations (EE) and viscous Navier-Stokes equations (NSE): the Taylor-Green vortex at Reynolds number Re = 3000 and initial Mach number Ma = 0.28, and an air- SF 6 -air gas curtain subjected to two shock waves at Ma = 1.2. The primary focus is on differences between NSE and EE predictions due to viscous effects. The outcome of the paper illustrates the advantages of utilizing NSE. In contrast to the EE, where the effective viscosity decreases upon grid refinement, NSE predictions can be assessed for simulations of flows with transition to turbulence at prescribed constant Re. The NSE can achieve better agreement between solutions and reference data, and the results converge upon grid refinement. On the other hand, the EE predictions do not converge with grid refinement, and can only exhibit similarities with the NSE results at coarse grid resolutions. We also investigate the effect of viscous effects on the dynamics of the coherent and turbulent fields, as well as on the mechanisms contributing to the production and diffusion of vorticity. The results show that nominally inviscid calculations can exhibit significantly varying flow dynamics driven by changing effective resolution-dependent Reynolds number, and highlight the role of viscous processes affecting the vorticity field. These tendencies become more pronounced upon grid refinement. The discussion of the results concludes with the assessment of the computational cost of inviscid and viscous computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ 54 , and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ 54 binding loops. Based on the DnaB analysis, we hypothesize that DnaB ‘grabs’ ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

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

Temperature-dependent model of helium bubble expansion and bursting in tungsten

Continuum-scale models that can reliably predict the behavior of helium in tungsten are of interest to the fusion community due to the projected impact of these materials on fusion reactor operation. In this work, we perform molecular dynamics simulations of spherical helium bubbles of various sizes in tungsten at different temperatures and depths with the goal of determining a mathematical model of the pressure and density at which the bubbles initially expand or burst as a function of depth, size, temperature, and surface orientation. The bubbles are small enough that their loop-punching pressures cannot be accurately predicted with continuum mechanics, and their expansion behavior is important, as it appears to cause many of the features observed on helium-irradiated tungsten surfaces. We vary the temperature, bubble size, bubble depth, and surface orientation in each case, recording the bubble pressure and density that result in bubble expansion. An exponential function with three adjustable parameters is found to fit the results well; the parameters that best fit our results are provided.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Core energies of dislocations in bcc metals

Accurate methods and an efficient workflow for computing and documenting dislocation core energies are developed and applied to 1 2 ( 111 ) and ( 100 ) dislocations in five body-centered cubic (bcc) metals W, Ta, V, Mo, and α -Fe represented by 13 model interatomic potentials. For each dislocation type, dislocation core energies are extracted for a large number of dislocation characters thoroughly sampling the entire 2-space of crystallographic line orientations of the bcc lattice. Of particular interest, core energies of the 1 2 ( 111 ) { 110 } dislocations are found to be distinctly asymmetric with respect to the sign of the character angle, whereas core energies of ( 100 ) { 110 } junction dislocations exhibit marked cusps for line orientations vicinal to the closed-packed ( 111 ) directions. Overall, our findings furnish substantial insights for developing accurate models of dislocation core energies employed in mesoscale dislocation dynamics simulations of crystal plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular dynamics simulations of the secondary-binding site in disaccharide-modified glycopeptide antibiotics

Oritavancin is a semisynthetic glycopeptide antibiotic used to treat severe infections by multidrug-resistant Gram-positive pathogens. Oritavancin is known to be a thousand times more potent than vancomycin against Gram-positive bacteria due to the additional interactions with bacterial peptidoglycan (PG) facilitated by a secondary-binding site. The presence of this secondary-binding site is evident in desleucyl-oritavancin, an Edman degradation product of oritavancin, still retaining its potency against Gram-positive bacteria, whereas desleucyl-vancomycin is devoid of any antimicrobial activities. Herein, using explicit solvent molecular dynamics (MD) simulations, steered MD simulations, and umbrella sampling, we show evidence of a secondary-binding site mediated by the disaccharide-modified hydrophobic sidechain of oritavancin interactions with the pentaglycyl-bridge segment of the PG. The interactions were characterized through comparison to the interaction of PG with chloroeremomycin, vancomycin, and the desleucyl analogs of the glycopeptides. Our results show that the enhanced binding of oritavancin to PG over the binding of the other complexes studied is due to an increase in the hydrophobic effect, electrostatic and van der Waals interactions, and not the average number of hydrogen bonds. Our ranking of the binding interactions of the biomolecular complexes directly correlates with the order based on their experimental minimum inhibitory concentrations. The results of our simulations provide insight into the modification of glycopeptides to increase their antimicrobial activities or the design of novel antibiotics against pathogenic Gram-positive bacteria.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Identification of Small-Molecule Inhibitors of Fibroblast Growth Factor 23 Signaling via In Silico Hot Spot Prediction and Molecular Docking to α-Klotho

Fibroblast growth factor 23 (FGF23) is a therapeutic target for treating hereditary and acquired hypophosphatemic disorders, such as X-linked hypophosphatemic (XLH) rickets and tumor-induced osteomalacia (TIO), respectively. FGF23-induced hypophosphatemia is mediated by signaling through a ternary complex formed by FGF23, the FGF receptor (FGFR), and α-Klotho. Currently, disorders of excess FGF23 are treated with an FGF23-blocking antibody, burosumab. Small-molecule drugs that disrupt protein/protein interactions necessary for the ternary complex formation offer an alternative to disrupting FGF23 signaling. In this study, the FGF23:α-Klotho interface was targeted to identify small-molecule protein/protein interaction inhibitors since it was computationally predicted to have a large fraction of hot spots and two druggable residues on α-Klotho. We further identified Tyr433 on the KL1 domain of α-Klotho as a promising hot spot and α-Klotho as an appropriate drug-binding target at this interface. Subsequently, we performed in silico docking of ~5.5 million compounds from the ZINC database to the interface region of α-Klotho from the ternary crystal structure. Following docking, 24 and 20 compounds were in the final list based on the lowest binding free energies to α-Klotho and the largest number of contacts with Tyr433, respectively. Herein, five compounds were assessed experimentally by their FGF23-mediated extracellular signal-regulated kinase (ERK) activities in vitro, and two of these reduced activities significantly. Both these compounds were predicted to have favorable binding affinities to α-Klotho but not have a large number of contacts with the hot spot Tyr433. ZINC12409120 was found experimentally to disrupt FGF23:α-Klotho interaction to reduce FGF23-mediated ERK activities by 70% and have a half maximal inhibitory concentration (IC50) of 5.0 ± 0.23 μM. Molecular dynamics (MD) simulations of the ZINC12409120:α-Klotho complex starting from in silico docking poses reveal that the ligand exhibits contacts with residues on the KL1 domain, the KL1–KL2 linker, and the KL2 domain of α-Klotho simultaneously, thereby possibly disrupting the regular function of α-Klotho and impeding FGF23:α-Klotho interaction. ZINC12409120 is a candidate for lead optimization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

Folding and translocation of the undecamer of poly-L-leucine across the water-hexane interface. A molecular dynamics study

The undecamer of poly-L-leucine at the water-hexane interface is studied by molecular dynamics simulations. This represents a simple model relevant to folding and insertion of hydrophobic peptides into membranes. The peptide, initially placed in a random coil conformation on the aqueous side of the system, rapidly translocates toward the hexane phase and undergoes interfacial folding into an alpha-helix in the subsequent 36 ns. Folding is nonsequential and highly dynamic. The initially formed helical segment at the N-terminus of the undecamer becomes transiently broken and, subsequently, reforms before the remainder of the peptide folds from the C-terminus. The formation of intramolecular hydrogen bonds during the folding of the peptide is preceded by a dehydration of the participating polar groups, as they become immersed in hexane. Folding proceeds through a short-lived intermediate, a 3(10)-helix, which rapidly interconverts to an alpha-helix. Both helices contribute to the equilibrium ensemble of folded structures. The helical peptide is largely buried in hexane, yet remains adsorbed at the interface. Its preferred orientation is parallel to the interface, although the perpendicular arrangement with the N-terminus immersed in hexane is only slightly less favorable. In contrast, the reversed orientation is highly unfavorable, because it would require dehydration of C-terminus carbonyl groups that do not participate in intramolecular hydrogen bonding. For the same reason, the transfer of the undecamer from the interface to the bulk hexane is also unfavorable. The results suggest that hydrophobic peptides fold in the interfacial region and, simultaneously, translocate into the nonpolar side of the interface. It is further implied that peptide insertion into the membrane is accomplished by rotating from the parallel to the perpendicular orientation, most likely in such a way that the N-terminus penetrates the bilayer.

NASA Discipline Exobiology↗
DOE OSTI2023

Fluctuating hydrodynamics and the Rayleigh–Plateau instability

The Rayleigh–Plateau instability occurs when surface tension makes a fluid column become unstable to small perturbations. At nanometer scales, thermal fluctuations are comparable to interfacial energy densities. Consequently, at these scales, thermal fluctuations play a significant role in the dynamics of the instability. These microscopic effects have previously been investigated numerically using particle-based simulations, such as molecular dynamics (MD), and stochastic partial differential equation–based hydrodynamic models, such as stochastic lubrication theory. In this paper, we present an incompressible fluctuating hydrodynamics model with a diffuse-interface formulation for binary fluid mixtures designed for the study of stochastic interfacial phenomena. An efficient numerical algorithm is outlined and validated in numerical simulations of stable equilibrium interfaces. We present results from simulations of the Rayleigh–Plateau instability for long cylinders pinching into droplets for Ohnesorge numbers of Oh = 0.5 and 5.0. Both stochastic and perturbed deterministic simulations are analyzed and ensemble results show significant differences in the temporal evolution of the minimum radius near pinching. Short cylinders, with lengths less than their circumference, were also investigated. As previously observed in MD simulations, we find that thermal fluctuations cause these to pinch in cases where a perturbed cylinder would be stable deterministically. Finally, we show that the fluctuating hydrodynamics model can be applied to study a broader range of surface tension–driven phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Clustering Effects on the Structure of Ionomer Solutions: A Combined SANS and Simulations Study

Ionic assemblies, or clusters, determine the structure and dynamics of ionizable polymers and enable their many applications. Fundamental to attaining well-defined materials is controlling the balance between van der Waals interactions that govern the backbone behavior and the forces that drive the formation of ionic clusters. Here, using small angle neutron scattering and fully atomistic molecular dynamics simulations, the structure of a model ionomer, sulfonated polystyrene in toluene solutions, was investigated as the cluster cohesion was tweaked by the addition of ethanol. The static structure factor was measured by both techniques and correlated with the size of the ionic clusters as the polymer concentration was varied. The conjunction of SANS results and molecular insight from MD simulations enabled the determination of the structure in these inhomogeneous networks on multiple length scales. Additionally, we find that across the entire concentration range studied, a network driven by the formation of ionic clusters was formed, where the size of the clusters drives the inhomogeneity of these systems. Tweaking the ionic clusters through the addition of ethanol impacts the packing of the sulfonated groups, their shape, and their size distribution, which in turn, affects the structure of these networks.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High pressure suppression of plasticity due to an overabundance of shear embryo formation

Abstract High pressure shear band formation is a critical phenomenon in energetic materials due to its influence on both mechanical strength and mechanochemical activation. While shear banding is known to occur in a variety of these materials, the governing dynamics of the mechanisms are not well defined for molecular crystals. We conduct molecular dynamics simulations of shock wave induced shear band formation in the energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) to assess shear band nucleation processes. We find, that at high pressures, the initial formation sites for shear bands, “embryos”, form in excess and rapidly lower deviatoric stresses prior to shear band formation and growth. This results in the suppression of plastic deformation. A local cluster analysis is used to quantify and contrast this mechanism with a more typical shear banding seen at lower pressures. These results demonstrate a mechanism that is reversible in nature and that supersedes shear band formation at increased pressures. We anticipate that these results will have a broad impact on the modeling and development of high-strain rate application materials such as those for high explosives and hypersonic systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting permeation of compounds across the outer membrane of P. aeruginosa using molecular descriptors

The ability Gram-negative pathogens have at adapting and protecting themselves against antibiotics has increasingly become a public health threat. Data-driven models identifying molecular properties that correlate with outer membrane (OM) permeation and growth inhibition while avoiding efflux could guide the discovery of novel classes of antibiotics. Here we evaluate 174 molecular descriptors in 1260 antimicrobial compounds and study their correlations with antibacterial activity in Gram-negative Pseudomonas aeruginosa. The descriptors are derived from traditional approaches quantifying the compounds’ intrinsic physicochemical properties, together with, bacterium-specific from ensemble docking of compounds targeting specific MexB binding pockets, and all-atom molecular dynamics simulations in different subregions of the OM model. Using these descriptors and the measured inhibitory concentrations, we design a statistical protocol to identify predictors of OM permeation/inhibition. We find consistent rules across most of our data highlighting the role of the interaction between the compounds and the OM. An implementation of the rules uncovered in our study is shown, and it demonstrates the accuracy of our approach in a set of previously unseen compounds. Our analysis sheds new light on the key properties drug candidates need to effectively permeate/inhibit P. aeruginosa, and opens the gate to similar data-driven studies in other Gram-negative pathogens.

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

0BGRaman: Graph Network based Simulator for Forecasting Molecular Polarizability

This report presents the work performed under the GRaman project, sponsored by the PCSD LDRD Seed program. The project aimed at accelerating ab initio molecular dynamics simulation using Graph Networks. The Graph Network framework is a ML framework that has been successfully employed to simulate the dynamics of several physical systems: including water splashing in a container and flags moving with the wind. In this effort, we performed a data collection campaign for 3 different molecules of interest. We have built tools for preprocessing the trajectories obtained by simulating Raman Spectroscopy with NWChem and translating them into a suitable format for training. We have developed a training algorithm to train the Graph Network based simulators based on our data and developed a simulator that produces trajectories in the same NWChem format. While the tool has improved with each iteration of development and subsequent experiments, the current state of the tool does not allow to directly incorporate the technology within the NWChem framework because the trajectories produced by the tool are not yet accurate enough. However, the technology has proved to have good potential and it is certainly worth further research and development.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quantifying the atomistic free-volume morphology of materials with graph theory

Here, we introduce a new computational methodology for the identification and characterization of free volume within/around atomistic configurations. This scheme employs a three-stage workflow, by which spheres are iteratively grown inside of voxels, and ultimately converted to planar graphs, which are then characterized via a graph-based order parameter. Our approach is computationally efficient, physically intuitive, and universally transferable to any material system. Validation of our methodology is performed on several sets of materials problems: (1) classification of unique free volumes in various crystal phases, (2) autonomous detection and classification of complex surface defects during epitaxial growth simulations, (3) characterization of free volume defects in metals/alloys, and (4) quantification of the spatio-temporal behavior of nano-scale free volume morphologies as a function of both temperature and free-volume size. Our method accurately identifies and characterizes unique free volumes over a multitude of systems and length scales, indicating its potential for future use in understanding the relationship between free volume morphology and material properties under both static and dynamic conditions.

36 MATERIALS SCIENCE↗