Molecular dynamics of human methemoglobin: the transmission of conformational information between subunits in an alpha beta dimer

Atomistic Simulations of Heme Dissociation Pathways in Human Methemoglobins Reveal Hidden Intermediates

Heme dissociations disrupt function and structural integrity of human hemoglobin and trigger various cardiovascular complications. These events become significant in methemoglobins that have undergone autoxidation of ferrous into ferric heme. We have structurally characterized the heme disassociation pathways for adult tetrameric methemoglobins using all-atom molecular dynamics simulations. These reveal that bis-histidine hemichromes, characterized here by the coordination of heme iron to both the F8 (proximal) and E7 (distal) histidines, are seen as intermediates following dissociation of the water molecule distally bound to each heme iron. Later, the breaking of coordination between heme iron and proximal histidine disrupts the F helix and pushes it away from the heme cavity, enabling both bulk solvent penetration and disruption of tetramer interface interactions. The interactions inhibiting heme dissociation were then seen to be (i) either a direct or a water-molecule-mediated interaction between distal histidine and heme iron and (ii) stacking between heme and the αCE1/βCD1 phenylalanine residue. These interactions are less important in the β than in α subunits due to a more flexible β subunit CE loop region. The absence of a distal histidine interaction in the H(E7)L mutant and increased heme cavity volume in the V(E11)A mutant both promoted heme escape from the protein interior. Adult and fetal hemoglobins were seen to share a general heme disassociation pathway and intermediates due to the conservation of key heme pocket residues. The intermediates seen here are analyzed in light of experimental studies of heme dissociation and pathways of certain hemoglobinopathies.

Collective dynamics underlying allosteric transitions in hemoglobin

Hemoglobin is the prototypic allosteric protein. Still, its molecular allosteric mechanism is not fully understood. To elucidate the mechanism of cooperativity on an atomistic level, we developed a novel computational technique to analyse the coupling of tertiary and quaternary motions. From Molecular Dynamics simulations showing spontaneous quaternary transitions, we separated the transition trajectories into two orthogonal sets of motions: one consisting of intra-chain motions only (referred to as tertiary-only) and one consisting of global inter-chain motions only (referred to as quaternary-only). The two underlying subspaces are orthogonal by construction and their direct sum is the space of full motions. Using Functional Mode Analysis, we were able to identify a collective coordinate within the tertiary-only subspace that is correlated to the most dominant motion within the quaternary-only motions, hence providing direct insight into the allosteric coupling mechanism between tertiary and quaternary conformation changes. This coupling-motion is substantially different from tertiary structure changes between the crystallographic structures of the T- and R-state. We found that hemoglobin's allosteric mechanism of communication between subunits is equally based on hydrogen bonds and steric interactions. In addition, we were able to affect the T-to-R transition rates by choosing different histidine protonation states, thereby providing a possible atomistic explanation for the Bohr effect.

Spontaneous quaternary and tertiary T-R transitions of human hemoglobin in molecular dynamics simulation

We present molecular dynamics simulations of unliganded human hemoglobin (Hb) A under physiological conditions, starting from the R, R2, and T state. The simulations were carried out with protonated and deprotonated HC3 histidines His(β)146, and they sum up to a total length of 5.6µs. We observe spontaneous and reproducible T→R quaternary transitions of the Hb tetramer and tertiary transitions of the α and β subunits, as detected from principal component projections, from an RMSD measure, and from rigid body rotation analysis. The simulations reveal a marked asymmetry between the α and β subunits. Using the mutual information as correlation measure, we find that the β subunits are substantially more strongly linked to the quaternary transition than the α subunits. In addition, the tertiary populations of the α and β subunits differ substantially, with the β subunits showing a tendency towards R, and the α subunits showing a tendency towards T. Based on the simulation results, we present a transition pathway for coupled quaternary and tertiary transitions between the R and T conformations of Hb.

As the prototypic allosteric protein, human hemoglobin (Hb) has drawn extensive scientific efforts for many decades. Human Hb exists in two quaternary conformations, the low-affinity (or deoxy) T state, and the high-affinity (or oxy) R state, and the transition between the T and the R state is mainly characterized by a 15° rotation of the α1β1 dimer with respect to the α1β2 dimer. Subsequent binding of molecular oxygen to the four heme groups drives the Hb tetramer from the deoxy T to the oxy R state, rendering the T→R transitions the molecular process that underlies the well-known cooperativity of ligand binding. Despite the wealth of structural information available for Hb, the mechanistic coupling between the quaternary transition and the tertiary transitions in the individual subunits remains poorly understood. We report spontaneous and reproducible T-R transitions of Hb in molecular dynamics simulations, allowing us to study the mechanism underlying the transitions in atomistic detail. We pay special attention to the interplay between the tertiary and quaternary transitions of Hb, as well as to the tertiary t/r populations of the subunits in a particular (T or R) quaternary state. Interestingly, we observe a pronounced asymmetry between the α and β subunits with the β subunits being more strongly linked to the quaternary transitions than the α subunits. The simulations allow us to propose a pathway for coupled quaternary and tertiary transitions between the R and T conformations of Hb.

