Global allostery model of hemoglobin. Modulation of O(2) affinity, cooperativity, and Bohr effect by heterotropic allosteric effectors

Structures of haemoglobin from woolly mammoth in liganded and unliganded states

The haemoglobin (Hb) of the extinct woolly mammoth has been recreated using recombinant genes expressed in Escherichia coli. The globin gene sequences were previously determined using DNA recovered from frozen cadavers. Although highly similar to the Hb of existing elephants, the woolly mammoth protein shows rather different responses to chloride ions and temperature. In particular, the heat of oxygenation is found to be much lower in mammoth Hb, which appears to be an adaptation to the harsh high-latitude climates of the Pleistocene Ice Ages and has been linked to heightened sensitivity of the mammoth protein to protons, chloride ions and organic phosphates relative to that of Asian elephants. To elucidate the structural basis for the altered homotropic and heterotropic effects, the crystal structures of mammoth Hb have been determined in the deoxy, carbonmonoxy and aquo-met forms. These models, which are the first structures of Hb from an extinct species, show many features reminiscent of human Hb, but underline how the delicate control of oxygen affinity relies on much more than simple overall quaternary-structure changes.

Local structure and dynamics of hemeproteins by X-ray absorption near edge structure spectroscopy

X-ray absorption near edge structure (XANES) spectroscopy is a synchrotron radiation technique sensitive to the local structure and dynamics around the metal site of a heme containing protein. Advances in detection techniques and theoretical/computational platforms in the last 15 years allowed the use of XANES as a quantitative probe of the key structural determinants driving functional changes, both in a concerted way with protein crystallography and EXAFS (extended X-ray absorption fine structure), or as a stand-alone method to apply in the crystal state as well as in solution. Moreover, the local dynamics of the heme site has been deeply investigated, on one hand, coupling XANES to classical photolysis experiments at cryogenic temperatures; on the other hand, the intrinsic property of the synchrotron radiation to induce radiolysis events, has been exploited to investigate specific cryotrapped intermediates, using X-rays both as a pump and a probe. Insights on the XANES method and some specific examples are presented to illustrate these topics.

Heme reactivity is uncoupled from quaternary structure in gel-encapsulated hemoglobin: a resonance Raman spectroscopic study

Encapsulation of hemoglobin (Hb) in silica gel preserves structure and function but greatly slows protein motion, thereby providing access to intermediates along the allosteric pathway that are inaccessible in solution. Resonance Raman (RR) spectroscopy with visible and ultraviolet laser excitation provides probes of heme reactivity and of key tertiary and quaternary contacts. These probes were monitored in gels after deoxygenation of oxyHb and after CO binding to deoxyHb, which initiate conformational change in the R-T and T-R directions, respectively. The spectra establish that quaternary structure change in the gel takes a week or more but that the evolution of heme reactivity, as monitored by the Fe-histidine stretching vibration, ν(FeHis), is completed within two days, and is therefore uncoupled from the quaternary structure. Within each quaternary structure, the evolving ν(FeHis) frequencies span the full range of values between those previously associated with the high- and low-affinity end states, R and T. This result supports the tertiary two-state (TTS) model, in which the Hb subunits can adopt high- and low-affinity tertiary structures, r and t, within each quaternary state. The spectra also reveal different tertiary pathways, involving the breaking and reformation of E and F interhelical contacts in the R-T direction but not the T-R direction. In the latter, tertiary motions are restricted by the T quaternary contacts.

Steric factors moderate conformational fluidity and contribute to the high proton sensitivity of Root effect hemoglobins

The structural basis of the extreme pH dependence of oxygen binding to Root effect Hbs is a long-standing puzzle in the field of protein chemistry. A previously unappreciated role of steric factors in the Root effect was revealed by a comparison of pH effects on oxygenation and oxidation processes in human Hb relative to Spot (Leiostomus xanthurus) and Carp (Cyprinodon carpio) Hbs. The Root effect confers five-fold increased pH sensitivity to oxygenation of Spot and Carp Hbs relative to Hb A(0) in the absence of anionic effectors, and even larger relative elevations of pH sensitivity of oxygenation in the presence of 0.2M phosphate. Remarkably, the Root effect was not evident in the oxidation of the Root effect Hbs. This finding rules out pH-dependent alterations in the thermodynamic properties of the heme iron, measured in the anaerobic oxidation reaction, as the basis of the Root effect. The alternative explanation supported by these results is that the elevated pH sensitivity of oxygenation of Root effect Hbs is attributable to globin-dependent steric effects that alter oxygen affinity by constraining conformational fluidity, but which have little influence on electron exchange via the heme edge. This elegant mode of allosteric control can regulate oxygen affinity within a given quaternary state, in addition to modifying the T-R equilibrium. Evolution of Hb sequences that result in proton-linked steric barriers to heme oxygenation could provide a general mechanism to account for the appearance of the Root effect in the structurally diverse Hbs of many species.

