Methemoglobin--it's not just blue: a concise review

Hemoglobin has functions besides carrying oxygen to the tissues, and regulates vascular tone and inflammation via a redox couple with methemoglobin. Hemoglobin has iron in the reduced valance Fe(II) and methemoglobin has iron in the oxidized valance Fe (III), with a free energy capable of producing water from oxygen. In generating methemoglobin the couple functions as a nitrite reductase. The degree of oxidation of hemoglobin senses the oxygen level in the blood and uses its ability to produce nitric oxide from nitrite to control vascular tone, increasing blood flood when the proportion of oxygenated hemoglobin falls. Additional cardiovascular damage is produced by methemoglobin mediated oxidation of light density lipoproteins, accelerating arteriosclerosis. In addition, the release of heme from methemoglobin is an important factor in inflammation. These physiologic functions are paralleled by the well-described role in the oxidation of various drugs resulting in methemoglobinemia.     

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.     

Accessible surface areas as a measure of the thermodynamic parameters of hydration of peptides

A method is described for the inclusion of the effects of hydration in empirical conformational energy computations on polypeptides. The free energy of hydration is composed of additive contributions of various functional groups. The hydration of each group is assumed to be proportional to the accessible surface area of the group. The constants of proportionality, representing the free energy of hydration per unit area of accessible surface, have been evaluated for seven classes of groups (occurring in peptides) by least-squares fitting to experimental free energies of solution of small monofunctional aliphatic and aromatic molecules. The same method has also been applied to the modeling of the enthalpy and heat capacity of hydration, each of which is computed from the accessible surface area.     

Crystalline variational methods

A surface free energy function is defined to be crystalline if its Wulff shape (the equilibrium crystal shape) is a polyhedron. All the questions that one considers for the area functional, where the surface free energy per unit area is 1 for all normal directions, can be considered for crystalline surface free energies. Such questions are interesting for both mathematical and physical reasons. Methods from the geometric calculus of variations are useful for studying a number of such questions; a survey of some of the results is given.     

Resonance Raman Spectra of Cobalt-Substituted Hemoglobin: Cooperativity and Displacement of the Cobalt Atom upon Oxygenation

The resonance Raman spectra of oxy and deoxy cobalt-substituted hemoglobin (CoHb) are reported. Comparison of these spectra to those of hemoglobin, methemoglobin, cytochrome c, and model cobalt porphyrin complexes suggests that the displacement of the cobalt atom upon oxygenation of CoHb is no greater than the out-of-plane distance in five-coordinate Co(II) porphyrins, 0.15 A. Combining this distance with the expected contraction of the cobalt-histidine bond, Ibers has estimated a maximum displacement of 0.37 A for the proximal histidine with respect to the heme plane upon oxygenation, about one-third the corresponding distance estimated for iron hemoglobin. The free energy of cooperativity for cobalt hemoglobin is also estimated to be one-third that of iron hemoglobin. These results are therefore consistent with Hopfield's distributed energy model, which predicts proportionality between proximal histidine displacement and the free energy of cooperativity. By implication they support Perutz's trigger mechanism for cooperativity.     

Free energy of burying hydrophobic residues in the interface between protein subunits

We have obtained an experimental estimate of the free energy change associated with variations at the interface between protein subunits, a subject that has raised considerable interest since the concept of accessible surface area was introduced by Lee and Richards [Lee, B. & Richards, F. M. (1971) J. Mol. Biol. 55, 379–400]. We determined by analytical ultracentrifugation the dimer–tetramer equilibrium constant of five single and three double mutants of human Hb. One mutation is at the stationary α₁β₁ interface, and all of the others are at the sliding α₁β₂ interface where cleavage of the tetramer into dimers and ligand-linked allosteric changes are known to occur. A surprisingly good linear correlation between the change in the free energy of association of the mutants and the change in buried hydrophobic surface area was obtained, after corrections for the energetic cost of losing steric complementarity at the αβ dimer interface. The slope yields an interface stabilization free energy of −15 ± 1.2 cal/mol upon burial of 1 Ų of hydrophobic surface, in very good agreement with the theoretical estimate given by Eisenberg and McLachlan [Eisenberg, D. & McLachlan, A. D. (1986) Nature (London) 319, 199–203].     

