The Structure−Function Relationship of Hemoglobin in Solution at Atomic Resolution

Mechanisms of hemoglobin adaptation to high altitude hypoxia

Evidence from a number of vertebrate taxa suggests that modifications of hemoglobin (Hb) function may often play a key role in mediating an adaptive response to high altitude hypoxia. The respiratory functions of Hb are a product of the protein's intrinsic O(2)-binding affinity and its interactions with allosteric effectors such as protons, chloride ions, CO(2), and organic phosphates. Here we review several case studies involving high altitude vertebrates where it has been possible to identify specific mechanisms of Hb adaptation to hypoxia. In addition to comparative studies of Hbs from diverse animal species, functional studies of human Hb mutants also suggest that there is ample scope for evolutionary adjustments in Hb-O(2) affinity through alterations of the equilibrium constants of O(2) binding to deoxy- and oxyHb or through changes in the allosteric equilibrium constants for the transition between the deoxy- and oxyHb quaternary structures. It may be the case that certain evolutionary paths are followed more often than others simply because they are subject to less stringent pleiotropic constraints.

High-altitude adaptations in vertebrate hemoglobins

Vertebrates at high altitude are subjected to hypoxic conditions that challenge aerobic metabolism. O(2) transport from the respiratory surfaces to tissues requires matching between the O(2) loading and unloading tensions and the O(2)-affinity of blood, which is an integrated function of hemoglobin's intrinsic O(2)-affinity and its allosteric interaction with cellular effectors (organic phosphates, protons and chloride). Whereas short-term altitudinal adaptations predominantly involve adjustments in allosteric interactions, long-term, genetically-coded adaptations typically involve changes in the structure of the haemoglobin molecules. The latter commonly comprise substitutions of amino acid residues at the effector binding sites, the heme-protein contacts, or at intersubunit contacts that stabilize either the low-affinity ('Tense') or the high-affinity ('Relaxed') structures of the molecules. Molecular heterogeneity (multiple isoHbs with differentiated oxygenation properties) can further broaden the range of physico-chemical conditions where Hb functions under altitudinal hypoxia. This treatise reviews the molecular and cellular mechanisms that adapt haemoglobin-oxygen affinities in mammals, birds and ectothermic vertebrates at high altitude.

A convergent coordination chemistry-based approach to dissymmetric macrocyclic cofacial porphyrin complexes

We report a highly convergent and modular approach for the synthesis of dissymmetric cofacial porphyrin complexes, which is based upon the weak-link approach to supramolecular coordination chemistry. Specifically, we have utilized a halide-induced ligand rearrangement reaction, which is capable of providing heteroligated mixed-metal porphyrin complexes in quantitative yield. Significantly, the adoption of a coordination chemistry based approach for the synthesis of these complexes allows for facile in situ regulation of the porphyrin-porphyrin interactions through the addition of external chemical stimuli.

Insights into the solution structure of human deoxyhemoglobin in the absence and presence of an allosteric effector

We present a nuclear magnetic resonance (NMR) study in solution of the structures of human normal hemoglobin (Hb A) in the deoxy or unligated form in the absence and presence of an allosteric effector, inositol hexaphosphate (IHP), using 15N-1H residual dipolar coupling (RDC) measurements. There are several published crystal structures for deoxyhemoglobin A (deoxy-Hb A), and it has been reported that the functional properties of Hb A in single crystals are different from those in solution. Carbonmonoxyhemoglobin A (HbCO A) can also be crystallized in several structures. Our recent RDC studies of HbCO A in the absence and presence of IHP have shown that the solution structure of this Hb molecule is distinctly different from its classical crystal structures (R and R2). To have a better understanding of the structure-function relationship of Hb A under physiological conditions, we need to evaluate its structures in both ligated and unligated states in solution. Here, the intrinsic paramagnetic property of deoxy-Hb A has been exploited for the measurement of RDCs using the magnetic-field dependence of the apparent one-bond 1H-15N J couplings. Our RDC analysis suggests that the quaternary and tertiary structures of deoxy-Hb A in solution differ from its recently determined high-resolution crystal structures. Upon binding of IHP, structural changes in deoxy-Hb A are also observed, and these changes are largely within the alpha1beta1 (or alpha2beta2) dimer itself. These new structural findings allow us to gain a deeper insight into the structure-function relationship of this interesting allosteric protein.

