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.     

Structural fluctuations of myoglobin from normal-modes,M?ssbauer,Raman, and absorption spectroscopy.

A normal-mode analysis of carbon monoxymyoglobin (MbCO) and deoxymyoglobin (Mb) with 170 water molecules is performed for (54)Fe and (57)Fe. A projection is defined that extracts iron out-of-plane vibrational modes and is used to calculate spectra that can be compared with those from resonance Raman scattering. The calculated spectra and the isotopic shift (57)Fe versus (54)Fe agree with the experimental data. At low temperatures the average mean square fluctuations (MSFs) of the protein backbone atoms agree with molecular dynamics simulation. Below 180 K the MSFs of the heme iron agree with the data from Mossbauer spectroscopy. The MSFs of the iron atom relative to the heme are an order of magnitude smaller than the total MSFs of the iron atom. They agree with the data from optical absorption spectroscopy. Thus the MSFs of the iron atom as measured by Mossbauer spectroscopy can be used to probe the overall motion of the heme within the protein matrix, whereas the Gaussian thermal line broadening of the Soret band and the resonance Raman bands can be used to detect local intramolecular iron-porphyrin motions.     

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.     

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.     

Molecular dynamics simulations of Hydrogenobacter thermophilus cytochrome c552: comparisons of the wild-type protein, a b-type variant, and the apo state

Molecular dynamic simulations have been performed for wild-type Hydrogenobacter thermophilus cytochrome c(552), a b-type variant of the protein, and the apo state with the heme prosthetic group removed. In the b-type variant, Cys 10 and Cys 13 were mutated to alanine residues, and so the heme group was no longer covalently bound to the protein. Two 8-ns simulations have been performed for each system at 298 and 360 K. The simulations of the wild-type protein at 298 K show a very close agreement with experimental NMR data. A fluxional process involving the side chain of Met 59, which coordinates to the heme iron, is observed in accord with proposals from NMR studies. Overall, the structure and dynamical behavior of the protein during the simulations of the b-type variant is closely similar to that of the wild-type protein. However, side chains in the heme-binding site show larger fluctuations in the b-type variant simulation at 360 K. In addition, structural changes are seen for a number of residues close to the heme group, particularly Gly 22 and Ser 51. The simulations of the apo state show significant conformational changes for residues 50-59. These residues form a loop region, which packs over the heme group in the wild-type protein and hydrogen bonds to the heme propionate groups. In the absence of heme, in the apo state simulations, these residues form short but persistent regions of beta-sheet secondary structure. These could provide nucleation sites for the conversion to amyloid fibrils.     

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.     

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.     

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.     

Investigations of optical line shapes and kinetic hole burning in myoglobin.

We present the results of an extensive investigation of the optical line shapes of deoxymyoglobin (Mb), the ligand-bound form (MbCO), and the low-temperature photoproduct (Mb*). The thermal properties and the pH dependence of the Soret band and the near infrared band III (approximately 760 nm) are analyzed, taking into account the underlying vibrational properties of the absorption bands. The strong temperature dependence associated with the Soret band of MbCO and band III of Mb indicates significant coupling to low-frequency modes that may not be directly observed in the resonance Raman spectra. On the basis of analogous line-shape studies in a variety of heme systems, we assign the low-frequency coupling in MbCO to torsional motions of the CO molecule. The low-frequency mode coupled to band III (approximately 70 cm-1) is found to lie quite close to the value for the heme-doming motion (approximately 50 cm-1) calculated by using the kinetically determined value of the force constant (17 N/m). Significant inhomogeneous broadening in the Soret region of Mb and Mb* is found to be due to a "nonkinetic" coordinate that we associate with the orientation of the proximal histidine. A "kinetic" coordinate, associated with the equilibrium displacement of the iron atom from the porphyrin plane (a) is found to contribute to the inhomogeneous broadening of both the Soret band and band III. The relaxation of the heme as the system evolves from from Mb* to Mb is followed optically as a function of temperature, and a sharp transition temperature is found at 185 K. The blue shifts of the Soret band and band III as Mb* evolves to Mb are found to be nearly identical (delta v*ABS approximately 140 cm-1) and attributed to changes in the mean value of a between Mb* (a*0) and Mb (a0 = 0.45 A). A simple quadratic model for the coordinate coupling that simultaneously accounts for the observed shift, delta v*ABS, the low-temperature kinetics and the kinetic hole burning predicts a*0 = 0.2 +/- 0.05 A and EA = 16 +/- 2 kJ/mol for the room temperature Arrhenius barrier height at the heme. A simple quantitative method for the analysis of kinetic hole-burning experiments is also developed and applied to recent studies involving quaternary and subunit-specific hemoglobin structures.     