The role of hydrophobic interactions in initiation and propagation of protein folding

Globular proteins fold by minimizing the nonpolar surface that is exposed to water, while simultaneously providing hydrogen-bonding interactions for buried backbone groups, usually in the form of secondary structures such as alpha-helices, beta-sheets, and tight turns. A primary thermodynamic driving force for the formation of globular structure is thus the sequestration of nonpolar groups, but the correlation between the parts of proteins that are observed to fold first (termed folding initiation sites) and the "hydrophobicity" (as customarily defined) of the amino acids in these regions has been quite weak. It has previously been noted that many amino acid side chains contain considerable nonpolar sections, even if they also contain polar or charged groups. For example, a lysine side chain contains four methylenes, which may undergo hydrophobic interactions if the charged epsilon-NH(3)(+) group is salt-bridged or hydrogen-bonded. Folding initiation sites might therefore contain not only accepted "hydrophobic" amino acids, but also larger charged side chains. Recent experiments on the folding of mutant apomyoglobins provides corroboration for models based on the hypothesis that folding initiation sites arise from hydrophobic interactions. A near-perfect correlation was observed between the areas of the molecule that are present in the burst-phase kinetic intermediate and both the free energy of formation of hydrophobic initiation sites and the parameter "average area buried upon folding," which pinpoints large side chains, even those containing charged or polar portions. These results provide a putative mechanism for the control of protein-folding initiation and growth by polar/nonpolar sequence propensity alone.     

Catalysis of soluble hemoglobin oxidation by free iron on sickle red cell membranes

Abnormal deposition of hemichrome on the inner aspect of the sickle red cell membrane promotes premature cell demise. The steps proximate to hemichrome formation in these cells are poorly understood. To test the hypothesis that the pathologic deposits of free ferric iron located on the inner aspect of sickle cell membranes would be redox active and promote oxidation of soluble oxyhemoglobin, we incubated native versus iron-stripped sickle or normal ghost membranes with oxyhemoglobin S. We found that sickle membranes exerted an exaggerated effect on methemoglobin formation in solution, an effect completely accounted for by their abnormal content of free iron. This ability of sickle membranes to promote hemoglobin oxidation was not diminished by catalase or by presence of a high-affinity, iron-inactivating chelator that is unable to remove membrane iron. Examination of those membranes likewise revealed that their free iron content promoted deposition of additional heme-protein. These results establish that the potential redox couple formed by membrane-associated ferric iron and cytoplasmic oxyhemoglobin is promotive of hemoglobin oxidation and deposition of hemichrome on the membrane. This predicts that removal of pathologic membrane iron might help prevent the detrimental formation of methemoglobin and hemichrome in vivo, insofar as this is accelerated by transition metal.     

Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water

The hydrophobic interaction, the tendency for nonpolar molecules to aggregate in solution, is a major driving force in biology. In a direct approach to the physical basis of the hydrophobic effect, nanosecond molecular dynamics simulations were performed on increasing numbers of hydrocarbon solute molecules in water-filled boxes of different sizes. The intermittent formation of solute clusters gives a free energy that is proportional to the loss in exposed molecular surface area with a constant of proportionality of 45 +/- 6 cal/mol A(2). The molecular surface area is the envelope of the solute cluster that is impenetrable by solvent and is somewhat smaller than the more traditional solvent-accessible surface area, which is the area transcribed by the radius of a solvent molecule rolled over the surface of the cluster. When we apply a factor relating molecular surface area to solvent-accessible surface area, we obtain 24 cal/mol A(2). Ours is the first direct calculation, to our knowledge, of the hydrophobic interaction from molecular dynamics simulations; the excellent qualitative and quantitative agreement with experiment proves that simple van der Waals interactions and atomic point-charge electrostatics account for the most important driving force in biology.     

Prediction of protein-folding mechanisms from free-energy landscapes derived from native structures

