Ligand binding to heme proteins: II. Transitions in the heme pocket of myoglobin.

Phenomena occurring in the heme pocket after photolysis of carbonmonoxymyoglobin (MbCO) below about 100 K are investigated using temperature-derivative spectroscopy of the infrared absorption bands of CO. MbCO exists in three conformations (A substrates) that are distinguished by the stretch bands of the bound CO. We establish connections among the A substates and the substates of the photoproduct (B substates) using Fourier-transform infrared spectroscopy together with kinetic experiments on MbCO solution samples at different pH and on orthorhombic crystals. There is no one-to-one mapping between the A and B substates; in some cases, more than one B substate corresponds to a particular A substate. Rebinding is not simply a reversal of dissociation; transitions between B substates occur before rebinding. We measure the nonequilibrium populations of the B substates after photolysis below 25 K and determine the kinetics of B substate transitions leading to equilibrium. Transitions between B substates occur even at 4 K, whereas those between A substates have only been observed above about 160 K. The transitions between the B substates are nonexponential in time, providing evidence for a distribution of substates. The temperature dependence of the B substate transitions implies that they occur mainly by quantum-mechanical tunneling below 10 K. Taken together, the observations suggest that the transitions between the B substates within the same A substate reflect motions of the CO in the heme pocket and not conformational changes. Geminate rebinding of CO to Mb, monitored in the Soret band, depends on pH. Observation of geminate rebinding to the A substates in the infrared indicates that the pH dependence results from a population shift among the substates and not from a change of the rebinding to an individual A substate.     

Rebinding and relaxation in the myoglobin pocket

The infrared stretching bands of carboxymyoglobin (MbCO) and the rebinding of CO to Mb after photodissociation have been studied in the temperature range 10-300 K in a variety of solvents. Four stretching bands imply that MbCO can exist in four substates, A0-A3. The temperature dependences of the intensities of the four bands yield the relative binding enthalpies and and entropies. The integrated absorbances and pH dependences of the bands permit identification of the substates with the conformations observed in the X-ray data (Kuriyan et al., J. Mol. Biol. 192 (1986) 133). At low pH, A0 is hydrogen-bonded to His E7. The substates A0-A3 interconvert above about 180 K in a 75% glycerol/water solvent and above 270 K in buffered water. No major interconversion is seen at any temperature if MbCO is embedded in a solid polyvinyl alcohol matrix. The dependence of the transition on solvent characteristics is explained as a slaved glass transition. After photodissociation at low temperature the CO is in the heme pocket B. The resulting CO stretching bands which are identified as B substates are blue-shifted from those of the A substates. At 40 K, rebinding after flash photolysis has been studied in the Soret, the near-infrared, and the integrated A and B substates. All data lie on the same rebinding curve and demonstrate that rebinding is nonexponential in time from at least 100 ns to 100 ks. No evidence for discrete exponentials is found. Flash photolysis with monitoring in the infrared region shows four different pathways within the pocket B to the bound substates Ai. Rebinding in each of the four pathways B----A is nonexponential in time to at least 10 ks and the four pathways have different kinetics below 180 K. From the time and temperature dependence of the rebinding, activation enthalpy distributions g(HBA) and preexponentials ABA are extracted. No pumping from one A substate to another, or one B substate to another, is observed below the transition temperature of about 180 K. If MbCO is exposed to intense white light for 10-10(3) s before being fully photolyzed by a laser flash, the amplitude of the long-lived states increases. The effect is explained in terms of a hierarchy of substates and substate symmetry breaking. The characteristics of the CO stretching bands and of the rebinding processes in the heme pocket depend strongly on the external parameters of solvent, pH and pressure. This sensitivity suggests possible control mechanisms for protein reactions.     

Ligand binding to heme proteins. VI. Interconversion of taxonomic substates in carbonmonoxymyoglobin.