A molecular dynamics simulation of bacteriophage T4 lysozyme.

An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. The simulation was carried out with all hydrogen atoms modeled explicitly, the inclusion of all 152 crystallographic waters and at a temperature of 300 K. Temporal analysis of the trajectory versus energy, hydrogen bond stability, r.m.s. deviation from the starting crystal structure and radius of gyration, demonstrates that the simulation was both stable and representative of the average experimental structure. Average structural properties were calculated from the enzyme trajectory and compared with the crystal structure. The mean value of the C alpha displacements of the average simulated structure from the X-ray structure was 1.1 +/- 0.1 A; differences of the backbone phi and psi angles between the average simulated structure and the crystal structure were also examined. Thermal-B factors were calculated from the simulation for heavy and backbone atoms and both were in good agreement with experimental values. Relationships between protein secondary structure elements and internal motions were studied by examining the positional fluctuations of individual helix, sheet and turn structures. The structural integrity in the secondary structure units was preserved throughout the simulation; however, the A helix did show some unusually high atomic fluctuations. The largest backbone atom r.m.s. fluctuations were found in non-secondary structure regions; similar results were observed for r.m.s. fluctuations of non-secondary structure phi and psi angles. In general, the calculated values of r.m.s. fluctuations were quite small for the secondary structure elements. In contrast, surface loops and turns exhibited much larger values, being able to sample larger regions of conformational space. The C alpha difference distance matrix and super-positioning analyses comparing the X-ray structure with the average dynamics structure suggest that a 'hinge-bending' motion occurs between the N- and C-terminal domains.     

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.     

Spectroscopic studies of myoglobin at low pH: heme structure and ligation

We explore heme structure and ligation subsequent to a low-pH conformational transition in sperm whale myoglobin. Below pH 4.0, the iron-histidine bond breaks in metMb and deoxyMb. In MbCO, the majority of the iron-histidine bonds remain intact down to pH 2.6; however, the observation of a weak Fe-CO mode at 526 cm-1 indicates that a small fraction of the sample has the histidine replaced by a weak ligand, possibly water. The existence of a sterically hindered CO subpopulation in MbCO and the continued association of the four-coordinate heme with the protein in deoxyMb suggest that the heme pocket remains at least partially intact in the acid-induced conformation. The global pH-dependent conformational change described here is clearly distinguished from the local "closed" to "open" transition described previously in MbCO [Morikis et al. (1989) Biochemistry 28, 4791-4800]. Further observations of the four-coordinate heme state yield insights on the mechanism of heme photoreduction and the assignment of the 760-nm band in deoxyMb.     

Low-temperature formation of a distal histidine complex in hemoglobin: a probe for heme pocket flexibility

Pocket dynamics of horse deoxyhemoglobin and methemoglobin in the temperature range from 80 to 260 K is investigated. In both hemoglobins reversible conversion to a low-spin iron complex is observed at temperatures as low as 210 K. Electron spin resonance (ESR) and M?ssbauer data assigned this low-spin iron complex to the coordination of N tau-His-E7 as a sixth nitrogenous ligand. The bonding of this ligand located 4 A from the iron indicates the presence of a thermally available conformation that exhibits a high degree of flexibility in the heme pocket. In deoxyhemoglobin, the formation of the bis(histidine) complex was accompanied by excitations of conformational fluctuations manifested through the temperature dependence of the M?ssbauer-Lamb factor. The rate for the formation of this complex, with an associated energy barrier (greater than 60 KJ mol-1), is shown to serve as an index of heme pocket flexibility. Measurements performed on partially liganded (carbonmonoxy) hemoglobin indicate that partial ligation enhances conversion of the unliganded subunits to the bis(histidine) complex, suggesting that pocket dynamics is affected by subunit interactions.     

