Refined crystal structure of deoxyhemoglobin S. I. Restrained least-squares refinement at 3.0-A resolution

Cyclic Peptides as Aggregation Inhibitors for Sickle Cell Disease

Sickle cell disease is a missense genetic disorder characterized by the aggregation of deoxy-HbS into helical fibers that distort erythrocytes into a sickle-like shape. Herein, we investigate, through molecular dynamics, the effect of nine 5-mer cyclic peptides (CPs), tailor-designed to block key lateral contacts of the fibers. Our results show that the CPs bind orthogonally to the main HbS pocket involved in the latter contacts, with some revealing exceedingly long residence times. These CPs display moderate to high specificity, exhibiting molecular recognition events even at a HbS/CP (1:1) ratio. A much lower HbS-CP binding free energy, longer residence times, and higher specificity are also found relative to a previously reported CP with modest in vitro antisickling activity. These results indicate that some of these CPs have the potential to reduce the concentration of aggregation-competent deoxy-HbS, precluding or delaying the formation of lateral contact at the homogeneous nucleation stage.

Voxelotor does not inhibit sickle hemoglobin fiber formation upon complete deoxygenation

The drug voxelotor (commercially known as Oxbryta) has been approved by the US Food and Drug Administration for the treatment of sickle cell disease. It is known to reduce disease-causing sickling by inhibiting the transformation of the non-polymerizing, high-oxygen-affinity R quaternary structure of sickle hemoglobin into its polymerizing, low-affinity T quaternary structure. It has not been established whether the binding of the drug has anti-sickling effects beyond restricting the change of quaternary structure. By using a laser photolysis method that employs microscope optics, we have determined that fully deoxygenated sickle hemoglobin will assume the T structure. We show that the nucleation rates essential to generate the sickle fibers are not significantly affected by voxelotor. The method employed here should be useful for determining the mechanism of sickling inhibition for proposed drugs.

The Sickle-Cell Fiber Revisited

Sickle cell disease is the consequence of a single point mutation on the surface of the β chains of the hemoglobin molecule leading to the formation of rigid polymers that disrupt circulation. It has long been established that the polymers are comprised of seven pairs of double strands that are twisted replicas of the double strands found in crystals. Here, we review several newer developments that elaborate on that simple model and provide deeper insights into the process.

Molecular Dynamics of Hemoglobin Reveals Structural Alterations and Explains the Interactions Driving Sickle Cell Fibrillation

In sickle cell anemia, deoxyhemoglobin deforms RBCs by forming fibrils inside that disintegrate on oxygenation. We studied 100 ns long all-atom molecular dynamics (MD) for sickle and normal hemoglobin fibril models to understand this process, complemented by multiple 1 μs MD for a single tetramer of sickle and normal hemoglobin in deoxy and oxy states. We find that the presence of hydrophobic residues without a bulky side chain at β-6 in hemoglobin is the reason for the stability of the fibrils. Moreover, the free energy landscapes from MD of hemoglobin starting in the tensed (T) state capture the putative transition from T to relaxed (R) state, associated with oxygen binding. The three conformational wells in the landscapes are characterized by the quaternary changes where one αβ dimer rotates with respect to the other. The conformational changes from the oxygenation of sickle hemoglobin hinder the intermolecular contacts necessary for fibril formation.

Allosteric control of hemoglobin S fiber formation by oxygen and its relation to the pathophysiology of sickle cell disease

The pathology of sickle cell disease is caused by polymerization of the abnormal hemoglobin S upon deoxygenation in the tissues to form fibers in red cells, causing them to deform and occlude the circulation. Drugs that allosterically shift the quaternary equilibrium from the polymerizing T quaternary structure to the nonpolymerizing R quaternary structure are now being developed. Here we update our understanding on the allosteric control of fiber formation at equilibrium by showing how the simplest extension of the classic quaternary two-state allosteric model of Monod, Wyman, and Changeux to include tertiary conformational changes provides a better quantitative description. We also show that if fiber formation is at equilibrium in vivo, the vast majority of cells in most tissues would contain fibers, indicating that it is unlikely that the disease would be survivable once the nonpolymerizing fetal hemoglobin has been replaced by adult hemoglobin S at about 1 y after birth. Calculations of sickling times, based on a recently discovered universal relation between the delay time prior to fiber formation and supersaturation, show that in vivo fiber formation is very far from equilibrium. Our analysis indicates that patients survive because the delay period allows the majority of cells to escape the small vessels of the tissues before fibers form. The enormous sensitivity of the duration of the delay period to intracellular hemoglobin composition also explains why sickle trait, the heterozygous condition, and the compound heterozygous condition of hemoglobin S with pancellular hereditary persistence of fetal hemoglobin are both relatively benign conditions.

On the Nonaggregation of Normal Adult Hemoglobin and the Aggregation of Sickle Cell Hemoglobin

Sickle cell disease is a genetic disorder associated with a single mutation (Glu-β6 → Val-β6) in the β chains of hemoglobin, causing the polymerization of deoxygenated sickle cell hemoglobin (deoxy-HbS). The deoxy-HbS binding free energy was recently studied through molecular simulations, and a value of -14 ± 1 kcal mol^-1 was found. Here, we studied the binding free energy of normal adult hemoglobin (deoxy-HbA), which does not polymerize at normal physiological conditions, with the aim of elucidating the importance of the presence of Val-β6 and of the absence of Glu-β6 on the aggregation of deoxy-HbS. A binding free energy of -4.4 ± 0.5 kcal mol^-1 was found from a one-dimensional potential of mean force. Hydrophobic interactions are shown to represent less than 20% of the interactions in the contact interface, and despite similarly strong hydrogen-bonded ion pairs (i.e., salt bridges) and water bridged electrostatic interactions are found for deoxy-HbA and deoxy-HbS, a large repulsive potential energy is associated with Glu-β6, whereas a mild attractive potential energy is connected with Val-β6. Interestingly, Asp-β73 switches from forming a major electrostatic repulsive pair with Glu-β6 in deoxy-HbA, to forming a major attractive residue pair with Val-β6 in deoxy-HbS, consistent with the view that damping of electrostatic repulsions involving Glu-β6, namely, those associated with Asp-β73, could be responsible for the polymerization of deoxy-HbA at high potassium phosphate concentrations. Solvation analysis shows that functional groups forming salt bridges and water bridged interactions preserve a nearly intact first hydration sphere, avoiding a complete dewetting free energy penalty. These results support the view that the absence of Glu-β6 is more important than the presence of Val-β6, and that although hydrophobic effects, associated with the Val-β6 dehydration and interaction with the hydrophobic pocket in the neighbor tetramer, are important, electrostatic interactions are dominant, opposite to a picture where HbS association is driven by hydrophobic interactions.