Toxicological consequences of extracellular hemoglobin: biochemical and physiological perspectives

Under normal physiology, human red blood cells (RBCs) demonstrate a circulating lifespan of approximately 100-120 days with efficient removal of senescent RBCs taking place via the reticuloendothelial system, spleen, and bone marrow phagocytosis. Within this time frame, hemoglobin (Hb) is effectively protected by efficient RBC enzymatic systems designed to allow for interaction between Hb and diffusible ligands while preventing direct contact between Hb and the external environment. Under normal resting conditions, the concentration of extracellular Hb in circulation is therefore minimal and controlled by specific plasma and cellular (monocyte/macrophage) binding proteins (haptoglobin) and receptors (CD163), respectively. However, during pathological conditions leading to hemolysis, extracellular Hb concentrations exceed normal plasma and cellular binding capacities, allowing Hb to become a biologically relevant vasoactive and redox active protein within the circulation and at extravascular sites. Under conditions of genetic, drug-induced, and autoimmune hemolytic anemias, large quantities of Hb are introduced into the circulation and often lead to acute renal failure and vascular dysfunction. Interestingly, the study of chemically modified Hb for use as oxygen therapeutics has allowed for some basic understanding of extracellular Hb toxicity, particularly in the absence of functional clearance mechanisms and in circulatory antioxidant depleted states.

Comparison of nitric oxide-induced oxidation of recombinant oxyhemoglobin subunits using a competition experiment

A low reaction rate with nitric oxide (NO) is one of the important characteristics of hemoglobin (Hb)-based oxygen carriers. The reaction rate between oxyHb and NO is usually measured by stopped-flow spectrophotometry. However, the reported rates vary due to the difficulty of accurately determining the NO concentration and the limit of the instrument dead time. To circumvent these problems, we developed an experiment using oxymyoglobin (oxyMb) to compete with oxyHb for NO that is released from an NO donor. Determination of the rate constants in the competition experiment no longer depends on accurate measurement of time or NO concentration, since this approach instead measures the ratio of rate constants for the reaction of oxyHb and oxyMb with NO. For recombinant mutant Hb alpha(L29F)beta the rates for alpha(L29F) and beta are approximately 15- and 1.6-fold smaller than for wild-type Hb. In conclusion, the competition experiment provides an alternative method for determination of relative reaction rates of recombinant Hb subunits with NO.

Subunit-selective interrogation of CO recombination in carbonmonoxy hemoglobin by isotope-edited time-resolved resonance Raman spectroscopy

Hemoglobin (Hb) is an allosteric tetrameric protein made up of alphabeta heterodimers. The alpha and beta chains are similar, but are chemically and structurally distinct. To investigate dynamical differences between the chains, we have prepared tetramers in which the chains are isotopically distinguishable, via reconstitution with (15)N-heme. Ligand recombination and heme structural evolution, following HbCO dissociation, was monitored with chain selectivity by resonance Raman (RR) spectroscopy. For alpha but not for beta chains, the frequency of the nu(4) porphyrin breathing mode increased on the microsecond time scale. This increase is a manifestation of proximal tension in the Hb T-state, and its time course is parallel to the formation of T contacts, as determined previously by UVRR spectroscopy. Despite the localization of proximal constraint in the alpha chains, geminate recombination was found to be equally probable in the two chains, with yields of 39 +/- 2%. We discuss the possibility that this equivalence is coincidental, in the sense that it arises from the evolutionary pressure for cooperativity, or that it reflects mechanical coupling across the alphabeta interface, evidence for which has emerged from UVRR studies of site mutants.

Protein dynamics from time resolved UV Raman spectroscopy

Raman spectroscopy can provide unique information on the evolution of structure in proteins over a wide range of time scales; the picosecond to millisecond range can be accessed with pump-probe techniques. Specific parts of the molecule are interrogated by tuning the probe laser to a resonant electronic transition, including the UV transitions of aromatic residues and of the peptide bond. Advances in laser technology have enabled the characterization of transient species at an unprecedented level of structural detail. Applications to protein unfolding and allostery are reviewed.

Red blood cells induce hypoxic lung inflammation