Orientation of deoxyhemoglobin at high magnetic fields: structural insights from RDCs in solution

Human normal adult hemoglobin (Hb A) is a tetrameric protein molecule of ~64 kDa consisting of two identical -chains and two identical -chains of 141 and 146 amino acid residues each and four bound heme moieties. In the oxygen-free form of Hb A, also known as deoxyhemoglobin A (deoxy-Hb A), the hemes are paramagnetic with S = 2. We have measured the one-bond spin-spin couplings (1JNH + 1DNH) on (15N,2H)-labeled deoxy-Hb A in solution as a function of magnetic field strengths from 11.7 to 21.1 T and found that these couplings are linearly proportional to the square of the magnetic field. This field dependence provides an opportunity to extract the residual dipolar couplings (RDCs, 1DNH) and, thus, to compare predictions about the solution structure of deoxy-Hb A to crystal structures for this molecule. Such comparison is essential for our understanding of the structure, dynamics, and function of this allosteric protein under conditions close to the physiological state. This report illustrates the usefulness of using the magnetic-field dependent RDCs to determine the solution structure of a large paramagnetic protein. This method is especially valuable for those proteins whose structures must be determined in an oxygen-free environment.

Historical reconstructions of evolving physiological complexity:O2 secretion in the eye and swimbladder of fishes

The ability of some fishes to inflate their compressible swimbladder with almost pure oxygen to maintain neutral buoyancy, even against the high hydrostatic pressure several thousand metres below the water surface, has fascinated physiologists for more than 200 years. This review shows how evolutionary reconstruction of the components of such a complex physiological system on a phylogenetic tree can generate new and important insights into the origin of complex phenotypes that are difficult to obtain with a purely mechanistic approach alone. Thus, it is shown that oxygen secretion first evolved in the eyes of fishes, presumably for improved oxygen supply to an avascular, metabolically active retina. Evolution of this system was facilitated by prior changes in the pH dependence of oxygen-binding characteristics of haemoglobin (the Root effect) and in the specific buffer value of haemoglobin. These changes predisposed teleost fishes for the later evolution of swimbladder oxygen secretion, which occurred at least four times independently and can be associated with increased auditory sensitivity and invasion of the deep sea in some groups. It is proposed that the increasing availability of molecular phylogenetic trees for evolutionary reconstructions may be as important for understanding physiological diversity in the postgenomic era as the increase of genomic sequence information in single model species.

Intracellular pH sensors: design principles and functional significance

Changes in intracellular pH regulate many cell behaviors, including proliferation, migration, and transformation. However, our understanding of how physiological changes in pH affect protein conformations and macromolecular assemblies is limited. We present design principles, current modeling predictions, and examples of pH sensors or proteins that have activities or ligand-binding affinities that are regulated by changes in intracellular pH.

Quantitative vibrational dynamics of iron in carbonyl porphyrins

