Use of 1,10-phenanthroline in diiron dithiolate derivatives related to the [Fe-Fe] hydrogenase active site

Treatment of [Fe(2)(micro-pdt)(CO)(6)] (pdt = S(CH(2))(3)S) with 1,10-phenanthroline (phen) in refluxing toluene affords the asymmetric complex [Fe(2)(micro-pdt)(CO)(4)(phen)] (1); the protonation of with HBF(4).OEt(2) in CD(2)Cl(2) at 203 K has been monitored by (1)H NMR.     

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.     

(Alpha-diimine)chromium complexes: molecular and electronic structures; a combined experimental and density functional theoretical study

Dark brown crystals of [Cr( (1)L) 2] ( 1) were obtained from the reaction of [Cr (III)(acac) 3] (acac (-) = 2,4-pentanedionate) with 2 equiv of 2-methyl-1,4-bis(2,6-dimethylphenyl)-1,4-diaza-1,3-butadiene ( (1)L) and 3 equiv of sodium in tetrahydrofuran (thf) under an Ar atmosphere. Complex 1 possesses an S = 1 ground state, which is attained via intramolecular antiferromagnetic coupling between a high-spin Cr (II) ion ( S Cr = 2) and two anionic alpha-diiminato(1-) ligand pi radicals ( (1)L (*)) (1-). The molecular structure of 1 exhibits a distorted tetrahedral, nearly square-planar geometry. The average C-N imine bond length at 1.346 A is characteristic for the pi radical anion ( (1)L (*)) (1-), and therefore, the electronic structure of 1 is best described as [Cr (II)( (1)L (*)) 2]. This has been confirmed by broken symmetry density functional theoretical calculations BS(4,2) (DFT) at the B3LYP level. The reaction of [Cr (III)(acac) 3] with 1 equiv of 2,3-dimethyl-1,4-bis(2,6-diisopropylphenyl)-1,4-diaza-1,3-butadiene ( (2)L) and 1 equiv of Na in thf under Ar yields red-brown crystals of [Cr (III)( (2)L (*))(acac) 2] ( 2) ( S = 1). The oxidation of 2 with 1 equiv of Fc(PF 6) (Fc (+) = ferrocenium) in CH 2Cl 2 affords crystals of [Cr (III)( (2)L (ox))(acac) 2](PF 6) ( 3) ( S = (3)/ 2). The crystal structure determinations of 2 and 3 revealed that 2 contains a neutral, octahedral Cr (III) species [Cr (III)( (2)L (*))(acac) 2], whereas in 3 the ligand is oxidized, yielding an octahedral monocation [Cr (III)( (2)L (ox))(acac) 2] (+). These electronic structures have been confirmed by DFT calculations.     

Comparison between the geometric and electronic structures and reactivities of [FeNO]7 and [FeO2]8 complexes: a density functional theory study

In a previous study, we analyzed the electronic structure of S = 3/2 [FeNO](7) model complexes [Brown et al. J. Am. Chem. Soc. 1995, 117, 715-732]. The combined spectroscopic data and SCF-X alpha-SW electronic structure calculations are best described in terms of Fe(III) (S = 5/2) antiferromagnetically coupled to NO(-) (S = 1). Many nitrosyl derivatives of non-heme iron enzymes have spectroscopic properties similar to those of these model complexes. These NO derivatives can serve as stable analogues of highly labile oxygen intermediates. It is thus essential to establish a reliable density functional theory (DFT) methodology for the geometry and energetics of [FeNO](7) complexes, based on detailed experimental data. This methodology can then be extended to the study of [FeO(2)](8) complexes, followed by investigations into the reaction mechanisms of non-heme iron enzymes. Here, we have used the model complex Fe(Me(3)TACN)(NO)(N(3))(2) as an experimental marker and determined that a pure density functional BP86 with 10% hybrid character and a mixed triple-zeta/double-zeta basis set lead to agreement between experimental and computational data. This methodology is then applied to optimize the hypothetical Fe(Me(3)TACN)(O(2))(N(3))(2) complex, where the NO moiety is replaced by O(2). The main geometric differences are an elongated Fe[bond]O(2) and a steeper Fe[bond]O[bond]O angle in the [FeO(2)](8) complex. The electronic structure of [FeO(2)](8) corresponds to Fe(III) (S = 5/2) antiferromagnetically coupled to O(2)(-) (S = 1/2), and, consistent with the extended bond length, the [FeO(2)](8) unit has only one Fe(III)-O(2)(-) bonding interaction, while the [FeNO](7) unit has both sigma and pi type Fe(III)-NO(-) bonds. This is in agreement with experiment as NO forms a more stable Fe(III)-NO(-) adduct relative to O(2)(-). Although NO is, in fact, harder to reduce, the resultant NO(-) species forms a more stable bond to Fe(III) relative to O(2)(-) due to the different bonding interactions.     

