Computational studies of [NiFe] and [FeFe] hydrogenases

Computational Studies of Bioorganometallic Enzymes and Cofactors

Because of their complex geometric and electronic structures, the active sites and cofactors of bioorganometallic enzymes, which are characterized by their metal–carbon bonds, pose a major challenge for computational chemists. However, recent progress in computer technology and theoretical chemistry, along with insights gained from mechanistic, spectroscopic, and X-ray crystallographic studies, have established an excellent foundation for the successful completion of computational studies aimed at elucidating the electronic structures and catalytic cycles of these species. This chapter briefly reviews the most popular computational approaches employed in theoretical studies of bioorganometallic species and summarizes important information obtained from computational studies of (i) the enzymatic formation and cleavage of the Co–C bond of coenzyme B12; (ii) the catalytic cycle of methyl-coenzyme M reductase and its nickel-containing cofactor F430; (iii) the polynuclear active-site clusters of the bifunctional enzyme carbon monoxide dehydrogenase/acetyl-coenzyme A synthase; and (iv) the magnetic properties of the active-site cluster of Fe-only hydrogenases.

Small molecule mimics of hydrogenases: hydrides and redox

This tutorial review is aimed at chemical scientists interested in understanding and exploiting the remarkable catalytic behavior of the hydrogenases. The key structural features are analyzed for the active sites of the two most important hydrogenases. Reactivity is emphasized, focusing on mechanism and catalysis. Through this analysis, gaps are identified in the synthesis of functional replicas of these fascinating and potentially useful enzymes.

A structure-consistent mechanism for dioxygen formation in photosystem II

In recent DFT studies a new mechanism for O-O bond formation at the oxygen evolving center (OEC) in photosystem II has been suggested. With the structure of the S(4) state required for that mechanism, the structures of the lower S states are investigated herein by adding protons and electrons. A model was used including the full amino acids for the ones ligating the OEC, and in which the backbone positions were held fixed from the X-ray structure. The only charged second-shell ligand Arg357 was also included. An optimized structure for the S(1) state was reached with a large similarity to one of those suggested by EXAFS. A full catalytic cycle was derived which can rationalize the structural relaxation in the S(2) to S(3) transition, and the fact that only an electron leaves in the transition before. Water is suggested to bind to the OEC in the S(2) to S(3), and S(4) to S(0) transitions. A new possibility for water exchange is suggested from the final energy diagram. The optimal O-O bond formation occurs between an oxygen radical and an oxo ligand. The alternative mechanism, where the oxygen radical reacts with an external water, has a barrier about 20 kcal mol(-1) higher.

Nitrosyl derivatives of diiron(I) dithiolates mimic the structure and Lewis acidity of the [FeFe]-hydrogenase active site

This study probes the impact of electronic asymmetry of diiron(I) dithiolato carbonyls. Treatment of Fe2(S2C(n)H(2n))(CO)(6-x)(PMe3)x compounds (n = 2, 3; x = 1, 2, 3) with NOBF4 gave the derivatives [Fe2(S2C(n)H(2n))(CO)(5-x)(PMe3)x(NO)]BF4, which are electronically unsymmetrical because of the presence of a single NO(+) ligand. Whereas the monophosphine derivative is largely undistorted, the bis(PMe3) derivatives are distorted such that the CO ligand on the Fe(CO)(PMe3)(NO)(+) subunit is semibridging. Two isomers of [Fe2(S2C3H6)(CO)3(PMe3)2(NO)]BF4 were characterized spectroscopically and crystallographically. Each isomer features electron-rich Fe(CO)2PMe3 and electrophilic Fe(CO)(PMe3)(NO)(+) subunits. These species are in equilibrium with an unobserved isomer that reversibly binds CO (DeltaH = -35 kJ/mol, DeltaS = -139 J mol(-1) K(-1)) to give the symmetrical adduct [Fe2(S2C3H6)(mu-NO)(CO)4(PMe3)2]BF4. In contrast to Fe2(S2C3H6)(CO)4(PMe3)2, the bis(PMe3) nitrosyl complexes readily undergo CO substitution to give the (PMe3)3 derivatives. The nitrosyl complexes reduce at potentials that are approximately 1 V milder than their carbonyl counterparts. Results of density functional theory calculations, specifically natural bond orbital analysis, reinforce the electronic resemblance of the nitrosyl complexes to the corresponding mixed-valence diiron complexes. Unlike other diiron dithiolato carbonyls, these species undergo reversible reductions at mild potentials. The results show that the novel structural and chemical features associated with mixed-valence diiron dithiolates (the so-called H(ox) models) can be replicated in the absence of mixed-valency by the introduction of electronic asymmetry.

