Computational Studies of the H-Cluster of Fe-Only Hydrogenases: Geometric, Electronic, and Magnetic Properties and Their Dependence on the [Fe4S4] Cubane

Facile electrocatalytic proton reduction by a [Fe-Fe]-hydrogenase bio-inspired synthetic model bearing a terminal CN- ligand

An azadithiolate bridged CN- bound pentacarbonyl bis-iron complex, mimicking the active site of [Fe-Fe] H2ase is synthesized. The geometric and electronic structure of this complex is elucidated using a combination of EXAFS analysis, infrared and Mössbauer spectroscopy and DFT calculations. The electrochemical investigations show that complex 1 effectively reduces H+ to H2 between pH 0-3 at diffusion-controlled rates (1011 M-1 s-1) i.e. 108 s-1 at pH 3 with an overpotential of 140 mV. Electrochemical analysis and DFT calculations suggests that a CN- ligand increases the pKa of the cluster enabling hydrogen production from its Fe(i)-Fe(0) state at pHs much higher and overpotential much lower than its precursor bis-iron hexacarbonyl model which is active in its Fe(0)-Fe(0) state. The formation of a terminal Fe-H species, evidenced by spectroelectrochemistry in organic solvent, via a rate determining proton coupled electron transfer step and protonation of the adjacent azadithiolate, lowers the kinetic barrier leading to diffusion controlled rates of H2 evolution. The stereo-electronic factors enhance its catalytic rate by 3 order of magnitude relative to a bis-iron hexacarbonyl precursor at the same pH and potential.

A personal account on 25 years of scientific literature on [FeFe]-hydrogenase

[FeFe]-hydrogenases are gas-processing metalloenzymes that catalyze H₂ oxidation and proton reduction (H₂ release) in microorganisms. Their high turnover frequencies and lack of electrical overpotential in the hydrogen conversion reaction has inspired generations of biologists, chemists, and physicists to explore the inner workings of [FeFe]-hydrogenase. Here, we revisit 25 years of scientific literature on [FeFe]-hydrogenase and propose a personal account on 'must-read' research papers and review article that will allow interested scientists to follow the recent discussions on catalytic mechanism, O₂ sensitivity, and the in vivo synthesis of the active site cofactor with its biologically uncommon ligands carbon monoxide and cyanide. Focused on-but not restricted to-structural biology and molecular biophysics, we highlight future directions that may inspire young investigators to pursue a career in the exciting and competitive field of [FeFe]-hydrogenase research.

Electronic and Structural Properties of the Double Cubane Iron-Sulfur Cluster

The double-cubane cluster (DCC) refers to an [Fe8S9] iron-sulfur complex that is otherwise only known to exist in nitrogenases. Containing a bridging µ2-S ligand, the DCC in the DCC-containing protein (DCCP) is covalently linked to the protein scaffold via six coordinating cysteine residues. In this study, the nature of spin coupling and the effect of spin states on the cluster’s geometry are investigated computationally. Using density functional theory (DFT) and a broken symmetry (BS) approach to study the electronic ground state of the system, we computed the exchange interaction between the spin-coupled spins of the four FeFe dimers contained in the DCC. This treatment yields results that are in excellent agreement with both computed and experimentally determined exchange parameters for analogously coupled di-iron complexes. Hybrid quantum mechanical (QM)/molecular mechanical (MM) geometry optimizations show that cubane cluster A closest to charged amino acid side chains (Arg312, Glu140, Lys146) is less compact than cluster B, indicating that electrons of the same spin in a charged environment seek maximum separation. Overall, this study provides the community with a fundamental reference for subsequent studies of DCCP, as well as for investigations of other [Fe8S9]-containing enzymes.

Spin Polarization Reveals the Coordination Geometry of the [FeFe] Hydrogenase Active Site in Its CO-Inhibited State

