Computational investigation of [FeFe]-hydrogenase models: characterization of singly and doubly protonated intermediates and mechanistic insights

Hybrids of [FeFe]- and [NiFe]-H2ase Active Site Models

Complexes of the type (diphosphine)Ni(μ-SR)2Fe(CO)3 are investigated with azadithiolate (adt, HN(CH2S-)2) as the dithiolate. The resulting complexes are hybrid models for the active sites of the [NiFe]- and [FeFe]-hydrogenases. The key complex (dppv)Ni(μ-adt)Fe(CO)3 (3) was prepared from the complex Ni[(SCH2)2NCbz](dppv), which contains a Cbz-protected adt ligand (Cbz = C(O)OCH2Ph, dppv = cis-1,2-(Ph2P)2C2H2). This complex combines with Fe2(CO)9 to give (dppv)Ni[(μ-SCH2)2NCbz]Fe(CO)3, which is readily deprotected to give 3. Complex 3 undergoes protonation at both Fe and N to give successively [(dppv)Ni(μ-adt)FeH(CO)3]+ ([H3]+) and [(dppv)Ni(μ-adtH)FeH(CO)3]2+ ([H3H]2+). The redox properties and dynamics of these complexes resemble previously reported analogues with propanedithiolate. Solutions of [H3]+ readily degrade to [(dppv)Ni[(μ-SCH2)2NCH2]Fe(CO)3]+ ([4]+), which features a methylene group linking N and Fe. Complex [4]+ can be made in high yield by reaction of [H3]+ with CH2O, and this conversion was also demonstrated with 13CH2O. Complex [4]+ undergoes hydrogenolysis by photochemical reaction with H2 to give [(dppv)Ni[(μ-SCH2)2NMe]FeH(CO)3]+, the N-methylated analogue of [H3]+. Upon treatment ith Me3O+, [4]+ undergoes quaternization, giving [(dppv)Ni[(μ-SCH2)2N(Me)CH2]Fe(CO)3]2+. In contrast with the lability of [H3]+, the phosphine-substituted derivative [(dppv)Ni(μ-adt)FeH(CO)2(PPh3)]+ did not degrade. Most complexes were characterized by X-ray crystallography.

Catalytic Mechanism of Competing Proton Transfer Events from Water and Acetic Acid by [CoII(bpbH2)Cl2] for Water Splitting Processes

We performed first principles simulations to explore the water reduction process of the cobalt complex [Co^II(bpbH₂)Cl₂], where bpbH₂ = N,N'-bis(2'-pyridine carboxamide)-1,2-benzene. We considered the sequence steps of electron reduction followed by the proton addition process to observe the hydrogen evolution process. An experimental study of the catalyst showed that the increase in the acetic acid concentration triggers catalytic current and reduction of Co(II) to Co(I), and protonation occurred, yielding a Co(III)-H intermediate. Therefore, we used water and acetic acid as the proton sources. We compare the proton transfer kinetics from both the water and acetic acid. The reduction potentials and proton transfer kinetics from water or acetic acid to the reaction center were studied in a DMF solvent through the implicit solvent model. The first proton transfer from the acetic acid is more favorable, forming a Co^III-H complex and further reducing to Co^II-H. The second proton transfer from water to the Co^II-H moiety requires less free energy than acetic acid and is the rate-limiting step. The nature of the reduction process is also examined through the charge analysis, which reveals that the ligand becomes softer due to the C═O groups.

Linear Scaling Relationships to Predict pKa's and Reduction Potentials for Bioinspired Hydrogenase Catalysis

