Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conformational ensembles”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Extended experimental inferential structure determination method in determining the structural ensembles of disordered protein states

Proteins with intrinsic or unfolded state disorder comprise a new frontier in structural biology, requiring the characterization of diverse and dynamic structural ensembles. Here we introduce a comprehensive Bayesian framework, the Extended Experimental Inferential Structure Determination (X-EISD) method, which calculates the maximum log-likelihood of a disordered protein ensemble. X-EISD accounts for the uncertainties of a range of experimental data and back-calculation models from structures, including NMR chemical shifts, J-couplings, Nuclear Overhauser Effects (NOEs), paramagnetic relaxation enhancements (PREs), residual dipolar couplings (RDCs), hydrodynamic radii (R h ), single molecule fluorescence Förster resonance energy transfer (smFRET) and small angle X-ray scattering (SAXS). We apply X-EISD to the joint optimization against experimental data for the unfolded drkN SH3 domain and find that combining a local data type, such as chemical shifts or J-couplings, paired with long-ranged restraints such as NOEs, PREs or smFRET, yields structural ensembles in good agreement with all other data types if combined with representative IDP conformers.

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

Exploring fingerprints of ultrafast structural dynamics in molecular solutions with an X-ray laser

We apply ultrashort x-ray laser pulses to track optically excited structural dynamics of [Ir 2 (dimen) 4 ] 2+ molecules in solution. In our exploratory study we determine angular correlations in the scattered x-rays, which comprise a complex fingerprint of the ultrafast dynamics. Model-assisted analysis of the experimental correlation data allows us to elucidate various aspects of the photoinduced changes in the excited molecular ensembles. We unambiguously identify that in our experiment the photoinduced transition dipole moments in [Ir 2 (dimen) 4 ] 2+ molecules are oriented perpendicular to the Ir-Ir bond. The analysis also shows that the ground state conformer of [Ir 2 (dimen) 4 ] 2+ with a larger Ir–Ir distance is mostly responsible for the formation of the excited state. We also reveal that the ensemble of solute molecules can be characterized with a substantial structural heterogeneity due to solvent influence. In conclusion, the proposed x-ray correlation approach offers an alternative path for studies of ultrafast structural dynamics of molecular ensembles in the liquid and gas phases.

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

Calabi-Yau CFTs and random matrices

Using numerical methods for finding Ricci-flat metrics, we explore the spectrum of local operators in two-dimensional conformal field theories defined by sigma models on Calabi-Yau targets at large volume. Focusing on the examples of K3 and the quintic, we show that the spectrum, averaged over a region in complex structure moduli space, possesses the same statistical properties as the Gaussian orthogonal ensemble of random matrix theory.

superstring vacua↗
DOE OSTI2022

Robust and scalable uncertainty estimation with conformal prediction for machine-learned interatomic potentials

Abstract Uncertainty quantification (UQ) is important to machine learning (ML) force fields to assess the level of confidence during prediction, as ML models are not inherently physical and can therefore yield catastrophically incorrect predictions. Established a-posteriori UQ methods, including ensemble methods, the dropout method, the delta method, and various heuristic distance metrics, have limitations such as being computationally challenging for large models due to model re-training. In addition, the uncertainty estimates are often not rigorously calibrated. In this work, we propose combining the distribution-free UQ method, known as conformal prediction (CP), with the distances in the neural network’s latent space to estimate the uncertainty of energies predicted by neural network force fields. We evaluate this method (CP+latent) along with other UQ methods on two essential aspects, calibration, and sharpness, and find this method to be both calibrated and sharp under the assumption of independent and identically-distributed (i.i.d.) data. We show that the method is relatively insensitive to hyperparameters selected, and test the limitations of the method when the i.i.d. assumption is violated. Finally, we demonstrate that this method can be readily applied to trained neural network force fields with traditional and graph neural network architectures to obtain estimates of uncertainty with low computational costs on a training dataset of 1 million images to showcase its scalability and portability. Incorporating the CP method with latent distances offers a calibrated, sharp and efficient strategy to estimate the uncertainty of neural network force fields. In addition, the CP approach can also function as a promising strategy for calibrating uncertainty estimated by other approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

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

