Ultrafast anisotropic protein quake propagation after CO photodissociation in myoglobin

Glycerol-slaved 1H-1H NMR cross-relaxation in quasi-native lysozyme

1H-1H nuclear cross-relaxation experiments have been carried out with lysozyme in variable glycerol viscosity to study intramolecular motion, self-diffusion, and isotropic rigid-body rotational tumbling at 298 K, pH 3.8. Dynamics of intramolecular 1H-1H cross-relaxation rates, the increase in internuclear spatial distances, and lateral and rotational diffusion coefficients all show fractional viscosity dependence with a power law exponent κ in the 0.17-0.83 range. The diffusion coefficient of glycerol Ds with the bulk viscosity itself is non-Stokesian, having a fractional viscosity dependence on the medium viscosity (Ds ∼ η-κ, κ ≈ 0.71). The concurrence and close similarity of the fractional viscosity dependence of glycerol diffusion on the one hand, and diffusion and intramolecular cross-relaxation rates of the protein on the other lead to infer that relaxation of glycerol slaves protein relaxations. Glycerol-transformed native lysozyme to a quasi-native state does not affect the conclusion that both global and internal fluctuations are slaved to glycerol relaxation.

Tracking the structural dynamics of proteins with time-resolved X-ray solution scattering

Relevant events during protein function such as ligand binding/release and interaction with substrates or with light are often accompanied by out-of-equilibrium structural dynamics. Time-resolved experimental techniques have been developed to follow protein structural changes as they happen in real time after a given reaction-triggering event. Time-resolved X-ray solution scattering is a promising approach that bears structural sensitivity with temporal resolution in the femto-to-millisecond time range, depending on the X-ray source characteristics and the triggering method. Here we present the basic principles of the technique together with a description of the most relevant results recently published and a discussion on the computational methods currently developed to achieve a structural interpretation of the time-resolved X-ray solution scattering experimental data.

Time-resolved spectroscopic mapping of vibrational energy flow in proteins: Understanding thermal diffusion at the nanoscale

Vibrational energy exchange between various degrees of freedom is critical to barrier-crossing processes in proteins. Hemeproteins are well suited for studying vibrational energy exchange in proteins because the heme group is an efficient photothermal converter. The released energy by heme following photoexcitation shows migration in a protein moiety on a picosecond timescale, which is observed using time-resolved ultraviolet resonance Raman spectroscopy. The anti-Stokes ultraviolet resonance Raman intensity of a tryptophan residue is an excellent probe for the vibrational energy in proteins, allowing the mapping of energy flow with the spatial resolution of a single amino acid residue. This Perspective provides an overview of studies on vibrational energy flow in proteins, including future perspectives for both methodologies and applications.

Origin of Protein Quake: Energy Waves Conducted by a Precise Mechanical Machine

A long-standing challenge in protein biophysics is to understand protein quake in myoglobin─the structural dynamics responsible for redistributing the excess heme energy after photolysis. Despite extensive efforts, the molecular mechanism of this process remains elusive. Using the energy flow theory, we uncovered a fundamental new phenomenon: the heme energy is redistributed by sinusoidal waves with a ubiquitous fundamental frequency and two overtones. The energy waves emanate from the heme into the myoglobin backbone via a conduit of five consecutive dihedrals of the proximal histidine and then travel quickly along the backbone to reach sidechains across the protein. This mechanism is far more effective than the diffusion-based mechanism from previous studies because waves are systematic while diffusion is random. To propagate energy waves, coordinates must cooperate, resulting in collective modes that are singular vectors of the generalized work functional. These modes show task partitioning: a handful of high-energy modes generate large-scale breathing motion, which loosens up the protein matrix to enable hundreds of low-energy vibrational modes for energy transduction.

Enhanced Mobility during Diels-Alder Reaction: Results of Molecular Simulations

