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: 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.     

Myoglobin-CO conformational substate dynamics: 2D vibrational echoes and MD simulations

Two-dimensional (2D) infrared vibrational echoes were performed on horse heart carbonmonoxymyoglobin (MbCO) in water over a range of temperatures. The A(1) and A(3) conformational substates of MbCO are found to have different dephasing rates with different temperature dependences. A frequency-frequency correlation function derived from molecular dynamics simulations on MbCO at 298 K is used to calculate the vibrational echo decay. The calculated decay shows substantial agreement with the experimentally measured decays. The 2D vibrational echo probes protein dynamics and provides an observable that can be used to test structural assignments for the MbCO conformational substates.     

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.     

Molecular dynamics simulation of sucrose- and trehalose-coated carboxy-myoglobin

We performed a room temperature molecular dynamics (MD) simulation on a system containing 1 carboxy-myoglobin (MbCO) molecule in a sucrose-water matrix of identical composition (89% [sucrose/(sucrose + water)] w/w) as for a previous trehalose-water-MbCO simulation (Cottone et al., Biophys J 2001;80:931-938). Results show that, as for trehalose, the amplitude of protein atomic mean-square fluctuations, on the nanosecond timescale, is reduced with respect to aqueous solutions also in sucrose. A detailed comparison as a function of residue number evidences mobility differences along the protein backbone, which can be related to a different efficacy in bioprotection. Different heme pocket structures are observed in the 2 systems. The joint distribution of the magnitude of the electric field at the CO oxygen atom and of the angle between the field and the CO unit vector shows a secondary maximum in sucrose, absent in trehalose. This can explain the CO stretching band profile (A substates distribution) differences evidenced by infrared spectroscopy in sucrose- and trehalose-coated MbCO (Giuffrida et al., J Phys Chem B 2004;108:15415-15421), and in particular the appearance of a further substate in sucrose. Analysis of hydrogen bonds at the protein-solvent interface shows that the fraction of water molecules shared between the protein and the sugar is lower in sucrose than in trehalose, in spite of a larger number of water molecules bound to the protein in the former system, thus indicating a lower protein-matrix coupling, as recently observed by Fourier transform infrared (FTIR) experiments (Giuffrida et al., J Phys Chem B 2004;108:15415-15421).     

The distal residue-CO interaction in carbonmonoxy myoglobins: a molecular dynamics study of two distal histidine tautomers.

Four independent 90 ps molecular dynamics simulations of sperm-whale wild-type carbonmonoxy myoglobin (MbCO) have been calculated using a new AMBER force field for the haem prosthetic group. Two trajectories have the distal 64N delta nitrogen protonated, and two have the 64N epsilon nitrogen protonated; all water molecules within 16 A of the carbonyl O are included. In three trajectories, the distal residue remains part of the haem pocket, with the protonated distal nitrogen pointing into the active site. This is in contrast with the neutron diffraction crystal structure, but is consistent with the solution phase CO stretching frequencies (upsilon CO) of MbCO and various of its mutants. There are significant differences in the "closed" pocket structures found for each tautomer: the 64N epsilon H trajectories both show stable distal-CO interactions, whereas the 64N delta H tautomer) has a weaker interaction resulting in a more mobile distal side chain. One trajectory (a 64N delta H tautomer) has the distal histidine moving out into the "solvent", leaving the pocket in an "open" structure, with a large unhindered entrance to the active site. These trajectories suggest that the three upsilon CO frequencies observed for wild-type MbCO in solution, rather than representing significantly different Fe-C-O geometries as such, arise from three different haem pocket structures, each with different electric fields at the ligand. Each pocket structure corresponds to a different distal histidine conformer: the A3 band to the 64N epsilon H tautomer, the A1,2 band to the 64N delta H tautomer, and the A0 band to the absence of any significant interaction with the distal side chain.     

Dehydration and crystallization of trehalose and sucrose glasses containing carbonmonoxy-myoglobin

We report a study wherein we contemporarily measured 1) the dehydration process of trehalose or sucrose glasses embedding carbonmonoxy-myoglobin (MbCO) and 2) the evolution of the A substates in saccharide-coated MbCO. Our results indicate that microcrystallization processes, sizeably different in the two saccharides, take place during dehydration; moreover, the microcrystalline structure is maintained unless the dry samples are equilibrated with a humidity >/=75% (>/=60%) at 25 degrees C for the trehalose (sucrose) sample. The evolution of the parameters that characterize the A substates of MbCO indicates that 1) the effects of water withdrawal are analogous in samples dried in the presence or in the absence of sugars, although much larger effects are observed in the samples without sugar; 2) the distribution of A substates is determined by the overall matrix structure and not only by the sample water content; and 3) the population of A0 substate (i. e., the substate currently put in relation with MbCO molecules having the distal histidine out of the heme pocket) is largely enhanced during the dehydration process. However, after rehumidification its population is largely decreased with respect to the values obtained, at similar water content, during the first dehydration run.     

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.     

Molecular dynamics simulation of carboxy-myoglobin embedded in a trehalose-water matrix

