Solution NMR study of environmental effects on substrate seating in human heme oxygenase: influence of polypeptide truncation, substrate modification and axial ligand

Solution proton NMR has been used here to show that, as either the high-spin ferric, protohemin (PH) substrate complex at neutral pH, or the low-spin ferric, cyanide-inhibited PH substrate complex, the active site electronic and molecular structure of the 233- and 265-residue recombinant constructs of human heme oxygenase-1, hHO, are essentially indistinguishable. It is shown, moreover, that the equilibrium PH orientational isomerism about the alpha,gamma-meso axis is 1:1 in the water-ligated, resting-state complex, but changes to a 4:1 equilibrium ratio as the cyanide-inhibited complex, with the minor species in solution corresponding to the only one found in crystals. The introduction of significant PH orientational preference in the cyanide over the aquo complex is rationalized by the crystallographic observation for the same H2O and CN ligated complexes of rat heme oxygenase (rHO), where the steric tilt of the Fe-CN unit resulted in a approximately 1 A transition of PH into the hydrophobic interior, and stronger interaction of the vinyls with the HO matrix [M. Sugishima, H. Sakamoto, M. Noguchi, K. Fukugama, Biochemistry 42 (2003) 9898-9905]. 1H NMR spectra of the cyanide-inhibited PH complex are the most used, and most useful, for determining the distribution of orientational isomerism for PH in complexes of HO. Hence, it is imperative that the time-course of the spectra after sample preparation be considered in order to reach conclusions that relate isomeric seating of the heme with variable isomeric biliverdin products. The natural orientational isomerism of PH leads to spectral congestion that has prompted the use of a synthetic, twofold symmetric substrate, 2,4-dimethyldeuterohemin, DMDH. While the hyperfine shift pattern for non-ligated residues are very similar and are consistent with largely conserved molecular structure with the alternate substrates, the steric tilt of the Fe-CN vector towards the protein interior, as determined by the orientation of the major magnetic axes, is 2 degrees smaller for DMDH than PH, and is rationalized by the substrate translating even further into the hydrophobic interior in the cyanide complex when the bulky vinyl groups are replaced by methyl groups.     

Assignment of resonances in the 1H nuclear magnetic resonance spectrum of the carbon monoxide complex of human hemoglobin alpha-chains

Assignments are reported for a substantial number of heme and amino acid proton resonances in the 1H nuclear magnetic resonance spectrum of the carbon monoxide complex of isolated hemoglobin alpha-chains. These resonances provide information on the solution conformation of the protein, particularly in the vicinity of the heme. The heme pocket structure is generally similar to that of carbonmonoxymyoglobin; several conserved residues adopt virtually identical positions relative to the heme in the two proteins. The largest conformational differences involve residues surrounding the ligand-binding site, notably Val62 (E11) and His58 (E7). The chemical shifts of the proximal His87 (F8) resonances are very similar in spectra of the two proteins, indicating a highly conserved coordination geometry and similar hydrogen bonding to the backbone carbonyl of Leu83 (F4).     

NMR studies of the heme pocket conformations of monomeric hemoglobins from Glycera dibranchiata. Implications for ligand binding

Two-dimensional 1H-NMR methods have been used to assign side-chain resonances for the tryptophan residues and for several amino acids located in the heme pockets of the carbon monoxide complexes of the major monomeric hemoglobins from Glycera dibranchiata. The NMR spectra reveal a high degree of conservation of the heme pocket structure in the different hemoglobins. However some conformational differences are evident and residues at positions B10 and G8 on the distal side of the heme pocket are not conserved. From the present NMR studies it appears that the monomeric G. dibranchiata hemoglobin examined by X-ray crystallography [Padlan, E. A. & Love, W. (1974) J. Biol. Chem. 249, 4067-4078] corresponds to HbC. Except that the orientation of the heme in solution is the reverse of that reported in the crystal structure, there is a close correspondence between the heme pocket structure in the crystal and in solution. The proximal histidine coordination geometry is almost identical in the CO complexes of the three monomeric hemoglobins studied. Distal residues are strongly implicated in determining the observed kinetic differences in ligand binding reactions. In particular, steric crowding of the ligand binding site in hemoglobin A is probably a major factor in the slower kinetics of this component.     

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.     

Proton electron nuclear double resonance from nitrosyl horse heart myoglobin: the role of His-E7 and Val-E11