Biomimetic catalysts inspired by the active site of the [FeFe] hydrogenase enzyme can convert protons into molecular hydrogen. Minimizing the overpotential of the electrocatalytic process remains a major challenge for practical application of the catalyst. The catalytic cycle of the hydrogen production follows an ECEC mechanism (E represents an electron transfer step, and C refers to a chemical step), in which the electron and proton transfer steps can be either sequential or coupled (PCET). In this study, we have calculated the pKa's and the reduction potentials for a series of commonly used ligands (80 different complexes) using density functional theory. We establish that the required acid strength for protonation at the Fe-Fe site correlates with the standard reduction potential of the di-iron complexes with a linear energy relationship. These linear relationships allow for fast screening of ligands and tuning of the properties of the catalyst. Our study also suggests that bridgehead ligand properties, such as bulkiness and aromaticity, can be exploited to alter or even break the linear scaling relationships.

Electrocatalytic Behavior of Tetrathiafulvalene (TTF) and Extended Tetrathiafulvalene (exTTF) [FeFe] Hydrogenase Mimics

TTF- and exTTF-containing [(μ-S₂)Fe₂(CO)₆] complexes have been prepared by the photochemical reaction of TTF or exTTF and [(μ-S₂)Fe₂(CO)₆]. These complexes are able to interact with PAHs. In the absence of air and in acid media an electrocatalytic dihydrogen evolution reaction (HER) occurs, similarly to analogous [(μ-S₂)Fe₂(CO)₆] complexes. However, in the presence of air, the TTF and exTTF organic moieties strongly influence the electrochemistry of these systems. The reported data may be valuable in the design of [FeFe] hydrogenase mimics able to combine the HER properties of the [FeFe] cores with the unique TTF properties.

Successes, challenges, and opportunities for quantum chemistry in understanding metalloenzymes for solar fuels research

Overview of the rich and diverse contributions of quantum chemistry to understanding the structure and function of the biological archetypes for solar fuel research, photosystem II and hydrogenases.

Molecular Electrocatalysts for the Hydrogen Evolution Reaction: Input from Quantum Chemistry

In the pursuit of carbon-free fuels, hydrogen can be considered as an apt energy carrier. The design of molecular electrocatalysts for hydrogen production is important for the development of renewable energy sources that are abundant, inexpensive, and environmentally benign. Over the last 20 years, a large number of electrocatalysts have been developed, and considerable efforts have been directed toward the design of earth-abundant, first-row transition-metal complexes capable of promoting electrocatalytic hydrogen evolution reaction (HER). In this context, numerical approaches have emerged as powerful tools to study the catalytic performances of these complexes. This review covers some of the most significant theoretical mechanistic studies of biomimetic and bioinspired homogeneous HER catalysts. The approaches employed to study the free energy landscapes are discussed and methods used to obtain accurate estimates of relevant observables required to study the HER are presented. Furthermore, the structural and electronic parameters that govern the reactivity, and are necessary to achieve efficient hydrogen production, are discussed in view of future research directions.

Can carbene decorated [FeFe]-hydrogenase model complexes catalytically produce dihydrogen? An insight from theory

Cyclic alkyl amino carbene (CAAC) anchored [FeFe]-hydrogenase model complex featuring rotated conformation at one of the iron centers are found to be promising candidate for effective production of dihydrogen. A stepwise comparison of the complete mechanism using the CAAC stabilized model complex [1]⁰ has been performed with that of an experimentally isolated one ([2]⁰). Interestingly, the reduction events involved in the catalytic cycles are found to be more favorable than those previously reported for a similar experimentally known system. Furthermore, the computed ΔpK_a values indicate that the distal iron center with a vacant coordination site is more basic compared to the amino nitrogen atom of the azadithiolate bridge. We also made an attempt to determine the oxidation states of the iron centers for the intermediates involved in the catalytic cycles on the basis of the computed Mössbauer isomer shift and Mulliken spin density values.

Performance of enhanced DuBois type water reduction catalysts (WRC) in artificial photosynthesis - effects of various proton relays during catalysis

