Picosecond phase grating spectroscopy of hemoglobin and myoglobin: energetics and dynamics of global protein motion

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.

Sub-100-ps structural dynamics of horse heart myoglobin probed by time-resolved X-ray solution scattering

Here we report sub-100-ps structural dynamics of horse heart myoglobin revealed by time-resolved X-ray solution scattering. By applying the time-slicing scheme to the measurement and subsequent deconvolution, we investigate the protein structural dynamics that occur faster than the X-ray temporal pulse width of synchrotrons (~100 ps). The singular value decomposition analysis of the experimental data suggests that two structurally distinguishable intermediates are formed within 100 ps. In particular, the global structural change occurring on the time scale of 70 ps is identified.

Conformational Substates of Myoglobin Intermediate Resolved by Picosecond X-ray Solution Scattering

Conformational substates of proteins are generally considered to play important roles in regulating protein functions, but an understanding of how they influence the structural dynamics and functions of the proteins has been elusive. Here, we investigate the structural dynamics of sperm whale myoglobin associated with the conformational substates using picosecond X-ray solution scattering. By applying kinetic analysis considering all of the plausible candidate models, we establish a kinetic model for the entire cycle of the protein transition in a wide time range from 100 ps to 10 ms. Four structurally distinct intermediates are formed during the cycle, and most importantly, the transition from the first intermediate to the second one (B → C) occurs biphasically. We attribute the biphasic kinetics to the involvement of two conformational substates of the first intermediate, which are generated by the interplay between the distal histidine and the photodissociated CO.

Anisotropic Picosecond X-ray Solution Scattering from Photo-selectively Aligned Protein Molecules

Anisotropic X-ray scattering patterns of transiently aligned protein molecules in solution are measured by using pump-probe X-ray solution scattering. When a linearly polarized laser pulse interacts with an ensemble of molecules, the population of excited molecules is created with their transition dipoles preferentially aligned along the laser polarization direction. We measured the X-ray scattering from the myoglobin protein molecules excited by a linearly polarized, short laser pulse and obtained anisotropic scattering patterns on 100 ps time scale. An anisotropic scattering pattern contains higher structural information content than a typical isotropic pattern available from randomly oriented molecules. In addition, multiple independent diffraction patterns measured by using various laser polarization orientations will give substantially increased amount of structural information compared with a single isotropic pattern. By monitoring the temporal change of the anisotropic scattering pattern from 100 ps to 1 μs, we observed the orientational dynamics of photo-generated myoglobin with the rotational diffusion time of ∼15 ns.

Direct observation of myoglobin structural dynamics from 100 picoseconds to 1 microsecond with picosecond X-ray solution scattering

Here we report structural dynamics of equine myoglobin (Mb) in response to the CO photodissociation visualized by picosecond time-resolved X-ray solution scattering. The data clearly reveal new structural dynamics that occur in the timescale of ∼360 picoseconds (ps) and ∼9 nanoseconds (ns), which have not been clearly detected in previous studies.

Dynamics of conformational changes of Arabidopsis phototropin 1 LOV2 with the linker domain

Conformational changes of Arabidopsis phot1-LOV2 with the linker (phot1-LOV2-linker) were investigated from the viewpoint of the changes in molecular volume and molecular diffusion coefficient (D) by time-resolved transient grating (TG) and transient lens (TrL) methods. Although the absorption spectrum change completes within a few microseconds, the D-value detected by the TG method decreased drastically with a time constant of 1.0 ms from 9.2(+/-0.4)x10(-11) m(2)/s to 5.0(+/-0.3)x10(-11) m(2)/s. This time-dependent D was interpreted in terms of the unfolding of alpha-helices in the linker region. The change of the alpha-helices was confirmed by observing the recovery of the circular dichroism intensity. The TrL signal showed that the molecular volume decreases with two time constants; 300 micros and 1.0 ms. The former time constant is close to the previously observed photo-dissociation reaction rate of the phot1-LOV2 (without the linker) dimer, and the latter one agrees well with the rate of the D-change. Considering a similar time constant of the dissociation reaction of the LOV2 dimer, we interpreted these kinetics in terms of the dissociation step of the linker region from the LOV2 domain (T(390)(pre) state). After this step, the protein volume and D are decreased significantly with the lifetime of 1.0 ms. The D decrease indicates the increase of the intermolecular interaction between the protein and water molecules. On the basis of these observations, a two-step mechanism of the linker unfolding is proposed.

