What Can Be Learned from Nuclear Resonance Vibrational Spectroscopy: Vibrational Dynamics and Hemes.

Nuclear resonance vibrational spectroscopy (NRVS; also known as nuclear inelastic scattering, NIS) is a synchrotron-based method that reveals the full spectrum of vibrational dynamics for Mössbauer nuclei. Another major advantage, in addition to its completeness (no arbitrary optical selection rules), is the unique selectivity of NRVS. The basics of this recently developed technique are first introduced with descriptions of the experimental requirements and data analysis including the details of mode assignments. We discuss the use of NRVS to probe 57Fe at the center of heme and heme protein derivatives yielding the vibrational density of states for the iron. The application to derivatives with diatomic ligands (O2, NO, CO, CN-) shows the strong capabilities of identifying mode character. The availability of the complete vibrational spectrum of iron allows the identification of modes not available by other techniques. This permits the correlation of frequency with other physical properties. A significant example is the correlation we find between the Fe-Im stretch in six-coordinate Fe(XO) hemes and the trans Fe-N(Im) bond distance, not possible previously. NRVS also provides uniquely quantitative insight into the dynamics of the iron. For example, it provides a model-independent means of characterizing the strength of iron coordination. Prediction of the temperature-dependent mean-squared displacement from NRVS measurements yields a vibrational "baseline" for Fe dynamics that can be compared with results from techniques that probe longer time scales to yield quantitative insights into additional dynamical processes.

How Does a Heme Carbene Differ from Diatomic Ligated (NO, CO, and CN-) Analogues in the Axial Bond?

Compared to well studied diatomic ligands (NO, CN-, CO), the axial bonds of carbene hemes is much less known although its significance in biological chemistry. The unusually large quadrupole splitting (Δ EQ = +2.2 mm·s-1) and asymmetric parameter (η = 0.9) of the five-coordinate heme carbene [Fe(TTP)(CCl2)], which is the largest among all known low spin ferrohemes, has driven investigations by means of Mössbauer effect Nuclear Resonance Vibrational Spectroscopy (NRVS). Three distinct measurements on one single crystal (two in-plane and one out-of-plane) have demonstrated comprehensive vibrational structures including stretch (429) and bending modes (472 cm-1) of the axial Fe-CCl2, and revealed iron vibrational anisotropy in three orthogonal directions for the first time. Frontier orbital analysis especially comparisons with diatomic analogues (NO, CN-, CO) suggest that CCl2, similar to NO, has led to strong but anisotropic π bonding in a ligand-based "4C"-coordinate which induced the vibrational anisotropies and very large Mössbauer parameters. This is contrasted to CN- and CO complexes which possess a porphyrin-based "4N"-coordinate electronic and vibrational structures due to inherent on-axis linear ligation.

3D Motions of Iron in Six-Coordinate {FeNO}(7) Hemes by Nuclear Resonance Vibration Spectroscopy.

The vibrational spectrum of a six-coordinate nitrosyl iron porphyrinate, monoclinic [Fe(TpFPP)(1-MeIm)(NO)] (TpFPP=tetra-para-fluorophenylporphyrin; 1-MeIm=1-methylimidazole), has been studied by oriented single-crystal nuclear resonance vibrational spectroscopy (NRVS). The crystal was oriented to give spectra perpendicular to the porphyrin plane and two in-plane spectra perpendicular or parallel to the projection of the FeNO plane. These enable assignment of the FeNO bending and stretching modes. The measurements reveal that the two in-plane spectra have substantial differences that result from the strongly bonded axial NO ligand. The direction of the in-plane iron motion is found to be largely parallel and perpendicular to the projection of the bent FeNO on the porphyrin plane. The out-of-plane Fe-N-O stretching and bending modes are strongly mixed with each other, as well as with porphyrin ligand modes. The stretch is mixed with v50 as was also observed for dioxygen complexes. The frequency of the assigned stretching mode of eight Fe-X-O (X=N, C, and O) complexes is correlated with the Fe-XO bond lengths. The nature of highest frequency band at ≈560 cm(-1) has also been examined in two additional new derivatives. Previously assigned as the Fe-NO stretch (by resonance Raman), it is better described as the bend, as the motion of the central nitrogen atom of the FeNO group is very large. There is significant mixing of this mode. The results emphasize the importance of mode mixing; the extent of mixing must be related to the peripheral phenyl substituents.

