Temperature dependence of Q-band electron paramagnetic resonance spectra of nitrosyl heme proteins

The Q-band (35 GHz) electron paramagnetic resonance (EPR) spectra of nitrosyl hemoglobin (HbNO) and nitrosyl myoglobin (MbNO) were studied as a function of temperature between 19 K and 200 K. The spectra of both heme proteins show two classes of variations as a function of temperature. The first one has previously been associated with the existence of two paramagnetic species, one with rhombic and the other with axial symmetry. The second one manifests itself in changes in the g-factors and linewidths of each species. These changes are correlated with the conformational substates model and associate the variations of g-values with changes in the angle of the N(his)-Fe-N(NO) bond in the rhombic species and with changes in the distance between Fe and N of the proximal (F8) histidine in the axial species.     

FTIR and resonance Raman studies of nitric oxide binding to H93G cavity mutants of myoglobin.

Nitric oxide (NO) binds to the myoglobin (Mb) cavity mutant, H93G, forming either a five- or six-coordinate Fe-NO complex. The H93G mutation eliminates the covalent attachment between the protein and the proximal ligand, allowing NO to bind H93G possibly from the proximal side of the heme rather than the typical diatomic binding pocket on the distal side. The question of whether NO binds on the distal or proximal side was addressed by FTIR spectroscopy of the N-O vibrational frequency nuN(-O) for a set of Mb mutants that perturb the electrostatic environment of the heme pocket. Vibrational spectra of five- and six-coordinate MbNO complexes indicate that nu(N-O) shifts (by as much as 26 cm(-1)) to higher energies for the distal mutants H64V and H64V/H93G relative to the energies of wild-type and H93G MbNO, while nu(N-O) is not affected by the proximal side mutation S92A/H93G. This result suggests that NO binds on the distal side of heme in the five- and six-coordinate MbNO complexes of H93G. Additionally, values of the Fe-NO vibrational frequency nu(Fe-NO) as measured by resonance Raman spectroscopy are reported for the distal and proximal double mutants of H93G. These results suggest that nu(Fe-NO) is not very sensitive to mutations that perturb the electrostatic environment of the heme pocket, leading to the observation that nu(N-O) and nu(Fe-NO) are not quantitatively correlated for the MbNO complexes presented here. Furthermore, nu(N-O) and nu(Fe-NO) do not correlate well with equilibrium constants for imidazole binding to the five-coordinate MbNO complexes of the H93G double mutants. The data presented here do not appear to support the presence of pi-back-bonding or an inverse trans effect of NO binding in Mb mutants that alter the electrostatic environment of the heme pocket.     

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.     

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.     

An EPR study of myeloperoxidase in human granulocytes.

1. EPR spectra of human granulocytes (4 - 10(8) cells per ml) show an intense high-spin ferric heme signal with rhombic symmetry (gx = 6.90 and gy = 5.07) for the heme group. These g-values are identical to those of partially purified myeloperoxidase and thus the signal is derived from ferric myeloperoxidase. In chicken granulocytes, which contain little or no myeloperoxidase, only an axial type of heme iron signal, weak in intensity, can be detected at g = 6.0. 2. Upon phagocytosis of latex particles by human granulocytes the high-spin heme signal with rhombic symmetry is slowly converted into a signal with axial symmetry (gx = gy = 6.0), showing that the EPR signals of myeloperoxidase in the intact cell can be used to study the involvement of the enzyme in metabolic changes during phagocytosis.     

Crystal structures of ferrous horse heart myoglobin complexed with nitric oxide and nitrosoethane

