Crystallographic, molecular modeling, and biophysical characterization of the valine beta 67 (E11)-->threonine variant of hemoglobin

X-ray crystallography and sickle cell disease drug discovery-a tribute to Donald Abraham

X-ray crystallography and structure-based drug discovery have played a major role in the discovery of antisickling agents that target hemoglobin (Hb) for the treatment of sickle cell disease (SCD). Sickle cell disease, the most common inherited hematologic disorder, occurs as a result of a single point mutation of βGlu6 in normal human adult hemoglobin (HbA) to βVal6 in sickle hemoglobin (HbS). The disease is characterized by polymerization of HbS and sickling of red blood cells (RBCs), leading to several secondary pathophysiologies, including but not limited to vaso-occlusion, hemolytic anemia, oxidative stress, inflammation, stroke, pain crisis, and organ damage. Despite the fact that SCD was the first disease to have its molecular basis established, the development of therapies was for a very long time a challenge and took several decades to find therapeutic agents. The determination of the crystal structure of Hb by Max Perutz in the early 60s, and the pioneering X-ray crystallography research by Donald J. Abraham in the early 80s, which resulted in the first structures of Hb in complex with small molecule allosteric effectors of Hb, gave much hope that structure-based drug discovery (SBDD) could be used to accelerate development of antisickling drugs that target the primary pathophysiology of hypoxia-induced HbS polymerization to treat SCD. This article, which is dedicated to Donald J. Abraham, briefly reviews structural biology, X-ray crystallography and structure-based drug discovery from the perspective of Hb. The review also presents the impact of X-ray crystallography in SCD drug development using Hb as a target, emphasizing the major and important contributions by Don Abraham in this field.

Hydrogen bonding and crystal packing favor a nonplanarCo(III)porphyrin conformation and unusually weak axial ligation in [Co(TPP)(benzylamine)2](SCN): A crystallographic and density functional theory investigation

The single crystal X-ray structure of [ Co ( TPP )( BzNH2)2]( SCN ), compound 2, where TPP = 5,10,15,20-tetraphenylporphyrin dianion and BzNH2= benzylamine, reveals that the SCN-ion is hydrogen-bonded to one of the coordinated amino group hydrogen atoms via its sulfur atom. Furthermore, the N – H ⋯ SCN interaction is balanced by a stronger N – H ⋯ NCS hydrogen bonding interaction for the trans BzNH2ligand as a result of the multiple hydrogen bond accepting character of the thiocyanate ion. Analysis of the crystal packing shows that these two hydrogen bonds play a major role in fixing unusual orientations for the axial ligands relative to the porphyrin ring in this system. This, in turn, leads to the formation of a nonplanar porphyrin core conformation that is a mixture of ruffle- and saddle-type distortions. The intricate hydrogen bonding between the cations and anions in 2 results in an unusually long mean Co – Naminecoordination distance of 2.033(4) Å, some 0.05 Å longer than previously observed for other bis(primary amine) complexes of Co(III) porphyrins with comparable porphyrin ligands. Density functional theory (DFT) calculations at the B3LYP/LACVP* level of theory have been used to gauge the perturbation of the electronic structure of the [ Co ( TPP )( BzNH2)2]+cation caused by the N – H ⋯ SCN and N – H ⋯ NCS hydrogen-bonded SCN-ions. The calculations show that partial mixing of the anion MOs with those of the porphyrin cation lead to changes in the electron populations of the 3d orbitals of up to 0.42 e as well as more nearly tetragonal electronic symmetry for the Co(III) ion as a result of adjustments of the relative energies of the MOs with predominantly 3d character.

Bis[2-(3-carboxyphenoxy)carbonylethyl]phosphinic acid (m-BCCEP): a novel affinity cross-linking reagent for the beta-cleft modification of human hemoglobin

