Inhomogeneous broadening in spectral bands of carbonmonoxymyoglobin. The connection between spectral and functional heterogeneity.

The rebinding kinetics of CO to myoglobin after flash photolysis is nonexponential in time below approximately 180 K; the kinetics is governed by a distribution of enthalpic barriers. This distribution results from inhomogeneities in the protein conformation, referred to as conformational substates. Hole-burning experiments on the Soret and IR CO-stretch bands test the assumption that an inhomogeneous distribution of conformational substates results in inhomogeneously broadened spectra. CO was slowly photolyzed at different wavelengths in the Soret band at 10 K. Both the Soret band and the CO-stretch band A1, centered at 1,945 cm-1, shift during photolysis, demonstrating that different wavelengths excite different parts of the distributed population. We have also done kinetic hole-burning experiments by measuring peak shifts in the Soret and A1 bands as the CO molecules rebind. The shifts indicate that the spectral and enthalpic distributions are correlated. In the A1 band, the spectral and enthalpic distributions are highly correlated while in the Soret the correlation is weak. From the peak shifts in the spectral and kinetic hole-burning experiments the inhomogeneous broadening is estimated to be approximately 15% of the total width in the Soret band and approximately 60% in A1. We have previously measured the tilt angle alpha between the bound CO and the heme normal (Ormos, P., D. Braunstein, H. Frauenfelder, M. K. Hong, S.-L. Lin, T. B. Sauke, and R. D. Young. 1988. Proc. Natl. Acad. Sci. USA. 85:8492-8496) and observed a wave number dependence of the tilt angles within the CO-stretch A bands. Thus the spectral and enthalpic distributions of the A bands are coupled to a heterogeneity of the structure.     

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.     

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.     

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.     

Infrared absorption study of the heme pocket dynamics of carbonmonoxyheme proteins

The temperature dependencies of the infrared absorption CO bands of carboxy complexes of horseradish peroxidase (HRP(CO)) in glycerol/water mixture at pH 6.0 and 9.3 are interpreted using the theory of optical absorption bandshape. The bands' anharmonic behavior is explained assuming that there is a higher-energy set of conformational substates (CSS(h)), which are populated upon heating and correspond to the protein substates with disordered water molecules in the heme pocket. Analysis of the second moments of the CO bands of the carboxy complexes of myoglobin (Mb(CO)) and hemoglobin (Hb(CO)), and of HRP(CO) with benzohydroxamic acid (HRP(CO)+BHA), shows that the low energy CSS(h) exists also in the open conformation of Mb(CO), where the heme pocket is spacious enough to accommodate a water molecule. In the HRP(CO)+BHA and closed conformations of Mb(CO) and Hb(CO), the heme pocket is packed with BHA and different amino acids, the CSS(h) has much higher energy and is hardly populated even at the highest temperatures. Therefore only motions of these amino acids contribute to the band broadening. These motions are linked to the protein surface and frozen in the glassy matrix, whereas in the liquid solvent they are harmonic. Thus the second moment of the CO band is temperature-independent in glass and is proportional to the temperature in liquid. The temperature dependence of the second moment of the CO peak of HRP(CO) in the trehalose glass exhibits linear coupling to an oscillator. This oscillator can be a moving water molecule locked in the heme pocket in the whole interval of temperatures or a trehalose molecule located in the heme pocket.     

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.     

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.     

Quaternary structure dependence of kinetic hole burning and conformational substates interconversion in hemoglobin

Using a sol-gel encapsulation technique, we have prepared samples of CO saturated human adult hemoglobin locked in the R or T quaternary conformation. We report time-resolved spectra of these samples in the Soret region following flash photolysis, in the time interval ranging from 250 ns to 200 ms and in the temperature interval of 100-170 K. A suitable analysis of the measured difference spectra enables us to obtain the spectral contribution of deoxyHb and HbCO molecules as a function of time and/or of the fraction N(t) of deoxyHb molecules. In our experimental time window geminate CO rebinding to hemoglobin in the T quaternary conformation is about 2 orders of magnitude slower than to hemoglobin in the R conformation: this suggests that the barrier distribution for the CO rebinding, g(H), depends strongly on the protein quaternary structure. In our temperature interval, spectral shifts due to kinetic hole burning (KHB) are present: for HbCO the KHB effect is large in the R conformation and small in the T conformation. For deoxyHb the opposite is true. We attribute the observed behavior to the effect of interconversion between the relevant substates. This effect is stronger for HbCO molecules in the T conformation and for deoxyHb molecules in the R conformation; it confirms the quaternary structure dependence of the hemoglobin energy landscape and suggests enhanced dynamics of ligation intermediate species such as T-state HbCO or R-state deoxyHb.     