Hypoxia, which commonly associates with respiratory and cardiovascular diseases, provokes an acute inflammatory response. However, underlying mechanisms are not well understood. Here we report that red blood cells (RBCs) induce hypoxic inflammation by producing reactive oxygen species (ROS) that diffuse to endothelial cells of adjoining blood vessels. Real-time fluorescence imaging of rat and mouse lungs revealed that in the presence of RBC-containing vascular perfusion, hypoxia increased microvascular ROS, and cytosolic Ca(2+), leading to P-selectin-dependent leukocyte recruitment. However, in the presence of RBC-free perfusion, all hypoxia-induced responses were completely inhibited. Because hemoglobin (Hb) autoxidation causes RBC superoxide formation that readily dismutates to H(2)O(2), hypoxia-induced responses were lost when we inhibited Hb autoxidation with CO or nitrite, or when the H(2)O(2) inhibitor, catalase was added to the infusion to neutralize the RBC-derived ROS. By contrast, perfusion with RBCs from BERK-trait mice that are more susceptible to Hb autoxidation and to hypoxia-induced superoxide production enhanced the hypoxia-induced responses. We conclude that in hypoxia, increased Hb autoxidation augments superoxide production in RBCs. Consequently, RBCs release H(2)O(2) that diffuses to the lung microvascular endothelium, thereby initiating Ca(2+)-dependent leukocyte recruitment. These findings are the first evidence that RBCs contribute to hypoxia-induced inflammation.

Protein dynamics explain the allosteric behaviors of hemoglobin

Bohr, Hasselbalch, and Krogh discovered homotropic and heterotropic allosteric behaviors of hemoglobin (Hb) in 1903/1904. A chronological description since then of selected principal models of the allosteric mechanism of Hb, such as the Adair scheme, the MWC two-state concerted model, the KNF induced-fit sequential model, the Perutz stereochemical model, the tertiary two-state model, and the global allostery model (an expanded MWC models), is concisely presented, followed by analysis and discussion of their limitations and deficiencies. The determination of X-ray crystallographic structures of deoxy- and ligated-Hb and the structure-based stereochemical model by Perutz are an epoch-making event in this history. However, his assignment of low-affinity deoxy- and high-affinity oxy-quaternary structures of Hb to the T- and R-states, respectively, though apparently reasonable, and as well as his hypothesis that the T-/R-quaternary structural transition regulates the oxygen-affinity, have created confusions and side-tracked studies of Hb on the structure-function relationship. The differences in static molecular structures of Hb between T(deoxy)- and R(oxy)-quaternary states reported in detail by Perutz and others are ligation-linked structural changes, but not related to the control/regulation of the oxygen-affinity. The oxygen-affinity (K(T) and K(R)) of Hb has been shown to be regulated by the heterotropic effector-linked tertiary structural changes without involving the T/R-quaternary changes. However, a recent high-resolution crystallographic analysis of Hb with different oxygen-affinities shows that static molecular structures of Hb determined by crystallography can neither identify the nature of the T(low-affinity) functional state nor decipher the mechanism by which Hb stores free energy in the T(low-affinity) functional state. Molecular dynamics simulations show that fluctuations of helices of oxy-Hb are increased upon de-oxygenation and/or binding 2,3-biphosphoglycerate. These are known to lower the oxygen-affinity of Hb. It is proposed that the coordination mode of the heme Fe with proximal and distal His is modulated by these helical fluctuations, resulting in the modulation of the oxygen-affinity of Hb. Therefore, it is proposed that the oxygen-affinity of Hb is regulated by pentanary (the 5th-order time-dependent or dynamic) tertiary structural changes rather than the T-/R-quaternary structural transitions in Hb. Homotropic and heterotropic allosteric effects of Hb are oxygen- and effector-linked, conformational entropy-driven entropy-enthalpy compensation phenomena and not much to do with static structural changes. The dynamic allostery model, which integrates these observations, provides the structural basis for the global allostery model (an expanded MWC model).

Effector-induced structural fluctuation regulates the ligand affinity of an allosteric protein: binding of inositol hexaphosphate has distinct dynamic consequences for the T and R states of hemoglobin

