Computational studies of [NiFe] and [FeFe] hydrogenases

Fe-H/D stretching and bending modes in nuclear resonant vibrational, Raman and infrared spectroscopies: comparisons of density functional theory and experiment

Infrared, Raman, and nuclear resonant vibrational (NRVS) spectroscopies have been used to address the Fe-H bonding in trans-HFe(CO) iron hydride compound, HFe(CO)(dppe)2, dppe = 1,2-bis(diphenylphosphino)ethane. H and D isotopomers of the compound, with selective substitution at the metal-coordinated hydrogen, have been considered in order to address the Fe-H/D stretching and bending modes. Experimental results are compared to the normal mode analysis by density functional theory (DFT). The results are that (i) the IR spectrum does not clearly show Fe-H stretching or bending modes; (ii) Fe-H stretching modes are clear but weak in the Raman spectrum, and Fe-H bending modes are weak; (iii) NRVS 57Fe spectroscopy resolves Fe-H bending clearly, but Fe-H or Fe-D stretching is above its experimentally resolved frequency range. DFT calculations (with no scaling of frequencies) show intensities and peak locations that allow unambiguous correlations between observed and calculated features, with frequency errors generally less than 15 cm(-1). Prospects for using these techniques to unravel vibrational modes of protein active sites are discussed.

Role of the azadithiolate cofactor in models for [FeFe]-hydrogenase: novel structures and catalytic implications

