Evolutionarily consistent families in SCOP: sequence, structure and function

Background

HutZ is the sole heme storage protein identified in the pathogenic bacterium Vibrio cholerae and is required for optimal heme utilization. However, no heme oxygenase activity has been observed with this protein. Thus far, HutZ??s structure and heme-binding mechanism are unknown.

Results

We report the first crystal structure of HutZ in a homodimer determined at 2.0 ? resolution. The HutZ structure adopted a typical split-barrel fold. Through a docking study and site-directed mutagenesis, a heme-binding model for the HutZ dimer is proposed. Very interestingly, structural superimposition of HutZ and its homologous protein HugZ, a heme oxygenase from Helicobacter pylori, exhibited a structural mismatch of one amino acid residue in ??6 of HutZ, although residues involved in this region are highly conserved in both proteins. Derived homologous models of different single point variants with model evaluations suggested that Pro¹⁴⁰ of HutZ, corresponding to Phe²¹⁵ of HugZ, might have been the main contributor to the structural mismatch. This mismatch initiates more divergent structural characteristics towards their C-terminal regions, which are essential features for the heme-binding of HugZ as a heme oxygenase.

Conclusions

HutZ??s deficiency in heme oxygenase activity might derive from its residue shift relative to the heme oxygenase HugZ. This residue shift also emphasized a limitation of the traditional template selection criterion for homology modeling.     

Bis-methionyl coordination in the crystal structure of the heme-binding domain of the streptococcal cell surface protein Shp

Surface proteins Shr, Shp, and the ATP-binding cassette (ABC) transporter HtsABC are believed to make up the machinery for heme uptake in Streptococcus pyogenes. Shp transfers its heme to HtsA, the lipoprotein component of HtsABC, providing the only experimentally demonstrated example of direct heme transfer from a surface protein to an ABC transporter in Gram-positive bacteria. To understand the structural basis of heme transfer in this system, the heme-binding domain of Shp (Shp¹⁸⁰) was crystallized, and its structure determined to a resolution of 2.1 Å. Shp¹⁸⁰ exhibits an immunoglobulin-like β-sandwich fold that has been recently found in other pathogenic bacterial cell surface heme-binding proteins, suggesting that the mechanisms of heme acquisition are conserved. Shp shows minimal amino acid sequence identity to these heme-binding proteins and the structure of Shp¹⁸⁰ reveals a unique heme-iron coordination with the axial ligands being two methionine residues from the same Shp molecule. A negative electrostatic surface of protein structure surrounding the heme pocket may serve as a docking interface for heme transfer from the more basic outer cell wall heme receptor protein Shr. The crystal structure of Shp¹⁸⁰ reveals two exogenous, weakly bound hemins, which form a large interface between the two Shp¹⁸⁰ molecules in the asymmetric unit. These “extra” hemins form a stacked pair with a structure similar to that observed previously for free hemin dimers in aqueous solution. The propionates of the protein-bound heme coordinate to the iron atoms of the exogenous hemin dimer, contributing to the stability of the protein interface. Gel filtration and analytical ultracentrifugation studies indicate that both full-length Shp and Shp¹⁸⁰ are monomeric in dilute aqueous solution.     

Discovery and Characterization of HemQ: AN ESSENTIAL HEME BIOSYNTHETIC PATHWAY COMPONENT*

Here we identify a previously undescribed protein, HemQ, that is required for heme synthesis in Gram-positive bacteria. We have characterized HemQ from Bacillus subtilis and a number of Actinobacteria. HemQ is a multimeric heme-binding protein. Spectroscopic studies indicate that this heme is high spin ferric iron and is ligated by a conserved histidine with the sixth coordination site available for binding a small molecule. The presence of HemQ along with the terminal two pathway enzymes, protoporphyrinogen oxidase (HemY) and ferrochelatase, is required to synthesize heme in vivo and in vitro. Although the exact role played by HemQ remains to be characterized, to be fully functional in vitro it requires the presence of a bound heme. HemQ possesses minimal peroxidase activity, but as a catalase it has a turnover of over 10⁴ min⁻¹. We propose that this activity may be required to eliminate hydrogen peroxide that is generated by each turnover of HemY. Given the essential nature of heme synthesis and the restricted distribution of HemQ, this protein is a potential antimicrobial target for pathogens such as Mycobacterium tuberculosis.     

1H, 13C, 15N backbone and side chain NMR resonance assignments of the N-terminal NEAr iron transporter domain 1 (NEAT 1) of the hemoglobin receptor IsdB of Staphylococcus aureus