Electron nuclear double resonance (ENDOR) spectroscopy has been used to study protons in nitrosyl horse heart myoglobin (MbNO). (1)H ENDOR spectra were recorded for different settings of the magnetic field. Detailed analysis of the ENDOR powder spectra, using computer simulation, based on the "orientation-selection" principle, leads to the identification of the available protons in the heme pocket. We observe hyperfine interactions of the N(HisF8)-Fe(2+)-N(NO) complex with five protons in axial and with eight protons in the rhombic symmetry along different orientations, including those of the principal axes of the g-tensor. Protons from His-E7 and Val-E11 residues are identified in the two symmetries, rhombic and axial, exhibited by MbNO. Our results indicate that both residues are present inside the heme pocket and help to stabilize one particular conformation.     

NMR studies of the conformations of leghemoglobins from soybean and lupin

Phase-sensitive two-dimensional NMR methods have been used to obtain extensive proton resonance assignments for the carbon monoxide complexes of lupin leghemoglobins I and II and soybean leghemoglobin a. The assigned resonances provide information on the solution conformations of the proteins, particularly in the vicinity of the heme. The structure of the CO complex of lupin leghemoglobin II in solution is compared with the X-ray crystal structure of the cyanide complex by comparison of observed and calculated ring current shifts. The structures are generally very similar but significant differences are observed for the ligand contact residues, Phe30, His63 and Val67, and for the proximal His97 ligand. Certain residues are disordered and adopt two interconverting conformations in lupin leghemoglobin II in solution. The proximal heme pocket structure is closely conserved in the lupin leghemoglobins I and II but small differences in conformation in the distal heme pocket are apparent. Larger conformational differences are observed when comparisons are made with the CO complex of soybean leghemoglobin. Altered protein-heme packing is indicated on the proximal side of the heme and some conformational differences are evident in the distal heme pocket. The small conformational differences between the three leghemoglobins probably contribute to the known differences in their O2 and CO association and dissociation kinetics. The heme pocket conformations of the three leghemoglobins are more closely related to each other than to sperm whale myoglobin. The most notable differences between the leghemoglobins and myoglobin are: (a) reduced steric crowding of the ligand binding site in the leghemoglobins, (b) different orientations of the distal histidine, and (c) small but significant differences in proximal histidine coordination geometry. These changes probably contribute to the large differences in ligand binding kinetics between the leghemoglobins and myoglobin.     

Incorporation of Tyrosine and Glutamine Residues into the Soluble Guanylate Cyclase Heme Distal Pocket Alters NO and O2 Binding

Nitric oxide (NO) is the physiologically relevant activator of the mammalian hemoprotein soluble guanylate cyclase (sGC). The heme cofactor of α1β1 sGC has a high affinity for NO but has never been observed to form a complex with oxygen. Introduction of a key tyrosine residue in the sGC heme binding domain β1(1–385) is sufficient to produce an oxygen-binding protein, but this mutation in the full-length enzyme did not alter oxygen affinity. To evaluate ligand binding specificity in full-length sGC we mutated several conserved distal heme pocket residues (β1 Val-5, Phe-74, Ile-145, and Ile-149) to introduce a hydrogen bond donor in proximity to the heme ligand. We found that the NO coordination state, NO dissociation, and enzyme activation were significantly affected by the presence of a tyrosine in the distal heme pocket; however, the stability of the reduced porphyrin and the proteins affinity for oxygen were unaltered. Recently, an atypical sGC from Drosophila, Gyc-88E, was shown to form a stable complex with oxygen. Sequence analysis of this protein identified two residues in the predicted heme pocket (tyrosine and glutamine) that may function to stabilize oxygen binding in the atypical cyclase. The introduction of these residues into the rat β1 distal heme pocket (Ile-145 → Tyr and Ile-149 → Gln) resulted in an sGC construct that oxidized via an intermediate with an absorbance maximum at 417 nm. This absorbance maximum is consistent with globin Fe^II-O₂ complexes and is likely the first observation of a Fe^II-O₂ complex in the full-length α1β1 protein. Additionally, these data suggest that atypical sGCs stabilize O₂ binding by a hydrogen bonding network involving tyrosine and glutamine.     