We report on a molecular dynamics (MD) simulation of carboxy-myoglobin (MbCO) embedded in a water-trehalose system. The mean square fluctuations of protein atoms, calculated at different temperatures in the 100-300 K range, are compared with those from a previous MD simulation on an H2O-solvated MbCO and with experimental data from M?ssbauer spectroscopy and incoherent elastic neutron scattering on trehalose-coated MbCO. The results show that, for almost all the atomic classes, the amplitude of the nonharmonic motions stemming from the interconversion among the protein's conformational substates is reduced with respect to the H2O-solvated system, and their onset is shifted toward higher temperature. Moreover, our simulation shows that, at 300 K, the heme performs confined diffusive motions as a whole, leaving the underlying harmonic vibrations unaltered.     

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.     

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.     

Structural dynamics of myoglobin: spectroscopic and structural characterization of ligand docking sites in myoglobin mutant L29W

We have studied CO binding to the heme and CO migration among protein internal cavities after photodissociation in sperm whale carbonmonoxy myoglobin (MbCO) mutant L29W using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) and kinetic experiments at cryogenic temperatures. Photoproduct intermediates, characterized by CO at particular locations in the protein, were selectively enhanced by applying special laser illumination protocols. These studies were performed on the L29W mutant protein and a series of double mutants constructed so that bulky amino acid side chains block passageways between cavities or fill these sites. Binding of xenon was also employed as an alternative means of occluding cavities. All mutants exhibit two conformations, A(I) and A(II), with distinctly different photoproduct states and ligand binding properties. These differences arise mainly from different positions of the W29 and H64 side chains in the distal heme pocket [Ostermann, A., et al. (2000) Nature 404, 205-208]. The detailed knowledge of the interplay between protein structure, protein dynamics, and ligand migration at cryogenic temperatures allowed us to develop a dynamic model that explains the slow CO and O(2) bimolecular association observed after flash photolysis at ambient temperature.     

Ligand migration and protein fluctuations in myoglobin mutant L29W

We have determined eight X-ray structures of myoglobin mutant L29W at various experimental conditions. In addition, infrared spectroscopic experiments are presented, which are discussed in the light of the X-ray structures. Two distinct conformations of the CO-ligated protein were identified, giving rise to two stretching bands of heme-bound CO. If L29W MbCO crystals are illuminated around 180 K, a deoxy species is formed. The CO molecules migrate to the proximal side of the heme and remain trapped in the so-called Xe1 cavity upon temperature decrease to 105 K. The structure of this photoproduct is almost identical to the equilibrium high-temperature deoxy Mb structure. If the temperature is cycled to increasingly higher values, CO recombination is observed. Three intermediate structures have been determined during the rebinding process. Efficient recombination occurs only above 180 K, the characteristic temperature for the onset of protein dynamics. Rebinding is remarkably slow because bulky residues His64 and Trp29 block important migration pathways of the CO molecule.     

FTIR spectroscopic study of the dynamics of conformational substates in hydrated carbonyl-myoglobin films via temperature dependence of the CO stretching band parameters.

Two hydrated carbonyl myoglobin (MbCO) films, one containing (0.30 g water)/(g MbCO) from MbCO solution in water at pH 5.5 and the other (0.32 g water)/(gMbCO) from 0.1 M potassium phosphate buffer solution at pH 6.8, were studied by FTIR spectroscopy from 293 K to 78 K at selected temperatures on cooling and reheating. Above approximately 180 K the general trend in temperature dependence of half-bandwidths, peak maxima, and band area ratios of the A1 and A3 conformer bands is similar to those reported by Ansari et al. (1987. Biophys. J. 26:337) for MbCO in 75% glycerol/water solution, but abrupt changes in slopes at approximately 180-200 K and freezing-in of conformer populations, which could be taken as indicator for glass transition of the solvent or the protein, are absent for the hydrated MbCO films. This is interpreted in terms of an exceptionally broad distribution of relaxation times, and is in accord with conclusions from recent calorimetric annealing studies of hydrated protein powders (Sartor et al. 1994. Biophys. J. 66:249). Exchange between the three A conformers does not stop at approximately 180-200 K but occurs over the whole temperature region studied. These results are then discussed with respect to MbCO's behavior in the glass-->liquid transition region of glass-forming solvents, and it is concluded that, in analogy to the behavior of low-molecular-weight compounds with a distribution of rapidly interconverting conformers, freezing-in of MbCO's A conformer populations by the solvent should not be mistaken for a glass transition of MbCO.     

Ligand binding to heme proteins: III. FTIR studies of His-E7 and Val-E11 mutants of carbonmonoxymyoglobin.

