Allosteric interpretation of haemoglobin properties

Wandering about allostery

It was a Lucky Strike to be working with Eraldo Antonini on hemoglobin and myoglobin when Jeffries Wyman arrived in Rome in 1961. I found myself connected with a number of creative scientists when the concept of allosteric control was conceived and gifted to the life science community. In retrospect, this was a demonstration of the skill and imagination of a few intelligent scientists that I happened to be close to. Those talents demonstrated the power of creativity as pictured by the motto "Mens agitat molem"; a celebration of humanism and intellect that paved the way to novel discoveries in the field of structure function relationships in proteins. I have presented hereby some of the events and the people as emerged from my memory over three decades of exciting scientific life.

Proteome microarray identifies autoantibody biomarkers for diagnosis of hepatitis B-related hepatocellular carcinoma

BACKGROUND AND AIMS

Hepatocellular carcinoma (HCC) has the highest mortality rate among malignant tumors worldwide. This study aimed to analyze the biological characteristics of serum proteins in hepatitis B (HBV)-related liver diseases, identify diagnostic biomarkers for HBV-infected HCC, and provide a scientific basis for its prevention and treatment.

MATERIALS AND METHODS

We used HuProt arrays to identify candidate biomarkers for HBV-related liver diseases and verified the differential biomarkers by using an HCC-focused array. The biological characteristics of serum proteins were analyzed via bioinformatics. Serum biomarkers levels were validated by ELISA.

RESULTS

We identified 547 differentially expressed proteins from HBV-infected HCC in a screening cohort. After analyzing the biological characteristics of serum proteins, we identified 10 potential differential autoantibodies against tumor-associated antigens (TAAbs) and a candidate biomarker panel (APEX2, RCSD1, and TP53) for the diagnosis of HBV-associated HCC with 61.9% sensitivity and 81.7% specificity in an HCC-focused array validation cohort. Finally, the protein levels and diagnostic capability of the biomarker panel were confirmed in a large-sample validation cohort, and this panel was found to be superior to alpha-fetoprotein, the standard hallmark for the diagnosis of HCC.

CONCLUSION

The APEX2, RCSD1, and TP53 biomarker panels could be used for the diagnosis of HBV-associated HCC, providing a scientific basis for clinical practice.

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.

A retrospective on statistical mechanical models for hemoglobin allostery

Understanding allosteric interactions in proteins has become one of the major research areas in protein science. The original aim of the famous theoretical model of Monod, Wyman, and Changeux (MWC) was to explain the regulation of enzymatic activity in biochemical pathways. However, its first successful quantitative application was to explain cooperative oxygen binding by hemoglobin, often called the "hydrogen molecule of biology." The combination of its original application and the enormous amount of research on hemoglobin has made it the paradigm for studies of allostery, especially for multi-subunit proteins, and for the development of statistical mechanical models to describe how structure determines function. This article is a historical account of the development of statistical mechanical models for hemoglobin to explain both the cooperative binding of oxygen (called homotropic effects by MWC) and how oxygen binding is affected by ligands that bind distant from the heme oxygen binding site (called heterotropic allosteric effects by MWC). This account makes clear the many remaining challenges for describing the relationship of structure to function for hemoglobin in terms of a satisfactory statistical mechanical model.

Control of Oxygen Affinity in Mammalian Hemoglobins: Implications for a System Biology Description of the Respiratory Properties of the Red Blood Cell

Hemoglobin and myoglobin have been considered for a long time the paradigmatic model systems for protein function, to the point of being defined the "hydrogen atom[s] of biology". Given this privileged position and the huge amount of quantitative information available on these proteins, the red blood cell might appear as the model system and"hydrogen atom" of system biology. Indeed, since the red cell's main function is O2 transport by hemoglobin, the gap between the protein and the cell may appear quite small. Yet, a surprisingly large amount of detailed biochemical information is required for the modelization of the respiratory properties of the erythrocyte. This problem is compounded if modelization aims at uncovering or explaining evolutionarily selected functional properties of hemoglobin. The foremost difficulty lies in the fact that hemoglobins having different intrinsic properties and relatively ancient evolutionary divergence may behave similarly in the complex milieu of blood, whereas very similar hemoglobins sharing a substantial sequence similarity may present important functional differences because of the mutation of a few key residues. Thus, the functional properties of hemoglobin and blood may reflect more closely the recent environmental challenges than the remote evolutionary history of the animal. We summarize in this review the case of hemoglobins from mammals, in an attempt to provide a reasoned summary of their complexity that, we hope, may be of help to scientists interested in the quantitative exploration of the evolutionary physiology of respiration. Indeed the basis of a meaningful modelization of the red cell requires a large amount of information collected in painstaking and often forgotten studies of the biochemical properties of hemoglobin carried out over more than a century.

Direct observation of ligand migration within human hemoglobin at work

Human hemoglobin is the textbook example of the stereochemistry of an allosteric protein and of the exquisite control that a protein can exert over ligand binding. However, the fundamental basis by which the protein facilitates the ligand movement remains unknown. In this study, we used cryogenic X-ray crystallography and a high-repetition pulsed laser irradiation technique to elucidate the atomic details of ligand migration processes in hemoglobin after photolysis of the bound CO. Our data clarify the distinct CO migration pathways in the individual subunits of hemoglobin and unravel the functional roles of the internal cavities and neighboring amino acid residues in ligand exit and entry. Our results also demonstrate the high gas permeability and porosity of hemoglobin, facilitating O₂ delivery.

