Testing the two-state model: anomalous effector binding to human hemoglobin

MWC allosteric model explains unusual hemoglobin-oxygen binding curves from sickle cell drug binding

An oxygen-affinity-modifying drug, voxelotor, has very recently been approved by the FDA for treatment of sickle cell disease. The proposed mechanism of action is by preferential binding of the drug to the R quaternary conformation, which cannot copolymerize with the T conformation to form sickle fibers. Here, we report widely different oxygen dissociation and oxygen association curves for normal blood in the presence of voxelotor and interpret the results in terms of the allosteric model of Monod, Wyman, and Changeux with the addition of drug binding. The model does remarkably well in quantitatively explaining a complex data set with just the addition of drug binding and dissociation rates for the R and T conformations. Whereas slow dissociation of the drug from R results in time-independent dissociation curves, the changing association curves result from slow dissociation of the drug from T, as well as extremely slow binding of the drug to T. By calculating true equilibrium curves from the model parameters, we show that there would be a smaller decrease in oxygen delivery from the left shift in the dissociation curve caused by drug binding if drug binding and dissociation for both R and T were rapid. Our application of the Monod, Wyman, and Changeux model demonstrates once more its enormous power in explaining many different kinds of experimental results for hemoglobin. It should also be helpful in analyzing oxygen binding and in vivo delivery in future investigations of oxygen-affinity-modifying drugs for sickle cell disease.

Molecular level probing of preferential hydration and its modulation by osmolytes through the use of pyranine complexed to hemoglobin

Two spectroscopic probes are used to expose molecular level changes in hydration shell water interactions that directly relate to such issues as preferential hydration and protein stability. The major focus of the present study is on the use of pyranine (HPT) fluorescence to probe as a function of added osmolytes (PEG, urea, trehalose, and magnesium), the extent to which glycerol is preferentially excluded from the hydration shell of free HPT and HPT localized in the diphosphoglycerate (DPG) binding site of hemoglobin in both solution and in Sol-Gel matrices. The pyranine study is complemented by the use of vibronic side band luminescence from the gadolinium cation that directly exposes the changes in hydrogen bonding between first and second shell waters as a function of added osmolytes. Together the results form the basis for a water partitioning model that can account for both preferential hydration and water/osmolyte-mediated conformational changes in protein structure.

Control of the allosteric equilibrium of hemoglobin by cross-linking agents

The kinetics of ligand rebinding have been studied for modified or cross-linked hemoglobins (Hbs). Several compounds were tested that interact with alpha Val 1 or involve a cross-link between alpha Val 1 and alpha Lys 99 of the opposite dimer. By varying the length of certain cross-linking molecules, a wide range in the allosteric equilibrium could be obtained. Several of the mono-aldehyde modified Hbs show a shift toward the high affinity conformation of Hb. At the other extreme, for certain di-aldehyde cross-linked Hbs, the CO kinetics are typical of binding to deoxy Hb, even at low photodissociation levels, with which the dominant photoproduct is the triply liganded species; in these cases the hemoglobin does not switch from the low to high affinity state until after the fourth ligand is bound. Although each modified Hb shows only two distinct rates, the kinetic data as a function of dissociation level cannot be simulated with a simple two-state model. A critical length is observed for the maximum shift toward the low affinity T-state. Longer or shorter lengths of the cross-linker yielded more high affinity R-state. Unlike native Hb, which is in equilibrium with free dimers, the cross-linked Hbs maintain the fraction slow kinetics, which is unique to Hb tetramers, even at 0.5 microM (total heme). Addition of HbCN to unmodified HbCO solutions results in dimer exchange, which decreases the relative fraction of slow bimolecular kinetics; the cross-linked Hbs did not show such an effect, indicating that they do not participate in dimer exchange.

Heterotropic effectors exert more significant strain on monoligated than on unligated hemoglobin

