An Experimental and Density Functional Theoretical Investigation of Iron-57 Mössbauer Quadrupole Splittings in Organometallic and Heme-Model Compounds: Applications to Carbonmonoxy-Heme Protein Structure

In Crystallo O2 Cleavage at a Preorganized Triiron Cluster

In Nature, the four-electron reduction of O2 is catalyzed at preorganized multimetallic active sites. These complex active sites often feature low-coordinate, redox-active metal centers precisely positioned to facilitate rapid O2 activation processes that obviate the generation of toxic, partially reduced oxygen species. Very few biomimetic constructs simultaneously recapitulate the complexity and reactivity of these biological cofactors. Herein, we report solid-state O2 activation at a triiron(II) active site templated by phosphinimide ligands. Insight into the structure of the O2 reduction intermediates was obtained via in crystallo O2 dosing experiments in conjunction with spectroscopic, structural, magnetic, and computational studies. These data support the in situ formation of an Fe2IIIFeIV-dioxo intermediate upon exposure to O2 that participates in oxygen atom and hydrogen atom transfer reactivity with exogenous substrates to furnish a stable FeIIFe2III-oxo species. Combined, these studies provide an extraordinary level of detail into the dynamics of bond-forming and -breaking processes operative at complex multimetallic active sites.

Protoporphyrin IX iron(II) revisited. An overview of the Mössbauer spectroscopic parameters of low-spin porphyrin iron(II) complexes

Mössbauer parameters of low-spin six-coordinate [Fe(II)(Por)L2] complexes (where Por is a synthetic porphyrin; L is a nitrogenous aliphatic, an aromatic base or a heterocyclic ligand, a P-bonding ligand, CO or CN) and low-spin [Fe(Por)LX] complexes (where L and X are different ligands) are reported. A known point charge calculation approach was extended to investigate how the axial ligands and the four porphyrinato-N atoms generate the observed quadrupole splittings (ΔEQ) for the complexes. Partial quadrupole splitting (p.q.s.) and partial chemical shifts (p.c.s.) values were derived for all the axial ligands, and porphyrins reported in the literature. The values for each porphyrin are different emphasising the importance/uniqueness of the [Fe(PPIX)] moiety, (which is ubiquitous in nature). This new analysis enabled the construction of figures relating p.c.s and p.q.s values. The relationships presented in the figures indicates that strong field ligands such as CO can, and do change the sign of the electric field gradient in the [Fe(II)(Por)L2] complexes. The limiting p.q.s. value a ligand can have and still form a six-coordinate low-spin [Fe(II)(Por)L2] complex is established. It is shown that the control the porphyrin ligands exert on the low-spin Fe(II) atom limits its bonding to a defined range of axial ligands; outside this range the spin state of the iron is unstable and five-coordinate high-spin complexes are favoured. Amongst many conclusions, it was found that oxygen cannot form a stable low-spin [Fe(II)(Por)L(O2)] complex and that oxy-haemoglobin is best described as an [Fe(III)(Por)L(O2-)] complex, the iron is ferric bound to the superoxide molecule.

Studies on the binding of CO to low-spin [Fe(II)(Por)L2] complexes: an aid to understanding the binding of CO to haemoglobin and myoglobin

The visible and Mössbauer spectra of [Fe(II)(Por)L₂] and [Fe(II)(Por)L(CO)] complexes (where Por = protoporphyrin IX (PPIX) or tetra(p-sulfophenyl)porphyrin (TPPS) and L = an aliphatic or aromatic nitrogenous base) are reported and discussed. The results are compared to those of previously reported [Fe(II)(Por)L(CO)] complexes (where Por = PPIX, TPPS, PMXPP, TPP, OMTBP and OEP; L = a nitrogenous aromatic ligand) and HbCO (where Hb = haemoglobin) and MyCO (where My = myoglobin). A new approach, to extracting information from the Mössbauer parameters has been developed by plotting those of the [Fe(II)(Por)L₂] complexes against those of [Fe(II)(Por)L(CO)] complexes for the same ligands, has yielded a series of trend lines that show a significant dependence on both the nature of the porphyrin and also of the nitrogenous ligand. Different trend lines were found for aromatic nitrogenous ligands to aliphatic nitrogenous ligands showing that the porphyrins could donate different amounts of charge to the Fe(II) cations as the L ligand changed, and hence, they display electron sink properties. From the plots, it was shown that haemoglobin and myoglobin both bind CO very strongly compared to the model complexes studied herein. Using the reported structural and Mössbauer data for the [Fe(II)(Por)L₂] and [Fe(II)(Por)L(CO)] complexes, it proved possible and instructive to plot the Mössbauer parameters against a number of the bond lengths around the Fe(II) cations. The interpretation of the resulting trend lines both supported and facilitated the extension of our findings enabling further understanding of the geometry of the bonding in CO haemoglobin and CO myoglobin.

57Fe Mössbauer isomer shift of pure iron and iron oxides at high pressure-An experimental and theoretical study

The ⁵⁷Fe isomer shift (IS) of pure iron has been measured up to 100 GPa using synchrotron Mössbauer spectroscopy in the time domain. Apart from the expected discontinuity due to the α → ε structural and spin transitions, the IS decreases monotonically with increasing pressure. The absolute shifts were reproduced without semi-empirical calibrations by periodic density functional calculations employing extensive localized basis sets with several common density functionals. However, the best numerical agreement is obtained with the B1WC hybrid functional. Extension of the calculations to 350 GPa, a pressure corresponding to the Earth's inner core, predicted the IS range of 0.00 to -0.85 mm/s, covering the span from Fe(0) to Fe(VI) compounds measured at ambient pressure. The calculations also reproduced the pressure trend from polymorphs of prototypical iron oxide minerals, FeO and Fe₂O₃. Analysis of the electronic structure shows a strong donation of electrons from oxygen to iron at high pressure. The assignment of formal oxidation to the Fe atom becomes ambiguous under this condition.

