Multiple geminate ligand recombinations in human hemoglobin

Dynamics based clustering of globin family members

A methodology to cluster proteins based on their dynamics’ similarity is presented. For each pair of proteins from a dataset, the structures are superimposed, and the Anisotropic Network Model modes of motions are calculated. The twelve slowest modes from each protein are matched using a local mode alignment algorithm based on the local sequence alignment algorithm of Smith–Waterman. The dynamical similarity distance matrix is calculated based on the top scoring matches of each pair and the proteins are clustered using a hierarchical clustering algorithm. The utility of this method is exemplified on a dataset of protein chains from the globin family and a dataset of tetrameric hemoglobins. The results demonstrate the effect of the quaternary structure of globin members on their intrinsic dynamics and show good ability to distinguish between different states of hemoglobin, revealing the dynamical relations between them.

Role of Heme Pocket Water in Allosteric Regulation of Ligand Reactivity in Human Hemoglobin

Water molecules can enter the heme pockets of unliganded myoglobins and hemoglobins, hydrogen bond with the distal histidine, and introduce steric barriers to ligand binding. The spectrokinetics of photodissociated CO complexes of human hemoglobin and its isolated α and β chains were analyzed for the effect of heme hydration on ligand rebinding. A strong coupling was observed between heme hydration and quaternary state. This coupling may contribute significantly to the 20-60-fold difference between the R- and T-state bimolecular CO binding rate constants and thus to the modulation of ligand reactivity that is the hallmark of hemoglobin allostery. Heme hydration proceeded over the course of several kinetic phases in the tetramer, including the R to T quaternary transition. An initial 150 ns hydration phase increased the R-state distal pocket water occupancy, nw(R), to a level similar to that of the isolated α (∼60%) and β (∼10%) chains, resulting in a modest barrier to ligand binding. A subsequent phase, concurrent with the first step of the R → T transition, further increased the level of heme hydration, increasing the barrier. The final phase, concurrent with the final step of the allosteric transition, brought the water occupancy of the T-state tetramer, nw(T), even higher and close to full occupancy in both the α and β subunits (∼90%). This hydration level could present an even larger barrier to ligand binding and contribute significantly to the lower iron reactivity of the T state toward CO.

Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.

An atomistic view on human hemoglobin carbon monoxide migration processes

A significant amount of work has been devoted to obtaining a detailed atomistic knowledge of the human hemoglobin mechanism. Despite this impressive research, to date, the ligand diffusion processes remain unclear and controversial. Using recently developed computational techniques, PELE, we are capable of addressing the ligand migration processes. First, the methodology was tested on myoglobin's CO migration, and the results were compared with the wealth of theoretical and experimental studies. Then, we explored both hemoglobin tense and relaxed states and identified the differences between the α-and β-subunits. Our results indicate that the proximal site, equivalent to the Xe1 cavity in myoglobin, is never visited. Furthermore, strategically positioned residues alter the diffusion processes within hemoglobin's subunits and suggest that multiple pathways exist, especially diversified in the α-globins. A significant dependency of the ligand dynamics on the tertiary structure is also observed.

Unveiling the timescale of the R-T transition in human hemoglobin

Time-resolved wide-angle X-ray scattering, a recently developed technique allowing to probe global structural changes of proteins in solution, was used to investigate the kinetics of R-T quaternary transition in human hemoglobin and to systematically compare it to that obtained with time-resolved optical spectroscopy under nearly identical experimental conditions. Our data reveal that the main structural rearrangement associated with the R-T transition takes place approximately 2 mus after the photolysis of hemoglobin at room temperature and neutral pH. This finding suggests that the 20-mus step observed with time-resolved optical spectroscopy corresponds to a small and localized structural change.

Photo-dissociation quantum yields of mammalian oxyhemoglobin investigated by a nanosecond laser technique

The photo-dissociations of oxyhemoglobin of several mammals, such as human, bovine, pig, horse, and rabbit, have been studied. By means of optical pump-probe technique, the quantum yields for photo-dissociation of these oxyhemoglobin have been determined at pH 7 and 20 degrees C. A nanosecond laser at 532 nm is used as the pumping source, and a xenon lamp through a monochrometer provides a probe light at 432 nm. The experimental results show that the quantum yields of these mammalian oxyhemoglobin are different from each other, especially for that of rabbit. By analyzing the amino acid sequences and tetramer structures as well as the flexibility and hydrophobicity of the different hemoglobin, possible explanations for the differences are proposed.

Enthalpy and conformational volume changes of mammalian oxy-hemoglobins investigated by pulsed photoacoustic calorimetry

Enthalpy and conformational volume changes induced by laser photo-dissociation reactions of mammalian oxy-hemoglobin, such as human, bovine, pig, horse and rabbit oxy-hemoglobins, are investigated by pulsed photoacoustic calorimetry. Generally, the response time of the photoacoustic calorimetry is restricted by the width of the laser pulse and the bandwidth of the acoustic detector. Considering the time window of the experimental system, the enthalpy and conformational volume changes detected should be caused by the tertiary relaxation of the heme proteins. In order to calculate the enthalpy and conformational volume changes, the quantum yields of the photo-dissociation products of oxy-hemoglobins must be measured and taken into account. Finally, the enthalpy and conformational volume changes of the oxy-hemoglobins connected with the tertiary relaxation are obtained, which show that for all measured mammalian oxy-hemoglobins the values of enthalpy changes are in the range of 30.0-46.8 kcal/mol and volume changes are of 2.3-7.8 ml/mol although the quantum yields for the different species may have much bigger differences. A possible explanation of the results is presented.

