Complex nonexponential relaxation in myoglobin after photodissociation of MbCO: measurement and analysis from 2 ps to 56 υs

Prediction of Allosteric Communication Pathways in Proteins

MOTIVATION

Allostery in proteins is an essential phenomenon in biological processes. In this paper, we present a computational model to predict paths of maximum information transfer between active and allosteric sites. In this information theoretic study, we use mutual information as the measure of information transfer, where transition probability of information from one residue to its contacting neighbors is proportional to the magnitude of mutual information between the two residues. Starting from a given residue and using a Hidden Markov Model, we successively determine the neighboring residues that eventually lead to a path of optimum information transfer. The Gaussian approximation of mutual information between residue pairs is adopted. The limits of validity of this approximation are discussed in terms of a nonlinear theory of mutual information and its reduction to the Gaussian form.

RESULTS

Predictions of the model are tested on six widely studied cases, CheY Bacterial Chemotaxis, B-cell Lymphoma extra-large Bcl-xL, Human proline isomerase cyclophilin A (CypA), Dihydrofolate reductase DHFR, HRas GTPase, and Caspase-1. The communication transmission rendering the propagation of local fluctuations from the active sites throughout the structure in multiple paths correlate well with the known experimental data. Distinct paths originating from the active site may likely represent a multi functionality such as involving more than one allosteric site and/or preexistence of some other functional states. Our model is computationally fast and simple, and can give allosteric communication pathways, which are crucial for the understanding and control of protein functionality.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Information Flow and Allosteric Communication in Proteins

Based on Schreiber's work on transfer entropy, a molecular theory of nonlinear information transfer in proteins is developed. The joint distribution function for residue fluctuations is expressed in terms of tensor Hermite polynomials which conveniently separate harmonic and nonlinear contributions to information transfer. The harmonic part of information transfer is expressed as the difference between time dependent and independent mutual information. Third order nonlinearities are discussed in detail. Amount and speed of information transfer between residues, important for understanding allosteric activity in proteins, are discussed. While mutual information shows the maximum amount of information that may be transferred between two residues, it does not explain the actual amount of transfer nor the transfer rate of information. For this, dynamic equations of the system are needed. The solution of the Langevin equation and molecular dynamics trajectories are used in the present work for this purpose. Allosteric communication in Human NAD-dependent isocitrate dehydrogenase is studied as an example. Calculations show that several paths contribute collectively to information transfer. Important residues on these paths are identified. Time resolved information transfer between these residues, their amplitudes and transfer rates, which are in agreement with time resolved ultraviolet resonance Raman measurements in general, are estimated. Estimated transfer rates are in the order of 1-20 megabits per second. Information transfer from third order contributions are one to two orders of magnitude smaller than the harmonic terms, showing that harmonic analysis is a good approximation to information transfer.

Ruthenium hydrides encapsulated in sol–gel glasses exhibit new ultrafast vibrational dynamics

Vibrational dynamics were measured by IR pump–probe spectroscopy and two-dimensional IR spectroscopy for triruthenium dodecacarbonyl and the undecacarbonyl hydride that forms when it is encapsulated in an alumina sol–gel glass. For comparison, a triruthenium undecacarbonyl hydride salt was also synthesized and studied in neat solution to identify the potential influence of the confined solvent environment on the dynamics experienced by carbon monoxide ligands. The vibrational lifetime was found to be significantly decreased for both hydride species relative to the dodecacarbonyl compound. Conversely, spectral diffusion of the CO vibrations was measured to be faster for the parent compound. The most significant dynamic changes occurred upon transformation from the starting compound to the hydride, while only minor differences were observed between the dynamics of the freely dissolved and sol–gel encapsulated hydrides. The results suggest that the structural change to the hydride has the largest impact on the dynamics and that its improved catalytic properties likely do not originate from confined solvent effects.

Impact of Conformational Substates and Energy Landscapes on Understanding Hemoglobin Kinetics and Function