The active site of [FeFe] hydrogenase features a binuclear iron cofactor Fe₂ADT(CO)₃(CN)₂, where ADT represents the bridging ligand aza-propane-dithiolate. The terminal diatomic ligands all coordinate in a basal configuration, and one CO bridges the two irons leaving an open coordination site at which the hydrogen species and the competitive inhibitor CO bind. Externally supplied CO is expected to coordinate in an apical configuration. However, an alternative configuration has been proposed in which, due to ligand rotation, the CN^- bound to the distal Fe becomes apical. Using selective ¹³C isotope labeling of the CN^- and CO_ext ligands in combination with pulsed ¹³C electron-nuclear-nuclear triple resonance spectroscopy, spin polarization effects are revealed that, according to density functional theory calculations, are consistent with only the "unrotated" apical CO_ext configuration.

Asymmetry in the Ligand Coordination Sphere of the [FeFe] Hydrogenase Active Site Is Reflected in the Magnetic Spin Interactions of the Aza-propanedithiolate Ligand

[FeFe] hydrogenases are very active enzymes that catalyze the reversible conversion of molecular hydrogen into protons and electrons. Their active site, the H-cluster, contains a unique binuclear iron complex, [2Fe]_H, with CN^- and CO ligands as well as an aza-propane-dithiolate (ADT) moiety featuring a central amine functionality that mediates proton transfer during catalysis. We present a pulsed ¹³C-ENDOR investigation of the H-cluster in which the two methylene carbons of ADT are isotope labeled with ¹³C. We observed that the corresponding two ¹³C hyperfine interactions are of opposite sign and corroborated this finding using density functional theory calculations. The spin polarization in the ADT ligand is shown to be linked to the asymmetric coordination of the distal iron site with its terminal CN^- and CO ligands. We propose that this asymmetry is relevant for the enzyme reactivity and is related to the (optimal) stabilization of the iron-hydride intermediate in the catalytic cycle.

Electronic Structure of Two Catalytic States of the [FeFe] Hydrogenase H-Cluster As Probed by Pulse Electron Paramagnetic Resonance Spectroscopy

The active site of the [FeFe] hydrogenase (HydA1), the H-cluster, is a 6-Fe cofactor that contains CO and CN^- ligands. It undergoes several different oxidation and protonation state changes in its catalytic cycle to metabolize H₂. Among them, the well-known H_ox state and the recently identified H_hyd state are thought to be directly involved in H₂ activation and evolution, and they are both EPR active with net spin S = 1/2. Herein, we report the pulse electronic paramagnetic spectroscopic investigation of these two catalytic states in Chlamydomonas reinhardtii HydA1 ( CrHydA1). Using an in vitro biosynthetic maturation approach, we site-specifically installed ¹³C into the CO or CN^- ligands and ⁵⁷Fe into the [2Fe]_H subcluster of the H-cluster in order to measure the hyperfine couplings to these magnetic nuclei. For H_ox, we measured ¹³C hyperfine couplings (¹³CO a_iso of 25.5, 5.8, and 4.5 MHz) corresponding to all three CO ligands in the H-cluster. We also observed two ⁵⁷Fe hyperfine couplings (⁵⁷Fe a_iso of ∼17 and 5.7 MHz) arising from the two Fe atoms in the [2Fe]_H subcluster. For H_hyd, we only observed two distinct ¹³CO hyperfine interactions (¹³CO a_iso of 0.16 and 0.08 MHz) and only one for ¹³CN^- (¹³CN a_iso of 0.16 MHz); the couplings to the ¹³CO/¹³CN^- on the distal Fe of [2Fe]_H may be too small to detect. We also observed a very small (<0.3 MHz) ⁵⁷Fe HFI from the labeled [2Fe]_H subcluster and four ⁵⁷Fe HFI from the labeled [4Fe-4S]_H subcluster (⁵⁷Fe a_iso of 7.2, 16.6, 28.2, and 35.3 MHz). These hyperfine coupling constants are consistent with the previously proposed electronic structure of the H-cluster at both H_ox and H_hyd states and provide a basis for more detailed analysis.

Interaction of the H-Cluster of FeFe Hydrogenase with Halides