Guided by recent experimental results suggesting that protein-folding rates and mechanisms are determined largely by native-state topology, we develop a simple model for protein folding free-energy landscapes based on native-state structures. The configurations considered by the model contain one or two contiguous stretches of residues ordered as in the native structure with all other residues completely disordered; the free energy of each configuration is the difference between the entropic cost of ordering the residues, which depends on the total number of residues ordered and the length of the loop between the two ordered segments, and the favorable attractive interactions, which are taken to be proportional to the total surface area buried by the ordered residues in the native structure. Folding kinetics are modeled by allowing only one residue to become ordered/disordered at a time, and a rigorous and exact method is used to identify free-energy maxima on the lowest free-energy paths connecting the fully disordered and fully ordered configurations. The distribution of structure in these free-energy maxima, which comprise the transition-state ensemble in the model, are reasonably consistent with experimental data on the folding transition state for five of seven proteins studied. Thus, the model appears to capture, at least in part, the basic physics underlying protein folding and the aspects of native-state topology that determine protein-folding mechanisms.     

Dynamic hydration shell restores Kauzmann's 1959 explanation of how the hydrophobic factor drives protein folding

Kauzmann''s explanation of how the hydrophobic factor drives protein folding is reexamined. His explanation said that hydrocarbon hydration shells are formed, possibly of clathrate water, and they explain why hydrocarbons have uniquely low solubilities in water. His explanation was not universally accepted because of skepticism about the clathrate hydration shell. A revised version is given here in which a dynamic hydration shell is formed by van der Waals (vdw) attraction, as proposed in 1985 by Jorgensen et al. [Jorgensen WL, Gao J, Ravimohan C (1985) J Phys Chem 89:3470–3473]. The vdw hydration shell is implicit in theories of hydrophobicity that contain the vdw interaction between hydrocarbon C and water O atoms. To test the vdw shell model against the known hydration energetics of alkanes, the energetics should be based on the Ben-Naim standard state (solute transfer between fixed positions in the gas and liquid phases). Then the energetics are proportional to n, the number of water molecules correlated with an alkane by vdw attraction, given by the simulations of Jorgensen et al. The energetics show that the decrease in entropy upon hydration is the root cause of hydrophobicity; it probably results from extensive ordering of water molecules in the vdw shell. The puzzle of how hydrophobic free energy can be proportional to nonpolar surface area when the free energy is unfavorable and the only known interaction (the vdw attraction) is favorable, is resolved by finding that the unfavorable free energy is produced by the vdw shell.When Kauzmann reviewed in 1959 (1) the possible sources of the free energy needed to drive protein folding, he found that that the known factors are not sufficient. He asked what the missing factor could be, and he ruled out peptide H-bonds because they do not provide enough free energy, based on Schellman''s (2) analysis of the problem in 1955. Then he discovered the previously unknown hydrophobic factor after observing that known DG values for transfer of hydrocarbons out of water into other solvents could supply the missing free energy. Then he needed to assume that the interior of a folded protein is water-free and the nonpolar side chains are buried inside the protein as folding occurs. In 1960 the 2-Å structure of myoglobin by Kendrew et al. (3) confirmed these predictions.     

Interfacial free energy and the hydrophobic effect

Interfacial free energies demonstrate clearly that the antipathy between hydrocarbon and water rests on the strong attraction of water for itself. However, the unfavorable free energy associated with this antipathy, per unit area of contact between bulk hydrocarbon and water, is about 3-fold larger than a similar figure derived from solubility data per unit area of contact between a single dissolved hydrocarbon molecule and water. The discrepancy illustrates the difficulty in applying macroscopic concepts such as “interfacial surface” at the molecular level and can be formally resolved, at least qualitatively, by the predicted effect of surface curvature on surface tension.     

Conservation and relative importance of residues across protein-protein interfaces

A core region surrounded by a rim characterizes biological interfaces. We ascertain the importance of the core by showing the sequence entropies of the residues comprising the core to be smaller than those in the rim. Such a distinction is not seen in the 2-fold-related, nonphysiological interfaces formed in crystal lattices of monomeric proteins, thereby providing a procedure for characterizing the oligomeric state from crystal structures of protein molecules. This method is better than those that rely on the comparison of the sequence entropies in the interface and the rest of the protein surface, especially in cases where the surface harbors additional binding sites. To a good approximation there is a correlation between the accessible surface area lost because of complexation and DeltaDeltaG values obtained through alanine-scanning mutagenesis (26-38 cal per A(2) of the surface buried) for residues located in the core, a relationship that is not discernable for rim residues. If, however, a residue participates in hydrogen bonding across the interface, the extent of stabilization is 52 cal/mol per 1 A(2) of the nonpolar surface area buried by the residue. As opposed to an amino acid classification used earlier, an environment-based grouping of residues yields a better discrimination in the sequence entropy between the core and the rim.     

