©1996 by The American Society for Biochemistry and Molecular Biology, Inc.
Probing the Interface of Human Hemoglobin by Mutagenesis
ROLE OF THE FG-C CONTACT REGIONS (*)

(Received for publication, November 6, 1995; and in revised form, January 26, 1996)

Beatrice Vallone Andrea Bellelli Adriana E. Miele Maurizio Brunori (§) Giulio Fermi (1)

From the Department of Biochemical Sciences ``A. Rossi Fanelli'' and the Consiglio Nazionale delle Ricerche Centre of Molecular Biology, University of Rome ``La Sapienza,'' Piazza Aldo Moro, 5, 00185 Rome, Italy Laboratory of Molecular Biology, Medical Research Council, Hills Road, CB2 2QH Cambridge, United Kingdom

ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
FOOTNOTES
ACKNOWLEDGEMENTS
REFERENCES

ABSTRACT

The allosteric transition of hemoglobin involves an extensive reorganization of the alpha(1)beta(2) interface, in which two contact regions have been identified. This paper concerns the effect of two mutations located in the ``switch'' (alphaC3 Thr Trp) and the ``flexible joint'' (betaC3 Trp Thr). We have expressed and characterized one double and two single mutants: Hb alphaT38W/betaW37T, Hb betaW37T, and Hb alphaT38W, whose structure has been determined by crystallography.

We present data on: (i) the interface structure in the two contact regions, (ii) oxygen and CO binding kinetics and cooperativity, (iii) dissociation rates of deoxy tetramers and association rates of deoxy dimers, and (iv) the effect of NaI on deoxy tetramer dissociation rate constant.

All the mutants are tetrameric and T-state in the deoxygenated derivative. Reassociation of deoxygenated dimers is not modified by interface mutations. DeoxyHb alphaT38W dimerizes 30% slower than HbA; Hb betaW37T and Hb alphaT38W/betaW37T dissociate much faster. We propose a binding site for I at the switch region.

The single mutants bind O(2) cooperatively; the double one is almost non-cooperative, a feature confirmed by CO binding. The functional data, analyzed with the two-state model, indicate that these mutations reduce the value of the allosteric constant L(0).


INTRODUCTION

The three-dimensional structure of liganded and deoxyhemoglobin (Baldwin and Chothia, 1979; Shaanan, 1983; Perutz et al., 1987) indicates that the allosteric transition can be described topographically as a change in the relative orientation of the two dimers identified as alpha(1)beta(1) and alpha(2)beta(2), which rotate with respect to each other and slide along the alpha(1)beta(2) and alpha(2)beta(1) interfaces. Therefore, the amino acid residues that contribute to the alpha(1)beta(2) (and alpha(2)beta(1)) interface (Fig. 1) play a major role in controlling the relative stability of the allosteric states and, as a consequence, cooperativity. This interface contains residues of helix C and the FG corner of both chains. In total 17 residues establish interactions across the alphaFG-betaC and the alphaC-betaFG contacts with a pseudo-symmetric arrangement (Fig. 1). Other interactions of minor significance are present across the alphaFG-betaFG and alphaC-betaC contacts (Baldwin and Chothia, 1979).


Figure 1: Schematic structure of the alpha(1)beta(2) interface of HbA. Ribbon representation of alpha(1)beta(2) dimer. The regions involved in the contact are framed by the black square, and the residues represented in ball and stick mode are Thr alpha38 (C3) and Trp beta37 (C3).



The alphaFG-betaC contact is extensive and forms a network of weak bonds, which is largely maintained in the allosteric transition, in spite of some changes in the orientation of specific amino acid side chains. beta37 (C3) Trp in particular makes contacts with alpha92(FG4) Arg, alpha94(G1) Asp, and alpha95(G2) Pro in both oxy- and deoxyhemoglobin; nevertheless, it is a good probe of the allosteric transition because its optical spectrum is perturbed upon ligand binding (Briehl and Hobbs, 1970; Perutz et al., 1974). The amino acid residues at the alphaC-betaFG contact region of the alpha(1)beta(2) interface undergo an important reorganization in the course of the T-R transition. In deoxyhemoglobin the side chain of beta97(FG4) His lies between alpha44(CD2) Pro and alpha41(C6) Thr, while in oxyhemoglobin this same His settles between alpha41(C6) Thr and alpha38(C3) Thr. Thus, this region of the alpha(1)beta(2) interface, which seems to be compatible with only two states, plays a key role in determining the allosteric transition (see Fermi and Perutz(1981)).

Baldwin and Chothia(1979) summarize the function of the alpha(1)beta(2) interface with the definition of the alphaFG-betaC contact region as a ``flexible joint'' and the alphaC-betaFG as a ``switch.'' Hence the amino acid residue at position C3 has a different role in the alpha and beta chains; beta37(C3) Trp is involved in a series of contacts, which stabilize the tetramer in both the oxy and deoxy derivatives, while alpha38(C3) Thr participates in the reorganization of the alpha(1)beta(2) interface associated with the allosteric transition.

Human Hb has been expressed in Escherichia coli and yeast, and several interesting mutants have been prepared to test specific hypotheses (Nagai et al., 1985; Martin de Llano et al., 1993; Komiyama et al., 1995). In view of the considerations reported above, we have investigated the role of the residue at C3 in both chains of human hemoglobin by synthesizing the site-directed mutants bearing in the beta chains the residue found in the alpha chains and vice versa. The type of substitution to be inserted at C3 was chosen to probe the effect of changes in the number of atomic contacts (Schaad et al., 1993), avoiding the introduction of charged residues. Three hemoglobins were therefore expressed in E. coli: the single mutants alpha38 Thr Trp (alphaT38W) (^1)and beta37 Trp Thr (betaW37T) and the double mutant alpha38 Thr Trp/beta37 Trp Thr (alphaT38W/betaW37T, called herein the double mutant). It may be recalled that there are no known natural mutants of alpha38 Thr, while the two natural mutants of beta37 Trp, Hb Hirose and Hb Rothschild, have Ser and Arg, respectively (Yamaoka, 1971; Gacon et al., 1977; Sasaki et al., 1978) (see also Huisman(1992)).

