Comparison of the evolution of heme binding and NAD binding proteins.

Of the 85 three-dimensionally characterized residues of cytochrome b5, 51 are structurally and topologically equivalent to the globin fold. When these proteins have been superimposed, the heme irons are found to be less than 1.4 A separated and the heme normals are inclined by less than 9.5 degrees. Comparison of minimum base changes per codon between heme binding and NAD binding proteins are of the same order.     

Structural comparisons of heme binding proteins.

Of the 82 three dimensionally characterized residues of cytochrome c551, 49 are found to be structurally and topologically equivalent to the globin fold and 41 are equivalent to the cytochrome b5 fold, with a respective root mean square separation of 3.5 and 4.9 A between equivalenced Calpha atoms. The common fold represents a central heme binding core, corresponding to the middle exon of certain globin genes. After superposition of the protein folds, the heme irons are found to be separated by 5.4 and 1.6 A, while their heme normals are inclined by 6 degrees and 32 degrees, respectively. Furthermore, the heme "face", determined by the asymmetric attachment of the vinyl and propionyl side chains, is directed similarly in all three heme proteins. The heme itself is rotated by 72 degrees and 116 degrees about its normal, respectively. The minimum base change per codon for the three pairwise comparisons corresponds to the expected value of random sequence comparisons. While all three heme proteins may have diverged from a common ancestor, their similarity may have arisen from the requirements of heme binding or the utilization of a particularly stable fold. Known structures within commonly accepted divergent families were superimposed in order to discriminate better between convergence and divergence. Minimum base changes per codon, number of deletions and insertions, percentage of equivalenced residues, precision of heme superposition, and root mean square separation of equivalenced Calpha atoms were tested as measures of evolutionary relationships.     

Non-functional conserved residues in globins and their possible role as a folding nucleus.

Structure-based sequence alignment of 728 sequences of different globin subfamilies shows that in each subfamily there are two clusters of consensually conserved residues. The first is the well-known "functional" cluster which includes six heme-binding conserved residues (Phe CD1, His F8; aliphatic E11, FG5; hydrophobic F4, G5) and seven other conserved residues (Pro C2; aliphatic H19; hydrophobic B10, B13, B14, CD4, E4) that do not bind the heme but belong to its immediate neighborhood. The second cluster revealed here (aliphatic A8, G16, G12; aromatic A12; hydrophobic H8 and possibly H12) is distant from the heme. It is entirely non-polar and includes one turn (i, i+4 positions) from each of helices A, G, and H. It is known that A, G, and H helices formed at the earliest stage of apomyoglobin folding remain relatively stable in the equilibrium molten globule state, and are likely to be tightly packed with each other in this state. We have shown the existence of two similar conserved clusters in c -type cytochromes, heme-binding and distal from the heme. The second cluster in c -cytochromes includes one turn from each of the N and C-terminal alpha-helices. These N and C-terminal helices in cytochrome c are formed at the earliest stage of protein folding, remain relatively stable in the molten globule state, and are tightly packed with each other in this state, similar to the observed behavior of the globins. At least these two large protein families (c -type cytochromes and globins) have a close similarity in the existence and mutual positions of non-functional conserved residues. We assume that non-functional conserved residues are requisite for the fast and correct folding of both of these protein families into their stable 3D structures.     

The structure of cytochrome b562 from Escherichia coli at 2.5 A resolution.

The structure of cytochrome b562 from Escherichia coli has been determined at 2.5 A resolution by x-ray diffraction methods. Protein phases were computed by the single isomorphous replacement method with anomalous scattering measurements from the native and uranyl acetate-substituted crystals. The electron density was averaged about the noncrystallographic 2-fold axis relating 2 molecules in the triclinic unit cell. The protein consists of four nearly parallel alpha helices and represents a new class of cytochrome structure. The heme group is inserted between the helices near one end of the molecule with one heme face partially exposed to solvent. The two heme ligands are histidine and methionine. The 2 phenylalanines are packed internally near the heme group, and the 2 tyrosines are on the surface, also near the heme group. The folding of the protein resembles that of hemerythrin and tobacco mosaic virus protein and shows a different topology from that of cytochrome b5, cytochrome c, or the globins.     

Invariant features of globin primary structure and coding of their secondary structure