The present study reports distinct dynamic consequences for the T- and R-states of human normal adult hemoglobin (Hb A) due to the binding of a heterotropic allosteric effector, inositol hexaphosphate (IHP). A nuclear magnetic resonance (NMR) technique based on modified transverse relaxation optimized spectroscopy (TROSY) has been used to investigate the effect of conformational exchange of Hb A in both deoxy and CO forms, in the absence and presence of IHP, at 14.1 and 21.1 T, and at 37 degrees C. Our results show that the majority of the polypeptide backbone amino acid residues of deoxy- and carbonmonoxy-forms of Hb A in the absence of IHP is not mobile on the micros-ms time scale, with the exception of several amino acid residues, that is, beta109Val and beta132Lys in deoxy-Hb A, and alpha40Lys in HbCO A. The mobility of alpha40Lys in HbCO A can be explained by the crystallographic data showing that the H-bond between alpha40Lys and beta146His in deoxy-Hb A is absent in HbCO A. However, the conformational exchange of beta109Val, which is located in the intradimer (alpha 1beta 1 or alpha 2beta 2) interface, is not consistent with the crystallographic observations that show rigid packing at this site. IHP binding appears to rigidify alpha40Lys in HbCO A, but does not significantly affect the flexibility of beta109Val in deoxy-Hb A. In the presence of IHP, several amino acid residues, especially those at the interdimer (alpha 1beta 2 or alpha 2beta 1) interface of HbCO A, exhibit significant conformational exchange. The affected residues include the proximal beta92His in the beta-heme pocket, as well as some other residues located in the flexible joint (betaC helix-alphaFG corner) and switch (alphaC helix-betaFG corner) regions that play an important role in the dimer-dimer rotation of Hb during the oxygenation process. These findings suggest that, upon IHP binding, HbCO A undergoes a conformational fluctuation near the R-state but biased toward the T-state, apparently along the trajectory of its allosteric transition, accompanied by structural fluctuations in the heme pocket of the beta-chain. In contrast, no significant perturbation of the dynamic features on the ms-micros time scale has been observed upon IHP binding to deoxy-Hb A. We propose that the allosteric effector-induced quaternary structural fluctuation may contribute to the reduced ligand affinity of ligated hemoglobin. Conformational exchange mapping of the beta-chain of HbCO A observed at 21.1 T shows significantly increased scatter in the chemical exchange contribution to the transverse relaxation rate ( R ex) values, relative to those at lower fields, due to the enhanced effect of the local chemical shift anisotropy (CSA) fluctuation. A spring-on-scissors model is proposed to interpret the dynamic phenomena induced by the heterotropic effector, IHP.

Molecular dynamics simulations of hemoglobin A in different states and bound to DPG: effector-linked perturbation of tertiary conformations and HbA concerted dynamics

Recent functional studies reported on human adult hemoglobin (HbA) show that heterotropic effector-linked tertiary structural changes are primarily responsible for modulating the oxygen affinity of hemoglobin. We present the results of 6-ns molecular dynamics simulations performed to gain insights into the dynamical and structural details of these effector-linked tertiary changes. All-atom simulations were carried out on a series of models generated for T- and R-state HbA, and for 2,3-diphosphoglycerate-bound models. Cross-correlation analyses identify both intra- and intersubunit correlated motions that are perturbed by the presence of the effector. Principal components analysis was used to decompose the covariance matrix extracted from the simulations and reconstruct the trajectories along the principal coordinates representative of functionally important collective motions. It is found that HbA in both quaternary states exists as ensembles of tertiary conformations that introduce dynamic heterogeneity in the protein. 2,3-Diphosphoglycerate induces significant perturbations in the fluctuations of both HbA states that translate into the protein visiting different tertiary conformations within each quaternary state. The analysis reveals that the presence of the effector affects the most important components of HbA motions and that heterotropic effectors modify the overall dynamics of the quaternary equilibrium via tertiary changes occurring in regions where conserved functionally significant residues are located, namely in the loop regions between helices C and E, E and F, and F and G, and in concerted helix motions. The changes are not apparent when comparing the available x-ray crystal structures in the presence and absence of effector, but are striking when comparing the respective dynamic tertiary conformations of the R and T tetramers.

Structural characterization of ferric hemoglobins from three antarctic fish species of the suborder notothenioidei

Spontaneous autoxidation of tetrameric Hbs leads to the formation of Fe (III) forms, whose physiological role is not fully understood. Here we report structural characterization by EPR of the oxidized states of tetrameric Hbs isolated from the Antarctic fish species Trematomus bernacchii, Trematomus newnesi, and Gymnodraco acuticeps, as well as the x-ray crystal structure of oxidized Trematomus bernacchii Hb, redetermined at high resolution. The oxidation of these Hbs leads to formation of states that were not usually detected in previous analyses of tetrameric Hbs. In addition to the commonly found aquo-met and hydroxy-met species, EPR analyses show that two distinct hemichromes coexist at physiological pH, referred to as hemichromes I and II, respectively. Together with the high-resolution crystal structure (1.5 A) of T. bernacchii and a survey of data available for other heme proteins, hemichrome I was assigned by x-ray crystallography and by EPR as a bis-His complex with a distorted geometry, whereas hemichrome II is a less constrained (cytochrome b5-like) bis-His complex. In four of the five Antartic fish Hbs examined, hemichrome I is the major form. EPR shows that for HbCTn, the amount of hemichrome I is substantially reduced. In addition, the concomitant presence of a penta-coordinated high-spin Fe (III) species, to our knowledge never reported before for a wild-type tetrameric Hb, was detected. A molecular modeling investigation demonstrates that the presence of the bulkier Ile in position 67beta in HbCTn in place of Val as in the other four Hbs impairs the formation of hemichrome I, thus favoring the formation of the ferric penta-coordinated species. Altogether the data show that ferric states commonly associated with monomeric and dimeric Hbs are also found in tetrameric Hbs.

