Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conformation dynamics”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Molecular mechanism of biased signaling at the kappa opioid receptor

The κ-opioid receptor (KOR) has emerged as an attractive drug target for pain management without addiction, and biased signaling through particular pathways of KOR may be key to maintaining this benefit while minimizing side-effect liabilities. As for most G protein-coupled receptors (GPCRs), however, the molecular mechanisms of ligand-specific signaling at KOR have remained unclear. To better understand the molecular determinants of KOR signaling bias, we apply structure determination, atomic-level molecular dynamics (MD) simulations, and functional assays. We determine a crystal structure of KOR bound to the G protein-biased agonist nalfurafine, the first approved KOR-targeting drug. We also identify an arrestin-biased KOR agonist, WMS-X600. Using MD simulations of KOR bound to nalfurafine, WMS-X600, and a balanced agonist U50,488, we identify three active-state receptor conformations, including one that appears to favor arrestin signaling over G protein signaling and another that appears to favor G protein signaling over arrestin signaling. These results, combined with mutagenesis validation, provide a molecular explanation of how agonists achieve biased signaling at KOR.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Polymer principles behind solubilizing lignin with organic cosolvents for bioenergy

Lignin solubilization is key to a viable biorefinery because its removal leads to facile deconstruction of biomass and because the isolated lignin can serve to derive precursors of novel high-value materials. The mixing of organic solvents with water has been shown to improve biomass fractionation and lignin conversion reactions. However, generally-applicable solubilization strategies are lacking because of the remarkable variability of lignin across plant feedstocks. Here, to obtain a predictive understanding of lignin solvation, we perform molecular dynamics simulations of model lignin polymers in two mixtures of water with polar aprotic solvents: tetrahydrofuran (THF) : water and γ-valerolactone (GVL) : water. The model lignins include H-, G- and S-only homopolymers and lignins with an S : G 1 : 1 ratio. We find that a well-established theory of self-avoiding polymers in a “good” solvent describes accurately the physical conformations of all types of lignin in both solvents. As the degree of methoxy substitution increases in the homopolymers, the distributions of the lignin radius of gyration and the Flory exponent ν , which describes the lignin-solvent interactions, do not change in THF : water, while ν shifts to slightly higher values in GVL : water. We attribute this increase to the interaction between the methyl group of GVL with the lignin methoxy groups. We also find that the reduction in the lignin radius of gyration due to branching is accurately described by the Zimm–Stockmayer theory for both THF : water and GVL : water. The above findings validate the applicability of polymer physics concepts to lignin and suggest that GVL : water may have the most favorable interaction with S-lignin, whereas the interactions of THF : water with lignin are independent of lignin monomeric content.

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

Three-point energy correlators in hadronic Higgs boson decays

We present the analytic calculation of the leading-order three-point energy correlator (EEEC) in hadronic Higgs decays, including both the gluon-initiated channel H → g g + X and quark-initiated channel H → q q ¯ + X . The phase-space integration is evaluated directly using Mandelstam variables s i j = ( p i + p j ) 2 , and the appearing square roots can be rationalized by either conformal ratios or celestial coordinate variables. Throughout the calculation, we observe the same transcendental function space as in N = 4 super–Yang-Mills theory and e + e − → hadrons . Different infrared limits are also explored using the full analytic result, offering the fixed-order data for EEEC factorization and resummation. Given its nontrivial shape dependence, the EEEC presents an excellent opportunity to explore the dynamics of gluon jets originating from the H → g g decay channel at future lepton colliders. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2022

A simple technique to classify diffraction data from dynamic proteins according to individual polymorphs

