Search NASA⌕ Search

SEARCH · Search NASA

Results for “ensemble 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computing Free Energies with Fluctuation Relations on Quantum Computers

One of the most promising applications for quantum computers is the dynamic simulation of quantum materials. Current hardware, however, sets stringent limitations on how long such simulations can run before decoherence begins to corrupt results. The Jarzynski equality, a fluctuation theorem that allows for the computation of equilibrium free energy differences from an ensemble of short, non-equilibrium dynamics simulations, can make use of such short-time simulations on quantum computers. Here, we present a quantum algorithm based on the Jarzynski equality for computing free energies of quantum materials. We demonstrate our algorithm using the transverse field Ising model on both a quantum simulator and real quantum hardware. As the free energy is a central thermodynamic property that allows one to compute virtually any equilibrium property of a physical system, the ability to perform this algorithm for larger quantum systems in the future has implications for a wide range of applications including the construction of phase diagrams, prediction of transport properties and reaction constants, and computer-aided drug design.

Bassman, Lindsay↗
DOE OSTI2025

In Situ/Operando Probing of Dynamic Phase Structures of Alumina‐Supported Ultrasmall Copper‐Gold Alloy Nanoparticles Under Reaction Conditions

Abstract The ability to control phase structures and surface sites of ultrasmall alloy nanoparticles under reaction conditions is essential for preparing catalysts by design. This is, however, challenging due to limited understanding of the atomic‐scale phases and their correlation with the ensemble‐averaged structures and activities of catalysts during catalytic reactions. We reveal here a dynamic structural stability of alumina‐supported ultrasmall and equiatomic copper‐gold alloy nanoparticles under reaction conditions as a model system in the in situ/operando study. In situ atomic‐scale morphological tracking under oxygen reveals temperature‐dependent dynamic crystalline‐amorphous dual‐phase structures, showing dynamic stability over an elevated temperature range. This atomic‐scale dynamic phase stability coincides with a “conversion plateau” observed for carbon monoxide oxidation on the catalyst. It is substantiated by the stable lattice ordering/disordering structures and surface sites with oscillatory characteristics shown by operando ensemble‐average structural tracking of the catalyst during the oxidation reaction. The understanding of the atomic‐scale dynamic phase structures in correlation with the ensemble‐average dynamic ordering/disordering phase structures and surface sites provides fresh insights into the unique synergy of the supported alloy nanoparticles. This understanding has implications for the design and structural tuning of active and stable ultrasmall alloy catalysts under elevated temperatures.

Chemistry↗
DOE OSTI2023

Robust projection of East Asian summer monsoon rainfall based on dynamical modes of variability

Abstract The Asian monsoon provides the freshwater that a large population in Asia depends on, but how anthropogenic climate warming may alter this key water source remains unclear. This is partly due to the prevailing point-wise assessment of climate projections, even though climate change patterns are inherently organized by dynamics intrinsic to the climate system. Here, we assess the future changes in the East Asian summer monsoon precipitation by projecting the precipitation from several large ensemble simulations and CMIP6 simulations onto the two leading dynamical modes of internal variability. The result shows a remarkable agreement among the ensembles on the increasing trends and the increasing daily variability in both dynamical modes, with the projection pattern emerging as early as the late 2030 s. The increase of the daily variability of the modes heralds more monsoon-related hydrological extremes over some identifiable East Asian regions in the coming decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effect of interfacial structures on phonon transport across atomically precise Si/Al heterojunctions