FeFe hydrogenases catalyze H₂ oxidation and production using an "H-cluster", where two Fe ions are bound by an aza-dithiolate (adt) ligand. Various hypotheses have been proposed (by us and others) to explain that the enzyme reversibly inactivates under oxidizing, anaerobic conditions: intramolecular binding of the N atom of adt, formation of the so-called "Hox/inact" state or nonproductive binding of H₂ to isomers of the H-cluster. Here, we show that none of the above explains the new finding that the anaerobic, oxidative, H₂-dependent reversible inactivation is strictly dependent on the presence of Cl^- or Br^-. We provide experimental evidence that chloride uncompetitively inhibits the enzyme: it reversibly binds to catalytic intermediates of H₂ oxidation (but not to the resting "Hox" state), after which oxidation locks the active site into a stable, saturated, inactive form, the structure of which is proposed here based on DFT calculations. The halides interact with the amine group of the H-cluster but do not directly bind to iron. It should be possible to stabilize the inhibited state in amounts compatible with spectroscopic investigations to explore further this unexpected reactivity of the H-cluster of hydrogenase.

Structural and Kinetic Studies of Intermediates of a Biomimetic Diiron Proton-Reduction Catalyst

One-electron reduction and subsequent protonation of a biomimetic proton-reduction catalyst [FeFe(μ-pdt)(CO)₆] (pdt = propanedithiolate), 1, were investigated by UV-vis and IR spectroscopy on a nano- to microsecond time scale. The study aimed to provide further insight into the proton-reduction cycle of this [FeFe]-hydrogenase model complex, which with its prototypical alkyldithiolate-bridged diiron core is widely employed as a molecular, precious metal-free catalyst for sustainable H₂ generation. The one-electron-reduced catalyst was obtained transiently by electron transfer from photogenerated [Ru(dmb)₃]⁺ in the absence of proton sources or in the presence of acids (dichloro- or trichloroacetic acid or tosylic acid). The reduced catalyst and its protonation product were observed in real time by UV-vis and IR spectroscopy, leading to their structural characterization and providing kinetic data on the electron and proton transfer reactions. 1 features an intact (μ²,κ²-pdt)(μ-H)Fe₂ core in the reduced, 1^-, and reduced-protonated states, 1H, in contrast to the Fe-S bond cleavage upon the reduction of [FeFe(bdt)(CO)₆], 2, with a benzenedithiolate bridge. The driving-force dependence of the rate constants for the protonation of 1^- (k_pt = 7.0 × 10⁵, 1.3 × 10⁷, and 7.0 × 10⁷ M^-1 s^-1 for the three acids used in this study) suggests a reorganization energy >1 eV and indicates that hydride complex 1H is formed by direct protonation of the Fe-Fe bond. The protonation of 1^- is sufficiently fast even with the weaker acids, which excludes a rate-limiting role in light-driven H₂ formation under typical conditions.

The Mössbauer Parameters of the Proximal Cluster of Membrane-Bound Hydrogenase Revisited: A Density Functional Theory Study

An unprecedented [4Fe-3S] cluster proximal to the regular [NiFe] active site has recently been found to be responsible for the ability of membrane-bound hydrogenases (MBHs) to oxidize dihydrogen in the presence of ambient levels of oxygen. Starting from proximal cluster models of a recent DFT study on the redox-dependent structural transformation of the [4Fe-3S] cluster, ⁵⁷Fe Mössbauer parameters (electric field gradients, isomer shifts, and nuclear hyperfine couplings) were calculated using DFT. Our results revise the previously reported correspondence of Mössbauer signals and iron centers in the [4Fe-3S]³⁺ reduced-state proximal cluster. Similar conflicting assignments are also resolved for the [4Fe-3S]⁵⁺ superoxidized state with particular regard to spin-coupling in the broken-symmetry DFT calculations. Calculated ⁵⁷Fe hyperfine coupling (HFC) tensors expose discrepancies in the experimental set of HFC tensors and substantiate the need for additional experimental work on the magnetic properties of the MBH proximal cluster in its reduced and superoxidized redox states.