T-quaternary structure of oxy human adult hemoglobin in the presence of two allosteric effectors, L35 and IHP

The cooperative O(2)-binding of hemoglobin (Hb) have been assumed to correlate to change in the quaternary structures of Hb: T(deoxy)- and R(oxy)-quaternary structures, having low and high O(2)-affinities, respectively. Heterotropic allosteric effectors have been shown to interact not only with deoxy- but also oxy-Hbs causing significant reduction in their O(2)-affinities and the modulation of cooperativity. In the presence of two potent effectors, L35 and inositol hexaphosphate (IHP) at pH 6.6, Hb exhibits extremely low O(2)-affinities (K(T)=0.0085mmHg(-1) and K(R)=0.011mmHg(-1)) and thus a very low cooperativity (K(R)/K(T)=1.3 and L(0)=2.4). (1)H-NMR spectra of human adult Hb with these two effectors were examined in order to determine the quaternary state of Hb in solution and to clarify the correlation between the O(2)-affinities and the structural change of Hb caused by the heterotropic effectors. At pH 6.9, (1)H-NMR spectrum of deoxy-Hb in the presence of L35 and IHP showed a marker of the T-quaternary structure (the T-marker) at 14ppm, originated from inter- dimeric α(1)β(2)- (or α(2)β(1)-) hydrogen-bonds, and hyperfine-shifted (hfs) signals around 15-25ppm, caused by high-spin heme-Fe(II)s. Upon addition of O(2), the hfs signals disappeared, reflecting that the heme-Fe(II)s are ligated with O(2), but the T-marker signals still remained, although slightly shifted and broadened, under the partial pressure of O(2) (P(O2)) of 760mmHg. These NMR results accompanying with visible absorption spectroscopy and visible resonance Raman spectroscopy reveal that oxy-Hb in the presence of L35 and IHP below pH 7 takes the ligated T-quaternary structure under the P(O2) of 760mmHg. The L35-concentration dependence of the T-marker in the presence of IHP indicates that there are more than one kind of L35-binding sites in the ligated T-quaternary structure. The stronger binding sites are probably intra-dimeric binding sites between α(1)G- and β(1)G-helices, and the other weaker binding site causes the R→T transition without release of O(2). The fluctuation of the tertiary structure of Hb seems to be caused by both the structural perturbation of α(1)β(1) (or α(2)β(2)) intra-dimeric interface, where the stronger L35-binding sites exist, and by the IHP-binding to the α(1)α(2)- (or β(1)β(2)-) cavity. The tertiary structural fluctuation induced by the allosteric effectors may contribute to the significant reduction of the O(2)-affinity of oxy-Hb, which little depends on the quaternary structures. Therefore, the widely held assumptions of the structure-function correlation of Hb - [the deoxy-state]=[the T-quaternary structure]=[the low O(2)-affinity state] and [the oxy-state]=[the R-quaternary structure]=[the high O(2)-affinity state] and the O(2)-affiny of Hb being regulated by the T/R-quaternary structural transition - are no longer sustainable. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Hemoglobin allostery: variations on the theme

The two-state allosteric model of Monod, Wyman and Changeux (1965) offers a simple and elegant, yet very powerful and comprehensive, description of the functional behavior of hemoglobin. Although the extensive body of structural and functional information available is by-and-large consistent with this conceptual framework, some discrepancies between theory and experiment have been extensively discussed and considered to demand modifications of the original hypothesis. More recently the role of tertiary structural changes has been re-analyzed leading to extended kinetic models or indicating that powerful heterotropic effectors may be of paramount importance in controlling the function of human hemoglobin. The aim of this review is to analyze, and possibly reconcile, some discrepancies. We always felt that by looking at hemoglobins other than human HbA, the relative role of tertiary and quaternary allosteric effects may be better understood. The model systems illustrated below are the different hemoglobins from trout's blood, since they are characterized by the most striking variability of heterotropic effects, ranging from totally absent to very extreme with dominant contributions of tertiary effects. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Hemoglobin-ligand binding: understanding Hb function and allostery on atomic level

The major physiological function of hemoglobin (Hb) is to bind oxygen in the lungs and deliver it to the tissues. This function is regulated and/or made efficient by endogenous heterotropic effectors. A number of synthetic molecules also bind to Hb to alter its allosteric activity. Our purpose is to review the current state of Hb structure and function that involves ensemble of tense and relaxed hemoglobin states and the dynamic equilibrium of the multistate due to the binding of endogenous heterotropic or synthetic allosteric effectors. The review also discusses the atomic interactions of synthetic ligands with the function or altered allosteric function of Hb that could be potentially harnessed for the treatment of diseases. This article is part of a Special Issue entitled: Protein Structure and Function in the Crystalline State.

An investigation of the distal histidyl hydrogen bonds in oxyhemoglobin: effects of temperature, pH, and inositol hexaphosphate