Inspired by natural photosynthesis, features such as proton relays have been integrated into water reduction catalysts (WRC) for effective production of hydrogen. Research by DuBois et al. showed the crucial influence of these relays, largely in the form of pendant amine functions. In this work catalysts are presented containing innovative diphosphinoamine ligands: [M(ii)Cl2(PNP-C1)], [M(ii)(MeCN)2(PNP-C1)]2+, [M(ii)(PNP-C1)2]2+, and [M(ii)Cl(PNP-C2)]+ (M = Pt2+, Pd2+, Ni2+, Co2+; PNP-C1 = N,N-bis{(di(2-methoxyphenyl)phosphino)methyl}-N-alkylamine, PNP-C2 = N,N-bis{(di(2-methoxyphenyl)phosphino)ethyl}-N-alkylamine and alkyl = Me, Et, iso-Pr, Bz). Synthetic strategies and detailed characterisation are covered, including 1H-, 13C-, and 31P-NMR analysis, mass spectroscopy and single crystal X-ray diffractometry (XRD). The catalytic properties have been explored by changing the pendant amines and auxiliary methoxy coordination sites, as well as enlarging the ligand backbone. Moreover, confirmed by density functional theory (DFT) calculations based on XRD data in vacuo and solvent environment, two very different catalytic cycles are proposed. PNP-C1 shows a classical proton relay, whereas PNP-C2 allows an additional coordination of nitrogen, acting optionally like a pincer. Through new insights into efficiency and stability-increasing influences of proton relays in general, their number per metal centre, an enlarged ligand backbone and the use of solvato instead of halogenido complexes, substantial improvements have been made in catalytic performance over the DuBois et al. catalysts and recently self-made WRCs. The turnover number (TON) related to the single site of cost-efficient nickel WRCs is increased from 11.4 to 637, whereas a corresponding palladium catalyst gives TON as high as 2289.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Models of the iron-only hydrogenase enzyme: structure, electrochemistry and catalytic activity of Fe2(CO)3(μ-dithiolate)(μ,κ1,κ2-triphos)

A series of Fe₂(triphos)(CO)₃(μ-dithiolate) complexes have been prepared and studied as models of the diiron centre in [FeFe]-hydrogenases.

Crystallographic and spectroscopic assignment of the proton transfer pathway in [FeFe]-hydrogenases

The unmatched catalytic turnover rates of [FeFe]-hydrogenases require an exceptionally efficient proton-transfer (PT) pathway to shuttle protons as substrates or products between bulk water and catalytic center. For clostridial [FeFe]-hydrogenase CpI such a pathway has been proposed and analyzed, but mainly on a theoretical basis. Here, eleven enzyme variants of two different [FeFe]-hydrogenases (CpI and HydA1) with substitutions in the presumptive PT-pathway are examined kinetically, spectroscopically, and crystallographically to provide solid experimental proof for its role in hydrogen-turnover. Targeting key residues of the PT-pathway by site directed mutagenesis significantly alters the pH-activity profile of these variants and in presence of H₂ their cofactor is trapped in an intermediate state indicative of precluded proton-transfer. Furthermore, crystal structures coherently explain the individual levels of residual activity, demonstrating e.g. how trapped H₂O molecules rescue the interrupted PT-pathway. These features provide conclusive evidence that the targeted positions are indeed vital for catalytic proton-transfer.

[FeFe]-hydrogenases catalyze H₂-evolution and -oxidation at very high turnover-rates. Here the authors provide experimental evidence for the proposed proton-transfer (PT) pathway by kinetically, spectroscopically, and crystallographically characterizing eleven mutants from the two [FeFe]-hydrogenases CpI and HydA1.

Reaction Coordinate Leading to H2 Production in [FeFe]-Hydrogenase Identified by Nuclear Resonance Vibrational Spectroscopy and Density Functional Theory