Hans Frauenfelder's discovery of conformational substates in studies of myoglobin carbon monoxide geminate rebinding kinetics at cryogenic temperatures (Austin RH, Beeson KW, Eisenstein L, Frauenfelder H, & Gunsalus IC (1975) Dynamics of Ligand Binding to Myoglobin. Biochemistry 14(24):5355-5373) followed by his introduction of energy landscape theory with Peter Wolynes (Frauenfelder H, Sligar SG, & Wolynes PG (1991) The Energy Landscapes and Motions of Proteins. Science 254(5038):1598-1603) marked the beginning of a new era in the physics and physical chemistry of proteins. Their work played a major role in demonstrating the power and importance of dynamics and of Kramers reaction rate theory for understanding protein function. The biggest impact of energy landscape theory has been in the protein folding field, which is well-known and has been documented in numerous articles and reviews, including a recent one of my own (Eaton WA (2021) Modern Kinetics and Mechanism of Protein Folding: a Retrospective. J. Phys. Chem. B. 125(14):3452-3467). Here I will describe the much less well-known impact of their modern view of proteins on both experimental and theoretical studies of hemoglobin kinetics and function. I will first describe how Frauenfelder's experiments motivated and influenced my own research on myoglobin, which were key ingredients to my work on understanding hemoglobin.

Temperature Dependence of Tryptophan Fluorescence Lifetime in Aqueous Glycerol and Trehalose Solutions

The temperature dependences of tryptophan fluorescence decay kinetics in aqueous glycerol and 1 M trehalose solutions were examined. The fluorescence decay kinetics were recorded in the spectral region of 292.5-417.5 nm with nanosecond time resolution. The kinetics curves were approximated by the sum of three exponential terms, and the spectral distribution (DAS) of these components was determined. An antisymbatic course of fluorescence decay times of two (fast and medium) components in the temperature range from -60 to +10°C was observed. The third (slow) component showed only slight temperature dependence. The antisymbatic behavior of fluorescence lifetimes of the fast and medium components was explained on the assumption that some of the excited tryptophan molecules are transferred from a short-wavelength B-form with short fluorescence lifetime to a long-wavelength R-form with an intermediate fluorescence lifetime. This transfer occurred in the indicated temperature range.

Enhanced vibrational solvatochromism and spectral diffusion by electron rich substituents on small molecule silanes

Fourier transform infrared and two-dimensional IR (2D-IR) spectroscopies were applied to two different silanes in three different solvents. The selected solutes exhibit different degrees of vibrational solvatochromism for the Si-H vibration. Density functional theory calculations confirm that this difference in sensitivity is the result of higher mode polarization with more electron withdrawing ligands. This mode sensitivity also affects the extent of spectral diffusion experienced by the silane vibration, offering a potential route to simultaneously optimize the sensitivity of vibrational probes in both steady-state and time-resolved measurements. Frequency-frequency correlation functions obtained by 2D-IR show that both solutes experience dynamics on similar time scales and are consistent with a picture in which weakly interacting solvents produce faster, more homogeneous fluctuations. Molecular dynamics simulations confirm that the frequency-frequency correlation function obtained by 2D-IR is sensitive to the presence of hydrogen bonding dynamics in the surrounding solvation shell.

Tryptophan-to-heme electron transfer in ferrous myoglobins

It was recently demonstrated that in ferric myoglobins (Mb) the fluorescence quenching of the photoexcited tryptophan 14 (*Trp(14)) residue is in part due to an electron transfer to the heme porphyrin (porph), turning it to the ferrous state. However, the invariance of *Trp decay times in ferric and ferrous Mbs raises the question as to whether electron transfer may also be operative in the latter. Using UV pump/visible probe transient absorption, we show that this is indeed the case for deoxy-Mb. We observe that the reduction generates (with a yield of about 30%) a low-valence Fe-porphyrin π [Fe(II)(porph(●-))] -anion radical, which we observe for the first time to our knowledge under physiological conditions. We suggest that the pathway for the electron transfer proceeds via the leucine 69 (Leu(69)) and valine 68 (Val(68)) residues. The results on ferric Mbs and the present ones highlight the generality of Trp-porphyrin electron transfer in heme proteins.

Dynamics of geminate rebinding of CO to cytochrome c in guanidine HCl probed by femtosecond vibrational spectroscopy

