Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Structural basis for severe pain caused by mutations in the voltage sensors of sodium channel NaV1.7

Voltage-gated sodium channels in peripheral nerves conduct nociceptive signals from nerve endings to the spinal cord. Mutations in voltage-gated sodium channel NaV1.7 are responsible for a number of severe inherited pain syndromes, including inherited erythromelalgia (IEM). Here, we describe the negative shifts in the voltage dependence of activation in the bacterial sodium channel NaVAb as a result of the incorporation of four different IEM mutations in the voltage sensor, which recapitulate the gain-of-function effects observed with these mutations in human NaV1.7. Crystal structures of NaVAb with these IEM mutations revealed that a mutation in the S1 segment of the voltage sensor facilitated the outward movement of S4 gating charges by widening the pathway for gating charge translocation. In contrast, mutations in the S4 segments modified hydrophobic interactions with surrounding amino acid side chains or membrane phospholipids that would enhance the outward movement of the gating charges. These results provide key structural insights into the mechanisms by which these IEM mutations in the voltage sensors can facilitate outward movements of the gating charges in the S4 segment and cause hyperexcitability and severe pain in IEM. Our work gives new insights into IEM pathogenesis at the near-atomic level and provides a molecular model for mutation-specific therapy of this debilitating disease.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

General framework for quantifying dissipation pathways in open quantum systems. III. Off-diagonal subsystem–bath couplings

This paper extends the previously reported theory of dissipation pathways [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111 (2024)] to incorporate off-diagonal subsystem–bath coupling, which is often required to model molecular systems where the environment directly influences transitions and couplings between subsystem states. We systematically derive master equations for both population transfer and dissipation into individual bath components, for which we also rigorously prove energy conservation and detailed balance. The approach is based on second-order perturbation theory with respect to the subsystem–bath couplings, whose form is not limited to any specific model. The accuracy of the developed method is tested by applying it to diverse model Hamiltonians involving linearly coupled harmonic oscillator baths and comparing the outcomes against the hierarchical equations of motion (HEOM) method. Overall, our method accurately quantifies the contributions of specific bath components to the overall dissipation while significantly reducing the computational cost compared to numerically exact methods such as HEOM, thus offering a path to examine how vibronic interactions steer non-adiabatic processes in realistic chemical systems.

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

Molecular mechanism of trehalose 6-phosphate inhibition of the plant metabolic sensor kinase SnRK1

SUCROSE-NON-FERMENTING1-RELATED PROTEIN KINASE1 (SnRK1), a central plant metabolic sensor kinase, phosphorylates its target proteins, triggering a global shift from anabolism to catabolism. Molecular modeling revealed that upon binding of KIN10 to GEMINIVIRUS REP-INTERACTING KINASE1 (GRIK1), KIN10’s activation T-loop reorients into GRIK1’s active site, enabling its phosphorylation and activation. Trehalose 6-phosphate (T6P) is a proxy for cellular sugar status and a potent inhibitor of SnRK1. T6P binds to KIN10, a SnRK1 catalytic subunit, weakening its affinity for GRIK1. Here, we investigate the molecular details of T6P inhibition of KIN10. Molecular dynamics simulations and in vitro phosphorylation assays identified and validated the T6P binding site on KIN10. Under high-sugar conditions, T6P binds to KIN10, blocking the reorientation of its activation loop and preventing its phosphorylation and activation by GRIK1. Under these conditions, SnRK1 maintains only basal activity levels, minimizing phosphorylation of its target proteins, thereby facilitating a general shift from catabolism to anabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Simulated Adhesive Response to Aging (CRADA 453 Final Report)

The goal of this project was to develop a suite of molecular models to investigate adhesive performance as a result of the presence of different levels of moisture. These different levels of moisture are accepted in the adhesives community to represent different states of aging. First, our team modeled the adhesive’s crosslinked structure. Next, we modeled adding water to the crosslinked structure and used high performance computing resources to simulate water diffusion through the bulk adhesive and possible sites for chain scission. Finally, we modeled the response of the simulated adhesive to tensile loads at various levels of moisture content. These simulations were validated by experimental characterization by the industry partner. We found that the Young’s modulus of the adhesive decreased by about 30% when the water content increased to 3 wt.%. Our simulations also identified the likely site where the polymer chain could break. Our integration of molecular level simulations and mechanical property predictions using high performance computing with experimental studies of a hydrated epoxy system has yielded fundamental insights that will help industry advance lightweight joining technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Correction to Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

