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

Crystal structure of carbonmonoxy sickle hemoglobin in R-state conformation

The fundamental pathophysiology of sickle cell disease is predicated by the polymerization of deoxygenated (T-state) sickle hemoglobin (Hb S) into fibers that distort red blood cells into the characteristic sickle shape. The crystal structure of deoxygenated Hb S (DeoxyHb S) and other studies suggest that the polymer is initiated by a primary interaction between the mutation βVal6 from one Hb S molecule, and a hydrophobic acceptor pocket formed by the residues βAla70, βPhe85 and βLeu88 of an adjacent located Hb S molecule. On the contrary, oxygenated or liganded Hb S does not polymerize or incorporate in the polymer. In this paper we present the crystal structure of carbonmonoxy-ligated sickle Hb (COHb S) in the quaternary classical R-state at 1.76Å. The overall structure and the pathological donor and acceptor environments of COHb S are similar to those of the isomorphous CO-ligated R-state normal Hb (COHb A), but differ significantly from DeoxyHb S as expected. More importantly, the packing of COHb S molecules does not show the typical pathological interaction between βVal6 and the βAla70, βPhe85 and βLeu88 hydrophobic acceptor pocket observed in DeoxyHb S crystal. The structural analysis of COHb S, COHb A and DeoxyHb S provides atomic level insight into why liganded hemoglobin does not form a polymer.

Design, Synthesis, and Investigation of Novel Nitric Oxide (NO)-Releasing Prodrugs as Drug Candidates for the Treatment of Ischemic Disorders: Insights into NO-Releasing Prodrug Biotransformation and Hemoglobin-NO Biochemistry

We have developed novel nitric oxide (NO)-releasing prodrugs of efaproxiral (RSR13) for their potential therapeutic applications in a variety of diseases with underlying ischemia. RSR13 is an allosteric effector of hemoglobin (Hb) that decreases the protein's affinity for oxygen, thereby increasing tissue oxygenation. NO, because of its vasodilatory property, in the form of ester prodrugs has been found to be useful in managing several cardiovascular diseases by increasing blood flow and oxygenation in ischemic tissues. We synthesized three NO-donor ester derivatives of RSR13 (DD-1, DD-2, and DD-3) by attaching the NO-releasing moieties nitrooxyethyl, nitrooxypropyl, and 1-(pyrrolidin-1-yl)diazen-1-ium-1,2-diolate, respectively, to the carboxylate of RSR13. In vitro studies demonstrated that the compounds released NO in a time-dependent manner upon being incubated with l-cysteine (1.8-9.3%) or human serum (2.3-52.5%) and also reduced the affinity of Hb for oxygen in whole blood (ΔP50 of 4.9-21.7 mmHg vs ΔP50 of 25.4-32.1 mmHg for RSR13). Crystallographic studies showed RSR13, the hydrolysis product of the reaction between DD-1 and deoxygenated Hb, bound to the central water cavity of Hb. Also, the hydrolysis product, NO, was observed exclusively bound to the two α hemes, the first such HbNO structure to be reported, capturing the previously proposed physiological bis-ligated nitrosylHb species. Finally, nitrate was observed bound to βHis97. Ultraperformance liquid chromatography-mass spectrometry analysis of the compounds incubated with matrices used for the various studies demonstrated the presence of the predicted reaction products. Our findings, beyond the potential therapeutic application, provide valuable insights into the biotransformation of NO-releasing prodrugs and their mechanism of action and into hemoglobin-NO biochemistry at the molecular level.

Identification of a novel class of covalent modifiers of hemoglobin as potential antisickling agents

Aromatic aldehydes and ethacrynic acid (ECA) exhibit antipolymerization properties that are beneficial for sickle cell disease therapy. Based on the ECA pharmacophore and its atomic interaction with hemoglobin, we designed and synthesized several compounds - designated as KAUS (imidazolylacryloyl derivatives) - that we hypothesized would bind covalently to βCys93 of hemoglobin and inhibit sickling. The compounds surprisingly showed weak allosteric and antisickling properties. X-ray studies of hemoglobin in complex with representative KAUS compounds revealed an unanticipated mode of Michael addition between the β-unsaturated carbon and the N-terminal αVal1 nitrogen at the α-cleft of hemoglobin, with no observable interaction with βCys93. Interestingly, the compounds exhibited almost no reactivity with the free amino acids, L-Val, L-His and L-Lys, but showed some reactivity with both glutathione and L-Cys. Our findings provide a molecular level explanation for the compounds biological activities and an important framework for targeted modifications that would yield novel potent antisickling agents.