Femtosecond vibrational spectroscopy was used to probe the rebinding dynamics of CO to cytochrome c (Cytc) in 1.8 and 7 M guanidine HCl (GdnHCl) after photodeligation of the corresponding CO-bound protein in D2O buffer (pD = 7.4) at 283 K. Geminate rebinding (GR) dynamics of CO to the folded Cytc in 1.8 M GdnHCl (nCytc) is similar to that to chemically modified cytochrome c (cCytc), suggesting that the overall conformations of nCytcCO and cCytcCO are similar. About 86% of the dissociated CO molecules were geminately rebound to nCytc nonexponentially within 1 ns. The efficient GR of CO to the folded Cytc can be attributed to the organized protein matrix near the active site of nCytc that provides an efficient trap for the diffusing CO ligand after photodissociation. Although the concentration of nCytc did not affect its GR yield of CO, GR yield of CO to the unfolded Cytc in 7 M GdnHCl (uCytc) increased from 5 to 30% as the protein concentration increased from 0.3 to 9 mM. Time-resolved spectra of the (13)CO dissociated from both 9 mM nCytc(13)CO and 9 mM uCytc(13)CO showed a growing band with a peak at 2090 cm(-1) on the picosecond time scale, which was assigned to (13)CO in D2O solvent. At 1 ns, the fraction of the CO band in the solvent was about 10% of the nascent photodeligated protein in nCytc and more than 50% in the concentrated uCytc. Whereas a small opening in the active site of nCytc is responsible for the ultrafast escape of CO to solution in the folded protein, a large fraction of the CO escape to the solvent in uCytc results from the denatured structure of the active site in the unfolded protein. The spectrum of the CO dissociated from the concentrated uCytcCO contained a band that decayed as efficiently as that for the folded protein, suggesting that some fraction of uCytcCO might form aggregates even in 7 M denaturant, such that the aggregate acts as an efficient trap for the diffusing CO after deligation. No hint of precipitate in the concentrated uCytcCO and protein refolding upon dilution of the GdnHCl indicate that the aggregate does not grow continuously but remains as a soluble oligomer. The delayed appearance of the solvated CO and the inefficient GR of CO in uCytcCO suggest that the monomeric unfolded CytcCO so loosely arranged that the protein matrix cannot trap CO efficiently but the bound CO is still buried within hydrophobic residues even under the harsh denaturation condition.

Proteins and Glasses

The structure, the energy landscape, and the dynamics of proteins and glasses are similar. Both types of systems display characteristic nonexponential time dependencies of relaxation phenomena. Experiments suggest that both, proteins and glasses, are heterogeneous and that this fact causes the observed time dependence. This result is discussed in terms of the rough energy landscape characteristic of complex systems.

Protein structural dynamics in solution unveiled via 100-ps time-resolved x-ray scattering

We have developed a time-resolved x-ray scattering diffractometer capable of probing structural dynamics of proteins in solution with 100-ps time resolution. This diffractometer, developed on the ID14B BioCARS (Consortium for Advanced Radiation Sources) beamline at the Advanced Photon Source, records x-ray scattering snapshots over a broad range of q spanning 0.02-2.5 A(-1), thereby providing simultaneous coverage of the small-angle x-ray scattering (SAXS) and wide-angle x-ray scattering (WAXS) regions. To demonstrate its capabilities, we have tracked structural changes in myoglobin as it undergoes a photolysis-induced transition from its carbon monoxy form (MbCO) to its deoxy form (Mb). Though the differences between the MbCO and Mb crystal structures are small (rmsd < 0.2 A), time-resolved x-ray scattering differences recorded over 8 decades of time from 100 ps to 10 ms are rich in structure, illustrating the sensitivity of this technique. A strong, negative-going feature in the SAXS region appears promptly and corresponds to a sudden > 22 A(3) volume expansion of the protein. The ensuing conformational relaxation causes the protein to contract to a volume approximately 2 A(3) larger than MbCO within approximately 10 ns. On the timescale for CO escape from the primary docking site, another change in the SAXS/WAXS fingerprint appears, demonstrating sensitivity to the location of the dissociated CO. Global analysis of the SAXS/WAXS patterns recovered time-independent scattering fingerprints for four intermediate states of Mb. These SAXS/WAXS fingerprints provide stringent constraints for putative models of conformational states and structural transitions between them.

A unified model of protein dynamics