Recent advances in biosynthetic modeling of nitric oxide reductases and insights gained from nuclear resonance vibrational and other spectroscopic studies.

This Forum Article focuses on recent advances in structural and spectroscopic studies of biosynthetic models of nitric oxide reductases (NORs). NORs are complex metalloenzymes found in the denitrification pathway of Earth's nitrogen cycle where they catalyze the proton-dependent two-electron reduction of nitric oxide (NO) to nitrous oxide (N2O). While much progress has been made in biochemical and biophysical studies of native NORs and their variants, a clear mechanistic understanding of this important metalloenzyme related to its function is still elusive. We report herein UV-vis and nuclear resonance vibrational spectroscopy (NRVS) studies of mononitrosylated intermediates of the NOR reaction of a biosynthetic model. The ability to selectively substitute metals at either heme or nonheme metal sites allows the introduction of independent (57)Fe probe atoms at either site, as well as allowing the preparation of analogues of stable reaction intermediates by replacing either metal with a redox inactive metal. Together with previous structural and spectroscopic results, we summarize insights gained from studying these biosynthetic models toward understanding structural features responsible for the NOR activity and its mechanism. The outlook on NOR modeling is also discussed, with an emphasis on the design of models capable of catalytic turnovers designed based on close mimics of the secondary coordination sphere of native NORs.

Deep proton tunneling in the electronically adiabatic and non-adiabatic limits: comparison of the quantum and classical treatment of donor-acceptor motion in a protein environment.

Analytical models describing the temperature dependence of the deep tunneling rate, useful for proton, hydrogen, or hydride transfer in proteins, are developed and compared. Electronically adiabatic and non-adiabatic expressions are presented where the donor-acceptor (D-A) motion is treated either as a quantized vibration or as a classical "gating" distribution. We stress the importance of fitting experimental data on an absolute scale in the electronically adiabatic limit, which normally applies to these reactions, and find that vibrationally enhanced deep tunneling takes place on sub-ns timescales at room temperature for typical H-bonding distances. As noted previously, a small room temperature kinetic isotope effect (KIE) does not eliminate deep tunneling as a major transport channel. The quantum approach focuses on the vibrational sub-space composed of the D-A and hydrogen atom motions, where hydrogen bonding and protein restoring forces quantize the D-A vibration. A Duschinsky rotation is mandated between the normal modes of the reactant and product states and the rotation angle depends on the tunneling particle mass. This tunnel-mass dependent rotation contributes substantially to the KIE and its temperature dependence. The effect of the Duschinsky rotation is solved exactly to find the rate in the electronically non-adiabatic limit and compared to the Born-Oppenheimer (B-O) approximation approach. The B-O approximation is employed to find the rate in the electronically adiabatic limit, where we explore both harmonic and quartic double-well potentials for the hydrogen atom bound states. Both the electronically adiabatic and non-adiabatic rates are found to diverge at high temperature unless the proton coupling includes the often neglected quadratic term in the D-A displacement from equilibrium. A new expression is presented for the electronically adiabatic tunnel rate in the classical limit for D-A motion that should be useful to experimentalists working near room temperature. This expression also holds when a broad protein conformational distribution of D-A equilibrium distances dominates the spread of the D-A vibrational wavefunction.

Comprehensive Fe-ligand vibration identification in {FeNO}6 hemes.