Staphylococcus aureus is an opportunistic pathogen that causes skin and severe infections in mammals. Critical to S. aureus growth is its ability to scavenge iron from host cells. To this effect, S. aureus has evolved a sophisticated pathway to acquire heme from hemoglobin (Hb) as a preferred iron source. The pathway is comprised of nine iron-regulated surface determinant (Isd) proteins involved in heme capture, transport, and degradation. A key protein of the heme acquisition pathway is the surface-anchored hemoglobin receptor protein IsdB, which is comprised of two NEAr transporter (NEAT) domains that act in concert to bind Hb and extract heme for subsequent transfer to downstream acquisition pathway proteins. Despite significant advances in the structural knowledge of other Isd proteins, the structural mechanisms and molecular basis of the IsdB-mediated heme acquisition process are not well understood. In order to provide more insights into the mode of function of IsdB, we have initiated NMR structural studies of the first NEAT domain of IsdB (IsdB^N1). Herein, we report the near complete ¹H, ¹³C and ¹⁵N resonance assignments of backbone and side chain atoms, and the secondary structural topology of the 148-residue IsdB NEAT 1 domain. The NMR results are consistent with the presence of eight β-strands and one α-helix characteristic of an immunoglobulin-like fold observed in other NEAT domain family proteins. This work provides a solid framework to obtain atomic-level insights toward understanding how IsdB mediates IsdB-Hb protein–protein interactions critical for heme capture and transfer.     

Coordination modes of tyrosinate-ligated catalase-type heme enzymes: magnetic circular dichroism studies of Plexaura homomalla allene oxide synthase, Mycobacterium avium ssp. paratuberculosis protein-2744c, and bovine liver catalase in their ferric and ferrous states

Bovine liver catalase (BLC), catalase-related allene oxide synthase (cAOS) from Plexaura homomalla, and a recently isolated protein from the cattle pathogen Mycobacterium avium ssp. paratuberculosis (MAP-2744c (MAP)) are all tyrosinate-ligated heme enzymes whose crystal structures have been reported. cAOS and MAP have low (< 20%) sequence similarity to, and significantly different catalytic functions from, BLC. cAOS transforms 8R-hydroperoxy-eicosatetraenoic acid to an allene epoxide, whereas the MAP protein is a putative organic peroxide-dependent peroxidase. To elucidate factors influencing the functions of these and related heme proteins, we have investigated the heme iron coordination properties of these tyrosinate-ligated heme enzymes in their ferric and ferrous states using magnetic circular dichroism and UV-visible absorption spectroscopy. The MAP protein shows remarkable spectral similarities to cAOS and BLC in its native Fe(III) state, but clear differences from ferric proximal heme ligand His93Tyr Mb (myoglobin) mutant, which may be attributed to the presence of an Arg⁺-N^ω-H···¯O-Tyr (proximal heme axial ligand) hydrogen bond in the first three heme proteins. Furthermore, the spectra of Fe(III)-CN¯, Fe(III)-NO, Fe(II)-NO (except for five-coordinate MAP), Fe(II)-CO, and Fe(II)-O₂ states of cAOS and MAP, but not H93Y Mb, are also similar to the corresponding six-coordinate complexes of BLC, suggesting that a tyrosinate (Tyr-O¯) is the heme axial ligand trans to the bound ligands in these complexes. The Arg⁺-N^ω-H to ¯O-Tyr hydrogen bond would be expected to modulate the donor properties of the proximal tyrosinate oxyanion and, combined with the subtle differences in the catalytic site structures, affect the activities of cAOS, MAP and BLC.     

Methionine-oxidized horse heart cytochromes c. II. Conformation and heme configuration