The design and synthesis of bis[2-(3-carboxyphenoxy)carbonylethyl]phosphinic acid (m-BCCEP, 1) as a site-directed affinity reagent for cross-linking human hemoglobin have been reported as part of our long-term goal to generate artificial blood for emergency transfusions. Molecular modeling techniques were used to design the reagent, employing crystal coordinates of human hemoglobin A(0) imported from the Protein Data Bank. It was synthesized in four steps commencing from 3-hydroxybenzoic acid. The reagent 1 was converted to its trisodium salt to allow effective cross-linking in an aqueous medium. The reagent 1, as its trisodium salt, was found to specifically cross-link stroma-free human hemoglobin A(0) in the beta-cleft under oxygenated reaction conditions at neutral pH. The SDS-PAGE analyses of the modified hemoglobin pointed to the molecular mass range of 32 kDa as anticipated. The HPLC analyses of the product suggested that the cross-link had formed between the beta(1)-beta(2) subunits. Molecular dynamics simulation studies on the reagent-HbA(0) complex suggested that the predominant amino acid residues involved in the cross-linking are N-terminus Val-1 or Lys-82 on one of the beta-subunits and Lys-144 on the other. These predictions were borne out by MALDI-TOF MS analyses data of the peptide fragments obtained from tryptic digestion of the cross-linked product. The data also suggested the presence of a minor cross-link between Val-1 and Lys-82 on the opposing subunits. The oxygen equilibrium measurements of the m-BCCEP-modified hemoglobin product at 37 degrees C showed oxygen affinity (P(50) = 25.8 Torr) comparable to that of the natural whole blood (P(50) = 27.0 Torr) and significantly lower than that of stroma-free hemoglobin (P(50) = 14.19 Torr) assayed under identical conditions. The measured Hill coefficient value of 1.91 of the m-BCCEP-modified Hb product points to the reasonable retainment of oxygen-binding cooperativity after the cross-link formation.

Design of recombinant hemoglobins for use in transfusion fluids

Molecular biology has been applied to the development of hemoglobin-based oxygen carrier (HBOC) proteins that can be expressed in bacteria or yeast. The transformation of the hemoglobin molecule into an HBOC requires a variety of modifications for rendering the acellular molecule of hemoglobin physiologically acceptable when transfused in circulation. Hemoglobins with different oxygen affinities can be obtained by introducing mutations at the heme pocket, the site of oxygen binding, or by introducing surface mutations that stabilize the hemoglobin molecule in the low-oxygen-affinity state. Modification of the size of the heme pocket is also used to hinder nitric oxide depletion and associated vasoconstriction. Introduction of cysteine residues on the hemoglobin surface allows formation of intermolecular bonds and formation of polymeric HBOCs. These polymers of recombinant hemoglobin have the characteristics of molecular size, molecular stability, and oxygen delivery to hypoxic tissue suitable for an HBOC.

The Oxygen Transport System in Three Species of the Boreal Fish Family Gadidae

The Arctic and Antarctic marine faunas differ by age and isolation. Fishes of the two polar regions have undergone different regional histories that have driven the physiological diversities. Antarctic fish are highly stenothermal, in keeping with stable water temperatures, whereas Arctic fish, being exposed to seasonal temperature variations, exhibit higher physiological plasticity. This study reports the characterization of the oxygen transport system of three Arctic species of the family Gadidae, namely the Arctic cod Arctogadus glacialis, the polar cod Boreogadus saida, and the Atlantic cod Gadus morhua. Unlike Antarctic notothenioids, the blood displays high multiplicity, i.e. it has three hemoglobins, similar to many other acanthomorph teleosts. In the most abundant hemoglobin, oxygen binding is modulated by heterotropic effectors, with marked Bohr and Root effects. Remarkably, in two species (A. glacialis and B. saida), the Hill coefficient is very close to one in the whole pH range, indicating the apparent absence of cooperativity. The amino acid sequences have been used to gain insight into the evolution history of globins of polar fish. The results indicate that Arctic and Antarctic globins have different phylogenies and lead us to suggest that the selective pressure of environment stability allows the phylogenetic signal to be maintained in the Antarctic sequences, whereas environmental variability would tend to disrupt this signal in the Gadidae sequences.

Structural dynamics of myoglobin: ligand migration and binding in valine 68 mutants