Recent measurements indicate enhanced mobility of solvent molecules during Diels-Alder (DA) and other common chemical reactions. We present results of molecular dynamics simulations of the last stages of the DA cycloaddition reaction, from the transition state configuration to product, of furfurylamine and maleimide in acetonitrile at reactant concentrations studied experimentally. We find enhanced mobility of solvent and reactant molecules up to at least a nanometer from the DA product over hundreds of picoseconds. Local heating is ruled out as a factor in the enhanced mobility observed in the simulations, which is instead found to be due to solvent relaxation following the formation of the DA product.

Impact of Ultrafast Electric Field Changes on Photoreceptor Protein Dynamics

Studies on photoreceptors provide a wealth of information on cofactor and protein dynamics on the microsecond to seconds time-scale. Up to now, ultrafast dynamics addresses mainly the cofactor or chromophore, but ultrafast protein dynamics are poorly understood. Increasing evidence show that protein responses can occur even faster than the cofactor dynamics. The causal reason for the ultrafast protein response cannot be explained by the localized cofactor excitation or its excited-state decay, alone. We propose a Coulomb interaction mechanism started by a shock wave and stabilized by a dipole moment change at least partially responsible for coherent oscillations in proteins, protonation changes, water dislocations, and protein changes prior to and beyond chromophore's excited-state decay. Photoexcitation changes the electron density distribution of the chromophore within a few femtoseconds: The Coulomb shock wave affects polar groups, hydrogen bonds, and protein bound water molecules. The process occurs on a time-scale even faster than excited-state decay of the chromophore. We discuss studies on selected photoreceptors in light of this mechanism and its impact on a detailed understanding of protein dynamics.

Reaction dynamics studied via femtosecond X-ray liquidography at X-ray free-electron lasers

X-ray free-electron lasers (XFELs) provide femtosecond X-ray pulses suitable for pump–probe time-resolved studies with a femtosecond time resolution. Since the advent of the first XFEL in 2009, recent years have witnessed a great number of applications with various pump–probe techniques at XFELs. Among these, time-resolved X-ray liquidography (TRXL) is a powerful method for visualizing structural dynamics in the liquid solution phase. Here, we classify various chemical and biological molecular systems studied via femtosecond TRXL (fs-TRXL) at XFELs, depending on the focus of the studied process, into (i) bond cleavage and formation, (ii) charge distribution and electron transfer, (iii) orientational dynamics, (iv) solvation dynamics, (v) coherent nuclear wavepacket dynamics, and (vi) protein structural dynamics, and provide a brief review on each category. We also lay out a plausible roadmap for future fs-TRXL studies for areas that have not been explored yet.

Femtosecond X-ray liquidography using X-ray free-electron lasers (XFELs) visualizes various aspects of reaction dynamics.

Advances and challenges in time-resolved macromolecular crystallography

Conformational changes within biological macromolecules control a vast array of chemical reactions in living cells. Time-resolved crystallography can reveal time-dependent structural changes that occur within protein crystals, yielding chemical insights in unparalleled detail. Serial crystallography approaches developed at x-ray free-electron lasers are now routinely used for time-resolved diffraction studies of macromolecules. These techniques are increasingly being applied at synchrotron radiation sources and to a growing diversity of macromolecules. Here, we review recent progress in the field, including visualizing ultrafast structural changes that guide the initial trajectories of light-driven reactions as well as capturing biologically important conformational changes on slower time scales, for which bacteriorhodopsin and photosystem II are presented as illustrative case studies.

Ultrafast coherent motion and helix rearrangement of homodimeric hemoglobin visualized with femtosecond X-ray solution scattering

Ultrafast motion of molecules, particularly the coherent motion, has been intensively investigated as a key factor guiding the reaction pathways. Recently, X-ray free-electron lasers (XFELs) have been utilized to elucidate the ultrafast motion of molecules. However, the studies on proteins using XFELs have been typically limited to the crystalline phase, and proteins in solution have rarely been investigated. Here we applied femtosecond time-resolved X-ray solution scattering (fs-TRXSS) and a structure refinement method to visualize the ultrafast motion of a protein. We succeeded in revealing detailed ultrafast structural changes of homodimeric hemoglobin involving the coherent motion. In addition to the motion of the protein itself, the time-dependent change of electron density of the hydration shell was tracked. Besides, the analysis on the fs-TRXSS data of myoglobin allows for observing the effect of the oligomeric state on the ultrafast coherent motion.