We use nuclear resonance vibrational spectroscopy and computational predictions based on density functional theory (DFT) to explore the vibrational dynamics of (57)Fe in porphyrins that mimic the active sites of histidine-ligated heme proteins complexed with carbon monoxide. Nuclear resonance vibrational spectroscopy yields the complete vibrational spectrum of a Mössbauer isotope, and provides a valuable probe that is not only selective for protein active sites but quantifies the mean-squared amplitude and direction of the motion of the probe nucleus, in addition to vibrational frequencies. Quantitative comparison of the experimental results with DFT calculations provides a detailed, rigorous test of the vibrational predictions, which in turn provide a reliable description of the observed vibrational features. In addition to the well-studied stretching vibration of the Fe-CO bond, vibrations involving the Fe-imidazole bond, and the Fe-N(pyr) bonds to the pyrrole nitrogens of the porphyrin contribute prominently to the observed experimental signal. All of these frequencies show structural sensitivity to the corresponding bond lengths, but previous studies have failed to identify the latter vibrations, presumably because the coupling to the electronic excitation is too small in resonance Raman measurements. We also observe the FeCO bending vibrations, which are not Raman active for these unhindered model compounds. The observed Fe amplitude is strongly inconsistent with three-body oscillator descriptions of the FeCO fragment, but agrees quantitatively with DFT predictions. Over the past decade, quantum chemical calculations have suggested revised estimates of the importance of steric distortion of the bound CO in preventing poisoning of heme proteins by carbon monoxide. Quantitative agreement with the predicted frequency, amplitude, and direction of Fe motion for the FeCO bending vibrations provides direct experimental support for the quantum chemical description of the energetics of the FeCO unit.

Structure-function relationship in a variant hemoglobin: a combined computational-experimental approach

Our study examines the functional and structural effects of amino acid substitution in the distal side of beta-chains of human Hb Duarte (alpha(2)beta(2)(62Ala-->Pro)). We have compared the functional properties of the purified Hb Duarte with those of HbA, and through proton NMR and molecular dynamics simulations we have investigated their tertiary and quaternary structures. The variant exhibits an increased oxygen affinity with a normal Hill coefficient and Bohr effect. The abnormal function of Hb Duarte is attributed to the presence of a proline residue at the beta62 position, since the functional properties of another Hb variant in the same position, Hb J-Europa (beta(62Ala-->Asp)), have been described as normal. Thereafter (1)H-NMR studies have shown that the beta62 Ala-->Pro substitution causes structural modifications of the tertiary structure of the beta globins, leaving the quaternary structure unaltered. These results have been confirmed by extensive all-atom molecular dynamics simulations. All these findings lead to the conclusion that the beta62 Ala-->Pro substitution produces a destabilization of the E-helix extending downward to the CD corner. Particularly, a cavity near the distal histidine of the beta-chains, connecting the heme pocket to the solvent, is affected, altering the functional properties of the protein molecule.

1.25 Å Resolution Crystal Structures of Human Haemoglobin in the Oxy, Deoxy and Carbonmonoxy Forms

The most recent refinement of the crystallographic structure of oxyhaemoglobin (oxyHb) was completed in 1983, and differences between this real-space refined model and later R state models have been interpreted as evidence of crystallisation artefacts, or numerous sub-states. We have refined models of deoxy, oxy and carbonmonoxy Hb to 1.25 A resolution each, and compare them with other Hb structures. It is shown that the older structures reflect the software used in refinement, and many differences with newer structures are unlikely to be physiologically relevant. The improved accuracy of our models clarifies the disagreement between NMR and X-ray studies of oxyHb, the NMR experiments suggesting a hydrogen bond to exist between the distal histidine and oxygen ligand of both the alpha and beta-subunits. The high-resolution crystal structure also reveals a hydrogen bond in both subunit types, but with subtly different geometry which may explain the very different behaviour when this residue is mutated to glycine in alpha or beta globin. We also propose a new set of relatively fixed residues to act as a frame of reference; this set contains a similar number of atoms to the well-known "BGH" frame yet shows a much smaller rmsd value between R and T state models of HbA.

Evolution of vertebrate haemoglobins: Histidine side chains, specific buffer value and Bohr effect