Invariant features of the primary structure of 67 globins are analysed. These features may be responsible for the formation of the secondary structure of these proteins at the first stage of self-organization (in the unfolded chain). It is shown that in primary structures of globins there are 11 sites or regions of one to four residues in which at least one of the residues Asn, Asp, His, Pro, Ser or Thr is located in every globin (haem-linking His residues are excluded from these sites). An unambiguous correlation exists between the position of these regions and the secondary structure of globins: all these regions (except one) are located near the ends of helices in globins whose three-dimensional structure is known and the ends of all helices (except for the helix F) are coded by such regions. A decrease in the set of residues listed above leads to a sharp drop in the number of regions invariantly occupied by the residues, while an addition of residues such as Tyr and Gly to this set does not eventually increase the number of invariant regions. Five residues (Asn, Asp, His, Ser and Thr) of the six that code the ends of helical regions have polar side groups with a small number of degrees of freedom capable of forming hydrogen bonds with atoms of the backbone with a relatively small loss of entropy. One residue (Pro) has no NH-group and, therefore, has less chance of participating in the formation of hydrogen bonds between atoms of the backbone. This corroborates the hypothesis that competition between hydrogen bonds of short polar side groups and hydrogen bonds in the backbone is essential for the formation of the secondary structure in unfolded protein chains. Amino acid replacements in hydrophobic cores of the 67 globins are considered in the Appendix.     

Solution 1H NMR study of the heme cavity and folding topology of the abbreviated chain 118-residue globin from the cyanobacterium Nostoc commune

The globin from the cyanobacterium Nostoc commune, abbreviated GlbN, which appears to serve as a part of a terminal oxidase rather than as a respiratory pigment, displays relatively normal O2 binding properties, despite the highly abbreviated polypeptide chain, (118 residues) relative to more conventional globins [Thorsteinsson, M. V. , Bevan, D. R., Potts, M., Dou, Y., Eich, R. F., Hargrove, M. S., Gibson, Q. H., and Olson, J. S. (1999) Biochemistry 38, 2117-2126]. The nature of the heme cavity and the general folding topology of this cyanoglobin were investigated by solution 1H NMR to establish the extent to which, and the manner in which, this compact globin adheres to the standard globin fold. This represents by far the smallest globin subjected to structural analysis. The paramagnetic cyanomet derivative was selected because its characteristically large magnetic anisotropy imparts significant dipolar shifts which both improve resolution to greatly facilitate assignments and serve as indicators of the folding topology of the globin. Identification of the axial His 70 and highly conserved Phe 35 (CD1) determined the absolute orientation of the heme and proximal His. Sequential assignments of four helical and one loop segments, which exhibit dipolar contacts to the heme and among each other, confirm the presence of well-conserved F, G, and H helices and the FG corner. The majority of the abbreviation of the chain relative to the more conventional length globins is accommodated in the A-D helices, of which the last is completely missing. The distal residue which provides a H-bond to bound ligand is identified as Gln 43, but the expected helical position E7 could not be confirmed. His 46, placed at position E10, is found to adopt alternate orientations into, and out of, the heme cavity depending on protonation state, suggesting the presence of a Bohr effect at low pH. It is shown that the dipolar shifts exhibited by backbone protons for the assigned residues conform well to those observed for other cyanomet globins and further support a conserved Mb fold. Perturbed medium-range dipolar contacts and the pH-independent backbone proton lability of the F helix are interpreted in terms of a holoprotein which is less stable than a conventional length globin.     

Sequence analysis of cytochrome bd oxidase suggests a revised topology for subunit I

Numerous sequences of the cytochrome bd quinol oxidase (cytochrome bd) have recently become available for analysis. The analysis has revealed a small number of conserved residues, a new topology for subunit I and a phylogenetic tree involving extensive horizontal gene transfer. There are 20 conserved residues in subunit I and two in subunit II. Algorithms utilizing multiple sequence alignments predicted a revised topology for cytochrome bd, adding two transmembrane helices to subunit I to the seven that were previously indicated by the analysis of the sequence of the oxidase from E. coli. This revised topology has the effect of relocating the N-terminus and C-terminus to the periplasmic and cytoplasmic sides of the membrane, respectively. The new topology repositions I-H19, the putative ligand for heme b595, close to the periplasmic edge of the membrane, which suggests that the heme b595/heme d active site of the oxidase is located near the outer (periplasmic) surface of the membrane. The most highly conserved region of the sequence of subunit I contains the sequence GRQPW and is located in a predicted periplasmic loop connecting the eighth and ninth transmembrane helices. The potential importance of this region of the protein was previously unsuspected, and it may participate in the binding of either quinol or heme d. There are two very highly conserved glutamates in subunit I, E99 and E107, within the third transmembrane helix (E. coli cytochrome bd-I numbering). It is speculated that these glutamates may be part of a proton channel leading from the cytoplasmic side of the membrane to the heme d oxygen-reactive site, now placed near the periplasmic surface. The revised topology and newly revealed conserved residues provide a clear basis for further experimental tests of these hypotheses. Phylogenetic analysis of the new sequences of cytochrome bd reveals considerable deviation from the 16sRNA tree, suggesting that a large amount of horizontal gene transfer has occurred in the evolution of cytochrome bd.     