Identification of a small molecule that increases hemoglobin oxygen affinity and reduces SS erythrocyte sickling

Small molecules that increase the oxygen affinity of human hemoglobin may reduce sickling of red blood cells in patients with sickle cell disease. We screened 38 700 compounds using small molecule microarrays and identified 427 molecules that bind to hemoglobin. We developed a high-throughput assay for evaluating the ability of the 427 small molecules to modulate the oxygen affinity of hemoglobin. We identified a novel allosteric effector of hemoglobin, di(5-(2,3-dihydro-1,4-benzodioxin-2-yl)-4H-1,2,4-triazol-3-yl)disulfide (TD-1). TD-1 induced a greater increase in oxygen affinity of human hemoglobin in solution and in red blood cells than did 5-hydroxymethyl-2-furfural (5-HMF), N-ethylmaleimide (NEM), or diformamidine disulfide. The three-dimensional structure of hemoglobin complexed with TD-1 revealed that monomeric units of TD-1 bound covalently to β-Cys93 and β-Cys112, as well as noncovalently to the central water cavity of the hemoglobin tetramer. The binding of TD-1 to hemoglobin stabilized the relaxed state (R3-state) of hemoglobin. TD-1 increased the oxygen affinity of sickle hemoglobin and inhibited in vitro hypoxia-induced sickling of red blood cells in patients with sickle cell disease without causing hemolysis. Our study indicates that TD-1 represents a novel lead molecule for the treatment of patients with sickle cell disease.

Differential control of heme reactivity in alpha and beta subunits of hemoglobin: a combined Raman spectroscopic and computational study

The use of hybrid hemoglobin (Hb), with mesoheme substituted for protoheme, allows separate monitoring of the α or β hemes along the allosteric pathway. Using resonance Raman (rR) spectroscopy in silica gel, which greatly slows protein motions, we have observed that the Fe-histidine stretching frequency, νFeHis, which is a monitor of heme reactivity, evolves between frequencies characteristic of the R and T states, for both α or β chains, prior to the quaternary R-T and T-R shifts. Computation of νFeHis, using QM/MM and the conformational search program PELE, produced remarkable agreement with experiment. Analysis of the PELE structures showed that the νFeHis shifts resulted from heme distortion and, in the α chain, Fe-His bond tilting. These results support the tertiary two-state model of ligand binding (Henry et al., Biophys. Chem. 2002, 98, 149). Experimentally, the νFeHis evolution is faster for β than for α chains, and pump-probe rR spectroscopy in solution reveals an inflection in the νFeHis time course at 3 μs for β but not for α hemes, an interval previously shown to be the first step in the R-T transition. In the α chain νFeHis dropped sharply at 20 μs, the final step in the R-T transition. The time courses are fully consistent with recent computational mapping of the R-T transition via conjugate peak refinement by Karplus and co-workers (Fischer et al., Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 5608). The effector molecule IHP was found to lower νFeHis selectively for α chains within the R state, and a binding site in the α1α2 cleft is suggested.

Structure of liganded T-state haemoglobin from cat (Felis silvestris catus), a low oxygen-affinity species, in two different crystal forms

Haemoglobin (Hb) is an iron-containing metalloprotein which plays a major role in the transportation of oxygen from the lungs to tissues and of carbon dioxide back to the lungs. Hb is in equilibrium between low-affinity tense (T) and high-affinity relaxed (R) states associated with its unliganded and liganded forms, respectively. Mammalian species can be classified into two groups on the basis of whether they express `high' or `low' oxygen-affinity Hbs. Although Hbs from the former group have been studied extensively, a more limited number of structural studies have been performed for low oxygen-affinity Hbs. Here, the crystal structure of low oxygen-affinity cat methaemoglobin (metHb) has been solved at 2.0 and 2.4 Å resolution in two different crystal forms. Even though both structures are fully liganded, they unusually adopt a T-state-like quaternary conformation but with several localized R-like tertiary-structural and quaternary-structural features. The study provides atomic-level insights into the ligand-binding properties of this Hb, including its low cooperativity, blunt response to allosteric effectors and low affinity for oxygen, as well as further contributing to the mechanism underlying Hb allostery.

Crystallographic characterization of the nitric oxide derivative of R-state human hemoglobin

