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

De novo characterization of the spleen transcriptome of the large yellow croaker (Pseudosciaena crocea) and analysis of the immune relevant genes and pathways involved in the antiviral response

The large yellow croaker (Pseudosciaena crocea) is an economically important marine fish in China. To understand the molecular basis for antiviral defense in this species, we used Illumia paired-end sequencing to characterize the spleen transcriptome of polyriboinosinic:polyribocytidylic acid [poly(I:C)]-induced large yellow croakers. The library produced 56,355,728 reads and assembled into 108,237 contigs. As a result, 15,192 unigenes were found from this transcriptome. Gene ontology analysis showed that 4,759 genes were involved in three major functional categories: biological process, cellular component, and molecular function. We further ascertained that numerous consensus sequences were homologous to known immune-relevant genes. Kyoto Encyclopedia of Genes and Genomes orthology mapping annotated 5,389 unigenes and identified numerous immune-relevant pathways. These immune-relevant genes and pathways revealed major antiviral immunity effectors, including but not limited to: pattern recognition receptors, adaptors and signal transducers, the interferons and interferon-stimulated genes, inflammatory cytokines and receptors, complement components, and B-cell and T-cell antigen activation molecules. Moreover, the partial genes of Toll-like receptor signaling pathway, RIG-I-like receptors signaling pathway, Janus kinase-Signal Transducer and Activator of Transcription (JAK-STAT) signaling pathway, and T-cell receptor (TCR) signaling pathway were found to be changed after poly(I:C) induction by real-time polymerase chain reaction (PCR) analysis, suggesting that these signaling pathways may be regulated by poly(I:C), a viral mimic. Overall, the antivirus-related genes and signaling pathways that were identified in response to poly(I:C) challenge provide valuable leads for further investigation of the antiviral defense mechanism in the large yellow croaker.

Kinetic trapping of a key hemoglobin intermediate

The complete binding cascade of human hemoglobin consists of a series of partially ligated intermediates. The individual intermediate binding constants cannot be distinguished in O(2) binding curves, however, each constant can be determined from the O(2)-induced change in assembly constant for the α(2)β(2) tetramer from its constituent αβ dimers. The characterization of these O(2) binding constants has shown the Hb cascade to be asymmetric in nature, with binding dependent upon the specific distribution of O(2) among the four hemesites. A stopped-flow approach to measuring the dissociation constant of a key doubly ligated intermediate, that in which one dimer is oxygenated and the other is not, is described. The intermediate is transiently formed in the absence of O(2) and then allowed to dissociate in the presence of O(2). The free dimers thus released are trapped by the plasma protein haptoglobin, the rate limiting step being that of tetramer dissociation. The kinetic constant observed for the dissociation of this intermediate confirms the value for its equilibrium O(2) binding constant, previously determined under equilibrium conditions by subzero isoelectric focusing.

Intrinsic regulation of hemoglobin expression by variable subunit interface strengths

The expression of the six types of human Hb subunits over time is currently considered to be regulated mainly by transcription factors that bind to upstream control regions of the gene (the 'extrinsic' component of regulation). Here, we describe how subunit pairing and further assembly to tetramers in the liganded state is influenced by the affinity of subunits for one another (the 'intrinsic' component of regulation). The adult Hb dimers have the strongest subunit interfaces and the embryonic Hbs the weakest, with fetal Hbs being of intermediate strength, corresponding to the temporal order of their expression. These variable subunit binding strengths and the attenuating effects of acetylation contribute to the differences with which these Hb types form functional O(2) -binding tetramers consistent with gene switching.

Structural basis of cooperative ligand binding by the glycine riboswitch

The glycine riboswitch regulates gene expression through the cooperative recognition of its amino acid ligand by a tandem pair of aptamers. A 3.6 Å crystal structure of the tandem riboswitch from the glycine permease operon of Fusobacterium nucleatum reveals the glycine binding sites and an extensive network of interactions, largely mediated by asymmetric A-minor contacts, that serve to communicate ligand binding status between the aptamers. These interactions provide a structural basis for how the glycine riboswitch cooperatively regulates gene expression.

Identification of a tertiary interaction important for cooperative ligand binding by the glycine riboswitch

The glycine riboswitch has a tandem dual aptamer configuration, where each aptamer is a separate ligand-binding domain, but the aptamers function together to bind glycine cooperatively. We sought to understand the molecular basis of glycine riboswitch cooperativity by comparing sites of tertiary contacts in a series of cooperative and noncooperative glycine riboswitch mutants using hydroxyl radical footprinting, in-line probing, and native gel-shift studies. The results illustrate the importance of a direct or indirect interaction between the P3b hairpin of aptamer 2 and the P1 helix of aptamer 1 in cooperative glycine binding. Furthermore, our data support a model in which glycine binding is sequential; where the binding of glycine to the second aptamer allows tertiary interactions to be made that facilitate binding of a second glycine molecule to the first aptamer. These results provide insight into cooperative ligand binding in RNA macromolecules.

Kinetic and spectroscopic studies of hemin acquisition in the hemophore HasAp from Pseudomonas aeruginosa

The extreme limitation of free iron has driven various pathogens to acquire iron from the host in the form of heme. Specifically, several Gram-negative pathogens secrete a heme binding protein known as HasA to scavenge heme from the extracellular environment and to transfer it to the receptor protein HasR for import into the bacterial cell. Structures of heme-bound and apo-HasA homologues show that the heme iron(III) ligands, His32 and Tyr75, reside on loops extending from the core of the protein and that a significant conformational change must occur at the His32 loop upon heme binding. Here, we investigate the kinetics of heme acquisition by HasA from Pseudomonas aeruginosa (HasAp). The rate of heme acquisition from human met-hemoglobin (met-Hb) closely matches that of heme dissociation which suggests a passive mode of heme uptake from this source. The binding of free hemin is characterized by an initial rapid phase forming an intermediate before further conversion to the final complex. Analysis of this same reaction using an H32A variant lacking the His heme ligand shows only the rapid phase to form a heme-protein complex spectroscopically equivalent to that of the wild-type intermediate. Further characterization of these reactions using electron paramagnetic resonance and resonance Raman spectroscopy of rapid freeze quench samples provides support for a model in which heme is initially bound by the Tyr75 to form a high-spin heme-protein complex before slower coordination of the His32 ligand upon closing of the His loop over the heme. The slow rate of this loop closure implies that the induced-fit mechanism of heme uptake in HasAp is not based on a rapid sampling of the H32 loop between open and closed configurations but, rather, that the H32 loop motions are triggered by the formation of the high-spin heme-HasAp intermediate complex.