[FeFe]-hydrogenases are metalloenzymes that reversibly reduce protons to molecular hydrogen at exceptionally high rates. We have characterized the catalytically competent hydride state (Hhyd) in the [FeFe]-hydrogenases from both Chlamydomonas reinhardtii and Desulfovibrio desulfuricans using 57Fe nuclear resonance vibrational spectroscopy (NRVS) and density functional theory (DFT). H/D exchange identified two Fe-H bending modes originating from the binuclear iron cofactor. DFT calculations show that these spectral features result from an iron-bound terminal hydride, and the Fe-H vibrational frequencies being highly dependent on interactions between the amine base of the catalytic cofactor with both hydride and the conserved cysteine terminating the proton transfer chain to the active site. The results indicate that Hhyd is the catalytic state one step prior to H2 formation. The observed vibrational spectrum, therefore, provides mechanistic insight into the reaction coordinate for H2 bond formation by [FeFe]-hydrogenases.

Interplay between Terminal and Bridging Diiron Hydrides in Neutral and Oxidized States

This study describes the structural, spectroscopic, and electrochemical properties of electronically unsymmetrical diiron hydrides. The terminal hydride Cp*Fe(pdt)Fe(dppe)(CO)H ([1(t-H)]⁰, Cp*^- = Me₅C₅^-, pdt^2- = CH₂(CH₂S^-)₂, dppe = Ph₂PC₂H₄PPh₂) was prepared by hydride reduction of [Cp*Fe(pdt)Fe(dppe)(CO)(NCMe)]⁺. As established by X-ray crystallography, [1(t-H)]⁰ features a terminal hydride ligand. Unlike previous examples of terminal diiron hydrides, [1(t-H)]⁰ does not isomerize to the bridging hydride [1(μ-H)]⁰. Oxidation of [1(t-H)]⁰ gives [1(t-H)]⁺, which was also characterized crystallographically as its BF₄^- salt. Density functional theory (DFT) calculations indicate that [1(t-H)]⁺ is best described as containing an Cp*Fe^III center. In solution, [1(t-H)]⁺ isomerizes to [1(μ-H)]⁺, as anticipated by DFT. Reduction of [1(μ-H)]⁺ by Cp₂Co afforded the diferrous bridging hydride [1(μ-H)]⁰. Electrochemical measurements and DFT calculations indicate that the couples [1(t-H)]^+/0 and [1(μ-H)]^+/0 differ by 210 mV. Qualitative measurements indicate that [1(t-H)]⁰ and [1(μ-H)]⁰ are close in free energy. Protonation of [1(t-H)]⁰ in MeCN solution affords H₂ even with weak acids via hydride transfer. In contrast, protonation of [1(μ-H)]⁰ yields 0.5 equiv of H₂ by a proposed protonation-induced electron transfer process. Isotopic labeling indicates that μ-H/D ligands are inert.

Direct Observation of an Iron-Bound Terminal Hydride in [FeFe]-Hydrogenase by Nuclear Resonance Vibrational Spectroscopy

[FeFe]-hydrogenases catalyze the reversible reduction of protons to molecular hydrogen with extremely high efficiency. The active site ("H-cluster") consists of a [4Fe-4S]H cluster linked through a bridging cysteine to a [2Fe]H subsite coordinated by CN- and CO ligands featuring a dithiol-amine moiety that serves as proton shuttle between the protein proton channel and the catalytic distal iron site (Fed). Although there is broad consensus that an iron-bound terminal hydride species must occur in the catalytic mechanism, such a species has never been directly observed experimentally. Here, we present FTIR and nuclear resonance vibrational spectroscopy (NRVS) experiments in conjunction with density functional theory (DFT) calculations on an [FeFe]-hydrogenase variant lacking the amine proton shuttle which is stabilizing a putative hydride state. The NRVS spectra unequivocally show the bending modes of the terminal Fe-H species fully consistent with widely accepted models of the catalytic cycle.

Proton Coupled Electronic Rearrangement within the H-Cluster as an Essential Step in the Catalytic Cycle of [FeFe] Hydrogenases