The Escape Process of Carbon Monoxide from Myoglobin to Solution at Physiological Temperature

The carbon monoxide (CO) docking sites involved in the ligand escape process from the iron atom in hem of myoglobin (Mb) to solution at physiological temperature were studied on the basis of the effect of xenon (Xe) on the ligand escape rate by the transient grating (TG) technique. The TG method provides a direct measurement of the changes in molecular volume. The apparent CO escaping rate and the volume contraction increase with increasing Xe pressure. The pressure dependence of the rate is consistent with that of the Xe population at the Xe(1) site. This result clearly shows that CO is trapped at the Xe(1) site before escaping to solvent in a Xe-free solution at room temperature. It is shown that only CO but not the trapped Xe is released by the photoexcitation of the Xe-trapped MbCO. A dissociation scheme is proposed to explain the enhancement of the escaping rate by the presence of Xe(1). There are two branches for the CO escaping pathway. The dominant part of the dissociated CO escapes to the solvent through the Xe(1) trapping site under the Xe-free condition, and there are at least three intermediate states along this pathway. When a Xe atom blocks the Xe(1) site, the CO escapes through another route.

Probing protein dynamics using temperature jump relaxation spectroscopy

There have been recent advances in initiating and perturbing chemical reactions on very fast timescales, as short as picoseconds, thus making it feasible to study a vast range of chemical kinetics problems that heretofore could not be studied. One such approach is the rapid heating of water solutions using laser excitation. Laser-induced temperature jump relaxation spectroscopy can be used to determine the dynamics of protein motion, an area largely unstudied for want of suitable experimental and theoretical probes, despite the obvious importance of dynamics to protein function. Coupled with suitable spectroscopic probes of structure, relaxation spectroscopy can follow the motion of protein atoms over an enormous time range, from picoseconds to minutes (or longer), and with substantial structural specificity.

A spectrally silent transformation in the photolysis of octopus rhodopsin: a protein conformational change without any accompanying change of the chromophore's absorption

A spectrally silent transformation in the photolysis of octopus rhodopsin was detected by the time-resolved transient grating method. Our results showed that at least two photointermediates, which share the same chromophore absorption spectrum, exist after the final absorption changes. Previously, mesorhodopsin was thought to decay to the final photoproduct, acid metarhodopsin with a lifetime of 38 micros at 15 degrees C, but the present results show that there is at least one intermediate species (called transient acid metarhodopsin) with a lifetime of 180 micros at 15 degrees C, before forming acid metarhodopsin. This indicates that the parts of the protein distant from the chromophore are still changing even after the changes in microenvironment around the chromophore are over. From the signal intensity detected by the transient grating method, the volume change of the spectrally silent transformation was found to be DeltaV = 13 ml/mol. The activation energy of the spectrally silent transformation is much lower than those of other transformations of octopus rhodopsin. Since stable acid metarhodopsin has not been shown to activate the G protein, this transient acid metarhodopsin may be responsible for G protein activation.

Volume and enthalpy profiles of CO binding to Fe(II) tetrakis-(4-sulfonatophenyl)porphyrin

The focus of this study is to examine volume and enthalpy profiles of ligand binding associated with CO-Fe(II) tetrakis-(4-sulfonato phenyl)-porphyrin (COFe(II)4SP) in aqueous solution. Temperature dependent photothermal beam deflection was employed to probe the overall enthalpy and volume changes associated with CO-photolysis and recombination. The analysis demonstrates that ligand recombination occurs with a pseudo first order rate constant of (2.5+/-0.2)x10(4) s(-1) (at 25 degrees C) with a corresponding volume decrease of 6+/-1 ml/mol. The activation enthalpy (DeltaH(double dagger)) and volume (DeltaV(double dagger)) change for CO recombination (determined from temperature/pressure dependent transient absorption spectroscopy) are found to be 3.9 kcal/mol and 8.2 ml/mol, respectively. These data are consistent with a mechanism in which photolysis yields a five-coordinate high spin (H(2)O)Fe(II)4SP complex that recombines in a single step to form the low spin (CO)(H(2)O)Fe(II)4SP complex. Base elimination, often associated with CO photolysis from hemes, is not observed in this system. The overall volume changes suggest a transition state with significant high spin character. Furthermore, these results demonstrate the utility of coupling photothermal techniques with variable pressure/temperature transient absorption spectroscopy to probe heme reaction dynamics.