History and current status of newborn screening for hemoglobinopathies

The impact of hemoglobinopathies on healthcare in the United States, particularly sickle cell disease (SCD), has been significant. Enactment of the Sickle Cell Anemia Control Act in 1972 significantly increased the federal interest in the SCDs and other hemoglobinopathies. Only since May 1, 2006, have all states required and provided universal newborn screening for SCD despite a national recommendation to this effect in 1987. In this article, we review the history of screening for SCD and other hemoglobinopathies, along with federal and state activities that have contributed to improved health outcomes for patients with SCD, as well as current newborn screening practices. We also chronicle the federal activities that have helped to shape and to refine laboratory screening and diagnostic proficiency. Finally, we review molecular testing strategies that have evolved and outline their possible future impacts on disease detection and outcome improvement.

Energetic differences at the subunit interfaces of normal human hemoglobins correlate with their developmental profile

A previously unrecognized function of normal human hemoglobins occurring during protein assembly is described, i.e. self-regulation of subunit pairings and their durations arising from the variable strengths of their subunit interactions. Although many mutant human hemoglobins are known to have altered subunit interface strengths, those of the normal embryonic, fetal, and adult human hemoglobins have not been considered to differ significantly. However, in a comprehensive study of both types of subunit interfaces of seven of the eight normal oxy human hemoglobins, we found that the strengths, i.e., the free energies of the tetramer-dimer interfaces, contrary to previous reports, differ by 3 orders of magnitude and display an undulating profile similar to the transitions ("switches") of various globin subunit types over time. The dimer interface strengths are also variable and correlate linearly with their developmental profile. Embryonic hemoglobins are the weakest; fetal hemoglobin is of intermediate strength, and adult hemoglobins are the strongest. The pattern also correlates generally with their different O(2) affinities and responses to allosteric regulatory molecules. Acetylation of fetal hemoglobin weakens its unusually strong subunit interactions and occurs progressively as its level of expression diminishes and adult hemoglobin A formation begins; a causal relationship is suggested. The relative contributions of globin gene order and competition among subunits due to differences in their interface strengths were found to be complementary and establish a connection among genetics, thermodynamics, and development.

Thrombin

Thrombin is a Na+-activated, allosteric serine protease that plays opposing functional roles in blood coagulation. Binding of Na+ is the major driving force behind the procoagulant, prothrombotic and signaling functions of the enzyme, but is dispensable for cleavage of the anticoagulant protein C. The anticoagulant function of thrombin is under the allosteric control of the cofactor thrombomodulin. Much has been learned on the mechanism of Na+ binding and recognition of natural substrates by thrombin. Recent structural advances have shed light on the remarkable molecular plasticity of this enzyme and the molecular underpinnings of thrombin allostery mediated by binding to exosite I and the Na+ site. This review summarizes our current understanding of the molecular basis of thrombin function and allosteric regulation. The basic information emerging from recent structural, mutagenesis and kinetic investigation of this important enzyme is that thrombin exists in three forms, E*, E and E:Na+, that interconvert under the influence of ligand binding to distinct domains. The transition between the Na+ -free slow from E and the Na+ -bound fast form E:Na+ involves the structure of the enzyme as a whole, and so does the interconversion between the two Na+ -free forms E* and E. E* is most likely an inactive form of thrombin, unable to interact with Na + and substrate. The complexity of thrombin function and regulation has gained this enzyme pre-eminence as the prototypic allosteric serine protease. Thrombin is now looked upon as a model system for the quantitative analysis of biologically important enzymes.

Chemical basis of glycine riboswitch cooperativity

The glycine binding riboswitch forms a unique tandem aptamer structure that binds glycine cooperatively. We employed nucleotide analog interference mapping (NAIM) and mutagenesis to explore the chemical basis of glycine riboswitch cooperativity. Based on the interference pattern, at least two sites appear to facilitate cooperative tertiary interactions, namely, the minor groove of the P1 helix from aptamer 1 and the major groove of the P3a helix from both aptamers. Mutation of these residues altered both the cooperativity and binding affinity of the riboswitch. The data support a model in which the P1 helix of the first aptamer participates in a tertiary interaction important for cooperativity, while nucleotides in the P1 helix of the second aptamer interface with the expression platform. These data have direct analogy to well-characterized mutations in hemoglobin, which provides a framework for considering cooperativity in this RNA-based system.

Human embryonic, fetal, and adult hemoglobins have different subunit interface strengths. Correlation with lifespan in the red cell