Nitric oxide (NO) is a signaling agent that is biosynthesized in vivo. NO binds to the heme center in human hemoglobin (Hb) to form the HbNO adduct. This reaction of NO with Hb has been studied for many decades. Of continued interest has been the effect that the bound NO ligand has on the geometrical parameters of the resulting heme-NO active site. Although the crystal structure of a T-state human HbNO complex has been published previously, that of the high affinity R-state HbNO derivative has not been reported to date. We have crystallized and solved the three-dimensional X-ray structure of R-state human HbNO to 1.90 Å resolution. The differences in the FeNO bond parameters and H-bonding patterns between the α and β subunits contribute to understanding of the observed enhanced stability of the α(FeNO) moieties relative to the β(FeNO) moieties in human R-state HbNO.

Non-site-specific allosteric effect of oxygen on human hemoglobin under high oxygen partial pressure

Protein allostery is essential for vital activities. Allosteric regulation of human hemoglobin (HbA) with two quaternary states T and R has been a paradigm of allosteric structural regulation of proteins. It is widely accepted that oxygen molecules (O₂) act as a “site-specific” homotropic effector, or the successive O₂ binding to the heme brings about the quaternary regulation. However, here we show that the site-specific allosteric effect is not necessarily only a unique mechanism of O₂ allostery. Our simulation results revealed that the solution environment of high O₂ partial pressure enhances the quaternary change from T to R without binding to the heme, suggesting an additional “non-site-specific” allosteric effect of O₂. The latter effect should play a complementary role in the quaternary change by affecting the intersubunit contacts. This analysis must become a milestone in comprehensive understanding of the allosteric regulation of HbA from the molecular point of view.

Therapeutic strategies to alter the oxygen affinity of sickle hemoglobin

The pathophysiology of sickle cell disease involves the polymerization of sickle hemoglobin in its T state, which develops under low oxygen saturation. One therapeutic strategy is to develop pharmacologic agents to stabilize the R state of hemoglobin, which has higher oxygen affinity and is expected to have slower kinetics of polymerization, potentially delaying the sickling of red cells during circulation. This strategy has stimulated the investigation of aromatic aldehydes, aspirin derivatives, thiols, and isothiocyanates that can stabilize the R state of hemoglobin in vitro. One representative aromatic aldehyde agent, 5-hydoxymethyl-2-furfural, protects sickle cell mice from the effects of hypoxia.

Structure of fully liganded Hb ζ2β2s trapped in a tense conformation

A variant Hb ζ2β2(s) that is formed from sickle hemoglobin (Hb S; α2β2(s)) by exchanging adult α-globin with embryonic ζ-globin subunits shows promise as a therapeutic agent for sickle-cell disease (SCD). Hb ζ2β2(s) inhibits the polymerization of deoxygenated Hb S in vitro and reverses characteristic features of SCD in vivo in mouse models of the disorder. When compared with either Hb S or with normal human adult Hb A (α2β2), Hb ζ2β2(s) exhibits atypical properties that include a high oxygen affinity, reduced cooperativity, a weak Bohr effect and blunted 2,3-diphosphoglycerate allostery. Here, the 1.95 Å resolution crystal structure of human Hb ζ2β2(s) that was expressed in complex transgenic knockout mice and purified from their erythrocytes is presented. When fully liganded with carbon monoxide, Hb ζ2β2(s) displays a central water cavity, a ζ1-β(s)2 (or ζ2-β(s)1) interface, intersubunit salt-bridge/hydrogen-bond interactions, C-terminal βHis146 salt-bridge interactions, and a β-cleft, that are highly unusual for a relaxed hemoglobin structure and are more typical of a tense conformation. These quaternary tense-like features contrast with the tertiary relaxed-like conformations of the ζ1β(s)1 dimer and the CD and FG corners, as well as the overall structures of the heme cavities. This crystallographic study provides insights into the altered oxygen-transport properties of Hb ζ2β2(s) and, moreover, decouples tertiary- and quaternary-structural events that are critical to Hb ligand binding and allosteric function.

Coarse-grained and all-atom modeling of structural states and transitions in hemoglobin