The two products from the reaction of horse heart ferricytochrome c with Chloramine-T, the F^III and F^II CT-cytochromes, contain modification of the methionines to methionine sulfoxides, but they are distinct in their physiological functions. Conformational and heme-configurational characterization of the two CT-cytochromes has been carried out by using absorption, circular dichroism, fluorescence, proton magnetic resonance, and resonance Raman spectroscopy. The pH-absorption spectroscopic behavior, thermal stability, and ionization of the phenolic hydroxyls have also been reported. Spectroscopic studies of the heme c fragment, H8, in the presence of dimethylsulfoxide, as a model for CT-cytochrome heme configuration, were also conducted. The ferric and the ferrous CT-cytochromes above pH 7.5 have similar, yet distinct, spectroscopic properties, absorption, CD, resonance Raman, and PMR spectra, typical of low-spin hexacoordinated hemes, but distinct from those of the unmodified protein. The ferric spectrum lacks the 695-nm band, and the reduced spectrum contains an additional inflection at about 400 nm, a feature also observed in the spectra of ferrous H8-DMSO systems. The CD, resonance Raman, and PMR spectra are typical of a cytochrome with a loosened heme crevice and altered coordination configuration. The Methionine-80 proton resonances are absent in the uupfield PMR spectra of both the CT-ferricytochromes. The ferrous spectra, on the other hand, contain all the Met-80 resonances, but with smaller upfield shifts than those of the native protein. Both CT-ferric cytochromes are less stable in the acid region and convert to high-spin forms with a two-step transition and with a distinct set of pK _a values. The overall conformation is nearly identical to that of the native protein, but it is less stable to thermal unfolding. All the factors differentiating the modified preparations from the unmodified protein are more pronunced in the case of F^II, with F^III being the closest to the unmodified form. The two functionally distinct CT-cytochromes are two conformational isomers; conformationally and heme configurationally, they are spectroscopically very similar, yet distinct. Both contain an altered heme iron coordination configuration. The sulfur of Met-80 is repalced by the oxygen of Met-80 sulfoxide of a different configuration, R or S. Both contain a loosened heme crevice and are conformationally less stable than the native protein, F^II CT-cytochrome c being the most deranged.     

Structure analysis and characterization of the cytochrome c-554 from thermophilic green sulfur photosynthetic bacterium Chlorobaculum tepidum

The cytochrome (Cyt) c-554 in thermophilic green photosynthetic bacterium Chlorobaculum tepidum serves as an intermediate electron carrier, transferring electrons to the membrane-bound Cyt c _z from various enzymes involved in the oxidations of sulfide, thiosulfate, and sulfite compounds. Spectroscopically, this protein exhibits an asymmetric α-absorption band for the reduced form and particularly large paramagnetic ¹H NMR shifts for the heme methyl groups with an unusual shift pattern in the oxidized form. The crystal structure of the Cyt c-554 has been determined at high resolution. The overall fold consists of four α-helices and is characterized by a remarkably long and flexible loop between the α3 and α4 helices. The axial ligand methionine has S-chirality at the sulfur atom with its C^εH₃ group pointing toward the heme pyrrole ring I. This configuration corresponds to an orientation of the lone-pair orbital of the sulfur atom directed at the pyrrole ring II and explains the lowest-field ¹H NMR shift arising from the 18¹ heme methyl protons. Differing from most other class I Cyts c, no hydrogen bond was formed between the methionine sulfur atom and polypeptide chain. Lack of this hydrogen bond may account for the observed large paramagnetic ¹H NMR shifts of the heme methyl protons. The surface-exposed heme pyrrole ring II edge is in a relatively hydrophobic environment surrounded by several electronically neutral residues. This portion is considered as an electron transfer gateway. The structure of the Cyt c-554 is compared with those of other Cyts c, and possible interactions of this protein with its electron transport partners are discussed.     

Crystal structure of HugZ, a novel heme oxygenase from Helicobacter pylori

The crystal structure of a heme oxygenase (HO) HugZ from Helicobacter pylori complexed with heme has been solved and refined at 1.8 Å resolution. HugZ is part of the iron acquisition mechanism of H. pylori, a major pathogen of human gastroenteric diseases. It is required for the adaptive colonization of H. pylori in hosts. Here, we report that HugZ is distinct from all other characterized HOs. It exists as a dimer in solution and in crystals, and the dimer adopts a split-barrel fold that is often found in FMN-binding proteins but has not been observed in hemoproteins. The heme is located at the intermonomer interface and is bound by both monomers. The heme iron is coordinated by the side chain of His²⁴⁵ and an azide molecule when it is present in crystallization conditions. Experiments show that Arg¹⁶⁶, which is involved in azide binding, is essential for HugZ enzymatic activity, whereas His²⁴⁵, surprisingly, is not, implying that HugZ has an enzymatic mechanism distinct from other HOs. The placement of the azide corroborates the observed γ-meso specificity for the heme degradation reaction, in contrast to most known HOs that have α-meso specificity. We demonstrate through sequence and structural comparisons that HugZ belongs to a new heme-binding protein family with a split-barrel fold. Members of this family are widespread in pathogenic bacteria and may play important roles in the iron acquisition of these bacteria.     

1H, 15N, and 13C resonance assignments of the 2/2 hemoglobin from the cyanobacterium Synechococcus sp. PCC 7002 in the ferric bis-histidine state

