TY - JOUR
T1 - A Density Functional Theory Investigation of Fe-N-O Bonding in Heme Proteins and Model Systems
AU - Zhang, Yong
AU - Gossman, William
AU - Oldfield, Eric
PY - 2003/12/31
Y1 - 2003/12/31
N2 - We report the results of a series of density functional theory (DFT) calculations of the Mössbauer quadrupole splittings and isomer shifts in NO heme model compounds, together with the results of calculations of the Mössbauer quadrupole splittings, isomer shifts, and electron paramagnetic resonance hyperfine coupling constants in a model Fe(II)(NO)(imidazole) complex as a function of Fe-NO bond length and Fe-N-O bond angle. The results of the Mössbauer quadrupole splitting and isomer shift calculations on the NO heme model compounds show good accord between theory and experiment, with the largest errors being observed for structures having the largest crystallographic R1 values. The results of the property surface calculations were then used to calculate Fe-NO bond length and Fe-N-O bond angle probability surfaces (Z-surfaces) for a nitrosyl hemoglobin, using, in addition, an energy filter. The results obtained yielded a most probable Fe-NO bond length (r) of 1.79 Å and an Fe-N-O bond angle (β) of 136°-137°. This bond length is somewhat longer than those observed in most model compounds but may be due, at least in part, to hydrogen bond formation with the distal His residue. Bond elongation was also observed in a geometry optimized Fe(II)(NO)(imidazole) complex hydrogen bonded to an imidazole residue, in which we find r = 1.76-1.78 Å and β = 137°-138°. The computed bond angles are close to the canonical ∼140° value found in most model systems. Highly bent Fe-N-O bond angles or very long Fe-NO bond lengths seem unlikely to occur in proteins, due to their high energies. We also investigated the molecular orbitals and spin densities in each of the six coordinate systems investigated and found the orbitals and spin densities to be generally similar those described previously for five coordinate systems. Taken together, these results show that Mössbauer quadrupole splittings and isomer shifts, in addition to electron paramagnetic resonance hyperfine coupling constants, can now be calculated for nitrosyl heme systems with relatively good accuracy and that the results so obtained can be used to determine Fe-N-O geometries in metalloproteins. The Z-surface approach is thus applicable to both diamagnetic (CO) and paramagnetic (NO) heme proteins with in both cases the metal-ligand binding geometries found in the proteins being very close to those seen in model systems.
AB - We report the results of a series of density functional theory (DFT) calculations of the Mössbauer quadrupole splittings and isomer shifts in NO heme model compounds, together with the results of calculations of the Mössbauer quadrupole splittings, isomer shifts, and electron paramagnetic resonance hyperfine coupling constants in a model Fe(II)(NO)(imidazole) complex as a function of Fe-NO bond length and Fe-N-O bond angle. The results of the Mössbauer quadrupole splitting and isomer shift calculations on the NO heme model compounds show good accord between theory and experiment, with the largest errors being observed for structures having the largest crystallographic R1 values. The results of the property surface calculations were then used to calculate Fe-NO bond length and Fe-N-O bond angle probability surfaces (Z-surfaces) for a nitrosyl hemoglobin, using, in addition, an energy filter. The results obtained yielded a most probable Fe-NO bond length (r) of 1.79 Å and an Fe-N-O bond angle (β) of 136°-137°. This bond length is somewhat longer than those observed in most model compounds but may be due, at least in part, to hydrogen bond formation with the distal His residue. Bond elongation was also observed in a geometry optimized Fe(II)(NO)(imidazole) complex hydrogen bonded to an imidazole residue, in which we find r = 1.76-1.78 Å and β = 137°-138°. The computed bond angles are close to the canonical ∼140° value found in most model systems. Highly bent Fe-N-O bond angles or very long Fe-NO bond lengths seem unlikely to occur in proteins, due to their high energies. We also investigated the molecular orbitals and spin densities in each of the six coordinate systems investigated and found the orbitals and spin densities to be generally similar those described previously for five coordinate systems. Taken together, these results show that Mössbauer quadrupole splittings and isomer shifts, in addition to electron paramagnetic resonance hyperfine coupling constants, can now be calculated for nitrosyl heme systems with relatively good accuracy and that the results so obtained can be used to determine Fe-N-O geometries in metalloproteins. The Z-surface approach is thus applicable to both diamagnetic (CO) and paramagnetic (NO) heme proteins with in both cases the metal-ligand binding geometries found in the proteins being very close to those seen in model systems.
UR - http://www.scopus.com/inward/record.url?scp=0345802951&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=0345802951&partnerID=8YFLogxK
U2 - 10.1021/ja030340v
DO - 10.1021/ja030340v
M3 - Article
C2 - 14692781
AN - SCOPUS:0345802951
SN - 0002-7863
VL - 125
SP - 16387
EP - 16396
JO - Journal of the American Chemical Society
JF - Journal of the American Chemical Society
IS - 52
ER -