Theoretical Design by First Principles Molecular Dynamics of a Bioinspired Electrode-Catalyst System for Electrocatalytic Hydrogen Production from Acidified Water

Bacterial di-iron hydrogenases produce hydrogen efficiently from water. Accordingly, we have studied by first-principles molecular-dynamics simulations (FPMD) electrocatalytic hydrogen production from acidified water by their common active site, the [FeFe]H cluster, extracted from the enzyme and linked directly to the (100) surface of a pyrite electrode. We found that the cluster could not be attached stably to the surface via a thiol link analogous to that which attaches it to the rest of the enzyme, despite the similarity of the (100) pyrite surface to the Fe4S4 cubane to which it is linked in the enzyme. We report here a systematic sequence of modifications of the structure and composition of the cluster devised to maintain the structural stability of the pyrite/cluster complex in water throughout its hydrogen production cycle, an example of the molecular design of a complex system by FPMD.

Density Functional Theory Calculation of Bonding and Charge Parameters for Molecular Dynamics Studies on [FeFe] Hydrogenases

We have developed and tested molecular mechanics parameters for [FeS] clusters found in known [FeFe] hydrogenases. Bond stretching, angle bending, dihedral and improper torsion parameters for models of the oxidized and reduced catalytic H-cluster, [4Fe4S](+,2+)Cys4, [4Fe4S](+,2+)Cys3His, and [2Fe2S](+,2+)Cys4, were calculated solely from Kohn-Sham density functional theory and Natural Population Analysis. Circumsphere analysis of the cubane clusters in the energy-minimized structure of the full Clostridium pasteurianum hydrogenase I showed the resulting metallocluster structures to be similar to known cubane structures. All clusters were additionally stable in molecular dynamics simulations over the course of 1.0 ns in the fully oxidized and fully reduced enzyme models. Normal modes calculated by quasiharmonic analysis from the dynamics data show unexpected couplings among internal coordinate motions, which may reflect the effects of the protein structure on metallocluster dynamics.

Artificially maturated [FeFe] hydrogenase from Chlamydomonas reinhardtii: a HYSCORE and ENDOR study of a non-natural H-cluster

Hydrogenases are enzymes that catalyze the oxidation of H2 as well as the reduction of protons to form H2. The active site of [FeFe] hydrogenase is referred to as the "H-cluster" and consists of a "classical" [4Fe-4S] cluster connected via a bridging cysteine thiol group to a unique [2Fe]H sub-cluster, containing CN(-) and CO ligands as well as a bidentate azadithiolate ligand. It has been recently shown that the biomimetic [Fe2(adt)(CO)4(CN)2](2-) (adt(2-) = azadithiolate) complex resembling the diiron sub-cluster can be inserted in vitro into the apo-protein of [FeFe] hydrogenase, which contains only the [4Fe-4S] part of the H-cluster, resulting in a fully active enzyme. This synthetic tool allows convenient incorporation of a variety of diiron mimics, thus generating hydrogenases with artificial active sites. [FeFe] hydrogenase from Chlamydomonas reinhardtii maturated with the biomimetic complex [Fe2(pdt)(CO)4(CN)2](2-) (pdt(2-) = propanedithiolate), in which the bridging adt(2-) ligand is replaced by pdt(2-), can be stabilized in a state strongly resembling the active oxidized (Hox) state of the native protein. This state is EPR active and the signal originates from the mixed valence Fe(I)Fe(II) state of the diiron sub-cluster. Taking advantage of the variant with (15)N and (13)C isotope labeled CN(-) ligands we performed HYSCORE and ENDOR studies on this hybrid protein. The (13)C hyperfine couplings originating from both CN(-) ligands were determined and assigned. Only the (15)N coupling from the CN(-) ligand bound to the terminal iron was observed. Detailed orientation selective ENDOR and HYSCORE experiments at multiple field positions enabled the extraction of accurate data for the relative orientations of the nitrogen and carbon hyperfine tensors. These data are consistent with the crystal structure assuming a g-tensor orientation following the local symmetry of the binuclear sub-cluster.

EPR Spectroscopic Studies of [FeFe]-Hydrogenase Maturation