The different types of naturally occurring, normal human hemoglobins vary in their tetramer-dimer subunit interface strengths (stabilities) by three orders of magnitude in the liganded (CO or oxy) state. The presence of embryonic zeta-subunits leads to an average 20-fold weakening of tetramer-dimer interfaces compared to corresponding hemoglobins containing adult alpha-subunits. The dimer-monomer interfaces of these hemoglobins differ by at least 500-fold in their strengths; such interfaces are weak if they contain zeta-subunits and exchange with added beta-subunits in the form of beta(4) (HbH) significantly faster than do those with alpha-subunits. Subunit exchange occurs at the level of the dimer, although tetramer formation reciprocally influences the amount of dimer available for exchange. Competition between subunit types occurs so that pairs of weak embryonic hemoglobins can exchange subunits to form the stronger fetal and adult hemoglobins. The dimer strengths increase in the order Hb Portland-2 (zeta(2)beta(2)) < Hb Portland-1 (zeta(2)gamma(2)) approximately equal Hb Gower-1 (zeta(2)epsilon(2)) < Hb Gower-2 (alpha(2)epsilon(2)) < HbF(1) < HbF (alpha(2)gamma(2)) < HbA(2) (alpha(2)delta(2)), i.e., from embryonic to fetal to adult types, representing maturation from weaker to stronger monomer-monomer subunit contacts. This increasing order recapitulates the developmental order in which globins are expressed (embryonic --> fetal --> adult), suggesting that the intrinsic binding properties of the subunits themselves regarding the strengths of interfaces they form with competing subunits play an important role in the dynamics of protein assemblies and networks.

Thrombin allostery

Thrombin is a Na(+)-activated, allosteric serine protease that plays multiple functional roles in blood pathophysiology. Binding of Na(+) is the major driving force behind the procoagulant, prothrombotic and signaling functions of the enzyme. This review summarizes our current understanding of the molecular basis of thrombin allostery with special emphasis on the kinetic aspects of Na(+) activation. The molecular mechanism of thrombin allostery is a remarkable example of long-range communication that offers a paradigm for many other biological systems.

Single residue modification of only one dimer within the hemoglobin tetramer reveals autonomous dimer function

The mechanism of cooperativity in the human hemoglobin tetramer (a dimer of alpha beta dimers) has historically been modeled as a simple two-state system in which a low-affinity structural form (T) switches, on ligation, to a high-affinity form (R), yielding a net loss of hydrogen bonds and salt bridges in the dimer-dimer interface. Modifications that weaken these cross-dimer contacts destabilize the quaternary T tetramer, leading to decreased cooperativity and enhanced ligand affinity, as demonstrated in many studies on symmetric double modifications, i.e., a residue site modified in both alpha- or both beta-subunits. In this work, hybrid tetramers have been prepared with only one modified residue, yielding molecules composed of a wild-type dimer and a modified dimer. It is observed that the cooperative free energy of ligation to the modified dimer is perturbed to the same extent whether in the hybrid tetramer or in the doubly modified tetramer. The cooperative free energy of ligation to the wild-type dimer is unperturbed, even in the hybrid tetramer, and despite the overall destabilization of the T tetramer by the modification. This asymmetric response by the two dimers within the same tetramer shows that loss of dimer-dimer contacts is not communicated across the dimer-dimer interface, but is transmitted through the dimer that bears the modified residue. These observations are interpreted in terms of a previously proposed dimer-based model of cooperativity with an additional quaternary (T/R) component.

Electron transfer kinetics between hemoglobin subunits

The kinetics for electron transfer have been measured for samples of hemoglobin valency hybrids with initially one type of subunit, alpha or beta, in the oxidized state. Incubation of these samples under anaerobic conditions tends to randomize the type of subunit that is oxidized. With a time coefficient of a few hours at pH 7, 25 degrees C, the Hb solution (0.1 mm heme) approaches a form with about 60% of beta chains reduced, indicating a faster transfer rate in the direction alpha to beta. There was no observable electron transfer for samples saturated with oxygen. The electron transfer occurs predominantly between deoxy and aquo-met subunits, both high spin species. Furthermore, electron transfer does not depend on the quaternary state of hemoglobin. Incubation of oxidized cross-linked tetramer Hb A with deoxy Hb S also displayed electron transfer, implying a mechanism via inter-tetramer collisions. A dependence on the overall Hb concentration confirms this mechanism, although a small contribution of transfer between subunits of the same tetramer cannot be ruled out. These results suggest that in vivo collisions between the Hb tetramers will be involved in the relative distribution of the methemoglobin between subunits in association with the reductase system present in the erythrocyte.

Very empirical treatment of solvation and entropy: a force field derived from log Po/w

A non-covalent interaction force field model derived from the partition coefficient of 1-octanol/water solubility is described. This model, HINT for Hydropathic INTeractions, is shown to include, in very empirical and approximate terms, all components of biomolecular associations, including hydrogen bonding, Coulombic interactions, hydrophobic interactions, entropy and solvation/desolvation. Particular emphasis is placed on: (1) demonstrating the relationship between the total empirical HINT score and free energy of association, deltaGinteraction; (2) showing that the HINT hydrophobic-polar interaction sub-score represents the energy cost of desolvation upon binding for interacting biomolecules; and (3) a new methodology for treating constrained water molecules as discrete independent small ligands. An example calculation is reported for dihydrofolate reductase (DHFR) bound with methotrexate (MTX). In that case the observed very tight binding, deltaGinteraction < or = -13.6 kcal mol(-1), is largely due to ten hydrogen bonds between the ligand and enzyme with estimated strength ranging between -0.4 and -2.3 kcal mol(-1). Four water molecules bridging between DHFR and MTX contribute an additional -1.7 kcal mol(-1) stability to the complex. The HINT estimate of the cost of desolvation is +13.9 kcal mol(-1).

Allosteric free energy changes at the alpha 1 beta 2 interface of human hemoglobin probed by proton exchange of Trp beta 37

The energetic changes that occur on ligand binding in human hemoglobin have been investigated by measurements of the exchange rates of the indole proton of Trpbeta37(C3). The Trpbeta37 residues are located in helices C of the beta-subunits and are involved in contacts with the segments FG of the alpha-subunits at the interdimeric alpha1beta2 and alpha2beta1 interfaces of the hemoglobin tetramer. In the quaternary structure change that accompanies ligand binding to hemoglobin, these contacts undergo minimal changes in relative orientation and in packing, thereby acting as hinges, or flexible joints. The exchange rates of the indole proton of Trpbeta37(C3) were measured by nuclear magnetic resonance spectroscopy, in both deoxygenated and ligated hemoglobin. The results indicate that, at 15 degrees C, the exchange rate is increased from 9.0. 10(-6) to 3.3. 10(-4) s(-1) upon ligand binding to hemoglobin. This change suggests that the structural units at the hinge regions of the alpha1beta2/alpha2beta1 interfaces containing Trpbeta37(C3) are specifically stabilized in unligated hemoglobin, and experience a change in structural free energy of approximately 4 kcal/(mol tetramer) upon ligand binding. Therefore, the hinge regions of the alpha1beta2/alpha2beta1 interfaces could play a role in the transmission of free energy through the hemoglobin molecule during its allosteric transition.