The active site of [FeFe] hydrogenases, the H-cluster, consists of a [4Fe-4S] cluster connected via a bridging cysteine to a [2Fe] complex carrying CO and CN^- ligands as well as a bridging aza-dithiolate ligand (ADT) of which the amine moiety serves as a proton shuttle between the protein and the H-cluster. During the catalytic cycle, the two subclusters change oxidation states: [4Fe-4S]_H²⁺ ⇔ [4Fe-4S]_H⁺ and [Fe(I)Fe(II)]_H ⇔ [Fe(I)Fe(I)]_H thereby enabling the storage of the two electrons needed for the catalyzed reaction 2H⁺ + 2e^- ⇄ H₂. Using FTIR spectro-electrochemistry on the [FeFe] hydrogenase from Chlamydomonas reinhardtii (CrHydA1) at different pH values, we resolve the redox and protonation events in the catalytic cycle and determine their intrinsic thermodynamic parameters. We show that the singly reduced state H_red of the H-cluster actually consists of two species: H_red = [4Fe-4S]_H⁺ - [Fe(I)Fe(II)]_H and H_redH⁺ = [4Fe-4S]_H²⁺ - [Fe(I)Fe(I)]_H (H⁺) related by proton coupled electronic rearrangement. The two redox events in the catalytic cycle occur on the [4Fe-4S]_H subcluster at similar midpoint-potentials (-375 vs -418 mV); the protonation event (H_red/H_redH⁺) has a pK_a ≈ 7.2.

Hydrogenase Enzymes and Their Synthetic Models: The Role of Metal Hydrides

Hydrogenase enzymes efficiently process H2 and protons at organometallic FeFe, NiFe, or Fe active sites. Synthetic modeling of the many H2ase states has provided insight into H2ase structure and mechanism, as well as afforded catalysts for the H2 energy vector. Particularly important are hydride-bearing states, with synthetic hydride analogues now known for each hydrogenase class. These hydrides are typically prepared by protonation of low-valent cores. Examples of FeFe and NiFe hydrides derived from H2 have also been prepared. Such chemistry is more developed than mimicry of the redox-inactive monoFe enzyme, although functional models of the latter are now emerging. Advances in physical and theoretical characterization of H2ase enzymes and synthetic models have proven key to the study of hydrides in particular, and will guide modeling efforts toward more robust and active species optimized for practical applications.

Brønsted-Lowry Acid Strength of Metal Hydride and Dihydrogen Complexes

Transition metal hydride complexes are usually amphoteric, not only acting as hydride donors, but also as Brønsted-Lowry acids. A simple additive ligand acidity constant equation (LAC for short) allows the estimation of the acid dissociation constant Ka(LAC) of diamagnetic transition metal hydride and dihydrogen complexes. It is remarkably successful in systematizing diverse reports of over 450 reactions of acids with metal complexes and bases with metal hydrides and dihydrogen complexes, including catalytic cycles where these reactions are proposed or observed. There are links between pKa(LAC) and pKa(THF), pKa(DCM), pKa(MeCN) for neutral and cationic acids. For the groups from chromium to nickel, tables are provided that order the acidity of metal hydride and dihydrogen complexes from most acidic (pKa(LAC) -18) to least acidic (pKa(LAC) 50). Figures are constructed showing metal acids above the solvent pKa scales and organic acids below to summarize a large amount of information. Acid-base features are analyzed for catalysts from chromium to gold for ionic hydrogenations, bifunctional catalysts for hydrogen oxidation and evolution electrocatalysis, H/D exchange, olefin hydrogenation and isomerization, hydrogenation of ketones, aldehydes, imines, and carbon dioxide, hydrogenases and their model complexes, and palladium catalysts with hydride intermediates.

DFT study of the mechanism of hydrogen evolution catalysed by molecular Ni, Co and Fe catalysts containing a diamine–tripyridine ligand

Electrolysis of water to obtain hydrogen is a practical way to transform surplus electrical power into clean and sustainable hydrogen fuels.

