Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular conformation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Characterization of Heparin’s Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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

Polarized Resonant X-ray Scattering to Measure Molecular Orientation and Conformation in Organic Nanostructures

The ordering of organic (carbon-based) molecules into nanostructures drives useful materials properties for applications such as organic light emitting diodes, flexible batteries, and polymer nanocarriers. However, there is a current lack of techniques to measure this ordering. X-ray scattering techniques involving molecular resonances have shown promise in identifying molecular types, locations, and orientation. However, resonant soft X-ray scattering (RSoXS) techniques have not been developed even lacking basic optical models. The objectives achieved by this research were the development of methods, models, and analysis techniques that harness these resonant phenomena in X-ray scattering measurements to quantitatively probe ordering of organic molecules within devices and nanostructures. This award resulted in 11 peer reviewed publications with two others soon to be submitted. New physical models for X-ray scattering that accurately describe these complex molecular resonances enable measurements of how molecules assemble, orient and conform at the nanometer scale. These new techniques were applied to both organic electronic devices and aqueous nanoparticles being developed for hydrocarbon sequestration. This new measurement capability will enable characterization of structure-property relationships that lead to disruptive technologies of printable, stretchable, and biocompatible devices based on organic materials.

14 SOLAR ENERGY↗
NASA NTRS1990

Molecular arrangements and conformations of liquid n-tridecane chains confined between two hard walls

Monte Carlo simulations based on united CH2 atoms linked with fixed bond lengths and angles, and continuously varying torsional angles subject to appropriate interaction energies, were used to study liquid n-tridecane chains confined between two parallel hard walls with 5 mm separation. The molecules located in the central region are noted to be in very close agreement with ideal unperturbed chains. Nearly all the chain units in the first segmental layer adjacent to the walls belong to two-dimensional chains, and exhibit considerable orientational correlations between neighboring segments of different chains.

Vacatello, Michele↗
DOE OSTI2025

Functional protein mining with conformal guarantees

Molecular structure prediction and homology detection offer promising paths to discovering protein function and evolutionary relationships. However, current approaches lack statistical reliability assurances, limiting their practical utility for selecting proteins for further experimental and in-silico characterization. To address this challenge, we introduce a statistically principled approach to protein search leveraging principles from conformal prediction, offering a framework that ensures statistical guarantees with user-specified risk and provides calibrated probabilities (rather than raw ML scores) for any protein search model. Our method (1) lets users select many biologically-relevant loss metrics (i.e. false discovery rate) and assigns reliable functional probabilities for annotating genes of unknown function; (2) achieves state-of-the-art performance in enzyme classification without training new models; and (3) robustly and rapidly pre-filters proteins for computationally intensive structural alignment algorithms. Our framework enhances the reliability of protein homology detection and enables the discovery of uncharacterized proteins with likely desirable functional properties.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Human IgE monoclonal antibodies define two unusual epitopes trapping dog allergen Can f 1 in different conformations

Abstract Molecular analysis of interactions between IgE antibody and allergen allows the structural basis of IgE recognition to be defined. Human IgE (hIgE) epitopes of respiratory lipocalin allergens, including Can f 1, remain elusive due to a lack of IgE‐allergen complexes. This study aims to map the structure of allergenic epitopes on Can f 1. The fragment antigen‐binding (Fab) regions of Can f 1 specific human IgE monoclonal antibodies (hIgE mAb) were used to determine the structures of IgE epitopes. Epitope mutants were designed to target Can f 1 epitopes. Immunoassays and a human FcεRIαtransgenic mouse model of passive anaphylaxis in vivo were used to assess the functional activity of epitope mutants. Crystal structures of natural or recombinant Can f 1 complexed with two hIgE mAb 1J11 and 12F3 Fabs, respectively, were determined. The hIgE mAb bound to two partially overlapping epitopes and recognized two different Can f 1 conformations. The hIgE mAb 12F3 showed an unusual mode of binding by protruding its heavy chain CDR3 inside the Can f 1 calyx. Epitope mutants generated based on the structural analyses displayed a 64%–89% reduction in IgE antibody binding and failed to induce passive anaphylaxis in a human FcεRIαtransgenic mouse model. In summary, the structures of Can f 1‐hIgE Fab complexes revealed two unique and partially overlapping epitopes on Can f 1. The modification of the identified IgE epitopes provides a pathway for the design of hypoallergens to treat dog allergies.