Effect of the Electron Density of the Heme Fe Atom on the Nature of Fe-O2 Bonding in Oxy Myoglobin

Mössbauer spectroscopy has been used to characterize oxygenated myoglobins (oxy Mbs) reconstituted with native and chemically modified ⁵⁷Fe-enriched heme cofactors with different electron densities of the heme Fe atom (ρ_Fe) and to elucidate the effect of a change in the ρ_Fe on the nature of the bond between heme Fe and oxygen (O₂), i.e., the Fe-O₂ bond, in the protein. Quadrupole splitting (ΔE_Q) was found to decrease with decreasing ρ_Fe, and the observed ρ_Fe-dependent ΔE_Q confirmed an increase in the contribution of the ferric-superoxide (Fe³⁺-O₂^-) form to the resonance hybrid of the Fe-O₂ fragment with decreasing ρ_Fe. These observations explicitly accounted for the lowering of O₂ affinity of the protein due to an increase in the O₂ dissociation rate and a decrease in the autoxidation reaction rate of oxy Mb through decreasing H⁺ affinity of the bound ligand with decreasing ρ_Fe. Therefore, the present study demonstrated the mechanism underlying the electronic control of O₂ affinity and the autoxidation of the protein through the heme electronic structure. Carbon monoxide (CO) adducts of reconstituted Mbs (CO-Mbs) were similarly characterized, and we found that the resonance between the two canonical forms of the Fe-CO fragment was also affected by a change in ρ_Fe. Thus, the nature of the Fe-ligand bond in the protein was found to be affected by the ρ_Fe.

Reliable Estimation of Prediction Uncertainty for Physicochemical Property Models

One of the major challenges in computational science is to determine the uncertainty of a virtual measurement, that is the prediction of an observable based on calculations. As highly accurate first-principles calculations are in general unfeasible for most physical systems, one usually resorts to parameteric property models of observables, which require calibration by incorporating reference data. The resulting predictions and their uncertainties are sensitive to systematic errors such as inconsistent reference data, parametric model assumptions, or inadequate computational methods. Here, we discuss the calibration of property models in the light of bootstrapping, a sampling method that can be employed for identifying systematic errors and for reliable estimation of the prediction uncertainty. We apply bootstrapping to assess a linear property model linking the ⁵⁷Fe Mössbauer isomer shift to the contact electron density at the iron nucleus for a diverse set of 44 molecular iron compounds. The contact electron density is calculated with 12 density functionals across Jacob's ladder (PWLDA, BP86, BLYP, PW91, PBE, M06-L, TPSS, B3LYP, B3PW91, PBE0, M06, TPSSh). We provide systematic-error diagnostics and reliable, locally resolved uncertainties for isomer-shift predictions. Pure and hybrid density functionals yield average prediction uncertainties of 0.06-0.08 mm s^-1 and 0.04-0.05 mm s^-1, respectively, the latter being close to the average experimental uncertainty of 0.02 mm s^-1. Furthermore, we show that both model parameters and prediction uncertainty depend significantly on the composition and number of reference data points. Accordingly, we suggest that rankings of density functionals based on performance measures (e.g., the squared coefficient of correlation, r², or the root-mean-square error, RMSE) should not be inferred from a single data set. This study presents the first statistically rigorous calibration analysis for theoretical Mössbauer spectroscopy, which is of general applicability for physicochemical property models and not restricted to isomer-shift predictions. We provide the statistically meaningful reference data set MIS39 and a new calibration of the isomer shift based on the PBE0 functional.

Comment on “Calibration of 57Fe Mössbauer constants by first principles” Phys. Chem. Chem. Phys., 2016, 18, 10201–10206

The proportionality constant between the iron-57 Mössbauer isomer shift and the electron probablility density at the nucleus is reevaluated from the correct experimental isomer shifts.

On Predicting Mössbauer Parameters of Iron-Containing Molecules with Density-Functional Theory

The performance of six frequently used density functional theory (DFT) methods (RPBE, OLYP, TPSS, B3LYP, B3LYP*, and TPSSh) in the prediction of Mössbauer isomer shifts(δ) and quadrupole splittings (ΔEQ) is studied for an extended and diverse set of Fe complexes. In addition to the influence of the applied density functional and the type of the basis set, the effect of the environment of the molecule, approximated with the conducting-like screening solvation model (COSMO) on the computed Mössbauer parameters, is also investigated. For the isomer shifts the COSMO-B3LYP method is found to provide accurate δ values for all 66 investigated complexes, with a mean absolute error (MAE) of 0.05 mm s-1 and a maximum deviation of 0.12 mm s-1. Obtaining accurate ΔEQ values presents a bigger challenge; however, with the selection of an appropriate DFT method, a reasonable agreement can be achieved between experiment and theory. Identifying the various chemical classes of compounds that need different treatment allowed us to construct a recipe for ΔEQ calculations; the application of this approach yields a MAE of 0.12 mm s-1 (7% error) and a maximum deviation of 0.55 mm s-1 (17% error). This accuracy should be sufficient for most chemical problems that concern Fe complexes. Furthermore, the reliability of the DFT approach is verified by extending the investigation to chemically relevant case studies which include geometric isomerism, phase transitions induced by variations of the electronic structure (e.g., spin crossover and inversion of the orbital ground state), and the description of electronically degenerate triplet and quintet states. Finally, the immense and often unexploited potential of utilizing the sign of the ΔEQ in characterizing distortions or in identifying the appropriate electronic state at the assignment of the spectral lines is also shown.