The interactions of nitric oxide (NO) and organic nitroso compounds with heme proteins are biologically important, and adduct formation between NO-containing compounds and myoglobin (Mb) have served as prototypical systems for studies of these interactions. We have prepared crystals of horse heart (hh) MbNO from nitrosylation of aqua-metMb crystals, and we have determined the crystal structure of hh MbNO at a resolution of 1.9 A. The Fe-N-O angle of 147 degrees in hh MbNO is larger than the corresponding 112 degrees angle previously determined from the crystal structure of sperm whale MbNO (Brucker et al., Proteins 1998;30:352-356) but is similar to the 150 degrees angle determined from a MS XAFS study of a frozen solution of hh MbNO (Rich et al., J Am Chem Soc 1998;120:10827-10836). The Fe-N(O) bond length of 2.0 A (this work) is longer than the 1.75 A distance determined from the XAFS study and suggests distal pocket influences on FeNO geometry. The nitrosyl N atom is located 3.0 A from the imidazole N(epsilon) atom of the distal His64 residue, suggesting electrostatic stabilization of the FeNO moiety by His64. The crystal structure of the nitrosoethane adduct of ferrous hh Mb was determined at a resolution of 1.7 A. The nitroso O atom of the EtNO ligand is located 2.7 A from the imidazole N(epsilon) atom of His64, suggesting a hydrogen bond interaction between these groups. To the best of our knowledge, the crystal structure of hh Mb(EtNO) is the first such determination of a nitrosoalkane adduct of a heme protein.     

1H NMR studies of the effect of mutation at Valine45 on heme microenvironment of cytochrome b5

1D and 2D (1)H NMR were employed to probe the effects on the heme microenvironment of cytochrome b(5) caused by the mutation from Val45 to Tyr45, His45 and Glu45. Compared with wild type (WT) cytochrome b(5), in all mutants the heme ring are CCW rotated relative to the imidazole planes of axial ligands and the angles beta between two axial ligand imidazole planes are not changed, being in agreement with the temperature dependence of the shifts of the heme protons. The ratios of heme isomers (major to minor) are smaller than that in WT. The 4-vinyl group of the heme in V45Y assumes cis-orientation, being similar to that of WT, while in V45E and V45H, both cis and trans orientation are found. The relationships between the structure and biological function of the mutants are discussed in terms of the geometry of heme and axial ligands, the hydrophobicity of heme pocket and the electrostatic potential of the heme-exposed area.     

2D NMR of paramagnetic metalloenzymes: Cyanide-inhibited horseradish peroxidase

Summary Two-dimensional (2D) proton NMR correlation spectroscopy, COSY, and nuclear Overhauser spectroscopy, NOESY, have been used to explore the applicability of these methods for the moderately large (42 KDa), paramagnetic cyanide-inhibited derivative of horseradish peroxidase, HRP-CN. The target resonances are those in the active site of HRP-CN which experience substantial hyperfine shifts and paramagnetic relaxation. The magnitude COSY experiment was found to yield cross peaks for all known spin-coupled heme substituents, as well as for the majority of non-heme hyperfine shifted protons, in spite of line widths of the order of 100 Hz. Moreover, the rapid relaxation of the hyperfine-shifted resonances allows the extremely rapid collection of useful 2D NMR data sets without the loss of information. For the heme, the combination of COSY cross peaks for the vinyl and propionate substituents, and NOESY cross peaks among these substituent protons and heme methyls, allows assignment of heme resonances without recourse to deuterium labeling of the heme. A seven-proton coupled spin system was identified in the upfield region that is consistent with originating from the proposed catalytic Arg³⁸ residue in the distal heme pocket, with orientation relative to the heme similar to that found in cytochromec peroxidase. The upfield hyperfine-shifted methyl group in the substrate binding pocket previously proposed to arise from Leu²³⁷ is shown to arise instead from an as yet unidentified Ile. NOESY spectra collected at very short (3 ms) and intermediate (20 ms) mixing times indicate that build-up curves can be obtained that should yield estimates of distances in the heme cavity. It is concluded that 2D NMR studies should be able to provide the heme assignments, aid in identifying the catalytic residues, and provide information on the spatial disposition of such residues in the active site for cyanide complexes of a number of intermediate to large paramagnetic heme peroxidases, as well as for other paramagnetic metalloenzymes with line widths of 100 Hz. Moreover, paramagnetic-induced hyperfine shifts and linewidths to 100 Hz need not interfere with the complete solution structure determination of a small paramagnetic protein solely on the basis of 2D NMR data.     