Biochemistry & Molecular Biology↗
DOE OSTI2021

Amphiphile conformation impacts aggregate morphology and solution structure across multiple lengthscales

Although the self-assembly of amphiphiles is well-studied in aqueous solutions, much less is understood about the fundamental driving forces and structure property relationships in non-polar media. In recent work [Journal of Physical Chemistry B, 2020, 124, 10822.] the authors have studied a series of malonomide-based amphiphiles that are relevant to liquid-liquid extraction. That work demonstrated that aggregation is largely driven by local dipole-dipole interactions between molecules. Here, we build upon this observation to develop a more detailed understanding of how the balance of dipole-dipole interactions (controlled by conformation) and molecular architecture influences the morphology of the aggregates across lengthscales. Using constrained molecular dynamics about key degrees of freedom, we demonstrate that the conformation of N,N-dimethyl,N,N-dioctylhexylethoxy malonamide (DMDOHEMA) and N,N-dimethyl,N,N-dibutyltetradecyl malonamide (DMDBTDMA) has a significant impact upon self-association - where appropriate conformational sampling is essential. To quantify the aggregate morphology, several graph theoretic and persistent homology based properties are determined. The former examines the patterns of intermolecular interactions within clusters, while the latter examines the 3-dimensional spatial distribution. Based upon these analyses, we find that the morphology of aggregates, particularly at higher malonamide concentration, depends on a balance of dipole alignment and alkyl tail sterics. Dipole alignment encourages linear patterns of the intermolecular interactions within aggregates, while the alkyl tail steric interactions between the malonamide result in noticeably less linear aggregates for DMDOHEMA than DMDBTDMA. This is reflected in the spatial distribution, where more holes or voids exist between extractants within the DMDOHEMA that distribute within the solution in more of a "swiss cheese" arrangement as opposed to the more filamentous distribution of DMBDTDMA. This study links conformation and molecular structure to the morphology of amphiphile assemblies, and serves as a basis for ongoing study of multicomponent amphiphile solutions with polar and other solutes, and how these impact aggregation phenomena.

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

Multiscale Molecular Dynamics Simulations: Accelerating Conformational Sampling of Biomolecular Systems by Iterating All-Atom and Coarse-Grained Simulations

We developed the atomistic-coarse-grained multiscale MD simulation method in the OpenMM simulation package by iterating between the all-atom (AA) and coarse-grained (CG) MD simulations to enhance the sampling of biomolecular conformations. As the free energy surfaces are flattened during CG MD simulations, we can accelerate the transitions between different low-energy conformations. The AA-CG-AA cycles are repeated, facilitating the accelerated sampling of biomolecular conformations at a CG level, while the finer atomistic interactions are refined with AA simulators.

Do, Hung Nguyen↗
DOE OSTI2025

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗
DOE OSTI2025

Depletion of the Protein Hydration Shell with Increasing Temperature Observed by Small-Angle X-ray Scattering and Molecular Simulations

The hydration shell is an integral part of proteins since it plays key roles in conformational transitions, molecular recognition, and enzymatic activity. While the dynamics of the hydration shell have been described by spectroscopic techniques, the structure of the hydration shell remains less understood due to the lack of hydration shell-sensitive structural probes with high spatial resolution. We combined temperature-ramp small-angle X-ray scattering (T-ramp SAXS) from 255 to 335 K with molecular simulations to demonstrate that the hydration shells of the IgG-binding domain of Protein G (GB3) and the villin headpiece are remarkably temperature-sensitive. For proteins in the folded state, T-ramp SAXS data and explicit-solvent SAXS predictions consistently demonstrate decays of protein contrasts and radii of gyration with increasing temperature, which are shown to reflect predominantly temperature-sensitive, depleting hydration shells. The depletion is caused not merely by enhanced disorder within the hydration shells but also by partial displacements of surface-coordinated water molecules. Together, T-ramp SAXS and explicit-solvent SAXS calculations provide a novel structural view of the protein hydration shell, which underlies temperature-dependent processes such as cold denaturation, thermophoresis, or biomolecular phase separation.

electron density↗
DOE OSTI2022

Adaptive Conformer Sampling for Property Prediction Using the Conductor-like Screening Model for Real Solvents