This errata is to correct the acknowledgment in the original published article. G.B. (partial support), V.-A.G., M.-T.N., and R.R. were supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division. Computer resources were provided by Research Computing at Pacific Northwest National Laboratory and the National Energy Research Scientific Computing Center (NERSC), a U.S. Department of Energy Office of Science User Facility operated under Contract No. DE-AC02-05CH11231. PNNL is operated by Battelle for the U.S. Department of Energy under Contract DE-AC05-76RL01830. G.B. acknowledges the EPSRC for partial support in PhD funding. The authors acknowledge the use of the IRIDIS High Performance Computing Facility (IRIDIS 5) and associated support services at the University of Southampton in the completion of this work. We are grateful to the UK Materials and Molecular Modeling Hub (Thomas HPC) for computational resources, which is partially funded by EPSRC (EP/P020194/1).

continuum solvation, metal/water interfaces, ab in↗
DOE OSTI2023

Critical fluctuations in liquid–liquid extraction organic phases controlled by extractant and diluent molecular structure

Extractant aggregation in liquid–liquid extraction organic phases impacts extraction energetics and is related to the deleterious efficiency-limiting liquid–liquid phase transition known as third phase formation. Using small angle X-ray scattering, we find that structural heterogeneities across a wide range of compositions in binary mixtures of malonamide extractants and alkane diluents are well described by Ornstein–Zernike scattering. This suggests that structure in these simplified organic phases originates from the critical point associated with the liquid–liquid phase transition. To confirm this, we measure the temperature dependence of the organic phase structure, finding critical exponents consistent with the 3D Ising model. Molecular dynamics simulations were also consistent with this mechanism for extractant aggregation. Due to the absence of water or any other polar solutes required to form reverse-micellar-like nanostructures, these fluctuations are inherent to the binary extractant/diluent mixture. Further, we also show how the molecular structure of the extractant and diluent modulate these critical concentration fluctuations by shifting the critical temperature: critical fluctuations are suppressed by increasing extractant alkyl tail lengths or decreasing diluent alkyl chain lengths. This is consistent with how extractant and diluent molecular structure are known to impact metal and acid loading capacity in many-component LLE organic phases, suggesting phase behavior of practical systems may be effectively studied in simplified organic phases. Overall, the explicit connection between molecular structure, aggregation and phase behavior demonstrated here will enable the design of more efficient separations processes.

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

Heat Transport Hysteresis Generated Through Frequency Switching of a Time-Dependent Temperature Gradient

A stochastic energetics framework is applied to examine how periodically shifting the frequency of a time-dependent oscillating temperature gradient affects heat transport in a nanoscale molecular model. We specifically examine the effects that frequency switching, i.e., instantaneously changing the oscillation frequency of the temperature gradient, has on the shape of the heat transport hysteresis curves generated by a particle connected to two thermal baths, each with a temperature that is oscillating in time. Analytical expressions are derived for the energy fluxes in/out of the system and the baths, with excellent agreement observed between the analytical expressions and the results from nonequilibrium molecular dynamics simulations. We find that the shape of the heat transport hysteresis curves can be significantly altered by shifting the frequency between fast and slow oscillation regimes. We also observe the emergence of features in the hysteresis curves such as pinched loops and complex multi-loop patterns due to the frequency shifting. The presented results have implications in the design of thermal neuromorphic devices such as thermal memristors and thermal memcapacitors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phage K gp102 Drives Temperature-Sensitive Antibacterial Activity on USA300 MRSA