A comparative NMR study of the polypeptide backbone dynamics of hemoglobin in the deoxy and carbonmonoxy forms

Model-free-based NMR dynamics studies have been undertaken for polypeptide backbone amide N-H bond vectors for both the deoxy and carbonmonoxy forms of chain-specific, isotopically (15N and 2H) labeled tetrameric hemoglobin (Hb) using 15N-relaxation parameters [longitudinal relaxation rate (R1), transverse relaxation rate (R2), and heteronuclear nuclear Overhauser effect (NOE)] measured at two temperatures (29 and 34 degrees C) and two magnetic field strengths (11.7 and 14.1 T). In both deoxy and carbonmonoxy forms of human normal adult hemoglobin (Hb A), the amide N-H bonds of most amino acid residues are rigid on the fast time scale (nanosecond to picosecond), except for the loop regions and certain helix-helix connections. Although rigid in deoxy-Hb A, beta146His has been found to be free from restriction of its backbone motions in the CO form, presumably due to the rupture of its hydrogen bond/salt bridge network. We now have direct dynamics evidence for this structural transition of Hb in solution. While remarkably flexible in the deoxy state, alpha31Arg and beta123Thr, neighbors in the intradimer (alpha1beta1) interface, exhibit stiffening upon CO binding. These findings imply a role for alpha31Arg and beta123Thr in the intradimer communication but contradict the results from X-ray crystallography. We have also found that there is considerable flexibility in the intradimer (alpha1beta1) interface (i.e., B, G, and H helices and the GH corner) and possible involvement of several amino acid residues (e.g., alpha31Arg, beta3Leu, beta41Phe, beta123Thr, and beta146His) in the allosteric pathway. Several amino acid residues at the intradimer interfaces, such as beta109Val, appear to be involved in possible conformational exchange processes. The dynamic picture derived from the present study provides new insights into the traditional description of the stereochemical mechanism for the cooperative oxygenation of Hb A based on X-ray crystallographic results.

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.

Modeling heme proteins using atomistic simulations

Heme proteins are found in all living organisms, and perform a wide variety of tasks ranging from electron transport, to the oxidation of organic compounds, to the sensing and transport of small molecules. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of several relevant issues in heme proteins chemistry: (i) conformational analysis, ligand migration, and solvation effects studied using classical molecular dynamics simulations; (ii) electronic structure and spin state energetics of the active sites explored using quantum-mechanics (QM) methods; (iii) the interaction of heme proteins with small ligands studied through hybrid quantum mechanics-molecular mechanics (QM-MM) techniques; (iv) and finally chemical reactivity and catalysis tackled by a combination of quantum and classical tools.

Nature of the FeO2 bonding in myoglobin and hemoglobin: A new molecular paradigm

The iron(II)-dioxygen bond in myoglobin and hemoglobin is a subject of wide interest. Studies range from examinations of physical-chemical properties dependent on its electronic structure, to investigations of the stability as a function of oxygen supply. Among these, stability properties are of particular importance in vivo. Like all known dioxygen carriers synthesized so far with transition metals, the oxygenated forms of myoglobin and hemoglobin are known to be oxidized easily to their ferric met-forms, which cannot bind molecular oxygen and are therefore physiologically inactive. The mechanistic details of this autoxidation reaction, which are of clinical, as well as of physical-chemical, interest, have long been investigated by a number of authors, but a full understanding of the heme oxidation has not been reached so far. Recent kinetic and thermodynamic studies of the stability of oxymyoglobin (MbO2) and oxyhemoglobin (HbO2) have revealed new features in the FeO2 bonding. In vivo, the iron center is always subject to a nucleophilic attack of the water molecule or hydroxyl ion, which can enter the heme pocket from the surrounding solvent and thereby irreversibly displace the bound dioxygen from MbO2 or HbO2 in the form of O2- so that the iron is converted to the ferric met-form. Since the autoxidation reaction of MbO2 or HbO2 proceeds through a nucleophilic displacement following one-electron transfer from iron(II) to the bound O2, this reaction may be viewed as a meeting point of the stabilization and the activation of molecular oxygen performed by hemoproteins. Along with these lines of evidence, we finally discuss the stability property of human HbO2 and provide with the most recent state of hemoglobin research. The HbA molecule contains two types of alphabeta contacts and seems to differentiate them quite properly for its functional properties. The alpha1beta2 or alpha2beta1 contact is associated with the cooperative oxygen binding, whereas the alpha1beta1 or alpha2beta2 contact is used for controlling the stability of the bound O2. We can thus form a unified picture for hemoglobin function by closely integrating the cooperative and the stable binding of molecular oxygen with iron(II) in aqueous solvent. These new views on the nature of FeO2 bonding and the possible role of globin moiety in stabilizing MbO2 and HbO2 are of primary importance, not only for a full understanding of various hemoprotein reactions with O2, but also for planning new molecular designs for synthetic oxygen carriers which may be able to function in aqueous solvent and at physiological temperature.