Proton reduction and H2 oxidation are key elementary reactions for solar fuel production. Hydrogenases interconvert H+ and H2 with remarkable efficiency and have therefore received much attention in this context. For [FeFe]-hydrogenases, catalysis occurs at a unique cofactor called the H-cluster. In this article, we discuss ways in which EPR spectroscopy has elucidated aspects of the bioassembly of the H-cluster, with a focus on four case studies: EPR spectroscopic identification of a radical en route to the CO and CN- ligands of the H-cluster, tracing 57Fe from the maturase HydG into the H-cluster, characterization of the auxiliary Fe-S cluster in HydG, and isotopic labeling of the CN- ligands of HydA for electronic structure studies of its Hox state. Advances in cell-free maturation protocols have enabled several of these mechanistic studies, and understanding H-cluster maturation may in turn provide insights leading to improvements in hydrogenase production for biotechnological applications.

[FeFe]- and [NiFe]-hydrogenase diversity, mechanism, and maturation

The [FeFe]- and [NiFe]-hydrogenases catalyze the formal interconversion between hydrogen and protons and electrons, possess characteristic non-protein ligands at their catalytic sites and thus share common mechanistic features. Despite the similarities between these two types of hydrogenases, they clearly have distinct evolutionary origins and likely emerged from different selective pressures. [FeFe]-hydrogenases are widely distributed in fermentative anaerobic microorganisms and likely evolved under selective pressure to couple hydrogen production to the recycling of electron carriers that accumulate during anaerobic metabolism. In contrast, many [NiFe]-hydrogenases catalyze hydrogen oxidation as part of energy metabolism and were likely key enzymes in early life and arguably represent the predecessors of modern respiratory metabolism. Although the reversible combination of protons and electrons to generate hydrogen gas is the simplest of chemical reactions, the [FeFe]- and [NiFe]-hydrogenases have distinct mechanisms and differ in the fundamental chemistry associated with proton transfer and control of electron flow that also help to define catalytic bias. A unifying feature of these enzymes is that hydrogen activation itself has been restricted to one solution involving diatomic ligands (carbon monoxide and cyanide) bound to an Fe ion. On the other hand, and quite remarkably, the biosynthetic mechanisms to produce these ligands are exclusive to each type of enzyme. Furthermore, these mechanisms represent two independent solutions to the formation of complex bioinorganic active sites for catalyzing the simplest of chemical reactions, reversible hydrogen oxidation. As such, the [FeFe]- and [NiFe]-hydrogenases are arguably the most profound case of convergent evolution. This article is part of a Special Issue entitled: Fe/S proteins: Analysis, structure, function, biogenesis and diseases.

Investigations on the role of proton-coupled electron transfer in hydrogen activation by [FeFe]-hydrogenase

Proton-coupled electron transfer (PCET) is a fundamental process at the core of oxidation-reduction reactions for energy conversion. The [FeFe]-hydrogenases catalyze the reversible activation of molecular H2 through a unique metallocofactor, the H-cluster, which is finely tuned by the surrounding protein environment to undergo fast PCET transitions. The correlation of electronic and structural transitions at the H-cluster with proton-transfer (PT) steps has not been well-resolved experimentally. Here, we explore how modification of the conserved PT network via a Cys → Ser substitution at position 169 proximal to the H-cluster of Chlamydomonas reinhardtii [FeFe]-hydrogenase (CrHydA1) affects the H-cluster using electron paramagnetic resonance (EPR) and Fourier transform infrared (FTIR) spectroscopy. Despite a substantial decrease in catalytic activity, the EPR and FTIR spectra reveal different H-cluster catalytic states under reducing and oxidizing conditions. Under H2 or sodium dithionite reductive treatments, the EPR spectra show signals that are consistent with a reduced [4Fe-4S]H(+) subcluster. The FTIR spectra showed upshifts of νCO modes to energies that are consistent with an increase in oxidation state of the [2Fe]H subcluster, which was corroborated by DFT analysis. In contrast to the case for wild-type CrHydA1, spectra associated with Hred and Hsred states are less populated in the Cys → Ser variant, demonstrating that the exchange of -SH with -OH alters how the H-cluster equilibrates among different reduced states of the catalytic cycle under steady-state conditions.