Oriented single-crystal nuclear resonance vibrational spectroscopy (NRVS) has been used to obtain all iron vibrations in two {FeNO}(6) porphyrinate complexes, five-coordinate [Fe(OEP)(NO)]ClO4 and six-coordinate [Fe(OEP)(2-MeHIm)(NO)]ClO4. A new crystal structure was required for measurements of [Fe(OEP)(2-MeHIm)(NO)]ClO4, and the new structure is reported herein. Single crystals of both complexes were oriented to be either parallel or perpendicular to the porphyrin plane and/or axial imidazole ligand plane. Thus, the FeNO bending and stretching modes can now be unambiguously assigned; the pattern of shifts in frequency as a function of coordination number can also be determined. The pattern is quite distinct from those found for CO or {FeNO}(7) heme species. This is the result of unchanging Fe-N(NO) bonding interactions in the {FeNO}(6) species, in distinct contrast to the other diatomic ligand species. DFT calculations were also used to obtain detailed predictions of vibrational modes. Predictions were consistent with the intensity and character found in the experimental spectra. The NRVS data allow the assignment and observation of the challenging to obtain Fe-Im stretch in six-coordinate heme derivatives. NRVS data for this and related six-coordinate hemes with the diatomic ligands CO, NO, and O2 reveal a strong correlation between the Fe-Im stretch and Fe-N(Im) bond distance that is detailed for the first time.

The diagnostic vibrational signature of pentacoordination in heme carbonyls.

Heme-carbonyl complexes are widely exploited for the insight they provide into the structural basis of function in heme-based proteins, by revealing the nature of their bonded and nonbonded interactions with the protein. This report presents two novel results which clearly establish a FeCO vibrational signature for crystallographically verified pentacoordination. First, anisotropy in the NRVS density of states for ν(Fe-C) and δ(FeCO) in oriented single crystals of [Fe(OEP)(CO)] clearly reveals that the Fe-C stretch occurs at higher frequency than the FeCO bend and considerably higher than any previously reported heme carbonyl. Second, DFT calculations on a series of heme carbonyls reveal that the frequency crossover occurs near the weak trans O atom donor, furan. As ν(Fe-C) occurs at lower frequencies than δ(FeCO) in all heme protein carbonyls reported to date, the results reported herein suggest that they are all hexacoordinate.

Anisotropic iron motion in nitrosyl iron porphyrinates: natural and synthetic hemes.

The vibrational spectra of two five-coordinate nitrosyl iron porphyrinates, [Fe(OEP)(NO)] (OEP = dianion of 2,3,7,8,12,13,17,18-octaethylporphyrin) and [Fe(DPIX)(NO)] (DPIX = deuteroporphyrin IX), have been studied by oriented single-crystal nuclear resonance vibrational spectroscopy. Single crystals (both are in the triclinic crystal system) were oriented to give vibrational spectra perpendicular to the porphyrin plane. Additionally, two orthogonal in-plane measurements that were also either perpendicular or parallel to the projection of the FeNO plane onto the porphyrin plane yield the complete set of vibrations with iron motion. In addition to cleanly enabling the assignment of the FeNO bending and stretching modes, the measurements reveal that the two in-plane spectra from the parallel and perpendicular in-plane directions for both compounds have substantial differences. The assignment of these in-plane vibrations were aided by density functional theory predictions. The differences in the two in-plane directions result from the strongly bonded axial NO ligand. The direction of the in-plane iron motion is thus found to be largely parallel and perpendicular to the projection of the FeNO plane on the porphyrin plane. These axial ligand effects on the in-plane iron motion are related to the strength of the axial ligand-to-iron bond.

Spectroscopic and computational study of a nonheme iron nitrosyl center in a biosynthetic model of nitric oxide reductase.

A major barrier to understanding the mechanism of nitric oxide reductases (NORs) is the lack of a selective probe of NO binding to the nonheme FeB center. By replacing the heme in a biosynthetic model of NORs, which structurally and functionally mimics NORs, with isostructural ZnPP, the electronic structure and functional properties of the FeB nitrosyl complex was probed. This approach allowed observation of the first S=3/2 nonheme {FeNO}(7) complex in a protein-based model system of NOR. Detailed spectroscopic and computational studies show that the electronic state of the {FeNO}(7) complex is best described as a high spin ferrous iron (S=2) antiferromagnetically coupled to an NO radical (S=1/2) [Fe(2+)-NO(.)]. The radical nature of the FeB -bound NO would facilitate N-N bond formation by radical coupling with the heme-bound NO. This finding, therefore, supports the proposed trans mechanism of NO reduction by NORs.