Hemoglobin is one of the best-characterized proteins with respect to structure and function, but the internal ligand diffusion pathways remain obscure and controversial. Here we captured the CO migration processes in the tense (T), relaxed (R), and second relaxed (R2) quaternary structures of human hemoglobin by crystallography using a high-repetition pulsed laser technique at cryogenic temperatures. We found that in each quaternary structure, the photodissociated CO molecules migrate along distinct pathways in the α and β subunits by hopping between the internal cavities with correlated side chain motions of large nonpolar residues, such as α14Trp(A12), α105Leu(G12), β15Trp(A12), and β71Phe(E15). We also observe electron density evidence for the distal histidine [α58/β63His(E7)] swing-out motion regardless of the quaternary structure, although less evident in α subunits than in β subunits, suggesting that some CO molecules have escaped directly through the E7 gate. Remarkably, in T-state Fe(II)-Ni(II) hybrid hemoglobins in which either the α or β subunits contain Ni(II) heme that cannot bind CO, the photodissociated CO molecules not only dock at the cavities in the original Fe(II) subunit, but also escape from the protein matrix and enter the cavities in the adjacent Ni(II) subunit even at 95 K, demonstrating the high gas permeability and porosity of the hemoglobin molecule. Our results provide a comprehensive picture of ligand movements in hemoglobin and highlight the relevance of cavities, nonpolar residues, and distal histidines in facilitating the ligand migration.

Water, Ions, and Hemoglobin: Effects on Allostery and Polymerization

Proteins that function in aqueous solution can be perturbed by the solvent. Here we present experimental studies on two such interactions in the hemoglobin molecule. (1) Hemoglobin's oxygen binding is altered by introduction of crowding species or osmoticants, such as sucrose, through the linked binding of ions such as Cl or CO₂, but not otherwise. This rules out a significant role of buried surface in the allosteric energetics. (2) Sickle hemoglobin (HbS) polymerizes more readily in high concentrations of phosphate buffer. Such polymerization is analyzed quantitatively here for the first time in terms of the double nucleation mechanism. The changes in solubility are found to account for the increase in monomer addition rates and nucleation rates without requiring additional parameter adjustments. In the analysis, we also show how the analytical formulation of HbS nucleation may be adapted to include water that occupies the interstices between the assembled molecules. While such a "correction" has been applied to the equilibrium process, it has not previously been applied to the nucleation process.

Experiments on Hemoglobin in Single Crystals and Silica Gels Distinguish among Allosteric Models

Trapping quaternary structures of hemoglobin in single crystals or by encapsulation in silica gels has provided a demanding set of data to test statistical mechanical models of allostery. In this work, we compare the results of those experiments with predictions of the four major allosteric models for hemoglobin: the quaternary two-state model of Monod, Wyman, and Changeux; the tertiary two-state model of Henry et al., which is the simplest extension of the Monod-Wyman-Changeux model to include pre-equilibria of tertiary as well as quaternary conformations; the structure-based model of Szabo and Karplus; and the modification of the latter model by Lee and Karplus. We show that only the tertiary two-state model can provide a near quantitative explanation of the single-crystal and gel experimental results.

Experimental basis for a new allosteric model for multisubunit proteins

Monod, Wyman, and Changeux (MWC) explained allostery in multisubunit proteins with a widely applied theoretical model in which binding of small molecules, so-called allosteric effectors, affects reactivity by altering the equilibrium between more reactive (R) and less reactive (T) quaternary structures. In their model, each quaternary structure has a single reactivity. Here, we use silica gels to trap protein conformations and a new kind of laser photolysis experiment to show that hemoglobin, the paradigm of allostery, exhibits two ligand binding phases with the same fast and slow rates in both R and T quaternary structures. Allosteric effectors change the fraction of each phase but not the rates. These surprising results are readily explained by the simplest possible extension of the MWC model to include a preequilibrium between two tertiary conformations that have the same functional properties within each quaternary structure. They also have important implications for the long-standing question of a structural explanation for the difference in hemoglobin oxygen affinity of the two quaternary structures.

Dynamic prestress in a globular protein

A protein at equilibrium is commonly thought of as a fully relaxed structure, with the intra-molecular interactions showing fluctuations around their energy minimum. In contrast, here we find direct evidence for a protein as a molecular tensegrity structure, comprising a balance of tensed and compressed interactions, a concept that has been put forward for macroscopic structures. We quantified the distribution of inter-residue prestress in ubiquitin and immunoglobulin from all-atom molecular dynamics simulations. The network of highly fluctuating yet significant inter-residue forces in proteins is a consequence of the intrinsic frustration of a protein when sampling its rugged energy landscape. In beta sheets, this balance of forces is found to compress the intra-strand hydrogen bonds. We estimate that the observed magnitude of this pre-compression is enough to induce significant changes in the hydrogen bond lifetimes; thus, prestress, which can be as high as a few 100 pN, can be considered a key factor in determining the unfolding kinetics and pathway of proteins under force. Strong pre-tension in certain salt bridges on the other hand is connected to the thermodynamic stability of ubiquitin. Effective force profiles between some side-chains reveal the signature of multiple, distinct conformational states, and such static disorder could be one factor explaining the growing body of experiments revealing non-exponential unfolding kinetics of proteins. The design of prestress distributions in engineering proteins promises to be a new tool for tailoring the mechanical properties of made-to-order nanomaterials.