We have combined Fourier transform infrared/temperature derivative (FTIR-TDS) spectroscopy at cryogenic temperatures and flash photolysis at ambient temperature to examine the effects of polar and bulky amino acid replacements of the highly conserved distal valine 68 in sperm whale myoglobin. In FTIR-TDS experiments, the CO ligand can serve as an internal voltmeter that monitors the local electrostatic field not only at the active site but also at intermediate ligand docking sites. Mutations of residue 68 alter size, shape, and electric field of the distal pocket, especially in the vicinity of the primary docking site (state B). As a consequence, the infrared bands associated with the ligand at site B are shifted. The effect is most pronounced in mutants with large aromatic side chains. Polar side chains (threonine or serine) have only little effect on the peak frequencies. Ligands that migrate toward more remote sites C and D give rise to IR bands with altered frequencies. TDS experiments separate the photoproducts according to their recombination temperatures. The rates and extent of ligand migration among internal cavities at cryogenic temperatures can be used to interpret geminate and bimolecular O2 and CO recombination at room temperature. The kinetics of geminate recombination can be explained by steric arguments alone, whereas both the polarity and size of the position 68 side chain play major roles in regulating bimolecular ligand binding from the solvent.

Introduction of a new regulatory mechanism into human hemoglobin

Previous studies on bovine hemoglobin (HbBv) have suggested amino acid substitutions, which might introduce into human hemoglobin (HbA) functional characteristics of HbBv, namely a low intrinsic oxygen affinity regulated by Cl(-). Accordingly, we have constructed and characterized a multiple mutant, PB5, [beta(V1M + H2 Delta + T4I + P5A + A76K)] replacing four amino acid residues of HbA with those present at structurally analogous positions in HbBv, plus an additional substitution, beta T4I, which does not occur in either HbBv or HbA. This 'pseudobovine' hemoglobin has oxygen binding properties very similar to those of HbBv: the P(50) of HbA, PB5 and HbBv in the absence of Cl(-) are 1.6, 4.6 and 4.8 torr, respectively, and in 100 mM Cl(-) are 3.7, 10.5 and 12 torr, respectively. Moreover, PB5 has 3-fold slower autoxidation rate compared to HbA and HbBv. These are desirable characteristics for a human hemoglobin to be considered for use as a clinical artificial oxygen carrier. Although the functional properties of PB5 and HbBv are similar, van't Hoff plots indicate that the two hemoglobins interact differently with water, suggesting that factors regulating the R to T equilibrium are not the same in the two proteins. A further indication that PB5 is not a functional mimic of HbBv derives from PB5(control), a human hemoglobin with the same substitutions as PB5, except the beta T4I replacement. PB5(control) has a high oxygen affinity (P(50)=2.3 torr) in the absence of Cl(-), but retains the Cl(-) effect of PB5. The Cl(-) regulation of oxygen affinity in PB5 involves lysine residues at beta 8 and beta 76. PB4, which has the same substitutions as PB5 except beta A76K, and PB6, which has all the substitutions of PB5 plus beta K8Q, both have a low intrinsic oxygen affinity, like HbBv and PB5, but exhibit a decreased sensitivity to Cl(-). Since HbBv has lysine residues at both beta 8 and beta 76, these results imply that Cl(-) regulation in HbBv likewise involves these two residues. The mechanism responsible for the low intrinsic oxygen affinity of HbBv remains unclear. It is suggested that residues peculiar to HbBv at the alpha(1)beta(1) interface may play a role.

The Heme–Globin and Dimerization Equilibria of Recombinant Human Hemoglobins Carrying Site-Specific β Chains Mutations

The heme-globin and dimer-tetramer equilibria of ferric recombinant human hemoglobins with site-specific beta chain mutations at the heme pocket or at either the a1beta1 or the alpha1beta2 interfaces have been determined. The heme pocket mutation V67T leads to a marked stabilization of the beta chain heme and does not affect the dimer-tetramer association constant, K2,4. In the C112 mutants, the intrinsic rate of beta chain heme loss with respect to recombinant HbA (HbA-wt) is significantly increased only in C112G with some heme released also from the alpha chains. Gel filtration experiments indicate that the K2,4 value is essentially unaltered in C112G and C112L, but is increased in C112V and decreased in C112N. Substitution of cysteine 93 with A or M leads to a slight decrease of the rate of beta chain heme release, whereas the obvserved K2,4 value is similar to that obtained for HbA-wt. Modifications in oxygen affinity were observed in all the mutant hemoglobins with the exception of V67T, C93A, and C112G. The data indicate that there is no correlation between tetramer stability, beta chain heme affinity, and hemoglobin functionality and therefore point to a separate regulation of these properties.