The cyanide ligands of [FeFe] hydrogenase: pulse EPR studies of (13)C and (15)N-labeled H-cluster

The two cyanide ligands in the assembled cluster of [FeFe] hydrogenase originate from exogenous l-tyrosine. Using selectively labeled tyrosine substrates, the cyanides were isotopically labeled via a recently developed in vitro maturation procedure allowing advanced electron paramagnetic resonance techniques to probe the electronic structure of the catalytic core of the enzyme. The ratio of the isotropic ¹³C hyperfine interactions for the two CN^– ligands—a reporter of spin density on their respective coordinating iron ions—collapses from ≈5.8 for the H_ox form of hydrogenase to <2 for the CO-inhibited form. Additionally, when the maturation was carried out using [¹⁵N]-tyrosine, no features previously ascribed to the nitrogen of the bridging dithiolate ligand were observed suggesting that this bridge is not sourced from tyrosine.

Quantum mechanical methods for the investigation of metalloproteins and related bioinorganic compounds

It is well known that transition metal ions are often bound to proteins, conveying very specific functional properties. In fact, metalloproteins play crucial biological roles in the transport and activation of small molecules such as H2, O2, and N2, as well as in several other biochemical processes. However, even if the presence of transition metals in the active site of proteins allows a very rich biochemistry, the experimental disclosure of structure-activity relationships in metalloproteins is generally difficult exactly because of the presence of transition metals, which are intrinsically characterized by a very versatile and often elusive chemistry. For this reason, computational methods are becoming very popular tools in the characterization of metalloproteins. In particular, since computing power is becoming less and less expensive, due to the continuous technological development of CPUs, the computational tools suited to investigate metalloproteins are becoming more accessible and therefore more commonly used also in molecular biology and biochemistry laboratories. Here, we present the main procedures and computational methods based on quantum mechanics, which are commonly used to study the structural, electronic, and reactivity properties of metalloproteins and related bioinspired compounds, with a specific focus on the practical and technical aspects that must be generally tackled to properly study such biomolecular systems.

Inhibition of [FeFe]-hydrogenases by formaldehyde and wider mechanistic implications for biohydrogen activation

Formaldehyde-a rapid and reversible inhibitor of hydrogen evolution by [FeFe]-hydrogenases-binds with a strong potential dependence that is almost complementary to that of CO. Whereas exogenous CO binds tightly to the oxidized state known as H(ox) but very weakly to a state two electrons more reduced, formaldehyde interacts most strongly with the latter. Formaldehyde thus intercepts increasingly reduced states of the catalytic cycle, and density functional theory calculations support the proposal that it reacts with the H-cluster directly, most likely targeting an otherwise elusive and highly reactive Fe-hydrido (Fe-H) intermediate.

Mixed-valence nickel-iron dithiolate models of the [NiFe]-hydrogenase active site

