Functional and subunit assembly properties of hemoglobin Alberta (alpha 2 beta 2(101) Glu----Gly)

Computationally accessible method for estimating free energy changes resulting from site-specific mutations of biomolecules: systematic model building and structural/hydropathic analysis of deoxy and oxy hemoglobins

A practical computational method for the molecular modeling of free-energy changes associated with protein mutations is reported. The de novo generation, optimization, and thermodynamic analysis of a wide variety of deoxy and oxy hemoglobin mutants are described in detail. Hemoglobin is shown to be an ideal candidate protein for study because both the native deoxy and oxy states have been crystallographically determined, and a large and diverse population of its mutants has been thermodynamically characterized. Noncovalent interactions for all computationally generated hemoglobin mutants are quantitatively examined with the molecular modeling program HINT (Hydropathic INTeractions). HINT scores all biomolecular noncovalent interactions, including hydrogen bonding, acid-base, hydrophobic-hydrophobic, acid-acid, base-base, and hydrophobic-polar, to generate dimer-dimer interface "scores" that are translated into free-energy estimates. Analysis of 23 hemoglobin mutants, in both deoxy and oxy states, indicates that the effects of mutant residues on structurally bound waters (and visa versa) are important for generating accurate free-energy estimates. For several mutants, the addition/elimination of structural waters is key to understanding the thermodynamic consequences of residue mutation. Good agreement is found between calculated and experimental data for deoxy hemoglobin mutants (r = 0.79, slope = 0.78, standard error = 1.4 kcal mol(-1), n = 23). Less accurate estimates were initially obtained for oxy hemoglobin mutants (r = 0.48, slope = 0.47, standard error = 1.4 kcal mol(-1), n = 23). However, the elimination of three outliers from this data set results in a better correlation of r = 0.87 (slope = 0.72, standard error = 0.75, n = 20). These three mutations may significantly perturb the hemoglobin quaternary structure beyond the scope of our structural optimization procedure. The method described is also useful in the examination of residue ionization states in protein structures. Specifically, we find an acidic residue within the native deoxy hemoglobin dimer-dimer interface that may be protonated at physiological pH. The final analysis is a model design of novel hemoglobin mutants that modify cooperative free energy (deltaGc)--the energy barrier between the allosteric transition from deoxy to oxy hemoglobin.

Mutagenic dissection of hemoglobin cooperativity: effects of amino acid alteration on subunit assembly of oxy and deoxy tetramers

Free energies of oxygen-linked subunit assembly and cooperative interaction have been determined for 34 molecular species of human hemoglobin, which differ by amino acid alterations as a result of mutation or chemical modification at specific sites. These studies required the development of extensions to our earlier methodology. In combination with previous results they comprise a data base of 60 hemoglobin species, characterized under the same conditions. The data base was analyzed in terms of the five following issues. (1) Range and sensitivity to site modifications. Deoxy tetramers showed greater average energetic response to structural modifications than the oxy species, but the ranges are similar for the two ligation forms. (2) Structural localization of cooperative free energy. Difference free energies of dimer-tetramer assembly (oxy minus deoxy) yielded delta Gc for each hemoglobin, i.e., the free energy used for modulation of oxygen affinity over all four binding steps. A structure-energy map constructed from these results shows that the alpha 1 beta 2 interface is a unique structural location of the noncovalent bonding interactions that are energetically coupled to cooperativity. (3) Relationship of cooperativity to intrinsic binding. Oxygen binding energetics for dissociated dimers of mutants strongly indicates that cooperativity and intrinsic binding are completely decoupled by tetramer to dimer dissociation. (4) Additivity, site-site coupling and adventitious perturbations. All these are exhibited by individual-site modifications of this study. Large nonadditivity may be correlated with global (quaternary) structure change. (5) Residue position vs. chemical nature. Functional response is solely dictated by structural location for a subset of the sites, but varies with side-chain type at other sites. The current data base provides a unique framework for further analyses and modeling of fundamental issues in the structural chemistry of proteins and allosteric mechanisms.

Hb Rancho Mirage [beta 143(H21)His----Asp]; a variant in the 2,3-DPG binding site showing normal oxygen affinity at physiological pH

Hb Rancho Mirage was detected in a 17-year-old male in association with a mild anemia. Hemoglobin electrophoresis revealed the variant had a mobility between Hbs A and J on cellulose acetate (pH 8.6) and a mobility like Hb F on citrate agar (pH 6.4). A substitution of His----Asp was found at position 143 in the beta chain, a residue that contributes to the anionic 2,3-DPG binding site in Hb. This variant exhibited normal oxygen affinity at physiologic pH and reduced affinity at alkaline pH. This suggested a subtle shift in the allosteric equilibrium due most likely to the introduction of a negative charge that stabilized the 2,3-DPG pocket. Both homotrophic (heme-heme) and heterotropic (2,3-DPG and protons) effects were reduced; this might be a consequence of an alteration in the carboxyl terminal region of the beta-subunits. Although a His----Asp substitution would be considered to cause reasonable disruption of the 2,3-DPG and C-terminal conformation of the beta- subunits, the properties of Hb Rancho Mirage suggest that, in fact, there appear to be no major perturbation of the critical C-terminal residues.