NMR, IR, Mössbauer and quantum chemical investigations of metalloporphyrins and metalloproteins

We review contributions made towards the elucidation of CO and O2binding geometries in respiratory proteins. Nuclear magnetic resonance, infrared spectroscopy, Mössbauer spectroscopy, X-ray crystallography and quantum chemistry have all been used to investigate the Fe –ligand interactions. Early experimental results showed linear correlations between17O chemical shifts and the infrared stretching frequency (νCO) of the CO ligand in carbonmonoxyheme proteins and between the17O chemical shift and the13CO shift. These correlations led to early theoretical investigations of the vibrational frequency of carbon monoxide and of the13C and17O NMR chemical shifts in the presence of uniform and non-uniform electric fields. Early success in modeling these spectroscopic observables then led to the use of computational methods, in conjunction with experiment, to evaluate ligand-binding geometries in heme proteins. Density functional theory results are described which predict57Fe chemical shifts and Mössbauer electric field gradient tensors,17O NMR isotropic chemical shifts, chemical shift tensors and nuclear quadrupole coupling constants (e2qQ/h) as well as13C isotropic chemical shifts and chemical shift tensors in organometallic clusters, heme model metalloporphyrins and in metalloproteins. A principal result is that CO in most heme proteins has an essentially linear and untilted geometry (τ = 4 °, β = 7 °) which is in extremely good agreement with a recently published X-ray synchrotron structure. CO / O2discrimination is thus attributable to polar interactions with the distal histidine residue, rather than major Fe–C–O geometric distortions.

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.

Experimental and theoretical Mössbauer study of an extended family of [Fe8(μ4-O)4(μ-4-R-px)12X4] clusters

Six [Fe(8)(μ(4)-O)(4)(μ-4-R-pyrazolato)(12)X(4)] complexes containing an identical Fe(8)(μ(4)-O)(4) core have been structurally characterized and studied by Mössbauer spectroscopy. In each case, an inner μ(4)-O bridged Fe(III) cubane core is surrounded by four trigonal bipyramidal iron centers, the two distinct sites occurring in a 1:1 ratio. The Mössbauer spectrum of each of the clusters consists of two quadrupole doublets, which, with one exception (X = NCS, R = H), overlap to give three absorption lines. The systematic variation of X and R causes significant changes in the Mössbauer spectra. A comparison with values for the same clusters computed using density functional theory allows us to establish an unequivocal assignment of these peaks in terms of a nested model for the overlapping doublets. The changes in Mössbauer parameters (both experimental and computed) for the 1-electron reduced species [Fe(8)(μ(4)-O)(4)(μ-4-Cl-pyrazolato)(12)Cl(4)](-) are consistent with a redox event that is localized within the cubane core.

The prediction of Fe Mössbauer parameters by the density functional theory: a benchmark study

We report the performance of eight density functionals (B3LYP, BPW91, OLYP, O3LYP, M06, M06-2X, PBE, and SVWN5) in two Gaussian basis sets (Wachters and Partridge-1 on iron atoms; cc-pVDZ on the rest of atoms) for the prediction of the isomer shift (IS) and the quadrupole splitting (QS) parameters of Mössbauer spectroscopy. Two sources of geometry (density functional theory-optimized and X-ray) are used. Our data set consists of 31 iron-containing compounds (35 signals), the Mössbauer spectra of which were determined at liquid helium temperature and where the X-ray geometries are known. Our results indicate that the larger and uncontracted Partridge-1 basis set produces slightly more accurate linear correlations of electronic density used for the prediction of IS and noticeably more accurate results for the QS parameter. We confirm and discuss the earlier observation of Noodleman and co-workers that different oxidation states of iron produce different IS calibration lines. The B3LYP and O3LYP functionals have the lowest errors for either IS or QS. BPW91, OLYP, PBE, and M06 have a mixed success whereas SVWN5 and M06-2X demonstrate the worst performance. Finally, our calibrations and conclusions regarding the best functional to compute the Mössbauer characteristics are applied to candidate structures for the peroxo and Q intermediates of the enzyme methane monooxygenase hydroxylase (MMOH), and compared to experimental data in the literature.

Density functional theory calculations on Mössbauer parameters of nonheme iron nitrosyls

Density Functional Theory (DFT) calculations on transition metal nitrosyls often reveal unusual spin density profiles, involving substantial spatial separation of majority and minority spin densities. Against this context, there is a significant lack of studies where DFT calculations have been quantitatively calibrated against experimental spectroscopic properties. Reported herein are DFT calculations of Mössbauer isomer shifts and quadrupole splittings for 21 nonheme iron complexes (26 distinct iron sites) including 9 iron nitrosyls. Low- (S = 1/2) and high-spin (S = 3/2) {FeNO}(7) complexes, S = 1/2 {Fe(NO)(2)}(9) species, and polynuclear iron nitrosyls are all represented within the set of compounds examined. The general conclusion with respect to isomer shifts is that DFT (OLYP/STO-TZP) performs comparably well for iron nitrosyls and for iron complexes in general. However, quadrupole splittings are less accurately reproduced for nitrosyl complexes.

Resilience of the iron environment in heme proteins

