Quantitative vibrational dynamics of iron in carbonyl porphyrins

The Biologically Relevant Coordination Chemistry of Iron and Nitric Oxide: Electronic Structure and Reactivity

Nitric oxide (NO) is an important signaling molecule that is involved in a wide range of physiological and pathological events in biology. Metal coordination chemistry, especially with iron, is at the heart of many biological transformations involving NO. A series of heme proteins, nitric oxide synthases (NOS), soluble guanylate cyclase (sGC), and nitrophorins, are responsible for the biosynthesis, sensing, and transport of NO. Alternatively, NO can be generated from nitrite by heme- and copper-containing nitrite reductases (NIRs). The NO-bearing small molecules such as nitrosothiols and dinitrosyl iron complexes (DNICs) can serve as an alternative vehicle for NO storage and transport. Once NO is formed, the rich reaction chemistry of NO leads to a wide variety of biological activities including reduction of NO by heme or non-heme iron-containing NO reductases and protein post-translational modifications by DNICs. Much of our understanding of the reactivity of metal sites in biology with NO and the mechanisms of these transformations has come from the elucidation of the geometric and electronic structures and chemical reactivity of synthetic model systems, in synergy with biochemical and biophysical studies on the relevant proteins themselves. This review focuses on recent advancements from studies on proteins and model complexes that not only have improved our understanding of the biological roles of NO but also have provided foundations for biomedical research and for bio-inspired catalyst design in energy science.

SciPhon: a data analysis software for nuclear resonant inelastic X-ray scattering with applications to Fe, Kr, Sn, Eu and Dy