The kinetic properties of the three taxonomic A substates of sperm whale carbonmonoxy myoglobin in 75% glycerol/buffer are studied by flash photolysis with monitoring in the infrared stretch bands of bound CO at nu(A0) approximately 1967 cm-1, nu(A1) approximately 1947 cm-1, and nu(A3) approximately 1929 cm-1 between 60 and 300 K. Below 160 K the photodissociated CO rebinds from the heme pocket, no interconversion among the A substates is observed, and rebinding in each A substate is nonexponential in time and described by a different temperature-independent distribution of enthalpy barriers with a different preexponential. Measurements in the electronic bands, e.g., the Soret, contain contributions of all three A substates and can, therefore, be only approximately modeled with a single enthalpy distribution and a single preexponential. The bond formation step at the heme is fastest for the A0 substate, intermediate for the A1 substate, and slowest for A3. Rebinding between 200 and 300 K displays several processes, including geminate rebinding, rebinding after ligand escape to the solvent, and interconversion among the A substates. Different kinetics are measured in each of the A bands for times shorter than the characteristic time of fluctuations among the A substates. At longer times, fluctuational averaging yields the same kinetics in all three A substates. The interconversion rates between A1 and A3 are determined from the time when the scaled kinetic traces of the two substates merge. Fluctuations between A1 and A3 are much faster than those between A0 and either A1 or A3, so A1 and A3 appear as one kinetic species in the exchange with A0. The maximum-entropy method is used to extract the distribution of rate coefficients for the interconversion process A0 <--> A1 + A3 from the flash photolysis data. The temperature dependencies of the A substate interconversion processes are fitted with a non-Arrhenius expression similar to that used to describe relaxation processes in glasses. At 300 K the interconversion time for A0 <--> A1 + A3 is 10 microseconds, and extrapolation yields approximately 1 ns for A1 <--> A3. The pronounced kinetic differences imply different structural rearrangements. Crystallographic data support this conclusion: They show that formation of the A0 substate involves a major change of the protein structure; the distal histidine rotates about the C(alpha)-C(beta) bond, and its imidazole sidechain swings out of the heme pocket into the solvent, whereas it remains in the heme pocket in the A1 <--> A3 interconversion. The fast A1 <--> A3 exchange is inconsistent with structural models that involve differences in the protonation between A1 and A3.     

Ligand binding to heme proteins: connection between dynamics and function

Ligand binding to heme proteins is studied by using flash photolysis over wide ranges in time (100 ns-1 ks) and temperature (10-320 K). Below about 200 K in 75% glycerol/water solvent, ligand rebinding occurs from the heme pocket and is nonexponential in time. The kinetics is explained by a distribution, g(H), of the enthalpic barrier of height H between the pocket and the bound state. Above 170 K rebinding slows markedly. Previously we interpreted the slowing as a "matrix process" resulting from the ligand entering the protein matrix before rebinding. Experiments on band III, an inhomogeneously broadened charge-transfer band near 760 nm (approximately 13,000 cm-1) in the photolyzed state (Mb*) of (carbonmonoxy)myoglobin (MbCO), force us to reinterpret the data. Kinetic hole-burning measurements on band III in Mb* establish a relation between the position of a homogeneous component of band III and the barrier H. Since band III is red-shifted by 116 cm-1 in Mb* compared with Mb, the relation implies that the barrier in relaxed Mb is 12 kJ/mol higher than in Mb*. The slowing of the rebinding kinetics above 170 K hence is caused by the relaxation Mb*----Mb, as suggested by Agmon and Hopfield [(1983) J. Chem. Phys. 79, 2042-2053]. This conclusion is supported by a fit to the rebinding data between 160 and 290 K which indicates that the entire distribution g(H) shifts. Above about 200 K, equilibrium fluctuations among conformational substates open pathways for the ligands through the protein matrix and also narrow the rate distribution. The protein relaxations and fluctuations are nonexponential in time and non-Arrhenius in temperature, suggesting a collective nature for these protein motions. The relaxation Mb*----Mb is essentially independent of the solvent viscosity, implying that this motion involves internal parts of the protein. The protein fluctuations responsible for the opening of the pathways, however, depend strongly on the solvent viscosity, suggesting that a large part of the protein participates. While the detailed studies concern MbCO, similar data have been obtained for MbO2 and CO binding to the beta chains of human hemoglobin and hemoglobin Zürich. The results show that protein dynamics is essential for protein function and that the association coefficient for binding from the solvent at physiological temperatures in all these heme proteins is governed by the barrier at the heme.     