This review highlights the use of analytical tools, recently developed in the comparative method of evolutionary biology, for the study of haemoglobin (Hb) adaptation. It focuses on the functional consequences of a previously largely ignored structural feature of Hb, namely the degree and positional specificity of histidine (His) substitution in Hb chains. The importance of His side chains for hydrogen ion buffering, blood CO(2) transport capacity and the molecular mechanism of the Bohr effect in vertebrate Hbs is discussed. Using phylogenetically independent contrasts, a significant correlation between the specific buffer value of Hb and the number of predicted physiological buffer groups from Hb sequence data is shown. In a new result, the evolution of the number of physiological buffer groups in 77 vertebrate species is reconstructed on a phylogenetic tree. The analysis predicts that teleost fishes, passeriform birds and some snakes have independently evolved a much-reduced specific buffer value of Hb, possibly for enhancing the efficiency of an acid load to change oxygen affinity via the Bohr effect. This analysis demonstrates how in comparative physiology analysis of genetic databases in an evolutionary framework can identify candidate species for further experimental in vitro and whole animal studies.

Dynamics of Allostery in Hemoglobin: Roles of the Penultimate Tyrosine H bonds

The tyrosine residues adjacent to the C termini of the hemoglobin (Hb) subunits, alphaY140 and betaY145, are expected to play important structural roles, because the C termini are the loci of T-state quaternary salt-bridges, and because the tyrosine side-chains bridge the H and F helices via H bonds to the alphaV93 and betaV98 carbonyl groups. These roles have been investigated via measurements of oxygen binding, (1)H NMR spectra, resonance Raman (RR) spectra, and time-resolved resonance Raman (TR(3)) spectra on site mutants in which the Hcdots, three dots, centeredF H bonds are eliminated by replacing the tyrosine residues with phenylalanine. The TR(3) spectra confirm the hypothesis, based on TR(3) studies of wild-type Hb, that the Hcdots, three dots, centeredF H bonds break and then re-form during the sub-microsecond phase of the R-T quaternary transition. The TR(3) spectra support the inference from other mutational studies that the alphabeta dimers act as single dynamic units in this early phase, motions of the E and F helices being coupled tightly across the dimer interface. Formation of T quaternary contacts occurs at about the same rate in the mutants as in HbA. However, these contacts are weakened substantially by the Y/F substitutions. Equilibrium perturbations are apparent also, especially for the alpha-subunits, in which relaxation of the Fe-His bond, strengthening of the Acdots, three dots, centeredE interhelical H bond, and weakening of the "switch" quaternary contact in deoxyHb are all apparent. Structural effects are less marked for the beta-chain Y/F replacement, but the Bohr effect is reduced by 25%, indicating that the salt-bridge and H bond interactions of the adjacent C terminus are loosened. The alpha-chain replacement reduces the Bohr effect much more, consistent with the global perturbations detected by the structure probes.

Correlation between computed gas-phase and experimentally determined solution-phase infrared spectra: Models of the iron–iron hydrogenase enzyme active site

Gas-phase density functional theory calculations (B3LYP, double zeta plus polarization basis sets) are used to predict the solution-phase infrared spectra for a series of CO- and CN-containing iron complexes. It is shown that simple linear scaling of the computed C--O and C--N stretching frequencies yields accurate predictions of the the experimentally determined nu(CO) and nu(CN) values for a variety of complexes of different charges and in solvents of varying polarity. As examples of the technique, the resulting correlation is used to assign structures to spectroscopically observed but structurally ambiguous species in two different systems. For the (mu-SCH2CH2CH2S)[Fe(CO)3]2 complex in tetrahydrofuran solution, our calculations show that the initial electrochemical reduction process leads to a simple one-electron reduced product with a structure very similar to the (mu-SCH2CH2CH2S)[Fe(CO)3]2 parent complex. For the iron-iron hydrogenase enzyme active site, our computations show that the absence or presence of a water molecule near the distal iron center (the iron center further from the [4Fe4S] cluster and protein backbone) has very little effect on the predicted infrared spectra.

R-state haemoglobin with low oxygen affinity: crystal structures of deoxy human and carbonmonoxy horse haemoglobin bound to the effector molecule L35