The effect of allosteric effectors, such as inositol hexakisphosphate and/or bezafibrate, has been investigated on the unliganded human adult hemoglobin both spectroscopically (employing electronic absorption, circular dichroism, resonance Raman, and x-ray absorption near-edge spectroscopies) and functionally (following the kinetics of the first CO binding step up to a final 4% ligand saturation degree). All data indicate that the unliganded T-state is not perturbed by the interaction with either one or both effectors, suggesting that their functional influence is only exerted when a ligand molecule is bound to the heme. This is confirmed by the observation that CO dissociation from partially liganded hemoglobin ( </= 0.04) is strongly altered by the presence of either effector, and the effect is enhanced whenever the two effectors are simultaneously present. Altogether, these data are a direct demonstration of the occurrence of a strain induced by the presence of a ligand molecule bound to the heme, and for the first time there is a clear indication that the expression of the functional heterotropic effect by these non-heme ligands requires this strain, which is not present in the unliganded molecule.

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.

A two-state analysis of co-operative oxygen binding in the three human embryonic haemoglobins

The binding of oxygen to the three human embryonic haemoglobins, at pH 7.4, has been shown to occur as a co-operative process. Analysis of oxygen-binding curves obtained in the absence of organic phosphate allosteric effectors shows that the process can be described quite accurately by the two-state model of allosteric action. In the presence of organic phosphates, the binding affinity for oxygen to the T-state of the alpha 2 epsilon 2 and zeta 2 epsilon 2 haemoglobins is significantly lowered. The values of the best-fit two-state parameters determined for each of the embryonic haemoglobins together with the temperature-dependence of the overall equilibrium binding process are discussed in terms of oxygen transfer from the maternal blood supply. Fast-reaction studies have been used to determine the rate constants of the oxygen association and dissociation processes occurring in the R-state and the rate of the allosteric R > T conformational transition. Analysis of these data suggests a likely reason for the high affinity and low co-operativity of the embryonic proteins and identifies the origins of the inability of equilibrium measurements to identify chain non-equivalence in the R-state.

Probing the hemoglobin central cavity by direct quantification of effector binding using fluorescence lifetime methods

Time-resolved fluorescence methods have been used to show that 8-hydroxy-1,3,6-pyrenetrisulfonate (HPT), a fluorescent analog of 2,3-diphosphoglycerate, binds to the central cavity of carboxyhemoglobin A (HbACO) at pH 6.35. A direct quantitative approach, based on the distinctive free and bound HPT fluorescent lifetimes of 5.6 ns and approximately 27 ps, respectively, was developed to measure the binding affinity of this probe. HPT binds to a single site and is displaced by inositol hexaphosphate at a 1:1 mol ratio, indicating that binding occurs at the 2,3-diphosphoglycerate site in the central cavity. Furthermore, the results imply that low pH HbACO exists as an altered R state and not an equilibrium mixture of R and T states. The probe was also used to monitor competitive effector binding and to compare the affinity of the binding site in several cross-bridged HbA derivatives.

Quaternary structure dynamics and carbon monoxide binding kinetics of hemoglobin valency hybrids

The kinetics of CO binding and changes in quaternary structure for symmetric valency hybrids of human hemoglobin have been extensively studied by laser photolysis techniques. Both alpha+beta and alpha beta+ hybrids were studied with five different ferric ligands, over a broad range of CO concentrations and photolysis levels. After full CO photolysis, the hybrid tetramers switch extensively and rapidly (< 200 microseconds) to the T quaternary structure. Both R --> T and T --> R transition rates for valency hybrid tetramers with 0 and 1 bound CO have been obtained, as well as the CO association rates for alpha and beta subunits in the R and T states. The results reveal submillisecond R reversible T interconversion, and, for the first time, the changes in quaternary rates and equilibria due to binding a single CO per tetramer have been resolved. The data also show significant alpha-beta differences in quaternary dynamics and equilibria. The allosteric constants do not vary with the spin states of the ferric subunits as predicted by the Perutz stereochemical model. For the alpha beta+CN hybrid the kinetics are heterogeneous and imply partial conversion to a T-like state with very low (seconds) R reversible T interconversion.

Oxygen and CO binding to triply NO and asymmetric NO/CO hemoglobin hybrids