How cytochromes with different folds control heme redox potentials

The electrochemical midpoint potentials (E(m)'s) of 13 cytochromes, in globin (c, c(2), c(551), c(553)), four-helix bundle (c', b(562)), alpha beta roll (b(5)), and beta sandwich (f) motifs, with E(m)'s spanning 450 mV were calculated with multiconformation continuum electrostatics (MCCE). MCCE calculates changes in oxidation free energy when a heme-axial ligand complex is moved from water into protein. Calculated and experimental E(m)'s are in good agreement for cytochromes with His-Met and bis-His ligated hemes, where microperoxidases provide reference E(m)'s. In all cytochromes, E(m)'s are raised by 130-260 mV relative to solvated hemes by the loss of reaction field (solvation) energy. However, there is no correlation between E(m) and heme surface exposure. Backbone amide dipoles in loops or helix termini near the axial ligands raise E(m)'s, but amides in helix bundles contribute little. Heme propionates lower E(m)'s. If the propionic acids are partially protonated in the reduced cytochrome, protons are released on heme oxidation, contributing to the pH dependence of the E(m). In all cytochromes studied except b(5)'s and low potential globins, buried side chains raise E(m)'s. MCCE samples ionizable group protonation states, heme redox states, and side chain rotamers simultaneously. Globins show the largest structural changes on heme oxidation and four-helix bundles the least. Given the calculated protein-induced E(m) shift and measured cytochrome E(m) the five-coordinate, His heme in c' is predicted to have a solution E(m) between that of isolated bis-His and His-Met hemes, while the reference E(m) for His-Ntr ligands in cytochrome f should be near that of His-Met hemes.     

Structural basis of human cytoglobin for ligand binding

Cytoglobin (Cgb), a newly discovered member of the vertebrate globin family, binds O(2) reversibly via its heme, as is the case for other mammalian globins (hemoglobin (Hb), myoglobin (Mb) and neuroglobin (Ngb)). While Cgb is expressed in various tissues, its physiological role is not clearly understood. Here, the X-ray crystal structure of wild type human Cgb in the ferric state at 2.4A resolution is reported. In the crystal structure, ferric Cgb is dimerized through two intermolecular disulfide bonds between Cys38(B2) and Cys83(E9), and the dimerization interface is similar to that of lamprey Hb and Ngb. The overall backbone structure of the Cgb monomer exhibits a traditional globin fold with a three-over-three alpha-helical sandwich, in which the arrangement of helices is basically the same among all globins studied to date. A detailed comparison reveals that the backbone structure of the CD corner to D helix region, the N terminus of the E-helix and the F-helix of Cgb resembles more closely those of pentacoordinated globins (Mb, lamprey Hb), rather than hexacoordinated globins (Ngb, rice Hb). However, the His81(E7) imidazole group coordinates directly to the heme iron as a sixth axial ligand to form a hexcoordinated heme, like Ngb and rice Hb. The position and orientation of the highly conserved residues in the heme pocket (Phe(CD1), Val(E11), distal His(E7) and proximal His(F8)) are similar to those of other globin proteins. Two alternative conformations of the Arg84(E10) guanidium group were observed, suggesting that it participates in ligand binding to Cgb, as is the case for Arg(E10) of Aplysia Mb and Lys(E10) of Ngb. The structural diversities and similarities among globin proteins are discussed with relevance to molecular evolutionary relationships.     

Variation of folded polypeptide surface area with probe size

