Structural Dynamics-Us最新文献

The use of laser-plasma-based x-ray sources is discussed, with a view to carrying out time-resolved x-ray absorption spectroscopy measurements, down to the femtosecond timescale. A review of recent experiments performed by our team is presented. They concern the study of the nonequilibrium transition of metals from solid to the warm dense regime, which imposes specific constraints (the sample being destroyed after each shot). Particular attention is paid to the description of experimental devices and methodologies. Two main types of x-ray sources are compared, respectively, based on the emission of a hot plasma, and on the betatron radiation from relativistic electrons accelerated by laser.

The direct observation of the structure of micrometer-sized vapor-deposited ice is performed at Pohang Accelerator Laboratory x-ray free electron laser (PAL-XFEL). The formation of micrometer-sized ice crystals and their structure is important in various fields, including atmospheric science, cryobiology, and astrophysics, but understanding the structure of micrometer-sized ice crystals remains challenging due to the lack of direct observation. Using intense x-ray diffraction from PAL-XFEL, we could observe the structure of micrometer-sized vapor-deposited ice below 150 K with a thickness of 2-50 μm grown in an ultrahigh vacuum chamber. The structure of the ice grown comprises cubic and hexagonal sequences that are randomly arranged to produce a stacking-disordered ice. We observed that ice with a high cubicity of more than 80% was transformed to partially oriented hexagonal ice when the thickness of the ice deposition grew beyond 5 μm. This suggests that precise temperature control and clean deposition conditions allow μm-thick ice films with high cubicity to be grown on hydrophilic Si₃N₄ membranes. The low influence of impurities could enable in situ diffraction experiments of ice nucleation and growth from interfacial layers to bulk ice.

Mechanistic studies of Geobacillus stearothermophilus tryptophanyl-tRNA synthetase (TrpRS) afford an unusually detailed description-the escapement mechanism-for the distinct steps coupling catalysis to domain motion, efficiently converting the free energy of ATP hydrolysis into biologically useful alternative forms of information and work. Further elucidation of the escapement mechanism requires understanding thermodynamic linkages between domain configuration and conformational stability. To that end, we compare experimental thermal melting of fully liganded and apo TrpRS with a computational simulation of the melting of its fully liganded form. The simulation also provides important structural cameos at successively higher temperatures, enabling more confident interpretation. Experimental and simulated melting both proceed through a succession of three transitions at successively higher temperature. The low-temperature transition occurs at approximately the growth temperature of the organism and so may be functionally relevant but remains too subtle to characterize structurally. Structural metrics from the simulation imply that the two higher-temperature transitions entail forming a molten globular state followed by unfolding of secondary structures. Ligands that stabilize the enzyme in a pre-transition (PreTS) state compress the temperature range over which these transitions occur and sharpen the transitions to the molten globule and fully denatured states, while broadening the low-temperature transition. The experimental enthalpy changes provide a key parameter necessary to convert changes in melting temperature of combinatorial mutants into mutationally induced conformational free energy changes. The TrpRS urzyme, an excerpted model representing an early ancestral form, containing virtually the entire catalytic apparatus, remains largely intact at the highest simulated temperatures.

The important role of structural dynamics in protein function is widely recognized. Thermal or B-factors and their anisotropy, seen in x-ray analysis of protein structures, report on the presence of atomic coordinate heterogeneity that can be attributed to motion. However, their quantitative evaluation in terms of protein dynamics by x-ray ensemble refinement remains challenging. NMR spectroscopy provides quantitative information on the amplitudes and time scales of motional processes. Unfortunately, with a few exceptions, the NMR data do not provide direct insights into the atomic details of dynamic trajectories. Residual dipolar couplings, measured by solution NMR, are very precise parameters reporting on the time-averaged bond-vector orientations and may offer the opportunity to derive correctly weighted dynamic ensembles of structures for cases where multiple high-resolution x-ray structures are available. Applications to the SARS-CoV-2 main protease, M^pro, and ubiquitin highlight this complementarity of NMR and crystallography for quantitative assessment of internal motions.

The field of time-resolved macromolecular crystallography has been expanding rapidly after free electron lasers for hard x rays (XFELs) became available. Techniques to collect and process data from XFELs spread to synchrotron light sources. Although time-scales and data collection modalities can differ substantially between these types of light sources, the analysis of the resulting x-ray data proceeds essentially along the same pathway. At the base of a successful time-resolved experiment is a difference electron density (DED) map that contains chemically meaningful signal. If such a difference map cannot be obtained, the experiment has failed. Here, a practical approach is presented to calculate DED maps and use them to determine structural models.