The synchrotron radiation technique of nuclear resonant inelastic X-ray scattering (NRIXS), also known as nuclear resonance vibrational spectroscopy or nuclear inelastic scattering, provides a wealth of information on the vibrational properties of solids. It has found applications in studies of lattice dynamics and elasticity, superconductivity, heme biochemistry, seismology, isotope geochemistry and many other fields. It involves probing the vibrational modes of solids by using the nuclear resonance of Mössbauer isotopes such as ⁵⁷Fe, ⁸³Kr, ¹¹⁹Sn, ¹⁵¹Eu and ¹⁶¹Dy. After data reduction, it provides the partial phonon density of states of the Mössbauer isotope that is investigated, as well as many other derived quantities such as the mean force constant of the chemical bonds and the Debye velocity. The data reduction is, however, not straightforward and involves removal of the elastic peak, normalization and Fourier-Log transformation. Furthermore, some of the quantities derived are highly sensitive to details in the baseline correction. A software package and several novel procedures to streamline and hopefully improve the reduction of the NRIXS data generated at sector 3ID of the Advanced Photon Source have been developed. The graphical user interface software is named SciPhon and runs as a Mathematica package. It is easily portable to other platforms and can be easily adapted for reducing data generated at other beamlines. Several tests and comparisons are presented that demonstrate the usefulness of this software, whose results have already been used in several publications. Here, the SciPhon software is used to reduce Kr, Sn, Eu and Dy NRIXS data, and potential implications for interpreting natural isotopic variations in those systems are discussed.

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 ⁵⁷Fe 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 (O₂, 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.

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 (N₂O). 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 ⁵⁷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.

New insights into nitric oxide reductases (NORs) are obtained. Using nuclear resonance vibrational spectroscopy, we probe both iron atoms in mononitrosylated intermediates of the NOR reaction in a biosynthetic protein model that reveal new insights into the structural and electronic features responsible for the NOR activity and its likely mechanism.

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.

Quantitative vibrational dynamics of the metal site in a tin porphyrin: an IR, NRVS, and DFT study

We used a newer, synchrotron-based, spectroscopic technique (nuclear resonance vibrational spectroscopy, NRVS) in combination with a more traditional one (infrared absorption, IR) to obtain a complete, quantitative picture of the metal center vibrational dynamics in a six-coordinated tin porphyrin. From the NRVS (119)Sn site-selectivity and the sensitivity of the IR signal to (112)Sn/(119)Sn isotope substitution, we identified the frequency of the antisymmetric stretching of the axial bonds (290 cm(-1)) and all the other vibrations involving Sn. Experimentally authenticated density functional theory (DFT) calculations aid the data interpretation by providing detailed normal mode descriptions for each observed vibration. These results may represent a starting point toward the characterization of the local vibrational dynamics of the metallic site in tin porphyrins and compounds with related structures. The quantitative complementariness between IR, NRVS, and DFT is emphasized.

Observation of the Fe-CN and Fe-CO vibrations in the active site of [NiFe] hydrogenase by nuclear resonance vibrational spectroscopy

Nuclear inelastic scattering of (57)Fe labeled [NiFe] hydrogenase is shown to give information on different states of the enzyme. It was thus possible to detect and assign Fe-CO and Fe-CN bending and stretching vibrations of the active site outside the spectral range of the Fe-S cluster normal modes.

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 spectroscopy (NRVS) of rubredoxin and MoFe protein crystals

We have applied ⁵⁷Fe nuclear resonance vibrational spectroscopy (NRVS) for the first time to study the dynamics of Fe centers in Fe-S protein crystals, including oxidized wild type rubredoxin crystals from Pyrococcus furiosus, and the MoFe protein of nitrogenase from Azotobacter vinelandii. Thanks to the NRVS selection rule, selectively probed vibrational modes have been observed in both oriented rubredoxin and MoFe protein crystals. The NRVS work was complemented by extended X-ray absorption fine structure spectroscopy (EXAFS) measurements on oxidized wild type rubredoxin crystals from Pyrococcus furiosus. The EXAFS spectra revealed the Fe-S bond length difference in oxidized Pf Rd protein, which is qualitatively consistent with the X-ray crystal structure.

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.

Cis-dicarbonyl binding at cobalt and iron porphyrins with saddle-shape conformation

Diatomic molecules attached to complexed iron or cobalt centres are important in many biological processes. In natural systems, metallotetrapyrrole units carry respiratory gases or provide sensing and catalytic functions. Conceiving synthetic model systems strongly helps to determine the pertinent chemical foundations for such processes, with recent work highlighting the importance of the prosthetic groups' conformational flexibility as an intricate variable affecting their functional properties. Here, we present simple model systems to investigate, at the single molecule level, the interaction of carbon monoxide with saddle-shaped iron- and cobalt-porphyrin conformers, which have been stabilized as two-dimensional arrays on well-defined surfaces. Using scanning tunnelling microscopy we identified a novel bonding scheme expressed in tilted monocarbonyl and cis-dicarbonyl configurations at the functional metal-macrocycle unit. Modelling with density functional theory revealed that the weakly bonded diatomic carbonyl adduct can effectively bridge specific pyrrole groups with the metal atom as a result of the pronounced saddle-shape conformation of the porphyrin cage.

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.

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.

Nuclear resonance vibrational spectroscopy applied to [Fe(OEP)(NO)]: the vibrational assignments of five-coordinate ferrous heme-nitrosyls and implications for electronic structure

This study presents Nuclear Resonance Vibrational Spectroscopy (NRVS) data on the five-coordinate (5C) ferrous heme-nitrosyl complex [Fe(OEP)(NO)] (1, OEP(2-) = octaethylporphyrinato dianion) and the corresponding (15)N(18)O labeled complex. The obtained spectra identify two isotope sensitive features at 522 and 388 cm(-1), which shift to 508 and 381 cm(-1), respectively, upon isotope labeling. These features are assigned to the Fe-NO stretch nu(Fe-NO) and the in-plane Fe-N-O bending mode delta(ip)(Fe-N-O), the latter has been unambiguously assigned for the first time for 1. The obtained NRVS data were simulated using our quantum chemistry centered normal coordinate analysis (QCC-NCA). Since complex 1 can potentially exist in 12 different conformations involving the FeNO and peripheral ethyl orientations, extended density functional theory (DFT) calculations and QCC-NCA simulations were performed to determine how these conformations affect the NRVS properties of [Fe(OEP)NO]. These results show that the properties and force constants of the FeNO unit are hardly affected by the conformational changes involving the ethyl substituents. On the other hand, the NRVS-active porphyrin-based vibrations around 340-360, 300-320, and 250-270 cm(-1) are sensitive to the conformational changes. The spectroscopic changes observed in these regions are due to selective mechanical couplings of one component of E(u)-type (in ideal D(4h) symmetry) porphyrin-based vibrations with the in-plane Fe-N-O bending mode. This leads to the observed variations in Fe(OEP) core mode energies and NRVS intensities without affecting the properties of the FeNO unit. The QCC-NCA simulated NRVS spectra of 1 show excellent agreement with experiment, and indicate that conformer F is likely present in the samples of this complex investigated here. The observed porphyrin-based vibrations in the NRVS spectra of 1 are also assigned based on the QCC-NCA results. The obtained force constants of the Fe-NO and N-O bonds are 2.83-2.94 (based on the DFT functional applied) and about 12.15 mdyn/A, respectively. The electronic structures of 5C ferrous heme-nitrosyls in different model complexes are then analyzed, and variations in their properties based on different porphyrin substituents are explained. Finally, the shortcomings of different DFT functionals in describing the axial FeNO subunit in heme-nitrosyls are elucidated.

Vibrational dynamics of iron in cytochrome C

Nuclear resonance vibrational spectroscopy (NRVS) and Raman spectroscopy on (54)Fe- and (57)Fe-enriched cytochrome c (cyt c) identify multiple bands involving vibrations of the heme Fe. Comparison with predictions from Fe isotope shifts reveals that 70% of the NRVS signal in the 300-450 cm(-1) frequency range corresponds to vibrations resolved in Soret-enhanced Raman spectra. This frequency range dominates the "stiffness", an effective force constant determined by the Fe vibrational density of states (VDOS), which measures the strength of nearest-neighbor interactions with Fe. The stiffness of the low-spin Fe environment in both oxidation states of cyt c significantly exceeds that for the high-spin Fe in deoxymyoglobin, where the 200-300 cm(-1) frequency range dominates the VDOS. This situation is reflected in the shorter Fe-ligand bond lengths in the former with respect to the latter. The longer Fe-S(Met80) in oxidized cyt c with respect to reduced cyt c leads to a decrease in the stiffness of the iron environment upon oxidation. Comparison with NRVS measurements allows us to assess assignments for vibrational modes resolved in this region of the heme Raman spectrum. We consider the possibility that the 372 cm(-1) band in reduced cyt c involves the Fe-S(Met80) bond.

Quantum mechanical investigations of heme structure and vibrational spectra: effects of conformation, oxidation state, and electric field

Three conformers of a complex consisting of a Fe-porphyrin ring with axial imidazole groups FeP(Im)(2) that differed by orientation of axial imidazole groups were studied with quantum mechanical calculations in both reduced (neutral) and oxidized (cationic) states. In the reduced state, all three conformers correspond to local minima. On the other hand, in the oxidized state the conformer with the imidazole groups perpendicular to one another does not have a local minimum, which suggests that there could be changes in the structure of heme groups depending on oxidation state. The relative energy differences between conformers are small as well the differences of orbital energies. However, the populations of Fe 3d orbitals, and hence charge distributions, are predicted to change significantly between conformations. The orientation of heme groups can affect the kinetics of interheme electron transfer, so this electronic population redistribution will affect electron transfer kinetics. A comparison of calculated Raman spectra with measured surface-enhanced Raman spectra (Biju et al. Langmuir 2007, 23, 1333) shows excellent agreement for frequencies, but correlations of Raman intensities are less satisfactory. Possible explanations of the observed discrepancies could arise due to a problem in assignments of vibration bands in the experiment, a dependence of Raman spectra on heme complex on oxidation state, long-range protein structure, or membrane electric field are discussed. In particular, it was shown that the oxidation state does not dramatically alter the calculated Raman spectrum. However, the cell membrane electric field can significantly modify Raman spectra. An interpretation of the experimental oxidized cytochrome surface-enhanced Raman spectra is discussed.

Intermolecular dynamics in crystalline iron octaethylporphyrin (FeOEP)

The new technique of nuclear resonance vibrational spectroscopy (NRVS) has increased the range and quality of dynamical data from Fe-containing molecules that when combined with Raman and infrared spectroscopies impose stricter constraints on normal mode simulations, especially at lower frequencies. Going beyond the usual single molecule approximation, a classical normal-mode analysis that includes intermolecular coupling and the full crystalline symmetry is found to produce a better fit with fewer free parameters for the heme compound iron octaethylporphyrin (FeOEP), using NRVS data from polycrystalline material. Off-diagonal force constants were completely unnecessary, indicating that their role in previous single molecule fits was just to emulate intermolecular coupling. Sound velocities deduced from the calculated phonon dispersion curves are compared to NRVS measurements to further constrain the intermolecular force constants. The NRVS data by themselves are insufficient to rigorously determine all unknown force constants for molecules of this size, but the improved crystal model fit indicates the necessity of including intermolecular interactions for normal-mode analyses.

Reversible NO motion in crystalline [Fe(Porph)(1-MeIm)(NO)] derivatives

The synthesis, characterization, and X-ray structures of three low-spin (nitrosyl)iron(II) tetraarylporphyrinates, [Fe(TpXPP)(NO)(1-MeIm)], where X = F (in a triclinic and a monoclinic form) and OCH(3) are reported. All three molecules, at 100 K, have a single orientation of NO. These structures are the first examples of ordered NO's in [Fe(Porph)(NO)(1-MeIm)] complexes. The three new derivatives have similar structural features including a previously unnoted "bowing" of the N(NO)-Fe-N(Im) angle caused by a concerted tilting of the axial Fe-N(NO) and Fe-N(Im) bonds. Structural features such as the displacement of Fe out of the mean porphyrin plane toward NO, tilting of the Fe-N(NO) bond off the heme normal, and the asymmetry of the Fe-N(por) bonds further strengthen and confirm observations from earlier studies. The [Fe(TpXPP)(NO)(1-MeIm)] complexes were also studied at temperatures between 125 and 350 K to investigate temperature-dependent variations and trends in the coordination group geometry. At varying temperatures (above 150 K), all three derivatives display a second orientation of the NO ligand. The population and depopulation of this second orientation are thermally driven, with no apparent hysteresis. Crystal packing appears to be the significant feature in defining the order/disorder of the NO ligand. The length of the bond trans to NO, Fe-N(Im), was also found to be sensitive to temperature variation. The Fe-N(Im) bond length increases with increased temperature, whereas no other bonds change appreciably. The temperature-dependent Fe-N(Im) bond length change and cell volume changes are consistent with a "soft" Fe-N(Im) bond. Variable-temperature measurements show that the N-O stretching frequency changes with the Fe-N(Im) bond length. Temperature-dependent changes in the Fe-NIm bond length and N-O stretching frequency were also found to be completely reversible with no apparent hysteresis.

Explorations in metalloporphyrin stereochemistry, physical properties and beyond

A review of selected portions of our work in the area of porphyrin structure and physical characterization is presented. Topics covered include early work on periodic trends in first row transtion metalloporphyrins, a survey of electronic structure of iron derivatives including spin-state trends, ligand orientation effects and the elucidtion of unusual low-spin states for iron(II). A discussion of the different tlypes of high-spin iron(II) complexes and the effects of hydrogen bonding is given. A survey of nitric oxide (NO) derivatives is presented as well as a brief introduction into the use of nuclear resonance vibrational spectroscopy for the study of iron porphyrins and heme proteins.