There is widespread interest in using obligately lytic bacteriophages (“phages”) to treat human bacterial infections. Among Staphylococcus aureus infections, the USA300 lineage is a frequent cause of invasive disease. We observed that phage K, a model S. aureus myophage, exhibits temperature-sensitive growth on USA300 strains, with the wild-type phage providing poorer growth suppression in broth and forming smaller and fainter plaques at 37 °C vs. 30 °C. We isolated 65 mutants of phage K that had improved plaquing characteristics at 37 °C when compared to the parental phage. In all 65 mutants, this phenotype was attributable to loss-of-function (LoF) mutations in gp102, which encodes a protein of unknown function that has homologs only among the Herelleviridae (SPO1-like myophages infecting gram-positive bacteria). Additional experiments with representative mutants consistently showed that the temperature-sensitive plaque phenotype was specific to USA300 MRSA strains and that Gp102 disruption was correlated with improved suppression of bacterial growth in broth and improved antibacterial activity in a mouse model of upper respiratory tract infection. The same genotype and in vitro phenotypes could be replicated in close relatives of phage K. Gp102 disruption did not have a detectable effect on adsorption but did delay cell culture lysis relative to wild-type under permissive infection conditions, suggesting that gp102 conservation might be maintained by selective pressure for more rapid replication. Expression of gp102 on a plasmid was toxic to both an MSSA and a USA300 MRSA strain. Molecular modeling predicts a protein with two helix-turn-helix domains that displays some similarity to DNA-binding proteins such as transcription factors. While its function remains unclear, gp102 is a conserved gene that is important to the infection process of Kayvirus phages, and it appears that the manner in which USA300 strains defend against them at 37 °C can be overcome by gp102 LoF mutations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Avoiding pitfalls in molecular simulation of vapor sorption: Example of propane and isobutane in metal–organic frameworks for adsorption cooling applications

This study introduces recommendations for conducting molecular simulations of vapor adsorption, with an emphasis on enhancing the accuracy, reproducibility, and comparability of results. The first aspect we address is consistency in the implementation of some details of typical molecular models, including tail corrections and cutoff distances, due to their significant influence on generated data. We highlight the importance of explicitly calculating the saturation pressures at relevant temperatures using methods such as Gibbs ensemble Monte Carlo simulations and illustrate some pitfalls in extrapolating saturation pressures using this method. For grand canonical Monte Carlo (GCMC) simulations, the input fugacity is usually calculated using an equation of state, which often requires the critical parameters of the fluid. We show the importance of using critical parameters derived from the simulation with the same model to ensure internal consistency between the simulated explicit adsorbate phase and the implicit bulk phase in GCMC. We show the advantages of presenting isotherms on a relative pressure scale to facilitate easier comparison among models and with experiment. Extending these guidelines to a practical case study, we evaluate the performance of various isoreticular metal–organic frameworks (MOFs) in adsorption cooling applications. This includes examining the advantages of using propane and isobutane as working fluids and identifying MOFs with a superior performance.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating the Atomic Structures of the Gel Layer Formed during Aluminoborosilicate Glass Dissolution: An Integrated Experimental and Simulation Study

The altered glasses produced during aqueous dissolution of silicate and borosilicate glasses are among the most complex structures to understand at the atomic level due to their amorphous nature, random porosity and various levels of hydration. In this study, we gained insights of the complex atomic structures of altered aluminoborosilicate glasses by combining a range of experimental and computational approaches. The altered glasses were prepared by dissolution of three glasses with varying level of alumina in acid for 7 days. A comprehensive set of experimental [elemental analysis, high-energy X-ray diffraction, 29 Si and 27 Al solid-state nuclear magnetic resonance (NMR), and O 1s X-ray photoelectron spectroscopy (XPS)] and modeling (molecular dynamics (MD) simulations using non-reactive and reactive force fields) were used to study the atomic structures of these altered glasses. Elemental analysis showed that most of the B in the pristine glasses was leached into the solution and was not contained in the altered glass. The 29 Si and 27 Al solid-state NMR spectra revealed that the altered glasses have more polymerized silicate network as compared to those in the pristine glasses due to the reformation of linkages among Si and Al oxygen polyhedral in the altered glasses. The bridging and non-bridging (or hydroxyl O) atoms in the altered glasses were also quantified from their O 1s XPS spectra. Atomic structure models of the altered glasses were constructed by using MD simulations using the reactive force field based on the compositional information obtained from experiments. Various pore structures were generated using the charge-scaling (CS) method by using different initial densities and CS temperatures; the best CS parameters of each alternated glass were then determined by comparing with the experimental structure factors obtained from high energy X-ray diffraction. Pore and the atomic structures and vibrational properties around these pore surfaces were analyzed. Furthermore, these results from this comprehensive study thus provides a realistic insight into the pore morphology, atomic structure, and vibrational properties of altered glasses.

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

Preparation and Characterization of Inclusion Complexes of β-Cyclodextrin and Phenolics from Wheat Bran by Combination of Experimental and Computational Techniques