The hemoglobin from the cyanobacterium Synechococcus sp. PCC 7002 is a monomeric 123-residue Group I 2/2 hemoglobin. Here, we report ¹H, ¹⁵N, and ¹³C assignments for the ferric (low-spin, S = ½) protein with a b heme cofactor and after post-translational modification leading to a c-like heme.     

Identification of a Thiol/Disulfide Redox Switch in the Human BK Channel That Controls Its Affinity for Heme and CO

Heme is a required prosthetic group in many electron transfer proteins and redox enzymes. The human BK channel, which is a large-conductance Ca²⁺ and voltage-activated K⁺ channel, is involved in the hypoxic response in the carotid body. The BK channel has been shown to bind and undergo inhibition by heme and activation by CO. Furthermore, evidence suggests that human heme oxygenase-2 (HO2) acts as an oxygen sensor and CO donor that can form a protein complex with the BK channel. Here we describe a thiol/disulfide redox switch in the human BK channel and biochemical experiments of heme, CO, and HO2 binding to a 134-residue region within the cytoplasmic domain of the channel. This region, called the heme binding domain (HBD) forms a linker segment between two Ca²⁺-sensing domains (called RCK1 and RCK2) of the BK channel. The HBD includes a CXXCH motif in which histidine serves as the axial heme ligand and the two cysteine residues can form a reversible thiol/disulfide redox switch that regulates affinity of the HBD for heme. The reduced dithiol state binds heme (K_d = 210 nm) 14-fold more tightly than the oxidized disulfide state. Furthermore, the HBD is shown to tightly bind CO (K_d = 50 nm) with the Cys residues in the CXXCH motif regulating affinity of the HBD for CO. This HBD is also shown to interact with heme oxygenase-2. We propose that the thiol/disulfide switch in the HBD is a mechanism by which activity of the BK channel can respond quickly and reversibly to changes in the redox state of the cell, especially as it switches between hypoxic and normoxic conditions.     

Dissecting the Functional Role of Key Residues in Triheme Cytochrome PpcA: A Path to Rational Design of G. sulfurreducens Strains with Enhanced Electron Transfer Capabilities

PpcA is the most abundant member of a family of five triheme cytochromes c ₇ in the bacterium Geobacter sulfurreducens (Gs) and is the most likely carrier of electrons destined for outer surface during respiration on solid metal oxides, a process that requires extracellular electron transfer. This cytochrome has the highest content of lysine residues (24%) among the family, and it was suggested to be involved in e⁻/H⁺ energy transduction processes. In the present work, we investigated the functional role of lysine residues strategically located in the vicinity of each heme group. Each lysine was replaced by glutamine or glutamic acid to evaluate the effects of a neutral or negatively charged residue in each position. The results showed that replacing Lys⁹ (located near heme IV), Lys¹⁸ (near heme I) or Lys²² (between hemes I and III) has essentially no effect on the redox properties of the heme groups and are probably involved in redox partner recognition. On the other hand, Lys⁴³ (near heme IV), Lys⁵² (between hemes III and IV) and Lys⁶⁰ (near heme III) are crucial in the regulation of the functional mechanism of PpcA, namely in the selection of microstates that allow the protein to establish preferential e⁻/H⁺ transfer pathways. The results showed that the preferred e⁻/H⁺ transfer pathways are only established when heme III is the last heme to oxidize, a feature reinforced by a higher difference between its reduction potential and that of its predecessor in the order of oxidation. We also showed that K43 and K52 mutants keep the mechanistic features of PpcA by establishing preferential e⁻/H⁺ transfer pathways at lower reduction potential values than the wild-type protein, a property that can enable rational design of Gs strains with optimized extracellular electron transfer capabilities.     

ShmR Is Essential for Utilization of Heme as a Nutritional Iron Source in Sinorhizobium meliloti

The bacterium Sinorhizobium meliloti is able to use heme as a nutritional iron source. Here, we show that the iron-regulated shmR gene encodes an outer membrane protein required for growth on heme. Furthermore, an shmR mutant is resistant to the toxic heme analog gallium protoporphyrin. Thus, the receptor protein of the heme transport system has been identified in S. meliloti.     

Spectroscopic characterization of (57)Fe-enriched cytochrome c