The bimolecular and geminate CO recombination kinetics have been measured for hemoglobin (Hb) with over 90% of the ligand binding sites occupied by NO. Since Hb(NO)4 with inositol hexaphosphate (IHP) at pH below 7 is thought to take on the low affinity (deoxy) conformation, the goal of the experiments was to determine whether the species IHPHb-(NO)3(CO) also exists in this quaternary structure, which would allow ligand binding studies to tetramers in the deoxy conformation. For samples at pH 6.6 in the presence of IHP, the bimolecular kinetics show only a slow phase with rate 7 x 10(4) M-1 s-1, characteristic of CO binding to deoxy Hb, indicating that the triply NO tetramers are in the deoxy conformation. Unlike Hb(CO)4, the fraction recombination occurring during the geminate phase is low (< 1%) in aqueous solutions, suggesting that the IHPHb(NO)3(CO) hybrid is also essentially in the deoxy conformation. By mixing stock solutions of HbCO and HbNO, the initial exchange of dimers produces asymmetric (alpha NO beta NO/alpha CO beta CO) hybrids. At low pH in the presence of IHP, this hybrid also displays a high bimolecular quantum yield and a large fraction of slow (deoxy-like) CO recombination; the slow bimolecular kinetics show components of equal amplitude with rates 7 and 20 x 10(4) M-1 s-1, probably reflecting the differences in the alpha and beta chains. Samples of symmetric hybrids (a2NOI32Co or a2Co922NO) showed a lower (R-like) bimolecular yield and less slow phase for the CO bimolecular recombination, relative to the asymmetric hybrid or the triply NO species. The slower (T state) bimolecular rate of 7 x 104 M-1 s-1 was observed for CO rebinding to a chain.While oxygen equilibrium studies with 'HPHb(NO)3 were hampered by a high oxidation rate, it was possible to perform experiments with samples equilibrated with a mixed CO/oxygen atmosphere. Photodissociation of CO allows a temporary exposure of the binding sites to oxygen. The results confirm that IHPHb(NO)3 has a low oxygen affinity.

Allosteric kinetics and equilibria of triligated, cross-linked hemoglobin

Using modulated excitation, we have measured the forward and reverse rates of the allosteric transition between relaxed (R) and tense (T) quaternary structures for triply ligated hemoglobin (Hb), cross-linked between the alpha chains at Lys 99. Oxygen, carbon monoxide, and water were used as ligands and were studied in phosphate and low Cl- bis-Tris buffers at neutral pH. Since the cross-link prohibits disproportionation, triply ligated aquomet Hb species with ferrous beta chains were specifically isolated by isoelectric focusing. Modulated excitation provides rate pairs and therefore gives equilibrium constants between quaternary structures. To coordinate with that information, oxygen binding curves of fully ferrous and tri-aquomet Hb were also measured. L3, the equilibrium constant between three liganded R and T structures, is determined by modulated excitation to be of order unity for O2 or CO (1.1 to 1.5 for 3O2 and 0.7 for 3CO bound), while with three aquomet subunits it is much greater (> or = 23). R-->T conversion rates are similar to those found for HbA, with weak sensitivity to changes in L3. The L3 values from HbXL O2 were used to obtain a unique allosteric decomposition of the ferrous O2 binding curve in terms of KT, KR, and L3. From these values and the O2 binding curve of tri-aquomet HbXL, L3 was calculated to be 2.7 for the tri-aquomet derivative. Consistency in L3 values between equilibrium and modulated excitation data for tri-aquomet-HbXL can be achieved if the equilibrium constant for O2 binding to the alpha chains is six times lower than that for binding to the beta chains in the R state, while the cooperative properties remain homogeneous. The results are in quantitative agreement with other studies, and suggest that the principal effect of the cross-link is to decrease the R state and T state affinity of the alpha subunits with almost no change in the affinity of the beta subunits, leaving the allosteric parameters L and c unchanged.

Photoselection in polarized photolysis experiments on heme proteins