Conformational flexibility is essential to the functional behavior of proteins. We use an effective force constant introduced by Zaccai, the resilience, to quantify this flexibility. Site-selective experimental and computational methods allow us to determine the resilience of heme protein active sites. The vibrational density of states of the heme Fe determined using nuclear resonance vibrational spectroscopy provides a direct experimental measure of the resilience of the Fe environment, which we compare quantitatively with values derived from the temperature dependence of atomic mean-squared displacements in molecular dynamics simulations. Vibrational normal modes in the THz frequency range dominate the resilience. Both experimental and computational methods find a higher resilience for cytochrome c than for myoglobin, which we attribute to the increased number of covalent links to the peptide in the former protein. For myoglobin, the resilience of the iron environment is larger than the average resilience previously determined for hydrogen sites using neutron scattering. Experimental results suggest a slightly reduced resilience for cytochrome c upon oxidation, although the change is smaller than reported in previous Mössbauer investigations on a bacterial cytochrome c, and is not reproduced by the simulations. Oxidation state also has no significant influence on the compressibility calculated for cyt c, although a slightly larger compressibility is predicted for myoglobin.

Influence of Hartree−Fock Exchange on the Calculated Mössbauer Isomer Shifts and Quadrupole Splittings in Ferrocene Derivatives Using Density Functional Theory

Influence of molecular geometry, type of exchange-correlation functional, and contraction scheme of basis set applied at the iron nuclei have been tested in the calculation of 57Fe Mössbauer isomer shifts and quadrupole splittings for a wide range of ligand types, as well as oxidation and spin states, in inorganic and organometallic systems. It has been found that uncontraction of the s-part of Wachter's full-electron basis set at the iron nuclei does not appreciably improve the calculated isomer shifts. The observed correlations for all tested sets of geometries are close to each other and predominantly depend on the employed exchange-correlation functional with B3LYP functional being slightly better as compared to BPW91. Both hybrid (B3LYP) and pure (BPW91) exchange-correlation functionals are suitable for the calculation of isomer shifts in organometallic compounds. Surprisingly, it has been found that the hybrid B3LYP exchange-correlation functional completely fails in accurate prediction of quadrupole splittings in ferrocenes, while performance of the pure BPW91 functional for the same systems was excellent. This observation has been explained on the basis of relationship between the amount of Hartree-Fock exchange involved in the applied exchange-correlation functional and the calculated HOMO-LUMO energy gap in ferrocenes. On the basis of this explanation, use of only pure exchange-correlation functionals has been suggested for accurate prediction of Mössbauer spectra parameters in ferrocenes.

Synthesis, Structure, and Properties of an Fe(II) Carbonyl [(PaPy3)Fe(CO)](ClO4): Insight into the Reactivity of Fe(II)−CO and Fe(II)−NO Moieties in Non-Heme Iron Chelates of N-Donor Ligands

An Fe(II) carbonyl complex [(PaPy3)Fe(CO)](ClO4) (1) of the pentadentate ligand N,N-bis(2-pyridylmethyl)amine-N-ethyl-2-pyridine-2-carboxamide (PaPy3H, H is the dissociable amide proton) has been synthesized and structurally characterized. This Fe(II) carbonyl exhibits its nu(CO) at 1972 cm(-1), and its 1H NMR spectrum in degassed CD3CN confirms its S = 0 ground state. The bound CO in 1 is not photolabile. Reaction of 1 with an equimolar amount of NO results in the formation of the {Fe-NO}7 nitrosyl [(PaPy3)Fe(NO)](ClO4) (2), while excess NO affords the iron(III) nitro complex [(PaPy3)Fe(NO2)](ClO4) (5). In the presence of [Fe(Cp)2]+ and excess NO, 1 forms the {Fe-NO}6 nitrosyl [(PaPy3)Fe(NO)](ClO4)2 (3). Complex 1 also reacts with dioxygen to afford the iron(III) mu-oxo species [{(PaPy3)Fe}2O](ClO4)2 (4). Comparison of the metric and spectral parameters of 1 with those of the previously reported {Fe-NO}6,7 nitrosyls 3 and 2 provides insight into the electronic distributions in the Fe(II)-CO, Fe(II)-NO, and Fe(II)-NO+ bonds in the isostructural series of complexes 1-3 derived from a non-heme polypyridine ligand with one carboxamide group.

Active site structure of class I ribonucleotide reductase intermediate X: a density functional theory analysis of structure, energetics, and spectroscopy

Several models for the active site structure of class I ribonucleotide reductase (RNR) intermediate X have been studied in the work described in this paper, using broken-symmetry density functional theory (DFT) incorporated with the conductor-like screening (COSMO) solvation model. The calculated properties, including geometries, spin states, 57Fe, 1H, and 17O hyperfine tensors, Mössbauer isomer shifts, and quadrupole splittings, and the estimation of the Fe(IV) d-d transition energies have been compared with the available experimental values. On the basis of the detailed analysis and comparisons, we propose a definite form for the active site structure of class I RNR intermediate X, which contains an Fe1(III)Fe2(IV) center (where Fe1 is the iron site closer to Tyr122, and the two iron sites are high-spin antiferromagnetically coupled to give a total 1/2 net spin), two mu-oxo bridges, one terminal water which binds to Fe1(III) and also H-bonds to both side chains of Asp84 and Glu238, and one bidentate carboxylate group from the side chain of Glu115.

DFT calculations of57Fe Mössbauer isomer shifts and quadrupole splittings for iron complexes in polar dielectric media: Applications to methane monooxygenase and ribonucleotide reductase