Hemoglobin (Hb), an oxygen-binding protein composed of four subunits (α1, α2, β1, and β2), is a well-known example of allosteric proteins that are capable of cooperative ligand binding. Despite decades of studies, the structural basis of its cooperativity remains controversial. In this study, we have integrated coarse-grained (CG) modeling, all-atom simulation, and structural data from X-ray crystallography and wide-angle X-ray scattering (WAXS), aiming to probe dynamic properties of the two structural states of Hb (T and R state) and the transitions between them. First, by analyzing the WAXS data of unliganded and liganded Hb, we have found that the structural ensemble of T or R state is dominated by one crystal structure of Hb with small contributions from other crystal structures of Hb. Second, we have used normal mode analysis to identify two distinct quaternary rotations between the α1β1 and α2β2 dimer, which drive the transitions between T and R state. We have also identified the hot-spot residues whose mutations are predicted to greatly change these quaternary motions. Third, we have generated a CG transition pathway between T and R state, which predicts a clear order of quaternary and tertiary changes involving α and β subunits in Hb. Fourth, we have used the accelerated molecular dynamics to perform an all-atom simulation starting from the T state of Hb, and we have observed a transition toward the R state of Hb. Further analysis of crystal structural data and the all-atom simulation trajectory has corroborated the order of quaternary and tertiary changes predicted by CG modeling.

Unlocking the binding and reaction mechanism of hydroxyurea substrates as biological nitric oxide donors

Hydroxyurea is the only FDA approved treatment of sickle cell disease. It is believed that the primary mechanism of action is associated with the pharmacological elevation of nitric oxide in the blood; however, the exact details of this are still unclear. In the current work, we investigate the atomic level details of this process using a combination of flexible-ligand/flexible-receptor virtual screening coupled with energetic analysis that decomposes interaction energies. Utilizing these methods, we were able to elucidate the previously unknown substrate binding modes of a series of hydroxyurea analogs to hemoglobin and the concomitant structural changes of the enzyme. We identify a backbone carbonyl that forms a hydrogen bond with bound substrates. Our results are consistent with kinetic and electron paramagnetic resonance (EPR) measurements of hydroxyurea-hemoglobin reactions, and a full mechanism is proposed that offers new insights into possibly improving substrate binding and/or reactivity.

Crystallographic analysis of human hemoglobin elucidates the structural basis of the potent and dual antisickling activity of pyridyl derivatives of vanillin

Vanillin has previously been studied clinically as an antisickling agent to treat sickle-cell disease. In vitro investigations with pyridyl derivatives of vanillin, including INN-312 and INN-298, showed as much as a 90-fold increase in antisickling activity compared with vanillin. The compounds preferentially bind to and modify sickle hemoglobin (Hb S) to increase the affinity of Hb for oxygen. INN-312 also led to a considerable increase in the solubility of deoxygenated Hb S under completely deoxygenated conditions. Crystallographic studies of normal human Hb with INN-312 and INN-298 showed that the compounds form Schiff-base adducts with the N-terminus of the α-subunits to constrain the liganded (or relaxed-state) Hb conformation relative to the unliganded (or tense-state) Hb conformation. Interestingly, while INN-298 binds and directs its meta-positioned pyridine-methoxy moiety (relative to the aldehyde moiety) further down the central water cavity of the protein, that of INN-312, which is ortho to the aldehyde, extends towards the surface of the protein. These studies suggest that these compounds may act to prevent sickling of SS cells by increasing the fraction of the soluble high-affinity Hb S and/or by stereospecific inhibition of deoxygenated Hb S polymerization.

Structures and oxygen affinities of crystalline human hemoglobin C (β6 Glu->Lys) in the R and R2 quaternary structures

Recent crystallographic studies suggested that fully liganded human hemoglobin can adopt multiple quaternary conformations that include the two previously solved relaxed conformations, R and R2, whereas fully unliganded deoxyhemoglobin may adopt only one T (tense) quaternary conformation. An important unanswered question is whether R, R2, and other relaxed quaternary conformations represent different physiological states with different oxygen affinities. Here, we answer this question by showing the oxygen equilibrium curves of single crystals of human hemoglobin in the R and R2 state. In this study, we have used a naturally occurring mutant hemoglobin C (β6 Glu→Lys) to stabilize the R and R2 crystals. Additionally, we have refined the x-ray crystal structure of carbonmonoxyhemoglobin C, in the R and R2 state, to 1.4 and 1.8 Å resolution, respectively, to compare precisely the structures of both types of relaxed states. Despite the large quaternary structural difference between the R and R2 state, both crystals exhibit similar noncooperative oxygen equilibrium curves with a very high affinity for oxygen, comparable with the fourth oxygen equilibrium constant (K(4)) of human hemoglobin in solution. One small difference is that the R2 crystals have an oxygen affinity that is 2-3 times higher than that of the R crystals. These results demonstrate that the functional difference between the two typical relaxed quaternary conformations is small and physiologically less important, indicating that these relaxed conformations simply reflect a structural polymorphism of a high affinity relaxed state.