On the basis of X-ray crystal structures and electron paramagnetic resonance (EPR) measurements, it has been inferred that the O(2) binding to hemoglobin is stabilized by the hydrogen bonds between the oxygen ligands and the distal histidines. Our previous study by multinuclear nuclear magnetic resonance (NMR) spectroscopy has provided the first direct evidence of such H-bonds in human normal adult oxyhemoglobin (HbO(2) A) in solution. Here, the NMR spectra of uniformly (15)N-labeled recombinant human Hb A (rHb A) and five mutant rHbs in the oxy form have been studied under various experimental conditions of pH and temperature and also in the presence of an organic phosphate, inositol hexaphosphate (IHP). We have found significant effects of pH and temperature on the strength of the H-bond markers, i.e., the cross-peaks for the side chains of the two distal histidyl residues, α58His and β63His, which form H-bonds with the O(2) ligands. At lower pH and/or higher temperature, the side chains of the distal histidines appear to be more mobile, and the exchange with water molecules in the distal heme pockets is faster. These changes in the stability of the H-bonds with pH and temperature are consistent with the changes in the O(2) affinity of Hb as a function of pH and temperature and are clearly illustrated by our NMR experiments. Our NMR results have also confirmed that this H-bond in the β-chain is weaker than that in the α-chain and is more sensitive to changes in pH and temperature. IHP has only a minor effect on these H-bond markers compared to the effects of pH and temperature. These H-bonds are sensitive to mutations in the distal heme pockets but not affected directly by the mutations in the quaternary interfaces, i.e., α(1)β(1) and/or α(1)β(2) subunit interface. These findings provide new insights regarding the roles of temperature, hydrogen ion, and organic phosphate in modulating the structure and function of hemoglobin in solution.

Protonation states of histidine and other key residues in deoxy normal human adult hemoglobin by neutron protein crystallography

The protonation states of the histidine residues key to the function of deoxy (T-state) human hemoglobin have been investigated using neutron protein crystallography. These residues can reversibly bind protons, thereby regulating the oxygen affinity of hemoglobin. By examining the OMIT F(o)-F(c) and 2F(o)-F(c) neutron scattering maps, the protonation states of 35 of the 38 His residues were directly determined. The remaining three residues were found to be disordered. Surprisingly, seven pairs of His residues from equivalent α or β chains, αHis20, αHis50, αHis58, αHis89, βHis63, βHis143 and βHis146, have different protonation states. The protonation of distal His residues in the α(1)β(1) heterodimer and the protonation of αHis103 in both subunits demonstrates that these residues may participate in buffering hydrogen ions and may influence the oxygen binding. The observed protonation states of His residues are compared with their ΔpK(a) between the deoxy and oxy states. Examination of inter-subunit interfaces provided evidence for interactions that are essential for the stability of the deoxy tertiary structure.

Exploring and exploiting allostery: Models, evolution, and drug targeting

The concept of allostery was elaborated almost 50years ago by Monod and coworkers to provide a framework for interpreting experimental studies on the regulation of protein function. In essence, binding of a ligand at an allosteric site affects the function at a distant site exploiting protein flexibility and reshaping protein energy landscape. Both monomeric and oligomeric proteins can be allosteric. In the past decades, the behavior of allosteric systems has been analyzed in many investigations while general theoretical models and variations thereof have been steadily proposed to interpret the experimental data. Allostery has been established as a fundamental mechanism of regulation in all organisms, governing a variety of processes that range from metabolic control to receptor function and from ligand transport to cell motility. A number of studies have shed light on how evolutionary pressures have favored and molded the development of allosteric features in specific macromolecular systems. The widespread occurrence of allostery has been recently exploited for the development and design of allosteric drugs that bind to either physiological or non-physiological allosteric sites leading to gain of function or loss of function. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Direct determination of protonation states of histidine residues in a 2 A neutron structure of deoxy-human normal adult hemoglobin and implications for the Bohr effect

We have investigated the protonation states of histidine residues (potential Bohr groups) in the deoxy form (T state) of human hemoglobin by direct determination of hydrogen (deuterium) positions with the neutron protein crystallography technique. The reversible binding of protons is key to the allosteric regulation of human hemoglobin. The protonation states of 35 of the 38 His residues were directly determined from neutron scattering omit maps, with 3 of the remaining residues being disordered. Protonation states of 5 equivalent His residues--alpha His20, alpha His50, alpha His89, beta His143, and beta His146--differ between the symmetry-related globin subunits. The distal His residues, alpha His58 and beta His63, are protonated in the alpha 1 beta 1 heterodimer and are neutral in alpha 2 beta 2. Buried residue alpha His103 is found to be protonated in both subunits. These distal and buried residues have the potential to act as Bohr groups. The observed protonation states of His residues are compared to changes in their pK(a) values during the transition from the T to the R state and the results provide some new insights into our understanding of the molecular mechanism of the Bohr effect.