Semihemoglobins, high oxygen affinity dimeric forms of human hemoglobin respond efficiently to allosteric effectors without forming tetramers

Significant reduction in oxygen affinity resulting from interactions between heterotropic allosteric effectors and hemoglobin in not only the unligated derivative but also the fully ligated form has been reported (Tsuneshige, A., Park, S. I., and Yonetani, T. (2002) Biophys. Chem. 98, 49-63; Yonetani, T., Park, S. I., Tsuneshige, A., Imai, K., and Kanaori, K. (2002) J. Biol. Chem. 277, 34508-34520). To further investigate this effect in more detail, alpha- and beta-semihemoglobins, namely, alpha(heme)beta(apo) and alpha(apo)beta(heme), respectively, were prepared and characterized with respect to the impact of allosteric effectors on both conformation and ligand binding properties. Semihemoglobins are dimers characterized by a high affinity for oxygen and lack of cooperativity. We found that, compared with stripped conditions, semihemoglobins responded to effectors (inositol hexaphosphate and L35) by decreasing the affinity for oxygen by 60- and 130-fold for alpha- and beta-semihemoglobins, respectively. 1H NMR and sedimentation velocity experiments carried out with their ligated and unligated forms in the absence and presence of effectors revealed that semihemoglobins always remain as single-heme-carrying dimers. Recombination kinetics of their photolyzed CO derivatives showed that effectors did indeed interact with their ligated forms. Measurements of the Fe-His stretching mode show that the semihemoglobins undergo a large ligand binding-induced conformational shift and that both ligand-free and ligand derivatives respond to the presence of effectors. Contradictions to the Monod-Wyman-Changeaux/Perutz allosteric model arise since 1) the modulation of ligand affinity is not achieved in semihemoglobins by the formation of a low affinity T conformation (quaternary effect) but by direct interaction with effectors, 2) effectors do interact significantly with ligated forms of high affinity semihemoglobins, and 3) modulation of the ligand affinity and the cooperativity are not necessarily linked but instead can be separated into two distinct phenomena that can be isolated.

Hemoglobin site-mutants reveal dynamical role of interhelical H-bonds in the allosteric pathway: time-resolved UV resonance Raman evidence for intra-dimer coupling

The dynamical effect of eliminating specific tertiary H-bonds in the hemoglobin (Hb) tetramer has been investigated by site-directed mutagenesis and time-resolved absorption and ultraviolet resonance Raman (UVRR) spectroscopy. The Trp alpha 14...Thr alpha 67 and Trp beta 15...Ser beta 72 H-bonds connect the A and E helices in the alpha and beta chains, and are proposed to break in the earliest protein intermediate (Rdeoxy) following photo-deligation of HbCO, along with a second pair of H-bonds involving tyrosine residues. Mutation of the acceptor residues Thr alpha 67 and Ser beta 72 to Val and Ala eliminates the A-E H-bonds, but has been shown to have no significant effect on ligand-binding affinity or cooperativity, or on spectroscopic markers of the T-state quaternary interactions. However, the mutations have profound and unexpected effects on the character of the Rdeoxy intermediate, and on the dynamics of the subsequent steps leading to the T state. Formation of the initial quaternary contact (RT intermediate) is accelerated, by an order of magnitude, but the locking-in of the T state is delayed by a factor of 2. These rate effects are essentially the same for either mutation, or for the double mutation, suggesting that the alpha beta dimer behaves as a mechanically coupled dynamical unit. Further evidence for intra-dimer coupling is provided by the Rdeoxy UVRR spectrum, in which either or both mutations eliminate the tyrosine difference intensity, although only tryptophan H-bonds are directly affected. A possible mechanism for mechanical coupling is outlined, involving transmission of forces through the alpha(1)beta(1) (and alpha(2)beta(2)) interface. The present observations establish that quaternary motions can occur on the approximately 100 ns time-scale. They show also that a full complement of interhelical H-bonds actually slows the initial quaternary motion in Hb, but accelerates the locking in of the T-contacts.

Human haemoglobin: a new paradigm for oxygen binding involving two types of alphabeta contacts