Protein functions require conformational motions. We show here that the dominant conformational motions are slaved by the hydration shell and the bulk solvent. The protein contributes the structure necessary for function. We formulate a model that is based on experiments, insights from the physics of glass-forming liquids, and the concepts of a hierarchically organized energy landscape. To explore the effect of external fluctuations on protein dynamics, we measure the fluctuations in the bulk solvent and the hydration shell with broadband dielectric spectroscopy and compare them with internal fluctuations measured with the Mössbauer effect and neutron scattering. The result is clear. Large-scale protein motions are slaved to the fluctuations in the bulk solvent. They are controlled by the solvent viscosity, and are absent in a solid environment. Internal protein motions are slaved to the beta fluctuations of the hydration shell, are controlled by hydration, and are absent in a dehydrated protein. The model quantitatively predicts the rapid increase of the mean-square displacement above approximately 200 K, shows that the external beta fluctuations determine the temperature- and time-dependence of the passage of carbon monoxide through myoglobin, and explains the nonexponential time dependence of the protein relaxation after photodissociation.

Ultrafast dynamics of myoglobin without the distal histidine: stimulated vibrational echo experiments and molecular dynamics simulations

Ultrafast protein dynamics of the CO adduct of a myoglobin mutant with the polar distal histidine replaced by a nonpolar valine (H64V) have been investigated by spectrally resolved infrared stimulated vibrational echo experiments and molecular dynamics (MD) simulations. In aqueous solution at room temperature, the vibrational dephasing rate of CO in the mutant is reduced by approximately 50% relative to the native protein. This finding confirms that the dephasing of the CO vibration in the native protein is sensitive to the interaction between the ligand and the distal histidine. The stimulated vibrational echo observable is calculated from MD simulations of H64V within a model in which vibrational dephasing is driven by electrostatic forces. In agreement with experiment, calculated vibrational echoes show slower dephasing for the mutant than for the native protein. However, vibrational echoes calculated for H64V do not show the quantitative agreement with measurements demonstrated previously for the native protein.

Hydrophobic Distal Pocket Affects NO−Heme Geminate Recombination Dynamics in Dehaloperoxidase and H64V Myoglobin

The recombination dynamics of NO with dehaloperoxidase (DHP) from Amphitrite ornata following photolysis were measured by femtosecond time-resolved absorption spectroscopy. Singular value decomposition (SVD) analysis reveals two important basis spectra. The first SVD basis spectrum reports on the population of photolyzed NO molecules and has the appearance of the equilibrium difference spectrum between the deoxy and NO forms of DHP. The first basis time course has two kinetic components with time constants of tau(11) approximately 9 ps and tau(12) approximately 50 ps that correspond to geminate recombination. The fast geminate process tau(11) arises from a contact pair with the heme iron in a bound state with S = 3/2 spin. The slow geminate process tau(12) corresponds to the recombination from a more remote docking site >3 A from the heme iron with the greater barrier corresponding to a S = 5/2 spin state. The second SVD basis spectrum represents a time-dependent Soret band shift indicative of heme photophysical processes and protein relaxation with time constants of tau(21) approximately 3 ps and tau(22) approximately 17 ps, respectively. A comparison between the more rapid rate constant of the slow geminate phase in DHP-NO and horse heart myoglobin (HHMbNO) or sperm whale myoglobin (SWMbNO) suggests that protein interactions with photolyzed NO are weaker in DHP than in the wild-type MbNOs, consistent with the hydrophobic distal pocket of DHP. The slower protein relaxation rate tau(22) in DHP-NO relative to HHMbNO implies less effective trapping in the docking site of the distal pocket and is consistent with a greater yield for the fast geminate process. The trends observed for DHP-NO also hold for the H64V mutant of SWMb (H64V MbNO), consistent with a more hydrophobic distal pocket for that protein as well. We examine the influence of solution viscosity on NO recombination by varying the glycerol content in the range from 0% to 90% (v/v). The dominant effect of increasing viscosity is the increase of the rate of the slow geminate process, tau(12), coupled with a population decrease of the slow geminate component. Both phenomena are similar to the effect of viscosity on wild-type Mb due to slowing of protein relaxation resulting from an increased solution viscosity and protein surface dehydration.

Extended subnanosecond structural dynamics of myoglobin revealed by Laue crystallography