Time-resolved methods in Biophysics. 2. Monitoring haem proteins at work with nanosecond laser flash photolysis

Haem proteins have long been the most studied proteins in biophysics, and have become paradigms for the characterization of fundamental biomolecular processes as ligand binding and regulatory conformational transitions. The presence of the haem prosthetic group, the absorbance spectrum of which has a ligation sensitive region conveniently located in the UV-visible range, has offered a powerful and sensitive tool for the investigation of molecular functions. The central Fe atom is capable of reversibly binding diatomic ligands, including O(2), CO, and NO. The Fe-ligand bond is photolabile, and a reactive unligated state can be transiently generated with a pulsed laser. The photodissociated ligands quickly rebind to the haem and the process can be monitored by transient absorbance methods. The ligand rebinding kinetics reflects protein dynamics and ligand migration within the protein inner cavities. The characterization of these processes was done in the past mainly by low temperature experiments. The use of silica gels to trap proteins allows the characterization of internal ligand dynamics at room temperature. In order to show the potential of the laser flash photolysis techniques, combined with modern numerical analysis methods, we report experiments conducted on two non-symbiotic haemoglobins from Arabidopsis thaliana. The comparison between time courses recorded on haemoglobins in solution and encapsulated in silica gels allows for the highlighting of different interplays between protein dynamics and ligand migration.

Geminate rebinding in R-state hemoglobin: kinetic and computational evidence for multiple hydrophobic pockets

Biphasic geminate rebinding of CO to myoglobin upon flash photolysis has been associated to ligand distribution in hydrophobic cavities, structurally detected by time-resolved crystallography, xenon occupancy, and molecular simulations. We show that the time course of CO rebinding to human hemoglobin also exhibits a biphasic geminate rebinding when the protein is entrapped in wet nanoporous silica gel. A simple branched kinetic scheme, involving the bound state A, the primary docking site C, and a secondary binding site B was used to calculate the microscopic rates and the time-dependent population of the intermediate species. The activation enthalpies of the associated transitions were determined in the absence and presence of 80% glycerol. Potential hydrophobic docking cavities within the alpha and beta chains of hemoglobin were identified by computational modeling using xenon as a probe. A hydrophobic pocket on the distal side of the heme, corresponding to Xe4 in Mb, and a nearby site that does not have a correspondence in Mb were detected. Neither potential xenon sites on the proximal side nor a migration channel from the distal to proximal site was located. The small enthalpic barriers between states B and C are in very good agreement with the location of the xenon sites on the distal side. Furthermore, the connection between the two xenon sites is relatively open, explaining why the decreased mobility of the protein with viscosity only slightly perturbs the energetics of ligand migration between the two sites.

Evidence for Two Geminate Rebinding States Following Laser Photolysis of R State Hemoglobin Encapsulated in Wet Silica Gels

In this letter we report the first experimental evidence for CO rebinding to human hemoglobin from multiple geminate states. The analysis of the rebinding kinetics using a maximum entropy method allowed the identification of two distinct rebinding states within the protein matrix, which become populated under conditions of increased viscosity in a silica gel at high glycerol concentration. Our findings suggest the presence of at least two distinct docking sites for the photolyzed ligand. Assuming a minimal four-state model, we estimate the microscopic rates and the activation energies for the elementary processes.

Spectroscopic and Functional Characterization of T State Hemoglobin Conformations Encapsulated in Silica Gels

Oxygen binding curves of sol-gel-encapsulated deoxy human adult hemoglobin (HbA) have previously revealed two distinct noncooperative populations with oxygen binding affinities approximately 1000 and 100 times lower than that of the high-affinity R state. The two populations which have been termed the low-affinity (LA) and high-affinity (HA) T states can be selectively stabilized using two different encapsulation protocols for deoxy-HbA. The present study seeks to understand the factors giving rise to these different affinity states. Visible and UV resonance Raman spectroscopies are used to characterize the conformational properties of both the deoxy and deoxy-turned-carbonmonoxy (CO) derivatives of HbA derived from the two encapsulation protocols. The geminate and bimolecular recombination of CO to the photodissociated CO derivatives is used to characterize the functional properties of the slowly evolving encapsulated populations. The results show that the initial deoxy-HbA populations are conformationally indistinguishable with respect to encapsulation protocol. The addition of CO to sol-gel-encapsulated deoxy-HbA triggers a detectable progression of conformational and functional changes. Visible resonance Raman spectra of the CO photoproduct reveal a progression of changes of the iron-proximal histidine stretching frequencies: 215, 222, 227, and 230 cm(-1). The low and high values correspond to the initial deoxy T state and liganded R (R(2)) state species, respectively. The 222 and 227 cm(-1) species are generated using encapsulation protocols that give rise to what are termed the LA and HA T states, respectively. The UV resonance Raman spectra of these and related species indicate that the progression from deoxy T to LA to HA is associated with a progressive loosening of T state constraints within the hinge and switch regions of the alpha(1)beta(2) interface. The time scale for the progression is determined by a balance between the ligation-initiated evolution toward high-affinity conformations and factors such as allosteric effectors, gel matrix, and added glycerol that slow ligand-binding-induced relaxation. Thus, it appears that the encapsulation protocol-dependent rate of ligand-binding-induced relaxation determines the functional properties of the initially encapsulated deoxy-HbA population.