This review summarizes the most recent state of haemoglobin (Hb) research based on the literature and our own results. In particular, an attempt is made to form a unified picture for haemoglobin function by reconciling the cooperative oxygen binding with the stabilization of the bound dioxygen in aqueous solvent. The HbA molecule contains two types of alphabeta contacts. One type is the alpha1beta2 or alpha2beta1 contacts, called sliding contacts, and these are strongly associated with the cooperative binding of O2 to the alpha2beta2 tetramer. The other type is the alpha1beta1 or alpha2beta2 contacts, called packing contacts, but whose role in Hb function was not clear until quite recently. However, detailed pH-dependence studies of the autoxidation rate of HbO2 have revealed that the alpha1beta1 and alpha2beta2 interfaces are used for controlling the stability of the bound O2. When the alpha1beta1 or alpha2beta2 contact is formed, the beta chain is subjected to a conformational constraint which causes the distal (E7) histidine to be tilted slightly away from the bound dioxygen, preventing the proton-catalysed nucleophilic displacement of O2- from the FeO2 by an entering water molecule. This is one of the most characteristic features of HbO2 stability. Finally we discuss the role of the alpha1beta1 or alpha2beta2 contacts by providing some examples of unstable haemoglobin mutants. These pathological mutations are found mostly on the beta chain, especially in the alpha1beta1 contact regions. In this way, HbA seems to differentiate two types of alphabeta contacts for its functional properties.

A pH-dependent aquomet-to-hemichrome transition in crystalline horse methemoglobin

In 1947, Perutz and co-workers reported that crystalline horse methemoglobin undergoes a large lattice transition as the pH is decreased from 7.1 to 5.4. We have determined the pH 7.1 and 5.4 crystal structures of horse methemoglobin at 1.6 and 2.1 A resolution, respectively, and find that this lattice transition involves a 23 A translation of adjacent hemoglobin tetramers as well as changes in alpha heme ligation and the tertiary structure of the alpha subunits. Specifically, when the pH is lowered from 7.1 to 5.4, the Fe(3+) alpha heme groups (but not the beta heme groups) are converted from the aquomet form, in which the proximal histidine [His87(F8)alpha] and a water molecule are the axial heme ligands, to the hemichrome (bishistidine) form, in which the proximal histidine and the distal histidine [His58(E7)alpha] are the axial heme ligands. Hemichrome formation is coupled to a large tertiary structure transition in the eight-residue segment Pro44(CD2)alpha-Gly51(D7)alpha that converts from an extended loop structure at pH 7.1 to a pi-like helix at pH 5.4. The formation of the pi helix forces Phe46(CD4)alpha out of the alpha heme pocket and into the interface between adjacent hemoglobin tetramers where it participates in crystal lattice contacts unique to the pH 5.4 structure. In addition, the transition from aquomet alpha subunits to bishistidine alpha subunits is accompanied by an approximately 1.2 A movement of the alpha heme groups to a more solvent-exposed position as well as the creation of a solvent channel from the interior of the alpha heme pocket to the outside of the tetramer. These changes and the extensive rearrangement of the crystal lattice structure allow the alpha heme group of one tetramer to make direct contact with an alpha heme group on an adjacent tetramer. These results suggest possible functional roles for hemichrome formation in vivo.

Cooperativity in Scapharca dimeric hemoglobin: simulation of binding intermediates and elucidation of the role of interfacial water

Cooperative binding of ligands to proteins can serve to increase their efficiency and to regulate their activity. Thus, understanding of the mechanism of cooperativity is one of the central concerns of molecular biology. For the tetrameric human hemoglobin (HbA), the cooperative mechanism involves a reasonably well understood combination of tertiary and quaternary changes that occur during the binding process. The dimeric hemoglobin of Scapharca (HbI), which is composed of subunits with the same fold as in HbA, is also highly cooperative but the structural changes on ligand binding are small. A re-orientation of Phe97 in the binding pocket and changes in the number of interfacial water molecules have been implicated in the cooperative mechanism. To explore the role of these factors, we have investigated models of partially liganded intermediate states of HbI with molecular dynamics simulation methods. Since, unlike HbA, no structures for intermediates are available, they were constructed by combining subunits from the unliganded and liganded dimers. Two structurally distinct intermediates were examined, and it was shown that the transition between the two intermediates is directly coupled to the number of interfacial water molecules. Further, it was found that there is a well-defined water channel that connects the interface between the subunits to bulk water. The bottleneck (gate) of the channel, which can be open or closed, is made of hydrophilic residues. The implication of the present results for the cooperative mechanism of HbI is discussed.

Single residue modification of only one dimer within the hemoglobin tetramer reveals autonomous dimer function