The Bohr effect of hemoglobin intermediates and the role of salt bridges in the tertiary/quaternary transitions

Understanding mechanisms in cooperative proteins requires the analysis of the intermediate ligation states. The release of hydrogen ions at the intermediate states of native and chemically modified hemoglobin, known as the Bohr effect, is an indicator of the protein tertiary/quaternary transitions, useful for testing models of cooperativity. The Bohr effects due to ligation of one subunit of a dimer and two subunits across the dimer interface are not additive. The reductions of the Bohr effect due to the chemical modification of a Bohr group of one and two alpha or beta subunits are additive. The Bohr effects of monoliganded chemically modified hemoglobins indicate the additivity of the effects of ligation and chemical modification with the possible exception of ligation and chemical modification of the alpha subunits. These observations suggest that ligation of a subunit brings about a tertiary structure change of hemoglobin in the T quaternary structure, which breaks some salt bridges, releases hydrogen ions, and is signaled across the dimer interface in such a way that ligation of a second subunit in the adjacent dimer promotes the switch from the T to the R quaternary structure. The rupture of the salt bridges per se does not drive the transition.

Confirmation of a unique intra-dimer cooperativity in the human hemoglobin alpha(1)beta(1)half-oxygenated intermediate supports the symmetry rule model of allosteric regulation

The contribution of the alpha(1)beta(1)half-oxygenated tetramer [alphabeta:alphaO(2)betaO(2)] (species 21) to human hemoglobin cooperativity was evaluated using cryogenic isoelectric focusing. The cooperative free energy of binding, reflecting O(2)-driven protein structure changes, was measured as (21)DeltaG(c) = 5.1 +/- 0. 3 kcal for the Zn/FeO(2) analog. For the Fe/FeCN analog, (21)DeltaG(c) was estimated as 4.0 kcal after correction for a CN ligand rearrangement artifact, demonstrating that ligand rearrangement does not invalidate previous conclusions regarding this species. In the context of the entire Hb cooperativity cascade, which includes eight intermediate species, the 21 tetramer is highly abundant relative to the other doubly-ligated species, providing strong support for the previously determined consensus partition function of O(2) binding and for the Symmetry Rule model of hemoglobin cooperativity (Ackers et al., Science 1992;255:54-63). Cooperativity of normal human hemoglobin is shown to depend on site-configuration, and not solely the number of O(2) bound, nor the occupancy of alpha vs. beta subunits. Verification of a unique contribution from the alpha(1)beta(1)doubly-oxygenated species to the equilibrium O(2) binding curve strongly reinforces the Symmetry Rule interpretation that the alpha(1)beta(1)dimer acts both as a structural and functional element in cooperative O(2) binding.

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.

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.

Hydrophobicity: is LogP(o/w) more than the sum of its parts?

The empirically calculated parameter LogP(o/w), the log(10) of the coefficient for solvent partitioning between 1-octanol and water, has been used to provide the key data for a unique non-covalent interaction force field called HINT (Hydropathic INTeractions). This experimentally-derived force field encodes entropic as well as enthalpic information and also includes some representation of solvation and desolvation energetics in biomolecular associations. The theoretical basis for the HINT model is discussed. This review includes: 1) discussion of calculational representation of the hydrophobic effect, 2) the rationale for describing the experimental LogP(o/w) based descriptors used in the HINT force field and model as free energy-like, 3) the relationship between hydrophobic fragment constants and partial group electrostatic charge, and 4) the implications of structurally-conserved water molecules on free energy of molecular association. Several recent applications of HINT in structure-based and ligand-based drug discovery are reviewed. Finally, future directions in the HINT model development are proposed.

Evolutionarily conserved pathways of energetic connectivity in protein families

For mapping energetic interactions in proteins, a technique was developed that uses evolutionary data for a protein family to measure statistical interactions between amino acid positions. For the PDZ domain family, this analysis predicted a set of energetically coupled positions for a binding site residue that includes unexpected long-range interactions. Mutational studies confirm these predictions, demonstrating that the statistical energy function is a good indicator of thermodynamic coupling in proteins. Sets of interacting residues form connected pathways through the protein fold that may be the basis for efficient energy conduction within proteins.

Additive effects of beta chain mutations in low oxygen affinity hemoglobin betaF41Y,K66T

In order to decrease significantly the oxygen affinity of human hemoglobin, we have associated the mutation betaF41Y with another point mutation also known to decrease the oxygen affinity of Hb. We have synthesized a recombinant Hb (rHb) with two mutations in the beta chains: rHb betaF41Y,K66T. In the absence of 2, 3-diphosphoglycerate, additive effects of the mutations are evident, since the doubly mutated Hb exhibits a larger decrease in oxygen affinity than for the individual single mutations. In the presence of 2,3-diphosphoglycerate, the second mutation did not significantly increase the P(50) value relative to the single mutations. However, the kinetics of CO binding still indicate combined effects on the allosteric equilibrium, as evidenced by more of the slow bimolecular phase characteristic of binding to the deoxy conformation. Dimer-tetramer equilibrium studies indicate an increase in stability of the mutants relative to rHb A; the double mutant rHb betaF41Y, K66T at pH 7.5 showed a K(4,2) value of 0.26 microM. Despite the lower oxygen affinity, the single mutant betaF41Y and double mutant betaF41Y,K66T show only a moderate increase of 20% in the autoxidation rate. These mutations are thus of interest in developing a Hb-based blood substitute.