Analytical approximation to the accessible surface area of proteins

We propose an analytical substitute to the geometrical construction that is commonly used in calculating the protein surface area that is accessible to the solvent. A statistical approach leads to an expression of accessible surface areas as a function of distances between pairs of atoms or of residues in the protein structure, assuming only that these atoms or residues are randomly distributed in space but not penetrating each other. This function gives good estimates of the accessible surface area and of the area buried in subunit contacts for a number of proteins. Its evaluation is very fast, and the function can be differentiated, which opens the way to new applications of accessibility measurements in the study of proteins. As an example, we show that the presence of domains is easily detected by an automatic procedure based on surface areas only.     

Methemoglobin Reduction: Studies of the Interaction between Cell Populations and of the Role of Methylene Blue

1. It has been shown that mixtures of normal and G-6-PD deficient nitrite-treated erythrocytes reduced methemoglobin at a rate considerably morerapid than that computed from their individual rates of methemoglobin reduction.

2. Differential agglutination studies demonstrated that when normal,methemoglobin-containing cells in such a mixture have reduced all of theirmethemoglobin, they facilitate methemoglobin reduction in G-6-PD deficienterythrocytes.

3. The same effect could be observed in other mixtures of cells, (e.g., normal with normal, G-6-PD deficient with G-6-PD deficient, etc.) and evenwith highly purified methemoglobin solutions.

4. This effect could be observed in the presence of methylene blue but notin the presence of another redox dye, Nile blue sulfate.

5. Lactate served as an effective substrate for methemoglobin reduction.Methemoglobin reduction by lactate was enhanced by methylene blue butnot by Nile blue sulfate.

6. In mixtures of normal and G-6-PD deficient erythrocytes, no deficit inthe rate of accumulation of lactate was found. This indicates that the mechanism of enhancement of methemoglobin reduction is not the diffusion oflactate from non-methemoglobin-containing cells to methemoglobin-containing cells.

7. It was demonstrated that leukomethylene blue could reduce highlypurified solutions of methemoglobin in the absence of the enzyme "methemoglobin reductase."

8. The possible mechanism by which non-methemoglobin-containing cellsmay reduce methemoglobin in methemoglobin-containing cells is discussed.It seems most probable that leukomethylene blue is the mediator of the effect.This implies, contrary to earlier suggestions, that "methemoglobin reductase"acts prior to the reduction of methylene blue in the electron transport chain.

Submitted on February 13, 1963 Accepted on April 3, 1963     

Endothelial-cell heme uptake from heme proteins: induction of sensitization and desensitization to oxidant damage.

Iron-derived reactive oxygen species are implicated in the pathogenesis of various vascular disorders including atherosclerosis, vasculitis, and reperfusion injury. The present studies examine whether heme, when liganded to physiologically relevant proteins as in hemoglobin, can provide potentially damaging iron to intact endothelium. We demonstrate that reduced ferrohemoglobin, while relatively innocuous to cultured endothelial cells, when oxidized to ferrihemoglobin (methemoglobin), greatly amplifies oxidant (H2O2)-mediated endothelial-cell injury. Drawing upon our previous observation that free heme similarly primes endothelium for oxidant damage, we posited that methemoglobin, but not ferrohemoglobin, releases its hemes that can then be incorporated into endothelial cells. In support, cultured endothelial cells exposed to methemoglobin--in contrast to exposure to ferrohemoglobin, cytochrome c, or metmyoglobin--rapidly increased their heme oxygenase mRNA and enzyme activity, thereby supporting heme uptake; ferritin production was also markedly increased after such exposure, thus attesting to eventual incorporation of Fe. These cellular methemoglobin effects were inhibited by the heme-scavenging protein hemopexin and by haptoglobin or cyanide, agents that strengthen the liganding between heme and globin. If the endothelium is exposed to methemoglobin for a more prolonged period (16 hr), it accumulates large amounts of ferritin; concomitantly, and presumably associated with iron sequestration by this protein, the endothelium converts from hypersusceptible to hyperresistant to oxidative damage. We conclude that when oxidation of hemoglobin facilitates release of its heme groups, catalytically active iron is provided to neighboring tissue environments. The effect of this relinquished heme on the vasculature is determined both by extracellular factors--i.e., plasma proteins, such as haptoglobin and hemopexin--as well as intracellular factors, including heme oxygenase and ferritin. Acutely, if both extra- and intracellular defenses are overwhelmed, cellular toxicity arises; chronically, when ferritin is induced, resistance to oxidative injury may supervene.     