Increased sequence hydrophobicity reduces conformational specificity: A mutational case study of the Arc repressor protein

The amino-acid sequences of soluble, globular proteins must have hydrophobic residues to form a stable core, but excess sequence hydrophobicity can lead to loss of native state conformational specificity and aggregation. Previous studies of polar-to-hydrophobic mutations in the β-sheet of the Arc repressor dimer showed that a single substitution at position 11 (N11L) leads to population of an alternate dimeric fold in which the β-sheet is replaced by helix. Two additional hydrophobic mutations at positions 9 and 13 (Q9V and R13V) lead to population of a differently folded octamer along with both dimeric folds. Here we conduct a comprehensive study of the sequence determinants of this progressive loss of fold specificity. We find that the alternate dimer-fold specifically results from the N11L substitution and is not promoted by other hydrophobic substitutions in the β-sheet. We also find that three highly hydrophobic substitutions at positions 9, 11, and 13 are necessary and sufficient for oligomer formation, but the oligomer size depends on the identity of the hydrophobic residue in question. The hydrophobic substitutions increase thermal stability, illustrating how increased hydrophobicity can increase folding stability even as it degrades conformational specificity. The oligomeric variants are predicted to be aggregation-prone but may be hindered from doing so by proline residues that flank the β-sheet region. Loss of conformational specificity due to increased hydrophobicity can manifest itself at any level of structure, depending upon the specific mutations and the context in which they occur.

Molecular insights of inhibition in sickle hemoglobin polymerization upon glutathionylation: hydrogen/deuterium exchange mass spectrometry and molecular dynamics simulation-based approach

In sickle cell anemia, polymerization of hemoglobin in its deoxy state leads to the formation of insoluble fibers that result in sickling of red blood cells. Stereo-specific binding of isopropyl group of βVal6, the mutated amino-acid residue of a tetrameric sickle hemoglobin molecule (HbS), with hydrophobic groove of another HbS tetramer initiates the polymerization. Glutathionylation of βCys93 in HbS was reported to inhibit the polymerization. However, the mechanism of inhibition in polymerization is unknown to date. In our study, the molecular insights of inhibition in polymerization were investigated by monitoring the conformational dynamics in solution phase using hydrogen/deuterium exchange-based mass spectrometry. The conformational rigidity imparted due to glutathionylation of HbS results in solvent shielding of βVal6 and perturbation in the conformation of hydrophobic groove of HbS. Additionally, molecular dynamics simulation trajectory showed that the stereo-specific localization of glutathione moiety in the hydrophobic groove across the globin subunit interface of tetrameric HbS might contribute to inhibition in polymerization. These conformational insights in the inhibition of HbS polymerization upon glutathionylation might be translated in the molecularly targeted therapeutic approaches for sickle cell anemia.

Targeting HbS Polymerization

The mutation of β6 from glu to val in hemoglobin is responsible for the polymer formation that leads to vaso-occlusion, and a range of severe consequences in sickle cell disease. The treatment of the disease can be addressed in many ways, but the prevention of polymer formation is one of the most fundamental approaches one can take. Such prevention includes affecting the polymer structure, or dilution of the fraction of polymerizable hemoglobin. The latter approach includes (1) induction of HbF, which does not itself, nor in hybrid form, join sickle polymers, or (2) restricting the allosteric change in hemoglobin that occurs in oxygen delivery, and which is required for polymer formation. These approaches will be critically reviewed, as well as the most recent developments that show the benefits of simply swelling the volume of the red cell.

Nitric oxide pathology and therapeutics in sickle cell disease

Sickle cell disease is caused by a mutant form of hemoglobin that polymerizes under hypoxic conditions which leads to red blood cell (RBC) distortion, calcium-influx mediated RBC dehydration, increased RBC adhesivity, reduced RBC deformability, increased RBC fragility, and hemolysis. These impairments in RBC structure and function result in multifaceted downstream pathology including inflammation, endothelial cell activation, platelet and leukocyte activation and adhesion, and thrombosis, all of which contribute vascular occlusion and substantial morbidity and mortality. Hemoglobin released upon RBC hemolysis scavenges nitric oxide (NO) and generates reactive oxygen species (ROS) and thereby decreases bioavailability of this important signaling molecule. As the endothelium-derived relaxing factor, NO acts as a vasodilator and also decreases platelet, leukocyte, and endothelial cell activation. Thus, low NO bioavailability contributes to pathology in sickle cell disease and its restoration could serve as an effective treatment. Despite its promise, clinical trials based on restoring NO bioavailability have so far been mainly disappointing. However, particular "NO donating" agents such as nitrite, which unlike some other NO donors can improve sickle RBC properties, may yet prove effective.

Wetting of nonconserved residue-backbones: A feature indicative of aggregation associated regions of proteins

Aggregation is an irreversible form of protein complexation and often toxic to cells. The process entails partial or major unfolding that is largely driven by hydration. We model the role of hydration in aggregation using "Dehydrons." "Dehydrons" are unsatisfied backbone hydrogen bonds in proteins that seek shielding from water molecules by associating with ligands or proteins. We find that the residues at aggregation interfaces have hydrated backbones, and in contrast to other forms of protein-protein interactions, are under less evolutionary pressure to be conserved. Combining evolutionary conservation of residues and extent of backbone hydration allows us to distinguish regions on proteins associated with aggregation (non-conserved dehydron-residues) from other interaction interfaces (conserved dehydron-residues). This novel feature can complement the existing strategies used to investigate protein aggregation/complexation.

Sickle Cell Hemoglobin with Mutation at αHis-50 Has Improved Solubility

The unliganded tetrameric Hb S has axial and lateral contacts with neighbors and can polymerize in solution. Novel recombinants of Hb S with single amino acid substitutions at the putative axial (recombinant Hb (rHb) (βE6V/αH20R) and rHb (βE6V/αH20Q)) or lateral (rHb (βE6V/αH50Q)) or double amino acid substitutions at both the putative axial and lateral (rHb (βE6V/αH20R/αH50Q) and rHb (βE6V/αH20Q/αH50Q)) contact sites were expressed in Escherichia coli and purified for structural and functional studies. The (1)H NMR spectra of the CO and deoxy forms of these mutants indicate that substitutions at either αHis-20 or αHis-50 do not change the subunit interfaces or the heme pockets of the proteins. The double mutants show only slight structural alteration in the β-heme pockets. All mutants have similar cooperativity (n50), alkaline Bohr effect, and autoxidation rate as Hb S. The oxygen binding affinity (P50) of the single mutants is comparable with that of Hb S. The double mutants bind oxygen with slightly higher affinity than Hb S under the acidic conditions. In high salt, rHb (βE6V/αH20R) is the only mutant that has a shorter delay time of polymerization and forms polymers more readily than Hb S with a dextran-Csat value of 1.86 ± 0.20 g/dl. Hb S, rHb (βE6V/αH20Q), rHb (βE6V/αH50Q), rHb (βE6V/αH20R/αH50Q), and rHb (βE6V/αH20Q/αH50Q) have dextran-Csat values of 2.95 ± 0.10, 3.04 ± 0.17, 11.78 ± 0.59, 7.11 ± 0.66, and 10.89 ± 0.83 g/dl, respectively. rHb (βE6V/αH20Q/αH50Q) is even more stable than Hb S under elevated temperature (60 °C).