Isocyanide binding kinetics to monomeric hemoproteins. A study on the ligand partition between solvent and heme pocket.

The kinetics of methyl-, ethyl-, iso-propyl-, and ter-butyl-isocyanide binding to Aplysia limacina myoglobin (distal His----Lys) and the isolated beta chains from hemoglobin Zurich (distal His----Arg) have been investigated by flash photolysis at various temperatures above 0 degrees C. Sperm whale (Physter catodon) myoglobin and the isolated beta chains from normal adult hemoglobin have been used as references. In most reaction systems investigated the apparent extent of photolysis increases with temperature. For sperm whale myoglobin and the normal beta chains the increase is of the same magnitude and not correlated to the type of ligand used. On the contrary, for the two proteins lacking the distal histidine, the phenomenon is dependent on the size of the alkyl side chain of the ligand. The results, analyzed on the basis of the multibarrier model (Austin, R.H., K.W. Beeson, L. Eisenstein, H. Frauenfelder, and I.C. Gunsalus, 1975, Biochemistry, 16:5355-5373), suggest that the partition of the ligand molecules between the solvent and the heme pocket, occurring during the photolysis process, is primarily determined by interactions between the ligand and residues in the heme cavity rather than by diffusion through the protein matrix.     

NMR studies of the heme pocket conformations of monomeric hemoglobins from Glycera dibranchiata. Implications for ligand binding

Two-dimensional 1H-NMR methods have been used to assign side-chain resonances for the tryptophan residues and for several amino acids located in the heme pockets of the carbon monoxide complexes of the major monomeric hemoglobins from Glycera dibranchiata. The NMR spectra reveal a high degree of conservation of the heme pocket structure in the different hemoglobins. However some conformational differences are evident and residues at positions B10 and G8 on the distal side of the heme pocket are not conserved. From the present NMR studies it appears that the monomeric G. dibranchiata hemoglobin examined by X-ray crystallography [Padlan, E. A. & Love, W. (1974) J. Biol. Chem. 249, 4067-4078] corresponds to HbC. Except that the orientation of the heme in solution is the reverse of that reported in the crystal structure, there is a close correspondence between the heme pocket structure in the crystal and in solution. The proximal histidine coordination geometry is almost identical in the CO complexes of the three monomeric hemoglobins studied. Distal residues are strongly implicated in determining the observed kinetic differences in ligand binding reactions. In particular, steric crowding of the ligand binding site in hemoglobin A is probably a major factor in the slower kinetics of this component.     

Contribution of protein conformation to stereochemistry and reactivity of the active center of heme proteins and enzymes. The existence of horseradish peroxidase conformations and their possible role in the catalysis mechanism

In the spectral region 350-800 nm at 4.2 K we measured magnetic circular dichroism (MCD) spectra of the pentacoordinated complex of protcheme with 2-methylimidazole, deoxyleghemoglobin, neutral and alkaline forms of reduced horseradish peroxidase in the equilibrium states, as well as in non-equilibrium states produced by low-temperature photolysis of their carbon monoxide derivatives. Earlier the corresponding results have been obtained for myoglobin, hemoglobin and cytochromes P-450 and P-420. The energies of Fe-N (proximal His) and Fe-N(pyrroles) bonds and their changes upon ligand binding in heme proteins and enzymes were compared with those in the model heme complex thus providing conformational contribution into stereochemistry of the active site. The examples of weak and strong conformational "pressure" on stereochemistry were analysed and observed. If conformational energy contribution into stereochemistry prevails the electronic one the heme stereochemistry remains unchanged on ligand binding as it was observed for leghemoglobin and alkaline horseradish peroxidase. The change of bond energies in myoglobin and hemoglobin on ligand binding are comparable with those in protein free pentacoordinated protoheme, giving an example of weak conformational contribution to heme stereochemistry. The role of protein conformation energy in the modulation of ligand binding properties of heme in leghemoglobin relative to those in myoglobins is discussed. The most striking result were obtained in the study of reduced horseradish peroxidase in the pH region of 6.0-10.2. It was found that such different perturbations as ligand binding and heme-linked ionization of the distal amino acid residue induce identical changes in heme stereochemistry. Neither heme-linked ionization in the carbon monoxide complex nor the geometry of Fe-Co bond affect the heme local structure of photoproducts. These and other findings suggest a very low conformation mobility of horseradish peroxidase whose protein constraints appear to allow only two preferable geometries of specific amino acid residues that form the heme pocket. The role of the two tertiary structure constraints on the heme in the mechanism of horseradish peroxidase function is discussed. It is supposed that one conformation produces a heme environment suitable for two-electron oxidation of the native enzyme to compound I by hydrogen peroxide while another conformation changes the heme stereochemistry in the direction favourable for back reduction of compound I by the substrate to the resting enzyme through two one-electron steps. The switch from one tertiary structure to another is expected to be induced by substrate bind     