Diffractive optics-based heterodyne-detected four-wave mixing signals of protein motion: from "protein quakes" to ligand escape for myoglobin

Ligand transport through myoglobin (Mb) has been observed by using optically heterodyne-detected transient grating spectroscopy. Experimental implementation using diffractive optics has provided unprecedented sensitivity for the study of protein motions by enabling the passive phase locking of the four beams that constitute the experiment, and an unambiguous separation of the Real and Imaginary parts of the signal. Ligand photodissociation of carboxymyoglobin (MbCO) induces a sequence of events involving the relaxation of the protein structure to accommodate ligand escape. These motions show up in the Real part of the signal. The ligand (CO) transport process involves an initial, small amplitude, change in volume, reflecting the transit time of the ligand through the protein, followed by a significantly larger volume change with ligand escape to the surrounding water. The latter process is well described by a single exponential process of 725 +/- 15 ns. at room temperature. The overall dynamics provide a distinctive signature that can be understood in the context of segmental protein fluctuations that aid ligand escape via a few specific cavities, and they suggest the existence of discrete escape pathways.

Dynamics of structure and energy of horse carboxymyoglobin after photodissociation of carbon monoxide

The energetics and structural volume changes after photodissociation of carboxymyoglobin are quantitatively investigated by laser-induced transient grating (TG) and photoacoustic calorimetric techniques. Various origins of the TG signal are distinguished: the phase grating signals due to temperature change, due to absorption spectrum change, and due to volume change. We found a new kinetics of approximately 700 ns (at room temperature), which was not observed by the flash photolysis technique. This kinetics should be attributed to the intermediate between the geminate pair and the fully dissociated state. The enthalpy of an intermediate species is determined to be 61 +/- 10 kJ/mol, which is smaller than the expected Fe-CO bond energy. The volume of MbCO slightly contracts (5 +/- 3 cm(3)/mol) during this process. CO is fully released from the protein by an exponential kinetics from 25 to -2 degrees C. During this escaping process, the volume expands by 14.7 +/- 2 cm(3)/mol at room temperature and 14 +/- 10 kJ/mol is released, which should represent the protein relaxation and the solvation of the CO (the enthalpy of this final state is 47 +/- 10 kJ/mol). A potential barrier between the intermediate and the fully dissociated state is DeltaH(*) = 41.3 kJ/mol and DeltaS(*) = 13.6 J mol(-1) K(-1). The TG experiment under a high wavenumber reveals that the volume expansion depends on the temperature from 25 to -2 degrees C. The volume changes and the energies of the intermediate species are discussed.

A multifrequency electron spin resonance study of T4 lysozyme dynamics

Electron spin resonance (ESR) spectroscopy at 250 GHz and 9 GHz is utilized to study the dynamics and local structural ordering of a nitroxide-labeled enzyme, T4 lysozyme (EC 3.2.1.17), in aqueous solution from 10 degrees C to 35 degrees C. Two separate derivatives, labeled at sites 44 and 69, were analyzed. The 250-GHz ESR spectra are well described by a microscopic ordering with macroscopic disordering (MOMD) model, which includes the influence of the tether connecting the probe to the protein. In the faster "time scale" of the 250-GHz ESR experiment, the overall rotational diffusion rate of the enzyme is too slow to significantly affect the spectrum, whereas for the 9-GHz ESR spectra, the overall rotational diffusion must be accounted for in the analysis. This is accomplished by using a slowly relaxing local structure model (SRLS) for the dynamics, wherein the tether motion and the overall motion are both included. In this way a simultaneous fit is successfully obtained for both the 250-GHz and 9-GHz ESR spectra. Two distinct motional/ordering modes of the probe are found for both lysozyme derivatives, indicating that the tether exists in two distinct conformations on the ESR time scale. The probe diffuses more rapidly about an axis perpendicular to its tether, which may result from fluctuations of the peptide backbone at the point of attachment of the spin probe.

Fast events in protein folding: the time evolution of primary processes

Most experimental studies on the dynamics of protein folding have been confined to timescales of 1 ms and longer. Yet it is obvious that many phenomena that are obligatory elements of the folding process occur on much faster timescales. For example, it is also now clear that the formation of secondary and tertiary structures can occur on nanosecond and microsecond times, respectively. Although fast events are essential to, and sometimes dominate, the overall folding process, with a few exceptions their experimental study has become possible only recently with the development of appropriate techniques. This review discusses new approaches that are capable of initiating and monitoring the fast events in protein folding with temporal resolution down to picoseconds. The first important results from those techniques, which have been obtained for the folding of some globular proteins and polypeptide models, are also discussed.