Rearrangement of the distal pocket accompanying E7 His----Gln substitution in elephant carbonmonoxy- and oxymyoglobin: 1H NMR identification of a new aromatic residue in the heme pocket

Two-dimensional 1H NMR methods have been used to assign side-chain resonances for the residues in the distal heme pocket of elephant carbonmonoxymyoglobin (MbCO) and oxymyoglobin (MbO2). It is shown that, while the other residues in the heme pocket are minimally perturbed, the Phe CD4 residue in elephant MbCO and MbO2 resonates considerably upfield compared to the corresponding residue in sperm whale MbCO. The new NOE connectivities to Val E11 and heme-induced ring current calculations indicate that Phe CD4 has been inserted into the distal heme pocket by reorienting the aromatic side chain and moving the CD corner closer to the heme. The C zeta H proton of the Phe CD4 was found to move toward the iron of the heme by approximately 4 A relative to the position of sperm whale MbCO, requiring minimally a 3-A movement of the CD helical backbone. The significantly altered distal conformation in elephant myoglobin, rather than the single distal E7 substitution, forms a plausible basis for its altered functional properties of lower autoxidation rate, higher redox potential, and increased affinity for CO ligand. These results demonstrate that one-to-one interpretation of amino acid residue substitution (E7 His----Gln) is oversimplified and that conformational changes of substituted proteins which are not readily predicted have to be considered for interpretation of their functional properties.     

Yeast iso-1-cytochrome c. A 2.8 A resolution three-dimensional structure determination

A molecular replacement approach, augmented with the results of predictive modeling procedures, solvent accessibility studies, packing analyses and translational coefficient searches, has been used to elucidate the 2.8 A (1 A = 0.1 nm) resolution structure of yeast iso-1-cytochrome c. An examination of the polypeptide chain folding of this protein shows it to have unique conformations in three regions, upon comparison with the structures of other eukaryotic cytochromes c. These include: residues -5 to +1 at the N-terminal end of the polypeptide chain, which are in an extended conformation and project in large part off the surface of the protein; residues 19 to 26, which form a surface beta-loop on the His18 ligand side of the central heme group; and, the C-terminal end of the helical segment composed of residues 49 to 56, which serves to form a part of the heme pocket. Structural studies also show that the highly reactive sulfhydryl group of Cys102 is buried within a hydrophobic region in the monomer form of yeast iso-1-cytochrome c. Dimerization of yeast iso-1-cytochrome c through disulfide bond formation between two such residues would require a substantial conformational change in the C-terminal helix of this protein. Another unique structural feature, the trimethylated side-chain of Lys72, is located on the surface of yeast iso-1-cytochrome c near the solvent-exposed edge of the bound heme prosthetic group. On the basis of the results of these and other structural studies, an analysis of the spatial conservation of structural features in the heme pocket of eukaryotic cytochromes c has been conducted. It was found that the residues involved could be divided into three general classes. The current structural analyses and additional modeling studies have also been used to explain the altered functional properties observed for mutant yeast iso-1-cytochrome c proteins.     

NMR studies of nitrophorin distal pocket side chain effects on the heme orientation and seating of NP2 as compared to NP1

The nitrophorins (NP) of the adult blood-sucking insect Rhodnius prolixus fall into two pairs based on sequence identity (NP1,4 (90%) and NP2,3 (79%)), which differ significantly in the size of side chains of residues which contact the heme. These residues include those in the distal pocket of NP2 (I120) and NP1 (T121) and the “belt” that surrounds the heme of NP2 (S40, F42), and NP1(A42, L44). To determine the importance of these residues and others conserved or very similar for the two pairs, including L122(123), L132(133), appropriate mutants of NP2 and NP1 have been prepared and studied by ¹H NMR spectroscopy. Wild-type NP2 has heme orientation ratio (A:B) of 1:8 at equilibrium, while wild-type NP1 has A:B ~ 1:1 at equilibrium. Another difference between NP2 and NP1 is in the heme seating with regard to His57(59). It is found that among the distal pocket residues investigated, the residue most responsible for heme orientation and seating is I120(T121). F42(L44) and L106(F107) may also be important, but must be investigated in greater detail.     

Proton NMR characterization of isomeric sulfmyoglobins: preparation, interconversion, reactivity patterns, and structural features