The N-terminal sequence affects distant helix interactions in hemoglobin. Implications for mutant proteins from studies on recombinant hemoglobin felix

The N-terminal 18-amino acid sequence of the beta-chain of hemoglobin, as far as the end of the A helix, has been replaced by the corresponding sequence of the gamma-chain of fetal hemoglobin with the remaining sequence of the beta-chain retained (helices B through H). The gamma-beta-chain had the correct mass, and its entire sequence was established by mass spectrometric analysis of its tryptic peptides; the alpha-chain also had the correct mass. This recombinant hemoglobin (named Hb Felix) retains cooperativity and has an oxygen affinity like that of HbA both in the presence and absence of the allosteric regulators, 2,3-diphosphoglycerate or chloride but differs from HbF in its 2,3-diphosphoglycerate response. However, Hb Felix has some features that resemble fetal hemoglobin, i.e. its significantly decreased tetramer-dimer dissociation and its circular dichroism spectrum, which measure the strength of the tetramer-dimer interface in the oxy conformation and its rearrangement to the deoxy conformation, respectively. Even though Hb Felix contains the HbA amino acids at its tetramer-dimer interface, which is located at a distance from the substitution sites, its interface properties resemble those of HbF. Therefore, the N-terminal sequence and not just those amino acids directly involved at the subunit interface contacts with alpha-chains must have a strong influence on this region of the molecule. The results reinforce the concept of fluid long range relationships among various parts of the hemoglobin tetramer (Dumoulin, A., Manning, L. R., Jenkins, W. T., Winslow, R. M., and Manning, J. M. (1997) J. Biol. Chem. 272, 31326-31332) and demonstrate the importance of the N-terminal sequence, especially in some mutant hemoglobins, in influencing its overall structure by affecting the relationship between helices.

Energetic components of the allosteric machinery in hemoglobin measured by hydrogen exchange

A hydrogen exchange (HX) functional labeling method was used to study allosterically active segments in human hemoglobin (Hb) at the alpha-chain N terminus and the beta-chain C terminus. Allosterically important interactions that contact these segments were removed one or more at a time by mutation (Hbs Cowtown, Bunbury, Barcelona, Kariya), proteolysis (desArg141alpha, desHis146beta), chemical modification (N-ethylsuccinimidyl-Cys93beta), and the withdrawal of extrinsic effectors (phosphate groups, chloride). The effects of each modification on HX rate at the local and the remote position were measured in the deoxy Hb T-state and translated into change in structural free energy at each position.The removal of individual salt links destabilizes local structure by 0.4 to 0.75 kcal/mol (pH 7.4, 0 degreesC, 0.35 M ionic strength) and often produces cross-subunit effects while hemoglobin remains in the T-state. In doubly modified hemoglobins, different changes that break the same links produce identical destabilization, changes that are structurally independent show energetic additivity, and changes that intersect show energetic overlap. For the overall T-state to R-state transition and for some but not all modifications within the T-state, the summed loss in stabilization free energy measured at the two chain termini matches the total loss in allosteric free energy measured by global methods. These observations illustrate the importance of evaluating the detailed energetics and the modes of energy transfer that define the allosteric machinery.

Signal transmission between subunits in the hemoglobin T-state

To study allosteric mechanism in hemoglobin, a hydrogen-exchange method was used to measure ligand-dependent changes in structural free energy at defined allosterically sensitive positions. When the two alpha-subunits are CN-met liganded, effects can be measured locally, within the alpha-subunit, and also remotely, within the beta-subunit, even though the quaternary structure remains in the T conformation. When the two beta-subunits are liganded, effects occur at the same positions. The effects seen are the same, independently of whether ligands occupy the alpha-chain hemes or the beta-chain hemes. Control experiments rule out modes of energy transfer other than programmed cross-subunit interaction within the T-state. Cross-subunit transfer may depend on pulling the heme trigger (moving the heme iron into the heme plane) rather than on liganding alone.

Tetramer-dimer equilibrium of oxyhemoglobin mutants determined from auto-oxidation rates

One of the main difficulties with blood substitutes based on hemoglobin (Hb) solutions is the auto-oxidation of the hemes, a problem aggravated by the dimerization of Hb tetramers. We have employed a method to study the oxyHb tetramer-dimer equilibrium based on the rate of auto-oxidation as a function of protein concentration. The 16-fold difference in dimer and tetramer auto-oxidation rates (in 20 mM phosphate buffer at pH 7.0, 37 degrees C) was exploited to determine the fraction dimer. The results show a transition of the auto-oxidation rate from low to high protein concentrations, allowing the determination of the tetramer-dimer dissociation coefficient K4,2 = [Dimer] 2/[Tetramer]. A 14-fold increase in K4,2 was observed for addition of 10 mM of the allosteric effector inositol hexaphosphate (IHP). Recombinant hemoglobins (rHb) were genetically engineered to obtain Hb with a lower oxygen affinity than native Hb (Hb A). The rHb alpha2beta2 [(C7) F41Y/(G4) N102Y] shows a fivefold increase in K4,2 at pH 7.0, 37 degrees C. An atmosphere of pure oxygen is necessary in this case to insure fully oxygenated Hb. When this condition is satisfied, this method provides an efficient technique to characterize both the tetramer-dimer equilibrium and the auto-oxidation rates of various oxyHb. For low oxygen affinity Hb equilibrated under air, the presence of deoxy subunits accelerates the auto-oxidation. Although a full analysis is complicated, the auto-oxidation studies for air equilibrated samples are more relevant to the development of a blood substitute based on Hb solutions. The double mutants, rHb alpha2beta2 [(C7) F41Y/(G4) N102A] and rHb alpha2beta2 [(C7) F41Y/(E10) K66T], show a lower oxygen affinity and a higher rate of oxidation than Hb A. Simulations of the auto-oxidation rate versus Hb concentration indicate that very high protein concentrations are required to observe the tetramer auto-oxidation rate. Because the dimers oxidize much more rapidly, even a small fraction dimer will influence the observed oxidation rate.