Although detailed crystal structures of haemoglobin (Hb) provide a clear understanding of the basic allosteric mechanism of the protein, and how this in turn controls oxygen affinity, recent experiments with artificial effector molecules have shown a far greater control of oxygen binding than with natural heterotropic effectors. Contrary to the established text-book view, these non-physiological compounds are able to reduce oxygen affinity very strongly without switching the protein to the T (tense) state. In an earlier paper we showed that bezafibrate (BZF) binds to a surface pocket on the alpha subunits of R state Hb, strongly reducing the oxygen affinity of this protein conformation. Here we report the crystallisation of Hb with L35, a related compound, and show that this binds to the central cavity of both R and T state Hb. The mechanism by which L35 reduces oxygen affinity is discussed, in relation to spectroscopic studies of effector binding.

Proton-linked subunit heterogeneity in ferrous nitrosylated human adult hemoglobin: an EPR study

The effect of pH on the X-band electron paramagnetic resonance (EPR) spectrum of ferrous nitrosylated human adult tetrameric hemoglobin (HbNO) as well as of ferrous nitrosylated monomeric alpha- and beta-chains has been investigated, at -163 degrees C. At pH 7.3, the X-band EPR spectrum of tetrameric HbNO and ferrous nitrosylated monomeric alpha- and beta-chains displays a rhombic shape. Lowering the pH from 7.3 to 3.0, tetrameric HbNO and ferrous nitrosylated monomeric alpha- and beta-chains undergo a transition towards a species characterized by a X-band EPR spectrum with a three-line splitting centered at 334mT. These pH-dependent spectroscopic changes may be taken as indicative of the cleavage, or the severe weakening, of the proximal HisF8-Fe bond. In tetrameric HbNO, the pH-dependent spectroscopic changes depend on the acid-base equilibrium of two apparent ionizing groups with pK(a) values of 5.8 and 3.8. By contrast, the pH-dependent spectroscopic changes occurring in ferrous nitrosylated monomeric alpha- and beta-chains depend on the acid-base equilibrium of one apparent ionizing group with pK(a) values of 4.8 and 4.7, respectively. The different pK(a) values for the proton-linked spectroscopic transition(s) of tetrameric HbNO and ferrous nitrosylated monomeric alpha- and beta-chains suggest that the quaternary assembly drastically affects the strength of the proximal HisF8-Fe bond in both subunits. This probably reflects a 'quaternary effect', i.e., structural changes in both subunits upon tetrameric assembly, which is associated to a relevant variation of functional properties (i.e., proton affinity).

Side-chain assignments of methyl-containing residues in a uniformly 13C-labeled hemoglobin in the carbonmonoxy form

Sequence-specific assignment of the methyl groups in large proteins can be obtained from an MQ-(H)CC(m)H(m)-TOCSY experiment on uniformly (13)C-labeled proteins without deuteration (Yang et al., 2004). Here the procedure is further demonstrated on a uniformly (13)C-labeled alpha-chain or beta-chain of human normal adult hemoglobin (65 kDa) in the carbonmonoxy form. In addition, a strategy is presented for assigning protons of methyl-containing residues of uniformly (13)C-labeled large proteins, on the basis of prior methyl assignments based on MQ-(H)CCH-TOCSY and H(C)C(m)H(m)-TOCSY experiments. Assignment of about 80% of the side-chain resonances of methyl-containing residues of carbonmonoxyhemoglobin has been obtained.

Red blood cell pH, the Bohr effect, and other oxygenation-linked phenomena in blood O2 and CO2 transport