To predict the isomer shifts of Fe complexes in different oxidation and spin states more accurately, we have performed linear regression between the measured isomer shifts (delta(exp)) and DFT (PW91 potential with all-electron triple-zeta plus polarization basis sets) calculated electron densities at Fe nuclei [rho(0)] for the Fe(2+,2.5+) and Fe(2.5+,3+,3.5+,4+) complexes separately. The geometries and electronic structures of all complexes in the training sets are optimized within the conductor like screening (COSMO) solvation model. Based on the linear correlation equation delta(exp) = alpha[rho(0) - 11884.0] + C, the best fitting for 17 Fe(2+,2.5+) complexes (totally 31 Fe sites) yields alpha = -0.405 +/- 0.042 and C = 0.735 +/- 0.047 mm s(-1). The correlation coefficient is r = -0.876 with a standard deviation of SD = 0.075 mm s(-1). In contrast, the linear fitting for 19 Fe(2.5+,3+,3.5+,4+) complexes (totally 30 Fe sites) yields alpha = -0.393 +/- 0.030 and C = 0.435 +/- 0.014 mm s(-1), with the correlation coefficient r = -0.929 and a standard deviation SD = 0.077 mm s(-1). We provide a physical rationale for separating the Fe(2+,2.5+) fit from the Fe(2.5+,3+,3.5+,4+) fit, which also is clearly justified on a statistical empirical basis. Quadrupole splittings have also been calculated for these systems. The correlation between the calculated (DeltaE(Q(cal))) and experimental (DeltaE(Q(exp))) quadrupole splittings based on |DeltaE(Q(exp))| = A |DeltaE(Q(cal))| + B yields slope A, which is almost the ideal value 1.0 (A = 1.002 +/- 0.030) and intercept B almost zero (B = 0.033 +/- 0.068 mm s(-1)). Further calculations on the reduced diferrous and oxidized diferric active sites of class-I ribonucleotide reductase (RNR) and the hydroxylase component of methane monooxygenase (MMOH), and on a mixed-valent [(tpb)Fe3+(mu-O)(mu-CH3CO2)Fe4+(Me3[9]aneN3)]2+ (S = 3/2) complex and its corresponding diferric state have been performed. Calculated results are in very good agreement with the experimental data.

Heme carbonyls: environmental effects on nu(C-O) and Fe-C/C-O bond length correlations

The synthesis and characterization of four low-spin (carbonyl)iron(II) tetraphenylporphyrinates, [Fe(TPP)(CO)(L)], where L = 1-methylimidazole, 2-methylimidazole, 1,2-dimethylimidazole (unsolvated), and 1,2-dimethylimidazole (toluene solvate) are reported. The complexes show nearly the same value of nu(C-O) in toluene solution (1969-72 cm(-1)) but a large range of CO stretching frequencies in the solid-state (1926-1968 cm(-1)). The large solid-state variation results from CO interactions in the solid state, as shown by an examination of the crystal structures of the four complexes. The high precision of the four structures obtained allows us to make a number of structural and spectroscopic correlations that describe the Fe-C-O and N(Im)-Fe-CO units. The values of nu(C-O) and the Fe-C and C-O bond distances are strongly correlated and provide a structural, as well as a spectroscopic, correlation of the pi back-bonding model. The interactions of CO described are closely related to the large range of CO stretching frequencies observed in heme proteins and specific interactions observed in carbonylmyoglobin (MbCO).

Variable pi-bonding in iron(II) porphyrinates with nitrite, CO, and tert-butyl isocyanide: characterization of [Fe(TpivPP)(NO2)(CO)]-

The addition of the strongly pi-bonding ligands CO or tert-butyl isocyanide to the low-spin five-coordinate iron(II) nitrite species [Fe(TpivPP)(NO2)]- (TpivPP = picket fence porphyrin) gives two new six-coordinate species [Fe(TpivPP)(NO2)(CO)]- and [Fe(TpivPP)(NO2)(t-BuNC)]-. These species have been characterized by single-crystal structure determinations and by UV-vis, IR, and Mössbauer spectroscopies. All evidence shows that in the mixed-ligand iron(II) porphyrin species, [Fe(TpivPP)(NO2)(CO)]-, the two trans, pi-accepting ligands CO and nitrite compete for pi density. The CO ligand however dominates the bonding. The Fe-N(NO2) bond lengths for the two independent anions in the unit cell at 2.006(4) and 2.009(4) A are lengthened compared to other nitrite species with either no trans ligands or non-pi-accepting trans ligands to nitrite. The Fe-C(CO) bond lengths are 1.782(4) A and 1.789(5) A for the two anions. The two Fe-C-O angles at 175.5(4) and 177.5(4) degrees are essentially linear in both anions. The quadrupole splitting for [Fe(TpivPP)(NO2)(CO)]- was determined to be 0.32 mm/s, and the isomer shift was 0.18 mm/s at room temperature in zero applied field. Both of the Mössbauer parameters are much smaller than those found for six-coordinate low-spin iron(II) porphyrinates with neutral nitrogen-donating ligands as well as iron(II) nitro complexes. However, the Mössbauer parameters are typical of other six-coordinate CO porphyrinates signifying that CO is the more dominant ligand. The CO stretching frequency of 1974 cm(-1) is shifted only slightly to higher energy compared to six-coordinate CO complexes with neutral nitrogen-donor ligands trans to CO. Crystal data for [K(222)][Fe(TpivPP)(NO2)(CO)].1/2C6H5Cl: monoclinic, space group P2(1)/c, Z = 8, a = 33.548(6) A, b = 18.8172(15) A, c = 27.187(2) A, beta = 95.240(7) degrees, V = 17091(4) A3.