The preparations of sulfmyoglobin (sulf-Mb) by standard procedures have been found heterogeneous by 1H NMR spectroscopy. Presented here are the results of a comprehensive study of the factors that influence the selection among the three dominant isomeric forms of sperm whale sulf-Mb and their resulting detailed optical and 1H NMR properties as related to their detectability and structural properties of the heme pocket. A single isomer is formed initially in the deoxy state; further treatment in any desired oxidation/ligation state can yield two other major isomers. Acid catalysis and chromatography facilitate formation of a second isomer, particularly in the high-spin state. At neutral pH, a third isomer is formed by a first-order process. The processes that alter oxidation/ligation state are found to be reversible and are judged to affect only the metal center, but the three isomeric sulf-Mbs are found to exhibit significantly different ligand affinity and chemical stability. The present results allow, for the first time, a rational approach for preparing a given isomeric sulf-Mb in an optimally pure state for subsequent characterization by other techniques. While optical spectroscopy can distinguish the alkaline forms, only 1H NMR clearly distinguishes all three ferric isomers. The ring current shifts in the carbonyl complexes of reduced sulf-Mb complexes support saturation for a pyrrole in each isomer. The hyperfine shift patterns in the various oxidation/spin states of sulf-Mbs indicate relatively small structural alteration, and the proximal and distal sides of the heme suggest that peripheral electronic effects are responsible for the differentially reduced ligand affinities for the three isomeric sulf-Mbs. The first 1H NMR spectra of sulfhemoglobins are presented, which indicate a structure similar to that of the initially formed sulf-Mb isomer but also suggest the presence of a similar molecular heterogeneity as found for sulf-Mb, albeit to a smaller extent.     

Restricting the ligand-linked heme movement in Scapharca dimeric hemoglobin reveals tight coupling between distal and proximal contributions to cooperativity.

Cooperative ligand binding in the dimeric hemoglobin from the blood clam Scapharca inaequivalvis results primarily from tertiary, rather than quaternary, structural changes. Ligand binding is coupled with conformational changes of key residues, including Phe 97, which is extruded from the proximal heme pocket, and the heme group, which moves deeper into the heme pocket. We have tested the role of the heme movement in cooperative function by mutating Ile 114, at the base of the heme pocket. Replacement of this residue with a Met did not disturb the hemoglobin structure or significantly alter equilibrium ligand binding properties. In contrast, substitution with a Phe at position 114 inhibits the ligand-linked movement of the heme group, and substantially reduces oxygen affinity and cooperativity. As the extent of heme movement to the normal position of the ligated state is diminished, Phe 97 is inhibited from its movement into the interface upon ligand binding. These results indicate a tight coupling between these two key cooperative transitions and suggest that the heme movement may be an obligatory trigger for expulsion of Phe 97 from the heme pocket.     

Crystal structure of heme oxygenase from the gram-negative pathogen Neisseria meningitidis and a comparison with mammalian heme oxygenase-1

We report the crystal structure of heme oxygenase from the pathogenic bacterium Neisseria meningitidis at 1.5 A and compare and contrast it with known structures of heme oxygenase-1 from mammalian sources. Both the bacterial and mammalian enzymes share the same overall fold, with a histidine contributing a ligand to the proximal side of the heme iron and a kinked alpha-helix defining the distal pocket. The distal helix differs noticeably in both sequence and conformation, and the distal pocket of the Neisseria enzyme is substantially smaller than in the mammalian enzyme. Key glycine residues provide the flexibility for the helical kink, allow close contact of the helix backbone with the heme, and may interact directly with heme ligands.     

A model for the misfolded bis-His intermediate of cytochrome c: the 1-56 N-fragment

We have characterized the ferric and ferrous forms of the heme-containing (1-56 residues) N-fragment of horse heart cytochrome c (cyt c) at different pH values and low ionic strength by UV-visible absorption and resonance Raman (RR) scattering. The results are compared with native cyt c in the same experimental conditions as this may provide a deeper insight into the cyt c unfolding-folding process. Folding of cyt c leads to a state having the heme iron coordinated to a histidine (His18) and a methionine (Met80) as axial ligands. At neutral pH the N-fragment (which lacks Met80) shows absorption and RR spectra that are consistent with the presence of a bis-His low spin heme, like several non-native forms of the parental protein. In particular, the optical spectra are identical to those of cyt c in the presence of a high concentration of denaturants; this renders the N-fragment a suitable model to study the heme pocket microenvironment of the misfolded (His-His) intermediate formed during folding of cyt c. Acid pH affects the ligation state in both cyt c and the N-fragment. Data obtained as a function of pH allow a correlation between the structural properties in the heme pocket of the N-fragment and those of non-native forms of cyt c. The results underline that the (57-104 residues) segment under native-like conditions imparts structural stability to the protein by impeding solvent access into the heme pocket.     