Oximetry with the NMR signals of hemoglobin Val E11 and Tyr C7

The NMR visibility of the signals from erythrocyte hemoglobin (Hb) presents an opportunity to assess the vascular PO₂ (partial pressure of oxygen) in vivo to gather insight into the regulation of O₂ transport, especially in contracting muscle tissue. Some concerns, however, have arisen about the validity of using the Val E11 signal as an indicator of PO₂, since its intensity depends on tertiary structural changes, in contrast to the quaternary structure changes associated with relaxed (R) and tense (T) transition during O₂ binding. We have examined the Val E11 and Tyr C7 signal intensity as a function of Hb saturation by developing an oximetry system, which permits the comparative analysis of the NMR and spectrophotometric measurements. The spectrophotometric assay defines the Hb saturation level at a given PO₂ and yields standard oxygen-binding curves. Under defined PO₂ and Hb saturation values, the NMR measurements have determined that the Val E11 signal, as well as the Tyr C7 signal, tracks closely Hb saturation and can therefore serve as a vascular oxygen biomarker.

Computing phenomenologic Adair-Klotz constants from microscopic MWC parameters

BACKGROUND

Modellers using the MWC allosteric framework have often found it difficult to validate their models. Indeed many experiments are not conducted with the notion of alternative conformations in mind and therefore do not (or cannot) measure relevant microscopic constant and parameters. Instead, experimentalists widely use the Adair-Klotz approach in order to describe their experimental data.

RESULTS

We propose a way of computing apparent Adair-Klotz constants from microscopic association constants and allosteric parameters of a generalised concerted model with two different states (R and T), with an arbitrary number of non-equivalent ligand binding sites. We apply this framework to compute Adair-Klotz constants from existing models of calmodulin and hemoglobin, two extreme cases of the general framework.

CONCLUSION

The validation of computational models requires methods to relate model parameters to experimentally observable quantities. We provide such a method for comparing generalised MWC allosteric models to experimentally determined Adair-Klotz constants.

Transition of haemoglobin between two tertiary conformations: Inositol hexakisphosphate increases the transition constant and the affinity of sheep haemoglobin for 5,5'-dithiobis(2-nitrobenzoate)

The equilibrium constant (K(equ)) for the reaction of 5,5'-dithiobis(2-nitrobenzoate) - DTNB - with the CysF9[93]beta sulphydryl group of the haemoglobins of the sheep decreases by about two orders of magnitude between pH approximately 5.6 and 9.2: from a mean of 7.2+/-1 to a mean of 0.044+/-0.01. Calculations from the pH dependence of K(equ) show that in the r-t tertiary conformational transition of haemoglobin the t isomer population is 50.7 and 61.8% for the major and minor haemoglobins, respectively. In the presence of inositol hexakisphosphate (inositol-P(6)), K(equ) increases for both haemoglobins by about an order of magnitude through most of the pH range. The t isomer population also increases to 82.1 and 79.6% for the major and minor haemoglobins, respectively. These results indicate that inositol-P(6) increases the affinity of the sulphydryl for DTNB by increasing the population of the t isomer. It is highly probable that a minimum four-state model that includes the r-t transition is required for a full understanding of haemoglobin function.

Structural basis for the heterotropic and homotropic interactions of invertebrate giant hemoglobin

The oxygen binding properties of extracellular giant hemoglobins (Hbs) in some annelids exhibit features significantly different from those of vertebrate tetrameric Hbs. Annelid giant Hbs show cooperative oxygen binding properties in the presence of inorganic cations, while the cooperativities of vertebrate Hbs are enhanced by small organic anions or chloride ions. To elucidate the structural basis for the cation-mediated cooperative mechanisms of these giant Hbs, we determined the crystal structures of Ca2+- and Mg2+-bound Hbs from Oligobrachia mashikoi at 1.6 and 1.7 A resolution, respectively. Both of the metal-bound structures were determined in the oxygenated state. Four Ca2+-binding sites and one Mg2+-binding site were identified in each tetramer subassembly. These cations are considered to stabilize the oxygenated form and increase affinity and cooperativity for oxygen binding, as almost all of the Ca2+ and Mg2+ cations were bound at the interface regions, forming either direct or hydrogen bond-mediated interactions with the neighboring subunits. A comparison of the structures of the oxygenated form and the partially unliganded form provides structural insight into proton-coupled cooperativity (Bohr effect) and ligand-induced transitions. Two histidine residues are assumed to be primarily associated with the Bohr effect. With regard to the ligand-induced cooperativity, a novel quaternary rotation mechanism is proposed to exist at the interface region of the dimer subassembly. Interactions among conserved residues Arg E10, His F3, Gln F7, and Val E11, together with the bending motion of the heme molecules, appear to be essential for quaternary rearrangement.

Oxidative effects of gemfibrozil on anion influx and metabolism in normal and Beta-thalassemic erythrocytes: physiological implications