The valorization of lignocellulose-derived bioproducts requires effective separation from excessive water. Liquid–liquid extraction is a promising low-energy separation technology, but effective extraction requires solvent selection based on the thermodynamic properties of the bioproduct and solvent components. We propose a computational framework for predicting such properties by developing an adaptive conformer selection approach for use with COSMO-RS (conductor-like screening model for real solvents) calculations. In this framework, molecular dynamics simulations are used to generate many molecular structures (conformers) at representative temperatures in varying solvent environments. Conformers are then clustered based on structural metrics in a low-dimensional space and selected using a mixed-integer quadratic programming problem to iteratively insert a sampled conformer. At each iteration, we determine bioproduct properties using COSMO-RS. Here, we demonstrate the capability of the proposed framework on representative bioproducts to show convergence of the adaptive sampling toward experimentally measured properties with fewer calculations than required by random conformer sampling, enabling the improved screening of solvent systems for liquid-phase separation.

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

Iron catalyst chemistry in modeling a high-pressure carbon monoxide nanotube reactor

The high-pressure carbon monoxide (HiPco) technique for producing single-wall carbon nanotubes (SWNTs) is analyzed with the use of a chemical reaction model coupled with flow properties calculated along streamlines, calculated by the FLUENT code for pure carbon monoxide. Cold iron pentacarbonyl, diluted in CO at about 30 atmospheres, is injected into a conical mixing zone, where hot CO is also introduced via three jets at 30 degrees with respect to the axis. Hot CO decomposes the Fe(CO)5 to release atomic Fe. Then iron nucleates and forms clusters that catalyze the formation of SWNTs by a disproportionation reaction (Boudouard) of CO on Fe-containing clusters. Alternative nucleation rates are estimated from the theory of hard sphere collision dynamics with an activation energy barrier. The rate coefficient for carbon nanotube growth is estimated from activation energies in the literature. The calculated growth was found be about an order of magnitude greater than measured, regardless of the nucleation rate. A study of cluster formation in an incubation zone prior to injection into the reactor shows that direct dimer formation from Fe atoms is not as important as formation via an exchange reaction of Fe with CO in FeCO.

Evaluation Studies↗
DOE OSTI2022

Bioimaging Science Program: 2022 Principal Investigator Meetings Proceedings

The mission of the U.S. Department of Energy’s (DOE) Biological and Environmental Research (BER) program’s Bioimaging Science Program (BSP) is to understand the translation of genomic information into the mechanisms that power living cells, communities of cells, and whole organisms. The goal of BSP is to develop new imaging and measurement technologies to visualize the spatial and temporal relationships of key metabolic processes governing phenotypic expression in plants and microbes. The extended goal of dynamic imaging is to functionally connect cellular components and interdependent organisms. Information about the time and place of chemical reactions in situ can identify causal relationships between biological activators and downstream effectors. BSP held its annual PI meeting virtually February 28–March 1. Contributing investigators are convened to review progress and current state-of-the-art bioimaging research. Holding the 2022 BSP meeting as part of the broader Genomic Science Program (GSP) PI meeting allowed researchers to interact with the extended GSP community. This convergence provided a platform for networking and exchange of ideas with experts in other technologies and in target BSP application areas, helping to forge new multidisciplinary collaborations among investigators from the sister programmatic areas within BER’s Biological Systems Science Division. An important highlight of the BSP meeting was the keynote presentation by Nobel Laureate Dr. Joachim Frank on Time-Resolved Macromolecular Imaging using Cryo-EM. He discussed microfluidic mixing and fast freezing to capture nonequilibrium intermediate states during molecular binding and conformational changes. The action of molecular machines can be captured at nanometer resolution and millisecond discrimination. BSP PIs made presentations describing their research focus and progress in plenary sessions on bioimaging science and on quantum-enabled bioimaging science research projects. BSP research at universities and DOE laboratories is presented in this report. A final discussion of the BSP was organized by meeting plenary session chairs, who prepared the following Executive Summary of current BSP research, research challenges, future opportunities, and potential ideas for expanding the BSP’s impact and interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1991

Ab initio study of the electrostatic multipole nature of torsional potentials in CH3SSCH3, CH3SSH, and HOOH