Polarized photolysis experiments have been performed on the carbon monoxide complex of myoglobin to assess the effects of photoselection on the kinetics of ligand rebinding and to investigate the reorientational dynamics of the heme plane. The results are analyzed in terms of the optical theory developed in the preceding paper by Ansari and Szabo. Changes in optical density arising from rotational diffusion of the photoselected population produce large deviations from the true geminate ligand rebinding curves if measurements are made with only a single polarization. The apparent ligand rebinding curves are significantly distorted even at photolysis levels greater than 90%. These deviations are eliminated by obtaining isotropically-averaged optical densities from measurements using both parallel and perpendicular polarizations of the probe pulse. These experiments also yield the optical anisotropy, which gives a novel method for accurately determining the degree of photolysis, as well as important information on the reorientational dynamics of the heme plane. The correlation time for the overall rotational diffusion of the molecule is obtained from the decay of the anisotropy. The anisotropy prior to rotational diffusion is lower than that predicted for a rigidly attached, perfectly circular absorber, corresponding to an apparent order parameter of S = 0.95 +/- 0.02. Polarized absorption data on single crystals suggest that the decreased anisotropy results more from internal motions of the heme plane which take place on time scales shorter than the duration of the laser pulse (10 ns) than from out-of-plane polarized transitions.

Speed of intersubunit communication in proteins

To determine the speed of communication between protein subunits, time-resolved absorption spectra were measured following partial photodissociation of the carbon monoxide complex of hemoglobin. The experiments were carried out using linearly polarized, 10-ns laser pulses, with the polarization of the excitation pulse both parallel and perpendicular to the polarization of the probe pulse. The substantial contribution to the observed spectra from photoselection effects was eliminated by isotropically averaging the polarized spectra, allowing a detailed comparison of the kinetics as a function of the degree of photolysis. These results show that prior to 1 microsecond both geminate ligand rebinding and conformational relaxation are independent of the number of ligands dissociated from the hemoglobin tetramer, as expected for a two-state allosteric model. After this time the kinetics depend on the ligation state of the tetramer. The conformational relaxation at 10 microseconds can be interpreted in terms of the two-state allosteric model as arising from the R to T quaternary conformational change of both unliganded and singly liganded molecules. These results suggest that communication between subunits requires about 1 microsecond and that the mechanism of the communication which occurs after this time is via the R to T conformational change. The optical anisotropy provides a novel means of accurately determining the extinction coefficients of the transient photoproduct. The decay in the optical anisotropy, moreover, provides an accurate determination of the rotational correlation time of 36 +/- 3 ns.

Application of linear free energy relations to protein conformational changes: the quaternary structural change of hemoglobin

The transition state for the R in equilibrium with T quaternary conformational change of hemoglobin has thermodynamic properties much closer to those of the R conformation than to those of the T conformation. This finding is based on a comparison of activation and equilibrium enthalpy and entropy changes and on the observation of a linear free energy relationship between quaternary rate and equilibrium constants. A previous theoretical study [Janin, J. & Wodak, S. J. (1985) Biopolymers 24, 509-526], using a highly simplified energy function, suggests that the R-like transition state is the result of a reaction pathway with the maximum buried surface area between alpha beta dimers.

Fitting abnormal oxygen equilibrium curves of hemoglobin

A growing number of oxygen equilibrium curves for hemoglobin (Hb) mutants, post-translational modifications, or the binding of potent new effectors of Hb cannot be fitted adequately with the two-state model. Examples are curves showing double maxima in the derivative of the Hill plot, or slopes of less than unity. We present such examples of modified hemoglobins and strong effectors in this study and calculate at which substate level the two-state model differs from the data. Analysis of hemoglobin oxygen equilibrium curves is reconsidered using the two-state model extended to allow variation of the individual substate probabilities. In this way the effect on the equilibrium due to perturbations in energy of each substate can be studied as a diagnostic tool.

Structure dynamics of the hemoglobin mutants Hb Hôtel Dieu, HbG Philadelphia, HbJ Mexico, Hb St. Mandé and Hb San Diego, studied by nanosecond-laser-flash photolysis