This paper summarizes studies on the redox behavior of synthetic models for the [FeFe]-hydrogenases, consisting of diiron dithiolato carbonyl complexes bearing the amine cofactor and its N-benzyl derivative. Of specific interest are the causes of the low reactivity of oxidized models toward H(2), which contrasts with the high activity of these enzymes for H(2) oxidation. The redox and acid-base properties of the model complexes [Fe(2)[(SCH(2))(2)NR](CO)(3)(dppv)(PMe(3))](+) ([2](+) for R = H and [2'](+) for R = CH(2)C(6)H(5), dppv = cis-1,2-bis(diphenylphosphino)ethylene)) indicate that addition of H(2) followed by deprotonation are (i) endothermic for the mixed valence (Fe(II)Fe(I)) state and (ii) exothermic for the diferrous (Fe(II)Fe(II)) state. The diferrous state is shown to be unstable with respect to coordination of the amine to Fe, a derivative of which was characterized crystallographically. The redox and acid-base properties for the mixed valence models differ strongly for those containing the amine cofactor versus those derived from propanedithiolate. Protonation of [2'](+) induces disproportionation to a 1:1 mixture of the ammonium [H2'](+) (Fe(I)Fe(I)) and the dication [2'](2+) (Fe(II)Fe(II)). This effect is consistent with substantial enhancement of the basicity of the amine in the Fe(I)Fe(I) state vs the Fe(II)Fe(I) state. The Fe(I)Fe(I) ammonium compounds are rapid and efficient H-atom donors toward the nitroxyl compound TEMPO. The atom transfer is proposed to proceed via the hydride. Collectively, the results suggest that proton-coupled electron-transfer pathways should be considered for H(2) activation by the [FeFe]-hydrogenases.

Potential hydrogen bottleneck in nickel-iron hydrogenase

The role of two-state reactivity at the enzyme active site with respect to binding of molecular H(2) for the high- and low-spin of [NiFe] hydrogenase (Ni-SI forms) is examined by density functional theory. In addition to examination of a single H(2) molecule binding at either the Ni or Fe of the active site, the possibility that H(2) binds simultaneously at each metal center in the active site of this enzyme is examined. The concurrent binding of two molecules of H(2) suggests a potential hydrogen bottleneck in which high concentrations might lead to a decrease in the rate of hydrogen oxidation.

A unifying structural and electronic concept for Hmd and [FeFe] hydrogenase active sites

The hydrogenases [FeFe] and Hmd feature at first sight rather different active sites. A closer inspection reveals striking similarities, which allow us to define swapped ligand spheres in such a way that the single active iron center of Hmd functions in a first-shell ligand environment resembling the reacting iron atom in [FeFe] hydrogenase and vice versa. These redesigned ligand environments can be conveniently studied with quantum chemical methods and point to general reactivity principles for iron centers with hydrogenase activity.

[3+2+1] cycloaddition involving alkynes, CO and bridging vinyliminium ligands in diiron complexes: a dinuclear version of the Dötz reaction?

The vinyliminium complex [Fe2{mu-eta(1):eta(3)-C(SiMe3)=CHC=N(Me)2}(mu-CO)(CO)(Cp)2][SO3CF3] reacts with HC[triple bond]CR (R = CPh2OH), affording a mixture of the 2,4,6-trisubstituted oxo-eta(5)-cyclohexadienyl complex [Fe{eta(5)-C6H2O(NMe2)(SiMe3)(R)}(Cp)], the 2,6-disubstituted phenol C6H3R(NMe2)OH (R = CHPh2) and 1,2,4-trisubstituted ferrocene [1-NMe2-2-R-4-SiMe3-Fc]. The corresponding reaction with HC[triple bond]CR (R = CMe2OH) yields analogous products: [Fe{eta(5)-C6H2O(NMe2)(SiMe3)(R')}(Cp)] (R' = CMe=CH2), the phenol C6H3R'(NMe2)OH together with [1-NMe2-2-R-4-SiMe3-Fc].

Isomerization of the hydride complexes [HFe2(SR)2(PR3)(x)(CO)(6-x)]+ (x = 2, 3, 4) relevant to the active site models for the [FeFe]-hydrogenases

The stepwise formation of bridging (mu-) hydrides of diiron dithiolates is discussed with attention on the pathway for protonation and subsequent isomerizations. Our evidence is consistent with protonations occurring at a single Fe center, followed by isomerization to a series of mu-hydrides. Protonation of Fe(2)(edt)(CO)(4)(dppv) (1) gave a single mu-hydride with dppv spanning apical and basal sites, which isomerized at higher temperatures to place the dppv into a dibasal position. Protonation of Fe(2)(pdt)(CO)(4)(dppv) (2) followed an isomerization pathway similar to that for [1H](+), except that a pair of isomeric terminal hydrides were observed initially, resulting from protonation at the Fe(CO)(3) or Fe(CO)(dppv) site. The first observable product from low temperature protonation of the tris-phosphine Fe(2)(edt)(CO)(3)(PMe(3))(dppv) (3) was a single mu-hydride wherein PMe(3) is apical and the dppv ligand spans apical and basal sites. Upon warming, this isomer converted fully but in a stepwise manner to a mixture of three other isomeric hydrides. Protonation of Fe(2)(pdt)(CO)(3)(PMe(3))(dppv) (4) proceeded similarly to the edt analogue 3, however a terminal hydride was observed, albeit only briefly and at very low temperatures (-90 degrees C). Low-temperature protonation of the bis-chelates Fe(2)(xdt)(CO)(2)(dppv)(2) produced exclusively the terminal hydrides [HFe(2)(xdt)(mu-CO)(CO)(dppv)(2)](+) (xdt = edt and pdt), which subsequently isomerized to a pair of mu-hydrides. At room temperature these (dppv)(2) derivatives convert to an equilibrium of two isomers, one C(2)-symmetric and the other C(s)-symmetric. The stability of the terminal hydrides correlates with the (C(2)-isomer)/(C(s)-isomer) equilibrium ratio, which reflects the size of the dithiolate. The isomerization was found to be unaffected by the presence of excess acid, by solvent polarity, and the presence of D(2)O. This isomerization mechanism is proposed to be intramolecular, involving a 120 degrees rotation of the HFeL(3) subunit to an unobserved terminal basal hydride as the rate-determining step. The observed stability of the hydrides was supported by DFT calculations, which also highlight the instability of the basal terminal hydrides. Isomerization of the mu-hydride isomers occurs on alternating FeL(3) via 120 degree rotations without generating D(2)O-exchangeable intermediates.

The useful properties of H2O as a ligand of a hydrogenase mimic

This paper investigates the required properties of Ru-coordinated ligands of a Ni-Ru based hydrogenase mimic. A series of ligands, including MeCN, pyridine, H(2)O and OH(-) were coordinated to Ru, with H(2)O being the only ligand to promote H(2)-activation. In addition, a tethered pyridyl moiety was synthesised and found to completely inhibit H(2)-activation. We conclude, therefore, that H(2)O is the ideal ligand for this mimic as a result of both its mild basicity and the availability of two lone pairs for simultaneous binding to Ru and H(2).

Comparative analysis of the performance of commonly available density functionals in the determination of geometrical parameters for zinc complexes

A set of 44 Zinc-ligand bond-lengths and of 60 ligand-metal-ligand bond angles from 10 diverse transition-metal complexes, representative of the coordination spheres of typical biological Zn systems, were used to evaluate the performance of a total of 18 commonly available density functionals in geometry determination. Five different basis sets were considered for each density functional, namely two all-electron basis sets (a double-zeta and triple-zeta formulation) and three basis sets including popular types of effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden. The results show that there are presently several better alternatives to the popular B3LYP density functional for the determination of Zn-ligand bond-lengths and angles. BB1K, MPWB1K, MPW1K, B97-2 and TPSS are suggested as the strongest alternatives for this effect presently available in most computational chemistry software packages. In addition, the results show that the use of effective-core potentials (in particular Stuttgart-Dresden) has a very limited impact, in terms of accuracy, in the determination of metal-ligand bond-lengths and angles in Zinc-complexes, and is a good and safe alternative to the use of an all-electron basis set such as 6-31G(d) or 6-311G(d,p).

Concerto catalysis--harmonising [NiFe]hydrogenase and NiRu model catalysts

This communication reports the successful merging of the chemical properties of a natural [NiFe]hydrogenase (Desulfovibrio vulgaris Miyazaki F) and our previously reported [NiRu] hydrogenase-mimic. The catalytic activity of both the natural enzyme and the mimic is almost identical, with the exception of working pH ranges, and this allows us to use them simultaneously in the same reaction flask. In such a manner, isotope exchange between D(2) and H(2)O could be conducted over an extended pH range (about 2-10) in one pot under mild conditions at ambient temperature and pressure.

Spectroelectrochemical characterization of the active site of the [FeFe] hydrogenase HydA1 from Chlamydomonas reinhardtii

Hydrogenases catalyze the reversible oxidation of molecular hydrogen. The active site of the [FeFe] hydrogenases (H-cluster) contains a catalytically active binuclear subcluster ([2Fe](H)) connected to a "cubane" [4Fe4S](H) subcluster. Here we present an IR spectroelectrochemical study of the [FeFe] hydrogenase HydA1 isolated from the green alga Chlamydomonas reinhardtii. The enzyme shows IR bands similar to those observed for bacterial [FeFe] hydrogenases. They are assigned to the stretching vibrations of the CN(-) and CO ligands on both irons of the [2Fe](H) subcluster. By following changes in frequencies of the IR bands during electrochemical titrations, two one-electron redox processes of the active enzyme could be distinguished. The reduction of the oxidized state (H(ox)) occurred at a midpoint potential of -400 mV vs NHE (H(ox)/H(red) transition) and relates to a change of the formal oxidation state of the binuclear subcluster. A subsequent reduction (H(red)/H(sred) transition) was determined to have a midpoint potential of -460 mV vs NHE. On the basis of the IR spectra, it is suggested that the oxidation state of the binuclear subcluster does not change in this transition. Tentatively, a reduction of the [4Fe4S](H) cluster has been proposed. In contrast to the bacterial [FeFe] hydrogenases, where the bridging CO ligand becomes terminal when going from H(ox) to H(red), in HydA1 the bridging CO is present in both the H(ox) and H(red) state. The removal of the bridging CO moiety has been observed in the H(red) to H(sred) transition. The significance of this result for the hydrogen conversion mechanism of this class of enzymes is discussed.

Sulfur oxygenates of biomimetics of the diiron subsite of the [FeFe]-hydrogenase active site: properties and oxygen damage repair possibilities

This study explores the site specificity (sulfur vs the Fe-Fe bond) of oxygenation of diiron (Fe(I)Fe(I) and Fe(II)Fe(II)) organometallics that model the 2-iron subsite in the active site of [FeFe]-hydrogenase: (mu-pdt)[Fe(CO)(2)L][Fe(CO)(2)L'] (L = L' = CO (1); L = PPh(3), L' = CO (2); L = L' = PMe(3) (4)) and (mu-pdt)(mu-H)[Fe(CO)(2)PMe(3)](2) (5). DFT computations find that the Fe-Fe bond in the Fe(I)Fe(I) diiron models is thermodynamically favored to produce the mu-oxo or oxidative addition product, Fe(II)-O-Fe(II); nevertheless, the sulfur-based HOMO-1 accounts for the experimentally observed mono- and bis-O-atom adducts at sulfur, i.e., (mu-pst)[Fe(CO)(2)L][Fe(CO)(2)L'] (pst = -S(CH(2))(3)S(O)-, 1,3-propanesulfenatothiolate; L = L' = CO (1-O); L = PPh(3), L' = CO (2-O); L = L' = PMe(3) (4-O)) and (mu-pds)[Fe(CO)(2)L][Fe(CO)(2)L'] (pds = -(O)S(CH(2))(3)S(O)-, 1,3-propanedisulfenato; L = PPh(3), L' = CO (2-O(2))). The Fe(II)(mu-H)Fe(II) diiron model (5), for which the HOMO is largely of sulfur character, exclusively yields S-oxygenation. The depressing effect of such bridging ligand modification on the dynamic NMR properties arising from rotation of the Fe(CO)(3) correlates with higher barriers to the CO/PMe(3) exchange of (mu-pst)[Fe(CO)(3)](2) as compared to (mu-pdt)[Fe(CO)(3)](2). Five molecular structures are confirmed by X-ray diffraction: 1-O, 2-O, 2-O(2), 4-O, and 6. Deoxygenation with reclamation of the mu-pdt parent complex occurs in a proton/electron-coupled process. The possible biological relevance of oxygenation and deoxygenation studies is discussed.

Dithiolate-bridged Fe-Ni-Fe trinuclear complexes consisting of Fe(CO)(3-n)(CN)(n) (n = 0, 1) components relevant to the active site of [NiFe] hydrogenase

A dithiolate-bridged Fe-Ni-Fe trinuclear carbonyl complex [(CO)(3)Fe(mu-ndt)Ni(mu-ndt)Fe(CO)(3)] (1, ndt = norbornane-exo-2,3-dithiolate) has been synthesized from the reaction of [Fe(CO)(4)I(2)] and Li(2)[Ni(ndt)(2)]. This reaction was found to occur with concomitant formation of a tetranuclear cluster [Ni(3)(mu-ndt)(4)FeI] (2). Treatment of 1 with Na[N(SiMe(3))(2)] transforms some of the CO ligands into CN(-), and the monocyanide complex (PPh(4))[(CO)(2)(CN)Fe(mu-ndt)Ni(mu-ndt)Fe(CO)(3)] (3) and the dicyanide complex (PPh(4))(2)[(CO)(2)(CN)Fe(mu-ndt)Ni(mu-ndt)Fe(CO)(2)(CN)] (4) were isolated. X-ray structural analyses of the trinuclear complexes revealed a Fe-Ni-Fe array in which the metal centers are connected by the ndt sulfur bridges and direct Fe-Ni bonds. Hydrogen bonding between the CN ligand in 3 and cocrystallized ethanol was found in the solid-state structure. The monocyanide complex 3 and dicyanide complex 4 reacted with acids such as HOTf or HCl generating insoluble materials, whereas complex 1 did not react.

(14)N HYSCORE investigation of the H-cluster of [FeFe] hydrogenase: evidence for a nitrogen in the dithiol bridge

Hydrogenases are enzymes catalyzing the reversible heterolytic splitting of molecular hydrogen. Despite extensive investigations of this class of enzymes its catalytic mechanism is not yet well understood. In this paper spectroscopic investigations of the active site of [FeFe] hydrogenase are presented. The so-called H-cluster consists of a bi-nuclear catalytically active subcluster connected to a [4Fe4S] ferredoxin-like unit via a Cys-thiol bridge. An important feature of the H-cluster is that both irons in the bi-nuclear subcluster are coordinated by CN and CO ligands. The bi-nuclear site also carries a dithiol bridge, whose central atom has not yet been identified. Nitrogen and oxygen are the most probable candidates from a mechanistic point of view. Here we present a study of the (14)N nuclear quadrupole and hyperfine interactions of the active oxidized state of the H-cluster using advanced EPR methods. In total three (14)N nuclei with quadrupole couplings of 0.95 MHz, 0.35 MHz and 1.23 MHz were detected using hyperfine sublevel correlation spectroscopy (HYSCORE). The assignment of the signals is based on their (14)N quadrupole couplings in combination with DFT calculations. One signal is assigned to the CN ligand of the distal iron, one to a Lys side chain nitrogen and one to the putative nitrogen of the dithiol bridge. Hence, these results provide the first experimental evidence for a di-(thiomethyl)amine ligand (-S-CH(2)-NH-CH(2)-S-) in the bi-nuclear subcluster. This finding is important for understanding the mechanism of [FeFe] hydrogenases, since the nitrogen is likely to act as an internal base facilitating the heterolytic splitting/formation of H(2).