Hydrogen bonding affects the [NiFe] active site of Desulfovibrio vulgaris Miyazaki F hydrogenase: a hyperfine sublevel correlation spectroscopy and density functional theory study

Pulse electron paramagnetic resonance and hyperfine sublevel correlation spectroscopy have been used to investigate nitrogen coordination of the active site of [NiFe] hydrogenase of Desulfovibrio vulgaris Miyazaki F in its oxidized "ready" state. The obtained (14)N hyperfine (A = [+1.32, +1.32, +2.07] MHz) and nuclear quadrupole (e(2)qQ/h = -1.9 MHz, eta = 0.37) coupling constants were assigned to the N(epsilon) of a highly conserved histidine (His88) by studying a hydrogenase preparation in which the histidines were (15)N labeled. The histidine is hydrogen-bonded via its N(epsilon)-H to the nickel-coordinating sulfur of a cysteine (Cys549) that carries an appreciable amount of spin density. Through the hydrogen bond a small fraction of the spin density ( approximately 1%) is delocalized onto the histidine ring giving rise to an isotropic (14)N hyperfine coupling constant of about 1.6 MHz. These conclusions are supported by density functional calculations. The measured (14)N quadrupole coupling constants are related to the polarization of the N(epsilon)-H bond, and the respective hydrogen bond can be classified as being weak.     

Protonation, electrochemical properties and molecular structures of halogen-functionalized diiron azadithiolate complexes related to the active site of iron-only hydrogenases

Diiron complexes [{(micro-SCH2)2NCH2C6H4X}{Fe(CO)2L}2] (L = CO, X = 2-Br, 1; 2-F, 2; 3-Br, 3; L = PMe(3), X = 2-Br, 4) were prepared as biomimetic models of the iron-only hydrogenase active site. The N-protonated species [(NH)]+ClO(4)(-), [(NH)](+)ClO(4)(-) and the micro-hydride diiron complex [4(FeHFe)]+PF(6)(-) were obtained in the presence of proton acids and well characterized. The protonation process of 4 was studied by in-situ IR and NMR spectroscopy, which suggests the formation of the diprotonated species [4(NH)(FeHFe)](2+) in the presence of an excess of proton acid. The molecular structures of 1, [(NH)]+ClO(4)(-), 4 and [4(FeHFe)]+PF(6)(-) were determined by X-ray crystallography. The single-crystal X-ray analysis reveals that an intramolecular H...Br contact (2.82 A) in the crystalline state of [1(NH)]+ClO(4)(-). In the presence of 1-6 equiv of the stronger acid HOTf, complex 1 is readily protonated on the bridged-N atom and can electrochemically catalyze the proton reduction at a relatively mild potential (ca.-1.0 V). Complex 4 is also electrocatalytic active at -1.4 V in the presence of HOTf with formation of the micro-hydride diiron species.     

Modeling [FeFe] hydrogenase: Synthesis and protonation of a diiron dithiolate complex containing a phosphine-N-heterocyclic-carbene ligand