The origin of torsional potentials in H3CSSCH3, H3CSSH, and HOOH and the anisotropy of the local charge distribution has been analyzed in terms of atomic multipoles calculated from the ab initio LCAO-MO-SCF wave function in the 6-31G* basis set. The results indicate that for longer -S-S-bonds the major contribution to these torsional barriers are electrostatic interactions of the atomic multipoles located on two atoms forming the rotated bond. This finding demonstrates the important role of electrostatic 1-2 interatomic interactions, usually neglected in conformational studies. It also opens the possibility to derive directly from accurate ab initio wave functions a simple nonempirical torsional potential involving atomic multipoles of two bonded atoms defining the torsional angle. For shorter -O-O- bonds, use of more precise models and inclusion of 1-3 interactions seems to be necessary.

NASA Program Exobiology↗
DOE OSTI2021

Cryo-EM structures of full-length Tetrahymena ribozyme at 3.1 Å resolution

Single-particle cryogenic electron microscopy (cryo-EM) has become a standard technique for determining protein structures at atomic resolution1,2,3. However, cryo-EM studies of protein-free RNA are in their early days. The Tetrahymena thermophila group I self-splicing intron was the first ribozyme to be discovered and has been a prominent model system for the study of RNA catalysis and structure–function relationships4, but its full structure remains unknown. Here we report cryo-EM structures of the full-length Tetrahymena ribozyme in substrate-free and bound states at a resolution of 3.1 Å. Additionally, newly resolved peripheral regions form two coaxially stacked helices; these are interconnected by two kissing loop pseudoknots that wrap around the catalytic core and include two previously unforeseen (to our knowledge) tertiary interactions. The global architecture is nearly identical in both states; only the internal guide sequence and guanosine binding site undergo a large conformational change and a localized shift, respectively, upon binding of RNA substrates. These results provide a long-sought structural view of a paradigmatic RNA enzyme and signal a new era for the cryo-EM-based study of structure–function relationships in ribozymes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A unifying Bayesian framework for merging X-ray diffraction data

Novel X-ray methods are transforming the study of the functional dynamics of biomolecules. Key to this revolution is detection of often subtle conformational changes from diffraction data. Diffraction data contain patterns of bright spots known as reflections. To compute the electron density of a molecule, the intensity of each reflection must be estimated, and redundant observations reduced to consensus intensities. Systematic effects, however, lead to the measurement of equivalent reflections on different scales, corrupting observation of changes in electron density. Here, we present a modern Bayesian solution to this problem, which uses deep learning and variational inference to simultaneously rescale and merge reflection observations. We successfully apply this method to monochromatic and polychromatic single-crystal diffraction data, as well as serial femtosecond crystallography data. We find that this approach is applicable to the analysis of many types of diffraction experiments, while accurately and sensitively detecting subtle dynamics and anomalous scattering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

LoTToR: An Algorithm for Missing-Wedge Correction of the Low-Tilt Tomographic 3D Reconstruction of a Single-Molecule Structure

A single-molecule three-dimensional (3D) structure is essential for understanding the thermal vibrations and dynamics as well as the conformational changes during the chemical reaction of macromolecules. Individual-particle electron tomography (IPET) is an approach for obtaining a snap-shot 3D structure of an individual macromolecule particle by aligning the tilt series of electron tomographic (ET) images of a targeted particle through a focused iterative 3D reconstruction method. The method can reduce the influence on the 3D reconstruction from large-scale image distortion and deformation. Due to the mechanical tilt limitation, 3D reconstruction often contains missing-wedge artifacts, presented as elongation and an anisotropic resolution. Here, we report a post-processing method to correct the missing-wedge artifact. This low-tilt tomographic reconstruction (LoTToR) method contains a model-free iteration process under a set of constraints in real and reciprocal spaces. A proof of concept is conducted by using the LoTToR on a phantom, i.e., a simulated 3D reconstruction from a low-tilt series of images, including that within a tilt range of ±15°. The method is validated by using both negative-staining (NS) and cryo-electron tomography (cryo-ET) experimental data. A significantly reduced missing-wedge artifact verifies the capability of LoTToR, suggesting a new tool to support the future study of macromolecular dynamics, fluctuation and chemical activity from the viewpoint of single-molecule 3D structure determination.

97 MATHEMATICS AND COMPUTING↗