Structural dynamics of myoglobin: ligand migration among protein cavities studied by Fourier transform infrared/temperature derivative spectroscopy.

Fourier transform infrared (FTIR) spectroscopy in the CO stretch bands combined with temperature derivative spectroscopy (TDS) was used to characterize intermediate states obtained by photolysis of two sperm whale mutant myoglobins, YQR (L29(B10)Y, H64(E7)Q, T67(E10)R) and YQRF (with an additional I107(G8)F replacement). Both mutants assume two different bound-state conformations, A(0) and A(3), which can be distinguished by their different CO bands near 1965 and 1933 cm(-1). They most likely originate from different conformations of the Gln-64 side chain. Within each A substate, a number of photoproduct states have been characterized on the basis of the temperature dependence of recombination in TDS experiments. Different locations and orientations of the ligand within the protein can be distinguished by the infrared spectra of the photolyzed CO. Recombination from the primary docking site, B, near the heme dominates below 50 K. Above 60 K, ligand rebinding occurs predominantly from a secondary docking site, C', in which the CO is trapped in the Xe4 cavity on the distal side, as shown by crystallography of photolyzed YQR and L29W myoglobin CO. Another kinetic state (C") has been identified from which rebinding occurs around 130 K. Moreover, a population appearing above the solvent glass transition at approximately 180 K (D state) is assigned to rebinding from the Xe1 cavity, as suggested by the photoproduct structure of the L29W sperm whale myoglobin mutant. For both the YQR and YQRF mutants, rebinding from the B sites near the heme differs for the two A substates, supporting the view that the return of the ligand from the C', C", and D states is not governed by the recombination barrier at the heme iron but rather by migration to the active site. Comparison of YQR and YQRF shows that access to the Xe4 site (C') is severely restricted by introduction of the bulky Phe side chain at position 107.     

Dynamics of carbon monoxide binding with neuronal nitric oxide synthase.

The dynamics of CO rebinding with neuronal NO synthase (nNOS) following laser flash photolysis have been investigated from 293 to 77 K in the absence and presence of its substrate L-arginine. The distribution functions of the rate parameters P(k) and of the activation enthalpy P(H) were determined using the maximum entropy method. In a fluid solvent near room temperature, bimolecular rebinding is biphasic, as previously reported by several groups. However, measurement of the rotational correlation time shows that the apparent biphasic rebinding is not relevant to the genuine dynamics of NOS. In addition to native dimeric nNOS, another species (possibly aggregated or partially unfolded conformation) with different hydrodynamic characteristics is responsible for the faster rebinding process. In a rigid environment at low temperature, the geminate internal rebinding is not affected by the presence of the nonnative species. nNOS exhibits a bimodal distribution of CO activation enthalpy with P(H) consisting of two distinct bands with temperature-dependent amplitudes down to 77 K. The similarity of these findings with those recently reported for cytochromes P-450 suggests a common hierarchical organization of conformational substates, with a splitting of each conformational substate into a doublet. Thus, thiolate-coordinated heme proteins are in clear contrast to histidine-coordinated oxygen-transport heme proteins. The present results with nNOS provide additional support to previous arguments incriminating the thiolate ligand as responsible for the splitting of conformational substates.     

Viscosity-dependent energy barriers and equilibrium conformational fluctuations in oxygen recombination with hemerythrin

The recombination kinetics of photo-dissociated oxyhemerythrin (Sipunculus nudus) have been investigated between 298 K and 90 K. Fast geminate recombinations compete with oxygen escape into the solvent, from which a subsequent slower bimolecular rebinding takes place. In phosphate buffer (pH 7.7) at 278 K, the fast and slow processes are exponential and have comparable amplitudes. Whereas the oxygen escape rate rapidly decreases upon increasing the viscosity, the inward rate from the solvent is found to be independent of viscosity, up to about 50 cP (50 mPa.s). The data suggest that a Brownian-motion-driven displacement of one or several side-chain residues is implied in oxygen escape from within the protein and also that hemerythrin undergoes a conformational change in the deoxy state. At higher viscosities and lower temperature only the geminate phase is observed and the kinetics progressively depart from an exponential. Below about 130 K, the kinetics resemble those reported in the literature for heme proteins. They are consistent with a temperature-independent non-equilibrium frozen distribution of conformational substates. However, between 190 K and 130 K, the profile of the kinetics is invariant on a log/log plot and the results simply differ by a translation along the log t axis. It is shown that this property is expected only for a temperature-dependent distribution of substates in a Boltzmann equilibrium. From room temperature, where rebinding is exponential, down to the 'freezing' temperature, the geminate recombinations display a variety of kinetic laws. It can be shown, however, that for a broad class of substate distributions, the initial slope of the kinetic plot follows an Arrhenius relationship. The activation energy is equal to that of the exponential rate constant measured at high temperature. This result establishes the conditions under which protein data obtained from low-temperature kinetics can be extrapolated to physiological temperature.     