A tensegrity structure is one composed of members that are permanently under either tension or compression, and the balance of these tensile and compressive forces provides the structure with its mechanical stability. Macroscale tensegrity structures, which include Buckminster Fuller's geodesic domes, achieve exceptional structural integrity with a minimal use of resources. The question we address in this work is whether nature makes use of molecular-scale tensegrity in the design of proteins. Using Molecular Dynamics simulations of the protein ubiquitin, we measure the network of pairwise forces connecting the amino acid residues and show that this network does indeed have the character of a tensegrity structure. Furthermore, we find that the arrangement of tensile and compressive forces is such that hydrogen bonds in the protein's beta sheet, which are crucial for bearing mechanical loads, are compressed. This pre-compression is enough to significantly lengthen the lifetime of a bond under a given force, and thus should be an important factor in determining the protein's mechanical strength. The rational design of molecular prestress networks promises to be a new avenue for the engineering of proteins with made-to-order mechanical properties, for applications in medicine, materials and nanotechnology.

Hemoglobin allostery: variations on the theme

The two-state allosteric model of Monod, Wyman and Changeux (1965) offers a simple and elegant, yet very powerful and comprehensive, description of the functional behavior of hemoglobin. Although the extensive body of structural and functional information available is by-and-large consistent with this conceptual framework, some discrepancies between theory and experiment have been extensively discussed and considered to demand modifications of the original hypothesis. More recently the role of tertiary structural changes has been re-analyzed leading to extended kinetic models or indicating that powerful heterotropic effectors may be of paramount importance in controlling the function of human hemoglobin. The aim of this review is to analyze, and possibly reconcile, some discrepancies. We always felt that by looking at hemoglobins other than human HbA, the relative role of tertiary and quaternary allosteric effects may be better understood. The model systems illustrated below are the different hemoglobins from trout's blood, since they are characterized by the most striking variability of heterotropic effects, ranging from totally absent to very extreme with dominant contributions of tertiary effects. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Resonance Raman spectra of photodissociated hemoglobins: implications on cooperative mechanisms

Resonance Raman spectra at cryogenic temperatures of photodissociated hemoglobins and the corresponding deoxygenated preparations are compared and significant differences are found in modes with contributions from peripheral substituents of the heme as well as in the iron-histidine stretching mode. These differences in heme vibrational spectra reflect differences in the tertiary structure of the heme pocket between deoxyhemoglobin and the CO-bound form. An analysis of the effects of cooperative energy storage on the tertiary structure around the heme is made and used to interpret this resonance Raman data. The differences between the spectra of the deoxygenated preparations and the photoproducts provide evidence that a fraction of the free energy of cooperativity, DeltaG, is located away from the heme. These data support models for cooperativity in which the cooperative energy is distributed over many bonds or is localized in protein bonds only, such as those at the subunit interface. In addition, the local changes in amino acid positions, which must occur following the change in the state of ligand binding, may drive the changes in the structural relationships of the subunits and hence form one of the initial steps for triggering the quaternary structure transition.

Oximetry with the NMR signals of hemoglobin Val E11 and Tyr C7

The NMR visibility of the signals from erythrocyte hemoglobin (Hb) presents an opportunity to assess the vascular PO₂ (partial pressure of oxygen) in vivo to gather insight into the regulation of O₂ transport, especially in contracting muscle tissue. Some concerns, however, have arisen about the validity of using the Val E11 signal as an indicator of PO₂, since its intensity depends on tertiary structural changes, in contrast to the quaternary structure changes associated with relaxed (R) and tense (T) transition during O₂ binding. We have examined the Val E11 and Tyr C7 signal intensity as a function of Hb saturation by developing an oximetry system, which permits the comparative analysis of the NMR and spectrophotometric measurements. The spectrophotometric assay defines the Hb saturation level at a given PO₂ and yields standard oxygen-binding curves. Under defined PO₂ and Hb saturation values, the NMR measurements have determined that the Val E11 signal, as well as the Tyr C7 signal, tracks closely Hb saturation and can therefore serve as a vascular oxygen biomarker.

Protein dynamics explain the allosteric behaviors of hemoglobin