Structural changes that occur upon photolysis of the Fe(II)(a3)-CO complex in the cytochrome ba(3)-oxidase of Thermus thermophilus: a combined X-ray crystallographic and infrared spectral study demonstrates CO binding to Cu(B).

The purpose of the work was to provide a crystallographic demonstration of the venerable idea that CO photolyzed from ferrous heme-a(3) moves to the nearby cuprous ion in the cytochrome c oxidases. Crystal structures of CO-bound cytochrome ba(3)-oxidase from Thermus thermophilus, determined at ~2.8-3.2Å resolution, reveal a Fe-C distance of ~2.0Å, a Cu-O distance of 2.4Å and a Fe-C-O angle of ~126°. Upon photodissociation at 100K, X-ray structures indicate loss of Fe(a3)-CO and appearance of Cu(B)-CO having a Cu-C distance of ~1.9Å and an O-Fe distance of ~2.3Å. Absolute FTIR spectra recorded from single crystals of reduced ba(3)-CO that had not been exposed to X-ray radiation, showed several peaks around 1975cm(-1); after photolysis at 100K, the absolute FTIR spectra also showed a significant peak at 2050cm(-1). Analysis of the 'light' minus 'dark' difference spectra showed four very sharp CO stretching bands at 1970cm(-1), 1977cm(-1), 1981cm(-1), and 1985cm(-1), previously assigned to the Fe(a3)-CO complex, and a significantly broader CO stretching band centered at ~2050cm(-1), previously assigned to the CO stretching frequency of Cu(B) bound CO. As expected for light propagating along the tetragonal axis of the P4(3)2(1)2 space group, the single crystal spectra exhibit negligible dichroism. Absolute FTIR spectrometry of a CO-laden ba(3) crystal, exposed to an amount of X-ray radiation required to obtain structural data sets before FTIR characterization, showed a significant signal due to photogenerated CO(2) at 2337cm(-1) and one from traces of CO at 2133cm(-1); while bands associated with CO bound to either Fe(a3) or to Cu(B) in "light" minus "dark" FTIR difference spectra shifted and broadened in response to X-ray exposure. In spite of considerable radiation damage to the crystals, both X-ray analysis at 2.8 and 3.2Å and FTIR spectra support the long-held position that photolysis of Fe(a3)-CO in cytochrome c oxidases leads to significant trapping of the CO on the Cu(B) atom; Fe(a3) and Cu(B) ligation, at the resolutions reported here, are otherwise unaltered.

Effects of imidazole deprotonation on vibrational spectra of high-spin iron(II) porphyrinates.

The effects of the deprotonation of coordinated imidazole on the vibrational dynamics of five-coordinate high-spin iron(II) porphyrinates have been investigated using nuclear resonance vibrational spectroscopy. Two complexes have been studied in detail with both powder and oriented single-crystal measurements. Changes in the vibrational spectra are clearly related to structural differences in the molecular structures that occur when imidazole is deprotonated. Most modes involving the simultaneous motion of iron and imidazolate are unresolved, but the one mode that is resolved is found at higher frequency in the imidazolates. These out-of-plane results are in accord with earlier resonance Raman studies of heme proteins. We also show the imidazole vs imidazolate differences in the in-plane vibrations that are not accessible to resonance Raman studies. The in-plane vibrations are at lower frequency in the imidazolate derivatives; the doming mode shifts are inconclusive. The stiffness, an experimentally determined force constant that averages the vibrational details to quantify the nearest-neighbor interactions, confirms that deprotonation inverts the relative strengths of axial and equatorial coordination.

Probing heme vibrational anisotropy: an imidazole orientation effect?