The mechanism of cooperativity in the human hemoglobin tetramer (a dimer of alpha beta dimers) has historically been modeled as a simple two-state system in which a low-affinity structural form (T) switches, on ligation, to a high-affinity form (R), yielding a net loss of hydrogen bonds and salt bridges in the dimer-dimer interface. Modifications that weaken these cross-dimer contacts destabilize the quaternary T tetramer, leading to decreased cooperativity and enhanced ligand affinity, as demonstrated in many studies on symmetric double modifications, i.e., a residue site modified in both alpha- or both beta-subunits. In this work, hybrid tetramers have been prepared with only one modified residue, yielding molecules composed of a wild-type dimer and a modified dimer. It is observed that the cooperative free energy of ligation to the modified dimer is perturbed to the same extent whether in the hybrid tetramer or in the doubly modified tetramer. The cooperative free energy of ligation to the wild-type dimer is unperturbed, even in the hybrid tetramer, and despite the overall destabilization of the T tetramer by the modification. This asymmetric response by the two dimers within the same tetramer shows that loss of dimer-dimer contacts is not communicated across the dimer-dimer interface, but is transmitted through the dimer that bears the modified residue. These observations are interpreted in terms of a previously proposed dimer-based model of cooperativity with an additional quaternary (T/R) component.

New insights into the allosteric mechanism of human hemoglobin from molecular dynamics simulations

It is still difficult to obtain a precise structural description of the transition between the deoxy T-state and oxy R-state conformations of human hemoglobin, despite a large number of experimental studies. We used molecular dynamics with the Path Exploration with Distance Constraints (PEDC) method to provide new insights into the allosteric mechanism at the atomic level, by simulating the T-to-R transition. The T-state molecule in the absence of ligands was seen to have a natural propensity for dimer rotation, which nevertheless would be hampered by steric hindrance in the "joint" region. The binding of a ligand to the alpha subunit would prevent such hindrance due to the coupling between this region and the alpha proximal histidine, and thus facilitate completion of the dimer rotation. Near the end of this quaternary transition, the "switch" region adopts the R conformation, resulting in a shift of the beta proximal histidine. This leads to a sliding of the beta-heme, the effect of which is to open the beta-heme's distal side, increasing the accessibility of the Fe atom and thereby the affinity of the protein. Our simulations are globally consistent with the Perutz strereochemical mechanism.

Effects of amino acid substitutions at beta 131 on the structure and properties of hemoglobin: evidence for communication between alpha 1 beta 1- and alpha 1 beta 2-subunit interfaces

Substitutions of Asn, Glu, and Leu for Gln at the beta131 position of the hemoglobin molecule result in recombinant hemoglobins (rHbs) with moderately lowered oxygen affinity and high cooperativity compared to human normal adult hemoglobin (Hb A). The mutation site affects the hydrogen bonds present at the alpha(1)beta(1)-subunit interface between alpha103His and beta131Gln as well as that between alpha122His and beta35Tyr. NMR spectroscopy shows that the hydrogen bonds are indeed perturbed; in the case of rHb (beta131Gln --> Asn) and rHb (beta131Gln --> Leu), the perturbations are propagated to the other alpha(1)beta(1)-interface H-bond involving alpha122His and beta35Tyr. Proton exchange measurements also detect faster exchange rates for both alpha(1)beta(1)-interface histidine side chains of the mutant rHbs in 0.1 M sodium phosphate buffer at pH 7.0 than for those of Hb A under the same conditions. In addition, the same measurements in 0.1 M Tris buffer at pH 7.0 show a much slower exchange rate for mutant rHbs and Hb A. One of the mutants, rHb (beta131Gln --> Asn), shows the conformational exchange of its interface histidines, and exchange rate measurements have been attempted. We have also conducted studies on the reactivity of the SH group of beta93Cys (a residue located in the region of the alpha(1)beta(2)-subunit interface) toward p-mercuribenzoate, and our results show that low-oxygen-affinity rHbs have a more reactive beta93Cys than Hb A in the CO form. Our results indicate that there is communication between the alpha(1)beta(1)- and alpha(1)beta(2)-subunit interfaces, and a possible communication pathway for the cooperative oxygenation of Hb A that allows the alpha(1)beta(1)-subunit interface to modulate the functional properties in conjunction with the alpha(1)beta(2) interface is proposed.

Microsecond to minute dynamics revealed by EX1-type hydrogen exchange at nearly every backbone hydrogen bond in a native protein

A previous comprehensive analysis of the pH dependence of native-state amide hydrogen (NH) exchange in turkey ovomucoid third domain (OMTKY3) yielded apparent opening and closing rate constants (k(op) and k(cl)) at 14 NH groups involved in global conformational changes. This analysis has been extended to 18 additional slowly exchanging NH groups. Quench-flow experiments were performed to monitor NH exchange in native OMTKY3 from neutral to very alkaline pH ( approximately 12) conditions. Above pH 10 the mechanism of exchange switched from one governed by a rapid equilibrium preceding the chemistry of exchange (i.e. EX2 exchange), to one where exchange was limited by the rate of opening (i.e. EX1 exchange). Kinetics of solvent exposure are now known for nearly all backbone NH groups in native OMTKY3, yielding rate constants that span five orders of magnitude, 0.004 to 200 s(-1).