Three types of polypeptide surface area (contact, accessible, and molecular) have been studied as a function of the radius of a probe sphere used to map the surface. The surfaces are: (1) three alpha-helices, the H-helix of myoglobin, the E-helix of leghemoglobin, and an artificial polyalanine helix, each with 26 residues; (2) two globins, myoglobin and leghemoglobin, each with 153 residues; and (3) a two-center model system for which the three types of surface area have been calculated analytically. The two globin helices have almost identical surface areas as a function of probe size as do the two globins. The polyalanine helix surface area is smaller but similar in shape to the globin helix areas. All three helix contact areas tend to the same limit as the probe size increases, and the globin contact areas behave similarly. Fractal dimensions were calculated for the helix and globin contact and molecular surfaces. All fractal dimensions showed strong dependence on probe size. The contact fractal dimension peaks at larger values for both the helices and globins. Most residues do not make contact with large probes (15 A).     

pH dependence of heme electrochemistry in cytochromes investigated by multiconformation continuum electrostatic calculations

Cytochromes belong to a diverse family of heme-containing redox proteins that function as intermediaries in electron transfer chains. They can be soluble, extrinsic, or intrinsic membrane proteins, and are found in different structural motifs (globin, 4-helix bundles, alpha beta roll, beta sandwich). Measured electrochemical midpoint potentials vary over a wide range even though the basic redox reaction at the heme is the same for all cytochromes. The perturbation of the heme electrochemistry is induced by the protein structure. Also, the pH dependence varies since it depends on the strength of interaction between the heme and surrounding residues as well as the ionization states of these groups. Multiconformation continuum electrostatics (MCCE) has been used to investigate the pH dependence of heme electrochemistry in cytochromes with different folds. Often propionates are the primary contributors for pH dependence especially if they are partially protonated in the reduced heme as it is shown for globin cytochrome c551 P. aeruginosa and cytochrome b5 R. norvegicus (alpha beta roll). However, if the propionates are already fully ionized at a certain pH they do not contribute to the pH dependence even if they have big interaction with the heme. At pH 7 there is no propionate contribution for cytochrome f C. reinhardtii (beta sandwich) and the 4-helix bundle c' R. palustris. Other residues can also change their ionization significantly during heme oxidation and therefore be involved in proton release and pH dependence. These residues have been identified for different cytochrome types.     

Common dynamics of globin family proteins

The recently discovered new members of the globin family, neurogobin and cytoglobin, are the object of sustained structural and functional studies aimed at understanding their physiological role and elucidating the impact of their bis-his heme hexacoordination. However, no studies have yet considered the dynamics of this protein family, an essential link between structure and function. In this communication, we present normal mode analysis results for neuroglobin, cytoglobin, hemoglobin and myoglobin to provide exploratory insights into globin characteristic motions. Our results show a clear correlation in the protein dynamics of this family. All four globins exhibit a high degree of correlated displacements involving residues in the C, E and F helices and link regions. They suggest that these motions play an important role in the reversible oxygen binding function of these proteins. Further, our results may help rationalize some functional features of the 6c-globins in that they alone exhibit correlated displacements of the G-helix region.     

Models for the complexes formed between cytochrome b5 and the subunits of methemoglobin

Computer graphics-generated models for the electron transfer complexes formed between cytochrome b5 and the subunits of methemoglobin are proposed. For both complexes, the orientation allowing optimal hydrogen bonding involves interaction between negatively charged residues on cytochrome b5 and positively charged residues on methemoglobin. In each complex, the heme groups of the interacting species are coplanar with the edges of the heme groups separated by 7-8 A and with the iron atoms 16 A apart. For the alpha-chain X cytochrome b5 complex, alpha-chain residues 56 (Lys), 60 (Lys), and 90 (Lys) interact with cytochrome b5 residues 44 (Glu), 43 (Glu), and 60 (Asp) respectively. A fourth hydrogen bond involves alpha-61 (Lys) bridging between a heme propionate from cytochrome b5 and a heme propionate from the alpha-chain. The contacts present in the beta-chain X cytochrome b5 complex involve hydrogen-bonding between beta-chain lysyl residues 59, 61, 65, and 95, and cytochrome b5 residues 48 (Glu), 44 (Glu), 43 (Glu), and 60 (Asp) respectively. An additional hydrogen bond can be formed by bridging of the epsilon-amino group of beta-66 (Lys) between a heme propionate from cytochrome b5 and a beta-chain heme propionate. In each complex, two nonionic interactions, one on each side of the heme groups, are also suggested. These interactions appear to effectively exclude external water molecules from the center of the protein-protein interaction domain. Comparison of the proposed binding loci for cytochrome b5 on the methemoglobin subunits with those proposed on cytochrome c reveals considerable structural homology between the cytochrome b5 binding sites.     