Exchange of subunit interfaces between recombinant adult and fetal hemoglobins. Evidence for a functional inter-relationship among regions of the tetramer

The inter-relationship between the interior subunit interfaces and the exterior diphosphoglycerate (DPG) binding region of the hemoglobin tetramer and the effects of a specific N-terminal acetylation on tetramer assembly have been evaluated. Tetrameric fetal hemoglobin F in the liganded state was found to dissociate to dimers much less than previously appreciated, i.e. about 70 times less than adult hemoglobin A (Kd = 0.01 microM and 0.68 microM, for HbF and HbA, at pH 7.5, respectively) over the pH range 6.2-7.5, whereas HbF1, in which the N termini of the gamma-chains are acetylated, dissociates like HbA. To determine whether this feature of HbF could be transferred to hemoglobin A, the single amino acid difference in their alpha1beta2/alpha1gamma2 interfaces and the 4 amino acid differences in their alpha1beta1/alpha1gamma1 interfaces have been substituted in HbA to those in HbF. This pentasubstituted recombinant HbA/F had the correct molecular weight as determined by mass spectrometry, the expected mobility on isoelectric focusing, the calculated amino acid composition, and normal circular dichroism properties, oxygen binding, and cooperativity. Although HbA/F has the same amino acid side chains that bind DPG as HbA, its diminished response to 2,3-DPG resembled that of HbF. However, its tetramer-dimer dissociation constant (Kd = 0.14 microM) was between that of HbA and HbF despite the fact that it was composed entirely of HbF subunit interfaces. The results indicate that regions of the tetramer distant from the tetramer-dimer interface influence its dissociation and, reciprocally, that the interfaces affect regions involved in the binding of allosteric regulators, suggesting flexible long range inter-relationships in hemoglobin.

Hydropathic analysis of the non-covalent interactions between molecular subunits of structurally characterized hemoglobins

The software program, HINT (Hydropathic INTeractions), which characterizes non-polar-non-polar, polar-polar, and non-polar-polar interactions, has been used to examine subunit interface associations involved in the hemoglobin allosteric transition at a residue and atomic level. HINT differs from many other computational programs in that it is based not on a statistical method or a force-field but employs parameters experimentally determined from solvent transfer experiments. The main focus of this study is to compare HINT scores that are based upon experimentally and thermodynamically derived measurements with experimentally determined thermodynamic results. The HINT analysis yields a good first-order approximation of experimentally measured energies for these interactions as determined by free energies of dimer-tetramer assembly for mutant hemoglobins. The results provide a framework for understanding subunit stabilities based upon individual atom interactions and repulsions. HINT, in agreement with previous analyses, indicates that: (1) the alpha1beta1 and alpha2beta2 subunit contacts are stabilized via several polar and many hydrophobic interactions with few repulsive contact areas in both the T (deoxyhemoglobin) and R (oxyhemoglobin) structures; (2) the alpha1alpha2 subunit contacts are primarily stabilized by polar salt bridge linkages in both T and R states; and (3) the alpha1beta2 and alpha2beta1 contacts have both strong positive and negative interactions in both T and R states with few hydrophobic interactions. The HINT scoring methodology provides a quantitative characterization of the major role of the alpha1beta2 and alpha2beta1 interfaces in the T-->R quaternary transition. HINT also confirms the stronger hydrogen bond formation in mutant Hb Rothschild (Trp 37beta-->Arg) with Asp94alpha1 that gives rise to a low-affinity (deoxy) hemoglobin. HINT shows that the stabilization of the alpha1beta2 interface with mutant Hb Ypsilanti (Asp99alpha-->Tyr) produces a high-affinity (oxy) hemoglobin by reducing hydrophobic-polar contacts in the R state. HINT interaction maps also identified specific sites for mutagenesis at the alpha1beta2 interface that can be explored to shift the allosteric equilibrium in either direction. In addition, the HINT program provides useful diagnostic data for checking the quality of refined crystallographic structures.

Detection of changes in pairwise interactions during allosteric transitions: coupling between local and global conformational changes in GroEL

A protein engineering approach for detecting and measuring local conformational changes that accompany allosteric transitions in proteins is described. Using this approach, we can identify interactions that are made or broken during allosteric transitions. The method is applied to probe for changes in pairwise interactions in the chaperonin GroEL during its ATP-induced allosteric transitions. Two pairwise interactions are investigated: one between subunits (Asp-41 with Thr-522) and the other within subunits (Glu-409 with Arg-501). We find that the intraring intersubunit interaction between Asp-41 and Thr-522 changes little during the allosteric transitions of GroEL, indicating that the hydrogen bond between these residues is maintained. In contrast, the intrasubunit salt bridge between Glu-409 and Arg-501 becomes significantly weaker during the ATP-induced allosteric transitions of GroEL. Our results are consistent with the electron microscopy observations of an ATP-induced hinge movement of the apical domains relative to the equatorial domains.

Effects of NaCl on the linkages between O2 binding and subunit assembly in human hemoglobin: titration of the quaternary enhancement effect