Bohr, Hasselbalch, and Krogh discovered homotropic and heterotropic allosteric behaviors of hemoglobin (Hb) in 1903/1904. A chronological description since then of selected principal models of the allosteric mechanism of Hb, such as the Adair scheme, the MWC two-state concerted model, the KNF induced-fit sequential model, the Perutz stereochemical model, the tertiary two-state model, and the global allostery model (an expanded MWC models), is concisely presented, followed by analysis and discussion of their limitations and deficiencies. The determination of X-ray crystallographic structures of deoxy- and ligated-Hb and the structure-based stereochemical model by Perutz are an epoch-making event in this history. However, his assignment of low-affinity deoxy- and high-affinity oxy-quaternary structures of Hb to the T- and R-states, respectively, though apparently reasonable, and as well as his hypothesis that the T-/R-quaternary structural transition regulates the oxygen-affinity, have created confusions and side-tracked studies of Hb on the structure-function relationship. The differences in static molecular structures of Hb between T(deoxy)- and R(oxy)-quaternary states reported in detail by Perutz and others are ligation-linked structural changes, but not related to the control/regulation of the oxygen-affinity. The oxygen-affinity (K(T) and K(R)) of Hb has been shown to be regulated by the heterotropic effector-linked tertiary structural changes without involving the T/R-quaternary changes. However, a recent high-resolution crystallographic analysis of Hb with different oxygen-affinities shows that static molecular structures of Hb determined by crystallography can neither identify the nature of the T(low-affinity) functional state nor decipher the mechanism by which Hb stores free energy in the T(low-affinity) functional state. Molecular dynamics simulations show that fluctuations of helices of oxy-Hb are increased upon de-oxygenation and/or binding 2,3-biphosphoglycerate. These are known to lower the oxygen-affinity of Hb. It is proposed that the coordination mode of the heme Fe with proximal and distal His is modulated by these helical fluctuations, resulting in the modulation of the oxygen-affinity of Hb. Therefore, it is proposed that the oxygen-affinity of Hb is regulated by pentanary (the 5th-order time-dependent or dynamic) tertiary structural changes rather than the T-/R-quaternary structural transitions in Hb. Homotropic and heterotropic allosteric effects of Hb are oxygen- and effector-linked, conformational entropy-driven entropy-enthalpy compensation phenomena and not much to do with static structural changes. The dynamic allostery model, which integrates these observations, provides the structural basis for the global allostery model (an expanded MWC model).

Band 3 (AE1, SLC4A1)-mediated transport of stilbenedisulfonates. III: Role of solute and protein structure in proton-activated stilbenedisulfonate influx

DBDS (4,4'-dibenzamido-2,2'-stilbenedisulfonate) influx into magnesium resealed ghosts (MRSG) occurs over the anion/proton co-transport pH range (pK approximately 5.0). Here, factors are studied which may influence the pH dependence of DBDS transport. Accumulation of various stilbenedisulfonate (SD) molecules was studied and found to be correlated with the hydrophobicity of the R-groups (Hansch factor), not protonation of the sulfonates. The role of proton binding to glutamate 681 was found not to be part of the rate-limiting step in DBDS uptake by MRSG. Finally, the pH dependence of changes in quaternary structure/conformational state was investigated using an assay involving photo-crosslinking of band 3 subunits in the presence of DASD (4,4'-diazido-2,2'-stilbenedisulfonate). Lowering the pH promoted intersubunit crosslinking by DASD, with a pK value of 4.75+/-1.0. This value is comparable to the pK for DBDS binding to the "second" class of sites on control band 3 (pK = 5.01+/-0.16), and to DBDS influx into control MRSG (pK values between 4.57+/-0.15 and 4.7+/-0.1). Susceptibility to photo-crosslinking was reversed by raising the pH prior to initiation of the reaction. Significantly, no photo-crosslinking was observed between pH 6.0 and 8.0, where band 3 subunits are known to exist as stable dimers and tetramers. We conclude that intersubunit photo-crosslinking does not simply involve random collision between photo-activated DASD and band 3. Rather, proton binding to band 3 either alters the conformation at the interface between subunits of pre-existing tetramers, or it promotes self-association of stable dimers to a "novel" tetrameric conformational state.

Quaternary structures of intermediately ligated human hemoglobin a and influences from strong allosteric effectors: resonance Raman investigation

The Fe-histidine stretching (nu(Fe-His)) frequency was determined for deoxy subunits of intermediately ligated human hemoglobin A in equilibrium and CO-photodissociated picosecond transient species in the presence and absence of strong allosteric effectors like inositol(hexakis)phosphate, bezafibrate, and 2,3-bisphosphoglycerate. The nu(Fe-His) frequency of deoxyHb A was unaltered by the effectors. The T-to-R transition occurred around m = 2-3 in the absence of effectors but m > 3.5 in their presence, where m is the average number of ligands bound to Hb and was determined from the intensity of the nu(4) band measured in the same experiment. The alpha1-beta2 subunit contacts revealed by ultraviolet resonance Raman spectra, which were distinctly different between the T and R states, remained unchanged by the effectors. This observation would solve the recent discrepancy that the strong effectors remove the cooperativity of oxygen binding in the low-affinity limit, whereas the (1)H NMR spectrum of fully ligated form exhibits the pattern of the R state.

The Molecular Code for Hemoglobin Allostery Revealed by Linking the Thermodynamics and Kinetics of Quaternary Structural Change. 1. Microstate Linear Free Energy Relations

A novel model linking the thermodynamics and kinetics of hemoglobin's allosteric (R --> T) and ligand binding reactions is applied to photolysis data for human HbCO. To describe hemoglobin's kinetics at the microscopic level of structural transitions and ligand-binding events for individual [ij]-ligation microstates ((ij)R --> (ij)T, (ij)R + CO --> ((i)(+1))(k)R, and (ij)T + CO --> ((i)(+1))(k)T), the model calculates activation energies, (ij)DeltaG(++), from previously measured cooperative free energies of the equilibrium microstates (Huang, Y., and Ackers, G. K. (1996) Biochemistry 35, 704-718) by using linear free energy relations ((ij)DeltaG(++) - (01)DeltaG(++) = alpha[(ij)DeltaG - (01)DeltaG], where the parameter alpha, describing the variation of activation energy with reaction energy perturbation, can depend on the natures of both the reaction and the perturbation). The alpha value measured here for the allosteric dynamics, 0.21 +/- 0.03, corresponds closely to values observed previously, strongly suggesting that the thermodynamic microstate energies directly underlie the allosteric kinetics (as opposed to the alpha((ij)DeltaG(RT)) serving merely as arbitrary fitting parameters). Besides systematizing the study of hemoglobin kinetics, the utility of the microstate linear free energy model lies in the ability to test microscopic aspects of allosteric dynamics such as the "symmetry rule" for quaternary change deduced previously from thermodynamic evidence (Ackers, G. K., et al. (1992) Science 255, 54-63). Reflecting a remarkably detailed correspondence between thermodynamics and kinetics, we find that a kinetic model that includes the large free energy splitting between doubly ligated T microstates implied by the symmetry rule fits the data significantly better than one that does not.

