Is there a toxicological advantage for non-hyperbolic kinetics in cytochrome P450 catalysis? Functional allostery from "distributive catalysis"

The cytochrome P450s (CYPs) are the major enzymatic detoxification and drug metabolism system. Recently, it has become clear that several CYP isoforms exhibit positive and negative homotropic cooperativity. However, the toxicological implications of allosteric kinetics have not been considered, nor understood. The allosteric kinetics are particularly enigmatic in several respects. In many cases, CYPs bioactivate substrates to more toxic products, thus making it difficult to rationalize a functional advantage for positive cooperativity. Also, CYPs exhibit cooperativity with many structurally diverse ligands, in marked contrast to the specificity observed with other allosteric systems. Here, kinetic simulations are used to compare the probabilistic time- and concentration-dependent integrated toxicity function during conversion of substrate to product for CYP models exhibiting Michaelis-Menten (non-cooperative) kinetics, positive cooperativity, or negative cooperativity. The results demonstrate that, at low substrate concentrations, the slower substrate turnover afforded by cooperative CYPs compared with Michaelis-Menten enzymes can be a significant toxicological advantage, when toxic thresholds exist. When present, the advantage results from enhanced "distribution" of toxin in two pools, substrate and product, for an extended period, thus minimizing the chance that either exceeds its toxic threshold. At intermediate concentrations, the allosteric kinetics can be a modest advantage or modest disadvantage, depending on the kinetic parameters. However, at high substrate concentrations associated with a high probability of toxicity, fast turnover is desirable, and this advantage is provided also by the cooperative enzymes. For the positive homotropic cooperativity, the allosteric kinetics minimize the probability of toxicity over the widest range of system parameters. Furthermore, this apparent functional cooperativity is achieved without specific molecular recognition that is the hallmark of "traditional" allostery.     

The contribution of individual interchain interactions to the stabilization of the T and R states of Escherichia coli aspartate transcarbamoylase

Stabilization of the T and R allosteric states of Escherichia coli aspartate transcarbamoylase is governed by specific intra- and interchain interactions. The six interchain interactions between Glu-239 in one catalytic chain of one catalytic trimer with both Lys-164 and Tyr-165 of a different catalytic chain in the other catalytic trimer have been shown to be involved in the stabilization of the T state. In this study a series of hybrid versions of aspartate transcarbamoylase was studied to determine the minimum number of these Glu-239 interactions necessary to maintain homotropic cooperativity and the T allosteric state. Hybrids with zero, one, and two Glu-239 stabilizing interactions do not exhibit cooperativity, whereas the hybrids with three or more Glu-239 stabilizing interactions exhibit cooperativity. The hybrid enzymes with one or more of the Glu-239 stabilizing interactions also exhibit heterotropic interactions. Two hybrids with three Glu-239 stabilizing interactions, in different geometric relationships, had identical properties. From this and previous studies, it is concluded that the 239 stabilizing interactions play a critical role in the manifestation of homotropic cooperativity in aspartate transcarbamoylase by the stabilization of the T state of the enzyme. As substrate binding energy is utilized, more and more of the T state stabilizing interactions are relaxed, and finally the enzyme shifts to the R state. In the case of the Glu-239 stabilizing interactions more than three of the interactions must be broken before the enzyme shifts to the R state. The interactions between the catalytic and regulatory chains and between the two catalytic trimers of aspartate transcarbamoylase provide a global set of interlocking interactions that stabilize the T and R states of the enzyme. The substrate-induced local conformational changes observed in the structure of the isolated catalytic subunit drive the quaternary T to R transition of aspartate transcarbamoylase and functionally induced homotropic cooperativity.     

Morpheeins--a new structural paradigm for allosteric regulation