Phonons are important carriers of energy and information in many cryogenic devices used for quantum information science and in fundamental physics experiments such as dark matter detectors. In these systems phonon behavior can be dominated by interfaces, hence there is increasing demand for a more detailed understanding of interfacial phonon transport in relevant material systems. Previous studies have focused on understanding thermal transport over the entire phonon spectrum at and above room temperature. At ultra-low temperatures, however, knowledge is missing regarding athermal phonon behavior due to the challenge in modeling the extreme conditions in micro-scale, heterogeneous cryogenic systems, as well as extracting single-phonon information from a large ensemble. In this work, we delineate interfacial phonon transport using a combination of classical molecular dynamics (MD) and phonon wave packet simulations, to illustrate the consistency and differences between the ensemble- and single-phonon dynamics. We consider three single-crystal Si surface reconstructions --- (1×1), (v3×v3) and (7×7) --- and model both experimentally observed Si(1×1)/Al interfaces and hypothesized Si(v3×v3)/Al and Si(7×7)/Al interfaces. The overall interfacial thermal conductance calculated from non-equilibrium MD shows that for the Si(1×1)/Al system, the presence of Al twin grain boundaries can hinder phonon transport and reduce thermal conductance by 2--12% relative to single-crystal Al; whereas the Si (v3×v3) and (7×7) reconstructions can enhance it by 6--19%. Normal mode decomposition reveals that both the increase and decrease in conductance are related to inelastic phonon scattering. Single-phonon wave packet simulations show results consistent with NEMD, while further suggesting that phonon polarization conversion is significant even when elastic transmission dominates, and that the interfacial structures have anisotropic impacts on atomic vibrations along different lattice directions. Our findings suggest alternative perspectives in achieving selective phonon transport through controlling interfacial structures of materials with atomically-precise fabrication techniques, and that the phonon wave packet formalism is a potentially powerful method for developing a detailed understanding of non-equilibrium phenomena in the low-temperature limit.

Lu, Zexi↗
DOE OSTI2025

CMIP6-based Dynamically Downscaled Hydroclimate Projection over the Conterminous US

This dataset presents a suite of downscaled hydro-climate projections over the conterminous United States (CONUS) based on multiple selected Global Climate Models (GCMs) from the Coupled Models Intercomparison Project phase 6 (CMIP6). The CMIP6 GCMs are downscaled dynamically using Regional Climate Model version 4 (RegCM4). Each ensemble member covers the 1980-2019 baseline and 2020-2059 near-future periods under the high-end (SSP585) emission scenario. This dataset is formulated to support the SECURE Water Act Section 9505 Assessment for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details on this dataset, please refer to Kao et al. (2022) and Rastogi et al. (2022).

13 HYDRO ENERGY↗
DOE OSTI2022

Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems

Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.

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

Distribution of centrality measures on undirected random networks via the cavity method

The Katz centrality of a node in a complex network is a measure of the node’s importance as far as the flow of information across the network is concerned. For ensembles of locally tree-like undirected random graphs, this observable is a random variable. Its full probability distribution is of interest but difficult to handle analytically because of its “global” character and its definition in terms of a matrix inverse. Leveraging a fast Gaussian Belief Propagation-Cavity algorithm to solve linear systems on tree-like structures, we show that i) the Katz centrality of a single instance can be computed recursively in a very fast way, and ii) the probability P ( K ) that a random node in the ensemble of undirected random graphs has centrality K satisfies a set of recursive distributional equations, which can be analytically characterized and efficiently solved using a population dynamics algorithm. We test our solution on ensembles of Erdős-Rényi and Scale Free networks in the locally tree-like regime, with excellent agreement. The analytical distribution of centrality for the configuration model conditioned on the degree of each node can be employed as a benchmark to identify nodes of empirical networks with over- and underexpressed centrality relative to a null baseline. We also provide an approximate formula based on a rank- 1 projection that works well if the network is not too sparse, and we argue that an extension of our method could be efficiently extended to tackle analytical distributions of other centrality measures such as PageRank for directed networks in a transparent and user-friendly way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Autodifferentiable Ensemble Kalman Filters

Data assimilation is concerned with sequentially estimating a temporally evolving state. This task, which arises in a wide range of scientific and engineering applications, is particularly challenging when the state is high-dimensional and the state-space dynamics are unknown. This paper introduces a machine learning framework for learning dynamical systems in data assimilation. Here, our auto-differentiable ensemble Kalman filters (AD-EnKFs) blend ensemble Kalman filters for state recovery with machine learning tools for learning the dynamics. In doing so, AD-EnKFs leverage the ability of ensemble Kalman filters to scale to high-dimensional states and the power of automatic differentiation to train high-dimensional surrogate models for the dynamics. Numerical results using the Lorenz-96 model show that AD-EnKFs outperform existing methods that use expectation-maximization or particle filters to merge data assimilation and machine learning. In addition, AD-EnKFs are easy to implement and require minimal tuning.

autodifferentiation↗
DOE OSTI2023