The Molecular Code for Hemoglobin Allostery Revealed by Linking the Thermodynamics and Kinetics of Quaternary Structural Change. 2. Cooperative Free Energies of (αFeCOβFe)2 and (αFeβFeCO)2 T-State Tetramers

Ligand photodissociation experiments are used to measure the prephotolysis equilibria between doubly liganded R and T quaternary conformers of the symmetric Fe-Co HbCO hybrids, (alpha(FeCO)beta(Co))(2) and (alpha(Co)beta(FeCO))(2). The free energies obtained from these data are used to calculate the cooperative free energies of the (alpha(FeCO)beta(Fe))(2) and (alpha(Fe)beta(FeCO))(2) intermediate CO-ligation states of normal hemoglobin in the T conformation, quantities important to the evaluation of current models of cooperativity. The symmetry rule model, incorporating sequential cooperativity of T-state ligand binding within an alphabeta dimer in addition to the traditional two-state cooperativity of the tetramer, predicts a larger free energy penalty for disturbing both dimers in a doubly liganded T tetramer than would be expected in the two-state model as currently formulated. (Cooperative energy penalties are simply proportional to the number of tetramer-bound ligands in the traditional two-state model.) The value found here for the energies of doubly liganded T microstates in which both dimers are perturbed, 7.9 +/- 0.3 kcal/mol, is consistent with the symmetry rule model but significantly higher than that expected (5-6 kcal/mol) in the two-state model of cooperativity.

Linus Pauling and sickle cell disease

The 1949 paper by Linus Pauling et al. [Science 110 (1949) 543-548] describing the discovery of sickle cell anemia as the first molecular disease had a major impact on biology and medicine. Inspired by the scholarly works of John Edsall on the history of hemoglobin research, I present a brief retrospective analysis of Pauling's paper. This is followed by some personal recollections of Edsall and Pauling.

A tertiary two-state allosteric model for hemoglobin

The two-state allosteric model of Monod, Wyman, and Changeux (MWC) provides an excellent description of homotropic effects in a vast array of equilibrium and kinetic measurements on cooperative ligand binding by hemoglobin. However, in contrast to experimental observations, the model does not allow for alteration of the ligand affinity of the T quaternary structure by allosteric effectors. This failure to explain heterotropic effects has been appreciated for over 30 years, and it has been generally assumed to result from tertiary conformational changes in the absence of a quaternary change. Here we explore a model that preserves the essential MWC idea that binding without a quaternary conformational change is non-cooperative, but where tertiary conformations of individual subunits play the primary role instead of the quaternary conformations. In this model, which we call the 'tertiary two-state (TTS) model', the two affinity states correspond to two tertiary conformations of individual subunits rather than the two quaternary conformations of the MWC two-state allosteric model. Ligation and the R quaternary structure bias the subunit population toward the high affinity tertiary conformation, while deligation and the T quaternary structure favor the low affinity tertiary conformation. We show that the model is successful in quantitatively explaining a demanding set of kinetic data from nanosecond carbon monoxide photodissociation experiments at times longer than approximately 1 micros. Better agreement between the model and the submicrosecond kinetic data may result from detailed considerations of the distribution and dynamics of conformational substates of the two tertiary conformations. The model is consistent with the results of solution, gel, and single crystal oxygen binding studies, but underestimates the population of doubly-liganded molecules determined in low-temperature electrophoresis experiments.

High and low oxygen affinity conformations of T state hemoglobin

To understand the interplay between tertiary and quaternary transitions associated with hemoglobin function and regulation, oxygen binding curves were obtained for hemoglobin A fixed in the T quaternary state by encapsulation in wet porous silica gels. At pH 7.0 and 15 degrees C, the oxygen pressure at half saturation (p50) was measured to be 12.4 +/- 0.2 and 139 +/- 4 torr for hemoglobin gels prepared in the absence and presence of the strong allosteric effectors inositol hexaphosphate and bezafibrate, respectively. Both values are in excellent agreement with those found for the binding of the first oxygen to hemoglobin in solution under similar experimental conditions. The corresponding Hill coefficients of hemoglobin gels were 0.94 +/- 0.02 and 0.93 +/- 0.03, indicating, in the frame of the Monod, Wyman, and Changeux model, that high and low oxygen-affinity tertiary T-state conformations have been isolated in a pure form. The values, slightly lower than unity, reflect the different oxygen affinity of alpha- and beta-hemes. Significantly, hemoglobin encapsulated in the presence of the weak effector phosphate led to gels that show intermediate oxygen affinity and Hill coefficients of 0.7 to 0.8. The heterogeneous oxygen binding results from the presence of a mixture of the high and low oxygen-affinity T states. The Bohr effect was measured for hemoglobin gels containing the pure conformations and found to be more pronounced for the high-affinity T state and almost absent for the low-affinity T state. These findings indicate that the functional properties of the T quaternary state result from the contribution of two distinct, interconverting conformations, characterized by a 10-fold difference in oxygen affinity and a different extent of tertiary Bohr effect. The very small degree of T-state cooperativity observed in solution and in the crystalline state might arise from a ligand-induced perturbation of the distribution between the high- and low-affinity T-state conformations.