One often observes small but measurable differences in the diffraction data measured from different crystals of a single protein. These differences might reflect structural differences in the protein and may reveal the natural dynamism of the molecule in solution. Partitioning these mixed-state data into single-state clusters is a critical step that could extract information about the dynamic behavior of proteins from hundreds or thousands of single-crystal data sets. Mixed-state data can be obtained deliberately (through intentional perturbation) or inadvertently (while attempting to measure highly redundant single-crystal data). To the extent that different states adopt different molecular structures, one expects to observe differences in the crystals; each of the polystates will create a polymorph of the crystals. After mixed-state diffraction data have been measured, deliberately or inadvertently, the challenge is to sort the data into clusters that may represent relevant biological polystates. Here, this problem is addressed using a simple multi-factor clustering approach that classifies each data set using independent observables, thereby assigning each data set to the correct location in conformational space. This procedure is illustrated using two independent observables, unit-cell parameters and intensities, to cluster mixed-state data from chymotrypsinogen (ChTg) crystals. It is observed that the data populate an arc of the reaction trajectory as ChTg is converted into chymotrypsin.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An Implementation of Replica Exchange with Dynamical Scaling for Efficient Large-Scale Simulations

An implementation of the replica exchange with dynamical scaling (REDS) method in the commonly used molecular dynamics program GROMACS is presented. REDS is a replica exchange method that requires fewer replicas than conventional replica exchange while still providing data over a range of temperatures and can be used in either constant volume or constant pressure ensembles. Details for running REDS simulations are given, and an application to the human islet amyloid polypeptide (hIAPP) 11-25 fragment shows that the model efficiently samples conformational space.

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

Spectral Element Applications in Complex Nuclear Reactor Geometries: Tet-to-Hex Meshing

The spectral element code Nek5000 is an open-source, higher-order computational fluid dynamics code developed at Argonne National Laboratory. It is designed to solve incompressible Navier-Stokes equations, but it also has a low-Mach-number approximation feature available. Large eddy simulation is approached by explicit filtering of the velocity field (and other fields) to mimic the effect of dissipation due to the unresolved scale. The computational domain is decomposed into second-order hexahedral elements that conform to the boundaries. However, generating a high-quality pure-hexahedral mesh can be challenging for some problems. For simple geometries, traditional blocking methods can be used to decompose the domain into smaller blocks to generate a so-called structural mesh. A structural mesh can maintain good orthogonality but can have a highly skewed mesh to conform to the geometry, as well as unnecessary refinement in the far field. Moreover, for geometries with relative complexity, blocking the geometry becomes impossible. To address these issues, we adopted a let-to-hex strategy to generate a pure hexahedral mesh for Nek5000. First, we generate a pure tetrahedral mesh for the geometry; then we divide one tetrahedral element into four hexahedral elements. A pure tetrahedral mesh could be easily generated for complex geometries by using many current meshing codes. In this paper, we use the commercial codes ANSYS meshing and ANSYSICEM to generate the pure tetrahedral mesh and then convert it to a pure hexahedral mesh. Boundary layers are extruded in ANSYSICEM to maintain near-wall resolution.

ANSYS ICEM↗
DOE OSTI2022

Binding Interactions between Receptor-Binding Domain of Spike Protein and Human Angiotensin Converting Enzyme-2 in Omicron Variant

The emergence of new SARS-CoV-2 Omicron variant of concern (OV) has exacerbated the COVID-19 pandemic because of a large number of mutations in the spike protein, particularly in the receptor-binding domain (RBD), resulting in highly contagious and/or vaccine-resistant strains. Herein, we present a systematic analysis based on detailed molecular dynamics (MD) simulations in order to understand how the OV RBD mutations affect the ACE2 binding. We show that the OV RBD binds to ACE2 more efficiently and tightly predominantly because of strong electrostatic interactions, thereby promoting increased infectivity and transmissibility compared to other strains. Some of the OV RBD mutations are predicted to affect the antibody neutralization either through their role in the S-protein conformational changes, such as S371L, S373P, and S375F, or through changing its surface charge distribution, such as G339D, N440K, T478K, and E484A. Other mutations, such as K417N, G446S, and Y505H, decrease the ACE2 binding, whereas S447N, Q493R, G496S, Q498R, and N501Y tend to increase it.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

On the Rapid Calculation of Binding Affinities for Antigen and Antibody Design and Affinity Maturation Simulations