Ligand binding to heme proteins. V. Light-induced relaxation in proximal mutants L89I and H97F of carbonmonoxymyoglobin.

We have studied the proximal mutants L89I and H97F of MbCO with FTIR and temperature-derivative spectroscopy at temperatures between 10 and 160 K. The mutations give rise only to minor alterations of the stretch spectra of the bound and photodissociated CO ligand. The most pronounced difference is a larger population in the A3 substate at approximately 1930 cm-1 in the mutants. The barrier distributions, as determined by temperature-derivative spectroscopy, are very similar to native MbCO after short illumination. Extended illumination leads to substantial increases of the rebinding barriers in native MbCO and the proximal mutants. A larger fraction of light-relaxed states is found in the proximal mutants, implying that the conformational energy landscape has been modified to more easily allow light-induced transitions. These and other spectroscopic data imply that the large changes in the binding properties are brought about by a light-induced conformational relaxation involving the structure at the heme iron. Similarities with spectral hole-burning studies and physical models are discussed.     

Reactive line-shape narrowing in low-temperature inhomogeneous geminate recombination of CO to myoglobin

The temporal shift in the near-IR absorption peak of myoglobin (Mb) following flash photolysis of MbCO at cryogenic temperatures appears to be due largely to an inhomogeneous reactive process rather than to relaxation. This conclusion, which follows from a new analysis of the experimental data, is based on the following three points: First, at very low temperatures (60 K) a transient line-narrowing effect can be detected. Second, there is a universal, temperature-independent, correlation between spectral shift and survival probability in the rebinding kinetics, and third, the same quantitative model which accounts for rebinding accounts semiquantitatively for the temporal shift in the peak. A fit to the model indicates that the inhomogeneous broadening of the near-IR peak in myoglobin is 15-20% of the total width. The same rebinding process which governs the loss of intensity of this peak is therefore most likely responsible for the shift in its center wavelength.     

The temperature dependence of steady-state kinetics: What can be learned about pig liver esterase stereospecificity?

The steady-state kinetic parameters for pig liver carboxylesterase (PLE)-catalyzed hydrolysis of the prochiral substrate dimethyl phenylmalonate (DMPM) (product enantioselectivity) and the separate enantiomers of three chiral 2-phenylpropionic acid esters (substrate enantioselectivity) were measured at seven temperatures between 288 K and 312 K. Arrhenius plots of turnover numbers against the reciprocal of experimental temperatures yielded enthalpies and entropies of activation at enzyme saturation. (+)-(S)-methyl-2-phenylpropionate, (+)-(S)-4-nitrophenyl 2-phenylpropionate, and both enantiomers of phenyl 2-phenylpropionate showed very similar activation enthalpies and entropies (approximately 50 kJ mol^?1 and ?50 J mol^?1 K^?1, respectively), but differences were observed for (?)-(R)-methyl 2-phenylpropionate and (?)-(R)-4-nitrophenyl 2-phenylpropionate. Whereas the entropies of activation of all 2-phenylpropionates were negative, positive entropies of activation were measured in the formation of monomethyl phenylmalonate enantiomers from DMPM. Enthalpy–entropy compensation analysis of the data indicates a common mechanism of PLE substrate and product enantiospecificity in the reactions studied here. © 1994 Wiley-Liss, Inc.     

Dynamics of carbon monoxide recombination to fully reduced cytochrome c oxidase in plant mitochondria after low-temperature flash photolysis.