ENDOR from nitrogens and protons in low spin ferric heme and hemoprotein.

Electron nuclear double resonance (ENDOR) signals have been obtained from iron-linked nitrogens in frozen solutions of cytochrome c, metmyoglobin cyanide, and a low spin protohemin mercaptide complex. Hyperfine couplings from heme protons have also been obtained from metmyoglobin cyanide and from a low spin protohemin cyanide complex. Several of these proton resonances are assigned to specific heme protons.     

High-pressure 1H NMR study of pressure-induced structural changes in the heme environments of metcyanomyoglobins

The effect of pressure on the heme environment structure of sperm whale and horse heart metcyanomyoglobins was investigated up to 300 MPa by high-pressure (1)H NMR spectroscopy. Pressure-induced changes in the distances between the observed protons and the heme iron atom were estimated from changes in the dipolar shift due to the paramagnetic effect on the protons. The changes showed that the heme peripheral structure as a whole was compressed by pressure; the movements of the protons in the heme peripheral residues were in the range of +0.16 to -0.54 A/300 MPa. One-dimensional compressibilities for the protons, excluding the protons of the distal His residue, were in the range of 1.0 x 10(-4) to 6.1 x 10(-4)/MPa. The movements of the protons induced by pressure correlated well with the distance between the protons and cavities in the protein. The distal His residue (His 64) moved toward the outside of the heme pocket, but remained in the pocket even at 300 MPa. This movement was driven dominantly by a change in the dihedral angle around the C(alpha)-C(beta) rotational bond of the residue. Comparative work on horse heart metcyanomyoglobin implied that the conformational change of the His 64 imidazole ring was larger in the horse heart metcyanomyoglobin than in the sperm whale metcyanomyoglobin.     

Assignment of heme and distal amino acid resonances in the 1H-NMR spectra of the carbon monoxide and oxygen complexes of sperm whale myoglobin

Assignments of resonances of the heme and distal amino acid protons in spectra of the CO and O2 complexes of sperm whale myoglobin are reported. These resonances provide information on the conformation of the heme pocket. For oxymyoglobin, the assignments of the heme meso protons disagree with those proposed previously on the basis of partial deuteration experiments. Rapid ring flips about the C beta-C gamma bond are detected for Phe-CD1. Recent claims for two conformational substates of valine-E11 in carbonmonoxymyoglobin (Bradbury, J.H. and Carver, J.A. (1984) Biochemistry 23, 4905-4913) are shown to be in error. The pK of His-97 (FG3) in carbonmonoxymyoglobin has been determined (pK = 5.9). This residue appears to influence many spectroscopic properties of myoglobin. The distal His-E7 in carbonmonoxymyoglobin has pK less than 5.0. Differences in the heme pocket conformation in the CO complexes of myoglobin and leghemoglobin are discussed. These differences may be influential in O2 and CO association reactions.     

Electron nuclear double resonance and electron paramagnetic resonance study on the structure of the NO-ligated heme alpha 3 in cytochrome c oxidase

