**2. Material and methods**

Mb, Hb, inositol hexakisphosphate (IHP), and bezafibrate (BZF) were purchased from Sigma-Aldrich and used as received. Fe(III) tetrakis(4-sulfonatophenyl)porphine (Fe(III)4SP) was obtained from Frontier-Scientific Inc. Oxymyoglobin (O2-Mb) samples were prepared by dissolving the protein in 50 mM HEPES buffer pH 7.0. The sample was then purged with Ar for 10 min and reduced by addition of a freshly prepared solution of sodium dithionite. The quality of the deoxymyoglobin (deoxyMb) was verified by UV-visible spectroscopy. (O2-Mb) was obtained by bubbling air through deoxyMb sample. The CO bound hemoglobin sample was prepared by desolving Hb in 100 mM HEPES buffer pH 7.0 in a 0.5 x 1cm quartz cuvette. The concentration of allosteric effectors was 5 mM for BZF and 1 mM for IHP. The sample was then sealed with a septum cap and purged with Ar for 10 min, reduced with a small amount of sodium dithionite to prepare deoxyhemoglobin (deoxyHb), and subsequently bubbled with CO for approximately 1 min. Preparation of O2-Mb and CO-Hb aducts was checked by UV-vis spectroscopy (Cary50, Varian).

#### **2.1 Quantum yield determination**

The quantum yield () was determined as described previously (Belogortseva et al., 2007). All transient absorption measurements were carried out on 50 µM samples in 50 or 100 mM HEPES buffer, pH 7.0, placed in a 2 mm path quartz cell. The cell was placed into a temperature controlled holder (Quantum Northwest) and the ligand photo-dissociation was triggered using a 532 nm output from a Nd:YAG laser (Minilite II, Continuum). The probe beam, an output from the Xe arc lamp (200 W, Newport) was propagated through the center of the cell and then focused on the input of a monochromator (Yvon-Jovin ). The intensity of the probe beam was detected using an amplified photodiode (PDA 10A, Thornlabs) and subsequently digitized (Wave Surfer 42Xs, 400 MHz). The power of the pump beam was kept below 50 µJ to match the laser power used in photoacoustic measurements. The quantum yield was determined by comparing the change in the sample absorbance at 440 nm with that of the reference, CO bound myoglobin of known quantum yield (ref= 0.96 (Henry et al., 1983)) according to Eq 1:

$$
\Phi = \frac{\Delta \mathbf{A}\_{\text{sam}} \Delta \varepsilon\_{\text{ref}} \Phi\_{\text{ref}}}{\Delta \mathbf{A}\_{\text{ref}} \Delta \varepsilon\_{\text{sam}}} \tag{1}
$$

where ΔAsam and ΔAref are the absorbance change of the sample and reference at 440 nm, respectively, and Δsam and Δref are the change in the extinction coefficient between the CO bound and reduced form of the sample and the reference, respectively.

Time Resolved Thermodynamics Associated with Diatomic Ligand Dissociation from Globins 305

�� <sup>+</sup> ���� (�����)

the first and subsequent step of the reaction, respectively. To retrieve thermodynamic and kinetic parameters, the reference trace is convoluted with the H(t) function using estimated parameters (i and i ) and the calculated E(t)calc is compared with the E(t)obs. The i and <sup>i</sup> values are varied until a satisfactory fit is obtained in terms of 2 and autocorrelation function. In practice, the lifetime for the prompt process is fixed to 1 ns, whereas other

For processes that occur with a quantum yield that is temperature dependent in the temperature range used in PAC measurements, the thermodynamic parameters for the fast phase (<20ns) are determined by plotting [Ehν(1-)]/] versus (Cpρ/β) according to Eq. 7 and the volume and enthalpy changes for the subsequent steps are obtained by plotting

� � � ��� + �� ���

� � � ��� + ��(���

Ligand migration in heme proteins is often described using the sequential three-state model (Henry et al., 1983) shown in Scheme 2. Upon cleavage of the coordination bond between the ligand and heme iron, the ligand is temporarily trapped within the protein matrix and then it either directly rebinds back to the heme iron in the so called "geminate rebinding" or diffuses from the protein matrix into the surrounding solvent. The subsequent bimolecular ligand binding to heme iron occurs on significantly longer time scales, hundreds of microseconds to milliseconds. The quantum yield for the geminate rebinding and for bimolecular association is strongly dependent on the character of the ligand and the protein. For example, CO rebinds to Mb predominantly through a bimolecular reaction with quantum yield close to unity (bim = 0.96 )(Henry et al., 1983), whereas the quantum yield for bimolecular O2 rebinding to heme proteins is significantly lower (Carver et al., 1990; Walda et al., 1994), and NO rebinds predominantly through geminate rebinding (Ye et al., 2002). To determine the thermodynamic parameters from acoustic data, the quantum yields for CO and O2 bimolecular rebinding to Hb and Mb, respectively, have to be known. The quantum yield for O2 binding to Mb was measured in the temperature range from 5 C to 35 C (Fig. 2) and the values show a weak temperature dependence with the quantum yield decreasing with increasing temperature. At 20 C the quantum yield is 0.09 ± 0.01 that is within the range of values reported previously ( = 0.057 (Walda et al., 1994) and ( = 0.12 (Carver et al., 1990)). We have also measured the quantum yield for CO bimolecular rebinding to Hb, and to Hb in the presence of effector molecules (Fig. 2). The quantum yield increases linearly with temperature and at 20 C, CO binds to Hb with quantum yield of

where ΔH and ΔV correspond to the reaction enthalpy and volume change, respectively.

�� � �� � � �

�(�)��� � �(�)��(�) (6)

) term in Eq. 4 and the 1 and 2 are the lifetime for

�� (5)

� (7)

� ) (8)

�(�) ������ �

����

(Ehν/) versus (Cpρ/β) according to Eq. 8 (Peters et al., 1992).

���(���)

����

where 1 and 2 correspond to the (����

parameters are allowed to be varied.

**3. Results** 

### **2.2 Photoacoustic calorimetry**

The photo-acoustic set-up used in our lab was described previously (Miksovska&Larsen, 2003). Briefly, the sample in a quartz cell was placed in a temperature controlled holder (Quantum Northwest). The 532 nm output from a Nd:YAG laser (7 ns pulse width , < 50 µJ power) was shaped using a narrow slit (100 µm) and focused on the center of the quartz cell. An acoustic detector (1 MHz, RV103, Panametrix) was coupled to the side of a quartz cell using a thin layer of honey and the detector output was amplified using an ultrasonic preamplifier (Panametrics 5662). The signal was then stored in a digitizer (Wave Surfer 42Xs, 400 MHz). The data were analyzed using Sound Analysis software (Quantum Northwest).

### **2.3 Data analysis**

The excitation of the photocleavable iron-ligand bond in heme proteins generates at least two processes that contribute to the photoacoustic wave: the volume change due to the heat released during the reaction (Q), and the volume change (ΔV') due to the photo-triggered structural changes (including bond cleavage / formation, electrostriction, solvation, etc.). The amplitude of the sample acoustic wave (Asam) can be expressed as:

$$A\_{sam} = K E\_a (Q \frac{\beta}{\varrho c\_p} + \Delta V') \tag{2}$$

where K is the instrument response constant, Ea is number of Einsteins absorbed, β is the expansion coefficient, ρ is the density, and Cp is the heat capacity. For water, the (β/Cpρ) term strongly varies with temperature mainly due to the temperature dependence of the β term. To evaluate the instrument response constant, the photo-acoustic traces are measured for a reference compound under experimental conditions (laser power, temperature, etc.) identical to those for the sample measurements. We have used Fe(III)4SP as a reference since it is non-fluorescent and photo-chemically stable. The amplitude of the reference acoustic trace can be described as:

$$A\_{ref} = K E\_a E\_{h\nu} \frac{\beta}{\varrho c\_p} \tag{3}$$

where Ehν is the energy of a photon at 532 nm (Ehν= 53.7 kcal mol-1). The amount of heat deposited to the solvent and the non-thermal volume changes can then be determined by measuring the acoustic wave for the sample and the reference for several temperatures and plotting the ratio of the sample and reference acoustic wave () as a function of (Cpρ/β) according to Eq. 4:

$$\frac{A\_{sam}}{A\_{ref}} = \left. \phi \right| = Q + \Delta V' \frac{\rho c\_p}{\beta} \tag{4}$$

For a multi-step process that exhibits volume and enthalpy changes on the time-scale between ~ 20 ns to 5 µs, the thermodynamic parameters for each individual step and corresponding lifetimes can be determined due to the sensitivity of the acoustic detector to the temporal profile of the pressure change. The time dependent sample acoustic signal E(t)obs can be expressed as a convolution of the time dependent function describing the volume change H(t) with the instrument response T(t) function (the reference acoustic wave):

$$H(t) = \phi\_1 e^{-\frac{t}{\tau 1}} + \frac{\phi\_2 k\_2}{(k\_2 - k\_1)} e^{-\frac{t}{\tau 1}} - e^{\frac{t}{\tau 2}} \tag{5}$$

$$E\left(\mathbf{t}\right)\_{obs} = H(\mathbf{t}) \otimes T(\mathbf{t})\tag{6}$$

where 1 and 2 correspond to the (���� ���� ) term in Eq. 4 and the 1 and 2 are the lifetime for the first and subsequent step of the reaction, respectively. To retrieve thermodynamic and kinetic parameters, the reference trace is convoluted with the H(t) function using estimated parameters (i and i ) and the calculated E(t)calc is compared with the E(t)obs. The i and <sup>i</sup> values are varied until a satisfactory fit is obtained in terms of 2 and autocorrelation function. In practice, the lifetime for the prompt process is fixed to 1 ns, whereas other parameters are allowed to be varied.

For processes that occur with a quantum yield that is temperature dependent in the temperature range used in PAC measurements, the thermodynamic parameters for the fast phase (<20ns) are determined by plotting [Ehν(1-)]/] versus (Cpρ/β) according to Eq. 7 and the volume and enthalpy changes for the subsequent steps are obtained by plotting (Ehν/) versus (Cpρ/β) according to Eq. 8 (Peters et al., 1992).

$$\frac{E\_{h\nu}(1-\phi)}{\phi} = \ -\Delta H + \Delta V \frac{c\_p \varrho}{\beta} \tag{7}$$

$$\frac{\phi\_{Eh\nu}}{\Phi} = -\Delta H + \Delta V (\frac{c\_p \rho}{\beta}) \tag{8}$$

where ΔH and ΔV correspond to the reaction enthalpy and volume change, respectively.

## **3. Results**

304 Thermodynamics – Interaction Studies – Solids, Liquids and Gases

The photo-acoustic set-up used in our lab was described previously (Miksovska&Larsen, 2003). Briefly, the sample in a quartz cell was placed in a temperature controlled holder (Quantum Northwest). The 532 nm output from a Nd:YAG laser (7 ns pulse width , < 50 µJ power) was shaped using a narrow slit (100 µm) and focused on the center of the quartz cell. An acoustic detector (1 MHz, RV103, Panametrix) was coupled to the side of a quartz cell using a thin layer of honey and the detector output was amplified using an ultrasonic preamplifier (Panametrics 5662). The signal was then stored in a digitizer (Wave Surfer 42Xs, 400 MHz). The data were analyzed using Sound Analysis software (Quantum

The excitation of the photocleavable iron-ligand bond in heme proteins generates at least two processes that contribute to the photoacoustic wave: the volume change due to the heat released during the reaction (Q), and the volume change (ΔV') due to the photo-triggered structural changes (including bond cleavage / formation, electrostriction, solvation, etc.).

ద

ఉ ద

where K is the instrument response constant, Ea is number of Einsteins absorbed, β is the expansion coefficient, ρ is the density, and Cp is the heat capacity. For water, the (β/Cpρ) term strongly varies with temperature mainly due to the temperature dependence of the β term. To evaluate the instrument response constant, the photo-acoustic traces are measured for a reference compound under experimental conditions (laser power, temperature, etc.) identical to those for the sample measurements. We have used Fe(III)4SP as a reference since it is non-fluorescent and photo-chemically stable. The amplitude of the reference acoustic

ఔܧܧܭ ൌ ܣ

ೞೌ 

where Ehν is the energy of a photon at 532 nm (Ehν= 53.7 kcal mol-1). The amount of heat deposited to the solvent and the non-thermal volume changes can then be determined by measuring the acoustic wave for the sample and the reference for several temperatures and plotting the ratio of the sample and reference acoustic wave () as a function of (Cpρ/β)

ൌ ߶ ൌ ܳ ȟܸԢ ఘ

For a multi-step process that exhibits volume and enthalpy changes on the time-scale between ~ 20 ns to 5 µs, the thermodynamic parameters for each individual step and corresponding lifetimes can be determined due to the sensitivity of the acoustic detector to the temporal profile of the pressure change. The time dependent sample acoustic signal E(t)obs can be expressed as a convolution of the time dependent function describing the volume change H(t) with the instrument response T(t) function (the reference acoustic

ȟܸԢሻ (2)

(3)

ఉ (4)

ܣ௦ ൌ ܧܭሺܳ <sup>ఉ</sup>

The amplitude of the sample acoustic wave (Asam) can be expressed as:

**2.2 Photoacoustic calorimetry** 

Northwest).

**2.3 Data analysis** 

trace can be described as:

according to Eq. 4:

wave):

Ligand migration in heme proteins is often described using the sequential three-state model (Henry et al., 1983) shown in Scheme 2. Upon cleavage of the coordination bond between the ligand and heme iron, the ligand is temporarily trapped within the protein matrix and then it either directly rebinds back to the heme iron in the so called "geminate rebinding" or diffuses from the protein matrix into the surrounding solvent. The subsequent bimolecular ligand binding to heme iron occurs on significantly longer time scales, hundreds of microseconds to milliseconds. The quantum yield for the geminate rebinding and for bimolecular association is strongly dependent on the character of the ligand and the protein. For example, CO rebinds to Mb predominantly through a bimolecular reaction with quantum yield close to unity (bim = 0.96 )(Henry et al., 1983), whereas the quantum yield for bimolecular O2 rebinding to heme proteins is significantly lower (Carver et al., 1990; Walda et al., 1994), and NO rebinds predominantly through geminate rebinding (Ye et al., 2002). To determine the thermodynamic parameters from acoustic data, the quantum yields for CO and O2 bimolecular rebinding to Hb and Mb, respectively, have to be known. The quantum yield for O2 binding to Mb was measured in the temperature range from 5 C to 35 C (Fig. 2) and the values show a weak temperature dependence with the quantum yield decreasing with increasing temperature. At 20 C the quantum yield is 0.09 ± 0.01 that is within the range of values reported previously ( = 0.057 (Walda et al., 1994) and ( = 0.12 (Carver et al., 1990)). We have also measured the quantum yield for CO bimolecular rebinding to Hb, and to Hb in the presence of effector molecules (Fig. 2). The quantum yield increases linearly with temperature and at 20 C, CO binds to Hb with quantum yield of

Time Resolved Thermodynamics Associated with Diatomic Ligand Dissociation from Globins 307

from the transient "deoxy-Mb" intermediate into the surrounding solvent and is approximately 3 times faster than the rate of the CO escape (Westrick et al., 1990), which suggests that O2 diffuses from the protein matrix through a transient channel with a lower activation barrier than CO. This result is consistent with the transient absorption studies that estimated the rate for O2 release to be approximately two times faster than that for CO (Carver et al., 1990). Interestingly, a similar time-constant of 200 ns to 300 ns was determined for CO escape from Mb at pH 3.5 (Angeloni&Feis, 2003). At acidic pH Mb adopts an open conformation with His 64 displaced toward the solvent giving a direct access to the distal cavity. These data suggest that the reorientation of His 64 may not be a

> CO-Hb + benzafibrate CO-Hb + phytic acid

CO-Hb

Fig. 2. Quantum yield for bimolecular photo-dissociation of O2 from the O2-Mb complex (bottom) and CO from the CO- Hb complex (top) as a function of temperature. The error of

5 10 15 20

temperature (oC)

0 10 20 30 40

**oC**)

**temperature** (

quantum yield is ± 0.05. The solid line demonstrates the trend.

rate limiting step for the O2 escape.

0.3

0.09

0.10

0.11

0.12

0.4

0.5

0.6

0.7

0.8

0.9

0.68 and in the presence of IHP and BZF 0.62 and 0.46, respectively. A similar quantum yield for CO bimolecular rebinding to Hb was reported previously by Unno et al. (bim =0.7 at 20 C) (Unno et al., 1990) and by Saffran and Gibson (=0.7 for CO binding to Hb and = 0.73 for CO association to Hb in the presence of IHP at 40 C) (Saffran&Gibson, 1977).

#### Scheme 2.

The photo-acoustic traces for O2 dissociation from Fe(II)Mb at pH 7.0 are shown in Fig. 3. At low temperatures (6 C to 15 C), the sample photoacoustic traces show a phase shift with respect to the reference trace indicating the presence of thermodynamic process(es) that occurs between 50 ns and ~ 5 µs. The sample traces were deconvoluted as described in the Materials and Methods section and the i values were plotted as a function of the temperature dependent factor (Cpρ/β) (Fig. 4). The extrapolated volume and enthalpy changes are listed in Table 1. The photo-cleavage of the Fe-O2 bond is associated with a fast structural relaxation (< 20 ns) forming a transient "deoxy-Mb intermediate". This initial transition is endothermic (ΔH = 21 ± 9 kcal mol-1) and leads to a small volume contraction of – 3.0 ± 0.5 mL mol-1. This initial relaxation is followed by ~ 250 ns kinetics that exhibit a volume increase of 5.5 ± 0.4 mL mol-1 and a very small enthalpy change of -8.9 ± 8.0 kcal mol-1. We associate the initial process with the photo-cleavage of the Fe-O2 bond. A similar volume decrease of approximately -3 mL mol-1 has been observed previously for the photodissociation of Fe-CO bond in Mb (Westrick&Peters, 1990; Westrick et al., 1990). The observed volume contraction reflects a fast relaxation of the heme binding pocket including: i) cleavage of the hydrogen bond between the distal histidine and oxygen molecule (Phillips&Schoenborn, 1981) ii) reorientation of distal residues within the heme binding pocket (Olson et al., 2007), and iii) fast migration of the photo-released ligand into the primary docking site and then into the internal cavities (Xe4 or Xe1) (Hummer et al., 2004). Also, the positive enthalpy change is consistent with the photo-cleavage of Fe-O2 bond. The subsequent 250 ns kinetics may reflect either the nanosecond geminate rebinding of the O2 molecule or the ligand diffusion from the protein matrix into the surrounding solvent. The kinetics for the geminate O2 rebinding were studied on femtosecond timescale by

Petrich et al. (Petrich et al., 1988), and on picosecond and nanosecond timescales (Carver et al., 1990; Miller et al., 1996). These studies identified two distinct sub-states of the "deoxyMb" intermediate: a "barrier-less" and a "photolyzable" sub-state. In the "barrierless" sub-state, oxygen rebinds to heme iron on sub-picosecond timescale whereas the oxygen association to the "photolyzable" substate occurs on nanosecond and microsecond timescales. Carver et al. (Carver et al., 1990) have reported the time constant for O2 nanosecond geminate rebinding to be 52 ± 14 ns at room temperature. This kinetic step has a lifetime that is comparable to the time resolution of our PAC instrument ( ~ 50 ns) and therefore it was not resolved in this study. The 250 ns step thus corresponds to the O2 escape

0.68 and in the presence of IHP and BZF 0.62 and 0.46, respectively. A similar quantum yield for CO bimolecular rebinding to Hb was reported previously by Unno et al. (bim =0.7 at 20 C) (Unno et al., 1990) and by Saffran and Gibson (=0.7 for CO binding to Hb and = 0.73 for CO association to Hb in the presence of IHP at 40 C) (Saffran&Gibson, 1977).

The photo-acoustic traces for O2 dissociation from Fe(II)Mb at pH 7.0 are shown in Fig. 3. At low temperatures (6 C to 15 C), the sample photoacoustic traces show a phase shift with respect to the reference trace indicating the presence of thermodynamic process(es) that occurs between 50 ns and ~ 5 µs. The sample traces were deconvoluted as described in the Materials and Methods section and the i values were plotted as a function of the temperature dependent factor (Cpρ/β) (Fig. 4). The extrapolated volume and enthalpy changes are listed in Table 1. The photo-cleavage of the Fe-O2 bond is associated with a fast structural relaxation (< 20 ns) forming a transient "deoxy-Mb intermediate". This initial transition is endothermic (ΔH = 21 ± 9 kcal mol-1) and leads to a small volume contraction of – 3.0 ± 0.5 mL mol-1. This initial relaxation is followed by ~ 250 ns kinetics that exhibit a volume increase of 5.5 ± 0.4 mL mol-1 and a very small enthalpy change of -8.9 ± 8.0 kcal mol-1. We associate the initial process with the photo-cleavage of the Fe-O2 bond. A similar volume decrease of approximately -3 mL mol-1 has been observed previously for the photodissociation of Fe-CO bond in Mb (Westrick&Peters, 1990; Westrick et al., 1990). The observed volume contraction reflects a fast relaxation of the heme binding pocket including: i) cleavage of the hydrogen bond between the distal histidine and oxygen molecule (Phillips&Schoenborn, 1981) ii) reorientation of distal residues within the heme binding pocket (Olson et al., 2007), and iii) fast migration of the photo-released ligand into the primary docking site and then into the internal cavities (Xe4 or Xe1) (Hummer et al., 2004). Also, the positive enthalpy change is consistent with the photo-cleavage of Fe-O2 bond. The subsequent 250 ns kinetics may reflect either the nanosecond geminate rebinding of the O2 molecule or the ligand diffusion from the protein matrix into the surrounding solvent. The kinetics for the geminate O2 rebinding were studied on femtosecond timescale by Petrich et al. (Petrich et al., 1988), and on picosecond and nanosecond timescales (Carver et al., 1990; Miller et al., 1996). These studies identified two distinct sub-states of the "deoxyMb" intermediate: a "barrier-less" and a "photolyzable" sub-state. In the "barrierless" sub-state, oxygen rebinds to heme iron on sub-picosecond timescale whereas the oxygen association to the "photolyzable" substate occurs on nanosecond and microsecond timescales. Carver et al. (Carver et al., 1990) have reported the time constant for O2 nanosecond geminate rebinding to be 52 ± 14 ns at room temperature. This kinetic step has a lifetime that is comparable to the time resolution of our PAC instrument ( ~ 50 ns) and therefore it was not resolved in this study. The 250 ns step thus corresponds to the O2 escape