Density functional study of the electric hyperfine interactions and the redox-structural correlations in the cofactor of nitrogenase. Analysis of general trends in (57)Fe isomer shifts

The influence of the interstitial atom, X, discovered in a recent crystallographic study of the MoFe protein of nitrogenase, on the electric hyperfine interactions of (57)Fe has been investigated with density functional theory. A semiempirical theory for the isomer shift, delta, is formulated and applied to the cofactor. The values of delta for the relevant redox states of the cofactor are predicted to be higher in the presence of X than in its absence. The analysis strongly suggests a [Mo(4+)4Fe(2+)3Fe(3+)] oxidation state for the S = 3/2 state M(N). Among C(4-), N(3-), and O(2-), oxide is found to be the least likely candidate for X. The analysis suggests that X should be present in the cofactor states M(OX) and M(R) as well as in the alternative nitrogenases. The calculations of the electric field gradients (EFGs) indicate that the small values for DeltaE(Q) in M(N) result from an extensive cancellation between valence and ligand contributions. X emerges from the analysis of the hyperfine interactions as an ionically bonded species. Its major effect is on the asymmetry parameters for the EFGs at the six equatorial sites, Fe(Eq). A spin-coupling scheme is proposed for the state [Mo(4+)4Fe(2+)3Fe(3+)] that is consistent with the measured (57)Fe A-tensors and DeltaE(Q) values for M(N) and identifies the unique site exhibiting the small A value with the terminal Fe site, Fe(T). The optimized structure of a cofactor model has been calculated for several oxidation states. The study reveals a contraction in the average Fe-Fe distance upon increasing the number of electrons stored in the cluster, in accord with extended X-ray absorption fine structure studies. The reliability of the adopted methodology for predicting redox-structural correlations is tested for cuboidal [4Fe-4S] clusters. The calculations reveal a systematic increase in the S...S sulfide distances, in quantitative agreement with the available data. These trends are rationalized by a simple electrostatic model.

Five- to six-coordination in (nitrosyl)iron(II) porphyrinates: effects of binding the sixth ligand

We report structural and spectroscopic data for a series of six-coordinate (nitrosyl)iron(II) porphyrinates. The structures of three tetraphenylporphyrin complexes [Fe(TPP)(NO)(L)], where L = 4-(dimethylamino)pyridine, 1-methylimidazole, 4-methylpiperidine, are reported here to a high degree of precision and allow observation of several previously unobserved structural features. The tight range of bonding parameters for the [FeNO] moiety for these three complexes suggests a canonical representation for six-coordinate systems (Fe-N(p) = 2.007 A, Fe-N(NO) = 1.753 A, angle FeNO = 138.5 degrees ). Comparison of these data with those obtained previously for five-coordinate systems allows the precise determination of the structural effects of binding a sixth ligand. These include lengthening of the Fe-N(NO) bond and a decrease in the Fe-N-O angle. Several other aspects of the geometry of these systems are also discussed, including the first examples of off-axis tilting of a nitrosyl ligand in a six-coordinate [FeNO](7) heme system. We also report the first examples of Mössbauer studies for these complexes. Measurements have been made in several applied magnetic fields as well as in zero field. The spectra differ from those of their five-coordinate analogues. To obtain reasonable fits to applied magnetic field data, rotation of the electrical field gradient is required, consistent with differing g-tensor orientations in the five- vs six-coordinate species.

Effect of the sixth axial ligand in CS-ligated iron(II)octaethylporphyrinates: structural and Mössbauer studies

The effect of a sixth ligand in a series of low-spin thiocarbonyl-ligated iron(II)octaethylporphyrinates has been investigated. Six-coordinate complexes have been synthesized and characterized by Mössbauer and infrared spectroscopy and single-crystal X-ray structure determinations. The results are compared with the five-coordinate parent complex. The crystal structures of [Fe(OEP)(CS)(1-MeIm)] and [Fe(OEP)(CS)(Py)] are reported and discussed. The 1-methylimidazole and pyridine derivatives exhibit Fe-C(CS) bond distances of 1.703(4) and 1.706(2) A that are significantly longer than the 1.662(3) A reported for five-coordinate [Fe(OEP)(CS)] (Scheidt, W. R.; Geiger, D. K. Inorg. Chem. 1982, 21, 1208). The trans Fe-N(ligand) distances of 2.112(3) and 2.1550(15) A observed for the 1-methylimidazole and pyridine complex are approximately 0.13 A longer than those observed for analogous bis-ligated complexes and are consistent with a significant structural trans effect for the CS ligand. Mössbauer investigations carried out for five- and six-coordinate thiocarbonyl derivatives with several different sixth axial ligands reveal interesting features. All derivatives exhibit very small isomer shift values, consistent with a very strong interaction between iron and CS. The five-coordinate derivative has delta(Fe) = 0.08 mm/s, and the six-coordinate complexes exhibit delta(Fe) = 0.14 to 0.19 mm/s at 4.2 K. The five-coordinate complex shows a large quadrupole splitting (DeltaE(q) = 1.93 mm/s at 4.2 K) which is reduced on coordination of the sixth ligand (DeltaE(q) = 0.42-0.80 mm/s at 4.2 K). Addition of a sixth ligand also leads to a small decrease in the value of nu(CS). Correlations in structural, IR, and Mössbauer results suggest that the sixth ligand effect is primarily induced by changes in sigma-bonding. The structure of [Fe(OEP)(CS)(CH(3)OH)] is briefly reported. Crystal data: [Fe(OEP)(CS)(1-MeIm)] crystallizes in the monoclinic system, space group P2(1)/n, Z = 4, a = 9.5906(5) A, b = 16.704(4) A, c = 23.1417(6) A, beta = 100.453(7) degrees. [Fe(OEP)(CS)(Py)] crystallizes in the triclinic system, space group P1, Z = 5, a = 13.9073(6) A, b = 16.2624(7) A, c = 22.0709(9) A, alpha = 70.586(1) degrees, beta = 77.242(1) degrees, gamma = 77.959(1) degrees. [Fe(OEP)(CS)(CH(3)OH)] crystallizes in the triclinic system, space group P1, Z = 1, a = 9.0599(5) A, b = 9.4389(5) A, c = 11.0676(6) A, alpha = 90.261(1) degrees, beta = 100.362(1) degrees, gamma = 114.664(1) degrees.