The complete iron vibrational spectrum of the five-coordinate high-spin complex [Fe(OEP)(2-MeHIm)], where OEP = octaethylporphyrinato and 2-MeHIm = 2-methylimidazole, has been obtained by oriented single-crystal nuclear resonance vibrational spectroscopy (NRVS) data. Measurements have been made in three orthogonal directions, which provides quantitative information for all iron motion. These experimental data, buttressed by density functional theory (DFT) calculations, have been used to define the effects of the axial ligand orientation. Although the axial imidazole removes the degeneracy in the in-plane vibrations, the imidazole orientation does not appear to control the direction of the in-plane iron motion. This is in contrast to the effect of the imidazolate ligand, as defined by DFT calculations, which does have substantial effects on the direction of the in-plane iron motion. The axial NO ligand has been found to have the strongest orientational effect (Angew. Chem., Int. Ed., 2010, 49, 4400). Thus the strength of the directional properties are in the order NO > imidazolate > imidazole, consistent with the varying strength of the Fe-ligand bond.

Infrared protein crystallography.

We consider the application of infrared spectroscopy to protein crystals, with particular emphasis on exploiting molecular orientation through polarization measurements on oriented single crystals. Infrared microscopes enable transmission measurements on individual crystals using either thermal or nonthermal sources, and can accommodate flow cells, used to measure spectral changes induced by exposure to soluble ligands, and cryostreams, used for measurements of flash-cooled crystals. Comparison of unpolarized infrared measurements on crystals and solutions probes the effects of crystallization and can enhance the value of the structural models refined from X-ray diffraction data by establishing solution conditions under which they are most relevant. Results on several proteins are consistent with similar equilibrium conformational distributions in crystal and solutions. However, the rates of conformational change are often perturbed. Infrared measurements also detect products generated by X-ray exposure, including CO(2). Crystals with favorable symmetry exhibit infrared dichroism that enhances the synergy with X-ray crystallography. Polarized infrared measurements on crystals can distinguish spectral contributions from chemically similar sites, identify hydrogen bonding partners, and, in opportune situations, determine three-dimensional orientations of molecular groups. This article is part of a Special Issue entitled: Protein Structure and Function in the Crystalline State.

Spectroscopic characterization of (57)Fe-enriched cytochrome c.

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

Vibrational probes and determinants of the S = 0 ⇌ S = 2 spin crossover in five-coordinate [Fe(TPP)(CN)]-.

The low-frequency vibrational characterization of the spin-crossover complex, five-coordinate cyano(tetraphenylporphyrinato)iron(II), [Fe(TPP)(CN)](-), is reported. Nuclear resonance vibrational spectroscopy has been used to measure all low-frequency vibrations involving iron at several temperatures; this yields vibrational spectra of both the low- (S = 0) and high-spin (S = 2) states. Multitemperature oriented single-crystal measurements facilitate assignments of the vibrational character of all modes and are consistent with the DFT-predicted spectra. The availability of the entire iron vibrational spectrum allows for the complete correlation of the modes between the two spin states. These data demonstrate that not only do the frequencies of the vibrations shift to lower values for the high-spin species as would be expected owing to the weaker bonds in the high-spin state, but also the mixing of iron modes with ligand modes changes substantially. Diagrams illustrating the changing character of the modes and their correlation are given. The reduced iron-ligand frequencies are the primary factor in the entropic stabilization of the high-spin state responsible for the spin crossover.

Nuclear resonance vibrational spectra of five-coordinate imidazole-ligated iron(II) porphyrinates.

Nuclear resonance vibrational spectra have been obtained for six five-coordinate imidazole-ligated iron(II) porphyrinates, [Fe(Por)(L)] (Por = tetraphenylporphyrinate, octaethylporphyrinate, tetratolylporphyrinate, or protoporphyrinate IX and L = 2-methylimidazole or 1,2-dimethylimidazole). Measurements have been made on both powder and oriented crystal samples. The spectra are dominated by strong signals around 200-300 cm(-1). Although the in-plane and out-of-plane vibrations are seriously overlapped, oriented crystal spectra allow their deconvolution. Thus, oriented crystal experimental data, along with density functional theory (DFT) calculations, enable the assignment of key vibrations in the spectra. Molecular dynamics are also discussed. The nature of the Fe-N(Im) vibrations has been elaborated further than was possible from resonance Raman studies. Our study suggests that the Fe motions are coupled with the porphyrin core and peripheral groups motions. Both peripheral groups and their conformations have significant influence on the vibrational spectra (position and shape).