Investigation of the heme iron dynamics in cytochrome c with Mössbauer spectroscopy and especially nuclear resonance vibrational spectroscopy requires the replacement of the natural abundant heme iron with the ⁵⁷Fe isotope. For demetallization, we use a safer and milder ferrous sulfate–hydrochloric acid method in addition to the harsher commonly used hydrofluoric acid-based procedure. The structural integrity of the ⁵⁷Fe-reconstituted protein in both oxidation states is confirmed from absorption spectra and a detailed analysis of the rich resonance Raman spectra. These results reinforce the application of metal-substituted heme c proteins as reliable models for the native proteins.     

Heme-iron utilization by Leptospira interrogans requires a heme oxygenase and a plastidic-type ferredoxin-NADP reductase

Background

Heme oxygenase catalyzes the conversion of heme to iron, carbon monoxide and biliverdin employing oxygen and reducing equivalents. This enzyme is essential for heme-iron utilization and contributes to virulence in Leptospira interrogans.

Methods

A phylogenetic analysis was performed using heme oxygenases sequences from different organisms including saprophytic and pathogenic Leptospira species. L. interrogans heme oxygenase (LepHO) was cloned, overexpressed and purified. The structural and enzymatic properties of LepHO were analyzed by UV–vis spectrophotometry and ¹H NMR. Heme-degrading activity, ferrous iron release and biliverdin production were studied with different redox partners.

Results

A plastidic type, high efficiently ferredoxin-NADP⁺ reductase (LepFNR) provides the electrons for heme turnover by heme oxygenase in L. interrogans. This catalytic reaction does not require a ferredoxin. Moreover, LepFNR drives the heme degradation to completeness producing free iron and α-biliverdin as the final products. The phylogenetic divergence between heme oxygenases from saprophytic and pathogenic species supports the functional role of this enzyme in L. interrogans pathogenesis.

Conclusions

Heme-iron scavenging by LepHO in L. interrogans requires only LepFNR as redox partner. Thus, we report a new substrate of ferredoxin-NADP⁺ reductases different to ferredoxin and flavodoxin, the only recognized protein substrates of this flavoenzyme to date. The results presented here uncover a fundamental step of heme degradation in L. interrogans.

General significance

Our findings contribute to understand the heme-iron utilization pathway in Leptospira. Since iron is required for pathogen survival and infectivity, heme degradation pathway may be relevant for therapeutic applications.     

Correlated behavior of the EPR signal of cytochrome b-559 heme Fe(III) ligated by OH− and the multiline signal of the Mn cluster in PS-II membrane fragments

EPR signals of Cyt b-559 heme Fe(III) ligated by OH⁻ and the multiline signal of the Mn cluster in PS-II membrane fragments have been investigated. In 2,3-dicyano-5,6-dichloro-p-benzoquinone-oxidized PS-II membrane fragments the light-induced decrease of the EPR signal of the heme Fe(III)-OH⁻ is accompanied by the appearance of the EPR multiline signal of the Mn cluster. Addition of F⁻ ions, which act as a stronger ligand for heme Fe(III) than OH⁻, decreases to the same extent the dark- and light-induced signal of the heme Fe(III)-OH⁻ and the light-induced multiline signal of the Mn cluster. These results are discussed in terms of the light-induced formation of a bound OH′ radical shared between the Cyt b-559 heme Fe and the Mn cluster as a first step of water oxidation.     

Staphylococcus aureus FepA and FepB Proteins Drive Heme Iron Utilization in Escherichia coli

EfeUOB-like tripartite systems are widespread in bacteria and in many cases they are encoded by genes organized into iron-regulated operons. They consist of: EfeU, a protein similar to the yeast iron permease Ftrp1; EfeO, an extracytoplasmic protein of unknown function and EfeB, also an extracytoplasmic protein with heme peroxidase activity, belonging to the DyP family. Many bacterial EfeUOB systems have been implicated in iron uptake, but a prefential iron source remains undetermined. Nevertheless, in the case of Escherichia coli, the EfeUOB system has been shown to recognize heme and to allow extracytoplasmic heme iron extraction via a deferrochelation reaction. Given the high level of sequence conservations between EfeUOB orthologs, we hypothesized that heme might be the physiological iron substrate for the other orthologous systems. To test this hypothesis, we undertook characterization of the Staphylococcus aureus FepABC system. Results presented here indicate: i) that the S. aureus FepB protein binds both heme and PPIX with high affinity, like EfeB, the E. coli ortholog; ii) that it has low peroxidase activity, comparable to that of EfeB; iii) that both FepA and FepB drive heme iron utilization, and both are required for this activity and iv) that the E. coli FepA ortholog (EfeO) cannot replace FepA in FepB-driven iron release from heme indicating protein specificity in these activities. Our results show that the function in heme iron extraction is conserved in the two orthologous systems.     