Sol-gel trapping of functional intermediates of hemoglobin: geminate and bimolecular recombination studies

The encapsulation of proteins in porous sol-gels is a promising technique for generating, trapping, and probing functionally significant nonequilibrium protein species. An essential step needed in the pursuit of that goal is establishing the degree to which the sol-gel limits conformational change upon adding or removing substrates. In the present study, geminate recombination and solvent phase bimolecular recombination of CO to human adult hemoglobin (HbA) are used as sensitive probes of the degree of conformational constraint within the sol-gel. Two forms of CO saturated encapsulated HbA are generated. In one case, designated [COHbA], the equilibrium form of COHbA is directly encapsulated. In the second case, designated as [deoxyHbA] + CO, the equilibrium form of deoxyHbA is encapsulated and only after the sample has aged is CO introduced to the HbA through the porous sol-gel matrix. Three different preparative protocols are used to generate the sol-gels for each of the two forms of encapsulated COHbA. The kinetic traces obtained from these encapsulated samples allow for an easy evaluation of the extent to which the sol-gel is locking in the initial tertiary/quaternary structure. The results show that the sol-gel encapsulated samples can be used with pulsed laser sources and that one of the tested encapsulation protocols is far superior with respect to conformational locking. This protocol is used to trap and probe nonequilibrium forms such as the liganded T state of HbA, a species whose properties are needed to fully explore allostery in HbA.

Site-directed mutagenesis in hemoglobin: test of functional homology of the F9 amino acid residues of hemoglobin alpha and beta chains

The cysteine residue at F9(93) of the human hemoglobin (Hb A) beta chain, conserved in mammalian and avian hemoglobins, is located near the functionally important alpha1-beta2 interface and C-terminal region of the beta chain and is reactive to sulfhydryl reagents. The functional roles of this residue are still unclear, although regulation of local blood flow through allosteric S-nitrosylation of this residue is proposed. To clarify the role of this residue and its functional homology to F9(88) of the alpha chain, we measured oxygen equilibrium curves, UV-region derivative spectra, Soret-band absorption spectra, the number of titratable -SH groups with p-mercuribenzoate and the rate of reaction of these groups with 4, 4'-dipyridine disulfide for three recombinant mutant Hbs with single amino acid substitutions: Ala-->Cys at 88alpha (rHb A88alphaC), Cys-->Ala at 93beta (rHb C93betaA) and Cys-->Thr at 93beta (rHb C93betaT). These Hbs showed increased oxygen affinities and impaired allosteric effects. The spectral data indicated that the R to T transition upon deoxygenation was partially restricted in these Hbs. The number of titratable -SH groups of liganded form was 3.2-3.5 for rHb A88alphaC compared with 2.2 for Hb A, whereas those for rHb C93betaA and rHb C93betaT were negligibly small. The reduction of rate of reaction with 4,4'-dipyridine disulfide upon deoxygenation in rHb A88alphaC was smaller than that in Hb A. Our experimental data have shown that the residues at 88alpha and 93beta have definite roles but they have no functional homology. Structure-function relationships in our mutant Hbs are discussed.

Functional analysis of hemoglobin molecules locked in doubly liganded conformations

A controversy still exists over whether the molecular basis of hemoglobin cooperativity can be more appropriately explained by one of two classic allosteric models, the concerted and sequential models. To distinguish these two models from the viewpoint of their fundamental processes, namely, the presence or absence of conformational equilibria, we have trapped the conformations of nickel(II)-iron(II) hybrid hemoglobin molecules with two CO-bound, alpha2(Fe-CO)beta2(Ni) and alpha2(Ni)beta2(Fe-CO), by encapsulation in the water-filled pores of sol-gel-derived transparent silica-gels. In our experimental system, nickel(II) protoporphyrin binds neither O2 nor CO and mimics a fixed deoxyheme, and the gel matrix provides a means of inhibiting large-scale protein structural changes, thus enabling O2 equilibrium study of the hybrids still in their doubly liganded conformations. Results showed that two conformations of widely different O2 affinity exist together in each doubly liganded hemoglobin, providing a direct proof of the concerted mechanism versus the sequential mechanism.

Quaternary structure sensitive tyrosine residues in human hemoglobin: UV resonance raman studies of mutants at alpha140, beta35, and beta145 tyrosine