The discovery of the S-shaped O2 equilibrium curve and the Bohr effect in 1904 stimulated a fertile and continued research into respiratory functions of blood and allosteric mechanisms in haemoglobin (Hb). The Bohr effect (influence of pH/CO2 on Hb O2 affinity) and the reciprocal Haldane effect (influence of HbO2 saturation on H+/CO2 binding) originate in the Hb oxy-deoxy conformational change and allosteric interactions between O2 and H+/CO2 binding sites. In steady state, H+ is passively distributed across the vertebrate red blood cell (RBC) membrane, and intracellular pH (pHi) changes are related to changes in extracellular pH, Hb-O2 saturation and RBC organic phosphate content. As the Hb molecule shifts between the oxy and deoxy conformation in arterial-venous gas transport, it delivers O2 and takes up CO2 and H+ in tissue capillaries (elegantly aided by the Bohr effect). Concomitantly, the RBC may sense local O2 demand via the degree of Hb deoxygenation and release vasodilatory agents to match local blood flow with requirements. Three recent hypotheses suggest (1) release of NO from S-nitroso-Hb upon deoxygenation, (2) reduction of nitrite to vasoactive NO by deoxy haems, and (3) release of ATP. Inside RBCs, carbonic anhydrase (CA) provides fast hydration of metabolic CO2 and ensures that the Bohr shift occurs during capillary transit. The formed H+ is bound to Hb (Haldane effect) while HCO3- is shifted to plasma via the anion exchanger (AE1). The magnitude of the oxylabile H+ binding shows characteristic differences among vertebrates. Alternative strategies for CO2 transport include direct HCO3- binding to deoxyHb in crocodilians, and high intracellular free [HCO3-] (due to high pHi) in lampreys. At the RBC membrane, CA, AE1 and other proteins may associate into what appears to be an integrated gas exchange metabolon. Oxygenation-linked binding of Hb to the membrane may regulate glycolysis in mammals and perhaps also oxygen-sensitive ion transport involved in RBC volume and pHi regulation. Blood O2 transport shows several adaptive changes during exposure to environmental hypoxia. The Bohr effect is involved via the respiratory alkalosis induced by hyperventilation, and also via the pHi change that results from modulation of RBC organic phosphate content. In teleost fish, beta-adrenergic activation of Na+/H+ exchange rapidly elevates pHi and O2 affinity, particularly under low O2 conditions.

Solution 1H NMR study of the active site molecular structure and magnetic properties of the cyanomet complex of the isolated, tetrameric beta-chain from human adult hemoglobin

The solution molecular structure and the electronic and magnetic properties of the heme pocket of the cyanomet complex of the isolated beta-chain of human adult hemoglobin, HbA, have been investigated by homonuclear 2D (1)H NMR in order to assess the extent of assignments allowed by (1)H NMR of a homo-tetrameric 65-kDa protein, to guide the future assignments of the heterotetrameric complex of HbA, and to compare the structure of the beta-chain to the crystallographically characterized complexes that contains the beta-chain. The target residues are those that exhibit significant (>|0.2| ppm) dipolar shifts, as predicted by a "preliminary" set of magnetic axes determined from a small set of easily assigned active site residues. All 104 target residues ( approximately 70% of total) were assigned by taking advantage of the temperature dependence predicted by the "preliminary" magnetic axes for the polypeptide backbone; they include all residues proposed to play a significant role in modulating the ligand affinity in the tetramer HbA. Left unassigned are the A-helix, the end of the G-helix and the beginning of the H-helix where dipolar shifts are less than |0.2| ppm. These comprehensive assignments allow the determination of a robust set of orientation and anisotropies of the paramagnetic susceptibility tensor that leads to quantitative interpretation of the dipolar shifts of the beta-chain in terms of the crystal coordinates of the beta-subunit in ligated HbA which, in turn, confirms a largely conserved molecular structure of the isolated beta-chain relative to that in the intact R-state HbA. The major magnetic axis, which is correlated with the tilt of the Fe-CN unit, is tilted approximately 10 degrees from the heme normal so that the Fe-CN unit is tilted toward the beta-meso-H in a fashion remarkably similar to the Fe-CO tilt in the beta-subunit of HbCO. It is concluded that a set of "preliminary" magnetic axes and the use of variable temperature 2D NMR spectra are crucial to effective assignments in the tetrameric cyanomet beta-chain and that this approach should be similarly effective in HbA.