Coordination sphere flexibility of active-site models for Fe-only hydrogenase: studies in intra- and intermolecular diatomic ligand exchange

Desymmetrized Diiron Azadithiolato Carbonyls: A Step Toward Modeling the Iron-Only Hydrogenases

Condensation of Fe(2)(SH)(2)(CO)(6), acetaldehyde, and (NH(4))(2)CO(3) affords the methyl-substituted azadithiolate Fe(2)[(SCHMe)(2)NH](CO)(6) (1). The complex exists mainly (~95%) as the meso diastereomer, but the d,l diastereoisomers could be detected. DFT calculations predict that the meso isomer would be 2.5 kcal/mol more stable than the d,l isomer due to conventional nonbonding interactions between the methyl groups and the ring hydrogen atoms. Crystallographic analysis of meso-1 confirms that the two methyl groups are equatorial, constraining the diferraazadithiolate bicycle to a conformation that desymmetrizes the diiron center. The lowered symmetry is confirmed by the observation of two (13)C NMR signals in the FeCO region under conditions of fast turnstile rotation at the Fe(CO)(3) groups. The pK(a) value of the amine in 1 is 7.89 (all pK(a)'s determined in MeCN solution), which is similar to a redetermined value for Fe(2)[(SCH(2))(2)NH](CO)(6) (2, pK(a) = 7.98) and only slightly less basic than the tertiary amine Fe(2)[(SCH(2))(2)NMe](CO)(6) (pK(a) = 8.14). Substitution of 1 with PMe(3) proceeded via the intermediacy of two isomers of Fe(2)[(SCHMe)(2)NH](CO)(5)(PMe(3)), affording Fe(2)[(SCHMe)(2)NH](CO)(4)(PMe(3))(2) (3). (31)P NMR spectra confirm that the two PMe(3) ligands in 3 are nonequivalent, consistent with the desymmetrizing effect of the dithiolate. The pK(a) value of the amine in 3 was found to be 11.3. Using triphenylphosphine, we prepared Fe(2)[(SCHMe)(2)NH](CO)(5)(PPh(3)) as a single regioisomer.

Computational definition of a mixed valent Fe(II)Fe(I) model of the [FeFe]hydrogenase active site resting state

Density-functional calculations have been used to examine the electronic structure and bonding in the recently reported complex [(PMe(3))(CO)(2)Fe(mu-pdt)(mu-CO)Fe(CO)(IMes)](+) (1(+), IMes=1,3-bis(2,4,6-trimethylphenyl)-imidazol-2-ylidene). This mixed valent Fe(II)Fe(I) complex features a rotated geometry that places a carbonyl ligand in a semi-bridging position, which makes it an accurate model of the S =(1/2) resting state of the [FeFe]-hydrogenase active site. Calculations indicate that the unpaired electron in this complex lies almost entirely on the rotated iron center, implying that this iron remains in the Fe(I) oxidation state, while the unrotated iron has been oxidized to Fe(II). The frontier molecular orbitals in 1(+) are compared with those in the neutral Fe(I)Fe(I) precursor (PMe(3))(CO)(2)Fe(mu-pdt)(mu-CO)Fe(CO)(IMes) at both its optimized geometry (1) and constrained to a rotated geometry (1(rot)). These theoretical results are used to address the role of the bridging CO ligand in 1(+) and to predict reactivity patterns; they are related back to the intricate biological mechanism of [FeFe]-hydrogenase.

Lewis vs. Brønsted-basicities of diiron dithiolates: spectroscopic detection of the "rotated structure" and remarkable effects of ethane- vs. propanedithiolate

The new complexes Fe2(S2CnH2n)(CO)2(dppv)2 (n = 2, 3; dppv = cis-1,2-C2H2(PPh2)2) form adducts with AlBr3 and B(C6F5)3, which adopt the "rotated structure" proposed for the active site of the Fe-only hydrogenases--the propanedithiolate is significantly more Lewis basic due to nonbonded interactions between the dithiolate strap and the ligands on Fe.

Studies on the Condensation Pathway to and Properties of Diiron Azadithiolate Carbonyls

Reaction of Fe(2)(SH)(2)(CO)(6) and HCHO, which gives Fe(2)[(SCH(2))(2)NH](CO)(6) in the presence of NH(3), affords the possible intermediate Fe(2)(SCH(2)OH)(2)(CO)(6), which has been characterized crystallographically as its axial-equatorial isomer. Fe(2)(SCH(2)OH)(2)(CO)(6) was shown to react with ammonia and amines to give Fe(2)[(SCH(2))(2)NR](CO)(6) (R = H, alkyl). Related hemithioacetal intermediates were generated by treatment of Fe(2)(SH)(2)(CO)(6) with RC(O)C(O)R (R = H, Ph, 4-F-C(6)H(4)) to give cycloadducts. The benzil derivative Fe(2)[S(2)C(2)(OH)(2)Ph(2)](CO)(6), a C(2)-symmetric species, was also characterized crystallographically. The acylated azadithiolate Fe(2)[(SCH(2))(2)NAc](CO)(6) was prepared by reaction of Li(2)Fe(2)S(2)(CO)(6) with (ClCH(2))(2)NC(O)Me. DNMR experiments show that the free energies of activation for rotation of the amide bond are the same for Fe(2)[(SCH(2))(2)NAc](CO)(6) and Fe(2)[(SCH(2))(2)NAc](CO)(4)(PMe(3))(2), which implies that the ligands on the iron centers do not strongly affect the basicity of the nitrogen. As a control, we showed that the thioamide Fe(2)[(SCH(2))(2)NC(S)Me](CO)(6) does exhibit a significantly higher barrier to rotation, attributable to the increased double-bond character of the N-C(S) bond.

Nature of hydrogen interactions with Ni(II) complexes containing cyclic phosphine ligands with pendant nitrogen bases

Studies of the role of proton relays in molecular catalysts for the electrocatalytic production and oxidation of H(2) have been carried out. The electrochemical production of hydrogen from protonated DMF solutions catalyzed by [Ni(P(2)(Ph)N(2)(Ph))(2)(CH(3)CN)](BF(4))(2), 3a (where P(2)(Ph)N(2)(Ph) is 1,3,5,7-tetraphenyl-1,5-diaza-3,7-diphosphacyclooctane), permits a limiting value of the H(2) production rate to be determined. The turnover frequency of 350 s(-1) establishes that the rate of H(2) production for the mononuclear nickel catalyst 3a is comparable to those observed for Ni-Fe hydrogenase enzymes. In the electrochemical oxidation of hydrogen catalyzed by [Ni(P(2)(Cy)N(2)(Bz))(2)](BF(4))(2), 3b (where Cy is cyclohexyl and Bz is benzyl), the initial step is the reversible addition of hydrogen to 3b (K(eq) = 190 atm(-1) at 25 degrees C). The hydrogen addition product exists as three nearly isoenergetic isomers 4A-4C, which have been identified by a combination of one- and two-dimensional (1)H, (31)P, and (15)N NMR spectroscopies as Ni(0) complexes with a protonated amine in each cyclic ligand. The nature of the isomers, together with calculations, suggests a mode of hydrogen activation that involves a symmetrical interaction of a nickel dihydrogen ligand with two amine bases in the diphosphine ligands. Single deprotonation of 4 by an external base results in a rearrangement to [HNi(P(2)(Cy)N(2)(Bz))(2)](BF(4)), 5, and this reaction is reversed by the addition of a proton to the nickel hydride complex. The small energy differences associated with significantly different distributions in electron density and protons within these molecules may contribute to their high catalytic activity.

Chelate control of diiron(I) dithiolates relevant to the [Fe-Fe]- hydrogenase active site

The reaction of Fe2(S2C2H4)(CO)6 with cis-Ph2PCH=CHPPh2 (dppv) yields Fe2(S2C2H4)(CO)4(dppv), 1(CO)4, wherein the dppv ligand is chelated to a single iron center. NMR analysis indicates that in 1(CO)4, the dppv ligand spans axial and basal coordination sites. In addition to the axial-basal isomer, the 1,3-propanedithiolate and azadithiolate derivatives exist as dibasal isomers. Density functional theory (DFT) calculations indicate that the axial-basal isomer is destabilized by nonbonding interactions between the dppv and the central NH or CH2 of the larger dithiolates. The Fe(CO)3 subunit in 1(CO)4 undergoes substitution with PMe3 and cyanide to afford 1(CO)3(PMe3) and (Et4N)[1(CN)(CO)3], respectively. Kinetic studies show that 1(CO)4 reacts faster with donor ligands than does its parent Fe2(S2C2H4)(CO)6. The rate of reaction of 1(CO)4 with PMe3 was first order in each reactant, k = 3.1 x 10(-4) M(-1) s(-1). The activation parameters for this substitution reaction, DeltaH = 5.8(5) kcal/mol and DeltaS = -48(2) cal/deg.mol, indicate an associative pathway. DFT calculations suggest that, relative to Fe2(S2C2H4)(CO)6, the enhanced electrophilicity of 1(CO)4 arises from the stabilization of a "rotated" transition state, which is favored by the unsymmetrically disposed donor ligands. Oxidation of MeCN solutions of 1(CO)3(PMe3) with Cp2FePF6 yielded [Fe2(S2C2H4)(mu-CO)(CO)2(dppv)(PMe3)(NCMe)](PF6)2. Reaction of this compound with PMe3 yielded [Fe2(S2C2H4)(mu-CO)(CO)(dppv)(PMe3)2(NCMe)](PF6)2.

Preparation, structures and electrochemical property of phosphine substituted diiron azadithiolates relevant to the active site of Fe-only hydrogenases

Mono- and di-phosphine diiron azadithiolate complexes [{(mu-SCH(2))(2)N(4-NO(2)C(6)H(4))}Fe(2)(CO)(5)(PMe(3))] (2), [{(mu-SCH(2))(2)N(4-NO(2)C(6)H(4))}{Fe(CO)(2)L}(2)] (3, L=PMe(3); 4, PMe(2)Ph) and the mu-hydride diiron complex [3(FeHFe)](+)[PF(6)](-) were prepared as biomimetic models of the active site of Fe-only hydrogenases. The complexes 2-4 and [3(FeHFe)](+)[PF(6)](-) were characterized by IR, (31)P, (1)H and (13)C NMR spectra and their molecular structures were determined by single crystal X-ray analyses. The PMe(3) ligand in complex 2 lies on the basal position. The PMe(3)-disubstituted complex 3 exists as two configuration isomers, transoid basal/basal and apical/basal, in the crystalline state, while two PMe(2)Ph ligands of 4 are in an apical/basal orientation. The variable temperature (31)P NMR spectra of 2 and 3 were made to have an insight into the existence of the possible conformation isomers of 2 and 3 in solution. The [3(FeHFe)](+) cation possesses the sole transoid ba/ba geometry as other reported mu-hydride diiron analogues. The electrocatalytic property of {(mu-SCH(2))(2)NC(6)H(5)}[Fe(CO)(2)PMe(3)](2) (5) was studied for proton reduction in the presence of HOAc.

Models of the iron-only hydrogenase: Structural studies of chelating diphosphine complexes [Fe2(CO)4(µ-pdt)(κ2P,P′-diphosphine)]

Six chelating diphosphine complexes, [Fe2(CO)4(micro-pdt)(kappa2P,P'-diphosphine)], have been crystallographically characterised allowing differences between basal-apical and dibasal conformations to be analysed.

The active site of the [FeFe]-hydrogenase from Desulfovibrio desulfuricans. II. Redox properties, light sensitivity and CO-ligand exchange as observed by infrared spectroscopy

In [FeFe]-hydrogenases, the H cluster (hydrogen-activating cluster) contains a di-iron centre ([2Fe]H subcluster, a (L)(CO)(CN)Fe(mu-RS2)(mu-CO)Fe(CysS)(CO)(CN) group) covalently attached to a cubane iron-sulphur cluster ([4Fe-4S]H subcluster). The Cys-thiol functions as the link between one iron (called Fe1) of the [2Fe]H subcluster and one iron of the cubane subcluster. The other iron in the [2Fe]H subcluster is called Fe2. The light sensitivity of the Desulfovibrio desulfuricans enzyme in a variety of states has been studied with infrared (IR) spectroscopy. The aerobic inactive enzyme (H(inact) state) and the CO-inhibited active form (H(ox)-CO state) were stable in light. Illumination of the H(ox) state led to a kind of cannibalization; in some enzyme molecules the H cluster was destroyed and the released CO was captured by the H clusters in other molecules to form the light-stable H(ox)-CO state. Illumination of active enzyme under 13CO resulted in the complete exchange of the two intrinsic COs bound to Fe2. At cryogenic temperatures, light induced the photodissociation of the extrinsic CO and the bridging CO of the enzyme in the H(ox)-CO state. Electrochemical redox titrations showed that the enzyme in the H(inact) state converts to the transition state (H(trans)) in a reversible one-electron redox step (E (m, pH 7) = -75 mV). IR spectra demonstrate that the added redox equivalent not only affects the [4Fe-4S]H subcluster, but also the di-iron centre. Enzyme in the H(trans) state reacts with extrinsic CO, which binds to Fe2. The H(trans) state converts irreversibly into the H(ox) state in a redox-dependent reaction most likely involving two electrons (E (m, pH 7) = -261 mV). These electrons do not end up on any of the six Fe atoms of the H cluster; the possible destiny of the two redox equivalents is discussed. An additional reversible one-electron redox reaction leads to the H(red) state (E (m, pH 7) = -354 mV), where both Fe atoms of the [2Fe]H subcluster have the same formal oxidation state. The possible oxidation states of Fe1 and Fe2 in the various enzyme states are discussed. Low redox potentials (below -500 mV) lead to destruction of the [2Fe]H subcluster.

Probing the mechanism of hypoxia selectivity of copper bis(thiosemicarbazonato) complexes: DFT calculation of redox potentials and absolute acidities in solution

Density functional theory (DFT) calculations have been performed using the uB3LYP/6-31++G(d,p) model to calculate the solution phase one-electron reduction potentials (E(calc)) and absolute pKa values of a series of copper bis(thiosemicarbazonato) complexes. The effects of solvation in water and dimethylsulfoxide (DMSO) are incorporated as a self-consistent reaction field (SCRF) using the integral equation formalism polarisable continuum model (IEFPCM) and are found to be essential for quantitative agreement with an average error in E(calc) of -0.02 V compared to experiment. The bonding and spin densities are examined through the use of Natural Bond Order analysis and the results used to rationalise the calculated and observed reduction potentials. Calculated estimates of pKa values of several copper(II) species are presented and their implications for the mechanisms of transport and trapping within hypoxic cells are considered. Reduction is found to be a prerequisite for protonation of the complexes which suggests their transport in the blood stream as neutral species, and the mechanistic sequence is identified as a sequential electrochemical-chemical (EC) process. The complex equilibria of protonation, reoxidation and dissociation are discussed and the copper(I) diprotonated, cationic complex of diacetyl bis(4-methyl-3-thiosemicarbazonato)copper(II), Cu(I)ATSMH2(+), is identified as a possible candidate for the initial species trapped in hypoxic cells.

Electrocatalytic proton reduction by dithiolate-bridged diiron carbonyl complexes: a connection to the H-cluster?

Spectroscopic and electrochemical investigation of electrocatalytic proton reduction by Fe(2)(mu-pdt)(CO)(6), 1, have been interpreted in terms of a reaction scheme involving sequential electron-proton reactions to give a two-electron, two-proton product that undergoes rate-limiting dihydrogen elimination. Further reduction, at slightly higher negative potentials, gives a more reactive product and this process dominates reactions conducted at higher acid concentrations. Inhibition of the electrocatalytic reaction by CO is due to the more efficient loss of catalyst and this is best modelled by a reaction that is second order in terms of 1(-). During electrocatalytic proton reduction a new species is observed, which features a bridging CO group and the wavenumbers of the nu(CO) modes of the terminally bound carbonyl groups are similar to those of the carbonyl groups bound to the oxidized form of the H-cluster.

Proton Electroreduction Catalyzed by Cobaloximes: Functional Models for Hydrogenases

Cobaloximes have been examined as electrocatalysts for proton reduction in nonaqueous solvent in the presence of triethylammonium chloride. [Co(III)(dmgH)2pyCl], working at moderate potentials (-0.90 V/(Ag/AgCl/3 mol x L(-1) NaCl) and in neutral conditions, is a promising catalyst as compared to other first-row transition metal complexes which generally function at more negative potentials and/or at lower pH. More than 100 turnovers can be achieved during controlled-potential electrolysis without detectable degradation of the catalyst. Cyclic voltammograms simulation is consistent with a heterolytic catalytic mechanism and allowed us to extract related kinetic parameters. Introduction of an electron-donating (electron-withdrawing) substituent in the axial pyridine ligand significantly increases (decreases) the rate constant of the catalytic cycle determining step. This effect linearly correlates with the Hammet coefficients of the introduced substituents. The influence of the equatorial glyoxime ligand was also investigated and the capability of the stabilized BF2-bridged species [Co(dmgBF2)2(OH2)2] for electrocatalyzed hydrogen evolution confirmed.

Dissecting the Intimate Mechanism of Cyanation of {2Fe3S} Complexes Related to the Active Site of All-Iron Hydrogenases by DFT Analysis of Energetics, Transition States, Intermediates and Products in the Carbonyl Substitution Pathway

A bridging carbonyl intermediate with key structural elements of the diiron sub-site of all-iron hydrogenase has been experimentally observed in the CN/CO substitution pathway of the {2Fe3S} carbonyl precursor, [Fe(2)(CO)(5){MeSCH(2)C(Me)(CH(2)S)(2)}]. Herein we have used density functional theory (DFT) to dissect the overall substitution pathway in terms of the energetics and the structures of transition states, intermediates and products. We show that the formation of bridging CO transitions states is explicitly involved in the intimate mechanism of dicyanation. The enhanced rate of monocyanation of {2Fe3S} over the {2Fe2S} species [Fe(2)(CO)(6){CH(2)(CH(2)S)(2)}] is found to rest with the ability of the thioether ligand to both stabilise a mu-CO transition state and act as a good leaving group. In contrast, the second cyanation step of the {2Fe3S} species is kinetically slower than for the {2Fe2S} monocyanide because the Fe2 atom is deactivated by coordination of the electron-donating thioether group. In addition, hindered rotation and the reaction coordinate of the approaching CN(-) group, are other factors which explain reactivity differences in {2Fe2S} and {2Fe3S} systems. The intermediate species formed in the second cyanation step of {2Fe3S} species is a mu-CO species, confirming the structural assignment made on the basis of FT-IR data (S. J. George, Z. Cui, M. Razavet, C. J. Pickett, Chem. Eur. J. 2002, 8, 4037-4046). In support of this we find that computed and experimental IR frequencies of structurally characterised {2Fe3S} species and those of the bridging carbonyl intermediate are in excellent agreement. In a wider context, the study may provide some insight into the reactivity of dinuclear systems in which neighbouring group on-off coordination plays a role in substitution pathways.

Synthesis, Structures and Electrochemical Properties of Nitro- and Amino-Functionalized Diiron Azadithiolates as Active Site Models of Fe-Only Hydrogenases

Complex [[(mu-SCH2)2N(4-NO2C6H4)]Fe2(CO)6] (4) was prepared by the reaction of the dianionic intermediate [(mu-S)2Fe2(CO)6](2-) and N,N-bis(chloromethyl)-4-nitroaniline as a biomimetic model of the active site of Fe-only hydrogenase. The reduction of 4 by Pd-C/H2 under a neutral condition afforded complex [[(mu-SCH2)2N(4-NH2C6H4)]Fe2(CO)6] (5) in 67 % yield. Both complexes were characterized by IR, 1H and 13C NMR spectroscopy and MS spectrometry. The molecular structure of 4, as determined by X-ray analysis, has a butterfly 2Fe2S core and the aryl group on the bridged-N atom slants to the Fe(2) site. Cyclic voltammograms of 4 and 5 were studied to evaluate their redox properties. It was found that complex 4 catalyzed electrochemical proton reduction in the presence of acetic acid. A plausible mechanism of the electrocatalytic proton reduction is discussed.

Enzymatic mechanism of Fe-only hydrogenase: density functional study on H-H making/breaking at the diiron cluster with concerted proton and electron transfers

The mechanism of the enzymatic hydrogen bond forming/breaking (2H(+) + 2e<==>H(2)) and the plausible charge and spin states of the catalytic diiron subcluster [FeFe](H) of the H cluster in Fe-only hydrogenases are probed computationally by the density functional theory. It is found that the active center [FeFe](H) can be rationally simulated as [[H](CH(3)S)(CO)(CN(-))Fe(p)(CO(b))(mu-SRS)Fe(d)(CO)(CN(-))L], where the monovalence [H] stands for the [4Fe4S](H)(2+) subcluster bridged to the [FeFe](H) moiety, (CH(3)S) represents a Cys-S, and (CO(b)) represents a bridging CO. L could be a CO, H(2)O, H(-), H(2), or a vacant coordination site on Fe(d). Model structures of possible redox states are optimized and compared with the X-ray crystallographic structures and FTIR experimental data. On the basis of the optimal structures, we study the most favorable path of concerted proton transfer and electron transfer in H(2)-forming/breaking reactions at [FeFe](H). Previous mechanisms derived from quantum chemical computations of Fe-only hydrogenases (Cao, Z.; Hall, M. B. J. Am. Chem. Soc. 2001, 123, 3734; Fan, H.; Hall, M. B. J. Am. Chem. Soc. 2001, 123, 3828) involved an unidentified bridging residue (mu-SRS), which is either a propanedithiolate or dithiomethylamine. Our proposed mechanism, however, does not require such a ligand but makes use of a shuttle of oxidation states of the iron atoms and a reaction site between the two iron atoms. Therefore, the hydride H(b)(-) (bridged to Fe(p) and Fe(d)) and eta(2)-H(2) at Fe(p) or Fe(d) most possibly play key roles in the dihydrogen reversible oxidation at the [FeFe](H) active center. This suggested way of H(2) formation/splitting is reminiscent of the mechanism of [NiFe] hydrogenases and therefore would unify the mechanisms of the two related enzymes.

Density functional theory investigation of the active site of [Fe]-hydrogenases: effects of redox state and ligand characteristics on structural, electronic, and reactivity properties of complexes related to the [2Fe]H subcluster

The effects of redox state and ligand characteristics on structural, electronic, and reactivity properties of complexes related to the [2Fe](H) subcluster of [Fe]-hydrogenases have been investigated by DFT calculations and compared with experimental and theoretical data obtained investigating both the enzyme and synthetic model complexes. Our results show that Fe(II)Fe(II) species characterized by OH or H(2)O groups terminally coordinated to the iron atom distal to the terminal sulfur ligand (Fe(d)) are less stable than corresponding mu-OH or mu-H(2)O species, suggesting that the latter are destabilized or kinetically inaccessible in the enzyme. In addition, results obtained investigating Fe(I)Fe(I) and Fe(II)Fe(I) complexes show that structure and relative stability of species characterized by a mu-CO group are significantly affected by the electronic properties of the ligands coordinated to the iron atoms. The investigation of reaction pathways for H(2) activation confirms and extends a previous hypothesis indicating that H(2) can be cleaved on Fe(II)Fe(II) species. In particular, even though [Fe]-hydrogenases are proposed to bind and activate H(2) at a single iron center, the comparison of our data with experimental results obtained studying synthetic complexes (Zhao, X.; Georgakaki, I. P.; Miller, M. L.; Mejia-Rodriguez, R.; Chiang, C.-Y.; Darensbourg, M. Y. Inorg. Chem. 2002, 41, 3917) suggests that activation paths involving both metal ions are also possible. Moreover, mu-H Fe(II)Fe(I) complexes are predicted to correspond to stable species and might be formed in the enzyme catalytic cycle.