New perspectives on iron-ligand vibrations of oxyheme complexes.

We report our studies of the vibrational dynamics of iron for three imidazole-ligated oxyheme derivatives that mimic the active sites of histidine-ligated heme proteins complexed with dioxygen. The experimental vibrational data are obtained from nuclear resonance vibrational spectroscopy (NRVS) measurements conducted on both powder samples and oriented single crystals, and which includes several in-plane (ip) and out-of-plane (oop) measurements. Vibrational spectral assignments have been made through a combination of the oriented sample spectra and predictions based on density functional theory (DFT) calculations. The two Fe-O(2) modes that have been previously observed by resonance Raman spectroscopy in heme proteins are clearly shown to be very strongly mixed and are not simply either a bending or stretching mode. In addition, a third Fe-O(2) mode, not previously reported, has been identified. The long-sought Fe-Im stretch, not observed in resonance Raman spectra, has been identified and compared with the frequencies observed for the analogous CO and NO species. The studies also suggest that the in-plane iron motion is anisotropic and is controlled by the orientation of the Fe-O(2) group and not sensitive to the in-plane Fe-N(p) bonds and/or imidazole orientations.

Spectroscopic identification of reactive porphyrin motions.

Nuclear resonance vibrational spectroscopy (NRVS) reveals the vibrational dynamics of a Mössbauer probe nucleus. Here, (57)Fe NRVS measurements yield the complete spectrum of Fe vibrations in halide complexes of iron porphyrins. Iron porphine serves as a useful symmetric model for the more complex spectrum of asymmetric heme molecules that contribute to numerous essential biological processes. Quantitative comparison with the vibrational density of states (VDOS) predicted for the Fe atom by density functional theory calculations unambiguously identifies the correct sextet ground state in each case. These experimentally authenticated calculations then provide detailed normal mode descriptions for each observed vibration. All Fe-ligand vibrations are clearly identified despite the high symmetry of the Fe environment. Low frequency molecular distortions and acoustic lattice modes also contribute to the experimental signal. Correlation matrices compare vibrations between different molecules and yield a detailed picture of how heme vibrations evolve in response to (a) halide binding and (b) asymmetric placement of porphyrin side chains. The side chains strongly influence the energetics of heme doming motions that control Fe reactivity, which are easily observed in the experimental signal.

Electronic structure and dynamics of nitrosyl porphyrins.

Nitric oxide (NO) is a signaling molecule employed to regulate essential physiological processes. Thus, there is great interest in understanding the interaction of NO with heme, which is found at the active site of many proteins that recognize NO, as well as those involved in its creation and elimination. We summarize what we have learned from investigations of the structure, vibrational properties, and conformational dynamics of NO complexes with ferrous porphyrins, as well as computational investigations in support of these experimental studies. Multitemperature crystallographic data reveal variations in the orientational disorder of the nitrosyl ligand. In some cases, equilibria among NO orientations can be analyzed using the van't Hoff relationship and the free energy and enthalpy of the solid-state transitions evaluated experimentally. Density functional theory (DFT) calculations predict that intrinsic barriers to torsional rotation are smaller than thermal energies at physiological temperatures, and the coincidence of observed NO orientations with minima in molecular mechanics potentials indicates that nonbonded interactions with other chemical groups control the conformational freedom of the bound NO. In favorable cases, reduced disorder at low temperatures exposes subtle structural features including off-axis tilting of the Fe-NO bond and anisotropy of the equatorial Fe-N bonds. We also present the results of nuclear resonance vibrational spectroscopy measurements on oriented single crystals of [Fe(TPP)(NO)] and [Fe(TPP)(1-MeIm)(NO)]. These describe the anisotropic vibrational motion of iron in five- and six-coordinate heme-NO complexes and reveal vibrations of all Fe-ligand bonds as well as low-frequency molecular distortions associated with the doming of the heme upon ligand binding. A quantitative comparison with predicted frequencies, amplitudes, and directions facilitates identification of the vibrational modes but also suggests that commonly used DFT functionals are not fully successful at capturing the trans interaction between the axial NO and imidazole ligands. This supports previous conclusions that heme-NO complexes exhibit an unusual degree of variability with respect to the computational method, and we speculate that this variability hints at a genuine electronic instability that a protein can exploit to tune its reactivity. We anticipate that ongoing characterization of heme-NO complexes will deepen our understanding of their structure, dynamics, and reactivity.