Protein dynamics derived from clusters of crystal structures

A method is presented to mathematically extract concerted structural transitions in proteins from collections of crystal structures. The "essential dynamics" procedure is used to filter out small-amplitude fluctuations from such a set of structures; the remaining large conformational changes describe motions such as those important for the uptake/release of substrate/ligand and in catalytic reactions. The method is applied to sets of x-ray structures for a number of proteins, and the results are compared with the results from essential dynamics as applied to molecular dynamics simulations of those proteins. A significant degree of similarity is found, thereby providing a direct experimental basis for the application of such simulations to the description of large concerted motions in proteins.

Conformational relaxation and ligand binding in myoglobin

Absorption spectroscopy with nanosecond time resolution shows that myoglobin undergoes conformational relaxation on the same time scale as geminate rebinding of carbon monoxide. Ligand rebinding following photodissociation of the heme-CO complex was measured from the amplitude of the average difference spectrum, while conformational changes were measured from changes in the detailed shape of the Soret spectra of the deoxyhemes. Experiments in which the solvent viscosity was varied between 1 and 300 cP and the temperature between 268 and 308 K were analyzed by fitting the multiwavelength kinetic data with both empirical and molecular models. Novel numerical techniques were employed in fitting the data, including the use of singular value decomposition to remove the effects of temperature and solvent on the spectra and of a Monte Carlo method to overcome the multiple minimum problem in searching parameter space. The molecular model is the minimal model that incorporates all of the major features of myoglobin kinetics at ambient temperatures, including a fast and slow rebinding conformation and two geminate states for each conformation. The results of fitting the kinetic data with this model indicate that the geminate-rebinding rates for the two conformations differ by at least a factor of 100. The differences between the spectra of the two conformations generated from the fits are similar to the differences between those of the R and T conformations of hemoglobin. In modeling the data, the dependence of the rates on temperature and viscosity was parametrized using a modification of Kramers theory which includes the contributions of both protein and solvent to the friction. The rate of the transition from the fast to the slow rebinding conformation is found to be inversely proportional to the viscosity when the viscosity exceeds about 30 cP and nearly viscosity independent at low viscosity. The viscosity dependence at high viscosities suggests that the two conformations differ by the global displacement of protein atoms on the proximal side of the heme observed by X-ray crystallography. We suggest that the conformational change observed in our experiments corresponds to the final portion of the nonexponential conformational relaxation recently observed by Anfinrud and co-workers, which begins on a picosecond time scale. Furthermore, extrapolation of our data to temperatures near that of the solvent glass transition suggests that this conformational relaxation may very well be the one postulated by Frauenfelder and co-workers to explain the decrease in the rate of geminate rebinding with increasing temperature above 180 K.

A connected-cluster of hydration around myoglobin: correlation between molecular dynamics simulations and experiment

An analysis of a molecular dynamics simulation of metmyoglobin in an explicit solvent environment of 3,128 water molecules has been performed. Both statics and dynamics of the protein-solvent interface are addressed in a comparison with experiment. Three-dimensional density distributions, temperature factors, and occupancy weights are computed for the solvent by using the trajectory coordinates. Analysis of the hydration leads to the localization of more than 500 hydration sites distributed into multiple layers of solvation located between 2.6 and 6.8 A from the atomic protein surface. After locating the local solvent density maxima or hydration sites we conclude that water molecules of hydration positions and hydration sites are distinct concepts. Both global and detailed properties of the hydration cluster around myoglobin are compared with recent neutron and X-ray data on myoglobin. Questions arising from differences between X-ray and neutron data concerning the locations of the protein-bound water are investigated. Analysis of water site differences found from X-ray and neutron experiments compared with our simulation shows that the simulation gives a way to unify the hydration picture given by the two experiments.

A photoacoustic calorimetry study of horse carboxymyoglobin on the 10-nanosecond time scale

The development of a photoacoustic calorimeter with a time resolution of 10 ns is presented, and the dynamics of the enthalpy and volume changes found in the photodissociation of CO from horse carboxymyoglobin are examined. With this enhanced time resolution a new transient species, the lifetime of which is 29 ns at 20 degrees C, is observed in the ligand dissociation process.