The accurate and efficient calculation of protein-protein binding affinities is an essential component in antibody and antigen design and optimization, and in computer modeling of antibody affinity maturation. Such calculations remain challenging despite advances in computer hardware and algorithms, primarily because proteins are flexible molecules, and thus, require explicit or implicit incorporation of multiple conformational states into the computational procedure. The astronomical size of the amino acid sequence space further compounds the challenge by requiring predictions to be computed within a short time so that many sequence variants can be tested. In this study, we compare three classes of methods for antibody/antigen (Ab/Ag) binding affinity calculations: (i) a method that relies on the physical separation of the Ab/Ag complex in equilibrium molecular dynamics (MD) simulations, (ii) a collection of 18 scoring functions that act on an ensemble of structures created using homology modeling software, and (iii) methods based on the molecular mechanics-generalized Born surface area (MM-GBSA) energy decomposition, in which the individual contributions of the energy terms are scaled to optimize agreement with the experiment. When applied to a set of 49 antibody mutations in two Ab/HIV gp120 complexes, all of the methods are found to have modest accuracy, with the highest Pearson correlations reaching about 0.6. In particular, the most computationally intensive method, i.e., MD simulation, did not outperform several scoring functions. The optimized energy decomposition methods provided marginally higher accuracy, but at the expense of requiring experimental data for parametrization. Within each method class, we examined the effect of the number of independent computational replicates, i.e., modeled structures or reinitialized MD simulations, on the prediction accuracy. We suggest using about ten modeled structures for scoring methods, and about five simulation replicates for MD simulations as a rule of thumb for obtaining reasonable convergence. We anticipate that our study will be a useful resource for practitioners working to incorporate binding affinity calculations within their protein design and optimization process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Challenges of conventional iterative all-atom and coarse-grained multiscale molecular dynamics

In this work, we evaluate the biomolecular dynamics behaviors when conventionally iterating between all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations over multiple cycles. We implemented the workflow to iterate between AA and CG in OpenMM, namely the iterative multiscale MD (iMMD) simulation workflow. In particular, we aim to identify practical applications for iterating between AA and CG simulations in a conventional manner without any constraints or model modifications. We evaluate the iMMD workflow on four representative systems, spanning folding of two soluble proteins and protein-protein as well as protein-lipid interactions of two membrane proteins. We observe that iteration between AA and CG representations could help the soluble proteins exit undesirable metastable states to fold, resulting from random protein structural distortions due to cycling. Consequently, the most reliable use of iterative AA and CG simulations appears to be to accelerating complex lipid mixing for membrane-bound protein systems rather than sampling protein conformational space. Our work explores the practical usages and limitations for iterative AA and CG simulations using readily available AA and CG force fields. The evaluated iMMD workflow in OpenMM is made available at https://github.com/lanl/iMMD.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Structure and dynamics of small polyimide oligomers with silicon as a function of aging

The effect of UV curing and shearing on the structure and behavior of a polyimide (PI) binder as it disperses silicon particles in a battery electrode slurry was investigated. PI dispersant effectiveness increases with UV curing time, which controls the overall binder molecular weight. Here, the shear force during electrode casting causes higher molecular weight PI to agglomerate, resulting in battery anodes with poorly dispersed Si particles that do not cycle well. It is hypothesized that when PI binder is added above a critical amount, it conformally coats the silicon particles and greatly impedes Li ion transport. There is an “interzonal region” for binder loading where it disperses silicon well and provides a coverage that facilitates Li transport through the anode material and into the silicon particles. These results have implications in ensuring reproducible electrode manufacturing and increasing cell performance by optimizing the PI structure and coordination with the silicon precursor.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interfacial Response and Structural Adaptation of Structured Polyelectrolyte Thin Films