Role of the beta4Thr-beta73Asp hydrogen bond in HbS polymer and domain formation from multinucleate-containing clusters

Fiber formation and domain formation from deoxy-HbS as well as from beta4 and beta73 HbS variants were investigated after temperature jump using DIC microscopy to gain a basic understanding of the determinants involved. Oversaturated deoxy-HbS generated numerous 14-stranded fibers and formed ovoid-shaped, multispherulitic domains. Domain number increased linearly as a function of time. Oversaturated deoxy-alpha2beta2(E6V,T4S) also generated time-dependent, ovoid-shaped spherulitic domains like HbS and alpha 2beta2(E6V,D73H) in the deoxy form. In contrast, alpha 2beta2(E6V,T4Y) and HbC-Harlem (alpha2beta2(E6V,D73N)) in the deoxy form generated time-dependent, ball-shaped domains containing many straight, crystalline-like fibers without evidence of branching. Some of these domains formed large needlelike crystals after overnight incubation. The inhibitory effect on polymer formation by beta4Tyr in HbS was stronger than that by beta4Ser but weaker than that by beta73Asn or beta73Leu. In contrast, both deoxy- and oxy-alpha2beta2(E6V,T4V) promoted formation of tiny, disordered amorphous aggregates without a delay time like oxy-HbS, which is in contrast to formation after a delay time of needlelike fibers for alpha 2beta2(E6V,D73L). Solubilities for both deoxy- and oxy-alpha 2beta2(E6V,T4V) were similar to that of deoxy-alpha 2beta2(E6V,D73H) but approximately 10-fold lower than that of deoxy-HbS. These results suggest that the strength of the hydrogen bond between beta4Thr and beta73Asp and the balance between the hydrogen bond and beta6Val hydrophobic interactions in deoxy-HbS polymers control formation of different types of fibers in a single domain or lead to formation of disordered, non-nucleated amorphous aggregates. These results also lead to a model in which multinucleation rather than a single-nucleation event occurs in a single cluster to generate numerous fibers growing from a single domain.

HbS-Savaria: the anti-polymerization effect of a single mutation in human alpha-chains

Recombinant alpha-Savaria globin (alpha(S49R)) was assembled with beta(S) chains by the alloplex intermediate pathway to generate tetrameric rHbS-Sarvaria (alpha (2) (S49R) beta (2) (E6V) ) that exhibited normal O(2) affinity and co-operatively at pH 7.4. Allosteric effectors, 2,3-DPG, L35, and NaCl increased O(2) affinity by 15%. Bohr effects were similar for rHbS-Savaria and HbS (0.38 +/- 0.025 vs. 0.46 +/- 0.03, respectively). The C(SAT) of HbS increased from 16.7 +/- 0.8 to 27.0 +/- 1.0 g/dL. Co-polymerization demonstrated inhibition predominantly by the Cis-dimer. Molecular modeling indicated that the positive charge at alpha-49 generated a strong anion-binding site and reduced flexibility of the CD-region by restricting movement in the E and F helices. The molecular distance between Arg-49 and Asn-78 in the neighboring double strand decreased, and electrostatic repulsion between the inter-double strands increased, resulting in inhibition of polymerization. The Savaria mutation may be useful for the design of super-inhibitory alpha-chains and gene therapy of sickle cell anemia.

The role of beta93 Cys in the inhibition of Hb S fiber formation

Recent studies have suggested that nitric oxide (NO) binding to hemoglobin (Hb) may lead to the inhibition of sickle cell fiber formation and the dissolution of sickle cell fibers. NO can react with Hb in at least 3 ways: 1) formation of Hb(II)NO, 2) formation of methemoglobin, and 3) formation of S-nitrosohemoglobin, through nitrosylation of the beta93 Cys residue. In this study, the role of beta93 Cys in the mechanism of sickle cell fiber inhibition is investigated through chemical modification with N-ethylmaleimide. UV resonance Raman, FT-IR and electrospray ionization mass spectroscopic methods in conjunction with equilibrium solubility and kinetic studies are used to characterize the effect of beta93 Cys modification on Hb S fiber formation. Both FT-IR spectroscopy and electrospray mass spectrometry results demonstrate that modification can occur at both the beta93 and alpha104 Cys residues under relatively mild reaction conditions. Equilibrium solubility measurements reveal that singly-modified Hb at the beta93 position leads to increased amounts of fiber formation relative to unmodified or doubly-modified Hb S. Kinetic studies confirm that modification of only the beta93 residue leads to a faster onset of polymerization. UV resonance Raman results indicate that modification of the alpha104 residue in addition to the beta93 residue significantly perturbs the alpha(1)beta(2) interface, while modification of only beta93 does not. These results in conjunction with the equilibrium solubility and kinetic measurements are suggestive that modification of the alpha104 Cys residue and not the beta93 Cys residue leads to T-state destabilization and inhibition of fiber formation. These findings have implications for understanding the mechanism of NO binding to Hb and NO inhibition of Hb S fiber formation.

The Hb A variant (beta73 Asp-->Leu) disrupts Hb S polymerization by a novel mechanism

Polymerization of a 1:1 mixture of hemoglobin S (Hb S) and the artificial mutant HbAbeta73Leu produces a dramatic morphological change in the polymer domains in 1.0 M phosphate buffer that are a characteristic feature of polymer formation. Instead of feathery domains with quasi 2-fold symmetry that characterize polymerization of Hb S and all previously known mixtures such as Hb A/S and Hb F/S mixtures, these domains are compact structures of quasi-spherical symmetry. Solubility of Hb S/Abeta73Leu mixtures was similar to that of Hb S/F mixtures. Kinetics of polymerization indicated that homogeneous nucleation rates of Hb S/Abeta73Leu mixtures were the same as those of Hb S/F mixtures, while exponential polymer growth (B) of Hb S/Abeta73Leu mixtures were about three times slower than those of Hb S/F mixtures. Differential interference contrast (DIC) image analysis also showed that fibers in the mixture appear to elongate between three and five times more slowly than in equivalent Hb S/F mixtures by direct measurements of exponential growth of mass of polymer in a domain. We propose that these results of Hb S/Abeta73Leu mixtures arise from a non-productive binding of the hybrid species of this mixture to the end of the growing polymer. This "cap" prohibits growth of polymers, but by nature is temporary, so that the net effect is a lowered growth rate of polymers. Such a cap is consistent with known features of the structure of the Hb S polymer. Domains would be more spherulitic because slower growth provides more opportunity for fiber bending to spread domains from their initial 2-fold symmetry. Moreover, since monomer depletion proceeds more slowly in this mixture, more homogeneous nucleation events occur, and the resulting gel has a far more granular character than normally seen in mixtures of non-polymerizing hemoglobins with Hb S. This mixture is likely to be less stiff than polymerized mixtures of other hybrids such as Hb S with HbF, potentially providing a novel approach to therapy.