Oxygen binding by human hemoglobin (Hb) and the coupled reactions of dimer-tetramer assembly were studied over a range of NaCl concentrations (from 0.08 M to 1.4 M) at pH 7.4 and 21.5 degrees C. A strategy of multi-dimensional analysis was employed [G.K. Ackers and H.R. Halvorson, Proc. Natl. Acad. Sci. U.S.A., 91, (1974) 4312] to optimize the resolution of the contributions to cooperativity and their heterotropic salt linkages at each stoichiometric degree of O2 binding. A wide range of Hb concentration was utilized at each [NaCl] in which O2-linked subunit assembly reactions contributed significantly to the positions and shapes of the binding isotherms. Kinetic determinations yielded forward and reverse rate constants for assembly of the unligated species. Amplitudes for the assembly rate data had concentration dependences in agreement with the independently determined dimer-tetramer assembly constants of oxyhemoglobin. Concentration-dependent binding isotherms were analyzed, in combination with the kinetically determined equilibrium constants, to yield salt-linked components of cooperativity at the four stages of oxygenation. The principal results of this study were as follows. (i) Assembly of fully oxygenated Hb tetramers is opposed by NaCl: the dimer-to-tetramer equilibrium constant becomes two orders of magnitude less favorable over the [NaCl] range 0.08 M to 1.4 M. By contrast, for deoxy-Hb the assembly equilibrium constant is reduced only two-fold. (ii) Oxygen binding to dimers is non-cooperative over the entire salt range, whereas dimer affinity is slightly favored by increasing the NaCl concentration. (iii) Overall affinity of tetramers for O2 is opposed by NaCl, becoming an order of magnitude less favorable over the range employed. Most of this decrease occurs at the fourth binding step, which shows a large, salt-mediated quaternary enhancement effect; i.e., the assembly of dimers into tetramers at 0.08 M NaCl is accompanied by an eight-fold increase in O2 affinity. (iv) The quaternary enhancement effect at the last O2-binding step is titrated progressively by salt until it reaches a negligible value near the highest [NaCl] of this study. The lowest [NaCl] condition (0.08 M) elicits the greatest tetramer cooperativity with the largest maximal Hill coefficient and the greatest suppression of intermediates. Possible origins and mechanistic implications of these phenomena are considered.

Hb Bruxelles, deletion of Phebeta42, shows a low oxygen affinity and low cooperativity of ligand binding

Functional studies of partially purified hemoglobin (Hb) Bruxelles, Phebeta42 (CD1) --> 0 indicate a major shift in the allosteric equilibrium toward the deoxy (T state) conformation. While Hb A shows a roughly symmetrical oxygenation curve with maximum cooperativity near half-saturation, Hb Bruxelles shows mainly properties of the low affinity (T state) form. The oxygen equilibrium curves for purified (>80%) Hb Bruxelles show little cooperativity and a P50 (without 2,3-diphosphoglycerate) about twice that of Hb A. The low cooperativity for Hb Bruxelles is partially compensated by an increase in oxygen affinity of the deoxy conformation and a lower 2,3-diphosphoglycerate effect. The beta chains of normal Hb have consecutive phenylalanine residues at positions 41 and 42. DNA sequencing studies of Hb Bruxelles showed a deletion of the codon TTT, which corresponds to residue Phe42. The CO rebinding kinetics after flash photolysis show mainly the slow phase, characteristic of CO binding to the deoxy conformation. In phosphate buffer at pH 7, the slow phase dominates even at low photolysis levels, where the main reaction is ligand binding to the triply liganded form. This indicates a switchover point, from the deoxy to oxy conformation, occurring beyond three ligands for Hb Bruxelles. There are few natural mutants that show a change in the oxygen affinity and cooperativity as large as that observed for Hb Bruxelles.

The oxygen-binding intermediates of human hemoglobin: evaluation of their contributions to cooperativity using zinc-containing hybrids

Hemoglobin tetramers [Zn/FeO(2)] containing oxygenated subunits (FeO(2)), in combination with unligated subunits containing zinc-substituted hemes (Zn), were analyzed to determine their contributions to the cooperativity of oxygen binding at the Fe sites. Energetic consequences of possible perturbation by zinc substitution were evaluated in all combinations of unligated Zn/Fe hybrid tetramers. A general thermodynamic strategy that corrects for the energetic effects of substituting a second metal for Fe showed the perturbations of Zn substitution to be negligible. This permitted cooperativity parameters of the native Fe/FeO(2) intermediates to be calculated from data on the corresponding Zn/FeO(2) molecules. These parameters, determined explicitly for all eight oxygen-binding intermediates (Fe/FeO(2)), were found to be identical to those predicted earlier from analyzing the O(2) binding data of normal hemoglobin according to the "molecular code" of hemoglobin allostery. The cooperativity parameters determined for this system showed the same distribution pattern found previously for five other oxygen analog systems (Fe/FeCN, FE/Mn(3+), Co/FECO, Co/FeCN, and Fe/FeCO).

Heterometallic hybrids of homometallic human hemoglobins

Hybridization experiments between normal Hb tetramers (Fe2+ Hb) and those with four metal-substituted hemes (i.e., replacement of Fe2+ by Co2+, Mg2+, Mn2+, Mn3+, Ni2+, or Zn2+) have revealed unexpected behavior. These homometallic Hbs have previously served as models that mimic the deoxy or oxy properties of normal Fe2+ Hb. In this study, hybrids were composed of one alpha 1 beta 1 dimer that is metal-substituted at both hemes, in association with a second dimer alpha 2 beta 2 that has normal Fe2+ hemes. Both metal-substituted subunits are unligated, whereas the two Fe2+ subunits either are both unligated or both ligated with O2, CO, or CN. It was found that four of the metal-substituted Hbs (Mg2+ Hb, Mn2+ Hb, Ni2+ Hb, and Zn2+ Hb) did not form detectable amounts of heterometallic hybrids with normal Fe2+ Hb even though (i) their homometallic parents formed tight tetrameric complexes with stabilities similar to that of Fe2+ Hb and (ii) hybrids with metal substitution at both alpha sites or both beta sites are known to form readily. This striking positional effect was independent of whether the normal Fe2+ hemes were ligated and of which ligand was used. These findings indicate that surprisingly large changes in tetramer behavior can arise from small and subtle perturbations at the heme sites. Possible origins of these effects are considered.

Stability of the heme-globin linkage in alphabeta dimers and isolated chains of human hemoglobin. A study of the heme transfer reaction from the immobilized proteins to albumin