Functional interpretation and structural insights of Arabidopsis lyrata cytochrome P450 CYP71A13 involved in auxin synthesis

Cytochrome P450 CYP71A13 of Arabidopsis lyrata is a heme protein involved in biosynthesis of indole-3-acetonitrile which leads to the formation of indolyl-3-acetic acid. It catalyzes a unique reaction: formation of a carbon-nitrogen triple bond and dehydration of indolyl-3-acetaldoxime. Homology model of this 57 kDa polypeptide revealed that the heme existed between H-helix and J- helix in the hydrophobic pocket, although both helixes are involved in catalytic activity, where Gly305 and Thr308, 311 of H- helix were involved in its stabilization. The substrate indole-3-acetaldoxime was tightly fitted into the substrate pocket with the aromatic ring being surrounded by amino acid residues creating a hydrophobic environment. The smaller size of the substrate binding pocket in cytochrome P450 CYP71A13 was due to the bulkiness of the two amino acid residues Phe182 and Trp315 pointing into the substrate binding cavity. The apparent role of the heme in cytochrome P450 CYP71A13 was to tether the substrate in the catalysis by indole-3-acetaldoxime dehydratase. Since the crystal structure of cytochrome P450 CYP71A13 has not yet been solved, the modeled structure revealed mechanism of substrate recognition and catalysis.     

Novel X-ray sequences and crystal structures of Persian and Starry sturgeon methemoglobins: Highlighting the role of heme pocket waters in causing autoxidation

Although there is a high sequence similarity between mammalian and fish hemoglobin (Hb), the oxidation and heme loss rates can vary greatly between them such that fish Hbs oxidise much more rapidly than mammalian Hbs. There is to date no sequence or structural data for any sturgeon Hb to reveal the level of autoxidation in these fish. In this study, novel high resolution X-ray sequences and crystal structures of methemoglobin (Met-Hb) from two sturgeon fish including Persian sturgeon (Acipenser percisus) and Starry sturgeon (Acipenser stellatus) belonging to the Caspian sea has been determined. A comprehensive sequence and structure comparison between these sturgeon Met-Hbs and a number of non-sturgeon and normal and sickle cell anaemia human Hb in varying heme states has been carried out highlighting (i) the structural variability in the heme propionate groups; (ii) the existence of certain residues or their displacement and shift in the heme pocket allowing entry of water molecules into the heme pocket; (iii) the importance of the number of water molecules in the heme pocket; (iv) the hydrogen bonding between oxygens of A and D propionate groups and that of waters in the heme pocket; and (v) the role of heme binding waters causing oxidative stress and heme autoxidation.     

Ligand selectivity of soluble guanylyl cyclase: effect of the hydrogen-bonding tyrosine in the distal heme pocket on binding of oxygen, nitric oxide, and carbon monoxide

Although soluble guanylyl cyclase (sGC) functions in an environment in which O(2), NO, and CO are potential ligands for its heme moiety, the enzyme displays a high affinity for only its physiological ligand, NO, but has a limited affinity for CO and no affinity for O(2). Recent studies of a truncated version of the sGC beta(1)-subunit containing the heme-binding domain (Boon, E. M., Huang, S H., and Marletta, M. A. (2005) Nat. Chem. Biol., 1, 53-59) showed that introduction of the hydrogen-bonding tyrosine into the distal heme pocket changes the ligand specificity of the heme moiety and results in an oxygen-binding sGC. The hypothesis that the absence of hydrogen-bonding residues in the distal heme pocket is sufficient to provide oxygen discrimination by sGC was put forward. We tested this hypothesis in a context of a complete sGC heterodimer containing both the intact alpha(1)- and beta(1)-subunits. We found that the I145Y substitution in the full-length beta-subunit of the sGC heterodimer did not produce an oxygen-binding enzyme. However, this substitution impeded the association of NO and destabilized the NO.heme complex. The tyrosine in the distal heme pocket also impeded both the binding and dissociation of the CO ligand. We propose that the mechanism of oxygen exclusion by sGC not only involves the lack of hydrogen bonding in the distal heme pocket, but also depends on structural elements from other domains of sGC.     