The kinetics of the change from the carboxy to the deoxy conformation of the mutated hemoglobins mentioned in the title and of normal human adult hemoglobin were determined from measurements of light absorption changes occurring up to 50 microseconds after nanosecond-laser photodissociation of the corresponding CO complexes. The spectral evolution of the mutated hemoglobins was found to be similar in its main features to that of normal hemoglobin. The kinetics could be decomposed into two phases with rates 1.1-1.8 x 10(6) s-1 and 0.17-0.34 x 10(6) s-1 (except Hb St. Mandé which displayed only the faster phase). Study of the mutated subunits of HbJ Mexico (alpha subunit) and Hb Hôtel Dieu (beta subunit) showed that they convert exponentially to the stable deoxy state after photodeligation at the same rates as the corresponding subunits of normal Hb: 1.1 x 10(6) s-1 (alpha) and 0.3 x 10(6) s-1 (beta). The results indicate that there is no direct correlation between the kinetics of spectral relaxation in the time range studied and the oxygenation properties for these hemoglobins. However, there is some indication that the kinetics are dependent upon the region of mutation.

Rate of allosteric change in hemoglobin measured by modulated excitation using fluorescence detection

We have measured the forward and reverse rates of the allosteric transition of hemoglobin A with three CO molecules bound by using modulated excitation coupled with fluorescence quenching of the DPG analogue, PTS (8-hydroxy-1,3,6 pyrene trisulfonic acid). This dye is observed to bind to the T state with significantly larger affinity than to the R state, and thus provides an unequivocal marker for the molecule's conformational change. The allosteric rates obtained with the fluorescent dye (pH 7.0, bis-Tris buffer) are (3.4 +/- 1.0) x 10(3)s-1 for the R to T transition and (2.1 +/- 0.5) x 10(4)s-1 for the T to R transition. This gives an equilibrium constant L3 of 0.16 +/- 0.06. These results provide good agreement with modulated difference spectra calibrated from model compounds, arguing that there is little if any difference in the kinetics observed by the heme spectra and the kinetics of the full subunit motion. The equilibrium constant between structures (L3) is smaller in the absence of phosphates than observed in phosphate buffer (0.33). However, the rates of the allosteric transition increase in the absence of phosphates as compared with the corresponding rates in phosphate buffer of 1.0 x 10(3)s-1 and 3.0 x 10(3)s-1. The effects of inorganic phosphates on the equilibrium can be separated from the effects on kinetics. We find that phosphates also affect the dynamic behavior of hemoglobin, and the presence of 0.15 M phosphate can be viewed as raising the transition state energy between R and T conformations by approximately 0.5 kcal/mol exclusive of the T state stabilization. Dissociation constants for the dye were measured to be 104 +/- 25 microM for unligated T state and 930 +/- 300 microM for the fully ligated R state. The best fit equilibrium constant (125 +/- 40 microM) for three ligands bound does not differ significantly from that measured without ligands bound. Incidental to the measurement technique is the determination of the rates of binding and release of the dye. The association rate for binding to the T state is large, (at least 4 x 10(9) M-1 s-1) and may be diffusion limited, while the association and dissociation rates for R state binding, while not determined with precision, are clearly much smaller, of the scale of 10(5) M-1 s-1 for association.

Analysis of hemoglobin oxygen equilibrium curves. Are unique solutions possible?

As a hemoglobin tetramer is oxygenated, it converts from a form with low ligand affinity to a high-affinity structure. This allosteric transition occurs in partially liganded molecules, typically after two or three ligands are bound. As a result of the co-operative nature of the process, the populations of the partially liganded forms are low. The relative proportions and precise properties of these intermediate substrates are therefore difficult to measure and are subject to controversy. The problem is compounded by compensating effects; for example, over-estimation of the oxygen affinity of triply liganded forms will result in under-estimation of the proportion of these species. Specifically, published values for the oxygen affinity of the triply liganded species vary for identical conditions by a factor of more than 6. In analyses based on the highest affinity values, the triply liganded species virtually disappears. However, this affinity is usually calculated from the last few per cent of the oxygenation curve, and this part of the curve is extremely sensitive to the normalization of the data. We conclude that unique solutions for the ligand affinity and substrate populations may be impossible to achieve, and that unusual mechanisms based on particular combinations of parameters, therefore, should be viewed with caution.