A series of mixed-valence nickel-iron dithiolates is described. Oxidation of (diphosphine)Ni(dithiolate)Fe(CO)(3) complexes 1, 2, and 3 with ferrocenium salts affords the corresponding tricarbonyl cations [(dppe)Ni(pdt)Fe(CO)(3)](+) ([1](+)), [(dppe)Ni(edt)Fe(CO)(3)](+) ([2](+)) and [(dcpe)Ni(pdt)Fe(CO)(3)](+) ([3](+)), respectively, where dppe = Ph(2)PCH(2)CH(2)PPh(2), dcpe = Cy(2)PCH(2)CH(2)PCy(2), (Cy = cyclohexyl), pdtH(2) = HSCH(2)CH(2)CH(2)SH, and edtH(2) = HSCH(2)CH(2)SH. The cation [2](+) proved unstable, but the propanedithiolates are robust. IR and EPR spectroscopic measurements indicate that these species exist as C(s)-symmetric species. Crystallographic characterization of [3]BF(4) shows that Ni is square planar. Interaction of [1]BF(4) with P-donor ligands (L) afforded a series of substituted derivatives of type [(dppe)Ni(pdt)Fe(CO)(2)L]BF(4) for L = P(OPh)(3) ([4a]BF(4)), P(p-C(6)H(4)Cl)(3) ([4b]BF(4)), PPh(2)(2-py) ([4c]BF(4)), PPh(2)(OEt) ([4d]BF(4)), PPh(3) ([4e]BF(4)), PPh(2)(o-C(6)H(4)OMe) ([4f]BF(4)), PPh(2)(o-C(6)H(4)OCH(2)OMe) ([4g]BF(4)), P(p-tol)(3) ([4h]BF(4)), P(p-C(6)H(4)OMe)(3) ([4i]BF(4)), and PMePh(2) ([4j]BF(4)). EPR analysis indicates that ethanedithiolate [2](+) exists as a single species at 110 K, whereas the propanedithiolate cations exist as a mixture of two conformers, which are proposed to be related through a flip of the chelate ring. Mössbauer spectra of 1 and oxidized S = 1/2 [4e]BF(4) are both consistent with a low-spin Fe(I) state. The hyperfine coupling tensor of [4e]BF(4) has a small isotropic component and significant anisotropy. DFT calculations using the BP86, B3LYP, and PBE0 exchange-correlation functionals agree with the structural and spectroscopic data, suggesting that the SOMOs in complexes of the present type are localized in an Fe(I)-centered d(z(2)) orbital. The DFT calculations allow an assignment of oxidation states of the metals and rationalization of the conformers detected by EPR spectroscopy. Treatment of [1](+) with CN(-) and compact basic phosphines results in complex reactions. With dppe, [1](+) undergoes quasi-disproportionation to give 1 and the diamagnetic complex [(dppe)Ni(pdt)Fe(CO)(2)(dppe)](2+) ([5](2+)), which features square-planar Ni linked to an octahedral Fe center.

O2 reactions at the six-iron active site (H-cluster) in [FeFe]-hydrogenase

Irreversible inhibition by molecular oxygen (O(2)) complicates the use of [FeFe]-hydrogenases (HydA) for biotechnological hydrogen (H(2)) production. Modification by O(2) of the active site six-iron complex denoted as the H-cluster ([4Fe4S]-2Fe(H)) of HydA1 from the green alga Chlamydomonas reinhardtii was characterized by x-ray absorption spectroscopy at the iron K-edge. In a time-resolved approach, HydA1 protein samples were prepared after increasing O(2) exposure periods at 0 °C. A kinetic analysis of changes in their x-ray absorption near edge structure and extended X-ray absorption fine structure spectra revealed three phases of O(2) reactions. The first phase (τ(1) ≤ 4 s) is characterized by the formation of an increased number of Fe-O,C bonds, elongation of the Fe-Fe distance in the binuclear unit (2Fe(H)), and oxidation of one iron ion. The second phase (τ(2) ≈ 15 s) causes a ∼50% decrease of the number of ∼2.7-Å Fe-Fe distances in the [4Fe4S] subcluster and the oxidation of one more iron ion. The final phase (τ(3) ≤ 1000 s) leads to the disappearance of most Fe-Fe and Fe-S interactions and further iron oxidation. These results favor a reaction sequence, which involves 1) oxygenation at 2Fe(H(+)) leading to the formation of a reactive oxygen species-like superoxide (O(2)(-)), followed by 2) H-cluster inactivation and destabilization due to ROS attack on the [4Fe4S] cluster to convert it into an apparent [3Fe4S](+) unit, leading to 3) complete O(2)-induced degradation of the remainders of the H-cluster. This mechanism suggests that blocking of ROS diffusion paths and/or altering the redox potential of the [4Fe4S] cubane by genetic engineering may yield improved O(2) tolerance in [FeFe]-hydrogenase.

A theoretical study on the enhancement of functionally relevant electron transfers in biomimetic models of [FeFe]-hydrogenases