The Molecular Code for Hemoglobin Allostery Revealed by Linking the Thermodynamics and Kinetics of Quaternary Structural Change. 1. Microstate Linear Free Energy Relations

A novel model linking the thermodynamics and kinetics of hemoglobin's allosteric (R --> T) and ligand binding reactions is applied to photolysis data for human HbCO. To describe hemoglobin's kinetics at the microscopic level of structural transitions and ligand-binding events for individual [ij]-ligation microstates ((ij)R --> (ij)T, (ij)R + CO --> ((i)(+1))(k)R, and (ij)T + CO --> ((i)(+1))(k)T), the model calculates activation energies, (ij)DeltaG(++), from previously measured cooperative free energies of the equilibrium microstates (Huang, Y., and Ackers, G. K. (1996) Biochemistry 35, 704-718) by using linear free energy relations ((ij)DeltaG(++) - (01)DeltaG(++) = alpha[(ij)DeltaG - (01)DeltaG], where the parameter alpha, describing the variation of activation energy with reaction energy perturbation, can depend on the natures of both the reaction and the perturbation). The alpha value measured here for the allosteric dynamics, 0.21 +/- 0.03, corresponds closely to values observed previously, strongly suggesting that the thermodynamic microstate energies directly underlie the allosteric kinetics (as opposed to the alpha((ij)DeltaG(RT)) serving merely as arbitrary fitting parameters). Besides systematizing the study of hemoglobin kinetics, the utility of the microstate linear free energy model lies in the ability to test microscopic aspects of allosteric dynamics such as the "symmetry rule" for quaternary change deduced previously from thermodynamic evidence (Ackers, G. K., et al. (1992) Science 255, 54-63). Reflecting a remarkably detailed correspondence between thermodynamics and kinetics, we find that a kinetic model that includes the large free energy splitting between doubly ligated T microstates implied by the symmetry rule fits the data significantly better than one that does not.

Laser-induced time-resolved photoacoustic calorimetry study on photo-dissociation of human and bovine oxyhemoglobin

The dynamics of the enthalpy and volume changes related to the photo-dissociation of oxygen from human and bovine oxyhemoglobin are investigated by nanosecond time-resolved photoacoustic calorimetry (PAC). The values of enthalpy and volume change associated with the above process are deltaH = 37.8 +/- 3 kcal/mol, deltaV = 5.0 +/- 1 ml/mol for human HbO(2); and deltaH = 35.7 +/- 3.5 kcal/mol, deltaV = 4.8 +/- 1 ml/mol for bovine HbO(2), respectively. A possible explanation for the similar values between both human and bovine oxyhemoglobin is proposed. In addition, the PAC results for human HbO(2) and HbCO are compared and discussed.

Kinetics of nitric oxide binding to R-state hemoglobin

Despite earlier work indicating otherwise, some recent reports have suggested that nitric oxide (NO) binds to hemoglobin cooperatively. In particular, it has been suggested that, under physiological conditions, NO binds to the high-affinity R-state hemoglobin as much as 100 times faster than to the low-affinity T-state hemoglobin. This rapid NO binding could provide a means of preserving NO bioactivity. However, using a flash-flow photolysis technique, we have determined that the rate of NO binding to normal adult R-state hemoglobin is (2.1 +/- 0.1) x 10(7) (s(-1) M(-1)), which is essentially the same as that reported for T-state NO binding. (c)2002 Elsevier Science (USA).

The effect of water on the rate of conformational change in protein allostery

The influence of solvation on the rate of quaternary structural change is investigated in human hemoglobin, an allosteric protein in which reduced water activity destabilizes the R state relative to T. Nanosecond absorption spectroscopy of the heme Soret band was used to monitor protein relaxation after photodissociation of aqueous HbCO complex under osmotic stress induced by the nonbinding cosolute poly(ethylene glycol) (PEG). Photolysis data were analyzed globally for six exponential time constants and amplitudes as a function of osmotic stress and viscosity. Increases in time constants associated with geminate rebinding, tertiary relaxation, and quaternary relaxation were observed in the presence of PEG, along with a decrease in the fraction of hemes rebinding CO with the slow rate constant characteristic of the T state. An analysis of these results along with those obtained by others for small cosolutes showed that both osmotic stress and solvent viscosity are important determinants of the microscopic R --> T rate constant. The size and direction of the osmotic stress effect suggests that at least nine additional water molecules are required to solvate the allosteric transition state relative to the R-state hydration, implying that the transition state has a greater solvent-exposed area than either end state.