Work carried out over the last 30 years unveiled the role of structural dynamics in controlling protein function. Cavity networks modulate structural dynamics trajectories and are functionally relevant; in globins they have been assigned a role in ligand migration and docking. These findings raised renewed interest for time-resolved structural investigations of myoglobin (Mb), a simple heme protein displaying a photosensitive iron-ligand bond. Photodissociation of MbCO generates a nonequilibrium population of protein structures relaxing over a time range extending from picoseconds to milliseconds. This process triggers ligand migration to matrix cavities with clear-cut effects on the rate and yield of geminate rebinding. Here, we report subnanosecond time-resolved Laue diffraction data on the triple mutant YQR-Mb [Leu-29(B10)Tyr, His-64(E7)Gln, Thr-67(E10)Arg] that depict the sequence of structural events associated with heme and protein relaxation from 100 ps to 316 ns and above. The photodissociated ligand rapidly (<0.1 ns) populates the Xe-binding cavity distal to the heme. Moreover, the heme relaxation toward the deoxy configuration is heterogeneous, with a slower phase ( approximately ns) evident in these experiments. Damping of the heme response appears to result from a strain exerted by the E-helix via the CD-turn; Phe-43(CD1), in close contact with heme, opposes tilt until the strain is relieved. A comparison with crystallographic data on wild-type Mb and mutants Leu(29)Phe or Leu(29)Trp suggests that the internal structure controls the rate and amplitude of the relaxation events. A correlation between structural dynamics as unveiled by Laue crystallography and functional properties of Mb is presented.

Ligand migration pathway and protein dynamics in myoglobin: a time-resolved crystallographic study on L29W MbCO

By using time-resolved x-ray crystallography at room temperature, structural relaxations and ligand migration were examined in myoglobin (Mb) mutant L29W from nanoseconds to seconds after photodissociation of carbon monoxide (CO) from the heme iron by nanosecond laser pulses. The data were analyzed in terms of transient kinetics by fitting trial functions to integrated difference electron density values obtained from select structural moieties, thus allowing a quantitative description of the processes involved. The observed relaxations are linked to other investigations on protein dynamics. At the earliest times, the heme has already completely relaxed into its domed deoxy structure, and there is no photo-dissociated CO visible at the primary docking site. Initial relaxations of larger globin moieties are completed within several hundred nanoseconds. They influence the concomitant migration of photo-dissociated CO to the Xe1 site, where it appears at approximately 300 ns and leaves again at approximately 1.5 ms. The extremely long residence time in Xe1 as compared with wild-type MbCO implies that, in the latter protein, the CO exits the protein from Xe1 predominantly via the distal pocket. A well-defined deligated state is populated between approximately 2 micros and approximately 1 ms; its structure is very similar to the equilibrium deoxy structure. Between 1.5 and 20 ms, no CO is visible in the protein interior; it is either distributed among many sites within the protein or has escaped to the solvent. Finally, recombination at the heme iron occurs after >20 ms.

Observation of a power-law memory kernel for fluctuations within a single protein molecule

The fluctuation of the distance between a fluorescein-tyrosine pair within a single protein complex was directly monitored in real time by photoinduced electron transfer and found to be a stationary, time-reversible, and non-Markovian Gaussian process. Within the generalized Langevin equation formalism, we experimentally determine the memory kernel K(t), which is proportional to the autocorrelation function of the random fluctuating force. K(t) is a power-law decay, t(-0.51 +/- 0.07) in a broad range of time scales (10(-3)-10 s). Such a long-time memory effect could have implications for protein functions.

Dynamics of hemoglobin in human erythrocytes and in solution: influence of viscosity studied by ultrafast vibrational echo experiments

