Frédéric Gloaguen*a and Thomas B. Rauchfuss*b
aUMR CNRS 6521, Université de Bretagne Occidentale, CS 93837, 29238 Brest cedex 3, France. E-mail: fgloaguen@gmail.com
bSchool of Chemical Sciences, University of Illinois, Urbana, IL 61801, USA. E-mail: rauchfuz@uiuc.edu
First published on 31st October 2008
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.
![]() Frédéric Gloaguen | Frédéric Gloaguen received his undergraduate degree in Brest in 1990 and his PhD in Grenoble in 1994. He worked at the University of Poitiers and was hired by CNRS in 1996. After a sabbatical leave in 2001 at the University of Illinois, he joined the University of Brest in 2002. He is interested in all aspects of the catalysis of electrochemical reactions involved in energy storage and utilization. |
![]() Thomas Rauchfuss | Thomas Rauchfuss received his PhD degree at Washington State University in 1975. After a postdoc in Australian National University, he joined the faculty at the University of Illinois where he conducts research in synthetic inorganic and organometallic chemistry. In recent years he has focused on models for the active sites of the hydrogenases. In collaboration with Gloaguen, his group discovered the catalytic hydrogen evolution by diiron dithiolates. |
In an effort to understand the molecular mechanisms by which H2ases operate, much research has been aimed at mimicking the structures of their active sites. In recent years, these models have begun to yield biochemically significant insights, although gaps remain. Perhaps most perplexing are the high rates achieved by these enzymes, especially in view of the fact that they utilize first row metals that typically display diminished affinities for dihydrogen.2 Furthermore and still more challenging, the H2ases effect their reactions via apparent 1e− changes, which require odd-electron intermediates. The one-electron chemistry of metal hydrides and metal–dihydrogen complexes is lightly studied; thus, the biochemical mechanisms present opportunities for learning new organometallic chemistry relevant to dihydrogen.
The literature on the production of hydrogen in solution, homogeneous hydrogenogenesis, is not extensive,3,4 but the coordination chemistry of dihydrogen has been active for decades2 and is obviously relevant to biological processes. Although H2 itself exhibits neither redox nor any acid–base reactivity, its metal complexes exhibit both, i.e. dihydrogen complexes can be highly acidic and, the derived metal hydrides can be oxidized. It therefore makes sense that the active sites of both the [NiFe]- and the [FeFe]-H2ases feature metals. This generalization extends to the recently discovered hydrogen-transfer enzyme Hmd.5 Interestingly, the active sites of all three classes of “hydrogen-processing” enzymes feature thiolato iron carbonyl entities.
This review brings a particular focus on the reactivity of models for the H2ases. Reactivities of interest include protonation, binding of H2 and CO (a common inhibitor of H2ases), and redox.
![]() | ||
| Fig. 1 Structure of H-cluster from C. pasteurianum, a line drawing of the same, omitting the water situated near the distal (rightmost) Fe center, and the structure of the prototypical Fe2(pdt)(CO)6. | ||
Redox is a key attribute of the H2ases, and thus early studies attempted, fruitlessly, to oxidize the dicyanides to mixed-valence derivatives. Useful insights have however been obtained from the low-temperature oxidation of [Fe2[(SCH2)2C(Me)CH2SR](CN)2(CO)4]2−, wherein the dithiolato ligand bears a pendant, uncoordinated thioether group. At −40 °C, this species undergoes 1e− oxidation to give a mixed-valence intermediate, which is proposed to be stabilized by coordination of the thioether. Characterization by EPR and IR spectroscopies confirm the similarity of this mixed valence species to the CO-inhibited (HCOox) state of the enzyme.8
Oxidation of the diiron dicyanides typically affords intractable solids that appear to be polymers containing μ-CN linkages. Some details of the redox-induced aggregation of the diiron cyanides were elucidated in the case of [Fe2(edt)(CN)(CO)3(dppv)]−. This species undergoes complete consumption upon treatment with only one equivalent of Fc+ to give a diamagnetic Fe4 derivative [FeI2(S2C2H4)(CO)3(dppv)](μ-CN)[FeII2(S2C2H4)(μ-CO)(CN)(CO)2(dppv)]. The result is consistent with the disproportionation of a mixed-valence diiron intermediate into a diamagnetic mixed valence tetrairon species, i.e., 2FeIFeII→ FeIFeIFeIIFeII (eqn (1)).9
![]() | (1) |
Particularly significant advances in biomimetic modeling have resulted from the study of diiron compounds where the cyanides are replaced all or in part by tertiary phosphine ligands, as well as carbenes and isocyanides. Phosphine complexes of the type Fe2(SR)2(CO)6−x(PR3)x, especially for x≥ 2, protonate readily to give characterizable hydrides. Whereas the isomeric hydrides of [Fe2(SR)2(μ-H)(CN)2(CO)4]− are observable by NMR spectroscopy, hydrides containing monodentate phosphine ligands typically are isolated as single isomers.12 Mixed phosphine-cyanide complexes represent a compromise between the uncomplicated properties of phosphines and the biologically relevant, highly basic cyanide. In contrast, the parent hexacarbonyls protonate only with the strongest acids, e.g. triflic acid in CH2Cl2 solution.
A representative and well-studied “μ-hydride” is [Fe2(edt)(μ-H)(CO)4(PMe3)2]+, wherein CO ligands are trans to the hydride and the phosphines occupy trans basal sites. Such μ-hydrides are invoked as intermediates in some pathways for electrocatalytic hydrogenogenesis (see section 4). Protonation strongly affects the diiron center: in contrast to its unprotonated precursor, the hydrides are susceptible to substitution reactions involving replacement of CO ligands with donor ligands, such as cyanide and phosphines. Further indicating the lability of the CO ligands, upon illumination the hydrido complexes also catalyze isotopic exchange between D2 and H2O. H2ases characteristically catalyze such exchange, although photoactivation is not required.13
Since the time of the original crystallographic characterization, the enzymatic reactions have been assumed to proceed via the intermediacy of a hydride located on the apical site of the distal iron center. A model for such a terminal hydrido species has been characterized crystallographically in the form of [HFe2(edt)(CO)2(PMe3)4]BF4. The hydride is located at one of the two apical sites and one CO group semibridges the two Fe centers. Unlike the isomeric μ-hydrido complex [Fe2(edt)(μ-H)(CO)2(PMe3)4]BF4, the terminal hydride reacts with strong acids to release H2 (eqn (2)).14
![]() | (2) |
In situ analyses show that diiron dithiolato carbonyls initially protonate a single iron center.15,16 In most cases, the resulting terminal hydride complexes rapidly rearrange to give the bridged hydride isomer. Unsymmetrical complexes such as Fe2(pdt)(CO)4(dppe) protonate mainly at the Fe(CO)3 site, not the ostensibly more basic Fe(CO)(dppe) site. The terminal hydrides are more stable when the diiron dithiolato center is both electron-rich and has bulky ligands that hinder turnstile rotation. Such species can even be observed near room temperature and interrogated spectroscopically in detail (Scheme 1). Thus, protonation of Fe2(pdt)(CO)2(dppv)2 with HBF4 gives the terminal hydride [HFe2(pdt)(CO)2(dppv)2]+, which unimolecularly isomerizes to a series of isomeric μ-hydrides. The terminal hydride reduces at potentials ca. 200 mV more positive than the isomeric bridging hydride [Fe2(pdt)(μ-H)(CO)2(dppv)2]+. This lowered reduction potential points to a thermodynamic advantage for terminal hydrides as precursors to H2.17 Reduction of this diferrous hydride generates a mixed-valence hydride that is highly reactive toward protons to release H2. Overall, these studies indicate several important aspects: terminal hydride ligands form readily by protonation, they more easily undergo protonolysis to H2, and they reduce more easily than the isomeric bridging hydrides.
N-Protonation is favored thermodynamically for Fe2(adt)(CO)6 (with pKa’s ≈ 8 in MeCN solution). Replacement of the carbonyl ligands with alkyl phosphines (and presumably cyanides) enhances the basicity of the iron centers. Even in Fe2(adt)(CO)4(PMe3)2, where the FeFe bond is 1000× more basic than the amine, N-protonation is kinetically favored.21 In cases where the terminal hydride is more stabilized, equilibration between the hydride and ammonium derivative can be observed (Scheme 1). The emerging mechanistic picture is that N-protonation precedes and facilitates the formation and deprotonation of the terminal hydride.
![]() | ||
| Scheme 1 N- and Fe-protonation of a diiron azadithiolato complex and its 1e−-reduction. | ||
Even propanedithiolate noticeably affects the coordinating tendency of the Fe centers since the central methylene group projects over this apical site. This steric shielding is evidenced by the effect of the pdt vs. edt on the stereochemistry of Fe2(edt/pdt)(CO)4(diphosphine). For edt, the diphosphine exclusively spans apical and basal sites, whereas in the pdt derivative, one observes significant amounts of the dibasal isomer. The rates of carbonylation of the unsaturated Hox models containing pdt are also slower for the edt derivatives (see below).
A model for the entire H-cluster has been prepared with a 4Fe-4S subunit linked via thiolato bridges to one or more Fe2(SR)2(CO)5 centers (Fig. 2). In MeCN solution, the synthetic 6Fe ensemble is an electrocatalyst for hydrogen evolution using 3,5-dimethylpyridinium as a proton source. In addition to representing a significant synthetic achievement, this 6Fe cluster provides insights into the nature of the linkage between the 2Fe and the 4Fe subunits, insights potentially relevant to engineering biomimetic catalysts.22 First, attachment of the Fe2(SR)2(CO)5 fragment to one thiolate of Fe4S4(SR)42− shifts the reduction potential of the 4Fe cluster in the positive direction by 150 mV. It remains unclear how the redox properties of the 4Fe-4S cluster will be affected by more realistic diiron fragments such as [Fe2(SR)2(CN)2(CO)3]x−. Second, IR measurements in the νCO region suggest that the 4Fe-4S cluster-ligand is thioether-like in its Lewis basicity. This μ-SR linkage between the [4Fe-4S] and the Fe2(SR)2(CO)5 subunits breaks easily upon reduction of the dianionic 6Fe, whereas in the enzyme this bond is enforced by the protein.
Fe2(SR)2(CO)6 centers have also been modified with abiological redox auxiliaries.23,24 These assemblies are designed to couple the hydrogen producing properties of the diiron dithiolato carbonyls with sources of photogenerated reducing equivalents. To this end, diiron dithiolato carbonyls have been covalently linked to light-harvesting chromophores, such as [Ru(bipy)3]2+ and Zn(TPP), which serve as antennae. Excitation of the sensitizer produces a reducing equivalent that, when transferred to the diiron assembly, could potentially reduce protons to H2. Initial designs have focused on matching the properties of the photoreductant—lifetimes, quantum yields, and reducing properties—with the redox properties of the diiron center. A long-range goal of such work is the development of a catalyst for water splitting.
Highlighting the versatility and robustness of the Fe2(SR)2(CO)6 species, photosensitizers have been attached via several means, including linking to the azadithiolate and through phosphine ligands (Fig. 2). The diiron(I) dithiolato framework is compatible with numerous conjugation and assembly methods which will enable matching of the reducing potential of the photogenerated reducing equivalents with that of the underlying diiron center (ca.−1.2 V vs. SCE).23
![]() | ||
| Scheme 2 Generation of models for the Hox and HoxCO states. | ||
Modeling studies have yielded related coordinatively unsaturated diiron species exhibiting the expected magnetism, structure, and tendency to carbonylate.26 The mixed phosphine-carbene complex Fe2(pdt)(CO)4(PMe3)(IMes) undergoes oxidation at very mild potentials (−450 mV vs. Fc0/+) to give [Fe2(pdt)(CO)4(PMe3)(IMes)]+, isolated as its PF6− salt. The tetra- and trisubstituted derivatives of the type Fe2(edt/pdt)(CO)2(PR3)4 and Fe2(edt/pdt)(CO)3(PR3)3 can also be oxidized by ferrocenium to yield the corresponding monocations. The cations are somewhat thermally sensitive, typically decomposing within minutes near room temperature, but they are amenable to EPR, IR, and crystallographic analyses that show a strong similarity to the active site in the Hox state.
For these mixed valence species, EPR studies indicate that the rotated Fe center is Fe(I), thus oxidation causes “rotation without redox” at distal Fe and “redox without rotation” at the proximal Fe. In [Fe2(pdt)(CO)4(PMe3)(IMes)]+, the rotated (Fe(I)) site carries the bulky IMes ligand. DFT calculations suggest that the rotated Fe(I) site is more electrophilic, despite its lower oxidation state. Diamagnetic analogues of Hox, prepared using NO in place of CO, indicate that rotation of the distal Fe is caused by the electronic asymmetry that induces the shift of one CO ligand to a semi-bridging site.27
The unsaturated [Fe2(SR)2(CO)6−x(PR3)x]+ species reversibly bind CO, as observed also for the Hox state of the enzyme.28 IR studies of the 13CO-labeled adduct, a model for the HoxCO state, indicate that the exogenous CO binds to the vacant apical site. EPR results, in particular the 31P hyperfine values, indicate that the resulting 35e− adducts are valence-delocalized, in contrast to the Hox models. At 2.5 Å, the Fe–Fe distance in the Hox models is virtually unchanged relative to the 34e− precursor complexes,29 but upon binding CO, the Fe–Fe distance elongates by 0.1 Å. Electrochemical studies show that these HoxCO models are susceptible to oxidation to diamagnetic diferrous species, which are akin to the diamagnetic Hoxair state. CO is virtually unique in forming stable adducts with the Hox models, stronger donors result in disproportionation (eqn (3)).
| 2[Fe2(SR)2(μ-CO)L6]+⇌ [Fe2(SR)2(μ-CO)L6]2+ + Fe2(SR)2L6 + CO | (3) |
![]() | ||
| Fig. 3 Structure of the active site of the [NiFe]-H2ase and line drawing of the same. | ||
Although the [NiFe]-H2ases are pervasive and the structure announced several years ago, functional models have lagged structural models.31 The barriers to functional modeling result from several factors: (i) the heterometallic nature of the [NiFe] site is inherently more challenging synthetically than homodimetallic species, (ii) the Ni center is structurally unusual, (iii) the ambidentate character of both the cyanide (on Fe) and terminal thiolato (on Ni) ligands complicates the assembly of discrete species, and (iv) easily accessed models are cationic, whereas it is likely that the active site is anionic or charge-neutral. Underscoring the synthetic challenge is the corresponding complexity of the biological assembly pathway.30
The [NiFe]-H2ases are characteristically redox active. Mononuclear iron(II) carbonyls resist oxidation, thus it is logical that redox is centered on Ni, as implicated in several biophysical studies. Nickel shuttles between S = 0 Ni(II) and the otherwise uncommon S = 1/2 Ni(III) near −250 mV (vs. NHE). Model studies confirm that thiolates stabilize Ni(III).32 Well defined reactivity of these Ni(III) thiolates toward H2 has not been demonstrated.
![]() | ||
| Scheme 3 Selected models for the [NiFe]-H2ases active site featuring the Ni(SR)2 core and at least one biomimetic ligand.5,29,34,35 | ||
Prior to the structural characterization of the [NiFe]-H2ases, hydride and carbonyl derivatives of ferrocyanides were essentially unknown. The analogy between the coordination properties of FeII(CN)5 and RuII(NH3)5, noted by Taube,40 suggested the possibility that the cyanoferrous centers would bind H2, as has been seen for the Ru(II) ammines.2 Although the iron hexacyanides have long been known, the systematic chemistry of mixed ligand [Fe(CN)6−x(CO)x](4−x)− complexes has only recently been developed and four such species are now known. These studies show that the presence of even one CO ligand strongly stabilizes the ferrous state, arguing against a redox role for the iron.41 The emerging picture is that the FeII(CN)2(CO) center is well suited to serve as a Lewis acidic site for stabilizing the bridging hydride.
The chemistry of iron hydrido and dihydrogen complexes is well developed, largely due to the pioneering studies of Morris et al. on the [Fe(H2)L(diphosphine)2]2+ systems.42 The hydrido carbonyl cyanide [HFe(CN)2(CO)3]− has been prepared as a mixture of two major isomers. This anion not only features several of the ligands observed in active site, it releases H2 upon protonation.43 The unanswered question is whether CO ligands in this complex could be further substituted to accommodate thiolate ligands as seen in the active site.
Proton reduction catalysis is generally indicated by changes in cyclic voltammetry upon addition of acid HA. The main criterion is that the height of the reduction peak of the homogeneous catalyst, or its protonated form, increases with [HA] (Fig. 4). The increase in peak current is due to a catalytic cycle that produces H2 after electron transfers to the protonated and oxidized forms of the catalysts giving more current than in the absence of acid. At high [HA]/[catalyst] ratios, the catalytic current reaches a plateau value indicating that the current is limited by the rate of the catalytic cycle, not by proton diffusion. With respect to their catalytic properties, [FeFe]-H2ase models can be classified into two categories, depending on the sequence of protonation and reduction.
![]() | ||
| Fig. 4 Cyclic voltammograms of a solution 2.2 mM Fe2(bdt)(CO)6 in MeCN–Bu4NPF6 upon addition of an increasing amount of toluenesulfonic acid (HOTs) up to 33 mM. Inset: dependence of the catalytic peak current vs. [HOTs] (conditions: glassy carbon electrode 0.071 cm2 in surface area; scan rate 0.1 V s−1). | ||
![]() | (4) |
![]() | ||
| Scheme 4 | ||
The primary reduction of the propanedithiolate Fe2(pdt)(CO)6 occurs at E1/2,red = −1.63 V, 360 mV more negative than the bdt derivatives, indicating the superior donor ability of the alkyldithiolate.48 The chemical reversibility of the [Fe2(pdt)(CO)6]0/− couple depends on the solvent and is improved in the presence of CO. The charge passed upon electrolysis under CO atmosphere is consistent with an overall two-electron reduction. However the one-electron reduced species [Fe2(pdt)(CO)6]−, which is proposed to be structurally similar to its neutral precursor, is involved in several subsequent chemical reactions to give eventually an Fe4 product that is also catalytically active. Voltammetry of Fe2(pdt)(CO)6 in the presence of HOTs shows two different proton reduction waves. The first one (process I, Scheme 5) is associated with the reduction of the pdt complex. The second one (process II), which occurs at a potential more negative, appears in the presence of excess acid. Simulation of the cyclic voltammograms suggests that process I corresponds to H2 production through an ECEC mechanism with slow liberation of H2 from the doubly reduced intermediate, nominally H2Fe2(pdt)(CO)6. As a result, an additional electron transfer occurs (process II) leading to a 3-electron-2-proton intermediate that quickly releases H2.
![]() | ||
| Scheme 5 | ||
The preceding examples show that the nature of the dithiolate bridge strongly influences the reduction potential of the Fe2(SR)2(CO)6 complexes, and thus the potential at which H+/H2 catalysis occurs. In addition, at least two mechanisms must be considered depending on thermodynamic (pKa) and/or kinetic factors (H2 liberation rate). When the reduction process is not chemically reversible, as for the pdt derivative, products (dimers) formed in the absence of acid may also be involved in the electrocatalytic proton reduction process.49
Although monosubstituted derivatives Fe2(SR)2(CO)5L do not readily protonate, solutions of both Fe2(pdt)(CO)5(P(OMe)3) (Ep,red = −1.98 V) and Fe2(bdt)(CO)5(P(OMe)3) (Ep,red = −1.53 V) exhibit acid-responsive peaks at Ep,red = −1.65 and −1.18 V, respectively.16,52 A CE mechanism implicates the facile reduction of [(H)Fe2(SR)2(CO)5L]+. Simulations of voltammetric responses and electrolysis experiments suggest that hydrogen evolution catalysis may occur by a bimolecular reaction, i.e., 2HFe2(SR)2(CO)5L → 2Fe2(SR)2(CO)5L + H2, even in the presence of excess acid. Similarly, proton reduction by Fe2(bdt)(CO)5(P(OMe)3), which occurs at mild potentials, is also slow, due also to a bimolecular reaction.
, the potential of the proton/hydrogen couple and pKa,HA, the acid dissociation constant (eqn (5)).![]() | (5) |
and pKa,HA have been measured in a large number of solvents (Table 1).53
| Catalyst | Ecat/V | Acid | E0HA/Va | −η/V | Ref. |
|---|---|---|---|---|---|
a E0HA = − 0.059 pKa,HA, with = −0.14 and −0.77 V in MeCN and DMF, respectively. | |||||
| Fe2(bdt)(CO)6 | −1.31 | HOTs | −0.65 | 0.66 | 47 |
| Fe2(bdt)(CO)6 | −2.10 | HOAc | −1.46 | 0.64 | 46 |
| Fe2(pdt)(CO)6 (process I) | −1.59 | HOTs | −0.65 | 0.94 | 54 |
| Fe2(pdt)(CO)6 (process II) | −1.81 | HOTs | −0.65 | 1.16 | 48 |
| [Fe2(pdt)(μ-H)(CO)4(PMe3)(CNH)]+ | −1.43 | HOTs | −0.65 | 0.78 | 50 |
| Fe2(bdt)(CO)5P(OMe)3 | −1.18 | HOTs | −0.65 | 0.53 | 52 |
| Co(dmgBF2)2(MeCN)2 | −0.93 | CF3CO2H | −0.89 | 0.04 | 4 |
The possible role of the adt cofactor raises the possibility that an internal base in the diiron models would enhance the rate of H+/H2 catalysis, thereby lowering the overpotential. The rate-determining step for proton reduction catalysis by Fe2(adt)(CO)6 is, however, the slow H2 release from the 2-electron-2-proton intermediate,55 as seen for Fe2(pdt)(CO)6. Moreover, protonation of Fe2(adt)(CO)4(PMe3)2 gives the μ-hydride, which is too distant to be affected by the adt.12 Intramolecular proton–hydride coupling will be favored by minimizing the protonation-induced reorganization of the iron coordination sphere. A second condition for low overpotential is a matching of the pKa values of the H2-binding site with the protonated pendant base. These requirements have not been fulfilled in the models described to date. Progress is however imminent as protonation of Fe2(adt)(CO)2(dppv)2 gives terminal hydride [HFe2(adt)(CO)2(dppv)2]+ that is reactive towards protons upon reduction.17
By deviating from the Fe2(SR)2 and NiFe(SR)2 stoichiometries of the active sites, it is possible to generate new families of hydrogen-processing catalysts, those that are more bio-inspired than they are biomimetic. It appears that many combinations of nickel, iron carbonyls, and sulfide will activate H2 or reduce protons catalytically. Recent examples include catalysts with Fe4(SR)6,48 NiFe2(SR)2,33 and RuNi(SR)2 cores.35
Footnotes |
| † Part of the renewable energy theme issue. |
| ‡ Abbreviations: Fc+ = ferrocenium, pdt2− = 1,3-propanedithiolate, edt2− = 1,2-ethanedithiolate, adt2− = 2-aza-1,3-propanedithiolate, bdt2− = 1,2-benzenedithiolate, Imes = 1,3-bis(2,4,6-trimethylphenyl)imidazol-2-ylidene, dppv = 1,2-bis(diphenylphosphino)ethylene, dppe = 1,2-bis(diphenylphosphino)ethane, TPP2− = tetraphenylporphyrinate, bipy = 2,2′-bipyridine. |
| This journal is © The Royal Society of Chemistry 2009 |