Crystal structure of recombinant trypsin-solubilized fragment of cytochrome b(5) and the structural comparison with Val61His mutant

The crystal structure of the recombinant trypsin-solubilized fragment of the microsomal cytochrome b(5) from bovine liver has been determined at 1.9 A resolution and compared with the reported crystal structure of the lipase-solubilized fragment of the membrane protein cytochrome b(5). The two structures are similar to each other. However, some detailed structural differences are observed: the conformation of the segment Asn16-Ser20 is quite different, some helices around the heme and some segments between the helices are shifted slightly, the heme is rotated about the normal of the mean plane of heme, one of the propionates of the heme exhibits a different conformation. The average coordination distances between the iron and the two nitrogen atoms of the imidazole ligands are the same in the two structures. Most of the structural differences can be attributed to the different intermolecular interactions which result from the crystal packing. The wild-type protein structure is also compared with its Val61His mutant, showing that the heme binding and the main chain conformations are basically identical with each other except for the local area of the mutation site. However, when Val61 is mutated to histidine, the large side chain of His61 is forced to point away from the heme pocket toward the solvent region, disturbing the micro-environment of the heme pocket and influencing the stability and the redox potential of the protein.     

Helix movements and the reconstruction of the haem pocket during the evolution of the cytochrome c family

Analysis of cytochromes c (tuna), c2 (Rhodospirillum rubrum), c550 (Paracoccus denitrificans) and c551 (Pseudomonas aeruginosa) shows that they contain 48 residues identifiable as homologous from superposition of the structures. The other 34 to 64 residues are in loops that vary greatly in sequence, length and conformation, or in alpha-helices that are found in only some of the structures. Of the 48 homologous residues, 17 are in three segments which pack onto the haem faces. In all four structures, these segments have the same conformations, and the same locations relative to the haem. The other 31 residues are in three alpha-helices which are in contact with each other. These form the back and one side of the haem pocket. In cytochrome c551 the positions of the three alpha-helices have shifted and rotated, in comparison with cytochromes c and c2, by up to 5 A and 25 degrees relative to the haem. These shifts, facilitated by mutations at the helix-helix interfaces, are related to the reconstruction of the propionic acid side of the haem pocket described by Almassy & Dickerson (1978). Together these effects produce alternative structures for the haem pocket. This mechanism of adaptation to mutation contrasts with that observed in the globins. In the globins, mutations also produce changes in helix interfaces and shifts of packed helices, but in the globins these shifts are coupled to conserve the structure of the haem pocket.     

A modeling study of the interaction and electron transfer between cytochrome b 5 and some oxidized haemoglobins

Using Brownian motion simulations we have studied the formation of docked complexes of reduced cytochrome b5 and oxidized haemoglobin. Our results indicate that the presence of molecular electrostatic fields has a significant role to play in the formation of these complexes. In contrast to previous modeling studies on this system, we clearly identify electron transfer within an ensemble of similarly docked complexes rather than the formation of a single complex. Docking involves a number of acidic residues surrounding the exposed haem edge of cytochrome b5 and a set of basic residues surrounding the exposed haem edge of the globins. Although amino acids from the partner globin proteins are involved to a small extent in the binding of some of the complexes, the reactivity of any particular globin is essentially independent of the nature of its partner globin chain within the haemoglobin molecule. Comparison of results from adult and embryonic haemoglobins indicates a significant difference in complex formation. Application of electron tunneling analysis to the complexes allows us to predict the rates of electron transfer within each ensemble of complexes. These data provide a theoretical insight into the important process of re-reduction of oxidized haemoglobins as well as explaining the general inability to produce crystalline forms of many docked electron transfer complexes.     

The N-terminal region of the heme-regulated eIF2alpha kinase is an autonomous heme binding domain.

The N-terminal domain (NTD) of the heme-regulated eukaryotic initiation factor (eIF)2alpha kinase (HRI) was aligned to sequences in the NCBI data base using ENTREZ and a PAM250 matrix. Significant similarity was found between amino acids 11-118 in the NTD of rabbit HRI and amino acids 16-120 in mammalian alpha-globins. Several conserved amino acid residues present in globins are conserved in the NTD of HRI. His83 of HRI was predicted to be equivalent to the proximal heme ligand (HisF8) that is conserved in all globins. Molecular modeling of the NTD indicated that its amino acid sequence was compatible with the globin fold. Recombinant NTD (residues 1-159) was expressed in Escherichia coli. Spectral analysis of affinity purified recombinant NTD indicated that the NTD contained stably bound hemin. Mutational analysis indicated that His83 played a critical structural role in the stable binding of heme to the NTD, and was required to stabilize full length HRI synthesized de novo in the rabbit reticulocyte lysate. These results indicate that the NTD of HRI is an autonomous heme-binding domain, with His83 possibly serving as the proximal heme binding ligand.     