Ultrafast spectrally resolved stimulated vibrational echo experiments are used to measure the vibrational dephasing of the CO stretching mode of hemoglobin-CO (HbCO) inside living human erythrocytes (red blood cells), in liquid solutions, and in a glassy matrix. A method is presented to overcome the adverse impact on the vibrational echo signal from the strong light scattering caused by the cells. The results from the cytoplasmic HbCO are compared to experiments on aqueous HbCO samples prepared in different buffers, solutions containing low and high concentrations of glycerol, and in a solid trehalose matrix. Measurements are also presented that provide an accurate determination of the viscosity at the very high Hb concentration that is found inside the cells. It is demonstrated that the dynamics of the protein, as sensed by the CO ligand, are the same inside the erythrocytes and in aqueous solution and are independent of the viscosity. In solutions that are predominantly glycerol, the dynamics are modified somewhat but are still independent of viscosity. The experiments in trehalose give the dynamics at infinite viscosity and are used to separate the viscosity-dependent dynamics from the viscosity-independent dynamics. Although the HbCO dynamics are the same in the red blood cell and in the equivalent aqueous solutions, differences in the absorption spectra show that the distribution of a protein's equilibrium substates is sensitive to small pH differences.

Protein-ligand interaction probed by time-resolved crystallography

Time-resolved (TR) crystallography is a unique method for determining the structures of intermediates in biomolecular reactions. The technique reached its mature stage with the development of the powerful third-generation synchrotron X-ray sources, and the advances in data processing and analysis of time-resolved Laue crystallographic data. A time resolution of 100 ps has been achieved and relatively small structural changes can be detected even from only partial reaction initiation. The remaining challenge facing the application of this technique to a broad range of biological systems is to find an efficient and rapid, system-specific method for the reaction initiation in the crystal. Other frontiers for the technique involve the continued improvement in time resolution and further advances in methods for determining intermediate structures and reaction mechanisms. The time-resolved technique, combined with trapping methods and computational approaches, holds the promise for a complete structure-based description of biomolecular reactions.

Different relaxations in myoglobin after photolysis

To clarify the interplay of kinetic hole-burning (KHB), structural relaxation, and ligand migration in myoglobin (Mb), we measured time-resolved absorption spectra in the Soret region after photolysis of carbon monoxide Mb (MbCO) in the temperature interval 120-260 K and in the time window 350 ns to 200 ms. The spectral contributions of both photolyzed (Mb*) and liganded Mb (MbCO) have been analyzed by taking into account homogeneous bandwidth, coupling to vibrational modes, and static conformational heterogeneity. We succeeded in separating the "time-dependent" spectral changes, and this work provides possibilities to identify the events in the process of ligand rebinding. KHB is dominant at T <190 K in both the Mb* and the MbCO components. For MbCO, conformational substates interconversion at higher temperatures tends to average out the KHB effect. At 230-260 K, whereas almost no shift is observed in the MbCO spectrum, a shift of the order of approximately 80 cm(-1) is observed in Mb*. We attribute this shift to protein relaxation coupled to ligand migration. The time dependence of the Mb* spectral shift is interpreted with a model that enables us to calculate the highly nonexponential relaxation kinetics. Fits of stretched exponentials to this kinetics yield Kohlrausch parameter values of 0.25, confirming the analogy between proteins and glasses.

Orientational distribution of CO before and after photolysis of MbCO and HbCO: a determination using time-resolved polarized Mid-IR spectroscopy

The technique of time-resolved polarized mid-IR spectroscopy was used to probe the orientational distribution of carbon monoxide (CO) bound to and docked within horse myoglobin, sperm whale myoglobin, and human hemoglobin A in neutral pH solution at 283 K. An accurate determination of the orientation required that the experimentally measured polarization anisotropy be corrected for the effects of fractional photolysis in an optically thick sample. The experimental method measures the direction of the transition dipole, which is parallel to the CO bond axis when docked and nearly parallel when bound to the heme. The polarization anisotropy of bound CO is virtually the same for all protein systems investigated and is unchanging across its inhomogeneously broadened mid-IR absorption spectrum. From these results, it was concluded that the transition dipole moment of bound CO is oriented </=7 degrees from the heme plane normal. The polarized absorbance spectra of docked CO are similar for all protein systems investigated, but in stark contrast to bound CO, the polarization anisotropy is strongly correlated with vibrational frequency. The frequency-dependent anisotropy imposes severe constraints on the orientational probability distribution function of the transition dipole, which is well described as a dipole bathed in a Stark field whose out-of-plane motion is constrained by a simple double-well potential. The orientational and spatial constraints imposed on docked CO by the surrounding highly conserved amino acids serve to mediate ligand transport to and from the binding site and thereby control the rates and pathways for geminate ligand rebinding and ligand escape.