To further clarify some peculiar molecular mechanisms related to the physiology and pathophysiology of erythrocytes with respect to oxygen binding and release, metabolism and senescence, we investigated the oxidative effects of gemfibrozil in normal and beta-thalassemic red blood cells. Our results showed that the oxidative stress promoted by the drug, through a direct interaction with hemoglobin, may lead to activation of caspase 3, which in turn influences the band 3 anion flux and glucose metabolism. In a comparative context, we also evaluated the effect on band 3 and caspase 3 activation of orthovanadate (a phosphatase inhibitor) and t-butylhydroperoxide (a known oxidant). The results support the hypothesis that gemfibrozil influences band 3 function through several mechanisms of action, centered on oxidative stress, which induces significant alterations of glucose metabolism.

Interfacial and distal-heme pocket mutations exhibit additive effects on the structure and function of hemoglobin

Protein engineering strategies seek to develop a hemoglobin-based oxygen carrier with optimized functional properties, including (i) an appropriate O 2 affinity, (ii) high cooperativity, (iii) limited NO reactivity, and (iv) a diminished rate of auto-oxidation. The mutations alphaL29F, alphaL29W, alphaV96W and betaN108K individually impart some of these traits and in combinations produce hemoglobin molecules with interesting ligand-binding and allosteric properties. Studies of the ligand-binding properties and solution structures of single and multiple mutants have been performed. The aromatic side chains placed in the distal-heme pocket environment affect the intrinsic ligand-binding properties of the mutated subunit itself, beyond what can be explained by allostery, and these changes are accompanied by local structural perturbations. In contrast, hemoglobins with mutations in the alpha 1beta 1 and alpha 1beta 2 interfaces display functional properties of both "R"- and "T"-state tetramers because the equilibrium between quaternary states is altered. These mutations are accompanied by global structural perturbations, suggesting an indirect, allostery-driven cause for their effects. Combinations of the distal-heme pocket and interfacial mutations exhibit additive effects in both structural and functional properties, contribute to our understanding of allostery, and advance protein-engineering methods for manipulating the O 2 binding properties of the hemoglobin molecule.

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.

The model of Monod, Wyman, and Changeux generates two sets of parameters when applied to oxygen binding in hemoglobin

The model of Monod, Wyman and Changeux is applied to binding phenomena where the Mass Law and its expansion are employed. In this communication the model of Monod, Wyman and Changeux (MWC) is applied to analyze the oxygen binding reaction in hemoglobin. The symmetrical structure of the MWC model with its three parameters is such that two sets of these parameters, rather than one, fit experimental data for the binding of oxygen to hemoglobin.

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.

Conformational dependence of hemoglobin reactivity under high viscosity conditions: the role of solvent slaved dynamics

The concept of protein dynamic states is introduced. This concept is based on (i) protein dynamics being organized hierarchically with respect to solvent slaving and (ii) which tier of dynamics is operative over the time window of a given measurement. The protein dynamic state concept is used to analyze the kinetic phases derived from the recombination of carbon monoxide to sol-gel-encapsulated human adult hemoglobin (HbA) and select recombinant mutants. The temperature-dependent measurements are made under very high viscosity conditions obtained by bathing the samples in an excess of glycerol. The results are consistent with a given tier of solvent slaved dynamics becoming operative at a time delay (with respect to the onset of the measurement) that is primarily solvent- and temperature-dependent. However, the functional consequences of the dynamics are protein- and conformation-specific. The kinetic traces from both equilibrium populations and trapped allosteric intermediates show a consistent progression that exposes the role of both conformation and hydration in the control of reactivity. Iron-zinc symmetric hybrid forms of HbA are used to show the dramatic difference between the kinetic patterns for T state alpha and beta subunits. The overall results support a model for allostery in HbA in which the ligand-binding-induced transition from the deoxy T state to the high -affinity R state proceeds through a progression of T state intermediates.

Intermediates detected by visible spectroscopy during the reaction of nitrite with deoxyhemoglobin: the effect of nitrite concentration and diphosphoglycerate

The reaction of nitrite with deoxyhemoglobin (deoxyHb) results in the reduction of nitrite to NO, which binds unreacted deoxyHb forming Fe(II)-nitrosylhemoglobin (Hb(II)NO). The tight binding of NO to deoxyHb is, however, inconsistent with reports implicating this reaction with hypoxic vasodilation. This dilemma is resolved by the demonstration that metastable intermediates are formed in the course of the reaction of nitrite with deoxyHb. The level of intermediates is quantitated by the excess deoxyHb consumed over the concentrations of the final products formed. The dominant intermediate has a spectrum that does not correspond to that of Hb(III)NO formed when NO reacts with methemoglobin (MetHb), but is similar to metHb resulting in the spectroscopic determinations of elevated levels of metHb. It is a delocalized species involving the heme iron, the NO, and perhaps the beta-93 thiol. The putative role for red cell reacted nitrite on vasodilation is associated with reactions involving the intermediate. (1) The intermediate is less stable with a 10-fold excess of nitrite and is not detected with a 100-fold excess of nitrite. This observation is attributed to the reaction of nitrite with the intermediate producing N2O3. (2) The release of NO quantitated by the formation of Hb(II)NO is regulated by changes in the distal heme pocket as shown by the 4.5-fold decrease in the rate constant in the presence of 2,3-diphosphoglycerate. The regulated release of NO or N2O3 as well as the formation of the S-nitroso derivative of hemoglobin, which has also been reported to be formed from the intermediates generated during nitrite reduction, should be associated with any hypoxic vasodilation attributed to the RBC.

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.