Landscape descriptions provide a framework for identifying functionally significant dynamic linkages in proteins but cannot supply details. Rate measurements of combinatorial mutations can implicate dynamic linkages in catalysis. A major difficulty is filtering dynamic linkages from the vastly more numerous static interactions that stabilize domain folding. The Geobacillus stearothermophilus (TrpRS) D1 switch is such a dynamic packing motif; it links domain movement to catalysis and specificity. We describe Thermofluor and far UV circular dichroism melting curves for all 16 D1 switch variants to determine their higher-order impact on unliganded TrpRS stability. A prominent transition at intermediate temperatures in TrpRS thermal denaturation is molten globule formation. Combinatorial analysis of thermal melting transcends the protein landscape in four significant respects: (i) bioinformatic methods identify dynamic linkages from coordinates of multiple conformational states. (ii) Relative mutant melting temperatures, δT_M, are proportional to free energy changes. (iii) Structural analysis of thermal melting implicates unexpected coupling between the D1 switch packing and regions of high local frustration. Those segments develop molten globular characteristics at the point of greatest complementarity to the chemical transition state and are the first TrpRS structures to melt. (iv) Residue F37 stabilizes both native and molten globular states; its higher-order interactions modify the relative intrinsic impacts of mutations to other D1 switch residues from those estimated for single point mutants. The D1 switch is a central component of an escapement mechanism essential to free energy transduction. These conclusions begin to relate the escapement mechanism to differential TrpRS conformational stabilities.

The dynamics of hydration water (HW) in 1,2-dimyristoyl-sn-glycero-3-phosphoethanolamine (DMPE) was investigated by means of quasi-elastic neutron scattering (QENS) and compared with those observed in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). The headgroup dynamics of DMPE was investigated using a mixture of tail-deuterated DMPE and D₂O, and the QENS profiles were interpreted as consisting of three modes. The fast mode comprised the rotation of hydrogen atoms in -NH₃⁺ and -CH₂- groups in the headgroup of DMPE, the medium-speed mode comprised fluctuations in the entire DMPE molecule, and the slow mode comprised fluctuations in the membrane. These interpretations were confirmed using molecular dynamics (MD) simulations. The HW dynamics analysis was performed on a tail-deuterated DMPE and H₂O mixture. The QENS profiles were analyzed in terms of three modes: (1) a slow mode, identified as loosely bound HW in the DMPC membrane; (2) a medium-speed mode similar to free HW in the DMPC membrane; and (3) a fast mode, identified as rotational motion. The relaxation time for the fast mode was approximately six times shorter than that of rotational water in DMPC, consistent with the results of terahertz time-domain spectroscopy. The activation energy of medium-speed HW in DMPE differed from that of free HW in DMPC, suggesting the presence of different hydration states or hydrogen-bonded networks around the phosphocholine and phosphoethanolamine headgroups.

We present spatially resolved measurements of the temperature of a flat liquid water microjet for varying ambient pressures, from vacuum to 100% relative humidity. The entire jet surface is probed in a single shot by a high-resolution infrared camera. Obtained 2D images are substantially influenced by the temperature of the apparatus on the opposite side of the infrared camera; a protocol to correct for the thermal background radiation is presented. In vacuum, we observe cooling rates due to water evaporation on the order of 10⁵ K/s. For our system, this corresponds to a temperature decrease in approximately 15 K between upstream and downstream positions of the flowing leaf. Making reasonable assumptions on the absorption of the thermal background radiation in the flatjet, we can extend our analysis to infer a thickness map. For a reference system, our value for the thickness is in good agreement with the one reported from white light interferometry.

Low-pass spectral analysis (LPSA) is a recently developed dynamics retrieval algorithm showing excellent retrieval properties when applied to model data affected by extreme incompleteness and stochastic weighting. In this work, we apply LPSA to an experimental time-resolved serial femtosecond crystallography (TR-SFX) dataset from the membrane protein bacteriorhodopsin (bR) and analyze its parametric sensitivity. While most dynamical modes are contaminated by nonphysical high-frequency features, we identify two dominant modes, which are little affected by spurious frequencies. The dynamics retrieved using these modes shows an isomerization signal compatible with previous findings. We employ synthetic data with increasing timing uncertainty, increasing incompleteness level, pixel-dependent incompleteness, and photon counting errors to investigate the root cause of the high-frequency contamination of our TR-SFX modes. By testing a range of methods, we show that timing errors comparable to the dynamical periods to be retrieved produce a smearing of dynamical features, hampering dynamics retrieval, but with no introduction of spurious components in the solution, when convergence criteria are met. Using model data, we are able to attribute the high-frequency contamination of low-order dynamical modes to the high levels of noise present in the data. Finally, we propose a method to handle missing observations that produces a substantial dynamics retrieval improvement from synthetic data with a significant static component. Reprocessing of the bR TR-SFX data using the improved method yields dynamical movies with strong isomerization signals compatible with previous findings.

The time-resolved x-ray absorption spectrum of the BT-1T cation (BT-1T⁺) is theoretically simulated in order to investigate the charge transfer reaction of the system. We employ both trajectory surface hopping and quantum dynamics to simulate the structural evolution over time and the changes in the state populations. To compute the static x-ray absorption spectra (XAS) of the ground and excited states, we apply both the time-dependent density functional theory and the coupled cluster singles and doubles method. The results obtained are in good agreement between the methods. It is, furthermore, found that the small structural changes that occur during the reaction have little effect on the static XAS. Hence, the tr-XAS can be computed based on the state populations determined from a nuclear dynamics simulation and one set of static XAS calculations, utilizing the ground state optimized geometry. This approach can save considerable computational resources, as the static spectra need not to be calculated for all geometries. As BT-1T is a relatively rigid molecule, the outlined approach should only be considered when investigating non-radiative decay processes in the vicinity of the Franck-Condon point.