Studies on erythrocruorin. IV. Circular dichroism of earthworm erythrocruorin.

Circular dichroism spectra of Lumbricus erythrocruorin in the absence and in the presence of heme ligands have been analyzed under a variety of experimental conditions in view of the peculiarities in ligand binding displayed by this high molecular weight heme protein (M_r = 3 × 10⁶).The undisaociated molecule exists in a “metastable” form with high cooperativity in oxygen binding, which can be converted into a stable form with low co-operativity either by changes in pH or temperature; circular dichroism spectra of oxyerythrocruorin in the Soret region give direct evidence of a local alteration in the heme environment under the conditions which affect co-operativity in oxygen binding of the undissociated molecule. Similar, although more pronounced changes in the same spectral region are observed in the dissociated molecule of M_r = 270,000, which displays low co-operativity in oxygen binding.Deoxygenation is accompanied by an inversion in the double Soret-Cotton effect, which indicates a substantial rearrangement in the heme environment upon removal of the ligand.The double peak in the Soret region found in all erythrocruorin derivatives can be taken as an indication of a distinctive distribution of the aromatic side-chains interacting with the heme chromophore.     

Resonance Raman spectroscopic studies of hydroperoxo derivatives of cobalt-substituted myoglobin

Recent progress in generating and stabilizing reactive heme protein enzymatic intermediates by cryoradiolytic reduction has prompted application of a range of spectroscopic approaches to effectively interrogate these species. The impressive potential of resonance Raman spectroscopy for characterizing such samples has been recently demonstrated in a number of studies of peroxo- and hydroperoxo-intermediates. While it is anticipated that this approach can be productively applied to the wide range of heme proteins whose reaction cycles naturally involve these peroxo- and hydroperoxo-intermediates, one limitation that sometimes arises is the lack of enhancement of the key intraligand ν(O-O) stretching mode in the native systems. The present work was undertaken to explore the utility of cobalt substitution to enhance both the ν(Co-O) and ν(O-O) modes of the CoOOH fragments of hydroperoxo forms of heme proteins bearing a trans-axial histidine linkage. Thus, having recently completed RR studies of hydroperoxo myoglobin, attention is now turned to its cobalt-substituted analogue. Spectra are acquired for samples prepared with ¹⁶O₂ and ¹⁸O₂ to reveal the ν(M-O) and ν(O-O) modes, the latter indeed being observed only for the cobalt-substituted proteins. In addition, spectra of samples prepared in deuterated solvents were also acquired, providing definitive evidence for the presence of the hydroperoxo-species.     

State of heme in heme c systems: Cytochrome c and heme c models

Polarized resonance Raman spectra of horse heart ferricytochrome c as a function of pH in the range 1.0–12, in the presence of the extrinsic ligands imidazole, cyanide, and azide, and in 4 M urea, are reported, as are resonance Raman spectra of heme undecapeptide in the presence of imidazole, pH 6.8 and pH 2.0, and with cyanide at pH 6.8. The range of investigation is 140–1700 cm^?1, using the 5145-, 4880-, and 4579-Å excitations. The spectra have been analyzed in terms of complexity, sensitivity, and the conformation-heme energetics of the systems. The state of heme in various forms is analyzed with regard to heme energetics, core size, nature of planarity, and coordination configuration. All low-spin forms of heme c systems, cytochrome c, and heme models are concluded to be hexacoordinated, in-plane heme iron systems. The effect of the location of the heme in the protein environment is found to be a slight expansion of the porphyrin core, ～0.01 Å, while the covalent linkage of heme to protein and a mixed nature of axial coordination configuration seem to have little effect on the energetics of the heme group. Complex formation with extrinsic ligand, imidazole, cyanide, or azide, results in a slight contraction of the heme core. The formation of cytochrome c form IV, the alkaline form, is shown to follow a process with apK _a of about 8.4, and similarly, acidic form II is created following the prior formation of an intermediate form with apK _a of about 3.6. The precursor to form IV is interpreted as containing perturbation of the pyrrol rings, whereas the precursor to the acidic form seems to reflect alteration of the energetics of the C_αC_{m α} structures of the heme group. The acidic form of heme undecapeptide is a hexacoordinated high-spin heme with an estimated displacement of 0.25 Å from the heme plane. The pH 2 form of cytochrome c is also a hexacoordinated high-spin form with two weak axial ligands, but iron is in the plane of the porphyrin ring.