Hemocyanin from E. californicum encapsulated in silica gels: Oxygen binding and conformational states

Cooperativity depends on the existence of equilibria among functionally distinct conformational states that are affected by homo and heterotropic effectors. In order to isolate the quaternary conformations of hemocyanin from E. californicum, the 24-meric giant protein was encapsulated in wet, nanoporous silica gels, either in the absence or presence of oxygen. The deoxy- and oxy-hemocyanin gels exhibit a p50 for oxygen of 11 and 2.5 torr, respectively, values in close agreement with those for hemocyanin in solution. The observed Hill coefficients are lower than unity, indicating a conformational heterogeneity within each locked conformational state, a finding in agreement with the assumption that at least four conformational states are required to explain the oxygen binding properties of E. californicum hemocyanin in solution.

Inhibiting HIV-1 integrase by shifting its oligomerization equilibrium

Proteins are involved in various equilibria that play a major role in their activity or regulation. The design of molecules that shift such equilibria is of great therapeutic potential. This fact was demonstrated in the cases of allosteric inhibitors, which shift the equilibrium between active and inactive (R and T) states, and chemical chaperones, which shift folding equilibrium of proteins. Here, we expand these concepts and propose the shifting of oligomerization equilibrium of proteins as a general methodology for drug design. We present a strategy for inhibiting proteins by "shiftides": ligands that specifically bind to an inactive oligomeric state of a disease-related protein and modulate its activity by shifting the oligomerization equilibrium of the protein toward it. We demonstrate the feasibility of our approach for the inhibition of the HIV-1 integrase (IN) protein by using peptides derived from its cellular-binding protein, LEDGF/p75, which specifically inhibit IN activity by a noncompetitive mechanism. The peptides inhibit the DNA-binding of IN by shifting the IN oligomerization equilibrium from the active dimer toward the inactive tetramer, which is unable to catalyze the first integration step of 3' end processing. The LEDGF/p75-derived peptides inhibit the enzymatic activity of IN in vitro and consequently block HIV-1 replication in cells because of the lack of integration. These peptides are promising anti-HIV lead compounds that modulate oligomerization of IN via a previously uncharacterized mechanism, which bears advantages over the conventional interface dimerization inhibitors.

Method for Determination of Association and Dissociation Rate Constants of Reversible Bimolecular Reactions by Isothermal Titration Calorimeters

The rate law equation for reversible bimolecular reactions, which are describable by association and dissociation rate constants (k1 and k-1), is not solvable to a plain formula under stoichiometric reaction conditions. Therefore, it is a general technique to observe such reactions under pseudo first-order conditions, which make the reactions a single-exponential process, and enable us to determine k1 and k-1 without any complicated iterative computations needed to analyze the same reactions under stoichiometric reaction conditions. However, the accelerated reaction rates under pseudo first-order conditions are not always favorable to the physicochemical tools employing a slow or medium response time, such as thermal analysis instruments. In this study, we have developed a simple non-iterative analytical method to determine k1 and k-1 of reversible bimolecular reactions under stoichiometric conditions on the basis of experimental data of isothermal titration calorimetry (ITC), which is generally used to determine thermodynamic parameters rather than kinetic constants. Our method is principally based on the general principle of chemical bindings caused along with the titration processes, that is, the chemical relaxation kinetics, which had been hitherto considered in the analysis on the ITC data.

A 45-ns molecular dynamics simulation of hemoglobin in water by vectorizing and parallelizing COSMOS90 on the earth simulator: Dynamics of tertiary and quaternary structures

Molecular dynamics (MD) simulations of human adult hemoglobin (HbA) were carried out for 45 ns in water with all degrees of freedom including bond stretching and without any artificial constraints. To perform such large-scale simulations, one of the authors (M.S.) accelerated his own software COSMOS90 on the Earth Simulator by vectorization and parallelization. The dynamical features of HbA were investigated by evaluating root-mean-square deviations from the initial X-ray structure (an oxy T-state hemoglobin with PDB code: 1GZX) and root-mean-square fluctuations around the average structure from the simulation trajectories. The four subunits (alpha(1), alpha(2), beta(1), and beta(2)) of HbA maintained structures close to their respective X-ray structures during the simulations even though no constraints were applied to HbA in the simulations. Dimers alpha(1)beta(1) and alpha(2)beta(2) also maintained structures close to their respective X-ray structures while they moved relative to each other like two stacks of dumbbells. The distance between the two dimers (alpha(1)beta(1) and alpha(2)beta(2)) increased by 2 A (7.4%) in the initial 15 ns and stably fluctuated at the distance with the standard deviation 0.2 A. The relative orientation of the two dimers fluctuated between the initial X-ray angle -100 degrees and about -105 degrees with intervals of a few tens of nanoseconds.