Rebinding of CO to reduced cytochrome c oxidase in plant mitochondria has been monitored optically at 590-630 nm after flash photolysis at low temperature from 160 to 200 K. (1) Under 100%-CO saturation, CO rebinding exhibits a four-step mechanism. The thermodynamic parameters of the first phase have been determined; its activation energy, Ea1, is 38.9 kJ.mol-1 and its enthalpy, delta H+/-1, and entropy, delta S+/-1, of activation are respectively 37.5 kJ.mol-1 and -75.8J.mol-1.K-1. (2) When the CO concentration is decreased to 0.2%, rebinding still occurs according to a four-step mechanism. The rate constant of the first phase is CO-concentration-independent. Under non-saturating conditions there is only one CO molecule per occupied site. The rebinding mechanism does not require additional CO molecules to be present in the haem pocket. (3) Dual-wavelength scanning experiments failed to detect optical forms correlated with the resolved phases. (4) Results are discussed with respect to previous work related to CO rebinding to mammalian cytochrome c oxidase and myoglobin.     

Conformational relaxation in hemoproteins: the cytochrome P-450cam case

Photodissociation of (CO)P-450(cam)(substrate) complexes was found to trigger a conformational relaxation process that interferes with ligand rebinding at temperatures as low as 140 K even though the protein conformational substates (CS(1)) remain frozen. To analyze the rebinding and relaxation kinetics, we developed a model that takes the distribution of relaxation rates explicitly into account and in which rebinding and relaxation rates are connected by a linear free energy relation. In all complexes heme relaxation occurs first and is probably faster than 100 ns even at 77 K. This is the only process found in substrate-free P-450(cam). Above 140 K and in the presence of a substrate, this initial, fast rebinding state (P) progressively relaxes to another state (P degrees ) in which rebinding is slower. The relaxation rate is independent of solvent rigidity and is governed by the protein's internal dynamics. Rebinding enthalpies in P and P degrees as well as the enthalpy shift brought about by relaxation correlate with the substrate propensity to block access to the iron site. In P degrees the barrier is higher because the substrate is closer to the heme normal and exerts more steric repulsion for CO binding. The relaxation process implies the return of substrate and heme to their ligand-free positions in which access to the heme is reduced.     

Filtering artifacts from lifetime distributions when maximizing entropy using a bootstrapped model

The maximum entropy method (MEM) has been used in many studies to reliably recover effective lifetimes from kinetics, whether measured experimentally or simulated computationally. Here, recent claims made by Mulligan et al. regarding MEM analyses of kinetics (Anal. Biochem. 421 (2012) 181-190) are shown to be unfounded. Their assertion that their software allows "analysis of datasets too noisy to process by existing iterative search algorithms" is refuted with a MEM analysis of their triexponential test case with increased noise. In addition, it is shown that lifetime distributions recovered from noisy kinetics data with the MEM can be improved by using a simple filter when bootstrapping the prior model. When deriving the bootstrapped model from the lifetime distribution obtained using a uniform model, only the slower processes are represented as Gaussians in the bootstrapped model. Using this new approach, results are clearly superior to those of Mulligan et al. despite the presence of increased noise. In a second example, ambiguity in the interpretation of Poisson kinetics in the presence of scattered excitation light is resolved by filtering the prior model.     

Conformational substates and dynamic properties of carbonmonoxy hemoglobin