Anisotropy in Sickle Hemoglobin Fibers from Variations in Bending and Twist

We have studied the variations of twist and bend in sickle hemoglobin fibers. We find that these variations are consistent with an origin in equilibrium thermal fluctuations, which allows us to estimate the bending and torsional rigidities and effective corresponding material moduli. We measure bending by electron microscopy of frozen hydrated fibers and find that the bending persistence length, a measure of the length of fiber required before it starts to be significantly bent due to thermal fluctuations, is 130microm, somewhat shorter than that previously reported using light microscopy. The torsional persistence length, obtained by re-analysis of previously published experiments, is found to be only 2.5microm. Strikingly this means that the corresponding torsional rigidity of the fibers is only 6x10(-27)Jm, much less than their bending rigidity of 5x10(-25)Jm. For (normal) isotropic materials, one would instead expect these to be similar. Thus, we present the first quantitative evidence of a very significant material anisotropy in sickle hemoglobin fibers, as might arise from the difference between axial and lateral contacts within the fiber. We suggest that the relative softness of the fiber with respect to twist deformation contributes to the metastability of HbS fibers: HbS double strands are twisted in the fiber but not in the equilibrium crystalline state. Our measurements inform a theoretical model of the thermodynamic stability of fibers that takes account of both bending and extension/compression of hemoglobin (double) strands within the fiber.

Structure of the thermolabile mutant aldolase B, A149P: molecular basis of hereditary fructose intolerance

Hereditary fructose intolerance (HFI) is a potentially lethal inborn error in metabolism caused by mutations in the aldolase B gene, which is critical for gluconeogenesis and fructose metabolism. The most common mutation, which accounts for 53% of HFI alleles identified worldwide, results in substitution of Pro for Ala at position 149. Structural and functional investigations of human aldolase B with the A149P substitution (AP-aldolase) have shown that the mutation leads to losses in thermal stability, quaternary structure, and activity. X-ray crystallography is used to reveal the structural basis of these perturbations. Crystals of AP-aldolase are grown at two temperatures (4 degrees C and 18 degrees C), and the structure solved to 3.0 angstroms resolution, using the wild-type structure as the phasing model. The structures reveal that the single residue substitution, A149P, causes molecular disorder around the site of mutation (residues 148-159), which is propagated to three adjacent beta-strand and loop regions (residues 110-129, 189-199, 235-242). Disorder in the 110-129-loop region, which comprises one subunit-subunit interface, provides an explanation for the disrupted quaternary structure and thermal instability. Greater structural perturbation, particularly at a Glu189-Arg148 salt bridge in the active-site architecture, is observed in the structure determined at 18 degrees C, which could explain the temperature-dependent loss in activity. The disorder revealed in these structures is far greater than that predicted by homology modeling and underscores the difficulties in predicting perturbations of protein structure and function by homology modeling alone. The AP-aldolase structure reveals the molecular basis of a hereditary disease and represents one of only a few structures known for mutant proteins at the root of the thousands of other inherited disorders.

Linkage of interactions in sickle hemoglobin fiber assembly: inhibitory effect emanating from mutations in the AB region of the alpha-chain is annulled by a mutation at its EF corner

The AB and GH regions of the alpha-chain are located in spatial proximity and contain a cluster of intermolecular contact residues of the sickle hemoglobin (HbS) fiber. We have examined the role of dynamics of AB/GH region on HbS polymerization through simultaneous replacement of non-contact Ala(19) and Ala(21) of the AB corner with more flexible Gly or rigid alpha-aminoisobutyric acid (Aib) residues. The polymerization behavior of HbS with Aib substitutions was similar to the native HbS. In contrast, Gly substitutions inhibited HbS polymerization. Molecular dynamics simulation studies of alpha-chains indicated that coordinated motion of AB and GH region residues present in native (Ala) as well as in Aib mutant was disrupted in the Gly mutant. The inhibitory effect due to Gly substitutions was further explored in triple mutants that included mutation of an inter-doublestrand contact (alphaAsn(78) --> His or Gln) at the EF corner. Although the inhibitory effect of Gly substitutions in the triple mutant was unaffected in the presence of alphaGln(78), His at this site almost abrogated its inhibitory potential. The polymerization studies of point mutants (alphaGln(78) --> His) indicated that the inhibitory effect due to Gly substitutions in the triple mutant was synergistically compensated for by the polymerization-enhancing activity of His(78). Similar synergistic coupling, between alphaHis(78) and an intra-double-strand contact point (alpha16) mutation located in the AB region, was also observed. Thus, two conclusions are made: (i) Gly mutations at the AB corner inhibit HbS polymerization by perturbing the dynamics of the AB/GH region, and (ii) perturbations of AB region (through changes in dynamics of the AB/GH region or abolition of a specific fiber contact site) that influence HbS polymerization do so in concert with alpha78 site at the EF corner. The overall results provide insights about the interaction-linkage between distant regions of the HbS tetramer in fiber assembly.

Twisted protein aggregates and disease: the stability of sickle hemoglobin fibers

We describe how twist could play an essential role in stabilizing 20 nm diameter sickle hemoglobin fibers. Our theory successfully reproduces the observed variation of helical pitch length with fiber diameter. With no remaining adjustable parameters it also yields a prediction for the torsional rigidity of sickle hemoglobin fibers that is in good agreement with experiment and hence retains the striking feature that such fibers can be highly mechanically anisotropic, even with a ratio of bending to torsional rigidity of about 50. We discuss how our study might be relevant to the development of treatment strategies.

Sickle hemoglobin polymer stability probed by triple and quadruple mutant hybrids