Amino acid sequence of the monomer subunit of the giant extracellular hemoglobin of the aquatic oligochaete, Tubifex tubifex

The extracellular hemoglobin of the aquatic oligochaete Tubifex tubifex consists of four subunits: a monomer of 16.5 kDa, a disulfide-bonded trimer of about 50 kDa and at least two subunits of about 30 kDa. The complete amino acid sequence of the monomeric subunit was determined: it consists of 141 amino acid residues and has a molecular mass of 16,286 Da including a heme group. 39 residues (28%) were found to be identical with those in the corresponding positions in the monomeric globin chains from Lumbricus terrestris, Pheretima sieboldi, and Tylorrhynchus heterochaetus. Tubifex and Lumbricus are most similar, with 75 amino acid identities (53%). There are eight invariant residues amongst these monomeric globins and the intracellular monomeric globin of Glycera and the human beta-globin. The monomeric globin from Tubifex aligns best with those of group A, globins which have a Cys in their second position and an invariant Lys-Val-Lys at positions 9-11 [Gotoh et al. (1987) Biochem. J. 241, 441-445]. The two cysteine residues, at positions 2 and 131, appear to be disulfide-bonded.     

Solution structure and dynamics of the functional domain of Paracoccus denitrificans cytochrome c(552) in both redox states

A soluble and fully functional 10.5 kDa fragment of the 18.2 kDa membrane-bound cytochrome c(552) from Paracoccus denitrificans has been heterologously expressed and (13)C/(15)N-labeled to study the structural features of this protein in both redox states. Well-resolved solution structures of both the reduced and oxidized states have been determined using high-resolution heteronuclear NMR. The overall protein topology consists of two long terminal helices and three shorter helices surrounding the heme moiety. No significant redox-induced structural differences have been observed. (15)N relaxation rates and heteronuclear NOE values were determined at 500 and 600 MHz. Several residues located around the heme moiety display increased backbone mobility in both oxidation states, while helices I, III, and V as well as the two concatenated beta-turns between Leu30 and Arg36 apparently form a less flexible domain within the protein structure. Major redox-state-dependent differences of the internal backbone mobility on the picosecond-nanosecond time scale were not evident. Hydrogen exchange experiments demonstrated that the slow-exchanging amide proton resonances mainly belong to the helices and beta-turns, corresponding to the regions with high order parameters in the dynamics data. Despite this correlation, the backbone amide protons of the oxidized cytochrome c(552) exchange considerably faster with the solvent compared to the reduced protein. Using both differential scanning calorimetry as well as temperature-dependent NMR spectroscopy, a significant difference in the thermostabilities of the two redox states has been observed, with transition temperatures of 349.9 K (76.8 degrees C) for reduced and 307.5 K (34.4 degrees C) for oxidized cytochrome c(552). These results suggest a clearly distinct backbone stability between the two oxidation states.     

Structure of beef liver catalase

The three-dimensional structure of beef liver catalase has been determined to 2.5 å resolution by a combination of isomorphous and molecular replacement techniques. Heavy-atom positions were found using vector search and difference Fourier methods. The tetrameric catalase molecule has 222 symmetry with one of its dyads coincident with a crystallographic 2-fold axis. The known polypeptide sequence has been unambiguously fitted to the electron density map. The heme is well buried in a hydrophobic pocket, 20 Å below the surface of the molecule, and accessible through a hydrophobic channel. Residues that line the heme pocket belong to two different subunits. Tyr357 is the proximal heme ligand and the catalytically important residues on the distal side are residues His74 and Asnl47. The tertiary structure consists of four domains: an extended non-globular amino-terminal arm, which stabilizes the quaternary structure; an anti-parallel, eight-stranded β-barrel providing the residues on the distal side of the heme; a rather random “wrapping domain” around the subunit exterior including the proximal heme ligand; and a final λ-helical structure resembling the E, F, G and H helices of the globins.