Interplay between disorder, local relaxation, and collective behavior for an ensemble of emitters outside versus inside a cavity

The interplay between collective optical response and molecular static and dynamic disorder is studied using simple effective Hamiltonians for an ensemble of two-level emitters inside and outside a single-mode cavity. We model environmental disorder by randomly modulating the molecular transition frequencies and the coupling between the emitters and the electromagnetic field. We also consider the effects of intermolecular interactions and orientational disorder. We investigate how these effects lead to new features in the steady-state absorption (outside the cavity), transmission spectra (inside the cavity), and the yield of local molecular processes such as a unimolecular reaction. Outside the cavity, the collective behavior is manifested in the linewidth of the steady-state absorption, the emission spectrum, and the local chemical yield. Inside the cavity, however, the collective behavior primarily determines the Rabi splitting. The effects of intermolecular interactions under orientational disorder are also studied. For the most part, for all types of disorder, if we increase disorder, we find a reduction in the collective nature of the molecular response (smaller effective N) and therefore the Rabi splitting contraction occurs with orientational disorder. Furthermore, we find that static disorder is more destructive to collective behavior than dynamic disorder.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Real-time tracking of protein unfolding with time-resolved x-ray solution scattering

The correct folding of proteins is of paramount importance for their function, and protein misfolding is believed to be the primary cause of a wide range of diseases. Protein folding has been investigated with time-averaged methods and time-resolved spectroscopy, but observing the structural dynamics of the unfolding process in real-time is challenging. Here, we demonstrate an approach to directly reveal the structural changes in the unfolding reaction. We use nano- to millisecond time-resolved x-ray solution scattering to probe the unfolding of apomyoglobin. The unfolding reaction was triggered using a temperature jump, which was induced by a nanosecond laser pulse. We demonstrate a new strategy to interpret time-resolved x-ray solution scattering data, which evaluates ensembles of structures obtained from molecular dynamics simulations. We find that apomyoglobin passes three states when unfolding, which we characterize as native, molten globule, and unfolded. The molten globule dominates the population under the conditions investigated herein, whereas native and unfolded structures primarily contribute before the laser jump and 30 μs after it, respectively. The molten globule retains much of the native structure but shows a dynamic pattern of inter-residue contacts. Our study demonstrates a new strategy to directly observe structural changes over the cause of the unfolding reaction, providing time- and spatially resolved atomic details of the folding mechanism of globular proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Catalytic site ensembles: A context to reexamine the Langmuir-Hinshelwood kinetic description

The Langmuir-Hinshelwood formalism describes catalytic reactions of Langmuirian surface species under the assumption that all adsorbates are randomly-distributed – enabling adjacency of surface-bound intermediates to be determined solely by coverages of single-site occupants. We demonstrate herein that this approximation is inappropriate even for simple catalytic reactions (e.g. A + A → A 2 ) and manifestly neglects islanding of slowly-consumed species and partitioning of highly-reactive species inherently engendered by ≥ two-site elementary steps (e.g. A*–A* → A 2(g) + *–*). Rigorous description of kinetically-consequential islanding/partitioning phenomena requires explicit description of the coverage and chemical dynamics of all multi-site ensembles. Higher-order, ensemble-specific rate terms identify the particular microscopic events relevant to each ensemble, and, in doing so, reveal that each elementary step (e.g. A (g) adsorption) describes not one event (e.g. A (g) + * → A*) , but a sum over all ensemble-specific paths (e.g. A (g) + *–* → A*–* and A (g) + A*–* → A*–A*). De-convoluting each elementary step into its constituent multi-site paths proffers kinetic detail otherwise inaccessible – enabling (i) identification of rate- and selectivity-determining site ensembles, (ii) calculation of rates and degrees of rate control of ensemble-specific elementary steps, (iii) incorporation of adsorbate surface diffusion, (iv) incorporation of lateral adsorbate interactions, and (v) quantitative description of catalysis of multi-site-occupying intermediates (e.g. *C n H m * species in hydrocarbon (de-)hydrogenation and C-C bond coupling/cleavage reactions) which we demonstrate here is inaccessible to the Langmuir-Hinshelwood formalism even if adsorbate surface diffusion is infinitely-fast.

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