Scheme 2.

from the transient "deoxy-Mb" intermediate into the surrounding solvent and is approximately 3 times faster than the rate of the CO escape (Westrick et al., 1990), which suggests that O2 diffuses from the protein matrix through a transient channel with a lower activation barrier than CO. This result is consistent with the transient absorption studies that estimated the rate for O2 release to be approximately two times faster than that for CO (Carver et al., 1990). Interestingly, a similar time-constant of 200 ns to 300 ns was determined for CO escape from Mb at pH 3.5 (Angeloni&Feis, 2003). At acidic pH Mb adopts an open conformation with His 64 displaced toward the solvent giving a direct access to the distal cavity. These data suggest that the reorientation of His 64 may not be a rate limiting step for the O2 escape.

Fig. 2. Quantum yield for bimolecular photo-dissociation of O2 from the O2-Mb complex (bottom) and CO from the CO- Hb complex (top) as a function of temperature. The error of quantum yield is ± 0.05. The solid line demonstrates the trend.

Time Resolved Thermodynamics Associated with Diatomic Ligand Dissociation from Globins 309

reaction volume can be expressed as the difference between the partial molar volume of products and reactants according to: ΔVslow= VO2 + VdeoxyMb- VO2-Mb - VH2O , where VO2 is the partial molar volume of oxygen, VH2O is the partial molar volume of water, VdeoxyMb is the partial molar volume of transient "deoxyMb" intermediate and VO2-Mb is the partial molar volume of oxy-Mb. Using VO2 = 28 mL mol-1 (Projahn et al., 1990) and VH2O = 15 mL mol-1 (the partial molar volume of water scaled to the occupancy of water molecule hydrogen bound to distal histidine) (Belogortseva et al., 2007), we estimate that the O2 release from Mb results in a structural volume change (VdoxyMb- VO2-Mb) of - 7.5 mL mol-1. This value is very similar to that reported previously for CO escape from Mb (ΔVstructural= VdoxyMb- VCO-Mb = - 6 mL mol-1) (Vetromile, et al., 2011) demonstrating that the overall structural changes accompanying the ligand bound to ligand free transition in Mb are very similar for both ligands. This is in agreement with the close resemblance of the X-ray structure of both the CO-bound and O2-bound Mb (Yang&Phillips Jr, 1996). The small enthalpy change measured for the 250 ns relaxation (ΔH = -8.9 ± 8.0 kcal mol-1) includes the enthalpy change for O2 solvation (ΔHsolv = -2.9 kcal mol-1 (Mills et al., 1979)) and the enthalpy change associated with H2O binding to the heme binding pocket (ΔHsolv = -7 kcal mol-1 (Vetromile, et al., 2011) indicating that the structural relaxation coupled to the ligand

The overall enthalpy change for O2 dissociation from Mb was determined to be 11.6 ± 8.5 kcal mol-1 and this value is in agreement with the value of 10 kcal mol-1 reported previously (Projahn et al., 1990). The overall reaction volume change determined here (ΔVoverall = +2.5 mL mol-1) is somewhat larger than the reaction volume change determined from the measurement of the equilibrium constant as a function of pressure (ΔV= - 2.9 mL mol-1) (Hasinoff, 1974) and significantly smaller than the reaction volume change determined as a difference between the activation volume for oxygen binding and dissociation from Mb that was reported to be 18 mL mol-1 (Projahn et al., 1990). Unlike photoacoustic studies that allow for reaction volume determination at ambient pressure, the high pressure measurements of equilibrium constant and/or rate constants (to determine activation volumes) may cause a pressure induced protein denaturation and/or structural changes, which may influence the magnitude of reaction volume changes in high pressure studies.

Fast phase -3.0 ± 0.5 20.5 ± 8.5

Slow phase 5.5 ± 0.4 -8.9 ± 8.0

Table 1. Volume and enthalpy changes associated with O2 dissociation from Mb in the

We have also probed the thermodynamic parameters associated with the CO photodissociation from Hb and the impact of the binding of BZF and IHP on the thermodynamics associated with the ligand migration between the heme binding pocket and surrounding solvent. The photo-acoustic traces for CO photo-dissociation from Hb are shown in Fig. 5. The sample and the reference acoustic wave overlay in phase indicating that the observed thermodynamic processes take place within 50 ns upon photo-dissociation, which is consistent with the fast CO diffusion from the heme matrix into the surrounding solvent.

**ΔV (mL mol-1) ΔH (kcal mol-1)** 

escape from the protein is entropy driven.

temperature range 6 - 10C.

Fig. 3. PAC traces for O2 photo-dissociation from O2-Mb complex at 6 C. Conditions: 40 µM Mb dissolved in 50 mM Hepes buffer pH 7.0. The absorbance of the reference compound, Fe(III)4SP, at excitation wavelength of 532 nm was identical as that of O2-Mb.

Fig. 4. Plot of the ratio of the acoustic amplitude for the photo-dissociation of the O2-Mb complex and the reference compound as a function of (Cpρ/β) term. 1 values that correspond to the prompt phase are shown as solid circles and the 2 values corresponding to the slow phase are shown as open squares. The data were fit with a linear curve and the corresponding volume and enthalpy changes were determined using Eq. 6 and Eq. 7.

The reaction volume change observed for the slow phase includes several factors: i) volume change due to the O2 escape into the surrounding solvent, ii) volume change associated with the heme hydration in deoxyMb and iii) volume change due to the structural changes. The