Analysis of the dynamics of rhizomucor miehei lipase at different temperatures.

The dynamics of Rhizomucor miehei lipase has been studied by molecular dynamics simulations at temperatures ranging from 200-500K. Simulations carried out in periodic boundary conditions and using explicit water molecules were performed for 400 ps at each temperature. Our results indicate that conformational changes and internal motions in the protein are significantly influenced by the temperature increase. With increasing temperature, the number of internal hydrogen bonds decreases, while surface accessibility, radius of gyration and the number of residues in random coil conformation increase. In the temperature range studied, the motions can be described in a low dimensional subspace, whose dimensionality decreases with increasing temperature. Approximately 80% of the total motion is described by the first (i) 80 eigenvectors at T=200K, (ii) 30 eigenvectors at T=300K and (iii) 10 eigenvectors at T=400K. At high temperature, the alpha-helix covering the active site in the native Rhizomucor miehei lipase, the helix at which end the active site is located, and in particular, the loop (Gly35-Lys50) show extensive flexibility.     

Structural dynamics of liganded myoglobin.

X-ray crystallography can reveal the magnitudes and principal directions of the mean-square displacements of every atom in a protein. This structural information is complementary to the temporal information obtainable by spectroscopic techniques such as nuclear magnetic resonance. Determination of the temperature dependence of the mean-square displacements makes it possible to separate large conformational motions from simple thermal vibrations. The contribution of crystal lattice disorder to the overall apparent displacement can be estimated by Mössbauer spectroscopy. This technique has been applied to high resolution x-ray diffraction data from sperm whale myoglobin in its Met iron and oxy cobalt forms. Both crystal structures display regions of large conformational motions, particularly at the chain termini and in the region of the proximal histidine. Overall, the mean-square displacement increases with increasing distance from the center of gravity of the molecule. Some regions of the heme pocket in oxy cobalt myoglobin are more rigid than the corresponding regions in Met myoglobin.     

Infrared Study of Carbon Monoxide Migration among Internal Cavities of Myoglobin Mutant L29W

Myoglobin, a small globular heme protein that binds gaseous ligands such asO₂, CO and NO reversibly at the heme iron, provides an excellent modelsystem for studying structural and dynamic aspects of protein reactions. Flashphotolysis experiments, performed over wide ranges in time and temperature, reveal a complex ligand binding reaction with multiple kinetic intermediates, resulting from protein relaxation and movements of the ligand within the protein. Our recent studies of carbonmonoxy-myoglobin (MbCO) mutant L29W, using time-resolved infrared spectroscopy in combination with x-ray crystallography, have correlated kinetic intermediates with photoproduct structures that are characterized by the CO residing in different internal protein cavities, so-called xenon holes. Here we have used Fourier transform infrared temperature derivative spectroscopy (FTIR-TDS) to further examine the role of internal cavities in the dynamics. Different cavities can be accessed by the CO ligands at different temperatures, and characteristic infrared absorption spectra have been obtained for the different locations of the CO ligand within the protein, enabling us to monitor ligand migration through the protein as well as conformational changes of the protein.     

Dynamics of hemoglobin investigated by M?ssbauer spectroscopy

M?ssbauer Spectra of Fe enriched horse hemoglobin and sperm whale myoglobin were measured in the temperature range from 80 K to 260 K. An analysis of the temperature dependence of the recoiless fraction (the Lamb-M?ssbauer factor) shows it to be sensitive to conformational fluctuations which affect the mean square displacement of the iron. We have found that the protein conformation has a dramatic effect on these measurements. For hemoglobin greater conformational fluctuations at lower temperatures are observed for carbonmonoxyhemoglobin in the liganded conformation than for deoxyhemoglobin in the unliganded conformation. On the other hand, the Lamb-M?ssbauer factor is insensitive to the binding of ligands to myoglobin and shows conformational fluctuations similar to deoxyhemoglobin even in the liganded state. It is also shown that a reversible complex with the distal histidine is formed in frozen deoxyhemoglobin solution above 200 K where the Lamb-M?ssbauer factor shows the excitation of new modes of conformational fluctuations. This complex is not formed with carbonmonoxyhemoglobin which already has a sixth ligand and with deoxymyoglobin which appears to undergo much more limited conformational fluctuations. A possible relationship between the formation of the distal histidine complex and the cooperative ligand binding reaction is suggested by results with partially liganded hemoglobin which indicate increased formation of the distal histidine complex.     

