Ligand binding to heme proteins: connection between dynamics and function

Structural factors controlling ligand binding to myoglobin: a kinetic hole-burning study

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as "kinetic hole burning." The study of these correlations gives direct insights into the structure-function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Electron transfer and protein dynamics in the photosynthetic reaction center

We have measured the kinetics of electron transfer (ET) from the primary quinone (Q(A)) to the special pair (P) of the reaction center (RC) complex from Rhodobacter sphaeroides as a function of temperature (5-300 K), illumination protocol (cooled in the dark and under illumination from 110, 160, 180, and 280 K), and warming rate (1.3 and 13 mK/s). The nonexponential kinetics are interpreted with a quantum-mechanical ET model (Fermi's golden rule and the spin-boson model), in which heterogeneity of the protein ensemble, relaxations, and fluctuations are cast into a single coordinate that relaxes monotonically and is sensitive to all types of relaxations caused by ET. Our analysis shows that the structural changes that occur in response to ET decrease the free energy gap between donor and acceptor states by 120 meV and decrease the electronic coupling between donor and acceptor states from 2.7 x 10(-4) cm(-1) to 1.8 x 10(-4) cm(-1). At cryogenic temperatures, conformational changes can be slowed or completely arrested, allowing us to monitor relaxations on the annealing time scale (approximately 10(3)-10(4) s) as well as the time scale of ET (approximately 100 ms). The relaxations occur within four broad tiers of conformational substates with average apparent Arrhenius activation enthalpies of 17, 50, 78, and 110 kJ/mol and preexponential factors of 10(13), 10(15), 10(21), and 10(25) s(-1), respectively. The parameterization provides a prediction of the time course of relaxations at all temperatures. At 300 K, relaxations are expected to occur from 1 ps to 1 ms, whereas at lower temperatures, even broader distributions of relaxation times are expected. The weak dependence of the ET rate on both temperature and protein conformation, together with the possibility of modeling heterogeneity and dynamics with a single conformational coordinate, make RC a useful model system for probing the dynamics of conformational changes in proteins.

Locally accessible conformations of proteins: multiple molecular dynamics simulations of crambin

Multiple molecular dynamics (MD) simulations of crambin with different initial atomic velocities are used to sample conformations in the vicinity of the native structure. Individual trajectories of length up to 5 ns sample only a fraction of the conformational distribution generated by ten independent 120 ps trajectories at 300 K. The backbone atom conformational space distribution is analyzed using principal components analysis (PCA). Four different major conformational regions are found. In general, a trajectory samples only one region and few transitions between the regions are observed. Consequently, the averages of structural and dynamic properties over the ten trajectories differ significantly from those obtained from individual trajectories. The nature of the conformational sampling has important consequences for the utilization of MD simulations for a wide range of problems, such as comparisons with X-ray or NMR data. The overall average structure is significantly closer to the X-ray structure than any of the individual trajectory average structures. The high frequency (less than 10 ps) atomic fluctuations from the ten trajectories tend to be similar, but the lower frequency (100 ps) motions are different. To improve conformational sampling in molecular dynamics simulations of proteins, as in nucleic acids, multiple trajectories with different initial conditions should be used rather than a single long trajectory.

Comparison of simulated and experimentally determined dynamics for a variant of the Lacl DNA-binding domain, Nlac-P

Recent advances in the experimentally determined structures and dynamics of the domains within LacI provide a rare context for evaluating dynamics calculations. A 1500-ps trajectory was simulated for a variant of the LacI DNA-binding domain, which consists of the first three helices in LacI and the hinge helix of the homologous PurR. Order parameters derived from dynamics simulations are compared to those obtained for the LacI DNA-binding domain with 15N relaxation NMR spectroscopy (Slijper et al., 1997. Biochemistry. 36:249-254). The MD simulations suggest that the unstructured loop between helices II and III does not exist in a discrete state under the conditions of no salt and neutral pH, but occupies a continuum of states between the DNA-bound and free structures. Simulations also indicate that the unstructured region between helix III and the hinge helix is very mobile, rendering motions of the hinge helix essentially independent of the rest of the protein. Finally, the alpha-helical hydrogen bonds in the hinge helix are broken after 1250 ps, perhaps as a prelude to helix unfolding.

Dynamic properties of monomeric insect erythrocruorin III from Chironomus thummi-thummi: relationships between structural flexibility and functional complexity

We have investigated the kinetics of geminate carbon monoxide binding to the monomeric component III of Chironomus thummi-thummi erythrocruorin, a protein that undergoes pH-induced conformational changes linked to a pronounced Bohr effect. Measurements were performed from cryogenic temperatures to room temperature in 75% glycerol and either 0.1 M potassium phosphate (pH 7) or 0.1 potassium borate (pH 9) after nanosecond laser photolysis. The distributions of the low temperature activation enthalpy g(H) for geminate ligand binding derived from the kinetic traces are quite narrow and are influenced by temperature both below and above approximately 170 K, the glass transition temperature. The thermal evolution of the CO binding kinetics between approximately 50 K and approximately 170 K indicates the presence of some degree of structural relaxation, even in this temperature range. Above approximately 220 K the width of the g(H) progressively decreases, and at 280 K geminate CO binding becomes exponential in time. Based on a comparison with analogous investigations of the homodimeric hemoglobin from Scapharca inaequivalvis, we propose a link between dynamic properties and functional complexity.

Contribution of Electron Paramagnetic Resonance to the studies of hemoglobin: the nitrosylhemoglobin system

Since the initial work of Ingram (8,10) Electron Paramagnetic Resonance contributed considerably to research in hemoglobins. Now, 40 years later we review some of the results of the application of EPR to nitrosylhemoglobin (HbNO), as an example of the diversity of information which this technique can provide.

Pressure effects on the proximal heme pocket in myoglobin probed by Raman and near-infrared absorption spectroscopy

The influence of high pressure on the heme protein conformation of myoglobin in different ligation states is studied using Raman spectroscopy over the temperature range from 30 to 295 K. Photostationary experiments monitoring the oxidation state marker bands demonstrate the change of rebinding rate with pressure. While frequency changes of vibrational modes associated with rigid bonds of the porphyrin ring are <1 cm(-1), we investigate a significant shift of the iron-histidine mode to higher frequency with increasing pressure (approximately 3 cm(-1) for deltaP = 190 MPa in Mb). The observed frequency shift is interpreted structurally as a conformational change affecting the tilt angle between the heme plane and the proximal histidine and the out-of-plane iron position. Independent evidence for iron motion comes from measurements of the redshift of band III in the near-infrared with pressure. This suggests that at high pressure the proximal heme pocket and the protein are altered toward the bound state conformation, which contributes to the rate increase for CO binding. Raman spectra of Mb and photodissociated MbCO measured at low temperature and variable pressure further support changes in protein conformation and are consistent with glasslike properties of myoglobin below 160 K.

Spectroscopic effects of polarity and hydration in the distal heme pocket of deoxymyoglobin

Distal pocket mutations at the E7 position (His64) of sperm whale deoxymyoglobin (deoxyMb) are used as a probe of distal pocket polarity and hydration. Changes of two key spectroscopic markers, the Fe-His(F8) stretch in the resonance Raman spectrum and the position of band III in the absorption spectrum, are monitored as the His64Tyr, His64Phe, His64Leu, and His64Gly mutations alter the distal heme pocket environment. The Fe-His vibration for the Phe, Leu, and Gly mutants is shifted to a lower frequency by 1-2 cm-1 relative to the Tyr mutant, wild type (WT), and native deoxyMb. Band III shifts to the red by approximately 4 nm ( approximately 70 cm-1) relative to WT and native deoxyMb for all the His64 mutants examined in this work. We correlate the small shift in the Fe-His frequency to the local electrostatic environment directly above the heme iron, affected by the presence of a localized water molecule in the heme pocket that is hydrogen-bonded to the E7 residue. The position of band III is roughly correlated to the displacement of the iron from the heme plane; however, the relatively large scatter in this correlation, as well as its dependence on distal pocket mutations, suggests that the heme pocket environment, particularly the E7 residue, also affects the energy of this transition.

On the origin of heme absorption band shifts and associated protein structural relaxation in myoglobin following flash photolysis

The role of the protein structural change monitored by absorption band shifts following flash photolysis of CO from myoglobin is discussed in terms of structure-function relationships. Evidence is presented that the Soret band shift does not depend primarily on the covalent linkage of the heme iron to the protein by using the mutation H93G(L) in which the proximal histidine 93 is replaced by glycine and an exogenous ligand L, which coordinates the heme iron but is not covalently bound to the globin. While CO rebinding kinetics depend strongly on the nature of the exogenous ligand L in H93G(L), the magnitude and time evolution of the Soret band shift in a viscous buffer on the nanosecond time scale are hardly perturbed in all cases studied. Comparison of the Soret band and charge transfer Band III shifts demonstrates that both have a similar time dependence on the nanosecond to microsecond time scale following flash photolysis in viscous solvents. We conclude that the nonexponential kinetics of protein relaxation probed by transient absorption band position shifts involves primarily distal coordinates prior to ligand escape. This result agrees with earlier measurements of Soret band shifts in distal pocket mutants of myoglobin (1). We suggest that the band shifts are primarily a response to changes in the electrostatic field around the heme (a transient Stark shift) associated with changes in protein structure that occur following ligand photodissociation.

Protein machine model of enzymatic reactions gated by enzyme internal dynamics

The slow character of conformational transition dynamics in native proteins, recently becoming more and more apparent, makes conventional theories of chemical reactions inapplicable for the description of enzymatic reactions. Any contemporary statistical theory of biochemical processes has to be based on a possibly simple but realistic model of microscopic dynamics of participating biomolecules. In a model considered in this paper the dynamics of enzymatic protein is approximated by a quasi-continuous diffusive motion of its solid-like structural elements relative to each other. The enzymatic reaction is assumed to involve three steps (a covalent tranformation preceded and followed by association-dissociation processes with the substrate and the product), each step being gated by conformational diffusion. In general, the reaction proceeds in three stages: initial, transient and steady-state. Carefully approximated analytical formulae describing the kinetics in each stage are derived. In the limit of the fast internal dynamics of the enzyme, when compared to the local chemical transformations, the initial stage of reaction, dependent on the initial distribution of enzyme conformations, is absent and all the formulae describing the remaining two stages simplify to those provided by the classical theory of Haldane. However, following recent studies, the rule seems to be that it is the conformational dynamics of the enzyme, and not the details of chemical mechanism, that affects the rate of enzymatic reaction. Apart from the possibility of the initial inhomogeneous kinetics, the important result obtained in the limit of slow conformational dynamics is that the kinetic mechanisms of a reaction differ in general between the transient and steady-state stages. Possibilities of carrying out an experimentum crucis directly discrediting the conventional approach are considered.

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.

Heme geometry in the 10 K photoproduct from sperm whale carbonmonoxymyoglobin

We have measured the Soret band of the photoproduct obtained by complete photolysis of sperm whale carbonmonoxymyoglobin at 10 K. The experimental spectrum has been modeled with an analytical expression that takes into account the homogeneous bandwidth, the coupling of the electronic transition with both high and low frequency vibrational modes, and the effects of static conformational heterogeneity. The comparison with deoxymyoglobin at low temperature reveals three main differences. In the photoproduct, the Soret band is shifted to red. The band is less asymmetric, and an enhanced coupling to the heme vibrational mode at 674 cm-1 is observed. These differences reflect incomplete relaxation of the active site after ligand dissociation. The smaller band asymmetry of the photoproduct can be explained by a smaller displacement of the iron atom from the mean porphyrin plane, in quantitative agreement with the X-ray structure analysis. The enhanced vibrational coupling is attributed to a subtle heme distortion from the planar geometry that is barely detectable in the X-ray structure.

Using buried water molecules to explore the energy landscape of proteins

Buried water molecules constitute a highly conserved, integral part of nearly all known protein structures. Such water molecules exchange with external solvent as a result of protein conformational fluctuations. We report here the results of water (17)O and (2)H magnetic relaxation dispersion measurements on wild-type and mutant bovine pancreatic trypsin inhibitor in aqueous solution at 4-80 degrees C. These data lead to the first determination of the exchange rate of a water molecule buried in a protein. The strong temperature dependence of this rate is ascribed to large-scale conformational fluctuations in an energy landscape with a statistical ruggedness of approximately 10 kJ mol(-1).

Protein tracking and detection of protein motion using atomic force microscopy

Height fluctuations over three different proteins, immunoglobulin G, urease, and microtubules, have been measured using an atomic force microscope (AFM) operating in fluid tapping mode. This was achieved by using a protein-tracking system, where the AFM tip was periodically repositioned above a single protein molecule (or structure) as thermal drifting occurred. Height (z-piezo signal) data were taken in 1 - or 2-s time slices with the tip over the molecule and compared to data taken on the support. The measured fluctuations were consistently higher when the tip was positioned over the protein, as opposed to the support the protein was adsorbed on. Similar measurements over patches of an amphiphile, where the noise was identical to that on the support, suggest that the noise increase is due to some intrinsic property of proteins and is not a result of different tip-sample interactions over soft samples. The orientation of the adsorbed proteins in these preliminary studies was not known; thus it was not possible to make correlations between the observed motion and specific protein structure or protein function beyond noting that the observed height fluctuations were greater for an antibody (anti-bovine IgG) and an enzyme (urease) than for microtubules.

Two-dimensional distributions of activation enthalpy and entropy from kinetics by the maximum entropy method

The maximum entropy method (MEM) is used to numerically invert the kinetics of ligand rebinding at low temperatures to obtain the underlying two-dimensional distribution of activation enthalpies and entropies, g(H,S). A global analysis of the rebinding of carbon monoxide (CO) to myoglobin (Mb), monitored in the Soret band at temperatures from 60 to 150 K, is performed using a Newton-Raphson optimization algorithm. The MEM approach describes the data much better than traditional least-squares analyses, reducing chi 2 by an order of magnitude. The MEM resolves two barrier distributions suggestive of rebinding to different bound conformations of MbCO, the so-called A1 and A3 substates, whose activation barriers have been independently estimated from kinetics monitored in the infrared. The distribution corresponding to A3 possesses higher activation entropies, also consistent with infrared measurements. Within an A substate, correlations of S and H are recovered qualitatively from simulated data but can be difficult to obtain from experimental data. When the rebinding measured at 60 K is excluded from the inversion, two peaks are no longer clearly resolved. Thus, data of very high quality are required to unambiguously determine the kinetic resolvability of subpopulations and the shape of the barrier distribution for a single A substate.

Thermal access to amplified chemical potential and the determination of equilibrium constants in protein solutions at subfreezing temperatures

During rapid cooling of ferric heme protein solutions containing fluoride, locally concentrated ligand cannot fully equilibrate with heme before the temperature drops below 200 K and into the range where energy is insufficient for exchange with iron-bound water. When temperature is then jumped above 200 K, exchange of fluoride for bound water is activated. Between 200 and 240 K, further fluoride complex formation takes place over several minutes; its extent is measured along the kinetic curve by reimmersing the sample into liquid nitrogen and taking EPR spectra. Kinetic curves for replacement of iron-bound water by fluoride in horse aquo-ferrimyoglobin and human aquo-ferrihemoglobin, and corresponding equilibrium constants have been obtained at temperatures between 200 and 240 K. The reaction rates are affected by sucrose. Results indicate that the kinetics of exchange of fluoride for heme-bound water at subfreezing temperatures is protein specific and not diffusion-controlled, and is not affected by the phase transition of ice which takes place at subfreezing temperature. Free energy changes accompanying these reactions are largely continuous as the systems pass from above to below freezing.

Kinetic evidence for the existence of a rate-limiting step in the reaction of ferric hemoproteins with anionic ligands

The kinetics of azide and fluroide binding to various monomeric and tetrameric ferric hemoproteins (sperm whale Mb, isolated alpha and beta chains of human Hb reacted with p-chloromercuribenzoate, dromeday, ox and human Hb) has been investigated (at pH 6.5 and 20 degrees C over a large range (20 microM to 2 M) of ligand concentration. It has been observed that the pseuo-first-order rate constant for azide binding to the hemoproteins investigated does not increase linearly with ligand concentration, but tends to level off toward an asymptomatic concentration-independent value typical for each hemoprotein. This behavior, which has been detected only by an investigation covering an unusually large range of ligand concentrations appears to be independent of the ionic strength, and it underlies the existence of a rate-limiting step in the dynamic pathway of azide binding to ferric hemoproteins, which is detectable whenever the observed pseudo- first-order rate constant becomes faster than a given value characteristic of the specific hemoprotein. Such a behavior is not observed in the case of fluroide binding probably because the pesudo- first-order rate constant for this ligand is much slower and never attains a value faster than that of the rate-limiting step. In general terms, this feature should involve a conformational equilibrium between at least two forms (possibly related to the interaction of H2O with distal histidine and its exchange with the bulk solvent) which modulates the access of the anionic ligand into the heme pocket and its reaction with the ferric iron.

Ligand binding and protein dynamics in cupredoxins

Type 1 copper sites bind nitric oxide (NO) in a photolabile complex. We have studied the NO binding properties of the type 1 copper sites in two cupredoxins, azurin and halocyanin, by measuring the temperature dependence of the ligand binding equilibria and the kinetics of the association reaction after photodissociation over a wide range of temperature (80-280 K) and time (10(-6)-10(2) s). In both proteins, we find nonexponential kinetics below 200 K that do not depend on the NO concentration. Consequently, this process is interpreted as geminate recombination. In azurin, the rebinding can be modeled with the Arrhenius law using a single pre-exponential factor of 10(8.3) s-1 and a Gaussian distribution of enthalpy barriers centered at 22 kJ/mol with a width [full width at half-maximum (FWHM)] of 11 kJ/mol. In halocyanin, a more complex behavior is observed. About 97% of the rebinding population can also be characterized by a Gaussian distribution of enthalpy barriers at 12 kJ/mol with a width of 6.0 kJ/mol (FWHM). The pre-exponential of this population is 1.6 x 10(12) s-1 at 100 K. After the majority population has rebound, a power-law phase that can be modeled with a gamma-distribution of enthalpy barriers is observed. Between 120 and 180 K, an additional feature that can be interpreted as a relaxation of the barrier distribution toward higher barriers shows up in the kinetics. Above 200 K, a slower, exponential rebinding appears in both cupredoxins. Since the kinetics depend on the NO concentration, this process is identified as bimolecular rebinding.(ABSTRACT TRUNCATED AT 250 WORDS)

Protein reaction kinetics in a room-temperature glass

Protein reaction kinetics in aqueous solution at room temperature are often simplified by the thermal averaging of conformational substates. These substates exhibit widely varying reaction rates that are usually exposed by trapping in a glass at low temperature. Here, it is shown that the solvent viscosity, rather than the low temperature, is primarily responsible for the trapping. This was demonstrated by placement of myoglobin in a glass at room temperature and subsequent observation of inhomogeneous reaction kinetics. The high solvent viscosity slowed the rate of crossing the energy barriers that separated the substates and also suppressed any change in the average protein conformation after ligand dissociation.

Automatic identification of discrete substates in proteins: singular value decomposition analysis of time-averaged crystallographic refinements

The singular value decomposition (SVD) provides a method for decomposing a molecular dynamics trajectory into fundamental modes of atomic motion. The right singular vectors are projections of the protein conformations onto these modes showing the protein motion in a generalized low-dimensional basis. Statistical analysis of the right singular vectors can be used to classify discrete configurational substates in the protein. The configuration space portraits formed from the right singular vectors can also be used to visualize complex high-dimensional motion and to examine the extent of configuration space sampling by the simulation.

Thermal fluctuations between conformational substates of the Fe(2+)-HisF8 linkage in deoxymyoglobin probed by the Raman active Fe-N epsilon (HisF8) stretching vibration

We have measured the VFe-His Raman band of horse heart deoxymyoglobin dissolved in an aqueous solution as a function of temperature between 10 and 300 K. The minimal model to which these data can be fitted in a statistically significant and physically meaningful way comprises four different Lorentzian bands with frequencies at 197, 209, 218, and 226 cm-1, and a Gaussian band at 240 cm-1, which exhibit halfwidths between 10 and 12.5 cm-1. All these parameters were assumed to be independent of temperature. The temperature dependence of the apparent total band shape's frequency is attributed to an intensity redistribution of the subbands at omega 1 = 209 cm-1, omega 2 = 218 cm-1, and omega 3 = 226 cm-1, which are assigned to Fe-N epsilon (HisF8) stretching modes in different conformational substrates of the Fe-HisF8 linkage. They comprise different out-of-plane displacements of the heme iron. The two remaining bands at 197 and 240 cm-1 result from porphyrin modes. Their intensity ratio is nearly temperature independent. The intensity ratio I3/I2 of the vFe-His subbands exhibits a van't Hoff behavior between 150 and 300 K, bending over in a region between 150 and 80 K, and remains constant between 80 and 10 K, whereas I2/I1 shows a maximum at 170 K and approaches a constant value at 80 K. These data can be fitted by a modified van't Hoff expression, which accounts for the freezing into a non-equilibrium distribution of substates below a distinct temperature Tf and also for the linear temperature dependence of the specific heat of proteins. The latter leads to a temperature dependence of the entropic and enthalpic differences between conformational substates. The fits to the intensity ratios of the vFe-His subbands yield a freezing temperature of Tf = 117 K and a transition region of delta T = 55 K. In comparison we have utilized the above thermodynamic model to reanalyze earlier data on the temperature dependence of the ratio Ao/A1 of two subbands underlying the infrared absorption band of the CO stretching vibration in CO-ligated myoglobin (A. Ansari, J. Berendzen, D. Braunstein, B. R. Cowen, H. Frauenfelder, M. K. Kong, I. E. T. Iben, J. Johnson, P. Ormos, T. B. Sauke, R. Scholl, A. Schulte, P. J. Steinbach, R. D. Vittitow, and R. D. Young, 1987, Biophys. Chem. 26:237-335). This yields thermodynamic parameters, in particular the freezing temperature (Tf = 231 K) and the width of the transition region (AT =8 K), which are significantly different from the corresponding parameters obtained from the above vFe-His data, but very close to values describing the transition of protein bound water from a liquid into an amorphous state. These findings and earlier reported data on the temperature dependence exhibited by the Soret absorption bands of various deoxy and carbonmonoxymyoglobins led us to the conclusion that the fluctuations between conformational substates of the heme environment in carbonmonoxymyoglobin are strongly coupled to motions within the hydration shell, whereas the thermal motions between the substates of the Fe-HisF8 linkage in deoxymyoglobin proceed on an energy landscape that is mainly determined by the intrinsic properties of the protein. The latter differ from protein fluctuations monitored by Mossbauer experiments ondeoxymyoglobin crystals which exhibit a strong coupling to the protein bound water and most probably reflect a higher tier in the hierarchical arrangement of substates and equilibrium fluctuations.

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.

Mössbauer spectroscopy on nonequilibrium states of myoglobin: a study of r-t relaxation

A frozen solution of 57Fe-enriched metmyoglobin was irradiated by x rays at 77 K. Mössbauer spectra showed a reduction of Fe(III) high spin by thermalized electrons and a production of a metastable Fe(II) low spin myoglobin complex with H2O at its sixth coordination site. The relaxation of the intermediate was investigated by Mössbauer spectroscopy as a function of temperature and time. The relaxation process starts above 140 K and is fully completed at approximately 200 K. At temperatures between 140 and 200 K, the relaxation lasts for hours and is nonexponential in time. Up to 180 K, the process can be described satisfactorily by a gamma distribution of activation enthalpies with an Arrhenius relation for the rate coefficient. The temperature and time dependence of the Mössbauer parameters indicates structural changes in the active center of the protein as early as 109 K that continue for several hours at higher temperatures. Above 180 K, structural rearrangements involving the whole protein molecule lead to a shift and narrowing of the barrier height distribution.