Fig. 3. PAC traces for O2 photo-dissociation from O2-Mb complex at 6 C. Conditions: 40 µM Mb dissolved in 50 mM Hepes buffer pH 7.0. The absorbance of the reference compound,

time (s)

3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4

Fig. 4. Plot of the ratio of the acoustic amplitude for the photo-dissociation of the O2-Mb complex and the reference compound as a function of (Cpρ/β) term. 1 values that

Cp

correspond to the prompt phase are shown as solid circles and the 2 values corresponding to the slow phase are shown as open squares. The data were fit with a linear curve and the corresponding volume and enthalpy changes were determined using Eq. 6 and Eq. 7.

10 15 20 25 30 35

(kcal mL-1

)

50

100

150

200

2 Eh /

(kcal mol-1

)

250

The reaction volume change observed for the slow phase includes several factors: i) volume change due to the O2 escape into the surrounding solvent, ii) volume change associated with the heme hydration in deoxyMb and iii) volume change due to the structural changes. The

Fe(III)4SP, at excitation wavelength of 532 nm was identical as that of O2-Mb.

 reference sample




Eh(1-1)/(kcal mol-1

)



0.00

0.02

PAC signal (a.u.)

0.04

0.06

0.08

0.10

reaction volume can be expressed as the difference between the partial molar volume of products and reactants according to: ΔVslow= VO2 + VdeoxyMb- VO2-Mb - VH2O , where VO2 is the partial molar volume of oxygen, VH2O is the partial molar volume of water, VdeoxyMb is the partial molar volume of transient "deoxyMb" intermediate and VO2-Mb is the partial molar volume of oxy-Mb. Using VO2 = 28 mL mol-1 (Projahn et al., 1990) and VH2O = 15 mL mol-1 (the partial molar volume of water scaled to the occupancy of water molecule hydrogen bound to distal histidine) (Belogortseva et al., 2007), we estimate that the O2 release from Mb results in a structural volume change (VdoxyMb- VO2-Mb) of - 7.5 mL mol-1. This value is very similar to that reported previously for CO escape from Mb (ΔVstructural= VdoxyMb- VCO-Mb = - 6 mL mol-1) (Vetromile, et al., 2011) demonstrating that the overall structural changes accompanying the ligand bound to ligand free transition in Mb are very similar for both ligands. This is in agreement with the close resemblance of the X-ray structure of both the CO-bound and O2-bound Mb (Yang&Phillips Jr, 1996). The small enthalpy change measured for the 250 ns relaxation (ΔH = -8.9 ± 8.0 kcal mol-1) includes the enthalpy change for O2 solvation (ΔHsolv = -2.9 kcal mol-1 (Mills et al., 1979)) and the enthalpy change associated with H2O binding to the heme binding pocket (ΔHsolv = -7 kcal mol-1 (Vetromile, et al., 2011) indicating that the structural relaxation coupled to the ligand escape from the protein is entropy driven.

The overall enthalpy change for O2 dissociation from Mb was determined to be 11.6 ± 8.5 kcal mol-1 and this value is in agreement with the value of 10 kcal mol-1 reported previously (Projahn et al., 1990). The overall reaction volume change determined here (ΔVoverall = +2.5 mL mol-1) is somewhat larger than the reaction volume change determined from the measurement of the equilibrium constant as a function of pressure (ΔV= - 2.9 mL mol-1) (Hasinoff, 1974) and significantly smaller than the reaction volume change determined as a difference between the activation volume for oxygen binding and dissociation from Mb that was reported to be 18 mL mol-1 (Projahn et al., 1990). Unlike photoacoustic studies that allow for reaction volume determination at ambient pressure, the high pressure measurements of equilibrium constant and/or rate constants (to determine activation volumes) may cause a pressure induced protein denaturation and/or structural changes, which may influence the magnitude of reaction volume changes in high pressure studies.


Table 1. Volume and enthalpy changes associated with O2 dissociation from Mb in the temperature range 6 - 10C.

We have also probed the thermodynamic parameters associated with the CO photodissociation from Hb and the impact of the binding of BZF and IHP on the thermodynamics associated with the ligand migration between the heme binding pocket and surrounding solvent. The photo-acoustic traces for CO photo-dissociation from Hb are shown in Fig. 5. The sample and the reference acoustic wave overlay in phase indicating that the observed thermodynamic processes take place within 50 ns upon photo-dissociation, which is consistent with the fast CO diffusion from the heme matrix into the surrounding solvent.

Time Resolved Thermodynamics Associated with Diatomic Ligand Dissociation from Globins 311

and triple ligated (CO)3Hb. Consequently, the observed reaction enthalpy corresponds to the enthalpy change due to the cleavage of the Fe-CO bond (ΔHFe-CO=17.5 kcal mol-1 (Leung et al., 1987; Miksovska et al., 2005)), the enthalpy change due to the solvation of a CO molecule (ΔHsolv = 2.6 kcal mol-1 (Leung et al., 1987)), the enthalpy change of structural relaxation associated with the ligand release from the protein matrix, and enthalpy of the distal pocket hydration. The occupancy of water molecules in the distal pocket of deoxyHb was determined to be significantly lower than that in Mb (~0.64 for the Hb - chain and ~ 0.33 for the Hb β-chain (Esquerra et al., 2010)). Using an average occupancy of 0.48, we estimate that the distal pocket hydration contributes to the overall enthalpy change by ~ - 3 kcal mol-1 (Vetromile, et al., 2011). Therefore, the structural relaxation coupled to the CO dissociation and diffusion into the surrounding solvent is

Fig. 6. The plot of the ratio of the acoustic amplitude for the CO photo-dissociation from the CO-Hb complex and the reference compound as a function of the temperature dependent factor (Cpρ/β) term. The reaction volume and enthalpy changes were extrapolated

CO-Hb +benzafibrate

CO-Hb+IHP

3.0 3.5 4.0 4.5 5.0 5.5 6.0

Cp (kcal mL-1)

Analogous to O2 photo-release from Mb, the observed reaction volume change for CO photorelease from Hb , ΔV=21.5 mL mol-1, can be expressed as: ΔV = VCO + V(CO)3Hb - V(CO)4 Hb - VH2O, where VCO is the partial molar volume of CO and V(CO)3Hb and V(CO)4 Hb are the partial molar volume of (CO)3Hb and (CO)4Hb, respectively. Using VCO = 35 mL mol-1 (Projahn et al., 1990) and VH2O = 9 mL mol-1 (partial molar volume of water scaled by the average occupancy of the Hb chain), we estimate that upon release of one CO molecule per Hb, the protein undergoes a small contraction of -7 mL mol-1. The small volume change observed here is consistent with the minor structural changes due to deligation of Hb in the R-state as observed in the X-ray structure that are predominantly

accompanied by a small enthalpy change of 2 kcal mol-1.

100 CO-Hb

according to Eq. 5

40

50

[Eh

() ]/kcal mol-1

60

70

80

90

The fast ligand escape from the heme binding pocket was observed in the presence of effectors (data not shown). Previous transient absorption studies showed that the CO photo-release from the fully ligated R-state Hb is followed by three relaxations with lifetimes of 50 ns, 1 µs, and 20 µs that were assigned to the unimolecular geminate rebinding, the tertiary structural relaxation, and the RT quaternary change, respectively (Goldbeck et al., 1996). The geminate rebinding occurs too fast to be resolved by our PAC detector, whereas the 20 µs RT transition, which strongly depends on the extent of heme ligation, is too slow to be resolved in PAC measurements. The 1 µs relaxation is within the time-window accessible by our detection system, however we were unable to resolve this step. Since this relaxation was observed as a small perturbation of the deoxy-Soret band (Goldbeck et al., 1996), it may reflect the structural relaxation localized within the vicinity of the heme binding pocket, which does not lead to measurable volume and enthalpy changes.

The volume and enthalpy changes associated with the diffusion of the photo-dissociated ligand to the surrounding solvent can be determined from the plot of the ratio of the amplitude of the acoustic trace for CO photo-dissociation from Hb and the reference as a function of temperature according to Eq. 7 (Fig. 6). The extrapolated thermodynamic values

Fig. 5. PAC traces for the CO photo-dissociation from the CO-Hb complex and the reference compound Fe(III)4SP. Conditions: 40 µM Hb in 100 mM HEPES buffer pH 7.0 and 20 C. The absorbance of the reference compound matched the absorbance of the sample at 532 nm.

are shown in Table 2. The CO photo-release from Hb is associated with a positive volume change of 21.5 ± 0.9 mL mol-1 and enthalpy change of 19.4 ± 1.2 kcal mol-1. These results are in agreement with the thermodynamic parameters reported previously by Peters et al: ΔV = 23.4 ± 0.5 mL mol-1 and ΔH = 18.0 ± 2.9 kcal mol-1 (Peters et al., 1992). Since the laser power used in this study was kept below 50 µJ, the low level of photo-dissociation was achieved that corresponds to 1 CO molecule per hemoglobin photo-released. Thus the observed thermodynamic parameters reflect the transition between fully ligated (CO)4Hb

The fast ligand escape from the heme binding pocket was observed in the presence of effectors (data not shown). Previous transient absorption studies showed that the CO photo-release from the fully ligated R-state Hb is followed by three relaxations with lifetimes of 50 ns, 1 µs, and 20 µs that were assigned to the unimolecular geminate rebinding, the tertiary structural relaxation, and the RT quaternary change, respectively (Goldbeck et al., 1996). The geminate rebinding occurs too fast to be resolved by our PAC detector, whereas the 20 µs RT transition, which strongly depends on the extent of heme ligation, is too slow to be resolved in PAC measurements. The 1 µs relaxation is within the time-window accessible by our detection system, however we were unable to resolve this step. Since this relaxation was observed as a small perturbation of the deoxy-Soret band (Goldbeck et al., 1996), it may reflect the structural relaxation localized within the vicinity of the heme binding pocket, which does not lead to measurable volume and enthalpy

The volume and enthalpy changes associated with the diffusion of the photo-dissociated ligand to the surrounding solvent can be determined from the plot of the ratio of the amplitude of the acoustic trace for CO photo-dissociation from Hb and the reference as a function of temperature according to Eq. 7 (Fig. 6). The extrapolated thermodynamic values

> 4SP CO-Hb

Fig. 5. PAC traces for the CO photo-dissociation from the CO-Hb complex and the reference compound Fe(III)4SP. Conditions: 40 µM Hb in 100 mM HEPES buffer pH 7.0 and 20 C. The absorbance of the reference compound matched the absorbance of the sample at 532

3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 4.6

time (s)

are shown in Table 2. The CO photo-release from Hb is associated with a positive volume change of 21.5 ± 0.9 mL mol-1 and enthalpy change of 19.4 ± 1.2 kcal mol-1. These results are in agreement with the thermodynamic parameters reported previously by Peters et al: ΔV = 23.4 ± 0.5 mL mol-1 and ΔH = 18.0 ± 2.9 kcal mol-1 (Peters et al., 1992). Since the laser power used in this study was kept below 50 µJ, the low level of photo-dissociation was achieved that corresponds to 1 CO molecule per hemoglobin photo-released. Thus the observed thermodynamic parameters reflect the transition between fully ligated (CO)4Hb

changes.

nm.


PAC amplitude (a.u.)

and triple ligated (CO)3Hb. Consequently, the observed reaction enthalpy corresponds to the enthalpy change due to the cleavage of the Fe-CO bond (ΔHFe-CO=17.5 kcal mol-1 (Leung et al., 1987; Miksovska et al., 2005)), the enthalpy change due to the solvation of a CO molecule (ΔHsolv = 2.6 kcal mol-1 (Leung et al., 1987)), the enthalpy change of structural relaxation associated with the ligand release from the protein matrix, and enthalpy of the distal pocket hydration. The occupancy of water molecules in the distal pocket of deoxyHb was determined to be significantly lower than that in Mb (~0.64 for the Hb - chain and ~ 0.33 for the Hb β-chain (Esquerra et al., 2010)). Using an average occupancy of 0.48, we estimate that the distal pocket hydration contributes to the overall enthalpy change by ~ - 3 kcal mol-1 (Vetromile, et al., 2011). Therefore, the structural relaxation coupled to the CO dissociation and diffusion into the surrounding solvent is accompanied by a small enthalpy change of 2 kcal mol-1.

Fig. 6. The plot of the ratio of the acoustic amplitude for the CO photo-dissociation from the CO-Hb complex and the reference compound as a function of the temperature dependent factor (Cpρ/β) term. The reaction volume and enthalpy changes were extrapolated according to Eq. 5

Analogous to O2 photo-release from Mb, the observed reaction volume change for CO photorelease from Hb , ΔV=21.5 mL mol-1, can be expressed as: ΔV = VCO + V(CO)3Hb - V(CO)4 Hb - VH2O, where VCO is the partial molar volume of CO and V(CO)3Hb and V(CO)4 Hb are the partial molar volume of (CO)3Hb and (CO)4Hb, respectively. Using VCO = 35 mL mol-1 (Projahn et al., 1990) and VH2O = 9 mL mol-1 (partial molar volume of water scaled by the average occupancy of the Hb chain), we estimate that upon release of one CO molecule per Hb, the protein undergoes a small contraction of -7 mL mol-1. The small volume change observed here is consistent with the minor structural changes due to deligation of Hb in the R-state as observed in the X-ray structure that are predominantly

Time Resolved Thermodynamics Associated with Diatomic Ligand Dissociation from Globins 313

We have also determined volume and enthalpy changes associated with the CO photodissociation from Hb in the presence of heterogenous effectors BZT and IHP (Fig. 6) and the thermodynamic profiles for CO photo-dissociation from CO-Hb complex in the presence and absence of effectors are presented in Fig.7 . Both effectors bind to Hb in the T-state and R-state and modulate the interaction of Hb with diatomic ligands (Coletta et al., 1999b; Marden et al., 1990). For example, the binding of BZF or IHP to CO-Hb complex decreases the CO association rate approximately four or eight times, respectively (Marden et al., 1990), and decreases the affinity of R state deoxy-Hb for oxygen (Tsuneshige et al., 2002). Coletta et al (Coletta et al., 1999a) have reported that simultaneous binding or IHP and BZF effectors to Hb at ambient pressure leads to the Hb intermediate with tertiary T-like structure in the quaternary R- conformation. Recently, using NMR spectroscopy Song et al. have shown that binding of IHP to the fully ligated Hb increase the conformational

The photoacoustic data presented here show that BZF binding to CO-Hb complex does not impact the reaction volume and enthalpy changes associated with CO photo-release. The crystal structure of horse CO-Hb in complex with BZF indicates that the structural changes due to BZF association to fully ligated Hb are localized in the -subunits (Shibayama et al., 2002). BZF binds to the surface of each -chain E-helix and decreases the distance between the heme iron and distal His and its binding site is surrounded by hydrophobic residues such as Ala 65, Leu 68, Leu 80 and Leu 83 (Shibayama et al., 2002). Such minor structural changes caused by BZF association are unlikely to alter the overall structural volume and enthalpy changes associated with the CO photo-release. However, due to the lower solubility of BZF, the effector concentration used is this study was 5 mM that results in a Hb fractional saturation of 0.25 (using KD of 15 mM (Ascenzi et al., 1993)). Such lower fractional saturation may prevent detection of BZF induced changes in Hb

On the other hand, the binding of IHP has a significant impact on the observed volume and enthalpy changes (Table 2). The reaction volume decreases by 10 mL mol-1 and the enthalpy change is more exothermic by nearly 30 kcal mol-1 compared to the thermodynamic parameters determined in the absence of effector molecules. Such negative reaction volume and exothermic enthalpy change indicates that electrostriction of solvent molecules caused by reorganization of salt bridges or redistribution of charges on protein surface contributes to the overall reaction volume and enthalpy change associated with the CO photo-release. Indeed, IHP interacts with charged residues along the Hb central cavity. At the Hb T-state, the IHP binding site is located at the interface of the β-chains involving Val 1, His2, Lys 82 and His 141 from each chain (Riccio et al., 2001); whereas at the R-state Hb, the IHP molecule interacts with the charged residues Lys 99 and Arg 141 from each -chain (Laberge et al., 2005). In the absence of the X-ray structure of IHP bound fully ligated and partially photolyzed CO-Hb, it is difficult to point out the factors that contribute to the observed volume and enthalpy changes on the molecular level. Arg 141 forms a salt bridge with Asp 126 in the T-state deoxy Hb that is absent in the fully ligated R- state (Park et al., 2006). We speculate that the transition between the fully ligated (CO)4Hb and partially ligated (CO)3Hb may be associated with a repositioning of the Arg 141 side chain leading to a partial exposure of either the IHP molecule and/or the Arg 141 side chain to the surrounding solvent molecules. Also, the

fluctuation of the R-state in both the - and β-chain (Song et al., 2008).

conformational dynamics.

localized in the the -chain and include reposition of the F-helix and shift of the EF and CD corner (Wilson et al., 1996).


Table 2. Volume and enthalpy changes associated with CO photo-dissociation from Hb.

Fig. 7. The thermodynamic profile for CO photo-dissociation from Hb in the absence of effector and in the presence of BZF and IHP.

localized in the the -chain and include reposition of the F-helix and shift of the EF and

CO-Hb 19.4 ± 1.2 21. 5 ± 0.9

CO-Hb + BZF 21.7 ± 7.9 22.3 ± 1.7

CO-Hb + IHP -9.9 ± 6.1 11.4 ± 1.3

Table 2. Volume and enthalpy changes associated with CO photo-dissociation from Hb.

Fig. 7. The thermodynamic profile for CO photo-dissociation from Hb in the absence of

effector and in the presence of BZF and IHP.

**ΔHprompt (kcal mol-1) ΔVprompt (mL mol-1)** 

CD corner (Wilson et al., 1996).

We have also determined volume and enthalpy changes associated with the CO photodissociation from Hb in the presence of heterogenous effectors BZT and IHP (Fig. 6) and the thermodynamic profiles for CO photo-dissociation from CO-Hb complex in the presence and absence of effectors are presented in Fig.7 . Both effectors bind to Hb in the T-state and R-state and modulate the interaction of Hb with diatomic ligands (Coletta et al., 1999b; Marden et al., 1990). For example, the binding of BZF or IHP to CO-Hb complex decreases the CO association rate approximately four or eight times, respectively (Marden et al., 1990), and decreases the affinity of R state deoxy-Hb for oxygen (Tsuneshige et al., 2002). Coletta et al (Coletta et al., 1999a) have reported that simultaneous binding or IHP and BZF effectors to Hb at ambient pressure leads to the Hb intermediate with tertiary T-like structure in the quaternary R- conformation. Recently, using NMR spectroscopy Song et al. have shown that binding of IHP to the fully ligated Hb increase the conformational fluctuation of the R-state in both the - and β-chain (Song et al., 2008).

The photoacoustic data presented here show that BZF binding to CO-Hb complex does not impact the reaction volume and enthalpy changes associated with CO photo-release. The crystal structure of horse CO-Hb in complex with BZF indicates that the structural changes due to BZF association to fully ligated Hb are localized in the -subunits (Shibayama et al., 2002). BZF binds to the surface of each -chain E-helix and decreases the distance between the heme iron and distal His and its binding site is surrounded by hydrophobic residues such as Ala 65, Leu 68, Leu 80 and Leu 83 (Shibayama et al., 2002). Such minor structural changes caused by BZF association are unlikely to alter the overall structural volume and enthalpy changes associated with the CO photo-release. However, due to the lower solubility of BZF, the effector concentration used is this study was 5 mM that results in a Hb fractional saturation of 0.25 (using KD of 15 mM (Ascenzi et al., 1993)). Such lower fractional saturation may prevent detection of BZF induced changes in Hb conformational dynamics.

On the other hand, the binding of IHP has a significant impact on the observed volume and enthalpy changes (Table 2). The reaction volume decreases by 10 mL mol-1 and the enthalpy change is more exothermic by nearly 30 kcal mol-1 compared to the thermodynamic parameters determined in the absence of effector molecules. Such negative reaction volume and exothermic enthalpy change indicates that electrostriction of solvent molecules caused by reorganization of salt bridges or redistribution of charges on protein surface contributes to the overall reaction volume and enthalpy change associated with the CO photo-release. Indeed, IHP interacts with charged residues along the Hb central cavity. At the Hb T-state, the IHP binding site is located at the interface of the β-chains involving Val 1, His2, Lys 82 and His 141 from each chain (Riccio et al., 2001); whereas at the R-state Hb, the IHP molecule interacts with the charged residues Lys 99 and Arg 141 from each -chain (Laberge et al., 2005). In the absence of the X-ray structure of IHP bound fully ligated and partially photolyzed CO-Hb, it is difficult to point out the factors that contribute to the observed volume and enthalpy changes on the molecular level. Arg 141 forms a salt bridge with Asp 126 in the T-state deoxy Hb that is absent in the fully ligated R- state (Park et al., 2006). We speculate that the transition between the fully ligated (CO)4Hb and partially ligated (CO)3Hb may be associated with a repositioning of the Arg 141 side chain leading to a partial exposure of either the IHP molecule and/or the Arg 141 side chain to the surrounding solvent molecules. Also, the

Time Resolved Thermodynamics Associated with Diatomic Ligand Dissociation from Globins 315

Carver, T. E., Rohlfs, R. J., Olson, J. S., Gibson, Q. H., Blackmore, R. S., Springer, B.

Coletta, M., Angeletti, M., Ascenzi, P., Bertollini, A., Della Longa, S., De Sanctis, G., Priori,

Coletta, M., Angeletti, M., Ascone, I., Boumis, G., Castellano, A. C., Dell Ariccia, M., Della

Esquerra, R. M., Lopez-Pena, I., Tipgunlakant, P., Birukou, I., Nguyen, R. L., Soman, J.,

Gensch, T.&Viappiani, C. (2003). Time-resolved photothermal methods: accessing time-

Gibson, Q. H., Olson, J. S., McKinnie, R. E.&Rohlfs, R. J. (1986). A kinetic description of ligand binding to sperm whale myoglobin. *J. Biol. Chem.,* 261, 22, pp. 10228-10239, Goldbeck, R. A., Paquette, S. J., Björling, S. C.&Kliger, D. S. (1996). Allosteric Intermediates

Hasinoff, B. B. (1974). Kinetic activation volumes of the binding of oxygen and carbon

Henry, E. R., Sommer, J. H., Hofrichter, J., Eaton, W. A.&Gellert, M. (1983). Geminate

Henzler-Wildman, K. A., Lei, M., Thai, V., Kerns, S. J., Karplus, M.&Kern, D. (2007). A

Hummer, G., Schotte, F.&Anfinrud, P. A. (2004). Unveiling functional protein motions with

Laberge, M., Kövesi, I., Yonetani, T.&Fidy, J. (2005). R-state hemoglobin bound to

Leung, W. P., Cho, K. C., Chau, S. K.&Choy, C. L. (1987). Measurement of the protein-ligand

photolysis apparatus. *Biochemistry,* 13, 15, pp. 3111-3117, 0006-2960

Unligated Hemoglobin. *Biophys. J.,* 76, 3, pp. 1532-1536, 0006-3495

*Photoch. Photobiol. Sci.,* 2, 7, pp. 699-721, 1474-905X

*J. Biol. Chem.,* 265, 32, pp. 20007-20020,

*Biol. Chem.,* 274, 11, pp. 6865-6874,

pp. 10270-10278, 1463-9076

8628-8639, 0006-2960

450, 7171, pp. 913-916, 1476-4687

579, 3, pp. 627-632, 0014-5793

*Acad. Sci. U.S.A.,* 101, 43, pp. 15330-15334,

*Phys. Lett.,* 141, 3, pp. 220-224, 0009-2614

0022-2836

A.&Sligar, S. G. (1990). Analysis of the kinetic barriers for ligand binding to sperm whale myoglobin using site-directed mutagenesis and laser photolysis techniques.

A. M., Santucci, R.&Amiconi, G. (1999a). Coupling of the Oxygen-linked Interaction Energy for Inositol Hexakisphosphate and Bezafibrate Binding to Human HbA0. *J.* 

Longa, S., De Sanctis, G., Priori, A. M., Santucci, R., Feis, A.&Amiconi, G. (1999b). Heterotropic Effectors Exert More Significant Strain on Monoligated than on

Olson, J. S., Kliger, D. S.&Goldbeck, R. A. (2010). Kinetic spectroscopy of heme hydration and ligand binding in myoglobin and isolated hemoglobin chains: an optical window into heme pocket water dynamics. *Phys. Chem. Chem. Phys.,* 12, 35,

resolved thermodynamics of photoinduced processes in chemistry and biology.

in Hemoglobin. 2. Kinetic Modeling of HbCO Photolysis. *Biochemistry,* 35, 26, pp.

monoxide to hemoglobin and myoglobin studied on a high-pressure laser flash

recombination of carbon monoxide to myoglobin. *J. Mol. Biol.,* 166, 3, pp. 443-451,

hierarchy of timescales in protein dynamics is linked to enzyme catalysis. *Nature,*

picosecond x-ray crystallography and molecular dynamics simulations.*Proc. Natl.* 

heterotropic effectors: models of the DPG, IHP and RSR13 binding sites. *FEBS Lett.,*

bond energy of carboxymyoglobin by pulsed photoacoustic calorimetry. *Chem.* 

ligand photo-release may be associated with the repositioning of the IHP molecule within the Hb central cavity. Based on a molecular dynamics simulation of IHP binding sides in south polar skua deoxyHb, an IHP migration pathway connecting the binding site at the interface between the -chains and the second binding site located between the β-chains was proposed suggesting that IHP interactions with Hb are dynamic and involve numerous positively charged residues situated along the central cavity (Riccio et al., 2001). Therefore, CO photo-release may trigger relocation of IHP within the central cavity resulting in larger exposure of IHP phosphate groups and/or charged amino acid residues and concomitant electrostriction of solvent molecules.