As part of an effort to understand the interactions in HbS polymerization, we have produced and studied a recombinant triple mutant, D6A(alpha)/D75Y(alpha)/E121R(beta), and a quadruple mutant comprising the preceding mutation plus the natural genetic mutation of sickle hemoglobin, E6V(beta). These recombinant hemoglobins expressed in yeast were extensively characterized, and their structure and oxygen binding cooperativity were found to be normal. Their tetramer-dimer dissociation constants were within a factor of 2 of HbA and HbS. Polymerization of these mutants mixed with HbS was investigated by a micromethod based on volume exclusion by dextran. The elevated solubility of mixtures of HbS with HbA and HbF in dextran could be accurately predicted without any variable parameters. Relative to HbS, the copolymerization probability of the quadruple mutant/HbS hybrid was found to be 6.2, and the copolymerization probability for the triple mutant/HbS hybrid was 0.52. The pure quadruple mutant had a solubility slightly above that of its hybrid with HbS. One way to explain these results is to require significant cis-trans differences in the polymer and that HbA assemble above 42.5 g/dl. A second way to explain these data is by the modification of motional freedom, thereby changing vibrational entropy in the polymer.

A role for the alpha 113 (GH1) amino acid residue in the polymerization of sickle hemoglobin. Evaluation of its inhibitory strength and interaction linkage with two fiber contact sites (alpha 16/23) located in the AB region of the alpha-chain

A cluster of amino acid residues located in the AB-GH region of the alpha-chain are shown in intra-double strand axial interactions of the hemoglobin S (HbS) polymer. However, alphaLeu-113 (GH1) located in the periphery is not implicated in any interactions by either crystal structure or models of the fiber, and its role in HbS polymerization has not been explored by solution experiments. We have constructed HbS Twin Peaks (betaGlu-6-->Val, alphaLeu-113-->His) to ascertain the hitherto unknown role of the alpha113 site in the polymerization process. The structural and functional behavior of HbS Twin Peaks was comparable with HbS. HbS Twin Peaks polymerized with a slower rate compared with HbS, and its polymer solubility (C(sat)) was found to be about 1.8-fold higher than HbS. To further authenticate the participation of the alpha113 site in the polymerization process as well as to evaluate its relative inhibitory strength, we constructed HbS tetramers in which the alpha113 mutation was coupled individually with two established fiber contact sites (alpha16 and alpha23) located in the AB region of the alpha-chain: HbS(alphaLys-16-->Gln, alphaLeu-113-->His), HbS(alphaGlu-23-->Gln, alphaLeu-113-->His). The single mutants at alpha16/alpha23 sites were also engineered as controls. The C(sat) values of the HbS point mutants involving sites alpha16 or alpha23 were higher than HbS but markedly lower as compared with HbS Twin Peaks. In contrast, C(sat) values of both double mutants were comparable with or higher than that of HbS Twin Peaks. The demonstration of the inhibitory effect of alpha113 mutation alone or in combination with other sites, in quantitative terms, unequivocally establishes a role for this site in HbS gelation. These results have implications for development of a more accurate model of the fiber that could serve as a blueprint for therapeutic intervention.

Polymerization of deoxy-sickle cell hemoglobin in high-phosphate buffer

Deoxy-sicklecell hemoglobin (HbS) polymerizes in 0.05 M phosphate buffer to form long helical fibers. The reaction typically occurs when the concentration of HbS is about 165 mg/ml. Polymerization produces a variety of polymorphic forms. The structure of the fibers can be probed by using site-directed mutants to examine the effect of altering the residues involved in intermolecular interactions. Polymerization can also be induced in the presence of 1.5 M phosphate buffer. Under these conditions polymerization occurs at much lower concentrations (ca. 5 mg/ml), which is advantageous when site-directed mutants are being used because only small quantities of the mutants are available. We have characterized the structure of HbS polymers formed in 1.5 M phosphate to determine how their structures are related to the polymers formed under more physiological conditions. Under both sets of conditions fibers are the first species to form. At pHs between 6.7 and 7.3 fibers initially form bundles and then crystals. At lower pHs fibers form macrofibers and then crystals. Fourier transforms of micrographs of the polymers formed in 1.5 M phosphate display the 32- and 64-A(-1) periodicity characteristic of fibers formed in 0.05 M phosphate buffer. The 64-A(-1) layer line is less prominent in Fourier transforms of negatively stained fibers formed in 1.5 M phosphate possibly because salt interferes with staining of the fibers. However, micrographs and Fourier transforms of frozen hydrated fibers formed in high and low phosphate display the same periodicities. Under both sets of reaction conditions HbS polymers form crystals with the same unit cell parameters as Wishner-Love crystals (a = 64 A, b = 185 A, c = 53 A). Some of the polymerization intermediates were examined in the frozen-hydrated state in order to determine whether their structures were significantly perturbed by negative staining. We have also carried out reconstructions of the frozen-hydrated fibers in high and low phosphate to compare their molecular coordinates. The helical projection of the reconstructions in low phosphate shows the expected 14-strand structure. In high phosphate the 14-strand fibers are also formed and their molecular coordinates are the same (within experimental error) as those of fibers formed in 0.05 M phosphate. In addition, the reconstructions of high-phosphate fibers reveal a new minor variant of fiber containing 10 strands. The polymerization products in 1.5 M phosphate buffer were generally indistinguishable from those formed in 0.05 M phosphate buffer. Micrographs of frozen hydrated specimens have facilitated the interpretation of previously published micrographs using negative staining.

A model for the sickle hemoglobin fiber using both mutation sites

The standard molecular model of the fiber of the sickle hemoglobin (HbS: beta6 Glu-->Val) has been revised to allow both beta6 mutation sites to participate in intermolecular contacts, rather than only one beta6 site as previously thought, for four molecules per 14-molecule fiber cross section. This structure accurately predicts the copolymerization of hybridized mixtures of HbS with HbA or HbC (beta6 Glu-->Lys), which could not be reconciled with prior models in which only half the beta6 sites were required for assembly. This model suggests new contacts within the fiber and raises the question of whether these cross-linked double strands could possess added stability important in such processes as nucleation.

Solution-active structural alterations in liganded hemoglobins C (beta6 Glu --> Lys) and S (beta6 Glu --> Val)

Based upon existing crystallographic evidence, HbS, HbC, and HbA have essentially the same molecular structure. However, important areas of the molecule are not well defined crystallographically (e.g. the N-terminal nonhelical portion of the alpha and beta chains), and conformational constraints differ in solution and in the crystalline state. Over the years, our laboratory and others have provided evidence of conformational changes in HbS and, more recently, in HbC. We now present data based upon allosteric perturbation monitored by front-face fluorescence, ultraviolet resonance Raman spectroscopy, circular dichroism, and oxygen equilibrium studies that confirm and significantly expand previous findings suggesting solution-active structural differences in liganded forms of HbS and HbC distal to the site of mutation and involving the 2,3-diphosphoglycerate binding pocket. The liganded forms of these hemoglobins are of significant interest because HbC crystallizes in the erythrocyte in the oxy form, and oxy HbS exhibits increased mechanical precipitability and a high propensity to oxidize. Specific findings are as follows: 1) differences in the intrinsic fluorescence indicate that the Trp microenvironments are more hydrophobic for HbS > HbC > HbA, 2) ultraviolet resonance Raman spectroscopy detects alterations in Tyr hydrogen bonding, in Trp hydrophobicity at the alpha1beta2 interface (beta37), and in the A-helix (alpha14/beta15) of both chains, 3) displacement by inositol hexaphosphate of the Hb-bound 8-hydroxy-1,3,6-pyrenetrisulfonate (the fluorescent 2,3-diphosphoglycerate analog) follows the order HbA > HbS > HbC, and 4) oxygen equilibria measurements indicate a differential allosteric effect by inositol hexaphosphate for HbC approximately HbS > HbA.