Bitterness often associated with whole wheat products may be related to phenolics in the bran. Cyclodextrins (CDs) are known to form inclusion complexes. The objective was to form inclusion complexes between β-CD and wheat phenolics. Pure phenolic acids (trans-ferulic acid (FA), caffeic acid (CA), and p-coumaric acid (CO)) and phenolic acids from wheat bran were used to investigate complex formation potential. Complexes were characterized by spectroscopy techniques, and a computational and molecular modeling study was carried out. The relative amount of complex formation between β-CD and wheat bran extract was CA > CO > FA. The phenolic compounds formed inclusion complexes with β-CDs by non-covalent bonds. The quantum-mechanical calculations supported the experimental results. The most stable complex was CO/β-CD complex. The ΔH value for CO/β-CD complex was -11.72 kcal/mol and was about 3 kcal/mol more stable than the other complexes. The QSPR model showed good correlation between binding energy and 1H NMR shift for the H5 signal. This research shows that phenolics and β-CD inclusion complexes could be utilized to improve the perception of whole meal food products since inclusion complexes have the potential to mask the bitter flavor and enhance the stability of the phenolics in wheat bran.

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

The influence of molecular design on structure–property relationships of a supramolecular polymer prodrug

Supramolecular self-assemblies of hydrophilic macromolecules functionalized with hydrophobic, structure-directing components have long been used for drug delivery. In these systems, loading of poorly soluble compounds is typically achieved through physical encapsulation during or after formation of the supramolecular assembly, resulting in low encapsulation efficiencies and limited control over release kinetics, which are predominately governed by diffusion and carrier degradation. To overcome these limitations, amphiphilic prodrugs that leverage a hydrophobic drug as both the therapeutic and structure-directing component can be used to create supramolecular materials with higher loading and controlled-release kinetics using biodegradable or enzymatically cleavable linkers. Here, we report the design, synthesis, and characterization of a library of supramolecular polymer prodrugs based on poly(ethylene glycol) (PEG) and the proregenerative drug 1,4-dihydrophenonthrolin-4-one-3-carboxylic acid (DPCA). Structure–property relationships were elucidated through experimental characterization of prodrug behavior in both the wet and dry states using scattering techniques and electron microscopy and corroborated by coarse-grained modeling. Molecular architecture and the hydrophobic-to-hydrophilic ratio of PEG–DPCA conjugates strongly influenced their physical state in water, ranging from fully soluble to supramolecular spherical assemblies and nanofibers. Molecular design and supramolecular structure, in turn, were shown to dramatically alter hydrolytic and enzymatic release and cellular transport of DPCA. In addition to potentially expanding therapeutic options for DPCA through control of supramolecular assemblies, the design principles elaborated here may inform the development of other supramolecular prodrugs based on hydrophobic small-molecule compounds.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Novel viruses of the family Partitiviridae discovered in Saccharomyces cerevisiae

It has been 49 years since the last discovery of a new virus family in the model yeast Saccharomyces cerevisiae . A large-scale screen to determine the diversity of double-stranded RNA (dsRNA) viruses in S . cerevisiae has identified multiple novel viruses from the family Partitiviridae that have been previously shown to infect plants, fungi, protozoans, and insects. Most S . cerevisiae partitiviruses (ScPVs) are associated with strains of yeasts isolated from coffee and cacao beans. The presence of partitiviruses was confirmed by sequencing the viral dsRNAs and purifying and visualizing isometric, non-enveloped viral particles. ScPVs have a typical bipartite genome encoding an RNA-dependent RNA polymerase (RdRP) and a coat protein (CP). Phylogenetic analysis of ScPVs identified three species of ScPV, which are most closely related to viruses of the genus Cryspovirus from the mammalian pathogenic protozoan Cryptosporidium parvum . Molecular modeling of the ScPV RdRP revealed a conserved tertiary structure and catalytic site organization when compared to the RdRPs of the Picornaviridae . The ScPV CP is the smallest so far identified in the Partitiviridae and has structural homology with the CP of other partitiviruses but likely lacks a protrusion domain that is a conspicuous feature of other partitivirus particles. ScPVs were stably maintained during laboratory growth and were successfully transferred to haploid progeny after sporulation, which provides future opportunities to study partitivirus-host interactions using the powerful genetic tools available for the model organism S . cerevisiae .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Bayesian optimization for chemical products and functional materials