Properties of human hemoglobins with increased polarity in the alpha- or beta-heme pocket. Carbonmonoxy derivatives

The spectroscopic, conformational, and functional properties of mutant carbonmonoxy hemoglobins in which either the beta-globin Val67(E11) or the alpha-globin Val62(E11) is replaced by threonine have been investigated. The thermal evolution of the Soret absorption band and the stretching frequency of the bound CO were used to probe the stereodynamic properties of the heme pocket. The functional properties were investigated by kinetic measurements. The spectroscopic and functional data were related to the conformational properties through molecular analysis. The effects of this nonpolar-to-polar isosteric mutation are: (i) increase of heme pocket anharmonic motions, (ii) stabilization of the A0 conformer in the IR spectrum, (iii) increased CO dissociation rates. The spectroscopic data indicate that for the carbonmonoxy derivatives, the Val --> Thr mutation has a larger conformational effect on the beta-subunits than on the alpha-subunits. This is at variance with the deoxy derivatives where the conformational modification was larger in the heme pocket of the alpha-subunit (Cupane, A., Leone, M., Militello, V., Friedman, R. K., Koley, A. P., Vasquez, G. P., Brinigar, W. S., Karavitis, M., and Fronticelli, C. (1997) J. Biol. Chem. 272, 26271-26278). These effects are attributed to a different electrostatic interaction between Ogamma of Thr(E11) and the bound CO molecule. Molecular analysis indicates a more favorable interaction of the bound CO with Thr Ogamma in the beta-subunit heme pocket.

Bis[2-(4-carboxyphenoxy)carbonylethyl]phosphinic acid (BCCEP): a novel affinity reagent for the beta-cleft modification of human hemoglobin

The design, synthesis, and hemoglobin cross-linking studies of a novel organic reagent, bis[2-(4-carboxyphenoxy)carbonylethyl]phosphinic acid (BCCEP, 1) have been reported. The reagent was designed with the aid of molecular modeling, employing crystal coordinates of human hemoglobin A0. It was synthesized in three steps commencing from 4-t-butoxycarbonylphenol. The tri-sodium salt of 1 was employed to cross-link human oxyHb. While SDS-PAGE analyses of the modified hemoglobin product pointed to the molecular mass range of 32 kDa, the HPLC analyse suggested that the cross-link had formed between the beta 1-beta 2 subunits. The oxygen equilibrium measurements of the modified hemoglobin at 37 degrees C showed significantly reduced oxygen affinity (P50 = 31.3 Torr) as compared with that of cell-free hemoglobin (P50 = 6.6 Torr). The sigmoidal shape of O2 curves of the modified Hb pointed to reasonable retainment of oxygen-binding cooperativity after the cross-link formation. Molecular dynamics simulation studies on the reagent-HbA0 complex suggested that the most likely amino acid residues involved in the cross-linking are N-terminus Val-1 or Lys-82 on one of the-chains, and Lys-144 on the other. These predictions were consistent with the results of MALDI-MS analyses of the peptide fragments obtained from tryptic digestion of the cross-linked product.

Conformation of the sebacyl beta1Lys82-beta2Lys82 crosslink in T-state human hemoglobin

