Synthetic and structural studies on [Fe2(SR)2(CN)x(CO)6-x](x-) as active site models for Fe-only hydrogenases

Synthesis and characterization of diiron dithiolate complexes containing a quinoxaline bridge

A potential model complex for the hydrogenase active site, [Fe(2){(μ-CH(2)S)(2)R}(CO)(6)] (1) (R = quinoxaline), was synthesized by condensation of [(μ-LiS)(2)Fe(2)(CO)(6)] with 2,3-bis(bromomethyl)quinoxaline. Reactions of 1 with bis(diphenylphosphino)methane (dppm) under a range of conditions yielded substituted complexes [Fe(2){(μ-CH(2)S)(2)R}(CO)(5)(dppm)] (2), [Fe(2){(μ-CH(2)S)(2)R}(CO)(4)(k(2)-dppm)] (3) and [Fe(2){(μ-CH(2)S)(2)R}(CO)(4)(μ-dppm)] (4). X-ray crystallography confirms that in 2, the dppm is terminally bonded to an iron atom via one phosphorus atom, whereas in 3, it acts as a chelating ligand to coordinate to an iron center in a dibasal-substituted manner. In 4, the dppm bridges the two iron atoms in a cis basal/basal fashion with one phosphorus bonded to each iron atom. Treatment of 1 with various tertiary phosphines at room temperature in acetonitrile (MeCN) generates a range of mono-substituted products [Fe(2){(μ-CH(2)S)(2)R}(CO)(5)L] (5, L = PEt(3); 6, PMe(3); 7, PPh(3); 8, Me(2)PPh). With Bu(t)NC, mono- and di-substituted [Fe(2){(μ-CH(2)S)(2)R}(CO)(5)(Bu(t)NC)] (9) and [Fe(2){(μ-CH(2)S)(2)R}(CO)(4)(Bu(t)NC)(2)] (10) complexes are generated. All the complexes were characterized by elemental analysis, IR, MS and NMR spectroscopy. IR and NMR spectroscopic studies suggest that addition of excess HBF(4)·OEt(2) acid to 1-4 led to the protonation of quinoxaline nitrogen atoms. In contrast, 5-10 were not stable in acidic media. Electrochemistry of 1-4 was investigated in the acetonitrile medium (0.1 M Bu(4)NPF(6)). The electrochemical instability of the reduced ligand, quinoxaline, and the reduced forms of these complexes revealed from the electrochemical studies suggests that they do not provide ideal models of the hydrogenase active site.

Targeting intermediates of [FeFe]-hydrogenase by CO and CN vibrational signatures

In this work, we employ density functional theory to assign vibrational signatures of [FeFe]-hydrogenase intermediates to molecular structures. For this purpose, we perform an exhaustive analysis of structures and harmonic vibrations of a series of CN and CO containing model clusters of the [FeFe]-hydrogenase enzyme active site considering also different charges, counterions, and solvents. The pure density functional BP86 in combination with a triple-ζ polarized basis set produce reliable molecular structures as well as harmonic vibrations. Calculated CN and CO stretching vibrations are analyzed separately. Scaled vibrational frequencies are then applied to assign intermediates in [FeFe]-hydrogenase's reaction cycle. The results nicely complement the previous studies of Darensbourg and Hall, and Zilberman et al. The infrared spectrum of the H(ox) form is in very good agreement with the calculated spectrum of the Fe(I)Fe(II) model complex featuring a free coordination site at the distal Fe atom, as well as, with the calculated spectra of the complexes in which H(2) or H(2)O are coordinated at this site. The spectrum of H(red) measured from Desulfovibrio desulfuricans is compatible with a mixture of a Fe(I)Fe(I) species with all terminal COs, and a Fe(I)Fe(I) species with protonated dtma ligand, while the spectrum of H(red) recently measured from Chlamydomonas reinhardtii is compatible with a mixture of a Fe(I)Fe(I) species with a bridged CO, and a Fe(II)Fe(II) species with a terminal hydride bound to the Fe atom.

Favorable Protonation of the (μ-edt)[Fe(2)(PMe(3))(4)(CO)(2)(H-terminal)](+) Hydrogenase Model Complex Over Its Bridging μ-H Counterpart: A Spectroscopic and DFT Study

The mechanism of hydrogen production in [FeFe] hydrogenase remains elusive. However, a species featuring a terminal hydride bound to the distal Fe is thought to be the key intermediate leading to hydrogen production. In this study, density functional theory (DFT) calculations on the terminal (H-term) and bridging (μ-H) hydride isomers of (μ-edt)-[Fe(2)(PMe(3))(4)(CO)(2)H](+) are presented in order to understand the factors affecting their propensity for protonation. Relative to H-term, μ-H is 12.7 kcal/mol more stable, which contributes to its decreased reactivity towards an acid. Potential energy surface (PES) calculations for the reaction of the H-term isomer with 4-nitropyridinium, a proton source, further reveal a lower activation energy barrier (14.5 kcal/mol) for H-term than for μ-H (29 kcal/mol). Besides these energetic considerations, the H-term isomer displays a key molecular orbital (MO <139>) that has a relatively strong hydride (1s) contribution (23%), which is not present in the μ-H isomer. This indicates a potential orbital control of the reaction of the hydride complexes with acid. The lower activation energy barrier and this key MO together control the overall catalytic activity of (μ-edt)[Fe(2)(PMe(3))(4)(CO)(2)(H-term)](+). Lastly, Raman and IR spectroscopy were performed in order to probe the ν(Fe-H) stretching mode of the two isomers and their deuterated counterparts. A ν(Fe-H) stretching mode was observed for the μ-H complex at 1220 cm(-1). However, the corresponding mode is not observed for the less stable H-term isomer.