Evidence of Functional Protein Dynamics from X-Ray Crystallographic Ensembles

It is widely recognized that representing a protein as a single static conformation is inadequate to describe the dynamics essential to the performance of its biological function. We contrast the amino acid displacements below and above the protein dynamical transition temperature, T_D∼215K, of hen egg white lysozyme using X-ray crystallography ensembles that are analyzed by molecular dynamics simulations as a function of temperature. We show that measuring structural variations across an ensemble of X-ray derived models captures the activation of conformational states that are of functional importance just above T_D, and they remain virtually identical to structural motions measured at 300K. Our results highlight the ability to observe functional structural variations across an ensemble of X-ray crystallographic data, and that residue fluctuations measured in MD simulations at room temperature are in quantitative agreement with the experimental observable.     

Dynamics of Hemoglobin Investigated by Mössbauer Spectroscopy

Abstract

Mössbauer Spectra of ⁵⁷Fe enriched horse hemoglobin and sperm whale myoglobin were measured in the temperature range from 80 K to 260 K. An analysis of the temperature dependence of the recoiless fraction (the Lamb-Mössbauer factor) shows it to be sensitive to conformational fluctuations which affect the mean square displacement of the iron. We have found that the protein conformation has a dramatic effect on these measurements. For hemoglobin greater conformational fluctuations at lower temperatures are observed for carbonmonoxyhemoglobin in the liganded conformation than for deoxyhemoglobin in the unliganded conformation. On the other hand, the Lamb-Mössbauer factor is insensitive to the binding of ligands to myoglobin and shows conformational fluctuations similar to deoxyhemoglobin even in the liganded state. It is also shown that a reversible complex with the distal histidine is formed in frozen deoxyhemoglobin solutions above 200 K where the Lamb-Mössbauer factor shows the excitation of new modes of conformational fluctuations. This complex is not formed with carbonmonoxyhemoglobin which already has a sixth ligand and with deoxymyoglobin which appears to undergo much more limited conformational fluctuations. A possible relationship between the formation of the distal histidine complex and the cooperative ligand binding reaction is suggested by results with partially liganded hemoglobin which indicate increased formation of the distal histidine complex.     

Light-induced protein-matrix uncoupling and protein relaxation in dry samples of trehalose-coated MbCO at room temperature

In humid samples of trehalose-coated carboxy-myoglobin (MbCO), thermally driven conformational relaxation takes place after photodissociation of the carbon monoxide (CO) molecule at room temperature. In such samples, because of the extreme viscosity of the external matrix, photodissociated CO cannot diffuse out of the protein and explores the whole (proximal and distal side) heme pocket, experiencing averaged protein heme pocket structures, as a results of the presence of Brownian motions. At variance, in very dry samples, a lower portion of the photodissociated CO diffuses from the distal to the proximal heme pocket side probing in nonaveraged structures. We revisit here the flash photolysis data by Librizzi et al. (2002) and report on new, room temperature experiments in MbCO-trehalose samples, shortly illuminated prior the laser pulse. In dry samples, pre-illumination increased the diffusion of CO from the distal to the proximal heme pocket side, which resulted in less structure than in non-pre-illuminated samples. Such an effect, which is absent in humid samples, stems from a decoupling of the protein internal degrees of freedom from those of the external water-sugar matrix. We suggest that such a decoupling can be brought about by the continuous attempts performed by the protein during pre-illumination to undergo relaxation toward the photodissociated deoxy state. This, in turn, causes a collapse in the hydrogen bond network, which connects the protein surface to the water-sugar matrix, as reported by Cottone et al. (2002) and Giuffrida et al. (2003). In the conclusion section, we discuss the possible involvement of the processes invoked to rationalize the present data, in the function of macromolecules and interactions in living cells.