The techniques of EPR and electron nuclear double resonance (ENDOR) were used to probe structure and electronic distribution at the nitric oxide (NO)-ligated heme alpha 3 in the nitrosylferrocytochrome alpha 3 moiety of fully reduced cytochrome c oxidase. Hyperfine and quadrupole couplings to NO (in both 15NO and 14NO forms), to histidine nitrogens, and to protons near the heme site were obtained. Parallel studies were also performed on NO-ligated myoglobin and model NO-heme-imidazole systems. The major findings and interpretations on nitrosylferrocytochrome alpha 3 were: 1) compared to other NO-heme-imidazole systems, the nitrosylferrocytochrome alpha3 gave better resolution of EPR and ENDOR signals; 2) at the maximal g value (gx = 2.09), particularly well resolved NO nitrogen hyperfine and quadrupole couplings and mesoproton hyperfine couplings were seen. These hyperfine and quadrupole couplings gave information on the electronic distribution on the NO, on the orientation of the g tensor with respect to the heme, and possibly on the orientation of the FeNO plane; 3) a combination of experimental EPR-ENDOR results and EPR spectral simulations evidenced a rotation of the NO hyperfine tensor with respect to the electronic g tensor; this implied a bent Fe-NO bond; 4) ENDOR showed a unique proton not seen in the other NO heme systems studied. The magnitude of this proton's hyperfine coupling was consistent with this proton being part of a nearby protein side chain that perturbs an axial ligand like NO or O2.     

Denaturation and renaturation of myeloperoxidase. Consequences for the nature of the prosthetic group.

The effects of the chaotropic agent, guanidine HCl, on the chlorinating activity, optical absorption, EPR, and resonance Raman spectra of myeloperoxidase have been studied. In the presence of the agent the Soret optical absorption of the reduced enzyme (lambda max = 474 nm) is blue shifted to 448 nm, a position similar to heme alpha-containing enzymes. The chlorinating activity of the enzyme disappears, and EPR spectra show a loss of intensity of the rhombic high spin heme signals (gx = 6.9; gy = 5.4) and the appearance of a more axial high spin signal (gx = gy = 6.0). Surprisingly the effects of guanidine HCl are partly reversible. Upon decreasing the concentration of the chaotropic agents by dilution, both the chlorinating activity and the original optical spectrum of native reduced enzyme (lambda max = 474 nm) are partly restored. The resonance Raman spectra of denatured cyanomyeloperoxidase are less complicated than those of native myeloperoxidase, which have been interpreted previously to suggest an iron chlorin chromophore. The multiple lines in the oxidation state marker region are not seen in the spectra of the denatured species. The changes suggest that upon denaturation the macrocycle is converted into a more symmetric structure. Since the effects on the optical absorption spectrum are reversible we speculate that, in the native enzyme, an apparent porphyrin macrocycle undergoes a reversible interaction with amino acid residues in the protein which creates an asymmetry in the electronic distribution of the macrocycle. Comparison of the Raman spectra of denatured cyanomyeloperoxidase with those of analogous heme alpha model complexes suggests the presence of a formyl group in the denatured species; our data, however, demonstrate that the chromophore structure is not identical to heme alpha and may contain a different C beta substitution on the ring macrocycle.     

Crystal structure of Nitrosomonas europaea cytochrome c peroxidase and the structural basis for ligand switching in bacterial di-heme peroxidases.

The crystal structure of the fully oxidized di-heme peroxidase from Nitrosomonas europaea has been solved to a resolution of 1.80 A and compared to the closely related enzyme from Pseudomonas aeruginosa. Both enzymes catalyze the peroxide-dependent oxidation of a protein electron donor such as cytochrome c. Electrons enter the enzyme through the high-potential heme followed by electron transfer to the low-potential heme, the site of peroxide activation. Both enzymes form homodimers, each of which folds into two distinct heme domains. Each heme is held in place by thioether bonds between the heme vinyl groups and Cys residues. The high-potential heme in both enzymes has Met and His as axial heme ligands. In the Pseudomonas enzyme, the low-potential heme has two His residues as axial heme ligands [Fulop et al. (1995) Structure 3, 1225-1233]. Since the site of reaction with peroxide is the low-potential heme, then one His ligand must first dissociate. In sharp contrast, the low-potential heme in the Nitrosomonas enzyme already is in the "activated" state with only one His ligand and an open distal axial ligation position available for reaction with peroxide. A comparison between the two enzymes illustrates the range of conformational changes required to activate the Pseudomonas enzyme. This change involves a large motion of a loop containing the dissociable His ligand from the heme pocket to the molecular surface where it forms part of the dimer interface. Since the Nitrosomonas enzyme is in the active state, the structure provides some insights on residues involved in peroxide activation. Most importantly, a Glu residue situated near the peroxide binding site could possibly serve as an acid-base catalytic group required for cleavage of the peroxide O--O bond.     