Protein Conformational Dynamics Probed by Single-Molecule Electron Transfer

Electron transfer is used as a probe for angstrom-scale structural changes in single protein molecules. In a flavin reductase, the fluorescence of flavin is quenched by a nearby tyrosine residue by means of photo-induced electron transfer. By probing the fluorescence lifetime of the single flavin on a photon-by-photon basis, we were able to observe the variation of flavin-tyrosine distance over time. We could then determine the potential of mean force between the flavin and the tyrosine, and a correlation analysis revealed conformational fluctuation at multiple time scales spanning from hundreds of microseconds to seconds. This phenomenon suggests the existence of multiple interconverting conformers related to the fluctuating catalytic reactivity.

Complex landscape of protein structural dynamics unveiled by nanosecond Laue crystallography

Although conformational changes are essential for the function of proteins, little is known about their structural dynamics at atomic level resolution. Myoglobin (Mb) is the paradigm to investigate conformational dynamics because it is a simple globular heme protein displaying a photosensitivity of the iron-ligand bond. Upon laser photodissociation of carboxymyoglobin Mb a nonequilibrium population of protein structures is generated that relaxes over a broad time range extending from picoseconds to milliseconds. This process is associated with migration of the ligand to cavities in the matrix and with a reduction in the geminate rebinding rate by several orders of magnitude. Here we report nanosecond time-resolved Laue diffraction data to 1.55-A resolution on a Mb mutant, which depicts the sequence of structural events associated with this extended relaxation. Motions of the distal E-helix, including the mutated residue Gln-64(E7), and of the CD-turn are found to lag significantly (100-300 ns) behind local rearrangements around the heme such as heme tilting, iron motion out of the heme plane, and swinging of the mutated residue Tyr-29(B10), all of which occur promptly (< or =3 ns). Over the same delayed time range, CO is observed to migrate from a cavity distal to the heme known to bind xenon (called Xe4) to another such cavity proximal to the heme (Xe1). We propose that the extended relaxation of the globin moiety reflects reequilibration among conformational substates known to play an essential role in controlling protein function.

Watching a protein as it functions with 150-ps time-resolved x-ray crystallography

We report picosecond time-resolved x-ray diffraction from the myoglobin (Mb) mutant in which Leu29 is replaced by Phe (L29Fmutant). The frame-by-frame structural evolution, resolved to 1.8 angstroms, allows one to literally "watch" the protein as it executes its function. Time-resolved mid-infrared spectroscopy of flash-photolyzed L29F MbCO revealed a short-lived CO intermediate whose 140-ps lifetime is shorter than that found in wild-type protein by a factor of 1000. The electron density maps of the protein unveil transient conformational changes far more dramatic than the structural differences between the carboxy and deoxy states and depict the correlated side-chain motion responsible for rapidly sweeping CO away from its primary docking site.

Slaving: solvent fluctuations dominate protein dynamics and functions

Protein motions are essential for function. Comparing protein processes with the dielectric fluctuations of the surrounding solvent shows that they fall into two classes: nonslaved and slaved. Nonslaved processes are independent of the solvent motions; their rates are determined by the protein conformation and vibrational dynamics. Slaved processes are tightly coupled to the solvent; their rates have approximately the same temperature dependence as the rate of the solvent fluctuations, but they are smaller. Because the temperature dependence is determined by the activation enthalpy, we propose that the solvent is responsible for the activation enthalpy, whereas the protein and the hydration shell control the activation entropy through the energy landscape. Bond formation is the prototype of nonslaved processes; opening and closing of channels are quintessential slaved motions. The prevalence of slaved motions highlights the importance of the environment in cells and membranes for the function of proteins.

Proximal ligand motions in H93G myoglobin