On understanding proton transfer to the biocatalytic [Fe-Fe](H) sub-cluster in [Fe-Fe]H(2)ases: QM/MM MD simulations

Proton transfer to the [Fe-Fe](H) sub-cluster in the Desulfovibrio desulfuricans (DdH) and Clostridium pasteurianum (CpI) [Fe-Fe] hydrogenases was investigated by a combination of first principles and empirical molecular dynamics simulations. Pathways that can be inferred from the X-ray crystal structures of DdH and CpI, i.e., (Glu159→Ser198→Glu156→water460→Cys178→DTMA([Fe-Fe](H)) and (Glu282→Ser319→Glu279→water612→Cys299), respectively, were considered. Proton transfer from Cys178 to DTMA in the [Fe-Fe](H) sub-cluster in DdH was readily observed in our results, specifically when [Fe-Fe](H) was in the reduced state ([Fe(I)-Fe(I)]) or in the mixed valence state for the protonated distal iron Fe(d) ([Fe(I)-Fe(II)-H(-)](H)). A concerted mechanism is proposed, where proton transfer in DdH from Glu159 to Glu156 via Ser198 and Glu156 to Cys178 via water460 readily occurred, as well as from Glu282 to Glu279 via Ser319 and Glu279 to Cys299 via water612 in CpI. The theoretical prediction of the proton transfer characteristics is consistent with the assumed biocatalytic mechanism of the [Fe-Fe] hydrogenases in which the proton binds at Fe(d), providing confirmation that has not been explored so far. The computational results were qualitatively validated by the agreement with experimental hydrogen production activity data for mutated CpI enzymes, relative to the wild-type protein. Finally, the insight provided by the simulations, combined, in part, with experimental validation, are important for establishing an approach in future exploration of proton transfer to the active site in this class of enzymes, and possibly also for biomimetic analogs.

Vibrational analysis of the model complex (mu-edt)[Fe(CO)(3)](2) and comparison to iron-only hydrogenase: the activation scale of hydrogenase model systems

Research on simple [FeFe] hydrogenase model systems of type (mu-S(2)R)[Fe(CO)(3)](2) (R = C(2)H(4) (edt), C(3)H(6) (pdt)) which have been shown to function as robust electrocatalysts for proton reduction, provides a reference to understand the electronic and vibrational properties of the active site of [FeFe] hydrogenases and of more sophisticated model systems. In this study, the solution and solid state Raman spectra of (mu-edt)[Fe(CO)(3)](2) and of the corresponding (13)CO-labeled complex are presented and analyzed in detail, with focus on the nu(C=O) and nu(Fe-CO)/delta(Fe-C=O) vibrational regions. These regions are specifically important as vibrations involving CO ligands serve as probes for the "electron richness" of low-valent transition metal centers and the geometric structures of the complexes. The obtained vibrational spectra have been completely assigned in terms of the nu(C=O), nu(Fe-CO), and delta(Fe-C=O) modes, and the force constants of the important C=O and Fe-CO bonds have been determined using our Quantum Chemistry Centered Normal Coordinate Analysis (QCC-NCA). In the 400-650 cm(-1) region, fifteen mixed nu(Fe-CO)/delta(Fe-C=O) modes have been identified. The most prominent Raman peaks at 454, 456, and 483 cm(-1) correspond to a combination of nu(Fe-CO) stretching and delta(Fe-C=O) linear bending modes. The less intense peaks at 416 cm(-1) and 419 cm(-1) correspond to pure delta(Fe-C=O) linear bends. In the nu(C=O) region, the nu(C=O) normal modes at lower energy (1968 and 1964 cm(-1)) are almost pure equatorial (eq) nu(C=O)(eq) stretching vibrations, whereas the remaining four nu(C=O) normal modes show dominant (C=O)(eq) (2070 and 1961 cm(-1)) and (C=O)(ax) (2005 and 1979 cm(-1); ax = axial) contributions. Importantly, an inverse correlation between the f(C=O)(ax/eq) and f(Fe-CO)(ax/eq) force constants is obtained, in agreement with the idea that the Fe(I)-CO bond in these types of complexes is dominated by pi-backdonation. Compared to the reduced form of [FeFe] hydrogenase (H(red)), the nu(C=O) vibrational frequencies of (mu-edt)[Fe(CO)(3)](2) are higher in energy, indicating that the dinuclear iron core in (mu-edt)[Fe(CO)(3)](2) is less electron rich compared to H(red) in the actual enzyme. Finally, quantum yields for the photodecomposition of (mu-edt)[Fe(CO)(3)](2) have been determined.