Novel heme ligation in a c-type cytochrome involved in thiosulfate oxidation: EPR and MCD of SoxAX from Rhodovulum sulfidophilum

The SoxAX complex of the bacterium Rhodovulum sulfidophilum is a heterodimeric c-type cytochrome that plays an essential role in photosynthetic thiosulfate and sulfide oxidation. The three heme sites of SoxAX have been analyzed using electronic absorption, electron paramagnetic resonance, and magnetic circular dichroism spectroscopies. Heme-3 in the ferric state is characterized by a Large g(max) EPR signal and has histidine and methionine axial heme iron ligands which are retained on reduction to the ferrous state. Hemes-1 and -2 both have thiolate plus nitrogenous ligand sets in the ferric state and give rise to rhombic EPR spectra. Heme-1, whose ligands derive from cysteinate and histidine residues, remains ferric in the presence of dithionite ion. Ferric heme-2 exists with a preparation-dependent mixture of two different ligand sets, one being cysteinate/histidine, the other an unidentified pair with a weaker crystal-field strength. Upon reduction of the SoxAX complex with dithionite, a change occurs in the ligands of heme-2 in which the thiolate is either protonated or replaced by an unidentified ligand. Sequence analysis places the histidine/methionine-coordinated heme in SoxX and the thiolate-liganded hemes in SoxA. SoxAX is the first naturally occurring c-type cytochrome in which a thiolate-coordinated heme has been identified.     

Characterization of the heme environmental structure of cytoglobin,a fourth globin in humans

Cytoglobin (Cgb) represents a fourth member of the globin superfamily in mammals, but its function is unknown. Site-directed mutagenesis, in which six histidine residues were replaced with alanine, was carried out, and the results indicate that the imidazoles of His81 (E7) and His113 (F8) bind to the heme iron as axial ligands in the hexacoordinate and the low-spin state. The optical absorption, resonance Raman, and IR spectral results are consistent with this conclusion. The redox potential measurements revealed an E' of 20 mV (vs NHE) in the ferric/ferrous couple, indicating that the imidazole ligands of His81 and His113 are electronically neutral. On the basis of the nu(Fe-CO) and nu(C-O) values in the resonance Raman and infrared spectra of the ferrous-CO complexes of Cgb and its mutants, it was found that CO binds to the ferrous iron after the His81 imidazole is dissociated, and three conformers are present in the resultant CO coordination structure. Two are in closed conformations of the heme pocket, in which the bound CO ligand interacts with the dissociated His81 imidazole, while the third is in an open conformation. The nu(Fe-O2) in the resonance Raman spectra of oxy Cgb can be observed at 572 cm(-1), suggesting a polar heme environment. These structural properties of the heme pocket of Cgb are discussed with respect to its proposed in vivo oxygen storage function.     

Accessibility of the cytochrome a heme in cytochrome c oxidase to exchangeable protons

The comparison of the resonance Raman spectrum of cytochrome a2+ from cytochrome oxidase in deuterated buffers to that in protonated buffers reveals many lines that have different frequency or intensity. Some of the frequency differences are very large, e.g. on the order of 10 cm-1. From these differences in the Raman spectra, we infer that the heme pocket is readily accessible to protons and that labile groups are either on the heme or interact strongly with it. These data suggest the possibility of direct participation in proton translocation and/or oxygen protonation by the heme of cytochrome a.     

The binding of molybdate to uteroferrin. Hyperfine interactions of the binuclear center with 95Mo, 1H, and 2H