The stability of the heme-globin linkage in alphabeta dimers and in the isolated chains of human hemoglobin has been probed by studying the transfer of heme from the proteins immobilized onto CNBr-activated Sepharose 4B to human albumin. The kinetic and equilibrium features of the reaction have been measured spectrophotometrically given the stability of the heme donors and the ease with which heme donor and acceptor can be separated. Isolated alpha and beta chains transfer heme to albumin at similar rates (1 6 x 10(-2) s-1 at pH 9.0 and 20 degrees C) in the ferrous CO-bound and in the ferric state. In alpha beta dimers the heme-globin linkage is strengthened considerably, albeit to a different extent in the ferrous CO-bound and ferric met-aquo derivatives. Only in the latter heme is lost at a measurable rate, 0.065 +/- 0.011 x 10(-2) s-1 for alpha heme and 2.8 +/- 0.6 x 10(-2) s-1 for beta heme at pH 9.0 and 20 degrees C, which is very close to the rate measured with soluble met-aquo-hemoglobin at micromolar concentrations. These results indicate that in human hemoglobin the heme-globin linkage in the alpha chains is stabilized by interactions between unlike chains at the alpha1 beta1 interface, whereas heme binding to the beta chains is stabilized by interactions at the alpha1beta2 interface. These long range factors have to be taken into account in addition to the local factors at the heme pocket when evaluating the effect of point mutation and chemical modification.

Subunit dissociations in natural and recombinant hemoglobins

A precise and rapid procedure employing gel filtration on Superose-12 to measure the tetramer-dimer dissociation constants of some natural and recombinant hemoglobins in the oxy conformation is described. Natural sickle hemoglobin was chosen to verify the validity of the results by comparing the values with those reported using an independent method not based on gel filtration. Recombinant sickle hemoglobin, as well as a sickle double mutant with a substitution at the Val-6(beta) receptor site, had approximately the same dissociation constant as natural sickle hemoglobin. Of the two recombinant hemoglobins with amino acid replacements in the alpha 1 beta 2 subunit interface, one was found to be extensively dissociated and the other completely dissociated. In addition, the absence of an effect of the allosteric regulators DPG and IHP on the dissociation constant was demonstrated. Thus, a tetramer dissociation constant can now be determined readily and used together with other criteria for characterization of hemoglobins and their interaction with small regulatory molecules.

Coupled reactions in hemoglobin. Heme-globin and dimer-dimer association

Five different human hemoglobins were used to test the postulate that dissociation of hemoglobin (Hb) tetramers into alpha beta dimers and dissociation of heme from globin are linked reactions. Spectrophotometric measurements of the initial rate of heme transfer from Hb to serum albumin were made over a 3000-fold range of Hb concentration and yielded the heme-globin dissociation rate constant for tetramers and that for dimers. The tetramer-dimer dissociation constant (K4,2) could then be calculated from the rate constant at intermediate concentrations. The values obtained for the five hemoglobins, spanning a 250-fold range in K4,2, were in good agreement with those found by direct methods. The relation between this new linkage reaction of hemoglobin and the classical ones, such as the reciprocal relation between the binding of oxygen and protons, is discussed briefly.

FT-IR difference spectroscopy of hemoglobins A and Kempsey: evidence that a key quaternary interaction induces protonation of Asp beta 99

Fourier transform infrared difference spectra are reported for the CO adduct of human hemoglobin versus deoxyHb, in H2O and D2O. In addition to the well-known CO stretching and S-H(D) stretching bands, the difference spectra reveal numerous bands in the 1200-1700 cm-1 region, a number of which are assigned. Several amide modes are identified via their frequencies and D2O sensitivities. Bands arising from histidine protonation have also been found via comparison of the difference spectra at different pH(D) values, with the aid of aqueous histidine spectra. Of particular interest is the observation of a negative band at 1697 cm-1, which is assigned to the C = O stretch of carboxylic acid. This carboxylic acid is tentatively identified as the side chain of Asp beta 99, because it is missing in the difference spectrum of Hb Kempsey, a mutant in which Asp beta 99 is replaced by Asn. Asp beta 99 forms a critical contact with Tyr alpha 42 across the alpha 1 beta 2 interface in deoxyHb, which is broken upon ligation. Protonation of Asp beta 99 in deoxyHb is consistent with UV resonance Raman evidence that Tyr alpha 42 is the acceptor rather than the donor of the quaternary H-bond.

A recombinant human hemoglobin with asparagine-102(beta) substituted by alanine has a limiting low oxygen affinity, reduced marginally by chloride

A recombinant (r) mutant hemoglobin (Hb) with Asn-102(beta) replaced by an Ala (N102A(beta)) has been prepared by PCR amplification of a mutagenic DNA fragment and expression of the recombinant protein in yeast. The side chain of Asn-102(beta) is part of an important region of the alpha 1 beta 2 interface that undergoes large structural changes in the transition between the deoxy and oxy conformations. Three natural mutant Hbs with neutral substitutions of Thr, Ser, or Tyr at this site have low oxygen affinities because a hydrogen bond between Asn-102(beta) and Asp-94(alpha) in normal HbA was considered to be absent in these mutants, thereby destabilizing the oxy conformation in favor of the deoxy conformation. This proposal has been tested by expression of an rHb containing alanine at position 102(beta); alanine was chosen because its methyl side chain cannot participate in hydrogen bond formation, yet it is small enough not to disrupt the subunit interface. The nature of the desired replacement was established by sequencing the entire mutated beta-globin gene as well as the tryptic peptide containing the substitution. Further characterization by SDS-PAGE, isoelectric focusing, HPLC analysis, mass spectrometry, amino acid analysis, and sequencing of the mutant tryptic peptide confirmed the purity of the rHb. Its oxygen binding curve (2.4 mM in heme) in the absence of chloride showed that it had a very low oxygen affinity with a P50 of 42 mm Hg.(ABSTRACT TRUNCATED AT 250 WORDS)