Analysis of the stability of hemoglobin S double strands

The deoxyhemoglobin S (deoxy-HbS) double strand is the fundamental building block of both the crystals of deoxy-HbS and the physiologically relevant fibers present within sickle cells. To use the atomic-resolution detail of the hemoglobin-hemoglobin interaction known from the crystallography of HbS as a basis for understanding the interactions in the fibers, it is necessary to define precisely the relationship between the straight double strands in the crystal and the twisted, helical double strands in the fibers. The intermolecular contact conferring the stability of the double strand in both crystal and fiber is between the beta6 valine on one HbS molecule and residues near the EF corner of an adjacent molecule. Models for the helical double strands were constructed by a geometric transformation from crystal to fiber that preserves this critical interaction, minimizes distortion, and makes the transformation as smooth as possible. From these models, the energy of association was calculated over the range of all possible helical twists of the double strands and all possible distances of the double strands from the fiber axis. The calculated association energies reflect the fact that the axial interactions decrease as the distance between the double strand and the fiber axis increases, because of the increased length of the helical path taken by the double strand. The lateral interactions between HbS molecules in a double strand change relatively little between the crystal and possible helical double strands. If the twist of the fiber or the distance between the double strand and the fiber axis is too great, the lateral interaction is broken by intermolecular contacts in the region around the beta6 valine. Consequently, the geometry of the beta6 valine interaction and the residues surrounding it severely restricts the possible helical twist, radius, and handedness of helical aggregates constructed from the double strands. The limitations defined by this analysis establish the structural basis for the right-handed twist observed in HbS fibers and demonstrates that for a subunit twist of 8 degrees, the fiber diameter cannot be more than approximately 300 A, consistent with electron microscope observations. The energy of interaction among HbS molecules in a double strand is very slowly varying with helical pitch, explaining the variable pitch observed in HbS fibers. The analysis results in a model for the HbS double strand, for use in the analysis of interactions between double strands and for refinement of models of the HbS fibers against x-ray diffraction data.

Systematic Enhancement of Polymerization of Recombinant Sickle Hemoglobin Mutants: Implications for Transgenic Mouse Model for Sickle Cell Anemia

To provide quantitative information on the sites that promote polymerization of sickle hemoglobin (HbS) after formation of the initial hydrophobic bond involving Val-6(β) [E6V(β)] and also to provide hemoglobins with an enhanced polymerization that could be used in a mouse model for sickle cell anemia, we have expressed recombinant double, triple, and quadruple HbS mutants with substitutions on both the α- and β-chains, E6V(β)/E121R(β), D75Y(α)/E6V(β)/E121R(β) and D6A(α)/D75Y(α)/E6V(β)/E121R(β). These recombinant hemoglobins were extensively characterized by high-performance liquid chromatography analysis, sodium dodecyl sulfate-polyacrylamide gel electrophoresis, isoelectric focusing, amino acid analysis, and mass spectroscopy. They retained the functional properties of the Hb tetramer and polymerized in a linear manner at progressively lower Hb concentration as a function of the degree of substitution, suggesting that these remote sites (αD6A, αD75Y, and βE121R) on the α- and β-chains exhibit additive, enhanced polymerization properties. The quadruple mutant has a polymerization concentration close to that of the purified SAD hemoglobin from transgenic mouse red blood cells consisting of HbS, Hb Antilles, and Hb D-Punjab. Normal mouse Hb increases the polymerization concentration of each mutant. Thus, the general approach of using recombinant Hbs as described here should prove useful in elucidating the quantitative aspects of the mechanism of HbS polymerization and in identifying the contribution of individual sites to the overall process. The strategy described here demonstrates the feasibility of a systematic approach to achieve future recombinant HbS mutants that could provide a new generation of the transgenic mouse model for sickle cell anemia.

Systematic Enhancement of Polymerization of Recombinant Sickle Hemoglobin Mutants: Implications for Transgenic Mouse Model for Sickle Cell Anemia

To provide quantitative information on the sites that promote polymerization of sickle hemoglobin (HbS) after formation of the initial hydrophobic bond involving Val-6(β) [E6V(β)] and also to provide hemoglobins with an enhanced polymerization that could be used in a mouse model for sickle cell anemia, we have expressed recombinant double, triple, and quadruple HbS mutants with substitutions on both the α- and β-chains, E6V(β)/E121R(β), D75Y(α)/E6V(β)/E121R(β) and D6A(α)/D75Y(α)/E6V(β)/E121R(β). These recombinant hemoglobins were extensively characterized by high-performance liquid chromatography analysis, sodium dodecyl sulfate-polyacrylamide gel electrophoresis, isoelectric focusing, amino acid analysis, and mass spectroscopy. They retained the functional properties of the Hb tetramer and polymerized in a linear manner at progressively lower Hb concentration as a function of the degree of substitution, suggesting that these remote sites (αD6A, αD75Y, and βE121R) on the α- and β-chains exhibit additive, enhanced polymerization properties. The quadruple mutant has a polymerization concentration close to that of the purified SAD hemoglobin from transgenic mouse red blood cells consisting of HbS, Hb Antilles, and Hb D-Punjab. Normal mouse Hb increases the polymerization concentration of each mutant. Thus, the general approach of using recombinant Hbs as described here should prove useful in elucidating the quantitative aspects of the mechanism of HbS polymerization and in identifying the contribution of individual sites to the overall process. The strategy described here demonstrates the feasibility of a systematic approach to achieve future recombinant HbS mutants that could provide a new generation of the transgenic mouse model for sickle cell anemia.

Inhibition of sickle beta-chain (betaS)-dependent polymerization by nonhuman alpha-chains. A superinhibitory mouse-horse chimeric alpha-chain