Free energy profiles from single-molecule pulling experiments

Nonequilibrium pulling experiments provide detailed information about the thermodynamic and kinetic properties of molecules. We show that unperturbed free energy profiles as a function of molecular extension can be obtained rigorously from such experiments without using work-weighted position histograms. An inverse Weierstrass transform is used to relate the system free energy obtained from the Jarzynski equality directly to the underlying molecular free energy surface. An accurate approximation for the free energy surface is obtained by using the method of steepest descent to evaluate the inverse transform. The formalism is applied to simulated data obtained from a kinetic model of RNA folding, in which the dynamics consists of jumping between linker-dominated folded and unfolded free energy surfaces.     

Directed assembly of cell-laden microgels for fabrication of 3D tissue constructs

We present a bottom-up approach to direct the assembly of cell-laden microgels to generate tissue constructs with tunable microarchitecture and complexity. This assembly process is driven by the tendency of multiphase liquid-liquid systems to minimize the surface area and the resulting surface free energy between the phases. We demonstrate that shape-controlled microgels spontaneously assemble within multiphase reactor systems into predetermined geometric configurations. Furthermore, we characterize the parameters that influence the assembly process, such as external energy input, surface tension, and microgel dimensions. Finally, we show that multicomponent cell-laden constructs could be generated by assembling microgel building blocks and performing a secondary cross-linking reaction. This bottom-up approach for the directed assembly of cell-laden microgels provides a powerful and highly scalable approach to form biomimetic 3D tissue constructs and opens a paradigm for directing the assembly of mesoscale materials.     

The reduction of methemoglobin levels by antioxidants

Preventing the oxidation of hemoglobin in solution is one of the major requirements for the successful production and long-term storage of hemoglobin-based blood substitutes. To this end we have studied the effects of antioxidants on the rate of methemoglobin formation and disappearance in solutions of human and bovine hemoglobin at 4 degrees C and 37 degrees C. Ascorbate and desferal (5 mM) were observed to act as prooxidants, increasing the rate of methemoglobin formation at 37 degrees C. Trehalose, mannitol, glucose, and EDTA (5 mM) had no significant effect. Glutathione and NADH (10 mM) were the most effective antioxidants tested, causing a significant decrease in the rate of methemoglobin formation at 37 degrees C for periods of up to 50 hours. The combination of these antioxidants in bovine hemoglobin at 4 degrees C resulted in the reduction of methemoglobin levels to nearly undetectable levels in approximately 150 hours. In addition, NADH and glutathione were found to reduce methemoglobin levels to 10% over a period of 100 hours in a sample of human hemoglobin that had been stored at 4 degrees C for one year and had 60% methemoglobin. These results suggest that the prevention and reversal of methemoglobin formation during the long-term storage of hemoglobin solutions and hemoglobin-based blood substitutes may now be possible.     

Unnecessary delay of diagnosis of buried bumper syndrome resulting in surgery

Percutaneous endoscopic catheter gastrostomy (PEG) is a convenient way to supply enteral nutrition for patients with swallowing disorders. One rare complication of PEG is the buried bumper syndrome where gastric mucosa overgrows the internal bumper and prevents free flow of the feeding solution. As a consequence, the application of enteral feeding has to be stopped until a free outflow is re-established. We report a case of buried bumper where symptoms were misinterpreted for several months as PEG stoma infection by the homecare service. This led to a vastly delayed diagnosis and treatment. As endoscopic intervention was unsuccessful, surgical PEG removal was required. In consequence, we recommend early endoscopic exploration in cases with prolonged inflammatory signs at the PEG stoma site in order to avoid misdiagnosis of buried bumper syndrome and to allow timely endoscopic intervention.