Structural studies of yeast iso-1 cytochrome c mutants by resonance Raman spectroscopy.

The Ser82 and Phe82 variants of yeast iso-1 cytochrome c were studied by resonance Raman spectroscopy. In both oxidation states, distinct spectral changes were observed for some of those bands in the low-frequency region, which sensitively respond to conformational perturbations of the protein environment of the heme. These bands can be assigned to modes which include strong contributions of vibrations largely localized in the propionate-carrying pyrrole rings A and D. This indicates structural differences in the deeper part of the heme crevice, remote from the mutation site. This conclusion is in line with previous results from X-ray crystallography and NMR spectroscopy. No differences in the resonance-Raman spectra were observed which can be directly correlated with conformational changes of the heme pocket in the vicinity of the mutation site. Temperature-dependent resonance Raman experiments of the oxidized mutants revealed spectral changes which are closely related to those observed for cytochrome c upon adsorption to charged silver surfaces by surface-enhanced resonance Raman spectroscopy. These spectral changes can be attributed to an opening of the heme crevice accompanied by a weakening of the iron-methionine ligand bond. The temperature-dependent conformational transition occurs at approximately 30 degrees C for the Ser82 variant and at about 45 degrees C for the Phe82 variant, implying that the Phe----Ser substitution significantly lowers the thermal stability of the heme pocket. The reduced forms of both mutants are stable up to 65 degrees C.     

NH2-terminal spin labeling of human hemoglobin--spin states and temperature dependence

In order to explore fully how ligand- and temperature-induced alterations in the spin states of heme iron are related to protein readjustments, the spin label 4-isothiocyanate (I) was covalently attached at beta-93 cysteines and at NH2-terminal valines of various heme-iron ligand forms of human hemoglobin. It was found that the mobility of NH2-terminally bound spin labels depends on the magnetic moment of the heme iron. There is a an approximately linear relationship between the magnetic moment of the heme iron and the mobility of NH2-terminally bound spin labels. In accordance with our previous results, the temperature dependence of ESR spectra of spin-labeled hemoglobin suggests the temperature-induced protein conformational change in those heme-iron ligand forms that are characterized by the equilibrium of the spin states of the heme iron. The conformational change was sensed at both spin-label-binding sites: at beta-93 cysteines and at NH2-terminal valines.     

Ionic channels with conformational substates.

Recent studies of protein dynamics suggest that ionic channels can assume many conformational substates. Long-lived substates have been directly observed in single-channel current records. In many cases, however, the lifetimes of conformational states will be far below the theoretical limit of time resolution of single-channel experiments. The existence of such hidden substates may strongly influence the observable (time-averaged) properties of a channel, such as the concentration dependence of conductance. A channel exhibiting fast, voltage-dependent transitions between different conductance states may behave as an intrinsic rectifier. In the presence of more than one permeable ion species, coupling between ionic fluxes may occur, even when the channel has only a single ion-binding site. In special situations the rate of ion translocation becomes limited by the rate of conformational transitions, meaning that the channel approaches the kinetic behavior of a carrier. As a result of the strong coulombic interaction between an ion in a binding site and polar groups of the protein, rate constants of conformational transitions may depend on the occupancy of the binding site. Under this condition a nonequilibrium distribution of conformational states is created when ions are driven through the channel by an external force. This may lead to an apparent violation of microscopic reversibility, i.e., to a situation in which the frequency of transitions from state A to state B is no longer equal to the transition frequency from state B to state A.     

Sampling field heterogeneity at the heme of c-type cytochromes by spectral hole burning spectroscopy and electrostatic calculations