Collaborative Research: Advancing Arctic Climate Projection Capability at Seasonal to Decadal Scales (Final Technical Report)

The Regional Arctic System Model (RASM) at process resolving configurations has been used to (i) advance understanding of physical processes and feedbacks involved in Arctic amplification and (ii) understand and potentially reduce uncertainty in prediction of arctic climate change at seasonal to decadal scales. RASM consists the atmosphere (Weather and Research Forecasting model, WRF), ocean (Parallel Ocean Program, POP), sea ice (CICE), land hydrology (Variable Infiltration Capacity model, VIC), river routing scheme (RVIC), marine biogeochemistry components and the coupling framework (CPL7). Its domain is pan-Arctic, with the atmosphere and land components configured on a 50-km or 25-km grid and four configurations of the ocean and sea ice components: 1/12°(~9.3km) or 1/48°(~2.4km) and 45 or 60 vertical layers. These RASM configurations have been motivated by the emerging exascale capability for high performance computing to improve model fidelity. The dynamical downscaling of reanalysis allows comparison of RASM results with observations in place and time to: (i) advance system level understanding of physical processes and coupling involved in an event, (ii) optimize model parameter space, (iii) diagnose and reduce model biases and (iv) produce realistic and consistent across all the components initial conditions for predictions and predictability studies, which are all unique capabilities not available in global Earth System Models (ESMs). An evaluation of RASM 1.0 (Cassano et al. 2017) revealed that it had a cold bias over the oceans and a warm bias over land areas due largely to cloud and radiation biases in the model, with too little cloud cover simulated over land and too much cloud cover simulated over sub-polar oceans. This study has motivated an upgrade to WRF version 3.7.1 in RASM and allowed for the inclusion of the radiative impact of convective clouds. A variety of atmospheric physics parameterizations were evaluated against observations (e.g. data from the Arctic Clouds in Summer Experiment (ACSE); Sedlar et al. 2020) to identify an optimal suite of WRF physics options in RASM. The RASM with the optimized WRF physics were used to study the impact of strong mesoscale winds over the ocean around the southern tip of Greenland (DuVivier and Cassano 2016) and their impact on oceanic convection (DuVivier et al. 2017a). Data from the PolarWinds field campaign were used to evaluate WRF boundary layer physics and resolution impacts on the simulation of a Greenland barrier wind event (DuVivier et al. 2017b). The RVIC streamflow routing model has been implemented in RASM to realistically represent high-resolution streamflow processes (Hamman et al. 2017) and to couple the land buoyancy fluxes to the ocean. The RASM-RVIC high-resolution data set of all coastal freshwater fluxes for the Arctic drainage basin and surrounding areas for 1979-2014 was published as a separate product (https://doi.org/10.5281/zenodo.293037). The fidelity of atmospheric momentum transfer to and the response of polar marine Ekman layer in RASM and Community Earth System Model (CESM) was investigated by Roberts et al. (2015). The increased frequency of oceanic flux exchange in CESM, following the RASM guidance, caused a considerable increase in the median inertial ice speed across the Southern Ocean and parts of the Arctic. A comprehensive evaluation of the RASM1.0 atmosphere-ocean-sea ice-land interface was completed by Brunke et al. (2018). RASM was also demonstrated for its capability to simulate extreme events in agreement with observations in space and time (Lee et al. submitted). In particular, the development of three open water events, known as polynyas, have been simulated north of Greenland in February of 2011, 2017 and 2018, in agreement with satellite observations for the past four decades. The optimized RASM sea ice results have been favorably evaluated against satellite observations and a subset of eleven CMIP6 models (Watts et al. submitted). In a complementary project, Jin et al. (2018) have shown that RASM with higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice conditions, concentrations of nutrients and ice algae, in comparison to results from the coarse-resolution (1°) CESM. In two other complementary studies, RASM results were used (i) to explain areas of concentrated use by bowhead whales, the seasonal progression in the use, and the physical environment within those areas (Citta et al. 2015) and (ii) for a synthesis of fall bowhead whales distribution and migration in the Bering-Chukchi-Beaufort (BCB) Sea to investigate whale movements and feeding to the local ocean hydrography and currents (Citta et al. 2018). However, the multi-decadal output from the CESM Large Ensemble yielded unrealistic forcing. Instead, the shorter NCEP CFSv2 9-month forecasts were successfully tested and afforded an increased ensemble size (~30) to demonstrate gains of dynamical downscaling at sub-seasonal to intra-annual time scales. The improved model physics and coupling among RASM model components have yielded more realistic representation of the sea ice cover and consistent across all model components initial conditions. Consequently, RASM demonstrates significant gains compared to simulation of sea ice in the NCEP reanalysis. In addition, RASM 6-month ensemble forecasts yield very realistic sea ice distribution, which demonstrates both significant gains of dynamical downscaling and the continued impact of the initial conditions on forecasts out to 6 months (https://nps.edu/web/rasm/predictions). A follow up study using RASM for dynamical downscaling of the more realistic CESM initialized Decadal Prediction Large Ensemble output is currently ongoing as part of the DOE RGMA HiLAT-RASM project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Progress in Normalizing Flows for 4d Gauge Theories

Normalizing flows have arisen as a tool to accelerate Monte Carlo sampling for lattice field theories. This work reviews recent progress in applying normalizing flows to 4-dimensional nonabelian gauge theories, focusing on two advancements: an architectural improvement referred to as learned active loops, and the application of correlated ensemble methods to QCD with N f = 2 dynamical fermions.

Abbott, Ryan [Massachusetts Institute of Technolog↗
DOE OSTI2021

Real-time interfacial electron dynamics revealed through temporal correlations in x-ray photoelectron spectroscopy

We present a novel technique to monitor dynamics in interfacial systems through temporal correlations in x-ray photoelectron spectroscopy (XPS) signals. To date, the vast majority of time-resolved x-ray spectroscopy techniques rely on pump–probe schemes, in which the sample is excited out of equilibrium by a pump pulse, and the subsequent dynamics are monitored by probe pulses arriving at a series of well-defined delays relative to the excitation. By definition, this approach is restricted to processes that can either directly or indirectly be initiated by light. It cannot access spontaneous dynamics or the microscopic fluctuations of ensembles in chemical or thermal equilibrium. Enabling this capability requires measurements to be performed in real (laboratory) time with high temporal resolution and, ultimately, without the need for a well-defined trigger event. The time-correlation XPS technique presented here is a first step toward this goal. The correlation-based technique is implemented by extending an existing optical-laser pump/multiple x-ray probe setup by the capability to record the kinetic energy and absolute time of arrival of every detected photoelectron. The method is benchmarked by monitoring energy-dependent, periodic signal modulations in a prototypical time-resolved XPS experiment on photoinduced surface-photovoltage dynamics in silicon, using both conventional pump–probe data acquisition, and the new technique based on laboratory time. The two measurements lead to the same result. The findings provide a critical milestone toward the overarching goal of studying equilibrium dynamics at surfaces and interfaces through time correlation-based XPS measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A Quantum-Based Approach to Predict Primary Radiation Damage in Polymeric Networks

Initial atomistic-level radiation damage in chemically reactive materials is thought to induce reaction cascades that can result in undesirable degradation of macroscale properties. Ensembles of quantum-based molecular dynamics (QMD) simulations can accurately predict these cascades, but extracting chemical insights from the many underlying trajectories is a labor-intensive process that can require substantial a priori intuition. We develop here a general and automated graph-based approach to extract all chemically distinct structures sampled in QMD simulations and apply our approach to predict primary radiation damage of polydimethylsiloxane (PDMS), the main constituent of silicones. A postprocessing protocol is developed to identify underlying polymer backbone structures as connected components in QMD trajectories. Furthermore, these backbones form a repository of radiation-damaged structures. A scheme for extracting and updating a library of isomorphically distinct structures is proposed to identify the spanning set and aid chemical interpretation of the repository. The analyses are applied to ensembles of cascade QMD simulations in which the four element types in PDMS are selectively excited in primary knock-on atom events. Our approach reveals a much higher degree of combinatorial complexity in this system than was inferred through radiolysis experiments. Probabilities are extracted for radiation-induced network changes including formation of branch points, carbon linkages, cycles, bond scissions, and carbon uptake into the Si–O siloxane backbone network. The general analysis framework presented here is readily extendable to modeling chemical degradation of other polymers and molecular materials and provides a basis for future quantum-informed multiscale modeling of radiation damage.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The role of SAXS and molecular simulations in 3D structure elucidation of a DNA aptamer against lung cancer

Aptamers are short, single-stranded DNA or RNA oligonucleotide molecules that function as synthetic analogs of antibodies and bind to a target molecule with high specificity. Aptamer affinity entirely depends on its tertiary structure and charge distribution. Therefore, length and structure optimization are essential for increasing aptamer specificity and affinity. Here, we present a general optimization procedure for finding the most populated atomistic structures of DNA aptamers. Based on the existed aptamer LC-18 for lung adenocarcinoma, a new truncated LC-18 (LC-18t) aptamer LC-18t was developed. A three-dimensional (3D) shape of LC-18t was reported based on small-angle X-ray scattering (SAXS) experiments and molecular modeling by fragment molecular orbital or molecular dynamic methods. Molecular simulations revealed an ensemble of possible aptamer conformations in solution that were in close agreement with measured SAXS data. The aptamer LC-18t had stronger binding to cancerous cells in lung tumor tissues and shared the binding site with the original larger aptamer. The suggested approach reveals 3D shapes of aptamers and helps in designing better affinity probes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Ligand modulation of the conformational dynamics of the A 2A adenosine receptor revealed by single-molecule fluorescence

G protein-coupled receptors (GPCRs) are the largest class of transmembrane proteins, making them an important target for therapeutics. Activation of these receptors is modulated by orthosteric ligands, which stabilize one or several states within a complex conformational ensemble. The intra- and inter-state dynamics, however, is not well documented. Here, we used single-molecule fluorescence to measure ligand-modulated conformational dynamics of the adenosine A 2A receptor (A 2A R) on nanosecond to millisecond timescales. Experiments were performed on detergent-purified A 2R in either the ligand-free (apo) state, or when bound to an inverse, partial or full agonist ligand. Single-molecule Förster resonance energy transfer (smFRET) was performed on detergent-solubilized A 2A R to resolve active and inactive states via the separation between transmembrane (TM) helices 4 and 6. The ligand-dependent changes of the smFRET distributions are consistent with conformational selection and with inter-state exchange lifetimes ≥ 3 ms. Local conformational dynamics around residue 229 6.31 on TM6 was measured using fluorescence correlation spectroscopy (FCS), which captures dynamic quenching due to photoinduced electron transfer (PET) between a covalently-attached dye and proximal aromatic residues. Global analysis of PET-FCS data revealed fast (150–350 ns), intermediate (50–60 μs) and slow (200–300 μs) conformational dynamics in A 2A R, with lifetimes and amplitudes modulated by ligands and a G-protein mimetic (mini-G s ). Most notably, the agonist binding and the coupling to mini-G s accelerates and increases the relative contribution of the sub-microsecond phase. Molecular dynamics simulations identified three tyrosine residues (Y112, Y288 7.53 , and Y290 7.55 ) as being responsible for the dynamic quenching observed by PET-FCS and revealed associated helical motions around residue 229 6.31 on TM6. As a result, this study provides a quantitative description of conformational dynamics in A 2A R and supports the idea that ligands bias not only GPCR conformations but also the dynamics within and between distinct conformational states of the receptor.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Molecular-dynamics simulation methods for macromolecular crystallography

It is investigated whether molecular-dynamics (MD) simulations can be used to enhance macromolecular crystallography (MX) studies. Historically, protein crystal structures have been described using a single set of atomic coordinates. Because conformational variation is important for protein function, researchers now often build models that contain multiple structures. Methods for building such models can fail, however, in regions where the crystallographic density is difficult to interpret, for example at the protein–solvent interface. To address this limitation, a set of MD–MX methods that combine MD simulations of protein crystals with conventional modeling and refinement tools have been developed. In an application to a cyclic adenosine monophosphate-dependent protein kinase at room temperature, the procedure improved the interpretation of ambiguous density, yielding an alternative water model and a revised protein model including multiple conformations. The revised model provides mechanistic insights into the catalytic and regulatory interactions of the enzyme. The same methods may be used in other MX studies to seek mechanistic insights.

59 BASIC BIOLOGICAL SCIENCES↗