One- and two-dimensional NMR investigations of the heme pocket in free alpha(CO) chains from human hemoglobin

Two-dimensional nuclear magnetic resonance techniques were used to assign resonances corresponding to heme pocket residues of the isolated alpha(CO) subunits of the human adult hemoglobin (HbA). The assignment procedure was based on the partial identification of the amino acid spin system from the J-correlated (COSY) spectrum and on the nuclear Overhauser effect connectivities (from NOSEY spectra) with the heme substituents. We present here partial assignments corresponding to five amino acid residues: Leu86, Leu-91, Val-93, Leu-101 and Leu-136. Starting from the known crystallographic structure of the alpha subunit in the hemoglobin tetramer, we applied a dipolar model to compute the ring-current shift of the protons from fifteen amino acid residues in the heme pocket. Comparison of the predicted and observed chemical shifts suggests that there is a very close similarity between the heme pocket tertiary structure of the alpha(CO) subunits in crystals of HbA(CO) and of the free alpha(CO) chains. The one-dimensional NMR spectra were used to monitor the pH-induced structural changes, the effects of chemical modification and of ligand substitution. Upon increasing the pH from 5.6 to 9.0 the structure of the heme environment appears to be invariant with the exception of some residues in the CD corner. The structure is also largely conserved when p-chloromercuribenzoate is bound to Cys-104. In contrast, the substitution of CO by O2 as ligand induces many large changes in the heme cavity which can be partially characterized by NMR spectroscopy.     

Comparison of apo- and heme-bound crystal structures of a truncated human heme oxygenase-2

Heme oxygenase (HO) catalyzes the first step in the heme degradation pathway. The crystal structures of apo- and heme-bound truncated human HO-2 reveal a primarily alpha-helical architecture similar to that of human HO-1 and other known HOs. Proper orientation of heme in HO-2 is required for the regioselective oxidation of the alpha-mesocarbon. This is accomplished by interactions within the heme binding pocket, which is made up of two helices. The iron coordinating residue, His(45), resides on the proximal helix. The distal helix contains highly conserved glycine residues that allow the helix to flex and interact with the bound heme. Tyr(154), Lys(199), and Arg(203) orient the heme through direct interactions with the heme propionates. The rearrangements of side chains in heme-bound HO-2 compared with apoHO-2 further elucidate HO-2 heme interactions.     

Substrate flow in catalases deduced from the crystal structures of active site variants of HPII from Escherichia coli.

The active site of heme catalases is buried deep inside a structurally highly conserved homotetramer. Channels leading to the active site have been identified as potential routes for substrate flow and product release, although evidence in support of this model is limited. To investigate further the role of protein structure and molecular channels in catalysis, the crystal structures of four active site variants of catalase HPII from Escherichia coli (His128Ala, His128Asn, Asn201Ala, and Asn201His) have been determined at approximately 2.0-A resolution. The solvent organization shows major rearrangements with respect to native HPII, not only in the vicinity of the replaced residues but also in the main molecular channel leading to the heme distal pocket. In the two inactive His128 variants, continuous chains of hydrogen bonded water molecules extend from the molecular surface to the heme distal pocket filling the main channel. The differences in continuity of solvent molecules between the native and variant structures illustrate how sensitive the solvent matrix is to subtle changes in structure. It is hypothesized that the slightly larger H(2)O(2) passing through the channel of the native enzyme will promote the formation of a continuous chain of solvent and peroxide. The structure of the His128Asn variant complexed with hydrogen peroxide has also been determined at 2.3-A resolution, revealing the existence of hydrogen peroxide binding sites both in the heme distal pocket and in the main channel. Unexpectedly, the largest changes in protein structure resulting from peroxide binding are clustered on the heme proximal side and mainly involve residues in only two subunits, leading to a departure from the 222-point group symmetry of the native enzyme. An active role for channels in the selective flow of substrates through the catalase molecule is proposed as an integral feature of the catalytic mechanism. The Asn201His variant of HPII was found to contain unoxidized heme b in combination with the proximal side His-Tyr bond suggesting that the mechanistic pathways of the two reactions can be uncoupled.