Our three mutants were compared to wild type HbA with respect to the reaction with oxygen and carbon monoxide, the stability of liganded and unliganded tetramers, the kinetics of association of deoxy dimers and the dissociation of deoxy tetramers, the effect of sodium iodide, and the changes in the molecular contacts at the interface by x-ray crystallography of deoxy-Hb alpha38 Thr Trp or by modeling.

Structural analysis leads to simple predictions on the role of the alpha(1)beta(2) interface, suggesting that (a) semi-conservative mutations at topological position C3 affect tetramer stability and cooperativity in a simple manner, without loss of allosteric behavior; (b) substitutions at equivalent positions of the alpha and beta chains should yield different functional effects, in spite of the pseudo-symmetry of that interface. By-and-large these predictions are fulfilled by our results.


EXPERIMENTAL PROCEDURES

Materials

The site-directed mutant hemoglobins were produced in E. coli and purified as described previously (Nagai et al., 1985; Hoffman et al., 1990; Vallone et al., 1993). Reagents were of analytical grade.

Functional Properties

The oxygen binding isotherms were determined using the tonometric method (Rossi Fanelli and Antonini, 1958). The time course of CO binding to deoxyhemoglobin was followed using either the Applied Photophysics stopped flow apparatus (Leatherhead, United Kingdom) or the flash photolysis apparatus described by Brunori and Giacometti(1981).

The time course of dissociation of deoxyhemoglobin into dimers was monitored by recording the optical spectrum (in the Soret region) after mixing deoxyhemoglobin with a stoichiometric amount of haptoglobin (isoform 1.1 purchased from Sigma), as described by Ip et al.(1976). Since haptoglobin binds rapidly and tightly to the free dimers, the equilibrium is driven toward the dissociated state (Nagel and Gibson, 1972).

The time course of association of deoxyhemoglobin dimers was followed by mixing oxyhemoglobin (which at micromolar concentration contains a substantial amount of oxygenated dimers) with 50 mM sodium dithionite in the rapid scanning stopped flow spectrometer described by Bellelli et al.(1990). After rapid dissociation of oxygen, the dimers recombine in a slow, concentration-dependent process (Kellett and Gutfreund, 1970).

Transient spectra were collected into a matrix (A), each column being a difference spectrum and each row a time course, and analyzed with the singular value deconvolution algorithm (SVD; Golub and Reinsch, 1980; Henry and Hofrichter, 1992). This deconvolution yields the three matrices U, S, and V, whose product U times S times V^T approximates A. The product U times S may be envisaged as a matrix of extinction coefficients of the spectroscopic components detected by the algorithm, while V columns represent their time evolution. The SVD analysis allowed to magnify the small absorption changes coupled to the allosteric transition and to increase the signal to noise ratio.

Structural Properties

The deoxygenated derivative of the mutant alphaT38W was crystallized following Perutz(1968). X-ray data from a single crystal were collected and processed as described by Perutz et al.(1993). A total of 25,572 reflections were collected in the shell between 22.0 and 2.48 Å; these reduced with a merging R factor of 9.6% to 15,560 uniquely indexed reflections, or 76% of the unique reflections in the shell. Cell dimensions were the same as those of native deoxy-HbA, within experimental error. Mutant structure amplitudes were scaled to native with a scale factor and a temperature factor; the R factor between mutant and native was 13.5% for 15,442 reflections included in the map. A difference map was calculated from the differences of the scaled amplitudes with phases of the native model (Brookhaven Protein Data Bank code 2HHB; Fermi et al., 1984). The map was symmetry-averaged about the molecular dyad to improve accuracy.

Molecular modeling was carried out on a Silicon Graphics workstation (Silicon Graphics Inc., Mountain View, CA) by use of the Discover/Insight package. Two structures were from the Brookhaven Protein Data Bank: deoxygenated recombinant hemoglobin (Kavanaugh et al., 1992a) and wild type oxyhemoglobin (Shaanan, 1983). For mutant alphaT38W we used the coordinates of the deoxy derivative. For the other mutants (i.e. the two proteins containing the substitution betaW37T and the oxy form of Hb alphaW38T) simulation of the mutation in the corresponding wild type structure was performed by substituting the side chain in the starting position leading to no unfavorable contacts with neighboring side chains, followed by energy minimization allowing only amino acids in a sphere of 6 Å from the mutated amino acid to move.


RESULTS

Oxygen Equilibrium Experiments

The oxygen binding properties of the three site-directed mutants and HbA, depicted in Fig. 2, show that both single mutants are cooperative and display Hill coefficients larger than 2, while the double mutant is almost non-cooperative (n = 1.3).


Figure 2: Oxygen binding isotherms of HbA and mutant hemoglobins. Open squares, HbA; closed squares, Hb alphaT38W; open circles, Hb betaW37T; closed circles, Hb alphaT38W/betaW37T. Continuous lines represent the best fit to a two-state model with the parameters reported in Table 1. Conditions: 0.1 M Bis-Tris/HCl buffer, pH 7.0, T = 20 °C. The magnitude of the symbols corresponds to approximately ± 1 standard deviation (± 0.017 on the ordinate scale).





In spite of its limitations (see, for example, Ackers et al. (1992)), the analysis of the oxygen binding isotherms was carried out according to the two-state allosteric model (Monod et al., 1965). Fit of data is often beset by uncertainties in K(R) and L(0). Nonetheless, since kinetic data on O(2) and CO have provided independent support that K(R) is identical or very close to that of the isolated alpha and beta chains and of the alphabeta dimer (Edelstein and Gibson, 1987; Szabo and Karplus, 1972), we fitted the data in Fig. 2assuming K(R) to be the same for HbA and all the mutants, and thus similar to the oxygen affinity of the isolated chains. This assumption proved compatible with a good fit, with the parameters presented in Table 1. By contrast, very poor fits were obtained if K(T) of the mutants was imposed to be the same as that of HbA, which is not surprising and agrees with the theory of Szabo and Karplus(1972). Within the limitations of the two-state model, the numerical values of the allosteric parameters indicate that mutations at C3 (i) affect the equilibrium constant between the two quaternary states, reducing the difference in stability between T(0) and R(0) in unliganded Hb (see values of L(0)), and (ii) reduce to some extent the constraints imposed on the subunits by the quaternary assembly as deduced from the values of K(T).