Precursors to [FeFe]-hydrogenase models: syntheses of Fe2(SR)2(CO)6 from CO-free iron sources

This report describes routes to iron dithiolato carbonyls that do not require preformed iron carbonyls. The reaction of FeCl 2, Zn, and Q 2S 2C n H 2 n (Q (+) = Na (+), Et 3NH (+)) under an atmosphere of CO affords Fe 2(S 2C n H 2 n )(CO) 6 ( n = 2, 3) in yields >70%. The method was employed to prepare Fe 2(S 2C 2H 4)( (13)CO) 6. Treatment of these carbonylated mixtures with tertiary phosphines, instead of Zn, gave the ferrous species Fe 3(S 2C 3H 6) 3(CO) 4(PR 3) 2, for R = Et, Bu, and Ph. Like the related complex Fe 3(SPh) 6(CO) 6, these compounds consist of a linear arrangement of three conjoined face-shared octahedral centers. Omitting the phosphine but with an excess of dithiolate, we obtained the related mixed-valence triiron species [Fe 3(S 2C n H 2 n ) 4(CO) 4] (-). The highly reducing all-ferrous species [Fe 3(S 2C n H 2 n ) 4(CO) 4] (2-) is implicated as an intermediate in this transformation. Reactive forms of iron, prepared by the method of Rieke, also combined with dithiols under a CO atmosphere to give Fe 2(S 2C n H 2 n )(CO) 6 in modest yields under mild conditions. Studies on the order of addition indicate that ferrous thiolates are formed prior to the onset of carbonylation. Crystallographic characterization demonstrated that the complexes Fe 3(S 2C 3H 6) 3(CO) 4(PEt 3) 2 and PBnPh 3[Fe 3(S 2C 3H 6) 4(CO) 4] feature high-spin ferrous and low-spin ferric as the central metal, respectively.

Diiron dithiolato carbonyls related to the H(ox)CO state of [FeFe]-hydrogenase

Oxidation of the electron-rich (E(1/2) = -175 vs Ag/AgCl) ethanedithiolato complex Fe2(S2C2H4)(CO)2(dppv)2 (1) under a CO atmosphere yielded [Fe2(S2C2H4)(mu-CO)(CO)2(dppv)2](+) ([1(CO)](+)), a model for the H(ox)(CO) state of the [FeFe]-hydrogenases. This complex exists as two isomers: a kinetically favored unsymmetrical derivative, unsym-[1(CO)](+), and a thermodynamically favored isomer, sym-[1(CO)](+), wherein both diphosphines span apical and basal sites. Crystallographic characterization of sym-[1(CO)](+) confirmed a C2-symmetric structure with a bridging CO ligand and an elongated Fe-Fe bond of 2.7012(14) A, as predicted previously. Oxidation of sym-[1(CO)](+) and unsym-[1(CO)](+) again by 1e(-) oxidation afforded the respective diamagnetic diferrous derivatives where the relative stabilities of the sym and unsym isomers are reversed. DFT calculations indicate that the stabilities of sym and unsym isomers are affected differently by the oxidation state of the diiron unit: the mutually trans CO ligands in the sym isomer are more destabilizing in the mixed-valence state than in the diferrous state. EPR analysis of mixed-valence complexes revealed that, for [1](+), the unpaired spin is localized on a single iron center, whereas for unsym/sym-[1(CO)](+), the unpaired spin was delocalized over both iron centers, as indicated by the magnitude of the hyperfine coupling to the phosphine ligands trans to the Fe-Fe vector. Oxidation of 1 by 2 equiv of acetylferrocenium afforded the dication [1](2+), which, on the basis of low-temperature IR spectrum, is structurally similar to [1](+). Treatment of [1](2+) with CO gives unsym-[1(CO)](2+).

Terminal hydride in [FeFe]-hydrogenase model has lower potential for H2 production than the isomeric bridging hydride

Protonation of the symmetrical tetraphosphine complexes Fe2(S2CnH2n)(CO)2(dppv)2 afforded the corresponding terminal hydrides, establishing that even symmetrical diiron(I) dithiolates undergo protonation at terminal sites. The terminal hydride [HFe2(S2C3H6)(CO)2(dppv)2](+) was found to catalyze proton reduction at potentials 200 mV milder than the isomeric bridging hydride, thereby establishing a thermodynamic advantage for catalysis operating via terminal hydride. The azadithiolate protonates to afford, [Fe2[(SCH2)2NH2](CO)2(dppv)2](+), [HFe2[(SCH2)2NH](CO)2(dppv)2](+), and [HFe2[(SCH2)2NH2](CO)2(dppv)2](2+), depending on conditions.