Resonance Raman spectroscopy has been used to observe changes in the iron-ligand stretching frequency in photoproduct spectra of the proximal cavity mutant of myoglobin H93G. The measurements compare the deoxy ferrous state of the heme iron in H93G(L), where L is an exogenous imidazole ligand bound in the proximal cavity, to the photolyzed intermediate of H93G(L)*CO at 8 ns. There are significant differences in the frequencies of the iron-ligand axial out-of-plane mode nu(Fe-L) in the photoproduct spectra depending on the nature of L for a series of methyl-substituted imidazoles. Further comparison was made with the proximal cavity mutant of myoglobin in the absence of exogenous ligand (H93G) and the photoproduct of the carbonmonoxy adduct of H93G (H93G-*CO). For this case, it has been shown that H2O is the axial (fifth) ligand to the heme iron in the deoxy form of H93G. The photoproduct of H93G-*CO is consistent with a transiently bound ligand proposed to be a histidine. The data presented here further substantiate the conclusion that a conformationally driven ligand switch exists in photolyzed H93G-*CO. The results suggest that ligand conformational changes in response to dynamic motions of the globin on the nanosecond and longer time scales are a general feature of the H93G proximal cavity mutant.

Heme charge-transfer band III is vibronically coupled to the Soret band

A complete resonance Raman excitation profile of the heme charge-transfer band known as band III is presented. The data obtained throughout the near-infrared region show preresonance with the Q-band, but the data also clearly show the enhancement of a number of modes in the spectral region of band III. Only nontotally symmetric modes are observed to have resonance enhancement in the band III region. The observed resonance enhancements in modes of B(1g) symmetry are compared with the enhancements of those same modes in the excitation profiles of the Q-band of deoxy myoglobin, also presented here for this first time. The Q-band data agree well with the theory of vibronic coupling in metalloporphyrins (Shelnutt, J. A. J. Chem. Phys. 1981, 74, 6644-6657). The strong vibronic coupling of the Q-band of the deoxy form of hemes is discussed in terms of the enhancement of modes with both B(1g) and A(2g) symmetry. The comparison between the Q-band and band III reveals that, consistent with the theory, only modes of B(1g) symmetry are enhanced in the vicinity of band III. These results show that band III is vibronically coupled to the Soret band. The coupling of band III to modes with strong rhombic distortion of the heme macrocycle calls into question the hypothesis that the axial iron out-of-plane displacement is primarily responsible for the structure-dynamics correlations observed in myoglobin.

The effect of ligand dynamics on heme electronic transition band III in myoglobin

Band III is a near-infrared electronic transition at ~13,000 cm(-1) in heme proteins that has been studied extensively as a marker of protein conformational relaxation after photodissociation of the heme-bound ligand. To examine the influence of the heme pocket structure and ligand dynamics on band III, we have studied carbon monoxide recombination in a variety of myoglobin mutants after photolysis at 3 K using Fourier transform infrared temperature-derivative spectroscopy with monitoring in three spectral ranges, (1) band III, the mid-infrared region of (2) the heme-bound CO, and (3) the photodissociated CO. Here we present data on mutant myoglobins V68F and L29W, which both exhibit pronounced ligand movements at low temperature. From spectral and kinetic analyses in the mid-infrared, a small number of photoproduct populations can be distinguished, differing in their distal heme pocket conformations and/or CO locations. We have decomposed band III into its individual photoproduct contributions. Each photoproduct state exhibits a different "kinetic hole-burning" (KHB) effect, a coupling of the activation enthalpy for rebinding to the position of band III. The analysis reveals that the heme pocket structure and the photodissociated CO markedly affect the band III transition. A strong kinetic hole-burning effect results only when the CO ligand resides in the docking site on top of the heme group. Migration of CO away from the heme group leads to an overall blue shift of band III. Consequently, band III can be used as a sensitive tool to study ligand dynamics after photodissociation in heme proteins.

Geminate rebinding in trehalose-glass embedded myoglobins reveals residue-specific control of intramolecular trajectories

It is becoming increasingly apparent that hydrophobic cavities (also referred to as xenon cavities) within proteins have significant functional implications. The potential functional role of these cavities in modulating the internal dynamics of carbon monoxide in myoglobin (Mb) is explored in the present study by using glassy matrices derived from trehalose to limit protein dynamics and to eliminate ligand exchange between the solvent and the protein. By varying the temperature (-15 to 65 degrees C) and humidity for samples of carbonmonoxy myoglobin embedded in trehalose-glass, it is possible to observe a hierarchy of distinct geminate recombination phases that extend from nanosecond to almost seconds that can be directly associated with rebinding from specific hydrophobic cavities. The use of mutant forms of Mb reveals the role of key residues in modulating ligand access between these cavities and the distal hemepocket.