Characterization of electronic structure and properties of a Bis(histidine) heme model complex

Ferric and ferrous hemes, such as those present in electron transfer proteins, often have low-lying spin states that are very close in energy. To explore the relationship between spin state, geometry, and cytochrome electron transfer, we investigate, using density functional theory, the relative energies, electronic structure, and optimized geometries for a high- and low-spin ferric and ferrous heme model complex. Our model consists of an iron-porphyrin axially ligated by two imidazoles, which model the interaction of a heme with histidine residues. Using the B3LYP hybrid functional, we found that, in the ferric model heme complex, the doublet is lower in energy than the sextet by 8.4 kcal/mol and the singlet ferrous heme is 6.7 kcal/mol more stable than the quintet. The difference between the high-spin ferric and ferrous model heme energies yields an adiabatic electron affinity (AEA) of 5.24 eV, and the low-spin AEA is 5.17 eV. Both values are large enough to ensure electron trapping, and electronic structure analysis indicates that the iron d(pi) orbital is involved in the electron transfer between hemes. Mössbauer parameters calculated to verify the B3LYP electronic structure correlate very well with experimental values. Isotropic hyperfine coupling constants for the ligand nitrogen atoms were also evaluated. The optimized geometries of the ferric and ferrous hemes are consistent with structures from X-ray crystallography and reveal that the iron-imidazole distances are significantly longer in the high-spin hemes, which suggests that the protein environment, modeled here by the imidazoles, plays an important role in regulating the spin state. Iron-imidazole dissociation energies, force constants, and harmonic frequencies were calculated for the ferric and ferrous low-spin and high-spin hemes. In both the ferric and the ferrous cases, a single imidazole ligand is more easily dissociated from the high-spin hemes.

Quantum chemical calculations of spectroscopic properties of metalloproteins and model compounds: EPR and Mössbauer properties

Recently developed theoretical methods to predict EPR and Mössbauer parameters open the way for close interactions between theorists and experimentalists to elucidate the geometric and electronic structures of metalloenzymes and model complexes and to obtain insight into their reactive properties. Spectral calculations (g-values, hyperfine couplings, zero-field splittings, isomer shifts and quadrupole splittings) are also a means to validate theoretical approaches and therefore complement the prediction of geometries, reaction energies and transition states.

Nuclear magnetic resonance shifts in paramagnetic metalloporphyrins and metalloproteins

We report the first detailed investigation of the (1)H, (13)C, (15)N, and (19)F nuclear magnetic resonance (NMR) spectroscopic shifts in paramagnetic metalloprotein and metalloporphyrin systems. The >3500 ppm range in experimentally observed hyperfine shifts can be well predicted by using density functional theory (DFT) methods. Using spin-unrestricted methods together with large, locally dense basis sets, we obtain very good correlations between experimental and theoretical results: R(2) = 0.941 (N = 37, p < 0.0001) when using the pure BPW91 functional and R(2) = 0.981 (N = 37, p < 0.0001) when using the hybrid functional, B3LYP. The correlations are even better for C(alpha) and C(beta) shifts alone: C(alpha), R(2) = 0.996 (N = 8, p < 0.0001, B3LYP); C(beta), R(2) = 0.995 (N = 8, p < 0.0001, B3LYP), but are worse for C(meso), in part because of the small range in C(meso) shifts. The results of these theoretical calculations also lead to a revision of previous heme and proximal histidine residue (13)C NMR assignments in deoxymyoglobin which are confirmed by new quantitative NMR measurements. Molecular orbital (MO) analyses of the resulting wave functions provide a graphical representation of the spin density distribution in the [Fe(TPP)(CN)(2)](-) (TPP = 5,10,15,20-tetraphenylporphyrinato) system (S = (1)/(2)), where the spin density is shown to be localized primarily in the d(xz) (or d(yz)) orbital, together with an analysis of the frontier MOs in Fe(TPP)Cl (S = (5)/(2)), Mn(TPP)Cl (S = 2), and a deoxymyoglobin model (S = 2). The ability to now begin to predict essentially all heavy atom NMR hyperfine shifts in paramagnetic metalloporphyrins and metalloproteins using quantum chemical methods should open up new areas of research aimed at structure prediction and refinement in paramagnetic systems in much the same way that DFT methods have been used successfully in the past to predict/refine elements of diamagnetic heme protein structures.

Mössbauer quadrupole splittings and electronic structure in heme proteins and model systems: a density functional theory investigation