Kinetics of Combination with CO

As an additional test of the conclusions drawn from O(2) binding data, we have determined the CO combination rate constants by stopped-flow and compared them to the values obtained by partial flash photolysis. The time course of CO binding to HbA and Hb alphaT38W by stopped flow is autocatalytic (Fig. 3), which confirms cooperative ligand binding; for Hb betaW37T and Hb alphaT38W/betaW37T, the same reaction appears non-autocatalytic, although the apparent rate constant is still much lower than that of the R-state.


Figure 3: Time course of CO combination by stopped flow. Open squares, HbA; closed squares, Hb alphaT38W; open circles, Hb betaW37T; closed circles, Hb alphaT38W/betaW37T. Lines are drawn according to a sequential scheme (see text). CO was 50 µM after mixing; other experimental conditions were as described for Fig. 2.



The data reported in Fig. 3were fitted to four consecutive and irreversible pseudo-first order processes, as described by Hopfield et al.(1971). The values of L(0) and c (K(R)/K(T)) used to calculate the relative amount of T- and R-state Hb at each ligation state were taken from Table 1; the values of ^Rk` were obtained by partial flash photolysis, which is known to populate partially liganded species that recombine as fast as the isolated chains (Sawicki and Gibson, 1976; Vandegriff et al., 1991; Jones et al., 1992). The CO combination time courses in Fig. 3are satisfactorily fitted with this model, regardless of the presence of evident kinetic cooperativity, and the rate constants are given in Table 2. Regarding Hb betaW37T, an autocatalytic time course may have been expected. However, it is known that speeding up of the CO combination is very sensitive to the value of c, which may not be identical for O(2) and CO, even though this difference is sometimes very difficult to demonstrate. Under these conditions the estimate of ^Tk` is subject to some uncertainties, and the overall time course is essentially exponential.



Even complete photolysis of dilute solutions of HbCO usually follows a biphasic rebinding time course (Antonini and Brunori, 1971). This has been attributed to the presence of rapidly reacting dimers in equilibrium with the slowly reacting tetramers (Gibson and Antonini, 1967; Edelstein et al., 1970). Thus flash photolysis has also been used to probe the extent of dissociation of HbCO into dimers. Employing this technique, Vallone et al. (1993) have shown that Hb alphaT38W CO is a more stable tetramer than HbA CO, by approximately 0.6 kcal/mol. Similar experiments carried out as a function of Hb concentration (data not shown) have demonstrated that the CO derivative of the two betaW37T mutants is almost completely dissociated into dimers at [Hb] = 15 µM, and thus a lower limit of 40 µM can be set to the value of the tetramer-dimer dissociation constant of these mutants.

Since O(2) and CO exhibit nearly parallel equilibrium curves, the data in Table 1and Table 2may be considered together even though most of the cooperativity is expressed in the dissociation rate constants in the case of O(2) and in the association rate constants in the case of CO (Antonini and Brunori, 1971; Szabo, 1978).

A significant conclusion is that all our mutants display a ligand affinity higher than HbA in the T-state. However, this increase is smaller for the two single mutants, and substantial for the double mutant, in keeping with the reduced number of contacts at the alpha(1)beta(2) interface (see below). More interestingly, the value of L(0) decreases significantly in the mutants, being, in the double mutant, 15-fold smaller than in HbA.

Recombination of Dimers following Oxygen Dissociation

The time course of the recombination of deoxygenated alphabeta dimers from HbA and the mutant Hbs was studied by recording the slow optical transitions that follow the rapid deoxygenation by dithionite of a dilute solution of HbO(2) (Antonini et al., 1968; Kellett and Gutfreund, 1970), given that, at micromolar concentration, HbO(2) is an equilibrium mixture of tetramers and dimers. This is possible because (a) the optical spectrum of deoxyhemoglobin in the T-state differs from that in the R-state (the transient R(0) Hb, as described by Sawicki and Gibson(1976); the deoxy dimers and the deoxy chains, as reviewed by Bellelli and Brunori(1994)), and (b) the T(0)-R(0) difference spectrum is, within errors, the same for HbA and our three mutants.

The experiment is described in the following scheme ().(^2)

An example of this type of experiment for HbA and the three mutants is reported in Fig. 4. The plot shows the second column of V (from the SVD analysis), which contains most of the slow optical transition. The faster phase seen in this figure is the time course of O(2) dissociation (upon mixing with dithionite), and is described by steps 1 and 3 of . The slower phase (step 2 in the same scheme) relates to the time course of association of deoxy dimers. The relative amplitude of this phase is a measure of the fraction of oxy dimers in equilibrium with oxy tetramers. Since addition of 0.5 M NaI promotes complete dissociation of the liganded tetramers but does not prevent the reassociation of unliganded dimers (Kellett and Gutfreund, 1970), experiments carried out under these experimental conditions yield the full amplitude of the slow phase. Estimates for the rate constant of deoxy dimers association for HbA and the three mutants (in the absence and in the presence of 0.5 M NaI) are reported in Table 3.


Figure 4: Time course of the oxygen dissociation (rapid upward phase) and dimers reassociation (slower downward phase) for HbA (A), Hb alphaT38W (B), Hb betaW37T (C), and Hb alphaT38W/betaW37T (D). Each panel depicts the time evolution of the amplitude of column 2 of the V matrix obtained form the SVD of an experiment carried out at two hemoglobin concentrations, namely 2.5 (open symbols) and 1.25 (closed symbols) µM after mixing (heme basis). Conditions were as described for Fig. 2.





Rate of Dissociation of Unliganded Tetramers into Dimers