The Kinetics of Nucleation and Growth of Sickle Cell Hemoglobin Fibers

Polymerization of sickle cell hemoglobin (HbS) in deoxy state is one of the basic events in the pathophysiology of sickle cell anemia. For insight into the polymerization process, we monitor the kinetics of nucleation and growth of the HbS polymer fibers. We define a technique for the determination of the rates J and delay times theta of nucleation and the fiber growth rates R of deoxy-HbS fibers, based on photolysis of CO-HbS by laser illumination. We solve numerically time-dependent equations of heat conductance and CO transport, coupled with respective photo-chemical processes, during kinetics experiments under continuous illumination. After calibration with experimentally determined values, we define a regime of illumination ensuring uniform temperature and deoxy-HbS concentration, and fast (within <1 s) egress to steady conditions. With these procedures, data on the nucleation and growth kinetics have relative errors of <5% and are reproducible within 10% in independent experiments. The nucleation rates and delay times have steep, exponential dependencies on temperature. In contrast, the average fiber growth rates only weakly depend on temperature. The individual growth rates vary by up to 40% under identical conditions. These variations are attributed to instability of the coupled kinetics and diffusion towards the growing end of a fiber. The activation energy for incorporation of HbS molecules into a polymer is E(A)=50 kJ mol(-1), a low value indicating the significance of the hydrophobic contacts in the HbS polymer. More importantly, the contrast between the strong theta(T) and weak R(T) dependencies suggests that the homogenous nucleation of HbS polymers occurs within clusters of a precursor phase. This conclusion may have significant consequences for the understanding of the pathophysiology of sickle cell anemia and should be tested in further work.

Allosteric Effectors Influence the Tetramer Stability of Both R- and T-states of Hemoglobin A

The contribution of heterotropic effectors to hemoglobin allostery is still not completely understood. With the recently proposed global allostery model, this question acquires crucial significance, because it relates tertiary conformational changes to effector binding in both the R- and T-states. In this context, an important question is how far the induced conformational changes propagate from the binding site(s) of the allosteric effectors. We present a study in which we monitored the interdimeric interface when the effectors such as Cl-, 2,3-diphosphoglycerate, inositol hexaphosphate, and bezafibrate were bound. We studied oxy-Hb and a hybrid form (alphaFeO2)2-(betaZn)2 as the T-state analogue by monitoring heme absorption and Trp intrinsic fluorescence under hydrostatic pressure. We observed a pressure-dependent change in the intrinsic fluorescence, which we attribute to a pressure-induced tetramer to dimer transition with characteristic pressures in the 70-200-megapascal range. The transition is sensitive to the binding of allosteric effectors. We fitted the data with a simple model for the tetramer-dimer transition and determined the dissociation constants at atmospheric pressure. In the R-state, we observed a stabilizing effect by the allosteric effectors, although in the T-analogue a stronger destabilizing effect was seen. The order of efficiency was the same in both states, but with the opposite trend as inositol hexaphosphate > 2,3-diphosphoglycerate > Cl-. We detected intrinsic fluorescence from bound bezafibrate that introduced uncertainty in the comparison with other effectors. The results support the global allostery model by showing that conformational changes propagate from the effector binding site to the interdimeric interfaces in both quaternary states.

Circular dichroism spectroscopy of tertiary and quaternary conformations of human hemoglobin entrapped in wet silica gels

The relative contributions to changes in visible and near UV circular dichroism spectra of hemoglobin of heme ligation and tertiary and quaternary conformational transitions were separated by exploiting the slowing down of structural relaxations for proteins encapsulated in wet, nanoporous silica gels. Spectral signatures, previously assumed to be characteristic of T and R quaternary states, were demonstrated to be specific to different tertiary conformations. The results support the view that ligation and allosteric effectors can modulate the structural and functional properties of hemoglobin by regulating the equilibrium between the same tertiary species within both quaternary states.

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.

Quaternary structure of carbonmonoxyhemoglobins in solution: structural changes induced by the allosteric effector inositol hexaphosphate