Ionizable block copolymers with distinctive block characteristics display the diversity crucial for the design of macromolecules for targeted applications. In contrast to van der Waals copolymers, their interfaces, which are critical to their function, consist of nanodomains, each of a different nature and thus unique interfacial behavior. Here, the interfacial response of a symmetric block copolymer with a sulfonated polystyrene polyelectrolyte center, tethered to polyethylene-r-propylene and terminated by poly(t-butyl styrene) is probed as polymer films are exposed to three polar solvents, water, propanol, and tetrahydrofuran (THF), using molecular dynamics simulations. Each of the solvents captures a distinctive interaction with the individual blocks. We find that at the film boundary, the interfacial response is initially dominated by that of the hydrophobic blocks to all solvents. At later times, the solvent distribution among the blocks, where water molecules associate predominantly with the sulfonated groups and propanol and THF reside at multiple different sites, determines the chemical composition and the polymer conformation at the interface. Altogether, these simulations provide the first direct molecular insight into the interfacial response of ionizable copolymers.

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

Excited-state vibronic coherences and intramolecular charge transfer dynamics of the photoinactive cyanobacteriochrome NpF2164g5

We have characterized the excited-state structural dynamics that follow optical excitation of the photoinactive cyanobacteriochrome (CBCR) NpF2164g5 to determine the first events along the photoisomerization reaction coordinate of the phycocyanobilin (PCB) chromophore in red-to-green photoactive CBCRs. Within 25 fs of photoexcitation to the first excited singlet state, S1, a cross peak begins to develop below the diagonal of the broadband two-dimensional electronic spectrum (2DES) owing to the formation of a twisted conformation of the PCB chromophore with an enhanced intramolecular charge-transfer (ICT) character. Excited-state coherent wavepacket motions, including a torsion of the methine bridge between rings C and D and the carbon–hydrogen out-of-plane (HOOP) wagging vibration, are rapidly damped as the cross peak forms. This finding supports an assignment of the torsional and HOOP modes to the reaction coordinate of a coherent nonadiabatic mechanism. The ICT process is accompanied by an ultrafast Stokes shift and rapidly damped oscillations at the torsional mode’s frequency, which are sensed by measurements of the energy gap of the cross peak using its first moment with respect to the detection frequency axis of the 2DES spectrum. Furthermore, these results have the further implication that polar side chains of nearby amino acid residues in the binding site of the PCB chromophore can manipulate the barrier height for photoisomerization in red-to-green photoactive CBCRs and in phytochromes by redistributing the π-electron density along the methine bridge between the C and D rings.

Charge transfer↗
DOE OSTI2021

The solution structures of higher-order human telomere G-quadruplex multimers

Human telomeres contain the repeat DNA sequence 5'-d(TTAGGG), with duplex regions that are several kilobases long terminating in a 3' single-stranded overhang. The structure of the single-stranded overhang is not known with certainty, with disparate models proposed in the literature. We report here the results of an integrated structural biology approach that combines small-angle X-ray scattering, circular dichroism (CD), analytical ultracentrifugation, size-exclusion column chromatography and molecular dynamics simulations that provide the most detailed characterization to date of the structure of the telomeric overhang. We find that the single-stranded sequences 5'-d(TTAGGG)n, with n = 8, 12 and 16, fold into multimeric structures containing the maximal number (2, 3 and 4, respectively) of contiguous G4 units with no long gaps between units. The G4 units are a mixture of hybrid-1 and hybrid-2 conformers. In the multimeric structures, G4 units interact, at least transiently, at the interfaces between units to produce distinctive CD signatures. Global fitting of our hydrodynamic and scattering data to a worm-like chain (WLC) model indicates that these multimeric G4 structures are semi-flexible, with a persistence length of ~34 Å. Investigations of its flexibility using MD simulations reveal stacking, unstacking, and coiling movements, which yield unique sites for drug targeting.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Unraveling Morphology and Chemistry Dynamics in Fluoroethylene Carbonate Generated Silicon Anode Solid Electrolyte Interphase Across Delithiated and Lithiated States: Relative Cycling Stability Enabled by an Elastomeric Polymer Matrix