The dissociation of deoxyhemoglobin was followed by the absorbance changes at 430 nm after reaction with haptoglobin, recorded either statically or kinetically (Ip et al., 1976). This method relies on the essentially irreversible binding of Hp to free alphabeta dimers, which forces the tetramer-dimer equilibrium toward dissociation, the rate-limiting step being dimerization (see also Nagel and Gibson (1972)). In view of the results reported above, the effect of NaI has also been tested on unliganded Hb.

The static difference spectra obtained upon binding to Hp are very similar to those obtained by the kinetics of dimer recombination (data not shown). The total absorbance change for complete dissociation of deoxy tetramers into dimers was found to be 13% of the total absorbance of the sample at 430 nm. This is in good agreement with the value expected for the T(0)-R(0) spectral change (for a review, see Bellelli and Brunori(1994)).

Fig. 5shows the observed time courses of dissociation. This experimental approach suffers from the extremely slow time course, which, at lower temperatures, extends over several days; thus, for some proteins and under some conditions, the optical transition could not be followed to completion. Nonetheless, the time course was fitted in all cases to a single exponential following Ip et al.(1976). In addition, the asymptotic value was independently estimated from the total optical density change (see above).


Figure 5: Time course of dissociation of the deoxygenated tetramer of HbA and mutants induced by binding of haptoglobin, in the presence and absence of NaI. The percentage of the reaction, followed by the decrease of absorbance at 430 nm, is plotted versus time, and fitted using a monoexponential equation. Panel A, HbA; panel B, Hb alphaT38W; panel C, Hb betaW37T; panel D, Hb alphaT38W/betaW37T. Open symbols, buffer 0.1 M Bis-Tris pH 7.0; closed symbols, the same buffer as before plus 0.5 M NaI, temperature = 20 °C.