Discreteness and integrality in Conformal Field Theory

Various observables in compact CFTs are required to obey positivity, discreteness, and integrality. Positivity forms the crux of the conformal bootstrap, but understanding of the abstract implications of discreteness and integrality for the space of CFTs is lacking. We systematically study these constraints in two-dimensional, non-holomorphic CFTs, making use of two main mathematical results. First, we prove a theorem constraining the behavior near the cusp of integral, vector-valued modular functions. Second, we explicitly construct non-factorizable, non-holomorphic cuspidal functions satisfying discreteness and integrality, and prove the non-existence of such functions once positivity is added. Application of these results yields several bootstrap-type bounds on OPE data of both rational and irrational CFTs, including some powerful bounds for theories with conformal manifolds, as well as insights into questions of spectral determinacy. We prove that in rational CFT, the spectrum of operator twists t ≥ c/12 is uniquely determined by its complement. Likewise, we argue that in generic CFTs, the spectrum of operator dimensions Δ > c–1/12 is uniquely determined by its complement, absent fine-tuning in a sense we articulate. Finally, we discuss implications for black hole physics and the (non-)uniqueness of a possible ensemble interpretation of AdS 3 gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Investigating SU(3) with Nf=8 fundamental fermions at strong renormalized coupling

Lattice simulations have observed a novel strong coupling symmetric mass generation (SMG) phase for the SU(3) gauge system with $N_f=8$ fundamental fermions (represented by two sets of staggered fields) at very large renormalized coupling ($g^2_{GF} \gtrsim 25$). The results of Phys.Rev.D 106 (2022) 014513 suggest that the SMG phase is separated from the weak coupling, conformal phase by a continuous phase transition, implying that the SMG phase exists in the continuum limit. To scrutinize these findings, we are generating a set of large volume zero temperature ensembles using nHYP improved staggered fermions with additional Pauli-Villars fields to tame gauge field fluctuations. We consider the low-lying meson spectrum and verify the existence of the SMG phase. Based on a finite size scaling analysis we predict that the phase transition between the strong and weak coupling phases is likely governed by a merged fixed point that is ultraviolet in the strong coupling but infrared in the weak coupling side. This finding suggests that the SU(3) 8-flavor system sits at the opening of the conformal window

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

CryoFold: Determining protein structures and data-guided ensembles from cryo-EM density maps

Cryoelectron microscopy requires molecular modeling for refinement of structures. Ensemble models arrive at low free-energy molecular structures, but are computationally expensive and limited to resolving only small proteins. Here, we introduce CryoFold, a pipeline of molecular dynamics simulations that determines ensembles of protein structures by integrating density data of varying sparsity at 3–5 Å resolution with sequence information and coarse-grained topological knowledge of the protein folds. We present six examples, folding proteins between 72 and 2,000 residues, including large membrane and multi-domain systems, and results from two Electron Microscopy Data Bank (EMDB) competitions. Driven by data from a single state, CryoFold discovers ensembles of common low-energy models together with rare low-probability structures that capture the equilibrium distribution of proteins constrained by the density maps. Many of these conformations are experimentally validated and functionally relevant. We arrive at a set of best practices for data-guided protein folding that are controlled using a Python graphical user interface (GUI).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Stochastic equilibrium Raman spectroscopy (STERS)

In this manuscript, we propose a new method for cavity- and surface-enhanced Raman spectroscopy (SERS) with improved temporal resolution in the measurement of stochastic Raman spectral fluctuations. Our approach combines Fourier spectroscopy and photon correlation to decouple the integration time from the temporal resolution. Using statistical optics Monte Carlo simulations, we establish the relationship between time resolution and Raman signal strength, revealing that typical Raman spectral fluctuations, commensurate with molecular conformational dynamics, can theoretically be resolved on micro- to millisecond timescales. The method can further extract average single-molecule dynamics from small sub-ensembles, thereby potentially mitigating challenges in achieving strictly single-molecule isolation on SERS substrates.