Horse alpha-chain inhibits sickle beta-chain-dependent polymerization; however, its inhibitory potential is not as high as that of mouse alpha-chain. Horse alpha-(1-30) and alpha-(31-141) segments make, respectively, minor and major contributions to the inhibitory potential of horse alpha-chain. The sum of the inhibitory potential of the two segments does not account for the inhibitory potential of the full-length horse alpha-chain. Although the polymerization inhibitory potential of horse alpha-chain is lower than mouse alpha-chain, the inhibitory potential of horse alpha-(31-141) is comparable to that of mouse alpha-(31-141). When mouse alpha-(1-30) is stitched to horse alpha-(31-141), the product is a chimeric alpha-chain with an inhibitory potential greater than mouse alpha-chain. In contrast, the stitching of horse alpha-(1-30) with mouse alpha-(31-141) had no additional inhibitory potential. Molecular modeling studies of HbS containing the mouse-horse chimeric alpha-chain indicate altered side-chain interactions at the alpha1beta1 interface when compared with HbS. In addition, the AB/GH corner perturbations facilitate a different stereochemistry for the interaction of the epsilon-amino group of Lys-16(alpha) with the beta-carboxyl group of Asp-116(alpha), resulting in a decrease in the accessibility of the side chain of Lys-16(alpha) to the solvent. Based on molecular modeling, we speculate that these perturbations by themselves, or in synergy with the altered conformational aspects of the alpha1beta1 interactions, represent the molecular basis of the superinhibitory potential of the mouse-horse chimeric alpha-chains.

The high resolution crystal structure of deoxyhemoglobin S

We have refined the crystal structure of deoxyhemoglobin S (beta Glu6-->Val) at 2.05 A resolution to an R-factor of 16.5% (free R=21. 5%) using crystals isomorphous to those originally grown by Wishner and Love. A predominant feature of this crystal form is a double strand of hemoglobin tetramers that has been shown by a variety of techniques to be the fundamental building block of the intracellular sickle cell fiber. The double strand is stabilized by lateral contacts involving the mutant valine interacting with a pocket between the E and F helices on another tetramer. The new structure reveals some marked differences from the previously refined 3.0 A resolution structure, including several residues in the lateral contact which have shifted by as much as 3.5 A. The lateral contact includes, in addition to the hydrophobic interactions involving the mutant valine, hydrophilic interactions and bridging water molecules at the periphery of the contact. This structure provides further insights into hemoglobin polymerization and may be useful for the structure-based design of therapeutic agents to treat sickle cell disease.

Influence of the A helix structure on the polymerization of hemoglobin S

Hb S variants containing Lys-beta132 --> Ala or Asn substitutions were engineered to evaluate the consequences of the A helix destabilization in the polymerization process. Previous studies suggested that the loss of the Glu-beta7-Lys-beta132 salt bridge in the recombinant Hb betaE6V/E7A could be responsible for the destabilization of the A helix. The recombinant Hb (rHb) S/beta132 variants polymerized with an increased delay time as well as decreased maximum absorbance and Hb solubility values similar to that of Hb S. These data indicate that the strength of the donor-acceptor site interaction may be reduced due to an altered conformation of the A helix. The question arises whether this alteration leads to a true inhibition of the polymerization process or to qualitatively different polymers. The oxygen affinity of the beta132 mutated rHbs was similar to that of Hb A and S, whereas the cooperativity and effects of organic phosphates were reduced. This could be attributed to modifications in the central cavity due to loss of the positively charged lysine. Since Lys-beta132 is involved in the stabilization of the alpha1-beta1 interface, the loss of the beta132(H10)-beta128(H6) salt bridge may be responsible for the marked thermal instability of the beta132 mutated rHbs.

Recombinant Sickle Hemoglobin Containing a Lysine Substitution at Asp-85(α): Expression in Yeast, Functional Properties, and Participation in Gel Formation

Clinical modalities based on inhibition of gelation of HbS are hindered by the lack of quantitative information on the extent of participation of different amino acid residues in the aggregation process. One such site is Asp-85(α), which is involved in a parallel interdouble strand ionic interaction with Lys-144(β) according to the crystal structure of HbS, but electron microscopy does not specifically show Asp-85(α) as a contact site for fiber formation. Using a yeast recombinant system, we have substituted this site by Lys to abolish ion pairing and to make a quantitative determination of its participation in aggregation. The purified double mutant was shown to have the expected pI, the calculated molecular weight, correct amino acid composition, and peptide map. The recombinant double mutant has an oxygen affinity of 10 mm Hg, which is identical to that for HbA and HbS under the same conditions; it also has high cooperativity with an average n value of 2.7. The change in P50 in response to chloride ions was about 25% less than that for HbA or HbS and is ascribed to the introduction of a new positive charge near one of the major oxygen-linked chloride binding sites of hemoglobin. The gelation concentration of the double mutant was measured by a new procedure (Bookchin et al, 1994); the maximal amount of soluble hemoglobin (Csat ) in the presence of dextran indicated a decreased tendency for gelation with a Csat of 53 mg/mL compared with 34 mg/mL for HbS. This inhibitory effect is smaller than that of the E6V(β)/L88A(β) (Csat , 67 mg/mL) and the E6V(β)/K95I(β) (Csat , 90 mg/mL) recombinant hemoglobins. Thus, we would classify Asp-85(α) as a moderate contributor to the strength of the HbS aggregate. This wide range of gelation values demonstrates that some sites are more important than others in promoting HbS aggregation.

Free energy simulations of axial contacts in sickle-cell hemoglobin

Molecular dynamics simulations have been used to investigate the thermodynamic stability of axial contacts in sickle-cell hemoglobin (HbS). Free energy changes were evaluated for the point mutation beta 121 Glu --> Gln in the axial contact region of HbS crystals. The calculations predict a free energy change of-3.6 kcal/mol per contact for the mutation, which is in qualitative agreement with experimental observations of aggravated sickling found in the double mutant Hb D Los Angeles (beta 6 Glu --> Val. beta 121 Glu --> Gln) relative to HbS (beta 6 Glu --> Val). The beta 121 Glu is sequestered in a salt link with beta 17 Lys located on the same polypeptide chain, making the Glu interactions with its surroundings similar in aggregates and individual hemoglobins. Due to this cancellation of the large electrostatic Glu contributions, the weak nonspecific interactions between the Gln and the neighboring polypeptide chain are the main contributing factor to the enhanced aggregation of Hb D Los Angeles relative to HbS. Together with the previous study of the lateral contact [K. Kuczera et al. (1990) Proceedings of the National Academy of Science USA, Vol. 87, pp, 8481-8485], the present results provide a more complete picture of the forces driving the sickling aggregation. A comparison of different treatments of internal flexibility in free energy simulations and analysis of rate of convergence of the different calculated properties has also been performed.

Functional studies and polymerization of recombinant hemoglobin Glu-alpha2beta26(A3) --> Val/Glu-7(A4) --> Ala