Recent studies noted the contribution of alpha42Tyr to the T-R-dependent UV resonance Raman (UVRR) spectral changes of HbA [Nagai, M., et al. (1996) J. Mol. Struct. 379, 65-75; Huang, S., et al. (1997) Biochemistry 36, 6197-6206], but the observed UVRR changes of the Tyr residue cannot be fully interpreted with alpha42Tyr alone. To identify the remaining contributions, the 235 nm-excited UVRR spectra of Tyr mutant Hbs at alpha140, beta35, and beta145 were investigated here. The Fe-His stretching mode demonstrated that all of these mutant Hbs take the T structure in the deoxy form under these experimental conditions. The UVRR change of the Trp residue of these mutants upon the T-R transition was the same as that in HbA, indicating that the T-R-dependent UVRR change of beta37Trp is not due to stacking with Tyr residues but is due to the formation or destruction of a hydrogen bond. The recombinant Hbs beta35Tyr --> Phe and beta35Tyr --> Thr both exhibited UVRR spectra identical with that of HbA, meaning that beta35Tyr is not responsible. In the spectra of des(beta146His,beta145Tyr)Hb with inositol hexaphosphate, the frequency shift of the Tyr RR bands was the same as that in HbA but the intensity enhancement in the CO form was small, suggesting that beta145Tyr contributes to a part of the intensity change, but scarcely relates to the frequency shift. In the spectra of Hb Rouen (alpha140Tyr --> His), the frequency shifts of bands at 1617 (Y8a) and 1177 (Y9a) cm-1 following ligation were half of those in HbA, while the intensity enhancement was not detected. This result means that alpha140Tyr is responsible for both the frequency shift and the intensity changes. It is suggested that the frequency shift of the Tyr RR bands upon the T --> R transition is due to changes in the hydrogen bonding state of alpha42- and alpha140Tyr and that the intensity enhancement is due to changes in the environment of the penultimate Tyr in both alpha and beta subunits (alpha140 and beta145). These alterations in the vibrational spectra clearly demonstrate which tyrosine residues are involved in the T-R transition as a result of modification of their local environments.

Fluorescent effector as a probe of the allosteric equilibrium in methemoglobin

A fluorescent analogue of diphosphoglycerate (DPG), hydroxy-pyrenetrisulfonate (HPT), was used as a probe of the allosteric equilibrium of methemoglobin. Like DPG, HPT binds, one per tetramer, with a higher affinity to deoxyHb than to oxyHb. Once bound, the HPT fluorescence is quenched by energy transfer to the hemes. HPT can thus serve as a probe of the conformational state of the hemoglobin tetramer: a higher quenching indicates a stronger binding and therefore, more of the deoxy conformation. Since HPT binds to the same site as DPG, it can be displaced by DPG in order to determine the fluorescence intensity of the free HPT under the same conditions, to correct for the inner filter effect. The high spin ferric ligands, such as water and F, showed less fluorescence (more of the deoxy state) than low spin cyano-metHb. The aquo-metHb samples showed a reversion to the oxyHb conformation above pH 7, as expected due to the acid-alkaline transition forming hydroxy-metHb. Effectors such as bezafibrate, which do not bind to the same site as DPG, show an increase in the deoxy-like characteristics.

T state hemoglobin binds oxygen noncooperatively with allosteric effects of protons, inositol hexaphosphate, and chloride

Hemoglobin is the paradigm of allosteric proteins. Over the years, cooperative oxygen binding has been explained by different models predicting that the T state of hemoglobin binds oxygen either noncooperatively or with some degree of cooperativity or with strong cooperativity. Therefore, a critical test that discriminates among models is to determine the oxygen binding by the T state of hemoglobin. Fixation of hemoglobin in the T state has been achieved either by crystallization from polyethylene glycol solutions or by encapsulation in wet porous silica gels. Hemoglobin crystals bind oxygen noncooperatively with reduced affinity compared with solution, with no Bohr effect and with no influence of other allosteric effectors. In this study, we have determined accurate oxygen-binding curves to the T state of hemoglobin in silica gels with the same microspectrophotometric apparatus and multiwavelengths analysis used in crystal experiments. The T state of hemoglobin in silica gels binds oxygen noncooperatively with an affinity and a Bohr effect similar to those observed in solution for the binding of the first oxygen molecule. Other allosteric effectors such as inositol hexaphosphate, bezafibrate, and chloride significantly affect oxygen affinity. Therefore, T state hemoglobins that are characterized by strikingly different functional properties share the absence of cooperativity in the binding of oxygen. These findings are fully consistent with the Monod, Wyman, and Changeux model and with most features of Perutz's stereochemical model, but they are not consistent with models of both Koshland and Ackers.

Can a two-state MWC allosteric model explain hemoglobin kinetics?

We have analyzed the nanosecond-millisecond kinetics of ligand binding and conformational changes in hemoglobin. The kinetics were determined from measurements of precise time-resolved optical spectra following nanosecond photodissociation of the heme-carbon monoxide complex. To fit the data, it was necessary to extend the two-state allosteric model of Monod, Wyman, and Changeux (MWC) to include geminate ligand rebinding and nonexponential tertiary relaxation within the R quaternary structure. Considerable simplification of the model is obtained by using a linear free energy relation for the rates of quaternary transitions, and by incorporating concepts from recent studies on the physics of geminate rebinding and conformational changes in myoglobin. The model, described by 85 coupled differential equations, quantitatively explains a demanding set of complex kinetic data. Moreover, with the same set of kinetic parameters it simultaneously fits the equilibrium data on ligand binding and the distribution of ligation states. The present results, together with those from single-crystal oxygen binding studies, indicate that the two-state MWC allosteric model has survived its most critical tests.

Allosteric effectors do not alter the oxygen affinity of hemoglobin crystals