Oriented single-crystal nuclear resonance vibrational spectroscopy of [Fe(TPP)(MI)(NO)]: quantitative assessment of the trans effect of NO.

This paper presents oriented single-crystal Nuclear Resonance Vibrational Spectroscopy (NRVS) data for the six-coordinate (6C) ferrous heme-nitrosyl model complex [(57)Fe(TPP)(MI)(NO)] (1; TPP(2-) = tetraphenylporphyrin dianion; MI = 1-methylimidazole). The availability of these data enables for the first time the detailed simulation of the complete NRVS data, including the porphyrin-based vibrations, of a 6C ferrous heme-nitrosyl, using our quantum chemistry centered normal coordinate analysis (QCC-NCA). Importantly, the Fe-NO stretch is split by interaction with a porphyrin-based vibration into two features, observed at 437 and 472 cm(-1). The 437 cm(-1) feature is strongly out-of-plane (oop) polarized and shows a (15)N(18)O isotope shift of 8 cm(-1) and is therefore assigned to nu(Fe-NO). The admixture of Fe-N-O bending character is small. Main contributions to the Fe-N-O bend are observed in the 520-580 cm(-1) region, distributed over a number of in-plane (ip) polarized porphyrin-based vibrations. The main component, assigned to delta(ip)(Fe-N-O), is identified with the feature at 563 cm(-1). The Fe-N-O bend also shows strong mixing with the Fe-NO stretching internal coordinate, as evidenced by the oop NRVS intensity in the 520-580 cm(-1) region. Very accurate normal mode descriptions of nu(Fe-NO) and delta(ip)(Fe-N-O) have been obtained in this study. These results contradict previous interpretations of the vibrational spectra of 6C ferrous heme-nitrosyls where the higher energy feature at approximately 550 cm(-1) had usually been associated with nu(Fe-NO). Furthermore, these results provide key insight into NO binding to ferrous heme active sites in globins and other heme proteins, in particular with respect to (a) the effect of hydrogen bonding to the coordinated NO and (b) changes in heme dynamics upon NO coordination. [Fe(TPP)(MI)(NO)] constitutes an excellent model system for ferrous NO adducts of myoglobin (Mb) mutants where the distal histidine (His64) has been removed. Comparison to the reported vibrational data for wild-type (wt) Mb-NO then shows that the effect of H bonding to the coordinated NO is weak and mostly leads to a polarization of the pi/pi* orbitals of bound NO. In addition, the observation that delta(ip)(Fe-N-O) does not correlate well with nu(N-O) can be traced back to the very mixed nature of this mode. The Fe-N(imidazole) stretching frequency is observed at 149 cm(-1) in [Fe(TPP)(MI)(NO)], and spectral changes upon NO binding to five-coordinate ferrous heme active sites are discussed. The obtained high-quality force constants for the Fe-NO and N-O bonds of 2.57 and 11.55 mdyn/A can further be compared to those of corresponding 5C species, which allows for a quantitative analysis of the sigma trans interaction between the proximal imidazole (His) ligand and NO. This is key for the activation of the NO sensor soluble guanylate cyclase. Finally, DFT methods are calibrated against the experimentally determined vibrational properties of the Fe-N-O subunit in 1. DFT is in fact incapable of reproducing the vibrational energies and normal mode descriptions of the Fe-N-O unit well, and thus, DFT-based predictions of changes in vibrational properties upon heme modification or other perturbations of these 6C complexes have to be treated with caution.