Cobb-Bruno, Colburn [University of California, Ber↗
DOE OSTI2025

Conformational free energy landscape of a glutamate transporter and microscopic details of its transport mechanism

Removing glutamate from the synaptic cleft is vital for proper function of the brain. Excitatory amino acid transporters mediate this process by uptaking the neurotransmitter from the synaptic cleft back to the cell after its release. The archaeal homolog, Glt Ph , an aspartate transporter fromPyrococcus horikoshii, presents the best structurally characterized model for this family of transporters. In order to transport, Glt Ph undergoes elevator-like conformational changes between inward-facing (IF) and outward-facing (OF) states. Here, we characterize, at an atomic level, the OF⇌IF transition of Glt Ph in differentapo/bound states using a combination of ensemble-based enhanced sampling techniques, employing more than two thousand of coupled simulation replicas of membrane-embedded Glt Ph . The resulting free-energy profiles portray the transition ofapo/bound states as a complex four-stage process, while sodium binding alone locks the structure in one of its states. Along the transition, the transport domain (TD) disengages from the scaffold domain (SD), allowing it to move as a piston sliding vertically with respect to the membrane during the elevator-like motion of TD. Lipid interactions with residues comprising the SD–TD interface directly influence the large-scale conformational changes and, consequently, the energetics of transport. Structural intermediates formed during the transition leak water molecules and may correlate to the uncoupled Cl − ion conductance observed experimentally in both prokaryotic and mammalian glutamate transporters. Mechanistic insights obtained from our study provide a structural framework for better development of therapeutic for neurological disorders.

Science & Technology - Other Topics↗
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 OSTI2024

Tailoring Hierarchical Structure and Rare Earth Affinity of Compositionally Identical Polymers via Sequence Control

Macromolecule sequence, structure, and function are inherently intertwined. While well-established relationships exist in proteins, they are more challenging to define for synthetic polymer nanoparticles due to their molecular weight, sequence, and conformational dispersities. Furthermore, to explore the impact of sequence on nanoparticle structure, we synthesized a set of 16 compositionally identical, sequence-controlled polymers with distinct monomer patterning of dimethyl acrylamide and a bioinspired, structure-driving di(phenylalanine) acrylamide (FF). Sequence control was achieved through multiblock polymerizations, yielding unique ensembles of polymer sequences which were simulated by kinetic Monte Carlo simulations. Systematic analysis of the global (tertiary- and quaternary-like) structure in this amphiphilic copolymer series revealed the effect of multiple sequence descriptors: the number of domains, the hydropathy of terminal domains, and the patchiness (density) of FF within a domain, each of which impacted both chain collapse and the distribution of single- and multichain assemblies. Furthermore, both the conformational freedom of chain segments and local-scale, β-sheet-like interactions were sensitive to the patchiness of FF. To connect sequence, structure, and target function, we evaluated an additional series of nine sequence-controlled copolymers as sequestrants for rare earth elements (REEs) by incorporating a functional acrylic acid monomer into select polymer scaffolds. We identified key sequence variables that influence the binding affinity, capacity, and selectivity of the polymers for REEs. Collectively, these results highlight the potential of and boundaries of sequence control via multiblock polymerizations to drive primary sequence ensembles hierarchical structures, and ultimately the functionality of compositionally identical polymeric materials.

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

Steering Molecular Dynamics Simulations of Membrane-Associated Proteins with Neutron Reflection Results

Herein, we present a novel method to incorporate structural results from surface-sensitive scattering, such as X-ray or neutron reflectometry, into molecular dynamics simulations. While reflectometry techniques generally provide a means to determine the molecular-scale structures of organized interfacial films, they were recently shown to offer the capability to characterize the structures of protein-membrane complexes supported by a solid substrate. One-dimensional information inherent in the experimental results is used in the form of component volume occupancy (CVO) profiles, which describe the distribution of molecular components within an interfacial architecture, to construct real-space constraints in the form of a biasing potential for the simulation that vanishes when the simulated and experimental profiles agree. This approach improves the correspondence between simulation and experiment, as shown in the re-evaluation of an neutron-reflection-derived structure which was approximated by an independent molecular dynamics simulation in earlier work, and it also leads to faster equilibration of ensemble structures. We further show that time averaging the CVO profile that develops in the simulation while biasing with this approach permits fluctuations about the average that are necessary for conformational exploration of the system. This method is particularly valuable for studies of proteins at interfaces that contain disordered regions since the conformation of such regions is difficult to judge from the analysis of one-dimensional experimental profiles and may take prohibitively long to equilibrate in simulations.

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

Uncovering a membrane-distal conformation of KRAS available to recruit RAF to the plasma membrane

The small GTPase KRAS is localized at the plasma membrane where it functions as a molecular switch, coupling extracellular growth factor stimulation to intracellular signaling networks. In this process, KRAS recruits effectors, such as RAF kinase, to the plasma membrane where they are activated by a series of complex molecular steps. Defining the membrane-bound state of KRAS is fundamental to understanding the activation of RAF kinase and in evaluating novel therapeutic opportunities for the inhibition of oncogenic KRAS-mediated signaling. We combined multiple biophysical measurements and computational methodologies to generate a consensus model for authentically processed, membrane-anchored KRAS. In contrast to the two membrane-proximal conformations previously reported, we identify a third significantly populated state using a combination of neutron reflectivity, fast photochemical oxidation of proteins (FPOP), and NMR. In this highly populated state, which we refer to as “membrane-distal” and estimate to comprise ~90% of the ensemble, the G-domain does not directly contact the membrane but is tethered via its C-terminal hypervariable region and carboxymethylated farnesyl moiety, as shown by FPOP. Subsequent interaction of the RAF1 RAS binding domain with KRAS does not significantly change G-domain configurations on the membrane but affects their relative populations. Overall, our results are consistent with a directional fly-casting mechanism for KRAS, in which the membrane-distal state of the G-domain can effectively recruit RAF kinase from the cytoplasm for activation at the membrane.

KRAS↗
DOE OSTI2023

Nonperturbative Exploration of Near-Conformal Dynamics

This project had three major goals: Explore qheories of composite dark matter and composite Higgs bosons through the work of the Lattice Strong Dynamics (LSD) collaboration.; Extend lattice gauge theories to curved manifolds through the work of the Quantum Finite Elements (QFE) collaboration.; Explore normalizing flow neural networks for producing ensembles of lattice gauge theory configurations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Entangled spin-polarized excitons from singlet fission in a rigid dimer

Abstract Singlet fission, a process that splits a singlet exciton into a biexciton, has promise in quantum information. We report time-resolved electron paramagnetic resonance measurements on a conformationally well-defined acene dimer molecule, TIPS-BP1', designed to exhibit strongly state-selective relaxation to specific magnetic spin sublevels. The resulting optically pumped spin polarization is a nearly pure initial state from the ensemble. The long-lived spin coherences modulate the signal intrinsically, allowing a measurement scheme that substantially removes noise and uncertainty in the magnetic resonance spectra. A nonadiabatic transition theory with a minimal number of spectroscopic parameters allows the quantitative assignment and interpretation of the spectra. In this work, we show that the rigid dimer TIPS-BP1' supports persistent spin coherences at temperatures far higher than those used in conventional superconducting quantum hardware.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Early events in G-quadruplex folding captured by time-resolved small-angle X-ray scattering

Abstract Time-resolved small-angle X-ray experiments are reported here that capture and quantify a previously unknown rapid collapse of the unfolded oligonucleotide as an early step in the folding of hybrid 1 and hybrid 2 telomeric G-quadruplex structures. The rapid collapse, initiated by a pH jump, is characterized by an exponential decrease in the radius of gyration from 24.3 to 12.6 Å. The collapse is monophasic and is complete in <600 ms. Additional hand-mixing pH-jump kinetic studies show that slower kinetic steps follow the collapse. The folded and unfolded states at equilibrium were further characterized by SAXS studies and other biophysical tools, showing that G4 unfolding was complete at alkaline pH, but not in LiCl solution as is often claimed. The SAXS Ensemble Optimization Method analysis reveals models of the unfolded state as a dynamic ensemble of flexible oligonucleotide chains with a variety of transient hairpin structures. These results suggest a G4 folding pathway in which a rapid collapse, analogous to molten globule formation seen in proteins, is followed by a confined conformational search within the collapsed particle to form the native contacts ultimately found in the stable folded form.

Biochemistry & Molecular Biology↗