Heme pocket dynamics of human carbonmonoxy hemoglobin (HbCO) is studied by Fourier transform infrared spectroscopy. The CO stretching band at various temperatures in the interval 300-10 K is analyzed in terms of three taxonomic A substates; however, in HbCO the band attributed to the A(1) taxonomic substate accounts for approximately 90% of the total intensity in the pH range 8.8-4.5. Two different regimes as a function of temperature are observed: below 160 K, the peak frequency and the bandwidth of the A(1) band have constant values whereas, above this temperature, a linear temperature dependence is observed, suggesting the occurrence of transitions between statistical substates within the A(1) taxonomic substate in this protein. The relationship between the heme pocket dynamics (as monitored by the thermal behavior of the CO stretching band), the overall dynamic properties of the protein matrix (as monitored by the thermal behavior of Amide II and Amide I' bands) and the glass transition of the solvent (as monitored by the thermal behavior of the bending band of water) is also investigated. From this analysis, we derive the picture of a very soft heme pocket of hemoglobin characterized by rather large anharmonic terms and strongly coupled to the dynamic properties of the solvent.     

Electron paramagnetic resonance- (EPR-) resolved kinetics of cryogenic nitric oxide recombination to cytochrome c oxidase and myoglobin.

By the electron paramagnetic resonance (EPR) technique, recovery kinetics for nitric oxide (NO) to heme following cryogenic photolysis were studied for the nitrosylferrocytochrome a3 center in cytochrome c oxidase and for myoglobin. The recovery was nonexponential, as has been observed in previous cryogenic CO and O2 rebinding to heme systems. NO rebinding to heme a3 started near a temperature of 50 K and was related to a distribution of thermal activation energies. At the peak of the distribution the activation energy was 3.1 kcal/mol, and the preexponential in the recovery rate was 10(9.9) s-1. For recovery of NO back to the a3 heme, the activation energy was threefold less than that for CO where CO binds to nearby Cua3 following photolysis from heme a3, but was larger than the activation energy for CO, O2, and probably NO rebinding to myoglobin. NO ligand rebinding to myoglobin occurred at a temperature as low as 15 K and in a temperature regime where tunneling could occur. However, the rate of NO rebinding to myoglobin did increase with temperature in the 15-25 K range.     

Competition with Xenon Elicits Ligand Migration and Escape Pathways in Myoglobin

Evidence for ligand migration toward the xenon-binding cavities in myoglobin comes from a number of laser photolysis studies of MbO₂ including mutants and from cryo- and time-resolved crystallography of MbCO. To explore ligand migration in greater detail, we investigated the rebinding kinetics of both MbO₂ and MbCO under a xenon partial pressure ranging from 1 to 16 atm over the temperature range (293–77 K). Below 180 K xenon affects to a significant, but minor, extent the thermodynamic parameters for rebinding from the primary docking site in each Mb taxonomic substate. Above 200 K the ligand migrates to the proximal Xe1 site but when the latter is occupied by xenon a new kinetic process appears. It is attributed to rebinding from transient docking sites located on the path between the primary and the secondary docking site of both ligands. Ligand escape exhibits a more complicated pattern than expected. At room temperature O₂ and CO escape appears to take place exclusively from the primary site. In contrast, at T ≈ 250 K, roughly 50% of the CO molecules that have escaped from the protein originate from the Xe1 secondary site.     

Dynamics of ligand binding to myoglobin.

Myoglobin rebinding of carbon monoxide and dioxygen after photodissociation has been observed in the temperature range between 40 and 350 K. A system was constructed that records the change in optical absorption at 436 nm smoothly and without break between 2 musec and 1 ksec. Four different rebinding processes have been found. Between 40 and 160 K, a single process is observed. It is not exponential in time, but approximately given by N(t) = (1 + t/to)-n, where to and n are temperature-dependent, ligand-concentration independent, parameters. At about 170 K, a second and at 200 K, a third concentration-independent process emerge. At 210 K, a concentration-dependent process sets in. If myoglobin is embedded in a solid, only the first three can be seen, and they are all nonexponential. In a liquid glycerol-water solvent, rebinding is exponential. To interpret the data, a model is proposed in which the ligand molecule, on its way from the solvent to the binding site at the ferrous heme iron, encounters four barriers in succession. The barriers are tentatively identified with known features of myoglobin. By computer-solving the differential equation for the motion of a ligand molecule over four barriers, the rates for all important steps are obtained. The temperature dependences of the rates yield enthalpy, entropy, and free-energy changes at all barriers. The free-energy barriers at 310 K indicate how myoglobin achieves specificity and order. For carbon monoxide, the heights of these barriers increase toward the inside; carbon monoxide consequently is partially rejected at each of the four barriers. Dioxygen, in contrast, sees barriers of about equal height and moves smoothly toward the binding site. The entropy increases over the first two barriers, indicating a rupturing of bonds or displacement of residues, and then smoothly decreases, reaching a minimum at the binding site. The magnitude of the decrease over the innermost barrier implies participation of heme and/or protein. The nonexponential rebinding observed at low temperatures and in solid samples implies that the innermost barrier has a spectrum of activation energies. The shape of the spectrum has been determined; its existence can be explained by assuming the presence of many conformational states for myoglobin. In a liquid at temperatures above about 230 K, relaxation among conformational states occurs and rebinding becomes exponential.     

Kinetic evidence for three photolyzable taxonomic conformational substates in oxymyoglobin

The kinetics of oxygen geminate binding with the taxonomic substates of MbO2 are reported. The maximum entropy method was used to analyze the rebinding kinetics of MbCO and MbO2 monitored in the Soret. The resulting rate distributions were found to consist of a small number of overlapping bands. A global parametric fit of a series of rate distributions recorded at several temperatures was performed using a Gaussian basis set to resolve the individual enthalpy distributions P(H). This approach was first validated by showing that the well-documented taxonomic substates of MbCO could be recovered. The method was then applied to MbO2. Three taxonomic substates were identified at pH 4.8, whereas only two of them contribute to oxygen geminate rebinding at pH 7.0. These findings show that, similarly to MbCO, MbO2 also exists as three photolyzable and kinetically different taxonomic substates and suggest reconsidering the issue of the photolysis quantum yield of MbO2.     

A hierarchy of functionally important relaxations within myoglobin based on solvent effects, mutations and kinetic model

Geminate CO rebinding in myoglobin is studied for two viscous solvents, trehalose and sol-gel (bathed in 100% glycerol) at several temperatures. Mutations in key distal hemepocket residues are used to eliminate or enhance specific relaxation modes. The time-resolved data are analyzed with a modified Agmon-Hopfield model which is capable of providing excellent fits in cases where a single relaxation mode is dominant. Using this approach, we determine the relaxation rate constants of specific functionally important modes, obtaining also their Arrhenius activation energies. We find a hierarchy of distal pocket modes controlling the rebinding kinetics. The "heme access mode" (HAM) is responsible for the major slow-down in rebinding. It is a solvent-coupled cooperative mode which restricts ligand return from the xenon cavities. Bulky side-chains, like those His64 and Trp29 (in the L29W mutant), operate like overdamped pendulums which move over and block the binding site. They may be either unslaved (His64) or moderately slaved (Trp29) to the solvent. Small side-chain relaxations, most notably of leucines, are revealed in some mutants (V68L, V68A). They are conjectured to facilitate inter-cavity ligand motion. When all relaxations are arrested (H64L in trehalose), we observe pure inhomogeneous kinetics with no temperature dependence, suggesting that proximal relaxation is not a factor on the investigated timescale.     

Investigation of laser-induced long-lived states of photolyzed MbCO

We present evidence from resonance Raman and absorption measurements that the extended exposure of MbCO to CW laser light at low temperatures alters the CO rebinding kinetics and leads to a significantly increased population of very long lived states of photolyzed MbCO. This optical "pumping" process is observed for samples frozen in both aqueous buffer and glycerol/buffer and exhibits power law behavior with a very weak temperature dependence. A comparison of the nonexponential rebinding kinetics of CO molecules from the pumped states with the rebinding observed in flash photolysis experiments suggests that the pumped states are distinct geminate states, not observed in flash photolysis experiments. Thus, a four-state model, with two geminate states, is implicated for MbCO. Pumped states may represent "separated geminate pair" states with the CO molecule still in the heme pocket or possibly trapped within a cavity on its way through the protein matrix, consistent with molecular dynamics simulations. The possibility of significant deoxyheme relaxation from a less domed to a more domed configuration, as a result of the multiple photolysis events associated with the pumping process, is also explored. However, the small changes observed in the Soret band line shape and position subsequent to pumping at T less than 180 K tend to rule out this explanation for the pumping process. Since the yield for creating a pumped state is small (e.g., less than 10(-7) for T greater than 100 K), pumping can be observed only after extended illumination and is absent in flash photolysis measurements, even after multiple flashes. At higher temperatures (T greater than 180 K), the escape of the CO molecule to the solvent is observed. Our data are consistent with a "phase transition" of the protein that is coupled to the surrounding matrix. The protein fluctuations are quenched below approximately 185 K for a solvent composed of 70% glycerol and below approximately 260 K for aqueous buffer. We also present the first large amplitude measurements of CO rebinding from the protein exterior, observed below 200 K after freezing the sample under laser illumination.