Diiron azadithiolates as models for the [FeFe]-hydrogenase active site and paradigm for the role of the second coordination sphere

The [FeFe] hydrogenases (H2ases) catalyze the redox reaction that interconverts protons and H2. This area of biocatalysis has attracted attention because the metal-based chemistry is unusual, and the reactions have practical implications. The active site consists of a [4Fe-4S] cluster bridged to a [Fe2(μ-dithiolate)(CN)2(CO)3](z) center (z = 1- and 2-). The dithiolate cofactor is [HN(CH2S)2](2-), called the azadithiolate ([adt(H)](2-)). Although many derivatives of Fe2(SR)2(CO)6-xLx are electrocatalysts for the hydrogen evolution reaction (HER), most operate by slow nonbiomimetic pathways. Biomimetic hydrogenogenesis is thought to involve intermediates, wherein the hydride substrate is adjacent to the amine of the adt(H), being bonded to only one Fe center. Formation of terminal hydride complexes is favored when the diiron carbonyl models contain azadithiolate. Although unstable in the free state, the adt cofactor is stable once it is affixed to the Fe2 center. It can be prepared by alkylation of Fe2(SH)2(CO)6 with formaldehyde in the presence of ammonia (to give adt(H) derivatives) or amines (to give adt(R) derivatives). Weak acids protonate Fe2(adt(R))(CO)2(PR3)4 to give terminal hydrido (term-H) complexes. In contrast, protonation of the related 1,3-propanedithiolate (pdt(2-)) complexes Fe2(pdt)(CO)2(PR3)4 requires strong acids. The amine in the azadithiolate is a kinetically fast base, relaying protons to and from the iron, which is a kinetically slow base. The crystal structure of the doubly protonated model [(term-H)Fe2(Hadt(H))(CO)2(dppv)2](2+) confirms the presence of both ammonium and terminal hydrido centers, which interact through a dihydrogen bond (dppv = cis-C2H2(PPh2)2). DFT calculations indicate that this H---H interaction is sensitive to the counterions and is strengthened upon reduction of the diiron center. For the monoprotonated models, the hydride [(term-H)Fe2(adt(H))(CO)2(dppv)2](+) exists in equilibrium with the ammonium tautomer [Fe2(Hadt(H))(CO)2(dppv)2](+). Both [(term-H)Fe2(Hadt(H))(CO)2(dppv)2](2+) and [(term-H)Fe2(adt(H))(CO)2(dppv)2](+) are highly active electrocatalysts for HER. Catalysis is initiated by reduction of the diferrous center, which induces coupling of the protic ammonium center and the hydride ligand. In contrast, the propanedithiolate [(term-H)Fe2(pdt)(CO)2(dppv)2](+) is a poor electrocatalyst for HER. Oxidation of H2 has been demonstrated, starting with models for the oxidized state ("Hox"), for example, [Fe2(adt(H))(CO)3(dppv)(PMe3)](+). Featuring a distorted Fe(II)Fe(I) center, this Hox model reacts slowly with high pressures of H2 to give [(μ-H)Fe2(adt(H))(CO)3(dppv)(PMe3)](+). Highlighting the role of the proton relay, the propanedithiolate [Fe2(pdt)(CO)3(dppv)(PMe3)](+) is unreactive toward H2. The Hox-model + H2 reaction is accelerated in the presence of ferrocenium salts, which simulate the role of the attached [4Fe-4S] cluster. The redox-complemented complex [Fe2(adt(Bn))(CO)3(dppv)(FcP*)](n+) catalyzes both proton reduction and hydrogen oxidation (FcP* = (C5Me5)Fe(C5Me4CH2PEt2)).

Application of the energetic span model to the electrochemical catalysis of proton reduction by a diiron azadithiolate complex

A method for the computation of TOF of catalysis of electrochemical reaction as a function of the potential was developed.