In hemoglobin (Hb) S the hydrophobic mutated residue Val-beta6(A3) (donor site) closely interacts with the hydrophobic side groups of Phe-beta85(F1) and Leu-beta88(F4) (EF pocket, acceptor site) of a neighboring tetramer, resulting in decreased solubility and polymerization of the deoxy-Hb. The beta6(A3) residue is followed by two charged residues Glu-beta7(A4) and Lys-beta8(A5). This cluster has no attraction for the hydrophobic EF pocket. We have modified the beta7(A4) residue next to the donor site Val-beta6(A3), replacing the charged Glu by a hydrophobic Ala-(rHb betaE6V/E7A). The single mutant Glu-beta7 --> Ala-(rHb betaE7A) was also engineered. Both rHbs exhibit a heat instability and an increased oxygen affinity compared to Hb A and Hb S. There was a concentration dependence of the ligand binding properties (1-300 microM in heme) indicating an increased amount of dimers relative to Hb A. The deoxy form of rHb betaE6V/E7A polymerizes in vitro, with a decreased rate of polymer formation relative to Hb S, while the single mutant betaE7A does not polymerize in the same experimental conditions. The Glu-beta7(A4) --> Ala substitution does not increase the hydrophobic interaction between donor and acceptor site. We speculate that the loss of the normal saline bridge between Glu-beta7(A4) and Lys-beta132(H10) leads to an increased flexibility of the A helix and may account for the difference of the polymerization for this Hb S mutant.

Significance of structural changes in proteins: expected errors in refined protein structures

A quantitative expression key to evaluating significant structural differences or induced shifts between any two protein structures is derived. Because crystallography leads to reports of a single (or sometimes dual) position for each atom, the significance of any structural change based on comparison of two structures depends critically on knowing the expected precision of each median atomic position reported, and on extracting it for each atom, from the information provided in the Protein Data Bank and in the publication. The differences between structures of protein molecules that should be identical, and that are normally distributed, indicating that they are not affected by crystal contacts, were analyzed with respect to many potential indicators of structure precision, so as to extract, essentially by "machine learning" principles, a generally applicable expression involving the highest correlates. Eighteen refined crystal structures from the Protein Data Bank, in which there are multiple molecules in the crystallographic asymmetric unit, were selected and compared. The thermal B factor, the connectivity of the atom, and the ratio of the number of reflections to the number of atoms used in refinement correlate best with the magnitude of the positional differences between regions of the structures that otherwise would be expected to be the same. These results are embodied in a six-parameter equation that can be applied to any crystallographically refined structure to estimate the expected uncertainty in position of each atom. Structure change in a macromolecule can thus be referenced to the expected uncertainty in atomic position as reflected in the variance between otherwise identical structures with the observed values of correlated parameters.

Mutational analysis of phenylalanine beta 85 in the valine beta 6 acceptor pocket during hemoglobin S polymerization

Hemoglobin (Hb) S containing Leu, Ala, Thr, or Trp substitutions at beta 85 were made and expressed in yeast in an effort to evaluate the role of Phe-beta 85 in the acceptor pocket during polymerization of deoxy Hb S. The four Hb S variants have the same electrophoretic mobility as Hb S, and these beta 85 substitutions do not significantly affect heme-globin interactions and tetramer helix content. Hb S containing Trp-beta 85 had decreased oxygen affinity, whereas those with Leu-, Ala-, and Thr-beta 85 had increased oxygen affinity. All four supersaturated beta 85 variants polymerized with a delay time as does deoxy Hb S. This is in contrast to deoxy Hb S containing Phe-beta 88, Ala-beta 88, Glu-beta 88, or Glu-beta 85, which polymerized with no clear delay time (Adachi K, Konitzer P, Paulraj CG, Surrey S, 1994, J Biol Chem 269:17477-17480; Adachi K, Reddy LR, Surrey S, 1994, J Biol Chem 269:31563-31566). Leu substitution at beta 85 accelerated deoxy Hb S polymerization, whereas Ala, Thr, or Trp substitution inhibited polymerization. The length of the delay time and total polymer formed for these beta 85 Hb S variants depended on hemoglobin concentration in the same fashion as for deoxy Hb S: the higher the concentration, the shorter the delay time and the more polymer formed. Critical concentrations required for polymerization of deoxy Hb SF veta 85L, Hb SF beta 85A, Hb SF beta 85T, and Hb SF beta 85W are 0.65-, 2.2-, 2.5- and 3-fold higher, respectively, than Hb S.(ABSTRACT TRUNCATED AT 250 WORDS)

Participation and strength of interaction of lysine 95(beta) in the polymerization of hemoglobin S as determined by its site-directed substitution by isoleucine

The role of Lys-95(beta), which is on the exterior of the hemoglobin (HbS) tetramer, in the aggregation process has been addressed because there is a lack of agreement on its importance. The early studies on the aggregation of HbS in the presence of other mutant hemoglobins are consistent with the subsequent electron microscopic studies in demonstrating the participation of Lys-95(beta) in gelation; the results of the crystal structure do not agree with these conclusions. Therefore, with the objective of clarifying its role we have carried out site-directed substitution of Lys-95(beta) to an isoleucine residue. The mutation was introduced by polymerase chain reaction recombination methodology, and the absence of other mutations in the beta-globin gene was established by sequencing the gene in its entirety. The recombinant mutant hemoglobin was expressed in yeast and characterized by peptide mapping and sequencing, which demonstrated that the only different tryptic peptide had the Ile substitution at position 95(beta). The recombinant hemoglobin had the correct amino acid composition and molecular weight by mass spectrometric analysis. It was also pure as judged by isoelectric focusing. It was fully functional because it had an average Hill coefficient of 3.1 and responded normally to the allosteric regulators, chloride, 2,3-diphosphoglycerate, and inositol hexaphosphate. Of particular interest was the finding that this hemoglobin mutant aggregated at a concentration of about 40 g/dl, nearly twice that at which HbS itself aggregated (24 g/dl). Therefore, Lys-95(beta) has a very important role in the aggregation process and is a good candidate site for the design of a therapeutic agent for sickle cell anemia.

Recognition of distantly related proteins through energy calculations

A new method to detect remote relationships between protein sequences and known three-dimensional structures based on direct energy calculations and without reliance on statistics has been developed. The likelihood of a residue to occupy a given position on the structural template was represented by an estimate of the stabilization free energy made after explicit prediction of the substituted side chain conformation. The profile matrix derived from these energy values and modified by increasing the residue self-exchange values successfully predicted compatibility of heat-shock protein and globin sequences with the three-dimensional structures of actin and phycocyanin, respectively, from a full protein sequence databank search. The high sensitivity of the method makes it a unique tool for predicting the three-dimensional fold for the rapidly growing number of protein sequences.