Table 4reports the first order rate constants determined from the data in Fig. 5. HbA and Hb alphaT38W behave similarly, while both Hb betaW37T and the double mutant dissociate faster than HbA (in spite of being tetrameric in the deoxygenated derivative, as judged by the recovery of the expected optical transition). The effect of NaI on the dissociation of HbA and Hb betaW37T is the same (see the ratio k`/k in Table 4). However, quite unexpectedly, this effect is considerably reduced in the two mutants in which alpha38 (C3) is mutated to Trp; this point will be discussed below, with reference to the structure of the interface.



Crystallographic Structure of Deoxy-Hb alpha38 Thr Trp

Fig. 6(panel A) shows the difference map F(Hb alphaT38W)-F(HbA) in the deoxy state; in panel B only the negative values of the difference map (red contours) are shown, and the HbA structure is displayed. The excellent fit of the Trp side chain to the main positive peak may be seen. The largest negative peak represents the absence in this mutant of Thr alpha38 O-1 and a bound water molecule, seen in the structure of HbA; a second less intense negative peak may correspond to another water molecule. In the mutant, Trp C overlaps with the native Thr C-2 of HbA and thus no difference peak is observed. It may be clearly seen that there is no evidence of any change in position of neighboring residues of the alpha(1)beta(2) interface.


Figure 6: Symmetry averaged difference map F(Hb alphaT38W) - F(HbA). The map is contoured at +0.15 (green contours) and -0.15 (red contours), or approximately 3 times the root-mean-square density value of the unaveraged map. Panel a, the superimposed model is that of wild type HbA with alpha38 Thr replaced by Trp with 1 angle -60°. Residues of the alpha(1) chain are labeled with their ordinary position numbers and the letter A, those of the beta(2) chain with position numbers increased by 600 and the letter D. Panel b, only the negative values of the difference map are shown (red contours) and the HbA structure is displayed. The dashed lines represent hydrogen bonds. Labels as in panel a.



Molecular Modeling

Modeling of the T and R structures of our mutants (over and above the crystallographic structure of deoxy-Hb alphaT38W) was limited to simple geometric computations, such as the number of atoms in contact with the side chain at position C3, and the number of side chains interactions lost or established after amino acid substitutions. Nevertheless even this simple analysis yielded information useful for the interpretation of functional experiments.

As shown in Table 5, upon introduction of a Trp at alpha38 (C3) the increase in the number of interface contacts between alpha(1) and beta(2) is negligible in deoxy, but considerable in oxy-Hb. Conversely the substitution Trp Thr at beta37 (C3) leads to a substantial decrease in the number of contacts in both oxy- and deoxy-Hb, with loss of a hydrogen bond between beta37 Trp and alpha94 Asp in the deoxy-Hb (see Fig. 7, panel C). Hence one could deduce that compared to HbA, the alphaT38W mutant has a more compact and extensive alpha(1)beta(2) interface especially in the R-state; on the other hand, Hb betaW37T has a looser interface, especially in the T-state, with loss of a hydrogen bond.




Figure 7: Residues at position C3 (shown in black) and their neighboring amino acids in a sphere of 4 Å for oxy- and deoxy-HbA (Table 5). Residues belonging to the alpha and beta chains are shown in light gray. The labels are indicated only for the monomer opposite the C3 residue. The dots represent the water accessible surface (probe radius 1.4 Å) within the HbA tetramer. Hydrogen bonds involving the C3 side chain (either with other amino acid residues or with water molecules) are represented by a gray dotted line. Oxy-Hb: a, alpha38Thr from the crystallographic structure (Shaanan, 1983); b, alpha38Trp modeled starting from a. Deoxy-Hb: c, beta37Trp from the crystallographic structure (Kavanaugh et al., 1992b); d, beta37Thr modeled from c; e, alpha38Thr from the crystallographic structure (Kavanaugh et al., 1992b); f, alpha38Trp form the crystallographic structure (this paper). It has to be noticed that the number of residues in a sphere of 4 Å from the amino acid side chains of interest are fewer when the side chain is smaller and larger when this is greater; therefore, the amino acids displayed in panels a and b, c and d, and e and f of this figure are not the same.



Assuming that these effects are additive in the double mutant, one may expect that in this hemoglobin cooperativity will be severely impaired, since the R-state would be favored with a concomitant destabilization of the T-state; this prediction agrees well with the experimental data reported above.

In order to compare the data obtained on Hb betaW37T with observations reported in the literature on Hb Hirose (which bears the substitution betaW37S, Sasaki et al., 1978), we have measured the atoms in a sphere of 4 Å to a Ser in betaC3. As may be seen in Table 5, the number of contacts in Hb Hirose is reduced, in both oxy and deoxy state, compared to our mutant betaW37T, in agreement with the experimental data (see ``Discussion'').

Inspection of the water accessible surface in the proximity of position C3 allows evaluation of the cavities at the edge of the alpha(1)beta(2) interface, and how the shape of the surface is affected by the mutations. The crystallographic structure of deoxy-Hb alphaT38W shows that the Trp in alphaC3 lies flat on the edge of the interface, excluding beta145 Tyr and beta100 Pro from contact with the solvent, and thus seems to act as a ``hydrophobic plug'' (Fig. 7, panel F). Modeling the oxygenated state of the same mutant shows that the side chain of the Trp introduced in alphaC3 could also fit nicely in a cavity, in contact with beta145 Tyr.

The substitution betaW37T, on the other hand, seems to deepen a cavity at the edge of the opposite side of the alpha(1)beta(2) interface, possibly allowing to alpha141 Arg greater accessibility to the external medium.


DISCUSSION

Allosteric Properties of the Mutants

This paper presents a correlation between structural and functional properties of two single mutants and one double mutant in the alpha(1)beta(2) interface of human hemoglobin. The mutations are in the so called ``flexible joint'' and ``switch'' regions, which are pseudo-symmetric, alphaT38W being in the switch region, and betaW37T in the flexible joint; the double mutant contains both mutations. The first mutation (Thr Trp) might be expected to affect primarily the allosteric properties, while the second (Trp Thr) might also affect the stability of liganded and unliganded tetramers. Most, but not all, of the results presented above bear out these predictions. The properties of these new mutants are discussed in the framework of the two-state MWC model (Monod et al., 1965), in spite of the limitations that have been discussed in the literature (see, for example, Ackers et al.(1992)). The model provides a convenient analytical description of the data and allows discussion of the more interesting findings with a nomenclature known to most.

Hb alphaT38W, which bears two Trp residues in the alpha(1)beta(2) interface, is a slightly more stable tetramer than HbA, especially in the liganded form where the tetramer dissociation constant is decreased 6-fold (Vallone et al., 1993). The three-dimensional structure of this mutant in the deoxy state shows that the indole side chain of Trp alpha38 has been accommodated by expulsion of two water molecules hydrogen-bonded to alpha38 Thr in HbA, without significant changes of the adjacent residues (Fig. 6). Thus, the enhanced stability of the oxy tetramer is consistent with the increased hydrophobic character of the interface (which is not perturbed by the bulkier side chain), with the more extensive contacts of the indole (Table 5), and with the hypothesis that the most significant contribution to the stability of the tetramer resides in the pseudo symmetric alphaFG-betaC contact, which is unmodified. The slightly increased affinity of the T-state observed for Hb alphaT38W as compared to HbA is not easily accounted for, given that the constraints in the T-state are not decreased. On the other hand, the lower value of L(0) is consistent with the increased stability of R and thus with a (slightly) reduced cooperativity.

Hb betaW37T, which has no Trp residues in the alpha(1)beta(2) interface, dissociates into dimers more readily than HbA, in both the unliganded and liganded derivatives, as shown by results from analytical ultracentrifuge and flash photolysis to be published elsewhere, indicating that Trp at betaC3 contributes to the stability of the alpha(1)beta(2) interface in both quaternary structures. In spite of the destabilization of this interface, cooperativity of betaW37T is reduced but preserved, not only because the contacts in the alphaC-betaFG switch region (beta97 His, alpha41 Thr, alpha44 Pro, and alpha38 Thr) are preserved, but also because Thr at beta37(C3) can to some extent fulfill the role of Trp in the flexible joint, as indicated by modeling (Table 5). The remarkable but limited competence of Thr at beta37 residue is reflected in a small increase of the O(2) affinity of the T-state (Table 1) and a somewhat greater increase of the rate constant for CO binding to deoxy-Hb (Table 2); thus, the constraints of the T-state are partially released by mutation of betaC3 and essentially maintained by mutation of alphaC3.

The interesting behavior of the double mutant Hb alphaT38W/betaW37T may be understood qualitatively on the basis of the properties of the two single mutants. Although the double mutant, like HbA, has only one Trp residue at the alpha(1)beta(2) interface, ligand binding and tetramer dissociation are both remarkably different from wild type. This is the best evidence that the effect of mutations at topological position C3 is asymmetric and that perturbing the flexible joint and the switch indeed has different consequences. This double mutant dissociates more readily into dimers in both the liganded and unliganded derivatives, which we attribute largely to the flexible joint mutation betaW37T. However, cooperativity is also reduced, even when compared to that of the two single mutants, as shown by the increased O(2) affinity of the T-state (Table 1). These properties may be understood on the basis of two synergistic effects, i.e. a significant release of the interface constraints due to betaW37T, leading to a more relaxed T-state (hence the smaller K(T)), and an increase in the stability of the R-state tetramer, related to mutation alphaT38W. As a result the energy difference between the two allosteric states is reduced, as indicated by the value of the allosteric constant L(0), which is smaller than HbA by 15-fold.

As already stated, there is no known natural mutant of position alpha38 (Huisman, 1992); however, two site-directed mutants of this residue, namely Hb alphaT38S and alphaT38V, have been expressed and characterized by Hashimoto et al.(1993). The functional properties of both these Hbs are by and large similar to those of HbA. However, the substitutions in these case are semiconservative, and it is easily conceivable that the residues effectively replace the wild type Thr.

As to position beta37, our data may be compared with those of the natural mutants Hb Hirose and Hb Rothschild. Comparison of Hb betaW37T with Hb Rothschild (beta37 Trp Arg; Gacon et al.(1977)) is complex, because of arginine's positive charge. Nonetheless, consistent with the other beta37 mutants, Hb Rothschild displays reduced cooperativity and an increased tendency to dissociate into dimers. However, the structure of the deoxygenated derivative (Kavanaugh et al., 1992b) shows a novel and strong chloride binding site, which affects the functional properties of this Hb.

Hb Hirose (beta37 Trp Ser) displays extremely low (if any) cooperativity, very rapid CO binding by flow, and a high tendency to dissociate into dimers in both the presence and the absence of oxygen (Sasaki et al., 1978). Hb betaW37T, on the other hand, although extensively dissociated into dimers when liganded, is fully associated as the deoxy derivative (Fig. 5) and maintains a higher cooperativity than Hb Hirose, confirming that Thr but not Ser is a partially competent substitute for beta37 Trp at the flexible joint. The computer-simulated structures of oxy- and deoxy-Hb betaW37T (Fig. 7) indicate that the methyl group of Thr beta37 makes contacts with alpha95 Pro and alpha140 Tyr in the oxy derivative and with alpha140 Tyr in the deoxy derivative; these contacts are lost or less extensive with Ser (Table 5).

A site-directed mutant of position beta37, in which a Phe substitutes the Trp, has been obtained by Ishimori et al.(1992). In this case the Phe would be expected to provide a more effective replacement for Trp than either Hb Hirose or our mutant. Unfortunately the authors did not measure the tetramer-dimer equilibrium and kinetics on their mutant Hb. As with betaW37T, Hb betaW37F displays high oxygen affinity and low cooperativity, even though the ^1H NMR spectrum of this mutant is consistent with that typical of T-state HbA.

Rates of Dimer Reassociation and Tetramer Dissociation in Deoxyhemoglobin

Fitting the time course of association of deoxy alphabeta dimers to a second order reaction (as reported in Fig. 4) shows small but systematic deviations, which are more marked for the mutants than for HbA. Nevertheless, the time course of reassociation of deoxy alphabeta dimers of HbA, when fitted to the simplest possible scheme (see , under ``Results''), yields a rate constant at low salt concentration slightly larger than, but not inconsistent with, that reported in the literature (Gray, 1974; Wiedermann and Olson, 1975). Small differences may be due either to the experimental conditions or the analysis. Among the rate constants reported in Table 3, those determined for Hb betaW37T and for the double mutant are more reliable because these Hbs are completely dissociated into dimers as HbO(2) and thus the total signal is larger and the initial dimer concentration does not need to be fitted.

Interestingly, Table 3shows that the rate constant of reassociation of unliganded dimers and the effect of NaI are the same (within a factor of two) for all four hemoglobins. Thus, the effect of interface mutation(s) on the stability of the deoxygenated Hb tetramer is only very slightly (if at all) reflected in its rate of association, in agreement with Ackers et al.(1992); the same conclusion holds for the dissociating effect of NaI. These observations imply that formation of the productive dimer-dimer complex is marginally affected by one or two conservative changes at the contact interface. This is not wholly surprising in view of the large surface buried in the contact between two alphabeta dimers (1500 Å^2).

On the other hand, the effect of mutations on the rate of dissociation of the deoxy tetramer is more marked and follows an interpretable trend (Table 4). In deoxy-Hb alphaT38W there are 8 intersubunit contacts compared to 7 in deoxy-HbA, consistent with just a 30% decrease in the rate of tetramer dissociation (Table 4). The beta37 mutants, on the other hand, both dissociate faster than Hb A, which we attribute to the loss (i) of the contact between alpha141 Arg and beta37 Trp with perturbation of the network of salt bridges present in the T-state, and (ii) of a hydrogen bond between beta37 Trp and alpha94 Asp (Fermi and Perutz, 1981). Therefore, the reduced number of contacts (Table 5) agrees well with the increased rate of dissociation of the beta37 mutants (which can be as much as 25-fold).

However, it should be pointed out that the rate constants for dissociation of deoxy-Hb betaW37T and Hb alphaT38W/betaW37T are quite different. This implies that the effect of the two mutations is not additive; otherwise, the double mutant would dissociate as fast as Hb betaW37T. This deviation from simple additivity suggest that the interface may be also distorted by long range effects. This is not unreasonable, given that perturbation of a single residue even at the alpha(1)beta(1) interface leads to a large destabilization of the T-state and loss of cooperativity (see, for example, Amiconi et al.(1989), Weber et al.(1993), and Zhang et al.(1996)). A thorough interpretation of non additive effects necessarily implies taking into account long range interactions, which indeed are well documented. However, we would rather to defer this analysis until the crystallographic structure of these mutants (now in progress) is available.

Effect of Sodium Iodide

The iodide ion strongly promotes dissociation of the hemoglobin tetramer (Kellett and Gutfreund, 1970), more than other anions such as Cl (Antonini and Brunori, 1971). We found that I affects both the combination and dissociation rate constants ( Table 3and Table 4). The 25-50-fold slower reassociation of deoxy dimers may be due to binding of I to several positively charged residues of the alpha(1)beta(2) interface (alpha92 Arg, beta40 Arg, and alpha40 His). Nonetheless, the effect of I on the bimolecular association rate constant of the mutants is approximately the same as for HbA (Table 3).

On the other hand, the effect of I on the rate of dissociation of the deoxy tetramer depends in a rational way on the amino acid residue at position C3. This is evident comparing the ratio of the rate constants k` and k reported in Table 4. With a Trp in the switch region as in Hb alphaT38W, I accelerates dissociation only 3-fold compared to 13-fold in HbA. We interpret the reduced effect as resulting from the more hydrophobic character of the interface and from the position of Trp, which seems to act as a hydrophobic plug, shielding from contact with the solvent beta145 Tyr and beta100 Pro and forming a barrier to the penetration of this ion.

When Trp is removed from the flexible joint as in Hb betaW37T, I accelerates dissociation approximately 13-fold, just as in HbA. In deoxy-Hb, beta37 Trp is on the edge of the alpha(1)beta(2) interface ( Fig. 1and Fig. 7), and its orientation is such that alpha95 Pro and alpha140 Tyr are also exposed to the external medium; mutation of Trp with Thr therefore does not result in exposure to the solvent of side chains previously buried in the interface. Thus, the similarity in the effect of I on the dissociation rate constants of HbA and Hb betaW37T mutant is fully consistent with the proposal that access of the anion to the alpha(1)beta(2) interface is toward the switch region, near position alpha38, and therefore not affected by substitution in the flexible joint region. In support of this hypothesis, the effect of I on the double mutant is again smaller and similar to that of Hb alphaT38W, because of the hydrophobic plug effect. Thus, we conclude that the effect of I on the stability of the deoxy tetramer is linked to binding of this anion to the interface near position C3. This prediction may be checked by examination of the crystallographic structure of deoxy-HbA in the presence of NaI.


FOOTNOTES

*
This work was supported by a grant from the Ministero dell'Universitá e della Ricerca Scientifica e Tecnologica (40%, Liveprotein). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore by hereby marked ``advertisement'' in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

The atomic coordinates (1GLI) and structure factors (R1GLISF) have been deposited in the Protein Data Bank, Brookhaven National Laboratory, Upton, NY.

§
To whom correspondence should be addressed.

(^1)
The abbreviations used are: alphaT38W, mutant Hb alpha38 Thr Trp; betaW37T, mutant Hb beta37 Trp Thr; alphaT38W/betaW37T, mutant Hb alpha38 Thr Trp/beta37 Trp Thr; SVD, singular value deconvolution algorithm.

(^2)
Although has been widely applied to the description of this reaction (Kellett and Gutfreund, 1970; Wiedermann and Olson, 1975), the data obtained by recording the complete optical transition over the range 380 to 520 nm, and analyzed by SVD (see ``Experimental Procedures'') suggest that the kinetics may be more complex than previously thought. This observation is being further explored and will be presented elsewhere, since is sufficient for description of the experiments presented here.


ACKNOWLEDGEMENTS

We express our sincere thanks to Dr. K. Nagai (Medical Research Council (MRC), Cambridge, United Kingdom) for assistance in the preparation of site-directed mutants and to Dr. M. F. Perutz (MRC) for carefully reading a preliminary version of the manuscript. We also thank B. Volpi of the Istituto Superiore di Sanitá (Rome) for invaluable assistance in bacterial growth. We gratefully acknowledge Somatogen Inc. for permission to use the Hb expression system in E. coli for scientific purposes.


REFERENCES

  1. Ackers, G. K., Doyle, M. S., Myers, D., and Daugherty, M. A. (1992) Science 255, 54-63 [Medline] [Order article via Infotrieve]
  2. Amiconi, G., Ascoli, F., Barra, D., Bertollini, A., Matarese, R. M., Verzili, D., and Brunori, M. (1989) J. Biol. Chem. 264, 17745-17749 [Abstract/Free Full Text]
  3. Antonini, E., and Brunori, M. (1971) Hemoglobin and Myoglobin in Their Reactions with Ligands , North Holland, Amsterdam
  4. Antonini, E., Brunori, M., and Anderson, S. (1968) J. Biol. Chem. 243, 1816-1822 [Abstract/Free Full Text]
  5. Baldwin, J., and Chothia, C. (1979) J. Mol. Biol. 129, 175-220 [Medline] [Order article via Infotrieve]
  6. Bellelli, A., and Brunori, M. (1994) Methods Enzymol. 232, 56-71 [Medline] [Order article via Infotrieve]
  7. Bellelli, A., Antonini, G., Brunori, M., Springer, B. A., and Sligar, S. G. (1990) J. Biol. Chem. 265, 18898-18901 [Abstract/Free Full Text]
  8. Briehl, R. W., and Hobbs, J. F. (1970) J. Biol. Chem. 245, 544-554 [Abstract/Free Full Text]
  9. Brunori, M., and Giacometti, G. M. (1981) Methods Enzymol. 76, 582-595 [Medline] [Order article via Infotrieve]
  10. Edelstein, S., Rehmar, M. J., Olson, J. S., and Gibson, Q. H. (1970) J. Biol. Chem. 245, 4372-4381 [Abstract/Free Full Text]
  11. Fermi, G., and Perutz, M. F. (1981) Hemoglobin and Myoglobin: Atlas of Protein Sequence and Structure , Vol. 2, Oxford University Press, New York
  12. Fermi, G., Perutz, M. F., Shaanan, B., and Fourme, R. J. (1984) J. Mol. Biol. 175, 159-174 [Medline] [Order article via Infotrieve]
  13. Gacon, G., Belkhodja, O., Wajcoman, H., Labie, D., and Najman, A. (1977) FEBS Lett. 82, 243-246 [CrossRef][Medline] [Order article via Infotrieve]
  14. Gibson, Q. H., and Antonini, E. (1967) J. Biol. Chem. 242, 4678-4681 [Abstract/Free Full Text]
  15. Gibson, Q. H., and Edelstein, S. (1987) J. Biol. Chem. 262, 516-519 [Abstract/Free Full Text]
  16. Golub, G. H., and Reinsch, C. (1970) Numer. Methods 14, 403-420
  17. Gray, R. (1975) J. Biol. Chem. 249, 2879-2885 [Abstract/Free Full Text]
  18. Hashimoto, M., Ishimori, K., Imai, K., Miyazaki, G., Morimoto, H., Wada, Y., and Morishima, I. (1993) Biochemistry 32, 13688-13695 [Medline] [Order article via Infotrieve]
  19. Henry, E. H., and Hofrichter, J. (1992) Methods Enzymol. 210, 129-192
  20. Hoffman, S. J., Looker, D. L., Roerich, J. M., Cozart, P. E., Durfee, S. L., Tedesco, J. L., and Stetler, G. L. (1990) Proc. Natl. Acad. Sci. U. S. A. 87, 8521-8525 [Abstract]
  21. Hopfield, J. J., Schulman, R. G., and Ogawa, S. (1971) J. Mol. Biol. 61, 424-433
  22. Huisman, T. H. J. (ed) (1992) International Information Center 16, 127-236
  23. Ip, S. H. C., Johnson, M. L., and Ackers, G. K. (1976) Biochemistry 15, 654-660 [Medline] [Order article via Infotrieve]
  24. Ishimori, K., Imai, K., Miyazaki, G., Kitagawa, T., Wada, Y., Morimoto, H., and Morishima, I. (1992) Biochemistry 31, 3256-3264 [Medline] [Order article via Infotrieve]
  25. Jones, C. M., Ansari, A., Henry, E. R., Garrott, W. C., Hofrichter, J., and Eaton, W. A. (1992) Biochemistry 31, 6692-6702 [Medline] [Order article via Infotrieve]
  26. Kavanaugh, J. S., Rogers, P. H., and Arnone, A. (1992a) Biochemistry 31, 8640-8649 [Medline] [Order article via Infotrieve]
  27. Kavanaugh, J. S., Rogers, P. H., Case, D. A., and Arnone, A. (1992b) Biochemistry 31, 4111-4121 [Medline] [Order article via Infotrieve]
  28. Kellettt, G. L., and Gutfreund, H. (1970) Nature 227, 921-926 [Medline] [Order article via Infotrieve]
  29. Komiyama, N. H., Miyazaki, G., Tame, J., and Nagai, K. (1995) Nature 373, 244-246 [CrossRef][Medline] [Order article via Infotrieve]
  30. Martin de Llano, J. J, Jones, W., Schneider, K., Chait, B. T., Manning, J. M., Rodgers, G., Benjamin, L. J., and Weksler, B. (1993) J. Biol. Chem. 268, 27004-27011 [Abstract/Free Full Text]
  31. McGovern, P., Reisberg, P., and Olson, J. S. (1976) J. Biol. Chem. 251, 7871-7879 [Abstract]
  32. Monod, J., Wyman, J., and Changeux, J. P. (1965) J. Mol. Biol. 12, 88-118 [Medline] [Order article via Infotrieve]
  33. Nagai, K., Perutz, M. F., and Poyart, C. (1985) Proc. Natl. Acad. Sci. U. S. A. 82, 7252-7255 [Abstract]
  34. Nagel, R.l., and Gibson, Q. H. (1972) Biochem. Biophys. Res. Commun. 48, 959-966 [Medline] [Order article via Infotrieve]
  35. Perutz, M. F. (1968) J. Cryst. Growth 2, 54-56 [CrossRef]
  36. Perutz, M. F, Ladner, J. F., Simon, S. R., and Ho, C. (1974) Biochemistry 13, 2163-2169 [Medline] [Order article via Infotrieve]
  37. Perutz, M. F., Fermi, G., Luisi, B., Shaanan, B., and Liddington, R. C. (1987) Acc. Chem. Res. 20, 309-321
  38. Perutz, M. F., Fermi, G., Poyart, C., Pagnier, J., and Kister, J. (1993) J. Mol. Biol. 233, 536-545 [CrossRef][Medline] [Order article via Infotrieve]
  39. Rossi Fanelli, A., and Antonini, E. (1958) Arch. Biochem. Biophys. 77, 478-482
  40. Sasaki, J., Imamura, T., Yanase, T., Atha, D. H., Riggs, A., Bonaventura, J., and Bonaventura, C. (1978) J. Biol. Chem. 253, 87-94 [Medline] [Order article via Infotrieve]
  41. Sawicki, C. A., and Gibson Q. H., (1976) J. Biol. Chem. 251, 1533-1542 [Abstract]
  42. Schaad, O., Vallone, B., and Edelstein, S. J. (1993) C. R. Acad. Sci. Paris 316, 564-571 [Medline] [Order article via Infotrieve]
  43. Shaanan, B. (1983) J. Mol. Biol. 171, 31-59 [Medline] [Order article via Infotrieve]
  44. Szabo, A. (1978) Proc. Natl. Acad. Sci. U. S. A. 75, 2108-2111 [Abstract]
  45. Szabo, A., and Karplus, M. (1972) J. Mol. Biol. 251, 1533-1542
  46. Vallone, B., Vecchini, P., Cavalli, V., and Brunori, M. (1993) FEBS Lett. 324, 117-122 [CrossRef][Medline] [Order article via Infotrieve]
  47. Vandegriff, K. D., Le Tellier, Y. C., Winslow, R. M., Rohlf, R. S., and Olson, J. S. (1991) J. Biol. Chem. 266, 17049-17059 [Abstract/Free Full Text]
  48. Weber, R. E., Jessen, T. H., Malte, H., and Tame, J. (1993) J. Appl. Physiol. 75, 2646-2655 [Abstract]
  49. Wiedermann, B. L., and Olson, J. S. (1975) J. Biol. Chem. 250, 5273-5275 [Abstract]
  50. Yamaoka, K. (1971) Blood 30, 730-738
  51. Zhang, J., Hua, Z., Tame, J. R. H., Lu, G., Zhang, R., and Gu, X. (1996) J. Mol. Biol. , in press

©1996 by The American Society for Biochemistry and Molecular Biology, Inc.




This Article
Abstract
Full Text (PDF)
Purchase Article
View Shopping Cart
Alert me when this article is cited
Alert me if a correction is posted
Citation Map
Services
Email this article to a friend
Similar articles in this journal
Similar articles in PubMed
Alert me to new issues of the journal
Download to citation manager
Copyright Permissions
Google Scholar
Articles by Vallone, B.
Articles by Fermi, G.
Articles citing this Article
PubMed
PubMed Citation
Articles by Vallone, B.
Articles by Fermi, G.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 All ASBMB Journals   Molecular and Cellular Proteomics 
 Journal of Lipid Research   Biochemistry and Molecular Biology Education 
Copyright © 1996 by the American Society for Biochemistry and Molecular Biology.