The design of chemical-based products and functional materials is vital to modern technologies, yet remains expensive and slow. Artificial intelligence and machine learning offer new approaches to leverage data to overcome these challenges. Herein this review focuses on recent applications of Bayesian optimization (BO) to chemical products and materials including molecular design, drug discovery, molecular modeling, electrolyte design, and additive manufacturing. Numerous examples show how BO often requires an order of magnitude fewer experiments than Edisonian search. The essential equations for BO are introduced in a self-contained primer specifically written for chemical engineers and others new to the area. Finally, the review discusses four current research directions for BO and their relevance to product and materials design.

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

Atomistic Mechanics of Torn Back Folded Edges of Triangular Voids in Monolayer WS 2

Triangular nanovoids in 2D materials transition metal dichalcogenides have vertex points that cause stress concentration and lead to sharp crack propagation and failure. Here, the atomistic mechanics of back folding around triangular nanovoids in monolayer WS 2 sheets is examined. Combining atomic-resolution images from annular dark-field scanning transmission electron microscopy with reactive molecular modelling, it is revealed that the folding edge formation has statistical preferences under geometric conditions based on the orientation mismatch. It is further investigated how loading directions and strong interlayer friction, interplay with WS 2 lattice's crack preference, govern the deformation and fracture pattern around folding edges. Finally, these results provide fundamental insights into the combination of fracture and folding in flexible monolayer crystals and the resultant Moiré lattices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Polymorphism and Negative Linear Compressibility in Pyrazine-d4

Pyrazine (1,4-diazine) is a widely used linker in coordination polymers, molecular magnets, and metal–organic frameworks, yet its solid-state polymorphism is controlled by weak C–H···N interactions. Using variable-temperature and high-pressure neutron diffraction, a previously unrecognized low-temperature transition to phase IV is identified, showing that the form stabilized under pressure can also be accessed below ∼90 K at ambient pressure in the perdeuterated compound. The transition from phase III → IV involves a symmetry-breaking lattice distortion driven by the reorganization of C–H···N hydrogen bonds, accompanied by negative linear compressibility, which is primarily accommodated through the reduction in void-space within the structure. Structural analysis reveals that while the overall monoclinic framework of phase IV is preserved under both temperature and pressure routes, the evolution of hydrogen-bonded network is distinct depending on the pathway taken. These results establish pyrazine as a model molecular solid for studying pressure–temperature polymorphism and demonstrate how weak directional interactions can produce anomalous mechanical responses in organic crystals.

Ridley, Christopher J. [Oak Ridge National Laborat↗
DOE OSTI2023

Targeting Aspergillus allergen oryzin with a chemical probe at atomic precision

We report the molecular basis of Aspergillus fumigatus oryzin, allergen Asp f 13, or alkaline proteinase ALP1, containing the sequence motif His–Asp–Ser of the subtilisin family, structure, and function at atomic detail. Given the resolution of the data (1.06 Å), we use fragment molecular replacement with ideal polyalanine α-helices to determine the first crystal structure of oryzin. We probe the catalytic serine through formation of an irreversible bond to a small molecule compound, specifically labeling it, describing the amino acid residues performing the catalytic function. Defined by a self-processed pro-peptide, the active site architecture shapes up pocket-like subsites that bind to and unveil the S1′–S4′ substrate binding preferences. We use molecular modeling to dock a model of the pro-peptide in the S1–S4 region and to dock collagen along the active site cleft. Opposite to the face harboring the catalytic serine, the enzyme binds to a calcium ion in a binding site created by backbone flipping. We use thermal unfolding to show that this metal ion provides structural stability. With no known host inhibitor identified thus far, this structure may hasten the progress of developing new therapeutic agents for diseases caused by pathogenic fungi.

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

A structural underpinning of the lower critical solution temperature (LCST) behavior behind temperature-switchable liquids

In this work, we use state-of-the-art X-ray scattering and molecular dynamics to analyze amine-water mixtures that show the unusual lower critical solution temperature (LCST) behavior. The goal is to provide direct experimental evidence for the entropy-lowering molecular cluster formation hypothesized as necessary for LCST behavior. Differential wide-angle X-ray scattering and pair distribution analysis and small-angle X-ray scattering measurements were combined with molecular modeling and liquid-liquid equilibrium measurements, revealing direct experimental evidence for the hypothesis. However, the response of the amine phase to accommodating water is even more subtle than the simple hypothesis suggests, with the formation of robust nanoscale reverse micelles. The techniques developed in this paper can be expected to yield insights in the use of temperature-switchable liquids in solvent extraction and other separations, and the stabilization of organelles in living cells that do not have physical membranes but do require compositional gradients to operate.

36 MATERIALS SCIENCE↗