Uteroferrin, an acid phosphatase with a spin-coupled and redox-active binuclear iron center, is paramagnetic in its pink, enzymatically active, mixed-valence (S = 1/2) state. Phosphate, a product and inhibitor of the enzymatic activity of uteroferrin, converts the pink, EPR-active form of the protein to a purple, EPR-silent species. In contrast, molybdate, a tetrahedral oxyanion analog of phosphate, transforms the EPR spectrum of uteroferrin from a rhombic to an axial form. With both electron spin echo envelope modulation (ESEEM) and electron nuclear double resonance (ENDOR) spectroscopies, we observe a hyperfine interaction of [95Mo]molybdate with the S = 1/2, Fe(II)-Fe(III) center of the protein. A pair of 95Mo resonances centered at the 95Mo Larmor frequency at the applied magnetic field and separated by a hyperfine coupling constant of 1.2 MHz is evident. Therefore, a single monomeric species of molybdate is close to, and likely a ligand of, the binuclear cluster. 1H ENDOR studies on uteroferrin reveal at least six sets of lines mirrored about the 1H Larmor frequency. Two pairs of these lines become reduced in intensity when the protein is exchanged against D2O. Moreover, ESEEM and 2H ENDOR spectra display resonances at the 2H Larmor frequency. Therefore, the metal-binding region of the protein is accessible to solvent. Additional deuterium lines observable by ESEEM spectroscopy provide evidence for a population of strongly coupled, readily exchangeable protons associated with the binuclear center. The measured hyperfine coupling constants for these deuterons are orientation-dependent with splittings of nearly 4 MHz at g3 = 1.59 and less than 1 MHz at g1 = 1.94. In the presence of molybdate, ESEEM spectra of D2O-exchanged samples reveal a resonance at the 2H Larmor frequency, with no evidence of spectral components due to strongly coupled deuterons. 1H ENDOR studies of the uteroferrin-molybdate complex show at least seven pairs of lines, mirrored about the 1H Larmor frequency, of which one pair becomes attenuated in amplitude upon deuteration. The active site thus remains accessible to solvent in the presence of molybdate.     

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.     

Hydrogenases in the "active" state: determination of g-matrix axes and electron spin distribution at the active site by 1H ENDOR spectroscopy

Hydrons and electrons are substrates for the enzyme hydrogenase, but cannot be observed in X-ray crystal structures. High-resolution 1H electron nuclear double resonance (ENDOR) spectroscopy offers a means to detect the distribution of protons and unpaired electrons. ENDOR spectra were recorded from frozen solutions of the nickel-iron hydrogenases of Desulfovibrio gigas and Desulfomicrobium baculatum, in the "active" state ("Ni-C" EPR signal) and analyzed by orientationally selective simulation methods. The experimental spectra were fitted using a structural model of the nickel-iron centre based on crystallographic results, allowing for differences in electron spin distribution as well as the spatial orientation of the g-matrix ( g-tensor), and anisotropic and isotropic hyperfine couplings of the protons nearest to the nickel ion. ENDOR signals, detected after complete deuterium exchange, were assigned to six protons of the cysteines bound to nickel. The assignment took advantage of the substitution of a selenium for a sulfur ligand, which occurs naturally between the [NiFeSe] and [NiFe] hydrogenases from Dm. baculatum and D. gigas, respectively, and was found to affect just two signals. The four signals with the largest hyperfine couplings, including isotropic contributions from 4.5 to 13.5 MHz, were assigned to the beta-methylene protons of the two terminal cysteine ligands, one of which is substituted by seleno-cysteine in [NiFeSe] hydrogenase. The electron spin is delocalized onto the nickel (50%) and its sulfur ligands, with a higher proportion on the terminal than the bridging ligands. The g-matrix was found to align with the active site in such a way that the g1- g2 plane is nearly coplanar (18.3 degrees) with the plane defined by nickel and three sulfur atoms, and the g2 axis deviates by 22.9 degrees from the vector between nickel and iron. Significantly for the reaction of the enzyme, direct evidence for the binding of hydrons at the active site was obtained by the detection of H/D-exchangeable ENDOR signals.