Femtosecond time-resolved X-ray solution scattering (fs-TRXSS) measurements provide information on the structural dynamics of proteins in solution. Here, the authors present a structure refinement method for the analysis of fs-TRXSS data and use it to characterise the ultrafast structural changes of homodimeric haemoglobin.

Capturing photochemical and photophysical transformations in iron complexes with ultrafast X-ray spectroscopy and scattering

Light-driven chemical transformations provide a compelling approach to understanding chemical reactivity with the potential to use this understanding to advance solar energy and catalysis applications. Capturing the non-equilibrium trajectories of electronic excited states with precision, particularly for transition metal complexes, would provide a foundation for advancing both of these objectives. Of particular importance for 3d metal compounds is characterizing the population dynamics of charge-transfer (CT) and metal-centered (MC) electronic excited states and understanding how the inner coordination sphere structural dynamics mediate the interaction between these states. Recent advances in ultrafast X-ray laser science has enabled the electronic excited state dynamics in 3d metal complexes to be followed with unprecedented detail. This review will focus on simultaneous X-ray emission spectroscopy (XES) and X-ray solution scattering (XSS) studies of iron coordination and organometallic complexes. These simultaneous XES-XSS studies have provided detailed insight into the mechanism of light-induced spin crossover in iron coordination compounds, the interaction of CT and MC excited states in iron carbene photosensitizers, and the mechanism of Fe-S bond dissociation in cytochrome c.

Interplays of electron and nuclear motions along CO dissociation trajectory in myoglobin revealed by ultrafast X-rays and quantum dynamics calculations

Ultrafast structural dynamics with different spatial and temporal scales were investigated during photodissociation of carbon monoxide (CO) from iron(II)-heme in bovine myoglobin during the first 3 ps following laser excitation. We used simultaneous X-ray transient absorption (XTA) spectroscopy and X-ray transient solution scattering (XSS) at an X-ray free electron laser source with a time resolution of 80 fs. Kinetic traces at different characteristic X-ray energies were collected to give a global picture of the multistep pathway in the photodissociation of CO from heme. In order to extract the reaction coordinates along different directions of the CO departure, XTA data were collected with parallel and perpendicular relative polarizations of the laser pump and X-ray probe pulse to isolate the contributions of electronic spin state transition, bond breaking, and heme macrocycle nuclear relaxation. The time evolution of the iron K-edge X-ray absorption near edge structure (XANES) features along the two major photochemical reaction coordinates, i.e., the iron(II)-CO bond elongation and the heme macrocycle doming relaxation were modeled by time-dependent density functional theory calculations. Combined results from the experiments and computations reveal insight into interplays between the nuclear and electronic structural dynamics along the CO photodissociation trajectory. Time-resolved small-angle X-ray scattering data during the same process are also simultaneously collected, which show that the local CO dissociation causes a protein quake propagating on different spatial and temporal scales. These studies are important for understanding gas transport and protein deligation processes and shed light on the interplay of active site conformational changes and large-scale protein reorganization.

Short-lived metal-centered excited state initiates iron-methionine photodissociation in ferrous cytochrome c

The dynamics of photodissociation and recombination in heme proteins represent an archetypical photochemical reaction widely used to understand the interplay between chemical dynamics and reaction environment. We report a study of the photodissociation mechanism for the Fe(II)-S bond between the heme iron and methionine sulfur of ferrous cytochrome c. This bond dissociation is an essential step in the conversion of cytochrome c from an electron transfer protein to a peroxidase enzyme. We use ultrafast X-ray solution scattering to follow the dynamics of Fe(II)-S bond dissociation and 1s3p (Kβ) X-ray emission spectroscopy to follow the dynamics of the iron charge and spin multiplicity during bond dissociation. From these measurements, we conclude that the formation of a triplet metal-centered excited state with anti-bonding Fe(II)-S interactions triggers the bond dissociation and precedes the formation of the metastable Fe high-spin quintet state.

