A QM/MM investigation of the activation and catalytic mechanism of Fe-only hydrogenases

Fe-only hydrogenases are enzymes that catalyze dihydrogen production or oxidation, due to the presence of an unusual Fe(6)S(6) cluster (the so-called H-cluster) in their active site, which is composed of a Fe(2)S(2) subsite, directly involved in catalysis, and a classical Fe(4)S(4) cubane cluster. Here, we present a hybrid quantum mechanical and molecular mechanical (QM/MM) investigation of the Fe-only hydrogenase from Desulfovibrio desulfuricans, in order to unravel key issues regarding the activation of the enzyme from its completely oxidized inactive state (Hoxinact) and the influence of the protein environment on the structural and catalytic properties of the H-cluster. Our results show that the Fe(2)S(2) subcluster in the Fe(II)Fe(II) redox state - which is experimentally observed for the completely oxidized form of the enzyme - binds a water molecule to one of its metal centers. The computed QM/MM energy values for water binding to the diferrous subsite are in fact over 70 kJ mol(-1); however, the affinity toward water decreases by 1 order of magnitude after a one-electron reduction of H(ox)(inact), thus leading to the release of coordinated water from the H-cluster. The investigation of a catalytic cycle of the Fe-only hydrogenase that implies formation of a terminal hydride ion and a di(thiomethyl)amine (DTMA) molecule acting as an acid/base catalyst indicates that all steps have reasonable reaction energies and that the influence of the protein on the thermodynamic profile of H(2) production catalysis is not negligible. QM/MM results show that the interactions between the Fe(2)S(2) subsite and the protein environment could give place to structural rearrangements of the H-cluster functional for catalysis, provided that the bidentate ligand that bridges the iron atoms in the binuclear subsite is actually a DTMA residue.     

Probing the effects of one-electron reduction and protonation on the electronic properties of the Fe-S clusters in the active-ready form of [FeFe]-hydrogenases. A QM/MM investigation