We report the results of a series of density functional theory (DFT) calculations aimed at predicting the (57)Fe Mössbauer electric field gradient (EFG) tensors (quadrupole splittings and asymmetry parameters) and their orientations in S = 0, (1)/(2), 1, (3)/(2), 2, and (5)/(2) metalloproteins and/or model systems. Excellent results were found by using a Wachter's all electron basis set for iron, 6-311G for other heavy atoms, and 6-31G for hydrogen atoms, BPW91 and B3LYP exchange-correlation functionals, and spin-unrestricted methods for the paramagnetic systems. For the theory versus experiment correlation, we found R(2) = 0.975, slope = 0.99, intercept = -0.08 mm sec(-)(1), rmsd = 0.30 mm sec(-)(1) (N = 23 points) covering a DeltaE(Q) range of 5.63 mm s(-)(1) when using the BPW91 functional and R(2) = 0.978, slope = 1.12, intercept = -0.26 mm sec(-)(1), rmsd = 0.31 mm sec(-)(1) when using the B3LYP functional. DeltaE(Q) values in the following systems were successfully predicted: (1) ferric low-spin (S = (1)/(2)) systems, including one iron porphyrin with the usual (d(xy))(2)(d(xz)d(yz))(3) electronic configuration and two iron porphyrins with the more unusual (d(xz)d(yz))(4)(d(xy))(1) electronic configuration; (2) ferrous NO-heme model compounds (S = (1)/(2)); (3) ferrous intermediate spin (S = 1) tetraphenylporphinato iron(II); (4) a ferric intermediate spin (S = (3)/(2)) iron porphyrin; (5) ferrous high-spin (S = 2) deoxymyoglobin and deoxyhemoglobin; and (6) ferric high spin (S = (5)/(2)) metmyoglobin plus two five-coordinate and one six-coordinate iron porphyrins. In addition, seven diamagnetic (S = 0, d(6) and d(8)) systems studied previously were reinvestigated using the same functionals and basis set scheme as used for the paramagnetic systems. All computed asymmetry parameters were found to be in good agreement with the available experimental data as were the electric field gradient tensor orientations. In addition, we investigated the electronic structures of several systems, including the (d(xy))(2)(d(xz),d(yz))(3) and (d(xz),d(yz))(4)(d(xy))(1) [Fe(III)/porphyrinate](+) cations as well as the NO adduct of Fe(II)(octaethylporphinate), where interesting information on the spin density distributions can be readily obtained from the computed wave functions.

Characterization of the relevant excited states in the photodissociation of CO-ligated hemoglobin and myoglobin

The relevant excited states in the rapid photodissociation process of hemoglobin and myoglobin are examined by means of time-dependent density functional theory. Our calculations clearly show that the photodissociation is mediated by two repulsive states (5 A' ' and 3 A') which cross the lowest excited states (1 A' and 1 A' ') at an internuclear Fe-C distance of about 2 A. Electron detachment/attachment density plots nicely explain the repulsive nature of the 5 A' ' and 3 A' states.

Chemical shifts in amino acids, peptides, and proteins: from quantum chemistry to drug design

This chapter discusses recent progress in the investigation and use of (13)C, (15)N, and (19)F nuclear magnetic resonance (NMR) chemical shifts and chemical shift tensors in proteins and model systems primarily using quantum chemical (ab initio Hartree-Fock and density functional theory) techniques. Correlations between spectra and structure are made and the techniques applied to other spectroscopic and electrostatic properties as well, including hydrogen bonding, ligand binding to heme proteins, J-couplings, electric field gradients, and atoms-in-molecules theory, together with a brief review of the use of NMR chemical shifts in drug design.

(57)Fe Mössbauer isomer shifts of heme protein model systems: electronic structure calculations

We report the results of density functional theory (DFT) calculations of the (57)Fe Mössbauer isomer shifts (delta(Fe)) for a series of 24 inorganic, organometallic, and metalloprotein/metalloporphyrin model systems in S = 0, (1)/(2), 1, (3)/(2), 2, and (5)/(2) spin states. We find an excellent correlation between calculation and experiment over the entire 2.34 mm s(-1) range of isomer shifts: a 0.07-0.08 mm s(-1) rms deviation between calculation and experiment (corresponding to 3-4% of the total delta(Fe) range, depending on the functionals used) with R(2) values of 0.973 and 0.981 (p < 0.0001). The best results are obtained by using the hybrid exchange-correlation functional B3LYP, used previously for (57)Fe Mössbauer quadrupole splittings and (57)Fe NMR chemical shifts and chemical shielding anisotropies. The relativistically corrected value of alpha, alpha(rel), converges with the large basis set used in this work, but the exact values vary somewhat with the methods used: -0.253 a(0)(3) mm s(-1) (Hartree-Fock; HF); -0.316 a(0)(3) mm s(-1) (hybrid HF-DFT; B3LYP), or -0.367 a(0)(3) mm s(-1) (pure DFT; BPW91). Both normal and intermediate spin state isomer shifts are well reproduced by the calculations, as is the broad range of delta(Fe) values: from [Fe(VI)O(4)](2-) (-0.90 mm s(-1) expt; -1.01 mm s(-1) calc) to KFe(II)F(3) (1.44 mm s(-1) expt; 1.46 mm s(-1) calc). Molecular orbital analyses of all inorganic solids as well as all organometallic and metalloporphyrin systems studied reveal that there are three major core MO contributions to rho(tot)(0), the total charge density at the iron nucleus (and hence delta(Fe)), that do not vary with changes in chemistry, while the valence MO contributions are highly correlated with delta(Fe) (R(2) = 0.915-0.938, depending on the functionals used), and the correlation between the valence MO contributions and the total MO contribution is even better (R(2) = 0.965-0.976, depending on the functionals used). These results are of general interest since they demonstrate that DFT methods now enable the accurate prediction of delta(Fe) values in inorganic, organometallic, and metalloporphyrin systems in all spin states and over a very wide range of delta(Fe) values with a very small rms error.