Low-Frequency Protein Motions Coupled to Catalytic Sites

This review examines low-frequency vibrational modes of proteins and their coupling to enzyme catalytic sites. That protein motions are critical to enzyme function is clear, but the kinds of motions present in proteins and how they are involved in function remain unclear. Several models of enzyme-catalyzed reaction suggest that protein dynamics may be involved in the chemical step of the catalyzed reaction, but the evidence in support of such models is indirect. Spectroscopic studies of low-frequency protein vibrations consistently show that there are underdamped modes of the protein with frequencies in the tens of wavenumbers where overdamped behavior would be expected. Recent studies even show that such underdamped vibrations modulate enzyme active sites. These observations suggest that increasingly sophisticated spectroscopic methods will be able to unravel the link between low-frequency protein vibrations and enzyme function.

Energy Transfer across Nonpolar and Polar Contacts in Proteins: Role of Contact Fluctuations

Molecular dynamics simulations of the villin headpiece subdomain HP36 have been carried out to examine relations between rates of vibrational energy transfer across non-covalently bonded contacts and equilibrium structural fluctuations, with focus on van der Waals contacts. Rates of energy transfer across van der Waals contacts vary inversely with the variance of the contact length, with the same constant of proportionality for all nonpolar contacts of HP36. A similar relation is observed for hydrogen bonds, but the proportionality depends on contact pairs, with hydrogen bonds stabilizing the α-helices all exhibiting the same constant of proportionality, one that is distinct from those computed for other polar contacts. Rates of energy transfer across van der Waals contacts are found to be up to 2 orders of magnitude smaller than rates of energy transfer across polar contacts.

Role of atomic contacts in vibrational energy transfer in myoglobin

Heme proteins are ideal systems to investigate vibrational energy flow at the atomic level. Upon photoexcitation, a large amount of excess vibrational energy is selectively deposited on heme due to extremely fast internal conversion. This excess energy is redistributed to the surrounding protein moiety and then to water. Vibrational energy flow in myoglobin (Mb) was examined using picosecond time-resolved anti-Stokes ultraviolet resonance Raman (UVRR) spectroscopy. We used the Trp residue directly contacting the heme group as a selective probe for vibrationally excited populations. Trp residues were placed at different position close to the heme by site-directed mutagenesis. This technique allows us to monitor the excess energy on residue-to-residue basis. Anti-Stokes UVRR measurements for Mb mutants suggest that the dominant channel for energy transfer in Mb is the pathway through atomic contacts between heme and nearby amino acid residues as well as that between the protein and solvent water. It is found that energy flow through proteins is analogous to collisional exchange processes in solutions.

Ultrafast Structural Changes Decomposed from Serial Crystallographic Data

Direct visualization of electronic and molecular events during biochemical reactions is essential to mechanistic insights. This Letter presents an in-depth analysis of the serial crystallographic data sets collected by Barends and Schlichting et al. ( Science 2015 , 350 , 445 ) that probe the ligand photodissociation in carbonmonoxy myoglobin. This analysis reveals electron density changes caused by the formation of high-spin 3d atomic orbitals of the heme iron upon photolysis and their dynamic behaviors within the first few picoseconds. The heme iron is found popping out of and recoiling back into the heme plane in succession. These findings provide long-awaited visual validations for previous works using ultrafast spectroscopy and molecular dynamics simulations. Electron density variations are also found largely in the solvent during the first period of a low-frequency oscillation. This work demonstrates the importance of the analytical methods in detecting and isolating weak, transient signals of electronic changes arising from chemical reactions.

Topology, landscapes, and biomolecular energy transport

While ubiquitous, energy redistribution remains a poorly understood facet of the nonequilibrium thermodynamics of biomolecules. At the molecular level, finite-size effects, pronounced nonlinearities, and ballistic processes produce behavior that diverges from the macroscale. Here, we show that transient thermal transport reflects macromolecular energy landscape architecture through the topological characteristics of molecular contacts and the nonlinear processes that mediate dynamics. While the former determines transport pathways via pairwise interactions, the latter reflects frustration within the landscape for local conformational rearrangements. Unlike transport through small-molecule systems, such as alkanes, nonlinearity dominates over coherent processes at even quite short time- and length-scales. Our exhaustive all-atom simulations and novel local-in-time and space analysis, applicable to both theory and experiment, permit dissection of energy migration in biomolecules. The approach demonstrates that vibrational energy transport can probe otherwise inaccessible aspects of macromolecular dynamics and interactions that underly biological function.