Recent advances aimed at modeling the chemistry of the active site of [FeFe]-hydrogenases (the H-cluster, composed by a catalytic Fe(2)S(2) subcluster and an Fe(4)S(4) portion) have led to the synthesis of binuclear coordination compounds containing a noninnocent organophosphine ligand [2,3-bis(diphenylphosphino)maleic anhydride, bma] that is able to undergo monoelectron reduction, analogously to the tetranuclear Fe(4)S(4) subcluster portion of the H-cluster. However, such a synthetic model was shown to feature negligible electronic communication between the noninnocent ligand and the remaining portion of the cluster, at variance with the enzyme active site. Here, we report a theoretical investigation that shows why the electron transfer observed in the enzyme upon protonation of the catalytic Fe(2)S(2) subsite cannot take place in the bma-containing cluster. In addition, we show that targeted modifications of the bma ligand are sufficient to restore the electronic communication within the model, such that electron density can be more easily withdrawn from the noninnocent ligand, as a result of protonation of the iron centers. Similar results were also obtained with a ligand derived from cobaltocene. The relevance of our findings is discussed from the perspective of biomimetic reproduction of proton reduction to yield molecular hydrogen.

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.

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.

Spin coupling in Roussin's red and black salts

Although DFT calculations have provided a first-order electronic-structural description for Roussin's red and black salts, a detailed study of spin coupling in these species has yet to be reported. Such an analysis is presented here for the first time, based on broken-symmetry density functional theory (DFT, chiefly OLYP/STO-TZP) calculations. Both the Noodleman and Yamaguchi formulas were used to evaluate the Heisenberg coupling constants (J). Three nitrosylated binuclear clusters were studied: [Fe(2)(NO)(2)(Et-HPTB)(O(2)CPh)](2+) (1; Et-HPTB=N,N,N',N'-tetrakis-(N-ethyl-2-benzimidazolylmethyl)-2-hydroxy-1,3-diaminopropane), [Fe(NO)(2){Fe(NO)(NS(3))}-S,S'] (2), and Roussin's red salt anion [Fe(2)(NO)(4)(μ-S)(2)](2-) (3). Although the Heisenberg J for 1 is small (≈10(2) cm(-1)), 2 and 3 exhibit J values that are at least an order of magnitude higher (≈10(3) cm(-1)), where the J values refer to the following Heisenberg spin Hamiltonian: ℋ=JS(A)⋅S(B). For Roussin's black salt anion, [Fe(4)(NO)(7)(μ(3)-S)(3)](-) (4), the Heisenberg spin Hamiltonian describing spin coupling between the {FeNO}(7) unit (S(A)=3/2) and the three {Fe(NO)(2)}(9) units (S(B)=S(C)=S(D)=1/2) in [Fe(4)(NO)(7)(μ(3)-S)(3)](-) was assumed to have the form: ℋ=J(12)(S(A)⋅S(B)+S(A)⋅S(C)+S(A)⋅S(D))+J(22)(S(B)⋅S(C)+S(B)⋅S(D)+S(C)⋅S(D)), in which J(12) corresponds to the interaction between the apical iron and a basal iron, and J(22) refers to that between any two basal iron centers. Although the basal-basal coupling constant J(22) was found to be small (≈10(2) cm(-1)), the apical-basal coupling constant J(12) is some forty times higher (≈4000 cm(-1)). Thus, the nitrosylated iron-sulfur clusters feature some exceptionally high J values relative to the non-nitrosylated {2Fe2S} and {4Fe4S} clusters. An analysis of spin-dependent bonding energies shed light on this curious feature. In essence, the energy difference between the high-spin (i.e., ferromagnetically coupled iron sites) and low-spin (i.e., maximum spin coupling) states of Roussin's salts are indeed rather similar to those of analogous non-nitrosylated iron-sulfur clusters. However, the individual Fe(NO)(x) (x=1, 2) site spins are lower in the nitrosylated systems, resulting in a smaller denominator in both the Noodleman and Yamaguchi formulas for J, which in turn translates into the very high J values.

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.