We report on a comparative investigation of the heme pocket fields of two Zn-substituted c-type cytochromes-namely yeast and horse heart cytochromes c-using a combination of hole burning Stark spectroscopy and electrostatic calculations. The spectral hole burning experiments are consistent with different pocket fields experienced at the hemes of the respective cytochromes. In the case of horse heart Zn-cytochrome c, two distinguishable electronic origins with different electrostatic properties are observed. The yeast species, on the other hand, displays a single electronic origin. Electrostatic calculations and graphics modeling using the linearized finite-difference Poisson-Boltzmann equation performed at selected time intervals on nanosecond-molecular dynamics trajectories show that the hemes of the respective cytochromes sample different potentials as they explore conformational space. The electrostatic potentials generated by the protein matrix at the heme show different patterns in both cytochromes, and we suggest that the cytochromes differ by the number of "electrostatic substates" that they can sample, thus accounting for the different spectral populations observed in the two cytochromes.     

Structural heterogeneity in protein crystals

Extensive conformational heterogeneity is reported in highly refined crystallographic models for the proteins crambin, erabutoxin, myohemerythrin, and lamprey hemoglobin. From 6% to 13% of the amino acid side chains of these four proteins are seen in multiple, discrete conformations. Most common are flexible side chains on the molecular surface, but structural heterogeneity occasionally extends to buried side chains or to the polypeptide backbone. A few instances of sequence heterogeneity are also very clear. Numerous solvent sites are multiplets, and at high resolution, multiple, mutually exclusive solvent networks are observed. The proteins have been studied with X-ray diffraction data extending to spacings of from 0.945 to 2.0 A. The extensive heterogeneity observed here provides detailed, accurate structures for conformational substates of these molecules and sets a lower bound on the number of substates accessible to each protein molecule in solution. Electron density is missing or very weak for only a few side chains in these protein crystals, revealing a strong preference for discrete over continuous conformational perturbations. The results at very high resolution further suggest that even rather small conformational fluctuations produce discrete substates and that unresolved conformers are accommodated in increased atomic thermal parameters.     

Mass-weighted molecular dynamics simulation of the protein-ligand complex of rhizopuspepsin and inhibitor.

The mass-weighted molecular dynamics simulation method was developed previously for sampling the multidimensional conformational space of linear and cyclic polypeptides and studying their conformational flexibility. Herein results from molecular dynamics simulations of the protein-ligand complex of the aspartyl protease rhizopuspepsin and a polypeptide inhibitor are reported. The dihedral conformational space sampling for the linear peptide inhibitor in situ was found to be increased in the mass-weighted simulation as in other molecular systems previously studied. More significantly, the physical space of the enzyme binding pocket was also sampled efficiently in the simulations and multiple binding sites were identified for the inhibitor. These results suggest that it may be possible now to study, by computer simulations, the putative initial enzyme-inhibitor complex suggested experimentally from the time-dependent kinetics of enzyme inhibition by slow-binding inhibitors (Morrison, J. F., and C. T. Walsh. 1988. Adv. Enzymol. 61:201), and/or conformational substates in protein-ligand complexes suggested in the study of reassociation dynamics of myoglobin and carbon monoxide following photolysis (Austin, R. H., K. W. Beeson, L. Eisenstein, H. Frauenfelder, and I. C. Gunsalus. 1975. Biochemistry. 14:5355). Moreover, the intermediate binding steps and the molecular flexibility of the inhibitor shown in the MWMD simulation may have crucial roles in the ligand binding process.     

Nanosecond time-resolved absorption studies of human oxyhemoglobin photolysis intermediates.

The time-resolved spectra of photoproducts from ligand photodissociation of oxyhemoglobin are measured in the Soret spectral region for times from 10 ns to 320 microseconds after laser photolysis. Four processes are detected at a heme concentration of 80 microM: a 38-ns geminate recombination, a 137-ns tertiary relaxation, and two bimolecular processes for rebinding of molecular oxygen. The pseudo-first-order rate constants for rebinding to the alpha and beta subunits of hemoglobin are 3.2 x 10(4) s-1 (31 microseconds lifetime) and 9.4 x 10(4) s-1 (11 microseconds lifetime), respectively. The significance of kinetic measurements made at different heme concentrations is discussed in terms of the equilibrium compositions of hemoglobin tetramer and dimer mixtures. The rebinding rate constants for alpha and beta chains are observed to be about two times slower in the dimer than in the tetramer, a finding that appears to support the observation of quaternary enhancement in equilibrium ligand binding by hemoglobin tetramers.     

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.