We have applied the residual dipolar coupling (RDC) method to investigate the solution quaternary structures of (2)H- and (15)N-labeled human normal adult recombinant hemoglobin (rHb A) and a low-oxygen-affinity mutant recombinant hemoglobin, rHb(alpha96Val-->Trp), both in the carbonmonoxy form, in the absence and presence of an allosteric effector, inositol hexaphosphate (IHP), using a stretched polyacrylamide gel as the alignment medium. Our recent RDC results [Lukin, J. A., Kontaxis, G., Simplaceanu, V., Yuan, Y., Bax, A., and Ho, C. (2003) Proc. Natl. Acad. Sci. U.S.A. 100, 517-520] indicate that the quaternary structure of HbCO A in solution is a dynamic ensemble between two previously determined crystal structures, R (crystals grown under high-salt conditions) and R2 (crystals grown under low-salt conditions). On the basis of a comparison of the geometric coordinates of the T, R, and R2 structures, it has been suggested that the oxygenation of Hb A follows the transition pathway from T to R and then to R2, with R being the intermediate structure [Srinivasan, R., and Rose, G. D. (1994) Proc. Natl. Acad. Sci. U.S.A. 91, 11113-11117]. The results presented here suggest that IHP can shift the solution quaternary structure of HbCO A slightly toward the R structure. The solution quaternary structure of rHbCO(alpha96Val-->Trp) in the absence of IHP is similar to that of HbCO A in the presence of IHP, consistent with rHbCO(alpha96Val-->Trp) having an affinity for oxygen lower than that of Hb A. Moreover, IHP has a much stronger effect in shifting the solution quaternary structure of rHbCO(alpha96Val-->Trp) toward the R structure and toward the T structure, consistent with IHP causing a more pronounced decrease in its oxygen affinity. The results presented in this work, as well as other results recently reported in the literature, clearly indicate that there are multiple quaternary structures for the ligated form of hemoglobin. These results also provide new insights regarding the roles of allosteric effectors in regulating the structure and function of hemoglobin. The classical two-state/two-structure allosteric mechanism for the cooperative oxygenation of hemoglobin cannot account for the structural and functional properties of this protein and needs to be revised.

Modulation of reactivity and conformation within the T-quaternary state of human hemoglobin: the combined use of mutagenesis and sol-gel encapsulation

A range of conformationally distinct functional states within the T quaternary state of hemoglobin are accessed and probed using a combination of mutagenesis and sol-gel encapsulation that greatly slow or eliminate the T --> R transition. Visible and UV resonance Raman spectroscopy are used to probe the proximal strain at the heme and the status of the alpha(1)beta(2) interface, respectively, whereas CO geminate and bimolecular recombination traces in conjunction with MEM (maximum entropy method) analysis of kinetic populations are used to identify functionally distinct T-state populations. The mutants used in this study are Hb(Nbeta102A) and the alpha99-alpha99 cross-linked derivative of Hb(Wbeta37E). The former mutant, which binds oxygen noncooperatively with very low affinity, is used to access low-affinity ligated T-state conformations, whereas the latter mutant is used to access the high-affinity end of the distribution of T-state conformations. A pattern emerges within the T state in which ligand reactivity increases as both the proximal strain and the alpha(1)beta(2) interface interactions are progressively lessened after ligand binding to the deoxy T-state species. The ligation and effector-dependent interplay between the heme environment and the stability of the Trp beta37 cluster in the hinge region of the alpha(1)beta(2) interface appears to determine the distribution of the ligated T-state species generated upon ligand binding. A qualitative model is presented, suggesting that different T quaternary structures modulate the stability of different alphabeta dimer conformations within the tetramer.

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.

Electronic reorganization: Origin of sigma trans promotion effect

Binding of two ligands trans to each other by some transition metal complexes may be cooperative [Khoroshun et al., Mol Phys 2002, 100, 523]. Several interesting consequent effects include (i) inverse relationship between bond strength and binding affinity; (ii) smaller coordination barriers to formation of weaker bonds; (iii) enhancement of Lewis acidity with increased number of ligands. We describe a simple model, sigma trans promotion effect (TPE), which considers electronic reorganization between two Lewis structures, and predicts the above-mentioned effects. The applied result of present study is the unified perspective on several facts of heme chemistry. Particularly, we reiterate an important but often overlooked notion, developed previously within the spin pairing model [Drago and Corden, Acc Chem Res 1980, 13, 353], that, in hemoproteins, the proximal histidine and the distal ligand such as O2 or CO cooperate in promoting electronic reorganization. As a result, depopulation of dz2 orbital upon ligand binding contributes to the phenomenon of hemoglobin cooperativity. The presented density functional (B3LYP) calculations on realistic models, the processes of carbon monoxide binding by Fe(II) porphyrins and dinitrogen binding by triamido/triamidoamine Mo(III) complexes, particularly the evaluation of the coordination barriers due to spin-state change by location of the minima on seams of crossing, support the TPE model predictions. From a broader theoretical perspective, the present study would hopefully stimulate the development of much needed frameworks and tools for facile comparisons of wave functions and their properties between different geometries, species, and electronic states. Advancement of practical wave function comparisons may yield fresh qualitative perspectives on chemical reactivity, and promote better understanding of related concepts such as electronic reorganization.