Understanding vibrational energy transfer in macromolecules has been challenging to both theory and experiment. Here the authors use non-equilibrium molecular dynamics to reveal the relationship between heat transport in a model peptide, emergent nonlinearity, and the underlying free energy landscape.

Temperature-jump solution X-ray scattering reveals distinct motions in a dynamic enzyme

Correlated motions of proteins are critical to function, but these features are difficult to resolve using traditional structure determination techniques. Time-resolved X-ray methods hold promise for addressing this challenge, but have relied on the exploitation of exotic protein photoactivity, and are therefore not generalizable. Temperature jumps, through thermal excitation of the solvent, have been utilized to study protein dynamics using spectroscopic techniques, but their implementation in X-ray scattering experiments has been limited. Here, we perform temperature-jump small- and wide-angle X-ray scattering measurements on a dynamic enzyme, cyclophilin A, demonstrating that these experiments are able to capture functional intramolecular protein dynamics on the microsecond timescale. We show that cyclophilin A displays rich dynamics following a temperature jump, and use the resulting time-resolved signal to assess the kinetics of conformational changes. Two relaxation processes are resolved: a fast process is related to surface loop motions, and a slower process is related to motions in the core of the protein that are critical for catalytic turnover.

Atomistic characterization of collective protein–water–membrane dynamics

Water mediates correlated vibrations of atoms of protein and membrane bilayer surfaces.

Thermodynamic Scaling of Interfering Hemoglobin Strain Field Waves

Hemoglobin (Hgb) forms tetramers (dimerized α-β dimers), which enhance its globular stability and may also facilitate small gas molecule transport, as shown by recent all-atom Newtonian solvated simulations. Hydropathic bioinformatic thermodynamic scaling enables close comparisons of hemoglobin dimers with myoglobin and neuroglobin, and reveals many nonlocal wave-like features of strained Hgb structures at the coarse-grained amino acid level. The thermodynamic analysis employs two hydropathic scales, one describing abrupt first-order unfolding transitions, the other continuous second-order transitions. Small molecule exchange at hemes is a first-order process. Wave-like collective tetrameric features appropriate to ligand absorption and release, seen in optical experiments (short times), are identified thermodynamically at long times. Strain fields localized near hemes interfere with extended strain fields associated with dimer interfacial misfit, resulting in novel wavelength dependent dimer correlation function Fano antiresonances.

Dynamic multiple-scattering treatment of X-ray absorption: Parameterization of a new molecular dynamics force field for myoglobin

We present a detailed analysis of the X-ray absorption near-edge structure (XANES) data on the Fe K-edge of CO Myoglobin based on a combined procedure of Molecular Dynamics (MD) calculations and MXAN (Minuit XANes) data analysis that we call D-MXAN. The ability of performing quantitative XANES data analysis allows us to refine classical force field MD parameters, thus obtaining a reliable tool for the atomic investigation of this important model system for biological macromolecules. The iterative procedure here applied corrects the greatest part of the structural discrepancy between classical MD sampling and experimental determinations. Our procedure, moreover, is able to discriminate between different heme conformational basins visited during the MD simulation, thus demonstrating the necessity of a sampling on the order of tens of nanoseconds, even for an application such X-ray absorption spectroscopy data analysis.

Retinal isomerization in bacteriorhodopsin captured by a femtosecond x-ray laser

Ultrafast isomerization of retinal is the primary step in photoresponsive biological functions including vision in humans and ion transport across bacterial membranes. We used an x-ray laser to study the subpicosecond structural dynamics of retinal isomerization in the light-driven proton pump bacteriorhodopsin. A series of structural snapshots with near-atomic spatial resolution and temporal resolution in the femtosecond regime show how the excited all-trans retinal samples conformational states within the protein binding pocket before passing through a twisted geometry and emerging in the 13-cis conformation. Our findings suggest ultrafast collective motions of aspartic acid residues and functional water molecules in the proximity of the retinal Schiff base as a key facet of this stereoselective and efficient photochemical reaction.