Classic models for the allosteric regulation of protein function consider an equilibrium among protein structures of constant oligomeric multiplicity. The morpheein (mor-phee'-in) concept expands this model to include a dynamic equilibrium of protein structures wherein a protein monomer can exist in more than one conformation and each monomer conformation dictates a different quaternary structure of finite multiplicity and different functionality. The morpheein concept provides a new framework for understanding allosteric regulation, kinetic cooperativity and hysteresis. Porphobilinogen synthase constitutes a prototype morpheein ensemble comprising several interconverting quaternary structure isoforms; one monomer conformation dictates assembly of a high-activity octamer, whereas an alternative monomer conformation dictates assembly of a low-activity hexamer. It is proposed here that the behavior of some other allosteric enzymes reflect dynamic morpheein equilibrium systems and six candidate proteins are enumerated.     

The pAR5 mutation and the allosteric mechanism of Escherichia coli aspartate carbamoyltransferase.

Mutation pAR5 replaces residues 145'-153' at the C terminus of the regulatory (r) chains of Escherichia coli ATCase by a new sequence of six residues. The mutated enzyme has been shown to lack substrate cooperativity and inhibition by CTP. Solution X-ray scattering curves demonstrate that, in the absence of ligands, its structure is intermediate between the T form and the R form. In the presence of N-phosphonacetyl-L-aspartate, the mutant is similar to the wild type. An examination of the crystal structure of unligated ATCase reveals that the mutated site is at an interface between r and catalytic (c) chains, which exists only in the T allosteric form. A computer simulation by energy minimization suggests that the pAR5 mutation destabilizes this interface and induces minor changes in the tertiary structure of r chains. The resulting lower stability of the T form explains the loss of substrate cooperativity. The lack of allosteric inhibition may be related to a new electrostatic interaction made in mutant r chains between the C-terminal carboxylate and a lysine residue of the allosteric domain.     

Creation of an allosteric phosphofructokinase starting with a nonallosteric enzyme. The case of dictyostelium discoideum phosphofructokinase.

An allosteric phosphofructokinase (PFK) was created by sequence manipulation of the nonallosteric enzyme from the slime mold Dictyostelium discoideum (DdPFK). Most amino acid residues proposed as important for catalytic and allosteric sites are conserved in DdPFK except for a few of them, and their reversion did not modify its kinetic behavior. However, deletions at the unique C-terminal extension of this PFK produced a markedly allosteric enzyme. Thus, a mutant lacking the last 26 C-terminal residues exhibited hysteresis in the time course, intense cooperativity (n(H) = 3.8), and a 200-fold decrease in the apparent affinity for fructose 6-phosphate (S(0.5) = 4500 microm), strong activation by fructose 2,6-bisphosphate (K(act) = 0.1 microm) and fructose 1,6-bisphosphate (K(act) = 40 microm), dependence on enzyme concentration, proton inhibition, and subunit association-dissociation in response to fructose 6-phosphate versus the nonhysteretic and hyperbolic wild-type enzyme (n(H) = 1.0; K(m) = 22 microm) that remained as a stable tetramer. Systematic deletions and point mutations at the C-tail region of DdPFK identified the last C-terminal residue, Leu(834), as critical to produce a nonallosteric enzyme. All allosteric mutants were practically insensitive to MgATP inhibition, suggesting that this effect does not involve the same allosteric transition as that responsible for fructose 6-phosphate cooperativity and fructose bisphosphate activation.     

Review: allostery in chaperonins

Chaperonins mediate protein folding in an ATP-dependent manner. ATP binding and hydrolysis by chaperonins are subject to both homotropic and heterotropic allosteric regulation. In the case of GroEL and CCT, homotropic regulation by ATP is manifested in nested cooperativity, which involves positive intra-ring cooperativity and negative inter-ring cooperativity in ATP binding. Both types of cooperativity are modulated by various heterotropic allosteric effectors, which include nonfolded proteins, ADP, Mg2+, monovalent ions such as K+, and cochaperonins in the case of type I chaperonins such as GroEL. Here, the allosteric properties of chaperonins are reviewed and new results of ours are presented with regard to allosteric effects of ADP. The role of allostery in the reaction cycle and folding function of chaperonins is discussed.     

Delineation of the allosteric mechanism of a cytidylyltransferase exhibiting negative cooperativity.

The dimeric enzyme CTP:glycerol-3-phosphate cytidylyltransferase (GCT) displays strong negative cooperativity between the first and second binding of its substrate, CTP. Using NMR to study the allosteric mechanism of this enzyme, we observe widespread chemical shift changes for the individual CTP binding steps. Mapping these changes onto the molecular structure allowed the formulation of a detailed model of allosteric conformational change. Upon the second step of ligand binding, NMR experiments indicate an extensive loss of conformational exchange broadening of the backbone resonances of GCT. This suggests that a fraction of the free energy of negative cooperativity is entropic in origin.     

Allosteric effects of four stereoisomers of a fused indole ring system with 3H-N-methylscopolamine and acetylcholine at M1-M4 muscarinic receptors

We previously demonstrated that brucine and some analogues allosterically enhance the affinity of ACh at muscarinic receptor subtypes M1, M3 or M4. Here we describe allosteric effects at human M1-M4 receptors of four stereoisomers of a pentacyclic structure containing features of the ring structure of brucine. All compounds inhibited 3H-NMS dissociation almost completely at all subtypes with slopes of 1, with similar affinity values at the 3H-NMS-occupied receptor to those estimated from equilibrium assays, consistent with the ternary complex allosteric model. Compound 1a showed positive cooperativity with H-NMS and small negative or neutral cooperativity with ACh at all subtypes. Its stereoisomer, 1b, showed strong negative cooperativity with both 3H-NMS and ACh across the subtypes. Compound 2a was positive with 3H-NMS at M2 and M4 receptors, neutral at M3 and negative at M1 receptors; it was negatively cooperative with ACh at all subtypes. Its stereoisomer, 2b, was neutral with 3H-NMS at M1 receptors and positive at the other subtypes; 2b was negatively cooperative with ACh at M1, M3 and M4 receptors but showed 3-fold positive cooperativity with ACh at M2 receptors. This latter result was confirmed with further 3H-NMS and 3H-ACh radioligand binding assays and with functional assays of ACh-stimulated 35S-GTPgammaS binding. These results provide the first well characterised instance of a positive enhancer of ACh at M2 receptors, and illustrate the difficulty of predicting such an effect.     

Control of kinetics by cooperative interactions

Cooperative effects in ligand binding and dissociation kinetics are much less investigated than steady state kinetics or equilibrium binding. Nevertheless, cooperativity in ligand binding leads necessarily to characteristic properties with respect to kinetic properties of the system. In case of positive cooperativity as found in oxygen binding proteins, a typical property is an autocatalytic ligand dissociation behavior leading to a time dependent, apparent ligand dissociation rate. To follow systematically the influence of the various potentially involved parameters on this characteristic property, simulations based on the simple MWC model were performed which should be relevant for all types of models based on the concept of an allosteric unit. In cases where the initial conformational distribution is very much dominated by the R-state, the intrinsic kinetic properties of the T-state are of minor influence for the observed ligand dissociation rate. Even for fast conformational transition rates, the R-state properties together with the size of the allosteric unit and the allosteric equilibrium constant define the shape of the curve. In such a case, a simplified model of the MWC-scheme (the irreversible n-chain model) is a good approximation of the full scheme. However, if in the starting conformational distribution some liganded T-molecules are present (a few percent is enough), the average off-rates can be significantly altered. Thus, the assignment of the initial rates to R-state properties has to be done with great care. However, if the R-state strongly dominates initially it is even possible to get an estimation of the lower limit for the number of interacting subunits from kinetic data: similar to the Hill-coefficient for equilibrium conditions, a measure for "kinetic cooperativity" can be derived by comparing initial and final ligand dissociation rates.     

Mapping pathways of allosteric communication in GroEL by analysis of correlated mutations

An interesting example of an allosteric protein is the chaperonin GroEL. It undergoes adenosine 5'-triphosphate-induced conformational changes that are reflected in binding of adenosine 5'-triphosphate with positive cooperativity within rings and negative cooperativity between rings. Herein, correlated mutations in chaperonins are analyzed to unravel routes of allosteric communication in GroEL and in its complex with its co-chaperonin GroES. It is shown that analysis of correlated mutations in the chaperonin family can provide information about pathways of allosteric communication within GroEL and between GroEL and GroES. The results are discussed in the context of available structural, genetic, and biochemical data concerning short- and long-range interactions in the GroE system.     

Excitation-contraction coupling gain and cooperativity of the cardiac ryanodine receptor: a modeling approach

During calcium-induced calcium-release, the ryanodine receptor (RyR) opens and releases large amounts of calcium from the sarcoplasmic reticulum into the cytoplasm of the myocyte. Recent experiments have suggested that cooperativity between the four monomers comprising the RyR plays an important role in the dynamics of the overall receptor. Furthermore, this cooperativity can be affected by the binding of FK506 binding protein, and hence, modulated by adrenergic stimulation through the phosphorylating action of protein kinase A. This has important implications for heart failure, where it has been hypothesized that RyR hyperphosphorylation, resulting in a loss of cooperativity, can lead to a persistent leak and a reduced sarcoplasmic-reticula content. In this study, we construct a theoretical model that examines the cooperativity via the assumption of an allosteric interaction between the four subunits. We find that the level of cooperativity, regulated by the binding of FK506 binding-protein, can have a dramatic effect on the excitation-contraction coupling gain and that this gain exhibits a clear maximum. These findings are compared to currently available data from different species and allows for an evaluation of the aforementioned heart-failure scenario.     

Inhibition of E. coli CTP synthase by the "positive" allosteric effector GTP

Cytidine 5'-triphosphate (CTP) synthase catalyzes the ATP-dependent formation of CTP from UTP using either ammonia or l-glutamine as the source of nitrogen. When glutamine is the substrate, GTP is required as a positive allosteric effector to promote catalysis of glutamine hydrolysis. We show that at concentrations exceeding approximately 0.15 mM, GTP actually behaves as a negative allosteric effector of E. coli CTP synthase, inhibiting glutamine-dependent CTP formation. In addition, GTP inhibits NH(3)-dependent CTP formation in a concentration-dependent manner. However, GTP does not inhibit the enzyme's intrinsic glutaminase activity. Although the activation of CTP synthase by GTP does not display cooperative behavior, inhibition of both CTP synthase-catalyzed ammonia- and glutamine-dependent CTP synthesis by GTP do exhibit positive cooperativity. These results suggest that GTP binding affects CTP synthase catalysis in two ways: it activates enzyme-catalyzed glutamine hydrolysis and it inhibits the utilization of NH(3) as a substrate by the synthase domain.     

Crystal structure and allosteric regulation of the cytoplasmic Escherichia coli L-asparaginase I

AnsA is the cytoplasmic asparaginase from Escherichia coli involved in intracellular asparagine utilization. Analytical ultracentifugation and X-ray crystallography reveal that AnsA forms a tetrameric structure as a dimer of two intimate dimers. Kinetic analysis of the enzyme reveals that AnsA is positively cooperative, displaying a sigmoidal substrate dependence curve with an [S](0.5) of 1 mM L-asparagine and a Hill coefficient (n(H)) of 2.6. Binding of L-asparagine to an allosteric site was observed in the crystal structure concomitant with a reorganization of the quarternary structure, relative to the apo enzyme. The carboxyl group of the bound asparagine makes salt bridges and hydrogen bonds to Arg240, while the N(delta2) nitrogen interacts with Thr162. Mutation of Arg240 to Ala increases the [S](0.5) value to 5.9 mM, presumably by reducing the affinity of the site for L-asparagine, although the enzyme retains cooperativity. Mutation of Thr162 to Ala results in an active enzyme with no cooperativity. Transmission of the signal from the allosteric site to the active site appears to involve subtle interactions at the dimer-dimer interface and relocation of Gln118 into the vicinity of the active site to position the probable catalytic water molecule. These data define the structural basis for the cooperative regulation of the intracellular asparaginase that is required for proper functioning within the cell.     

Stabilization of the R allosteric structure of Escherichia coli aspartate transcarbamoylase by disulfide bond formation

Here we report the first use of disulfide bond formation to stabilize the R allosteric structure of Escherichia coli aspartate transcarbamoylase. In the R allosteric state, residues in the 240s loop from two catalytic chains of different subunits are close together, whereas in the T allosteric state they are far apart. By substitution of Ala-241 in the 240s loop of the catalytic chain with cysteine, a disulfide bond was formed between two catalytic chains of different subunits. The cross-linked enzyme did not exhibit cooperativity for aspartate. The maximal velocity was increased, and the concentration of aspartate required to obtain one-half the maximal velocity, [Asp](0.5), was reduced substantially. Furthermore, the allosteric effectors ATP and CTP did not alter the activity of the cross-linked enzyme. When the disulfide bonds were reduced by the addition of 1,4-dithio-dl-threitol the resulting enzyme had kinetic parameters very similar to those observed for the wild-type enzyme and regained the ability to be activated by ATP and inhibited by CTP. Small-angle x-ray scattering was used to verify that the cross-linked enzyme was structurally locked in the R state and that this enzyme after reduction with 1,4-dithio-dl-threitol could undergo an allosteric transition similar to that of the wild-type enzyme. The complete abolition of homotropic and heterotropic regulation from stabilizing the 240s loop in its closed position in the R state, which forms the catalytically competent active site, demonstrates the significance that the quaternary structural change and closure of the 240s loop has in the functional mechanism of aspartate transcarbamoylase.     

Pyruvate kinase “Göttingen1,2”: Congenital hemolytic anemia,evidence of double heterozygosity,and lack of enzyme cooperativity

Summary Double heterozygosity of pyruvate kinase (PK) deficiency associated with hereditary hemolytic anemia is emphasized by studies of a kindred harboring two distinct mutant forms of this enzyme. The hematologically unaffected parents exhibit slightly reduced PK activity, a normal Hill coefficient, and a normal thermodynamic dissociation constant for the overall reaction. The paternal enzyme is characterized by normal substrate affinities and decreased activities with the substrate analogues CDP and GDP, whereas the maternal enzyme shows normal affinity for PEP, but an increased affinity for ADP and low thermostability. It is assumed that the erythrocytes of the parents contain a mixture of normal PK and a functionally abnormal isoenzyme, the latter differing between the parents. The two children suffer from hereditary hemolytic anemia. Their PK must be a combination of the mutant paternal and maternal isoenzymes, and their activities are reduced to about 30%. These enzymes are characterized by an increased affinity for PEP and a decreased affinity for ADP, a Hill coefficient of about 1 (indicating lack of cooperativity due to a loss of its allosteric properties), a decreased overall catalytic activity, and a higher resistance to heat denaturation. Further differences are observed in the SDS-gel electrophoresis between the two patients' enzymes. From the enzymological point of view it is impossible to characterize true PK variants in such double heterozygous cases which contain a combination of two different isoenzymes. The cause of chronic hemolysis appears to depend mainly on the loss of the allosteric properties, i.e., the lack of enzyme cooperativity.     

A loop involving catalytic chain residues 230-245 is essential for the stabilization of both allosteric forms of Escherichia coli aspartate transcarbamylase

The allosteric transition of Escherichia coli aspartate transcarbamylase involves significant alterations in structure at both the quaternary and tertiary levels. On the tertiary level, the 240s loop (residues 230-245 of the catalytic chain) repositions, influencing the conformation of Arg-229, a residue near the aspartate binding site. In the T state, Arg-229 is bent out of the active site and may be stabilized in this position by an interaction with Glu-272. In the R state, the conformation of Arg-229 changes, allowing it to interact with the beta-carboxylate of aspartate, and is stabilized in this position by a specific interaction with Glu-233. In order to ascertain the function of Arg-229, Glu-233, and Glu-272 in the catalytic and cooperative interactions of the enzyme, three mutant enzymes were created by site-specific mutagenesis. Arg-229 was replaced by Ala, while both Glu-233 and Glu-272 were replaced by Ser. The Arg-229----Ala and Glu-233----Ser enzymes exhibit 10,000-fold and 80-fold decreases in maximal activity, respectively, and they both exhibit a 2-fold increase in the aspartate concentration at half the maximal observed velocity, [S]0.5. The Arg-229----Ala enzyme still exhibits substantial homotropic cooperativity, but all cooperativity is lost in the Glu-233----Ser enzyme. The Glu-233----Ser enzyme also shows a 4-fold decrease in the carbamyl phosphate [S]0.5, while the Arg-229----Ala enzyme shows no change in the carbamyl phosphate [S]0.5 compared to the wild-type enzyme. The Glu-272 to Ser mutation results in a slight reduction in maximal activity, an increase in [S]0.5 for both aspartate and carbamyl phosphate, and reduced cooperativity. Analysis of the isolated catalytic subunits from these three mutant enzymes reveals that in each case the changes in the kinetic properties of the isolated catalytic subunit are similar to the changes caused by the mutation in the holoenzyme. PALA was able to activate the Glu-233----Ser enzyme, at low aspartate concentrations, even though the mutant holoenzyme did not exhibit any cooperativity, indicating that cooperative interactions still exist between the active sites in this enzyme. It is proposed that Glu-233 of the 240s loop helps create the high-activity-high-affinity R state by positioning the side chain of Arg-229 for aspartate binding while Glu-272 helps stabilize the low-activity-low-affinity T state by positioning the side chain of Arg-229 so that it cannot interact with aspartate.(ABSTRACT TRUNCATED AT 400 WORDS)     

A quantitative model for allosteric control of purine reduction by murine ribonucleotide reductase

The reduction of purine nucleoside diphosphates by murine ribonucleotide reductase requires catalytic (R1) and free radical-containing (R2) enzyme subunits and deoxynucleoside triphosphate allosteric effectors. A quantitative 16 species model is presented, in which all pertinent equilibrium constants are evaluated, that accounts for the effects of the purine substrates ADP and GDP, the deoxynucleoside triphosphate allosteric effectors dGTP and dTTP, and the dimeric murine R2 subunit on both the quaternary structure of murine R1 subunit and the dependence of holoenzyme (R1(2)R2(2)) activity on substrate and effector concentrations. R1, monomeric in the absence of ligands, dimerizes in the presence of substrate, effectors, or R2(2) because each of these ligands binds R1(2) with higher affinity than R1 monomer. This leads to apparent positive heterotropic cooperativity between substrate and allosteric effector binding that is not observed when binding to the dimeric protein itself is evaluated. Allosteric activation results from an increase in k(cat) for substrate reduction upon binding of the correct effector, rather than from heterotropic cooperativity between effector and substrate. Neither the allosteric site nor the active site displays nucleotide base specificity: dissociation constants for dGTP and dTTP are nearly equivalent and K(m) and k(cat) values for both ADP and GDP are similar. R2(2) binding to R1(2) shows negative heterotropic cooperativity vis-à-vis effectors but positive heterotropic cooperativity vis-à-vis substrates. Binding of allosteric effectors to the holoenzyme shows homotropic cooperativity, suggestive of a conformational change induced by activator binding. This is consistent with kinetic results indicating full dimer activation upon binding a single equivalent of effector per R1(2)R2(2).     

Is Disrupted Nucleotide-Substrate Cooperativity a Common Trait for Cushing's Syndrome Driving Mutations of Protein Kinase A?

Somatic mutations in the PRKACA gene encoding the catalytic α subunit of protein kinase A (PKA-C) are responsible for cortisol-producing adrenocortical adenomas. These benign neoplasms contribute to the development of Cushing's syndrome. The majority of these mutations occur at the interface between the two lobes of PKA-C and interfere with the enzyme's ability to recognize substrates and regulatory (R) subunits, leading to aberrant phosphorylation patterns and activation. Rarely, patients with similar phenotypes carry an allosteric mutation, E31V, located at the C-terminal end of the αA-helix and adjacent to the αC-helix, but structurally distinct from the PKA-C/R subunit interface mutations. Using a combination of solution NMR, thermodynamics, kinetic assays, and molecular dynamics simulations, we show that the E31V allosteric mutation disrupts central communication nodes between the N- and C- lobes of the enzyme as well as nucleotide-substrate binding cooperativity, a hallmark for kinases' substrate fidelity and regulation. For both orthosteric (L205R and W196R) and allosteric (E31V) Cushing’s syndrome mutants, the loss of binding cooperativity is proportional to the density of the intramolecular allosteric network. This structure–activity relationship suggests a possible common mechanism for Cushing's syndrome driving mutations in which decreased nucleotide/substrate binding cooperativity is linked to loss in substrate fidelity and dysfunctional regulation.     

Recombinant hemoglobins with low oxygen affinity and high cooperativity

By introducing an additional H-bond in the alpha(1)beta(2) subunit interface or altering the charge properties of the amino acid residues in the alpha(1)beta(1) subunit interface of the hemoglobin molecule, we have designed and expressed recombinant hemoglobins (rHbs) with low oxygen affinity and high cooperativity. Oxygen-binding measurements of these rHbs under various experimental conditions show interesting properties in response to pH (Bohr effect) and allosteric effectors. Proton nuclear magnetic resonance studies show that these rHbs can switch from the oxy (or CO) quaternary structure (R) to the deoxy quaternary structure (T) without changing their ligation states upon addition of an allosteric effector, inositol hexaphosphate, and/or reduction of the ambient temperature. These results indicate that if we can provide extra stability to the T state of the hemoglobin molecule without perturbing its R state, we can produce hemoglobins with low oxygen affinity and high cooperativity. Some of these rHbs are also quite stable against autoxidation compared to many of the known abnormal hemoglobins with altered oxygen affinity and cooperativity. These results have provided new insights into the structure-function relationship in hemoglobin.