A QM/MM investigation of the active-ready (H(ox)) form of [FeFe]-hydrogenase from D. desulfuricans, in which the electronic properties of all Fe-S clusters (H, F and F') have been simultaneously described using DFT, was carried out with the aim of disclosing a possible interplay between the H-cluster and the accessory iron-sulfur clusters in the initial steps of the catalytic process leading to H(2) formation. It turned out that one-electron addition to the active-ready form leads to reduction of the F'-cluster and not of the H-cluster. Protonation of the H-cluster in H(ox) is unlikely, and in any case it would not trigger electron transfer from the accessory Fe(4)S(4) clusters to the active site. Instead, one-electron reduction and protonation of the active-ready form trigger electron transfer within the protein, a key event in the catalytic cycle. In particular, protonation of the H-cluster after one-electron reduction of the enzyme lowers the energy of the lowest unoccupied molecular orbitals localized on the H-cluster to such an extent that a long-range electron transfer from the F'-cluster towards the H-cluster itself is allowed.     

A theoretical study on the enhancement of functionally relevant electron transfers in biomimetic models of [FeFe]-hydrogenases

Recent advances aimed at modeling the chemistry of the active site of [FeFe]-hydrogenases (the H-cluster, composed by a catalytic Fe(2)S(2) subcluster and an Fe(4)S(4) portion) have led to the synthesis of binuclear coordination compounds containing a noninnocent organophosphine ligand [2,3-bis(diphenylphosphino)maleic anhydride, bma] that is able to undergo monoelectron reduction, analogously to the tetranuclear Fe(4)S(4) subcluster portion of the H-cluster. However, such a synthetic model was shown to feature negligible electronic communication between the noninnocent ligand and the remaining portion of the cluster, at variance with the enzyme active site. Here, we report a theoretical investigation that shows why the electron transfer observed in the enzyme upon protonation of the catalytic Fe(2)S(2) subsite cannot take place in the bma-containing cluster. In addition, we show that targeted modifications of the bma ligand are sufficient to restore the electronic communication within the model, such that electron density can be more easily withdrawn from the noninnocent ligand, as a result of protonation of the iron centers. Similar results were also obtained with a ligand derived from cobaltocene. The relevance of our findings is discussed from the perspective of biomimetic reproduction of proton reduction to yield molecular hydrogen.     

Evidence for the formation of terminal hydrides by protonation of an asymmetric iron hydrogenase active site mimic

Treatment of [Fe2(mu-pdt)(CO)6] [pdt=S(CH2)3S] with dppe (Ph2PCH2CH2PPh2) in refluxing toluene affords the asymmetric complex [Fe2(mu-pdt)(CO)4(dppe)] (1). Protonation of 1 with HBF4-Et2O in CH2Cl2 gives at room temperature the mu-hydrido derivative [Fe2(mu-pdt)(CO)4(dppe)(mu-H)](BF4) (2). Monitoring the reaction by 1H, 31P, and 13C NMR at low temperature reveals unambiguously that the process of the protonation of 1 implies terminal hydride intermediates.     

Density functional theory investigation of the active site of [Fe]-hydrogenases: effects of redox state and ligand characteristics on structural,electronic, and reactivity properties of complexes related to the [2Fe]H subcluster

The effects of redox state and ligand characteristics on structural, electronic, and reactivity properties of complexes related to the [2Fe](H) subcluster of [Fe]-hydrogenases have been investigated by DFT calculations and compared with experimental and theoretical data obtained investigating both the enzyme and synthetic model complexes. Our results show that Fe(II)Fe(II) species characterized by OH or H(2)O groups terminally coordinated to the iron atom distal to the terminal sulfur ligand (Fe(d)) are less stable than corresponding mu-OH or mu-H(2)O species, suggesting that the latter are destabilized or kinetically inaccessible in the enzyme. In addition, results obtained investigating Fe(I)Fe(I) and Fe(II)Fe(I) complexes show that structure and relative stability of species characterized by a mu-CO group are significantly affected by the electronic properties of the ligands coordinated to the iron atoms. The investigation of reaction pathways for H(2) activation confirms and extends a previous hypothesis indicating that H(2) can be cleaved on Fe(II)Fe(II) species. In particular, even though [Fe]-hydrogenases are proposed to bind and activate H(2) at a single iron center, the comparison of our data with experimental results obtained studying synthetic complexes (Zhao, X.; Georgakaki, I. P.; Miller, M. L.; Mejia-Rodriguez, R.; Chiang, C.-Y.; Darensbourg, M. Y. Inorg. Chem. 2002, 41, 3917) suggests that activation paths involving both metal ions are also possible. Moreover, mu-H Fe(II)Fe(I) complexes are predicted to correspond to stable species and might be formed in the enzyme catalytic cycle.     

Modeling [Fe-Fe] hydrogenase: evidence for bridging carbonyl and distal iron coordination vacancy in an electrocatalytically competent proton reduction by an iron thiolate assembly that operates through Fe(0)-Fe(II) levels

IR spectroelectrochemistry of Fe4{Me(CH2S)3}2(CO)8 (4Fe6S) in the nu(CO) region shows that the neutral and anion forms have all their CO groups terminally bound to the Fe atoms; however, for the dianion there is a switch of the coordination mode of at least one of the CO groups. The available structural and nu(CO) spectra are closely reproduced by density-functional theory calculations. The calculated structure of 4Fe6S2- closely mirrors that of the diiron subsite of the [Fe-Fe] hydrogenase H cluster with a bridging CO group and an open coordination site on the outer Fe atom of pairs of dithiolate-bridged Fe0FeII subunits connected by two bridging thiolates. Geometry optimization based on the all-terminal CO isomer of 4Fe6S2- does not give a stable structure but reveals a second-order saddle point ca. 11.53 kcal mol(-1) higher in energy than the CO-bridged form. Spectroelectrochemical studies of electrocatalytic proton reduction by 4Fe6S show that slow turnover from the primary reduction process (E1/2'=-0.71 V vs Ag/AgCl) involves rate-limiting protonation of 4Fe6S- followed by reduction to H:4Fe6S-. Rapid electrocatalytic proton reduction is obtained at potentials sufficient to access 4Fe6S2-, where the rate of dihydrogen elimination from the FeIIFeII core of 4Fe6S is ca. 500 times faster than that from the FeIFeI core of Fe2(mu-S(CH2)3S)(CO)6. The dramatically increased rate of electrocatalysis obtained from 4Fe6S over all previously identified model compounds appears to be related to the features uniquely common between it and the H-cluster, namely, that turnover involves the same formal redox states of the diiron unit (FeIFeII and Fe0FeII), the presence of an open site on the outer Fe atom of the Fe0FeII unit, and the thiolate-bridge to a second one-electron redox unit.     

Sulfur K-edge XAS and DFT calculations on nitrile hydratase: geometric and electronic structure of the non-heme iron active site

The geometric and electronic structure of the active site of the non-heme iron enzyme nitrile hydratase (NHase) is studied using sulfur K-edge XAS and DFT calculations. Using thiolate (RS(-))-, sulfenate (RSO(-))-, and sulfinate (RSO(2)(-))-ligated model complexes to provide benchmark spectral parameters, the results show that the S K-edge XAS is sensitive to the oxidation state of S-containing ligands and that the spectrum of the RSO(-) species changes upon protonation as the S-O bond is elongated (by approximately 0.1 A). These signature features are used to identify the three cysteine residues coordinated to the low-spin Fe(III) in the active site of NHase as CysS(-), CysSOH, and CysSO(2)(-) both in the NO-bound inactive form and in the photolyzed active form. These results are correlated to geometry-optimized DFT calculations. The pre-edge region of the X-ray absorption spectrum is sensitive to the Z(eff) of the Fe and reveals that the Fe in [FeNO](6) NHase species has a Z(eff) very similar to that of its photolyzed Fe(III) counterpart. DFT calculations reveal that this results from the strong pi back-bonding into the pi antibonding orbital of NO, which shifts significant charge from the formally t(2)(6) low-spin metal to the coordinated NO.     

Dihydrogen activation by a diruthenium analogue of the Fe-only hydrogenase active site

The photochemical reaction of Ru2(S2C3H6)(CO)4(PCy3)2 (1) and H2 gives the dihydride Ru2(S2C3H6)(mu-H)(H)(CO)3(PCy3)2 (2). NMR and crystallographic studies reveal mutually trans basal phosphine ligands and both bridging and terminal hydrides. Ru2(S2C2H4)(CO)4(PCy3)2 behaves similarly. Other HX substrates undergo photoaddition to 1, affording Ru2(S2C3H6)(mu-H)(X)(CO)3(PCy3)2 for X = OTs (3a), Cl (3b), and SPh (3c). Treatment of Ru2(S2C3H6)(mu-H)(H)(CO)3(PCy3)2 with [H(OEt2)]BArF4 (ArF = B(C6H3-3,5-(CF3)2) in CD2Cl2 gives [Ru2(S2C3H6)(mu-H)(CO)3(PCy3)2(H2)]+ (4), which catalyzes H2-D2 exchange. The reaction of 2 with [D(OEt2)]BArF4 gave [Ru2(S2C3H6)(mu-H)(CO)3(PCy3)2(HD)]+ (JH-D = 31 Hz). These studies provide the first models for the Fe-only hydrogenases that bear dihydrogen and terminal hydrido ligands.     

Photoreversible multiple additions of hydrogen to a highly unsaturated platinum-rhenium cluster complex

The compound Pt3Re2(CO)6(PBut3)3, 1, was obtained from the reaction of Re2(CO)10 with Pt(PBut3)2 in octane solvent at reflux. Compound 1 consists of a trigonal bipyramidal cluster of five metal atoms with three platinum atoms in the trigonal plane and the two rhenium atoms in the apical positions. The metal cluster is formally unsaturated by 10 electrons. Compound 1 sequentially adds 3 equiv of hydrogen at room temperature/1 atm to form the series of compounds Pt3Re2(CO)6(PBut3)3(mu-H)2, 2, Pt3Re2(CO)6(PBut3)3(mu-H)4, 3, and Pt3Re2(CO)6(PBut3)3(mu-H)6, 4. A small but significant kinetic isotope effect was observed, kH/kD = 1.3. The rate of addition of hydrogen is unaffected by the presence of a 20-fold excess of free PBut3 in solutions of 1. Compounds 2-4 each consist of a trigonal bipyramidal cluster of three platinum and two rhenium atoms similar to that of 1. The hydrido ligands in 2-4 bridge the platinum-rhenium bonds and are arranged to give structures having overall C2v symmetry for 2 and 3 and approximate D3h symmetry for 4. Some of the hydrido ligands were expelled from 4 in the form of hydrogen upon exposure of solutions to UV-vis irradiation to yield compound 3 and then 2 in reasonable yields, but the elimination of all hydrido ligands to yield 1 was achieved only under the most forcing UV irradiation and then only with a major loss of the complex due to decomposition. The electronic structures of 1-4 were investigated by DFT calculations. Additional DFT calculations have suggested some mechanisms for the activation of hydrogen at multicenter metal sites without ligand eliminations prior to the hydrogen additions.     

Proton reduction and dihydrogen oxidation on models of the [2Fe]H cluster of [Fe] hydrogenases. A density functional theory investigation

Density functional theory was used to compare reaction pathways for H2 formation and H+ reduction catalyzed by models of the binuclear cluster found in the active site of [Fe] hydrogenases. Terminal H+ binding to an Fe(I)-Fe(I) form, followed by monoelectron reduction and protonation of the di(thiomethyl)amine ligand, can conveniently lead to H2 formation and release, suggesting that this mechanism could be operative within the enzyme active site. However, a pathway that implies the initial formation of Fe(II)-Fe(II) mu-H species and release of H2 from an Fe(II)-Fe(I) form is characterized by only slightly less favored energy profiles. In both cases, H2 formation becomes less favored when taking into account the competition between CN and amine groups for H+ binding, an observation that can be relevant for the design of novel synthetic catalysts. H2 cleavage can take place on Fe(II)-Fe(II) redox species, in agreement with previous proposals [Fan, H.-J.; Hall, M. B. J. Am. Chem. Soc. 2001, 123, 3828] and, in complexes characterized by terminal CO groups, does not need the involvement of an external base. The step in H2 oxidation characterized by larger energy barriers corresponds to the second H+ extraction from the cluster, both considering Fe(II)-Fe(II) and Fe(II)-Fe(III) species. A comparison of the different reaction pathways reveals that H2 formation could involve only Fe(I)-Fe(I), Fe(II)-Fe(I), and Fe(II)-Fe(II) species, whereas Fe(III)-Fe(II) species might be relevant in H2 cleavage.     

Time resolved infrared spectroscopy: kinetic studies of weakly binding ligands in an iron-iron hydrogenase model compound

Solution photochemistry of (μ-pdt)[Fe(CO)(3)](2) (pdt = μ(2)-S(CH(2))(3)S), a precursor model of the 2-Fe subsite of the H-cluster of the hydrogenase enzyme, has been studied using time-resolved infrared spectroscopy. Following the loss of CO, solvation of the Fe center by the weakly binding ligands cyclohexene, 3-hexyne, THF, and 2,3-dihydrofuran (DHF) occurred. Subsequent ligand substitution of these weakly bound ligands by pyridine or cyclooctene to afford a more stable complex was found to take place via a dissociative mechanism on a seconds time scale with activation parameters consistent with such a pathway. That is, the ΔS(?) values were positive and the ΔH(?) parameters closely agreed with bond dissociation enthalpies (BDEs) obtained from DFT calculations. For example, for cyclohexene replacement by pyridine, experimental ΔH(?) and ΔS(?) values were determined to be 19.7 ± 0.6 kcal/mol (versus a theoretical prediction of 19.8 kcal/mol) and 15 ± 2 eu, respectively. The ambidentate ligand 2,3-DHF was shown to initially bind to the iron center via its oxygen atom followed by an intramolecular rearrangement to the more stable η(2)-olefin bound species. DFT calculations revealed a transition state structure with the iron atom almost equidistant from the oxygen and one edge of the olefinic bond. The computed ΔH(?) of 10.7 kcal/mol for this isomerization process was found to be in excellent agreement with the experimental value of 11.2 ± 0.3 kcal/mol.     

Reactivity of triruthenium thiophyne and furyne clusters: competitive S-C and P-C bond cleavage reactions and the generation of highly unsymmetrical alkyne ligands

The synthesis and reactivity of the thiophyne and furyne clusters [Ru3(CO)7(mu-dppm)(mu3-eta2-C4H2E)(mu-P(C4H3E)2)(mu-H)] (E = S, O) is reported. Addition of P(C4H3E)3 to [Ru3(CO)10(mu-dppm)] (1) at room temperature in the presence of Me3NO gives simple substitution products [Ru3(CO)9(mu-dppm)(P(C4H3E)3)] (E = S, 2; E = O, 3). Mild thermolysis in the presence of further Me3NO affords the thiophyne and furyne complexes [Ru3(CO)7(mu-dppm)(mu3-eta2-C4H2E)(mu-P(C4H3E)2)(mu-H)] (E = S, 4; E = O, 6) resulting from both carbon-hydrogen and carbon-phosphorus bond activation. In each the C4H2E (E = S, O) ligand donates 4-electrons to the cluster and the rings are tilted with respect to the mu-dppm and the phosphido-bridged open triruthenium unit. Heating 4 at 80 degrees C leads to the formation of the ring-opened cluster [Ru3(CO)5(mu-CO)(mu-dppm)(mu3-eta3-SC4H3)(mu-P(C4H3S)2)] (5) resulting from carbon-sulfur bond scission and carbon-hydrogen bond formation and containing a ring-opened mu3-eta3-1-thia-1,3-butadiene ligand. In contrast, a similar thermolysis of 3 affords the phosphinidene cluster [Ru3(CO)7(mu-dppm)(mu3-eta2-C4H2O)(mu3-P(C4H3O))] (7) resulting from a second phosphorus-carbon bond cleavage and (presumably) elimination of furan. Treatment of 4 and 6 with PPh3 affords the simple phosphine-substituted products [Ru3(CO)6(PPh3)(mu-dppm)(mu3-eta2-C4H2E)(mu-P(C4H3E)2)(mu-H)] (E = S, 8; E = O, 9). Both thiophyne and furyne clusters 4 and 6 readily react with hydrogen bromide to give [Ru3(CO)6Br(mu-Br)(mu-dppm)(mu3-eta2-eta1-C4H2E)(mu-P(C4H3E)2)(mu-H)] (E = S, 10; E = O, 11) containing both terminal and bridging bromides. Here the alkynes bind in a highly unsymmetrical manner with one carbon acting as a bridging alkylidene and the second as a terminally bonded Fisher carbene. As far as we are aware, this binding mode has only previously been noted in ynamine complexes or those with metals in different oxidation states. The crystal structures of seven of these new triruthenium clusters have been carried out, allowing a detailed analysis of the relative orientations of coordinated ligands.     

2,5-Dioxido-1,4-benzoquinonediimine (H2L2-), a hydrogen-bonding noninnocent bridging ligand related to aminated topaquinone: different oxidation state distributions in complexes [{(bpy)2Ru}2(mu-H2L)]n (n=0,+,2+,3+,4+) and [{(acac)2Ru}2(mu-H2L)]m (m=2-,-,0,+,2+)

The symmetrically dinuclear title compounds were isolated as diamagnetic [(bpy)2Ru(mu-H2L)Ru(bpy)2](ClO4)2 (1-(ClO4)2) and as paramagnetic [(acac)2Ru(mu-H2L)Ru(acac)2] (2) complexes (bpy=2,2'-bipyridine; acac- = acetylacetonate = 2,4-pentanedionato; H2L = 2,5-dioxido-1,4-benzoquinonediimine). The crystal structure of 22 H2O reveals an intricate hydrogen-bonding network: Two symmetry-related molecules 2 are closely connected through two NH(H2L2-)O(acac-) interactions, while the oxygen atoms of H2L2- of two such pairs are bridged by an (H2O)8 cluster at half-occupancy. The cluster consists of cyclic (H2O)6 arrangements with the remaining two exo-H2O molecules connecting two opposite sides of the cyclo-(H2O)6 cluster, and oxido oxygen atoms forming hydrogen bonds with the molecules of 2. Weak antiferromagnetic coupling of the two ruthenium(III) centers in 2 was established by using SQUID magnetometry and EPR spectroscopy. Geometry optimization by means of DFT calculations was carried out for 1(2+) and 2 in their singlet and triplet ground states, respectively. The nature of low-energy electronic transitions was explored by using time-dependent DFT methods. Five redox states were reversibly accessible for each of the complexes; all odd-electron intermediates exhibit comproportionation constants K(c)>10(8). UV-visible-NIR spectroelectrochemistry and EPR spectroscopy of the electrogenerated paramagnetic intermediates were used to ascertain the oxidation-state distribution. In general, the complexes 1n+ prefer the ruthenium(II) configuration with electron transfer occurring largely at the bridging ligand (mu-H2Ln-), as evident from radical-type EPR spectra for 13+ and (+. Higher metal oxidation states (iii, iv) appear to be favored by the complexes 2m; intense long-wavelength absorption bands and RuIII-type EPR signals suggest mixed-valent dimetal configurations of the paramagnetic intermediates 2+ and 2-.     

Terminal Hydride Species in [FeFe]‐Hydrogenases Are Vibrationally Coupled to the Active Site Environment

A combination of nuclear resonance vibrational spectroscopy (NRVS), FTIR spectroscopy, and DFT calculations was used to observe and characterize Fe?H/D bending modes in CrHydA1 [FeFe]‐hydrogenase Cys‐to‐Ser variant C169S. Mutagenesis of cysteine to serine at position 169 changes the functional group adjacent to the H‐cluster from a ‐SH to ‐OH, thus altering the proton transfer pathway. The catalytic activity of C169S is significantly reduced compared to that of native CrHydA1, presumably owing to less efficient proton transfer to the H‐cluster. This mutation enabled effective capture of a hydride/deuteride intermediate and facilitated direct detection of the Fe?H/D normal modes. We observed a significant shift to higher frequency in an Fe?H bending mode of the C169S variant, as compared to previous findings with reconstituted native and oxadithiolate (ODT)‐substituted CrHydA1. On the basis of DFT calculations, we propose that this shift is caused by the stronger interaction of the ‐OH group of C169S with the bridgehead ‐NH‐ moiety of the active site, as compared to that of the ‐SH group of C169 in the native enzyme.     

Bimetallic carbonyl thiolates as functional models for Fe-only hydrogenases

The anion [Fe(2)(S(2)C(3)H(6))(CN)(CO)(4)(PMe(3))](-) (2(-)) is protonated by sulfuric or toluenesulfonic acid to give HFe(2)(S(2)C(3)H(6))(CN)(CO)(4)(PMe(3)) (2H), the structure of which has the hydride bridging the Fe atoms with the PMe(3) and CN(-) trans to the same sulfur atom. (1)H, (13)C, and (31)P NMR spectroscopy revealed that HFe(2)(S(2)C(3)H(6))(CN)(CO)(4)(PMe(3)) is stereochemically rigid on the NMR time scale with four inequivalent carbonyl ligands. Treatment of 2(-) with (Me(3)O)BF(4) gave Fe(2)(S(2)C(3)H(6))(CNMe)(CO)(4)(PMe(3)) (2Me). The Et(4)NCN-induced reaction of Fe(2)(S(2)C(3)H(6))(CO)(6) with P(OMe)(3) gave [Fe(2)(S(2)C(3)H(6))(CN)(CO)(4)[P(OMe)(3)]](-) (4). Spectroscopic and electrochemical measurements indicate that 2H can be further protonated at nitrogen to give [HFe(2)(S(2)C(3)H(6))(CNH)(CO)(4)(PMe(3))](+) (2H(2)(+)). Electrochemical and analytical data show that reduction of 2H(2)(+) gives H(2) and 2(-). Parallel electrochemical studies on [HFe(2)(S(2)C(3)H(6))(CO)(4)(PMe(3))(2)](+) (3H(+)) in acidic solutions led also to catalytic proton reduction. The 3H(+)/3H couple is reversible, whereas the 2H(2)(+)/2H(2) couple is not, because of the efficiency of the latter as a proton reduction catalyst. Proton reduction is proposed to involve protonation of reduced diiron hydrides. DFT calculations establish that the regiochemistry of protonation is subtly dependent on the coligands but is more favorable to occur at the Fe-Fe bond for [Fe(2)(S(2)C(3)H(6))(CN)(CO)(4)(PMe(3))](-) than for [Fe(2)(S(2)C(3)H(6))(CN)(CO)(4)(PH(3))](-) or [Fe(2)(S(2)C(3)H(6))(CN)(CO)(4)[P(OMe)(3)]](-). The Fe(2)H unit stabilizes the conformer with eclipsed CN and PMe(3) because of an attractive electrostatic interaction between these ligands.     

FeMo cofactor of nitrogenase: a density functional study of states M(N), M(OX), M(R), and M(I).

The M(N) S = (3)/(2) resting state of the FeMo cofactor of nitrogenase has been proposed to have metal-ion valencies of either Mo(4+)6Fe(2+)Fe(3+) (derived from metal hyperfine interactions) or Mo(4+)4Fe(2+)3Fe(3+) (from M?ssbauer isomer shifts). Spin-polarized broken-symmetry (BS) density functional theory (DFT) calculations have been undertaken to determine which oxidation level best represents the M(N) state and to provide a framework for understanding its energetics and spectroscopy. For the Mo(4+)6Fe(2+)Fe(3+) oxidation state, the spin coupling pattern for several spin state alignments compatible with S = (3)/(2) were generated and assessed by energy and geometric criteria. The most likely BS spin state is composed of a Mo3Fe cluster with spin S(a) = 2 antiferromagnetically coupled to a 4Fe' cluster with spin S(b) = (7)/(2). This state has a low DFT energy for the isolated FeMoco cluster and the lowest energy when the interaction with the protein and solvent environment is included. This spin state also displays calculated metal hyperfine and M?ssbauer isomer shifts compatible with experiment, and optimized geometries that are in excellent agreement with the protein X-ray data. Our best model for the actual spin-coupled state within FeMoco alters this BS state by a slight canting of spins and is analogous in several respects to that found in the 8Fe P-cluster in the same protein. The spin-up and spin-down components of the LUMO contain atomic contributions from Mo(4+) and the homocitrate and from the central prismane Fe sites and muS(2) atoms, respectively. This qualitative picture of the accepting orbitals for M(N) is consistent with observations from M?ssbauer spectra of the one-electron reduced states. Similar calculations for the Mo(4+)4Fe(2+)3Fe(3+) oxidation state yield results that are in poorer agreement with experiment. Using the Mo(4+)6Fe(2+)Fe(3+) oxidation level as the most plausible resting state, the geometric, electronic and energetic properties of the one-electron redox transition to the oxidized state, M(OX), catalytically observed M(R) and radiolytically reduced M(I) states have also been explored.     

A counterintuitive structural effect of metal-metal bond protonation and its electronic underpinnings

Protonation across the metal-metal bond in the complexes [(CO)(2)M(mu-dppm)(mu-PtBu(2))(mu-H)M(CO)(2)] (M=Fe or Ru, dppm=Ph(2)PCH(2)PPh(2)) induces M-M bond shortening of up to about 0.05 A. DFT calculations on simplified iron models reproduce this trend well. Conversely, the computations show that the M-M distance in the dimer [{Cp*Ir(CO)}(2)] lengthens with two consecutive protonations, but there are no crystal structure determinations to highlight the effects on the Ir-Ir bond. DFT calculations and the analogous cobalt system confirm that the transformation of a two-electron, two-center (2e-2c) bond into a 2e-3c bond is accompanied by the predicted elongation. An MO analysis indicated similar nature and evolution of the M-M bonding these cases. In particular, the HOMOs of the mono-hydrido cations [Cp(CO)M(mu-H)M(CO)Cp](+) (M=Ir, Co) have evident M-M bent-bond character, and hence subsequent protonation invariably causes a decrease in the bond index. The Fe(2) and Co(2) systems have also been analyzed with the quantum theory of atoms in molecules (QTAIM) method, but in no case was an M-M bond critical point located unless an artificially shorter M-M distance was imposed. However, the trends for the atoms-in-molecules (AIM) bond delocalization indexes delta(M-M) confirm the overall M-M bond weakening on protonation. In conclusion, all the computational results for the iron system indicate that the paradigm of a direct correlation between bond strength and distance is not always applicable. This is attributable to a very flat potential energy surface and various competing effects imposed by the bridging ligands.     

Sulfur K-edge XAS and DFT calculations on [Fe4S4]2+ clusters: effects of H-bonding and structural distortion on covalency and spin topology

Sulfur K-edge X-ray absorption spectroscopy of a hydrogen-bonded elongated [Fe4S4]2+ cube is reported. The data show that this synthetic cube is less covalent than a normal compressed cube with no hydrogen bonding. DFT calculations reveal that the observed difference in electronic structure has significant contributions from both the cluster distortion and from hydrogen bonding. The elongated and compressed Fe4S4 structures are found to have different spin topologies (i.e., orientation of the delocalized Fe2S2 subclusters which are antiferromagnetically coupled to each other). It is suggested that the H-bonding interaction with the counterion does not contribute to the cluster elongation. A magneto-structural correlation is developed for the Fe4S4 cube that is used to identify the redox-active Fe2S2 subclusters in active sites of HiPIP and ferredoxin proteins involving these clusters.