Probing molecular forces in multi-component physiological membranes

Biological membranes are remarkably heterogeneous, composed of diverse lipid mixtures with distinct chemical structure and composition. By combining molecular dynamics simulations and the newly developed Lipid-Force Distribution Analysis (L-FDA), we explore force transmission in complex multi-component membrane models mimicking eukaryotic organelles. We found that the chemical-moiety based segmentation at membrane interfaces revealed a distinctive distribution of bonded and non-bonded forces in diverse membrane environment. Our molecular stress analysis could have far-reaching implications in describing the relationship between membrane mechanical properties and functional states of chemically distinct lipids.

Stress Propagation through Biological Lipid Bilayers in Silico

Membrane tension plays various critical roles in the cell. We here asked how fast and how far localized pulses of mechanical stress dynamically propagate through biological lipid bilayers. In both coarse-grained and all-atom molecular dynamics simulations of a dipalmitoylphosphatidylcholine lipid bilayer, we observed nanometer-wide stress pulses, propagating very efficiently longitudinally at a velocity of approximately 1.4 ± 0.5 nm/ps (km/s), in close agreement with the expected speed of sound from experiments. Remarkably, the predicted characteristic attenuation time of the pulses was in the order of tens of picoseconds, implying longitudinal stress propagation over length scales up to several tens of nanometers before damping. Furthermore, the computed dispersion relation leading to such damping was consistent with proposed continuum viscoelastic models of propagation. We suggest this mode of stress propagation as a potential ultrafast mechanism of signaling that may quickly couple mechanosensitive elements in crowded biological membranes.

Progress in small-angle scattering from biological solutions at high-brilliance synchrotrons

Small-angle X-ray scattering (SAXS) is an established technique that provides low-resolution structural information on macromolecular solutions. Recent decades have witnessed significant progress in both experimental facilities and in novel data-analysis approaches, making SAXS a mainstream method for structural biology. The technique is routinely applied to directly reconstruct low-resolution shapes of proteins and to generate atomistic models of macromolecular assemblies using hybrid approaches. Very importantly, SAXS is capable of yielding structural information on systems with size and conformational polydispersity, including highly flexible objects. In addition, utilizing high-flux synchrotron facilities, time-resolved SAXS allows analysis of kinetic processes over time ranges from microseconds to hours. Dedicated bioSAXS beamlines now offer fully automated data-collection and analysis pipelines, where analysis and modelling is conducted on the fly. This enables SAXS to be employed as a high-throughput method to rapidly screen various sample conditions and additives. The growing SAXS user community is supported by developments in data and model archiving and quality criteria. This review illustrates the latest developments in SAXS, in particular highlighting time-resolved applications aimed at flexible and evolving systems.

Reaction Path Averaging: Characterizing the Structural Response of the DNA Double Helix to Electron Transfer

A polarizable environment, prominently the solvent, responds to electronic changes in biomolecules rapidly. The knowledge of conformational relaxation of the biomolecule itself, however, may be scarce or missing. In this work, we describe in detail the structural changes in DNA undergoing electron transfer between two adjacent nucleobases. We employ an approach based on averaging of tens to hundreds of thousands of nonequilibrium trajectories generated with molecular dynamics simulation, and a reduction of dimensionality suitable for DNA. We show that the conformational response of the DNA proceeds along a single collective coordinate that represents the relative orientation of two consecutive base pairs, namely, a combination of helical parameters shift and tilt. The structure of DNA relaxes on time scales reaching nanoseconds, contributing marginally to the relaxation of energies, which is dominated by the modes of motion of the aqueous solvent. The concept of reaction path averaging (RPA), conveniently exploited in this context, makes it possible to filter out any undesirable noise from the nonequilibrium data, and is applicable to any chemical process in general.