The ligand exchange reaction of I_Me-(CH₂)₂-PPh₂ (I_Me = 1-methyimidazol-2-ylidene) and the hexacarbonyl complex [{Fe₂{μ-S(CH₂)₃S}(CO)₆] (1) resulted in the formation of the chelated complex [{Fe₂{μ-S(CH₂)₃S}(CO)₄(I_Me-(CH₂)₂-PPh₂)] (2). The molecular structure of 2 was confirmed by spectroscopic and X-ray analyses. This complex catalyzes proton reduction. Low temperature NMR studies on the protonation of 2 revealed the formation of a terminal hydride intermediate.     

Mechanism of oxidative shuttling for [2]rotaxane in a Stoddart-Heath molecular switch: density functional theory study with continuum-solvation model

The central component of the programmable molecular switch demonstrated recently by Stoddart and Heath is [2]rotaxane, which consists of a cyclobis-(paraquat-p-phenylene) ring-shaped shuttle [(CBPQT(4+))(PF(6)(-))(4)] encircling a finger and moving between two stations on the finger: tetrathiafulvalene (TTF) and 1,5-dioxynaphthalene (DNP). We report here a quantum mechanics (QM) study of the mechanism by which movement of the ring (and in turn the on-off switching) is controlled by the oxidation-reduction process. We use B3LYP density functional theory to describe how oxidation of the [2]rotaxane components (in using Poisson-Boltzmann continuum-solvation theory for acetonitrile solution) induces the motions associated with switching (translation of the ring). These calculations support the proposal that oxidation occurs on TTF, leading to repulsion between two positive charge centers (TTF(2+) and CBPQT(4+)) that drives the CBPQT(4+) ring from the TTF(2+) station toward the neutral DNP station. The theory also supports the experimental observation that the first and second oxidation potentials are nearly the same (separated by 0.09 eV in the QM). This excellent agreement between the QM and experiment suggests that QM can be useful in designing new systems.     

Hydrogen generation from weak acids: electrochemical and computational studies of a diiron hydrogenase mimic

Extended investigation of electrocatalytic generation of dihydrogen using [(mu-1,2-benzenedithiolato)][Fe(CO)3]2 has revealed that weak acids, such as acetic acid, can be used. The catalytic reduction producing dihydrogen occurs at approximately -2 V for several carboxylic acids and phenols resulting in overpotentials of only -0.44 to -0.71 V depending on the weak acid used. This unusual catalytic reduction occurs at a potential at which the starting material, in the absence of a proton source, does not show a reduction peak. The mechanism for this process and structures for the intermediates have been discerned by electrochemical and computational analysis. These studies reveal that the catalyst is the monoanion of the starting material and an ECEC mechanism occurs.     

Resolving the CO/CN ligand arrangement in CO-inactivated [FeFe] hydrogenase by first principles density functional theory calculations

The currently presumed assignment of CO/CN ligands in the structure of the active cluster in CO-inactivated [FeFe] hydrogenase is shown to be inconsistent with the available IR data in the enzyme from Clostridium pasteurianum I. A different arrangement has the correct qualitative and quantitative features, reproducing the observed line spacing and intensities and the observed line shift consequent to inactivation with labeled 13CO instead of 12CO. The new assignment is also consistent with the observed change from rhombic to axial symmetry of the electron paramagnetic resonance g tensor upon inactivation.     

Mechanistic study of the electrochemical oxygen reduction reaction on Pt(111) using density functional theory

Density functional theory (DFT) was used to study the electrolyte solution effects on the oxygen reduction reaction (ORR) on Pt(111). To model the acid electrolyte, an H(5)O(2)(+) cluster was used. The vibrational proton oscillation modes for adsorbed H(5)O(2)(+) computed at 1711 and 1010 cm(-1), in addition to OH stretching and H(2)O scissoring modes, agree with experimental vibrational spectra for proton formation on Pt surfaces in ultrahigh vacuum. Using the H(5)O(2)(+) model, protonation of adsorbed species was found to be facile and consistent with the activation barrier of proton transfer in solution. After protonation, OOH dissociates with an activation barrier of 0.22 eV, similar to the barrier for O(2) dissociation. Comparison of the two pathways suggests that O(2) protonation precedes dissociation in the oxygen reduction reaction. Additionally, an OH diffusion step following O protonation inhibits the reaction, which may lead to accumulation of oxygen on the electrode surface.     

Synthesis,characterization, and electrochemical properties of diiron azadithiolate complexes related to the active site of [FeFe]-hydrogenases

Treatment of [(μ-SCH₂)₂NPh]Fe₂(CO)₆ (A) with PPh₃ or PPh₂H in the presence of the decarbonylating agent Me₃NO·2H₂O afforded complexes [(μ-SCH₂)₂NPh]Fe₂(CO)₅(PPh₃) (1) and [(μ-SCH₂)₂NPh]Fe₂(CO)₅(PPh₂H) (2) in 87% and 74% yields, respectively. Complexes 1 and 2 were characterized by elemental analysis and various spectroscopic techniques. The molecular structures of 1 and 2 were further determined by X-ray crystallography. In both cases, the monophosphine ligand resides in an axial position of the square-pyramidal Fe atom and trans to the benzene ring of the azadithiolate ligand, in order to minimize steric repulsion. On the basis of electrochemical studies, all these complexes were found to catalyze proton reduction to H₂ in the presence of acetic acid.     

Electron delocalization from the fullerene attachment to the diiron core within the active-site mimics of [FeFe]hydrogenase

Attachment of the redox-active C(60)(H)PPh(2) group modulates the electronic structure of the Fe(2) core in [(μ-bdt)Fe(2)(CO)(5)(C(60)(H)PPh(2))]. The neutral complex is characterized by X-ray crystallography, IR, NMR spectroscopy, and cyclic voltammetry. When it is reduced by one electron, the spectroscopic and density functional theory results indicate that the Fe(2) core is partially spin-populated. In the doubly reduced species, extensive electron communication occurs between the reduced fullerene unit and the Fe(2) centers as displayed in the spin-density plot. The results suggest that the [4Fe4S] cluster within the H cluster provides an essential role in terms of the electronic factor.     

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.     

Ionization potentials of tantalum-carbide clusters: an experimental and density functional theory study

We have used photoionization efficiency spectroscopy to determine the ionization potentials (IP) of the tantalum-carbide clusters, Ta3Cn (n = 1-3) and Ta4Cn (n = 1-4). The ionization potentials follow an overall reduction as the number of carbon atoms increases; however, the trend is not steady as expected from a simple electrostatic argument. Instead, an oscillatory behavior is observed such that clusters with an odd number of carbon atoms have higher IPs and clusters with an even number of carbon atoms have lower IPs, with the Ta4C4 cluster exhibiting the lowest IP. Excellent agreement is found with relative IPs calculated using density functional theory for the lowest energy structures, which are consistent with the development of a 2 x 2 x 2 face-centered nanocrystal. This work shows that IPs may be used as a reliable validation for the geometries of metal-carbide clusters calculated by theory. The variation in IP can also be interpreted qualitatively with application of a simple model based upon isolobal frontier orbitals.     

Synthesis,characterization, and crystal structures of N-functionalized diiron azadithiolate complexes related to the active site of [FeFe]-hydrogenases

A series of N-functionalized diiron azadithiolate complexes, [(µ-SCH₂)₂NCH₂CO₂Me]Fe₂(CO)_5?L [L?=?CO (1); PPh₃ (2); Ph₂PCH₂PPh₂ (3)], as active site models of [FeFe]-hydrogenases has been prepared and characterized. While 1 was prepared by a sequential reaction of (µ-HS)₂Fe₂(CO)₆ with two equiv. of aqueous HCHO, followed by treatment of (µ-HOCH₂S)₂Fe₂(CO)₆ with one equiv. of H₂NCH₂CO₂Me in 46% yield; 2 and 3 were prepared by a carbonyl substitution reaction of 1 with PPh₃ or Ph₂PCH₂PPh₂ in the presence of Me₃NO?·?2H₂O in 90% and 85% yields, respectively. The crystal structures of 1 and 2 revealed that the substituent attached to the bridgehead nitrogen occupies an equatorial position and the PPh₃ ligand resides in an axial position of the square pyramid of Fe2.     

Tris(N-pyrrolidinyl)phosphine substituted diiron dithiolate related to iron-only hydrogenase active site: Synthesis, characterization and electrochemical properties

A tris(N-pyrrolidinyl)phosphine (P(NC₄H₈)₃) monosubstituted complex, [(μ-pdt)Fe₂(CO)₅P(NC₄H₈)₃] (2) was synthesized as a functional model of the hydrogen-producing capability of the iron hydrogenase active site. The structure was fully characterized by X-ray crystallography. IR and electrochemical studies have indicated that the P(NC₄H₈)₃ ligand has better electron-donating ability than that of those phosphine ligands, such as PMe₃, PTA (1,3,5-triaza-7-phosphaadamantane), PMe₂Ph PPh₃, and P(OEt)₃. The electrocatalytic activity of 2 was recorded in CH₃CN in the absence and presence of weak acid, HOAc. The cathodic shift of potential at −1.98 V and the dependence of current on acid concentration have indicated that complex 2 can catalyze the reduction of protons to hydrogen at its Fe⁰Fe^I level in the presence of HOAc. IR spectroelectrochemical experiments are conducted during the reduction of 2 under nitrogen and carbon monoxide, respectively. The formation of a bridging CO group during the reduction of 2 at −1.98 V has been identified using IR spectroelectrochemical techniques, and an electrocatalytic mechanism of 2 consistent with the spectroscopic and electrochemical results is proposed.     

Proton reduction and dihydrogen oxidation on models of the [2Fe]H cluster of [Fe] hydrogenases. A density functional theory investigation

Density functional theory was used to compare reaction pathways for H2 formation and H+ reduction catalyzed by models of the binuclear cluster found in the active site of [Fe] hydrogenases. Terminal H+ binding to an Fe(I)-Fe(I) form, followed by monoelectron reduction and protonation of the di(thiomethyl)amine ligand, can conveniently lead to H2 formation and release, suggesting that this mechanism could be operative within the enzyme active site. However, a pathway that implies the initial formation of Fe(II)-Fe(II) mu-H species and release of H2 from an Fe(II)-Fe(I) form is characterized by only slightly less favored energy profiles. In both cases, H2 formation becomes less favored when taking into account the competition between CN and amine groups for H+ binding, an observation that can be relevant for the design of novel synthetic catalysts. H2 cleavage can take place on Fe(II)-Fe(II) redox species, in agreement with previous proposals [Fan, H.-J.; Hall, M. B. J. Am. Chem. Soc. 2001, 123, 3828] and, in complexes characterized by terminal CO groups, does not need the involvement of an external base. The step in H2 oxidation characterized by larger energy barriers corresponds to the second H+ extraction from the cluster, both considering Fe(II)-Fe(II) and Fe(II)-Fe(III) species. A comparison of the different reaction pathways reveals that H2 formation could involve only Fe(I)-Fe(I), Fe(II)-Fe(I), and Fe(II)-Fe(II) species, whereas Fe(III)-Fe(II) species might be relevant in H2 cleavage.