The silicon solid electrolyte interphase (SEI) faces cyclical cracking and reconstruction due to the ∼350% volume expansion. Understanding the SEI dynamic morphology and chemistry evolution from delithiated to lithiated states is thereby paramount to engineering a stable Si anode. Fluoroethylene carbonate (FEC) is a preferred additive with widely demonstrated enhancement of the Si cycling. Thus, insights into the dynamics of the FEC-SEI may provide hints toward engineering the Si interface. Herein, complementary ATR-FTIR, AFM, tip IR, and XPS probing reveal the presence of an elastomeric polycarbonate-like matrix in the FEC-generated SEI which is absent from the FEC-free SEI. Adding FEC to the baseline 1 M LiPF 6 in EC:EMC (1:1) electrolyte promotes formation of a thinner and more conformal SEI, and subdues morphology and chemistry changes between consecutive half-cycles. From AFM, morphological stabilization of the FEC-SEI occurs earlier. Furthermore, conventional SEI biproducts such as Li 2 CO 3 and LiEDC appear in reduced quantities in the FEC-SEI implying a reduced quantity of Li-consuming species. The thin polymeric FEC-SEI enables deeper (de)lithiation of silicon. In conclusion, the enhanced mechanical compliance, chemical invariance, and reduced Li inventory consumption of the FEC-SEI are logically the key features underlying the Si cycling enhancement by FEC.

25 ENERGY STORAGE↗
DOE OSTI2025

Dynamic Interfacial Architectures: Cruciferin‐Stabilized Oil/Water Interfaces for Sustainable Emulsions

Stabilizing oil-water interfaces in emulsions by plant-based proteins provides sustainable and tunable ways for designing emulsions with specific properties, for food, healthcare, and pharmaceuticals. Cruciferin, a protein from rapeseed, has great potential as green emulsifier, but details about its structure and mobility at oil-water interfaces are largely unknown. Here, these properties are studied with small angle neutron and x-ray scattering, and neutron spin echo spectroscopy, analyzed by atomistic modelling of scattering curves and coarse-grained modelling, to gain insight into interface coverage, and molecular conformation and mobility at the interface. Cruciferin assumes trimeric conformations at the interface, as in solution, but with its protrusions from the central core of the subunits (“arms”) more compressed. Interfacial mobility is only marginally lower than in solution, indicating the arms still transiently extend and preserve a network, for the first time revealing the mechanism how cruciferin forms highly elastic 2d gel-like oil-water interfaces, as observed in macroscopic rheology. The high interfacial mobility may help in self-repairing non-stabilized interfacial fractions, reducing coalescence. These findings provide a deeper molecular level understanding of proteins at oil-water interfaces, which can stimulate development of new plant-based emulsion products, and contribute to the global protein transition.

cruciferin↗
DOE OSTI2025

Photoactivation Transition State and Dynamical Response of the Orange Carotenoid Protein

The orange carotenoid protein (OCP) regulates light harvesting in cyanobacteria by acting as a photoreceptor in its resting form, OCP O , and by effecting the quenching of bilin excitons upon binding to the core of the phycobilisome in its photoactivated red form, OCP R . We show herein using fluorescence anisotropy measurements and the action spectra for the rate constants of the two light-driven steps in the mechanism that the photoactivation of the OCP from Synechocystis sp. PCC 6803 is triggered by excited-state motions of the canthaxanthin chromophore that yield a twisted and bent conformation. Well-tempered metadynamics simulations reveal that a bicycle-pedal configuration, due to twisting of the two adjacent C=C bonds at the C13−C14 and C15−C15′ positions in the center of canthaxanthin’s π-conjugated isoprenoid backbone, can be accommodated by the binding site in the OCP, with the energy of a captured photon required to cross the local activation energy barriers from the dark equilibrium structure. The bicycle-pedal configuration breaks the conserved hydrogen-bonding interactions between the carbonyl substituent of the β-ionone end ring of canthaxanthin and the adjacent W288 and Y201 residues in the C-terminal domain. The action spectra are modulated by the vibronic excitation prepared by absorption transitions to the S 2 state, indicating that the photoactivation reactions are triggered by the canthaxanthin chromophore well prior to vibrational equilibration. These findings show that an ultrafast structural response of the OCP protein to the excited-state motions of the canthaxanthin chromophore controls the photoactivation yield and the sensing of blue light.

Fluorescence↗