Fouier-transform infrared (FTIR) difference spectra of several His-E7 and Val-E11 mutants of sperm whale carbonmonoxymyoglobin were obtained by photodissociation at cryogenic temperatures. The IR absorption of the CO ligand shows characteristic features for each of the mutants, both in the ligand-bound (A) state and in the photodissociated (B) state. For most of the mutants, a single A substate band is observed, which points to the crucial role of the His-E7 residue in determining the A substrate spectrum of the bound CO in the native structure. The fact that some of the mutants show more than one stretch band of the bound CO indicates that the appearance of multiple A substates is not exclusively connected to the presence of His-E7. In all but one mutant, multiple stretch bands of the CO in the photodissociated state are observed; these B substates are thought to arise from discrete positions and/or orientations of the photodissociated ligand in the heme pocket. The red shifts of the B bands with respect to the free-gas frequency indicate weak binding in the heme pocket. The observation of similar red shifts in microperoxidase (MP-8), where there is no residue on the distal side, suggests that the photodissociated ligand is still associated with the heme iron. Photoselection experiments were performed to determine the orientation of the bound ligand with respect to the heme normal by photolyzing small fractions of the sample with linearly polarized light at 540 nm.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

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.     

Influence of the heme pocket conformation on the structure and vibrations of the Fe-CO bond in myoglobin: a QM/MM density functional study.

The influence of the distal pocket conformation on the structure and vibrations of the heme-CO bond in carbonmonoxy myoglobin (MbCO) is investigated by means of hybrid QM/MM calculations based on density functional theory combined with a classical force field. It is shown that the heme-CO structure (QM treated) is quite rigid and not influenced by the distal pocket conformation (MM treated). This excludes any relation between FeCO distortions and the different CO absorptions observed in the infrared spectra of MbCO (A states). In contrast, both the CO stretch frequency and the strength of the CO...His64 interaction are very dependent on the orientation and tautomerization state of His64. Our calculations indicate that the CO...N(epsilon) type of approach does not contribute to the A states, whereas the CO...H-N(epsilon) interaction is the origin of the A(1) and A(3) states, the His64 residue being protonated at N(epsilon). The strength of the CO...His64 interaction is quantified, in comparison with the analogous O(2)...His64 interaction and with the observed changes in the CO stretch frequency. Additional aspects of the CO...His64 interaction and its biological implications are discussed.     

Kinetic, structural, and spectroscopic identification of geminate states of myoglobin: a ligand binding site on the reaction pathway

Elementary steps or geminate states in the reaction of gaseous ligands with transport proteins delineate the trajectory of the ligand and its rebinding to the heme. By use of kinetic studies of the 765-nm optical "conformation" band, three geminate states were identified for temperatures less than approximately 100 K. MbCO, which is accumulated by photolysis between 1.2 and approximately 10 K, was characterized by our previous optical and X-ray absorption studies [Chance, B., Fischetti, R., & Powers, L. (1983) Biochemistry 22, 3820-3829]. Between 10 and approximately 100 K, geminate states that are also identified that have recombination rates of approximately 10(3) s-1 and approximately 10(-5) s-1 (40 K). Thus, it is possible to maintain a steady-state nearly homogeneous population of the slowest recombining geminate state, Mb, by regulated continuous illumination (optical pumping). Both X-ray absorption and resonance Raman studies under similar conditions of optical pumping show that the heme structure around the iron in Mb is similar to that of MbCO. In both geminate states, the iron-proximal histidine distance remains unchanged (+/- 0.02 A) from that of MbCO while the iron to pyrrole nitrogen average distance has not fully relaxed to that of the deoxy state. In MbCO the CO remains close to iron but not bound, and the Fe...CO angle, which is bent in MbCO (127 +/- 4 degrees C), is decreased by approximately 15 degrees [Powers, L., Sessler, J. L., Woolery, G. L., & Chance, B. (1984) Biochemistry 23, 5519-5523]. The CO molecule in Mb, however, has moved approximately 0.7 A further from iron. Computer graphics modeling of the crystal structure of MbCO places the CO in a crevice in the heme pocket that is just large enough for the CO molecule end-on. Above approximately 100 K resonance Raman studies show that this structure relaxes to the deoxy state.     

Conformational substates in enzyme mechanism: the 120 K structure of alpha-lytic protease at 1.5 A resolution.

Insight into the dynamic properties of alpha-lytic protease (alpha LP) has been obtained through the use of low-temperature X-ray crystallography and multiple-conformation refinement. Previous studies of alpha LP have shown that the residues around the active site are able to move significantly to accommodate substrates of different sizes. Here we show a link between the ability to accommodate ligands and the dynamics of the binding pocket. Although the structure of alpha LP at 120 K has B-factors with a uniformly low value of 4.8 A2 for the main chain, four regions stand out as having significantly higher B-factors. Because thermal motion should be suppressed at cryogenic temperatures, the high B-factors are interpreted as the result of trapped conformational substates. The active site residues that are perturbed during accommodation of different substrates are precisely those showing conformational substates, implying that substrate binding selects a subset of conformations from the ensemble of accessible states. To better characterize the precise nature of these substates, a protein model consisting of 16 structures has been refined and evaluated. The model reveals a number of features that could not be well-described by conventional B-factors: for example, 40% of the main-chain residue conformations are distributed asymmetrically or in discrete clusters. Furthermore, these data demonstrate an unexpected correlation between motions on either side of the binding pocket that we suggest is a consequence of "dynamic close packing." These results provide strong evidence for the role of protein dynamics in substrate binding and are consistent with the results of dynamic studies of ligand binding in myoglobin and ribonuclease A.