The crystal structure of human T state hemoglobin crosslinked with bis(3,5-dibromo-salicyl) sebacate has been determined at 1.9 A resolution. The final crystallographic R factor is 0.168 with root-mean-square deviations (RMSD) from ideal bond distance of 0.018 A. The 10-carbon sebacyl residue found in the beta cleft covalently links the two betaLys82 residues. The sebacyl residue assumes a zigzag conformation with cis amide bonds formed by the NZ atoms of betaLys82's and the sebacyl carbonyl oxygens. The atoms of the crosslink have an occupancy factor of 1.0 with an average temperature factor for all atoms of 34 A2. An RMSD of 0.27 for all CA's of the tetramer is observed when the crosslinked deoxyhemoglobin is compared with deoxyhemoglobin refined by using a similar protocol, 2HHD [Fronticelli et al. J. Biol. Chem. 269: 23965-23969, 1994]. Thus, no significant perturbations in the tertiary or quaternary structure are introduced by the presence of the sebacyl residue. However, the sebacyl residue does displace seven water molecules in the beta cleft and the conformations of the beta1Lys82 and beta2Lys82 are altered because of the crosslinking. The carbonyl oxygen that is part of the amide bond formed with the NZ of beta2Lys82 forms a hydrogen bond with side chain of beta2Asn139 that is in turn hydrogen-bonded to the side chain of beta2Arg104. A comparison of the observed conformation with that modeled [Bucci et al. Biochemistry 35:3418-3425, 1996] shows significant differences. The differences in the structures can be rationalized in terms of compensating changes in the estimated free-energy balance, based on differences in exposed surface areas and the observed shift in the side-chain hydrogen-bonding pattern involving beta2Arg104, beta2Asn139, and the associated sebacyl carbonyl group.

Alpha-subunit oxidation in T-state crystals of a sebacyl cross-linked human hemoglobin with unusual autoxidation properties

In the crystal structure of human T-state hemoglobin with a sebacyl residue cross-linking the two beta-subunit Lys82,s (DecHb), the Fe atoms of the alpha-subunit hemes are found to be oxidized with a water molecule bound. The three-dimensional structure and heme geometries were compared to those of deoxyhemoglobin and other partially and fully oxidized hemoglobins [R. Liddington, Z. Derewenda, E. Dodson, R. Hubbard, G. Dodson, J. Mol. Biol. 228 (1992) 551]. The heme geometries of the alpha-subunits are consistent with those observed in oxidized structures. The proximal histidines of the alpha-subunits move toward the heme plane shifting the F-helix and FG-corner in a manner observed for partially oxidized human hemoglobin. This supports the hypothesis that these perturbations may precede the T- to R-state transition. Circular dichroism studies comparing DecHb and natural human hemoglobin in the deoxy and CO ligated forms confirm that the conformations of the deoxy forms are identical, but the ligated forms have slight differences in the solution structures. DecHb is found to be more resistant to autoxidation than natural hemoglobin. The time course of autoxidation of DecHb shows that it is virtually absent for the first 1500 min followed by a rapid increase. Thus, the discovery of the oxidation of the alpha-subunits in the deoxy-crystals is quite unexpected. The data confirm that ligation of the alpha-subunits precedes that of the beta-subunits. This may suggest a low ligand affinity of the alpha-diligated form of hemoglobin.

Modification of alpha-chain or beta-chain heme pocket polarity by Val(E11) --> thr substitution has different effects on the steric, dynamic, and functional properties of human recombinant hemoglobin. Deoxy derivatives

The dynamic and functional properties of mutant deoxyhemoglobins in which either the beta-globin Val67(E11) or the alpha-globin Val62(E11) is replaced by threonine have been investigated through the thermal evolution of the Soret absorption band in the temperature range 300 to 20 K and through the kinetics of CO rebinding after flash photolysis at room temperature. The conformational properties of the modified alpha chain and beta chain distal heme pockets were also studied through x-ray crystallography and molecular modeling. The data obtained with the various techniques consistently indicate that the polar isosteric mutation in the distal side of the alpha chain heme pocket has a larger effect on the investigated properties than the analogous mutation on the beta chain. We attribute the observed differences to the presence of a water molecule in the distal heme pocket of the modified alpha chains, interacting with the hydroxyl of the threonine side chain. This is indicated by molecular modeling which showed that the water molecule present in the alpha chain distal heme pocket can bridge by H bonding between Thr62(E11) and His58(E7) without introducing any unfavorable steric interactions. Consistent with the dynamic and functional data, the presence of a water molecule in the distal heme pocket of the modified beta chains is not observed by x-ray crystallography.