In solution, the oxygen affinity of hemoglobin in the T quaternary structure is decreased in the presence of allosteric effectors such as protons and organic phosphates. To explain these effects, as well as the absence of the Bohr effect and the lower oxygen affinity of T-state hemoglobin in the crystal compared to solution, Rivetti C et al. (1993a, Biochemistry 32:2888-2906) suggested that there are high- and low-affinity subunit conformations of T, associated with broken and unbroken intersubunit salt bridges. In this model, the crystal of T-state hemoglobin has the lowest possible oxygen affinity because the salt bridges remain intact upon oxygenation. Binding of allosteric effectors in the crystal should therefore not influence the oxygen affinity. To test this hypothesis, we used polarized absorption spectroscopy to measure oxygen binding curves of single crystals of hemoglobin in the T quaternary structure in the presence of the "strong" allosteric effectors, inositol hexaphosphate and bezafibrate. In solution, these effectors reduce the oxygen affinity of the T state by 10-30-fold. We find no change in affinity (< 10%) of the crystal. The crystal binding curve, moreover, is noncooperative, which is consistent with the essential feature of the two-state allosteric model of Monod J, Wyman J, and Changeux JP (1965, J Mol Biol 12:88-118) that cooperative binding requires a change in quaternary structure. Noncooperative binding by the crystal is not caused by cooperative interactions being masked by fortuitous compensation from a difference in the affinity of the alpha and beta subunits. This was shown by calculating the separate alpha and beta subunit binding curves from the two sets of polarized optical spectra using geometric factors from the X-ray structures of deoxygenated and fully oxygenated T-state molecules determined by Paoli M et al. (1996, J Mol Biol 256:775-792).

Allosteric proteins after thirty years: the binding and state functions of the neuronal alpha 7 nicotinic acetylcholine receptors

A key statement of the 1965 Monod-Wyman-Changeux (MWC) model for allosteric proteins concerns the distinction between the ligand-binding function (Y) and the relevant state function (R). Sequential models predict overlapping behavior of the two functions. In contrast, a straightforward experimental consequence of the MWC model is that for an oligomeric protein the parameters which characterize the two functions should differ significantly. Two situations, where R > Y and the system is hyper-responsive or where R < Y and the system is hypo-responsive, have been encountered. Indeed, the hyper-responsive pattern was first observed for the enzyme aspartate transcarbamoylase, by comparing Y with R monitored by a change in sedimentation. Extensions of the theory to ligand-gated channels led to the suggestion that, on the one hand, hyper-responsive properties also occur with high-affinity mutants. On the other hand, native channels of the acetylcholine neuronal alpha 7 receptor and low-affinity mutants of the glycine receptor can be interpreted in terms of the hypo-responsive pattern. For the ligand-gated channels, whereas R is detected directly by ion flux, ligand binding has rarely been measured and the formation of desensitized states may complicate the analysis. However, stochastic models incorporating both binding and channel opening for single molecules predict differences that should be measurable with new experimental approaches, particularly fluorescence correlation spectroscopy.

Allosteric intermediates in hemoglobin. 2. Kinetic modeling of HbCO photolysis

Nanosecond absorption spectra are measured in the Soret and near-UV spectral regions of human hemoglobin (Hb) after laser photolysis of the carbonyl adduct in order to study the dynamics of globin tertiary and quaternary conformational changes. Spectra and concentrations of physical intermediates, distinguished by extent of heme ligation and intraprotein relaxation, are obtained from a global analysis using a microscopic kinetic model that explicitly accounts for six observed relaxation and recombination processes. Three observed rate constants for CO rebinding and two intraprotein relaxation constants obtained are similar to constants determined by Hofrichter et al. [(1983) Proc. Natl. Acad. Sci. U.S.A. 80, 2235], the latter two comprising the 20-30-microseconds R --> T quaternary transition and a previously unassigned 1-microseconds intraprotein relaxation. On the basis of the modeled intermediate spectra, as well as UV circular dichroism results observed on this time scale [Björling, S.C., Goldbeck, R.A., Paquette, S.J., Milder, S.J., & Kliger, D.S. (1996) Biochemistry 35, 8619-8627], the 1-microsecond relaxation is assigned to heme conformational changes concomitant with a relaxation of protein conformation at the alpha 1 beta 2 interface corresponding to an initial step in a compound R --> T reaction path.

Modulated excitation of singly ligated carboxyhemoglobin

We have extended the method of modulated excitation, a small perturbation kinetic method, to study ligand binding and conformational change of hemoglobin tetramers with a single ligand bound. To restrict the excitation to the first ligand, only 1% of the hemes have bound CO, and the remainder are kept unliganded. A detailed theory is presented which agrees well with the experimental observations. This method of observing ligand recombination also provides a novel and simple method for determination of hemoglobin concentration. Additional relaxation processes are also observed. By fitting independently determined spectra to the spectra associated with the relaxations, these processes are assigned as thermal excitation and thermally driven protonation/deprotonation reactions. These added relaxations arise from the deoxy-Hb portion of the samples, and demonstrate that modulated excitation can be used effectively for temperature perturbation in the absence of photodissociation. The spectra observed are not well described by the spectra of allosteric change, however, and we conclude that there is no significant mixing of quaternary states at the first ligation step. In an appendix we present a derivation of the particular features seen in thermally modulated protonation reactions.