FeMo cofactor of nitrogenase: a density functional study of states M(N), M(OX), M(R), and M(I).

The M(N) S = (3)/(2) resting state of the FeMo cofactor of nitrogenase has been proposed to have metal-ion valencies of either Mo(4+)6Fe(2+)Fe(3+) (derived from metal hyperfine interactions) or Mo(4+)4Fe(2+)3Fe(3+) (from M?ssbauer isomer shifts). Spin-polarized broken-symmetry (BS) density functional theory (DFT) calculations have been undertaken to determine which oxidation level best represents the M(N) state and to provide a framework for understanding its energetics and spectroscopy. For the Mo(4+)6Fe(2+)Fe(3+) oxidation state, the spin coupling pattern for several spin state alignments compatible with S = (3)/(2) were generated and assessed by energy and geometric criteria. The most likely BS spin state is composed of a Mo3Fe cluster with spin S(a) = 2 antiferromagnetically coupled to a 4Fe' cluster with spin S(b) = (7)/(2). This state has a low DFT energy for the isolated FeMoco cluster and the lowest energy when the interaction with the protein and solvent environment is included. This spin state also displays calculated metal hyperfine and M?ssbauer isomer shifts compatible with experiment, and optimized geometries that are in excellent agreement with the protein X-ray data. Our best model for the actual spin-coupled state within FeMoco alters this BS state by a slight canting of spins and is analogous in several respects to that found in the 8Fe P-cluster in the same protein. The spin-up and spin-down components of the LUMO contain atomic contributions from Mo(4+) and the homocitrate and from the central prismane Fe sites and muS(2) atoms, respectively. This qualitative picture of the accepting orbitals for M(N) is consistent with observations from M?ssbauer spectra of the one-electron reduced states. Similar calculations for the Mo(4+)4Fe(2+)3Fe(3+) oxidation state yield results that are in poorer agreement with experiment. Using the Mo(4+)6Fe(2+)Fe(3+) oxidation level as the most plausible resting state, the geometric, electronic and energetic properties of the one-electron redox transition to the oxidized state, M(OX), catalytically observed M(R) and radiolytically reduced M(I) states have also been explored.     

Metal substitution in the active site of nitrogenase MFe(7)S(9) (M = Mo(4+), V(3+), Fe(3+))

The unifying view that molybdenum is the essential component in nitrogenase has changed over the past few years with the discovery of a vanadium-containing nitrogenase and an iron-only nitrogenase. The principal question that has arisen for the alternative nitrogenases concerns the structures of their corresponding cofactors and their metal-ion valence assignments and whether there are significant differences with that of the more widely known molybdenum-iron cofactor (FeMoco). Spin-polarized broken-symmetry (BS) density functional theory (DFT) calculations are used to assess which of the two possible metal-ion valence assignments (4Fe(2+)4Fe(3+) or 6Fe(2+)2Fe(3+)) for the iron-only cofactor (FeFeco) best represents the resting state. For the 6Fe(2+)2Fe(3+) oxidation state, the spin coupling pattern for several spin state alignments compatible with S = 0 were generated and assessed by energy criteria. The most likely BS spin state is composed of a 4Fe cluster with spin S(a) = (7)/(2) antiferromagnetically coupled to a 4Fe' cluster with spin S(b) = (7)/(2). This state has the lowest DFT energy for the isolated FeFeco cluster and displays calculated M?ssbauer isomer shifts consistent with experiment. Although the S = 0 resting state of FeFeco has recently been proposed to have metal-ion valencies of 4Fe(2+)4Fe(3+) (derived from experimental M?ssbauer isomer shifts), our isomer shift calculations for the 4Fe(2+)4Fe(3+) oxidation state are in poorer agreement with experiment. Using the Mo(4+)6Fe(2+)Fe(3+) oxidation level of the cofactor as a starting point, the structural consequences of replacement of molybdenum (Mo(4+)) with vanadium (V(3+)) or iron (Fe(3+)) in the cofactor have been investigated. The size of the cofactor cluster shows a dependency on the nature of the heterometal and increases in the order FeMoco < FeVco < FeFeco.     

High-nuclearity sulfide-rich molybdenum[bond]iron[bond]sulfur clusters: reevaluation and extension

High-nuclearity Mo[bond]Fe[bond]S clusters are of interest as potential synthetic precursors to the MoFe(7)S(9) cofactor cluster of nitrogenase. In this context, the synthesis and properties of previously reported but sparsely described trinuclear [(edt)(2)M(2)FeS(6)](3-) (M = Mo (2), W (3)) and hexanuclear [(edt)(2)Mo(2)Fe(4)S(9)](4-) (4, edt = ethane-1,2-dithiolate; Zhang, Z.; et al. Kexue Tongbao 1987, 32, 1405) have been reexamined and extended. More accurate structures of 2-4 that confirm earlier findings have been determined. Detailed preparations (not previously available) are given for 2 and 3, whose structures exhibit the C(2) arrangement [[(edt)M(S)(mu(2)-S)(2)](2)Fe(III)](3-) with square pyramidal Mo(V) and tetrahedral Fe(III). Oxidation states follow from (57)Fe M?ssbauer parameters and an S = (3)/(2) ground state from the EPR spectrum. The assembly system 2/3FeCl(3)/3Li(2)S/nNaSEt in methanol/acetonitrile (n = 4) affords (R(4)N)(4)[4] (R = Et, Bu; 70-80%). The structure of 4 contains the [Mo(2)Fe(4)(mu(2)-S)(6)(mu(3)-S)(2)(mu(4)-S)](0) core, with the same bridging pattern as the [Fe(6)S(9)](2-) core of [Fe(6)S(9)(SR)(2)](4-) (1), in overall C(2v) symmetry. Cluster 4 supports a reversible three-member electron transfer series 4-/3-/2- with E(1/2) = -0.76 and -0.30 V in Me(2)SO. Oxidation of (Et(4)N)(4)[4] in DMF with 1 equiv of tropylium ion gives [(edt)(2)Mo(2)Fe(4)S(9)](3-) (5) isolated as (Et(4)N)(3)[5].2DMF (75%). Alternatively, the assembly system (n = 3) gives the oxidized cluster directly as (Bu(4)N)(3)[5] (53%). Treatment of 5 with 1 equiv of [Cp(2)Fe](1+) in DMF did not result in one-electron oxidation but instead produced heptanuclear [(edt)(2)Mo(2)Fe(5)S(11)](3-) (6), isolated as the Bu(4)N(+)salt (38%). Cluster 6 features the previously unknown core Mo(2)Fe(5)(mu(2)-S)(7)(mu(3)-S)(4) in molecular C(2) symmetry. In 4-6, the (edt)MoS(3) sites are distorted trigonal bipramidal and the FeS(4) sites are distorted tetrahedral with all sulfide ligands bridging. M?ssbauer spectroscopic data for 2 and 4-6 are reported; (mean) iron oxidation states increase in the order 4 < 5 approximately 1 < 6 approximately 2. Redox and spectroscopic data attributed earlier to clusters 2 and 4 are largely in disagreement with those determined in this work. The only iron and molybdenum[bond]iron clusters with the same sulfide content as the iron[bond]molybdenum cofactor of nitrogenase are [Fe(6)S(9)(SR)(2)](4-) and [(edt)(2)Mo(2)Fe(4)S(9)](3-)(,4-).     

Density functional calculations on the binding of dinitrogen to the FeFe cofactor in Fe-only nitrogenase: FeFeco(mu 6-N2) as intermediate in nitrogen fixation

The geometries and stabilities of the FeFe cofactor at different oxidation states and its complexes with N(2) have been determined by density functional calculations. These calculations support an EPR-inactive resting state of the FeFe cofactor with four Fe(2+) and four Fe(3+) sites (4Fe(2+)4Fe(3+)). FeFeco(mu(6)-N(2)) with a central dinitrogen ligand is predicted to be the most stable complex of the FeFe cofactor with N(2). It is easily formed by penetration of N(2) into the trigonal Fe(6) prism of the FeFe cofactor with an approximate barrier of 4 kcal mol(-1). The present DFT results suggest that an FeFeco(mu(6)-N(2)) entity is a plausible intermediate in dinitrogen fixation by nitrogenase. CO is calculated to bind even more strongly than N(2) to the FeFe cofactor so that CO may inhibit the reduction of nitrogen by Fe-only nitrogenase.     

Synthesis of a molecular Mo2Fe6S9 cluster with the topology of the PN cluster of nitrogenase by rearrangement of an edge-bridged Mo2Fe6S8 double cubane

The structures of the P cluster and cofactor cluster of nitrogenase are well-defined crystallographically. They have been obtained only by biosynthesis; their chemical synthesis remains a challenge. Synthetic routes are sought to the P cluster in the P(N) state in which two cuboidal Fe(3)S(3) units are connected by a mu(6)-S atom and two Fe-(mu(2)-S(Cys))-Fe bridges. A reaction scheme affording a Mo(2)Fe(6)S(9) cluster in molecular form having the topology of the P(N) cluster has been devised. Reaction of the single cubane [(Tp)MoFe(3)S(4)Cl(3)](1)(-) with PEt(3) gives [(Tp)MoFe(3)S(4)(PEt(3))(3)](1+) (2), which upon reduction with BH(4)(-) affords the edge-bridged all-ferrous double cubane [(Tp)(2)Mo(2)Fe(6)S(8)(PEt(3))(4)] (4) (Tp = tris(pyrazolylhydroborate(1-)). Treatment of 4 with 3 equiv of HS(-) produces [(Tp)(2)Mo(2)Fe(6)S(9)(SH)(2)](3)(-) (7) as the Et(4)N(+) salt in 86% yield. The structure of 7 is built of two (Tp)MoFe(3)(mu(3)-S)(3) cuboidal fragments bridged by two mu(2)-S atoms and one mu(6)-S atom in an arrangement of idealized C(2) symmetry. The cluster undergoes three one-electron oxidation reactions and is oxidatively cleaved by p-tolylthiol to [(Tp)MoFe(3)S(4)(S-p-tol)(3)](2)(-) and by weak acids to [(Tp)MoFe(3)S(4)(SH)(3)](2-). The cluster core of 7 has the bridging pattern [Mo(2)Fe(6)(mu(2)-S)(2)(mu(3)-S)(6)(mu(6)-S)](1+) with the probable charge distribution [Mo(3+)(2)Fe(2+)(5)Fe(3+)S(9)](1+). Cluster 7 is a topological analogue of the P(N) cluster but differs in having two heteroatoms and two Fe-(mu(2)-S)-Fe instead of two Fe-(mu(2)-S(Cys))-Fe bridges. A best-fit superposition of the two cluster cores affords a weighted rms deviation in atom positions of 0.38 A. Cluster 7 is the first molecular topological analogue of the P(N) cluster. This structure had been prepared previously only as a fragment of complex high-nuclearity Mo-Fe-S clusters.     

簇离子Mo_3S_4M~（n＋）中M与Mo成键作用及电子光谱（M＝Fe、Ni，n＝4；M＝Cu，n＝5）

应用ＩＮＤＯ／Ｓ半经验量子化学方法，对簇合物离子Ｍｏ３Ｓ和Ｍｏ３Ｓ４Ｍｎ＋（Ｍ＝Ｆｅ、Ｎｉ，ｎ＝４；Ｍ＝Ｃｕ，ｎ＝５）分别进行分子轨道计算。根据计算得到的簇离子中的原子表观电荷和成键指标，说明Ｆｅ、Ｎｉ、Ｃｕ＋与Ｍｏ３Ｓ成键作用的相对强度依次是Ｆｅ－Ｍｏ＞Ｎｉ－Ｍｏ＞Ｃｕ＋－Ｍｏ。比较了用含组态作用的ＩＮＤＯ／Ｓ方法计算得到的电子跃迁能与实验得到的电子吸收光谱值，并讨论了吸收峰归属情况。对于Ｍ为Ｆｅ、Ｎｉ的簇离子Ｍｏ３Ｓ４Ｍ４＋，最低能量的电子跃迁吸收峰起源于异金属间电荷转移跃迁（ＭＭ’ＣＴ）；而Ｍｏ３Ｓ４Ｃｕ（５＋）簇离子观察到的吸收峰主要是Ｍｏ３Ｓ芯的局域内电荷转移跃迁。根据理论计算结果，由Ｃｕ＋离子到Ｍｏ３Ｓ的电荷转移跃迁谱线，大约在４６０００ｃｍ－１以上才能观察到吸收峰。从Ｍｏ３Ｓ４Ｆｅ４＋次低能量吸收峰的实验值１６６００ｃｍ－１和理论值１６５００ｃｍ－１与Ｍｏ３Ｓ的最低能量吸收峰的实验值１６６００ｃｍ－１和理论值１６９００ｃｍ－１比较，表明无论从理论上或实验上都能证实簇离子Ｍｏ３Ｓ４Ｆｅ４＋在能量为１６６００ｃｍ－１处的吸收峰是起因于Ｍｏ３Ｓ芯的局域内电荷转移跃迁。     

Single- and double-cubane clusters in the multiple oxidation states [VFe(3)S(4)](3+,2+,1+)

A new series of cubane-type [VFe(3)S(4)](z)() clusters (z = 1+, 2+, 3+) has been prepared as possible precursor species for clusters related to those present in vanadium-containing nitrogenase. Treatment of [(HBpz(3))VFe(3)S(4)Cl(3)](2)(-) (2, z = 2+), protected from further reaction at the vanadium site by the tris(pyrazolyl)hydroborate ligand, with ferrocenium ion affords the oxidized cluster [(HBpz(3))VFe(3)S(4)Cl(3)](1)(-) (3, z = 3+). Reaction of 2 with Et(3)P results in chloride substitution to give [(HBpz(3))VFe(3)S(4)(PEt(3))(3)](1+) (4, z = 2+). Reaction of 4 with cobaltocene reduced the cluster with formation of the edge-bridged double-cubane [(HBpz(3))(2)V(2)Fe(6)S(8)(PEt(3))(4)] (5, z = 1+, 1+), which with excess chloride underwent ligand substitution to afford [(HBpz(3))(2)V(2)Fe(6)S(8)Cl(4)](4)(-) (6, z = 1+, 1+). X-ray structures of (Me(4)N)[3], [4](PF(6)), 5, and (Et(4)N)(4)[6] x 2MeCN are described. Cluster 5 is isostructural with previously reported [(Cl(4)cat)(2)(Et(3)P)(2)Mo(2)Fe(6)S(8)(PEt(3))(4)] and contains two VFe(3)S(4) cubanes connected across edges by a Fe(2)S(2) rhomb in which the bridging Fe-S distances are shorter than intracubane Fe-S distances. M?ssbauer (2-5), magnetic (2-5), and EPR (2, 4) data are reported and demonstrate an S = 3/2 ground state for 2 and 4 and a diamagnetic ground state for 3. Analysis of (57)Fe isomer shifts based on an empirical correlation between shift and oxidation state and appropriate reference shifts results in two conclusions. (i) The oxidation 2 --> 3 + e(-) results in a change in electron density localized largely or completely on the Fe(3) subcluster and associated sulfur atoms. (ii) The most appropriate charge distributions are [V(3+)Fe(3+)Fe(2+)(2)S(4)](2+) (Fe(2.33+)) for 1, 2, and 4 and [V(3+)Fe(3+)(2)Fe(2+)S(4)](3+) (Fe(2.67+)) for 3 and [V(2)Fe(6)S(8)(SEt)(9)](3+). Conclusion i applies to every MFe(3)S(4) cubane-type cluster thus far examined in different redox states at parity of cluster ligation. The formalistic charge distributions are regarded as the best current approximations to electron distributions in these delocalized species. The isomer shifts require that iron atoms are mixed-valence in each cluster.     

Hydroxide-promoted core conversions of molybdenum-iron-sulfur edge-bridged double cubanes: oxygen-ligated topological PN clusters

The occurrence of a heteroatom X (C, N, or O) in the MoFe7S9X core of the iron-molybdenum cofactor of nitrogenase has encouraged synthetic attempts to prepare high-nuclearity M-Fe-S-X clusters containing such atoms. We have previously shown that reaction of the edge-bridged double cubane [(Tp)2Mo2Fe6S8(PEt3)4] (1) with nucleophiles HQ- affords the clusters [(Tp)2Mo2Fe6S8Q(QH)2](3-) (Q = S, Se) in which HQ- is a terminal ligand and Q(2-) is a mu2-bridging atom in the core. Reactions with OH- used as such or oxygen nucleophiles generated in acetonitrile from (Bu3Sn)2O or Me3SnOH and fluoride were examined. Reaction of 1 with Et4NOH in acetonitrile/water generates [(Tp)2Mo2Fe6S9(OH)2]3- (3), isolated as [(Tp)2Mo2Fe6S9(OH)(OC(=NH)Me)(H2O)](3-) and shown to have the [Mo2Fe6(mu2-S)2(mu3-S)6(mu6-S)] core topology very similar to the P(N) cluster of nitrogenase. The reaction system 1/Et4NOH in acetonitrile/methanol yields the P(N)-type cluster [(Tp)2Mo2Fe6S9(OMe)2(H2O)](3-) (5). The system 1/Me3SnOH/F- affords the oxo-bridged double P(N)-type cluster {[(Tp)2Mo2Fe6S9(mu2-O)]2}5- (7), convertible to the oxidized cluster {[(Tp)2Mo2Fe6S9(mu2-O)]2}4- (6), which is prepared independently from [(Tp)2Mo2Fe6S9F2(H2O)](3-)/(Bu3Sn)2O. In the preparations of 3-5 and 7, hydroxide liberates sulfide from 1 leading to the formation of P(N)-type clusters. Unlike reactions with HQ-, no oxygen atoms are integrated into the core structures of the products. However, the half-dimer composition [Mo2Fe6S9O] relates to the MoFe7S9 constitution of the putative native cluster with X = O. (Tp = hydrotris(pyrazolyl) borate(1-)).     

Interconversion and Reactivity of Two Heterometallic Tin-Containing Cuboidal Clusters from [Mo(3)S(4)(H(2)O)(9)](4+): X-ray Structure of the Single Cube with an Mo(3)SnS(4) Core

The Mo(3)SnS(4)(6+) single cube is obtained by direct addition of Sn(2+) to [Mo(3)S(4)(H(2)O)(9)](4+). UV-vis spectra of the product (0.13 mM) in 2.00 M HClO(4), Hpts, and HCl indicate a marked affinity of the Sn for Cl(-), with formation of the more strongly yellow [Mo(3)(SnCl(3))S(4)(H(2)O)(9)](3+) complex complete in as little as 0.050 M Cl(-). The X-ray crystal structure of (Me(2)NH(2))(6)[Mo(3)(SnCl(3))S(4)(NCS)(9)].0.5H(2)O has been determined and gives Mo-Mo (mean 2.730 ?) and Mo-Sn (mean 3.732 ?) distances, with a difference close to 1 ?. The red-purple double cube cation [Mo(6)SnS(8)(H(2)O)(18)](8+) is obtained by reacting Sn metal with [Mo(3)S(4)(H(2)O)(9)](4+). The double cube is also obtained in approximately 50% yield by BH(4)(-) reduction of a 1:1 mixture of [Mo(3)SnS(4)(H(2)O)(10)](6+) and [Mo(3)S(4)(H(2)O)(9)](4+). Conversely two-electron oxidation of [Mo(6)SnS(8)(H(2)O)(18)](8+) with [Co(dipic)(2)](-) or [Fe(H(2)O(6)](3+) gives the single cube [Mo(3)SnS(4)(H(2)O)(12)](6+) and [Mo(3)S(4)(H(2)O)(9)](4+) (up to 70% yield), followed by further two-electron oxidation to [Mo(3)S(4)(H(2)O)(9)](4+) and Sn(IV). The kinetics of the first stages have been studied using the stopped-flow method and give rate laws first order in [Mo(6)SnS(8)(H(2)O)(18)](8+) and the Co(III) or Fe(III) oxidant. The oxidation with [Co(dipic)(2)](-) has no [H(+)] dependence, [H(+)] = 0.50-2.00 M. With Fe(III) as oxidant, reaction steps involving [Fe(H(2)O)(6)](3+) and [Fe(H(2)O)(5)OH](2+) are implicated. At 25 degrees C and I = 2.00 M (Li(pts)) k(Co) is 14.9 M(-)(1) s(-)(1) and k(a) for the reaction of [Fe(H(2)O)(6)](3+) is 0.68 M(-)(1) s(-)(1) (both outer-sphere reactions). Reaction of Cu(2+) with the double but not the single cube is observed, yielding [Mo(3)CuS(4)(H(2)O)(10)](5+). A redox-controlled mechanism involving intermediate formation of Cu(+) and [Mo(3)S(4)(H(2)O)(9)](4+) accounts for the changes observed.     

Structural conversions of molybdenum-iron-sulfur edge-bridged double cubanes and P(n)-type clusters topologically related to the nitrogenase P-cluster

Edge-bridged Mo-Fe-S double cubanes are versatile precursors for the synthesis of other clusters of the same nuclearity. Thus, the double cubane [(Tp)(2)Mo(2)Fe(6)S(8)(PEt(3))(4)] sustains terminal ligand substitution with retention of the Mo(2)Fe(6)(micro(3)-S)(6)(micro(4)-S)(2) core structure and rearrangement to the Mo(2)Fe(6)(micro(2)-S)(2)(micro(3)-S)(6)(micro(6)-S) topology of the nitrogenase P(N) cluster upon reaction with certain nucleophiles. Four distinct processes for the conversion of double cubanes to P(N)-type clusters are documented, affording the products [(Tp)(2)Mo(2)Fe(6)S(9)(SR)(2)](3)(-), [(Tp)(2)Mo(2)Fe(6)S(8)(OMe)(3)](3)(-), and [(Tp)(2)Mo(2)Fe(6)S(7)(OMe)(4)](2)(-). In the latter clusters, two methoxides are terminal ligands and one or two are micro(2)-bridging ligands. The reverse transformation of a P(N)-type cluster to an edge-bridged double cubane has been demonstrated by the reaction of [(Tp)(2)Mo(2)Fe(6)S(8)(OMe)(3)](3)(-) with Me(3)SiX to afford [(Tp)(2)Mo(2)Fe(6)S(8)X(4)](2)(-) (X = Cl(-), Br(-)). Edge-bridged double cubanes have been obtained in the oxidation states [Mo(2)Fe(6)S(8)](2+,3+,4+). The stable oxidation state of P(N)-type clusters is [Mo(2)Fe(6)S(9)](+). Structures of five double cubanes and four P(N)-type clusters are reported. The P(N)-type clusters are synthetic representations of the biologically unique topology of the native P(N) cluster. Best-fit superpositions of the native and synthetic cluster cores gives weighted rms deviations in atom positions of 0.20-0.38 A. This study and an earlier investigation (Zhang, Y.; Holm, R. H. J. Am. Chem. Soc. 2003, 125, 3910-3920) provide a comprehensive account of the synthesis of structural analogues of the native P(N) cluster and provide the basis for continuing investigation of the synthesis of weak-field Mo-Fe-S clusters related to nitrogenase. (Tp = tris(pyrazolyl)hydroborate(1-).)     

Density functional theory calculations and exploration of a possible mechanism of N2 reduction by nitrogenase

Density functional theory (DFT) calculations have been performed on the nitrogenase cofactor, FeMoco. Issues that have been addressed concern the nature of M-M interactions and the identity and origin of the central light atom, revealed in a recent crystallographic study of the FeMo protein of nitrogenase (Einsle, O.; et al. Science 2002, 297, 871). Introduction of Se in place of the S atoms in the cofactor and energy minimization results in an optimized structure very similar to that in the native enzyme. The nearly identical, short, lengths of the Fe-Fe distances in the Se and S analogues are interpreted in terms of M-M weak bonding interactions. DFT calculations with O or N as the central atoms in the FeMoco marginally support the assignment of the central atom as N rather than O. The assumption was made that the central atom is the N atom, and steps of a catalytic cycle were calculated starting with either of two possible states for the cofactor and maintaining the same charge throughout (by addition of equal numbers of H(+) and e(-)) between steps. The states were [(Cl)Fe(II)(6)Fe(III)Mo(IV)S(9)(H(+))(3)N(3-)(Gl)(Im)](2-), [I-N-3H](2-), and [(Cl)Fe(II)(4)Fe(III)(3)Mo(IV)S(9)(H(+))(3)N(3-)(Gl)(Im)], [I-N-3H](0) (Gl = deprotonated glycol; Im = imidazole). These are the triply protonated ENDOR/ESEEM [I-N](5-) and M?ssbauer [I-N](3-) models, respectively. The proposed mechanism explores the possibilities that (a) redox-induced distortions facilitate insertion of N(2) and derivative substrates into the Fe(6) central unit of the cofactor, (b) the central atom in the cofactor is an exchangeable nitrogen, and (c) the individual steps are related by H(+)/e(-) additions (and reduction of substrate) or aquation/dehydration (and distortion of the Fe(6) center). The Delta E's associated with the individual steps of the proposed mechanism are small and either positive or negative. The largest positive Delta E is +121 kJ/mol. The largest negative Delta E of -333 kJ/mol is for the FeMoco with a N(3-) in the center (the isolated form) and an intermediate in the proposed mechanism.     

Ligand-bound S = 1/2 FeMo-cofactor of nitrogenase: hyperfine interaction analysis and implication for the central ligand X identity

Broken symmetry density functional theory (BS-DFT) has been used to address the hyperfine parameters of the single atom ligand X, proposed to be coordinated by six iron ions in the center of the paramagnetic FeMo-cofactor (FeMoco) of nitrogenase. Using the X = N alternative, we recently found that any hyperfine signal from X would be small (calculated A(iso)(X = (14)N) = 0.3 MHz) due to both structural and electronic symmetry properties of the [Mo-7Fe-9S- X] FeMoco core in its resting S = 3/2 state. Here, we extend our BS-DFT approach to the 2e(-) reduced S = 1/2 FeMoco state. Alternative substrates coordinated to this FeMoco state effectively perturb the electronic and/or structural symmetry properties of the cofactor. Using an example of an allyl alcohol (H2C=CH-CH2-OH) product ligand, we consider three different binding modes at single iron site and three different BS-DFT spin state structures and show that this binding would enhance the key hyperfine signal A(iso)(X) by at least 1 order of magnitude (3.8 < or = A(iso)(X = (14)N) < or = 14.7 MHz), and this result should not depend strongly on the exact identity of X (nitrogen, carbon, or oxygen). The interstitial atom, when the nucleus has a nonzero magnetic moment, should therefore be observable by ESR methods for some ligand-bound FeMoco states. In addition, our results illustrate structural details and likely spin-coupling patterns for models for early intermediates in the catalytic cycle.     

Electronic structure and spectroscopy of "superoxidized" iron centers in model systems: theoretical and experimental trends

Recent advances in synthetic chemistry have led to the discovery of "superoxidized" iron centers with valencies Fe(v) and Fe(vi) [K. Meyer et al., J. Am. Chem. Soc., 1999, 121, 4859-4876; J. F. Berry et al., Science, 2006, 312, 1937-1941; F. T. de Oliveira et al., Science, 2007, 315, 835-838.]. Furthermore, in recent years a number of high-valent Fe(iv) species have been found as reaction intermediates in metalloenzymes and have also been characterized in model systems [C. Krebs et al., Acc. Chem. Res., 2007, 40, 484-492; L. Que, Jr, Acc. Chem. Res., 2007, 40, 493-500.]. These species are almost invariably stabilized by a highly basic ligand X(n-) which is either O(2-) or N(3-). The differences in structure and bonding between oxo- and nitrido species as a function of oxidation state and their consequences on the observable spectroscopic properties have never been carefully assessed. Hence, fundamental differences between high-valent iron complexes having either Fe=O or Fe=N multiple bonds have been probed computationally in this work in a series of hypothetical trans-[FeO(NH(3))(4)OH](+/2+/3+) (1-3) and trans-[FeN(NH(3))(4)OH](0/+/2+) (4-6) complexes. All computational properties are permeated by the intrinsically more covalent character of the Fe=N multiple bond as compared to the Fe=O bond. This difference is likely due to differences in Z* between N and O that allow for better orbital overlap to occur in the case of the Fe=N multiple bond. Spin-state energetics were addressed using elaborate multireference ab initio computations that show that all species 1-6 have an intrinsic preference for the low-spin state, except in the case of 1 in which S=1 and S=2 states are very close in energy. In addition to M?ssbauer parameters, g-tensors, zero-field splitting and iron hyperfine couplings, X-ray absorption Fe K pre-edge spectra have been simulated using time-dependent DFT methods for the first time for a series of compounds spanning the high-valent states +4, +5, and +6 for iron. A remarkably good correlation of these simulated pre-edge features with experimental data on isolated high-valent intermediates has been found, allowing us to assign the main pre-edge features to excitations into the empty Fe d(z(2)) orbital, which is able to mix with Fe 4p(z), allowing an efficient mechanism for the intensification of pre-edge features.     

Mössbauer study of the three-coordinate planar Fe(II) thiolate complex [Fe(SR)(3)](-) (R = C(6)H(2)-2,4,6-tBu(3)): model for the trigonal iron sites of the MoFe(7)S(9):homocitrate cofactor of nitrogenase

The cofactor (M-center) of the MoFe protein of nitrogenase, a MoFe(7)S(9):homocitrate cluster, contains six Fe sites with a (distorted) trigonal sulfido coordination. These sites exhibit unusually small quadrupole splittings, Delta E(Q) approximately 0.7 mm/s, and isomer shifts, delta approximately 0.41 mm/s. M?ssbauer and ENDOR studies have provided the magnetic hyperfine tensors of all iron sites in the S = 3/2 state M(N). To assess the intrinsic zero-field splittings and hyperfine parameters of the cofactor sites, we have studied with M?ssbauer spectroscopy two salts of the three-coordinated Fe(II) thiolate complex [Fe(SR)(3)](-) (R = C(6)H(2)-2,4,6-tBu(3)). One of the salts, [Ph(4)P][Fe(SR)(3)] x 2MeCN x C(7)H(8), 1, has a planar geometry with idealized C(3h) symmetry. This S = 2 complex has an axial zero-field splitting with D = +10.2 cm(-1). The magnetic hyperfine tensor components A(x) = A(y) = -7.5 MHz and A(z) = -29.5 MHz reflect an orbital ground state with d(z(2)) symmetry. A(iso) = (A(x) +A(y) +A(z))/3 = -14.9 MHz, which includes the contact interaction (kappa P = -21.9 MHz) and an orbital contribution (+7 MHz), which is substantially smaller than A(iso) approximately -22 MHz of the tetrahedral Fe(II)(S-R)(4) sites of both rubredoxin and [PPh(4)](2)[Fe(II)(SPh)(4)]. The largest component of the electric field gradient (EFG) tensor is negative, as expected for a d(z(2)) orbital. However, Delta E(Q) = -0.83 mm/s, which is smaller than expected for a high-spin ferrous site. This reduction can be attributed to a ligand contribution, which in planar complexes provides a large positive EFG component perpendicular to the ligand plane. The isomer shift of 1, delta = 0.56 mm/s, approaches the delta-values reported for the six trigonal cofactor sites. The parameters of 1 and their importance for the cofactor cluster of nitrogenase are discussed.     

P(+) state of nitrogenase p-cluster exhibits electronic structure of a [fe(4)s(4)](+) cluster

Mo nitrogenase consists of two component proteins: the Fe protein, which contains a [Fe(4)S(4)] cluster, and the MoFe protein, which contains two different classes of metal cluster: P-cluster ([Fe(8)S(7)]) and FeMoco ([MoFe(7)S(9)C·homocitrate]). The P-cluster is believed to mediate the electron transfer between the Fe protein and the MoFe protein via interconversions between its various oxidation states, such as the all-ferrous state (P(N)) and the one- (P(+)) and two-electron (P(2+)) oxidized states. While the structural and electronic properties of P(N) and P(2+) states have been well characterized, little is known about the electronic structure of the P(+) state. Here, a mutant strain of Azotobacter vinelandii (DJ1193) was used to facilitate the characterization of the P(+) state of P-cluster. This strain expresses a MoFe protein variant (designated ΔnifB β-188(Cys) MoFe protein) that accumulates the P(+) form of P-cluster in the resting state. Magnetic circular dichroism (MCD) spectrum of the P-cluster in the oxidized ΔnifB β-188(Cys) MoFe protein closely resembles that of the P(2+) state in the oxidized wild-type MoFe protein, except for the absence of a major charge-transfer band centered at 823 nm. Moreover, magnetization curves of ΔnifB β-188(Cys) and wild-type MoFe proteins suggest that the P(2+) species in both proteins have the same spin state. MCD spectrum of the P(+) state in the ΔnifB β-188(Cys) MoFe protein, on the other hand, is associated with a classic [Fe(4)S(4)](+) cluster, suggesting that the P-cluster could be viewed as two coupled 4Fe clusters and that it could donate either one or two electrons to FeMoco by using one or both of its 4Fe halves. Such a mode of action of P-cluster could provide energetic and kinetic advantages to nitrogenase in the complex mechanism of N(2) reduction.     

DFT calculations of 57Fe Mössbauer isomer shifts and quadrupole splittings for iron complexes in polar dielectric media: applications to methane monooxygenase and ribonucleotide reductase

To predict the isomer shifts of Fe complexes in different oxidation and spin states more accurately, we have performed linear regression between the measured isomer shifts (delta(exp)) and DFT (PW91 potential with all-electron triple-zeta plus polarization basis sets) calculated electron densities at Fe nuclei [rho(0)] for the Fe(2+,2.5+) and Fe(2.5+,3+,3.5+,4+) complexes separately. The geometries and electronic structures of all complexes in the training sets are optimized within the conductor like screening (COSMO) solvation model. Based on the linear correlation equation delta(exp) = alpha[rho(0) - 11884.0] + C, the best fitting for 17 Fe(2+,2.5+) complexes (totally 31 Fe sites) yields alpha = -0.405 +/- 0.042 and C = 0.735 +/- 0.047 mm s(-1). The correlation coefficient is r = -0.876 with a standard deviation of SD = 0.075 mm s(-1). In contrast, the linear fitting for 19 Fe(2.5+,3+,3.5+,4+) complexes (totally 30 Fe sites) yields alpha = -0.393 +/- 0.030 and C = 0.435 +/- 0.014 mm s(-1), with the correlation coefficient r = -0.929 and a standard deviation SD = 0.077 mm s(-1). We provide a physical rationale for separating the Fe(2+,2.5+) fit from the Fe(2.5+,3+,3.5+,4+) fit, which also is clearly justified on a statistical empirical basis. Quadrupole splittings have also been calculated for these systems. The correlation between the calculated (DeltaE(Q(cal))) and experimental (DeltaE(Q(exp))) quadrupole splittings based on |DeltaE(Q(exp))| = A |DeltaE(Q(cal))| + B yields slope A, which is almost the ideal value 1.0 (A = 1.002 +/- 0.030) and intercept B almost zero (B = 0.033 +/- 0.068 mm s(-1)). Further calculations on the reduced diferrous and oxidized diferric active sites of class-I ribonucleotide reductase (RNR) and the hydroxylase component of methane monooxygenase (MMOH), and on a mixed-valent [(tpb)Fe3+(mu-O)(mu-CH3CO2)Fe4+(Me3[9]aneN3)]2+ (S = 3/2) complex and its corresponding diferric state have been performed. Calculated results are in very good agreement with the experimental data.     

ENDOR characterization of a synthetic diiron hydrazido complex as a model for nitrogenase intermediates

Molybdenum-dependent nitrogenase binds and reduces N2 at the [Fe7, Mo, S9, X, homocitrate] iron-molybdenum cofactor (FeMo-co). Kinetic and spectroscopic studies of nitrogenase variants indicate that a single Fe-S face is the most likely binding site. Recently, substantial progress has been made in determining the structures of nitrogenase intermediates formed during alkyne and N2 reduction through use of ENDOR spectroscopy. However, constraints derived from ENDOR studies of biomimetic complexes with known structure would powerfully contribute in turning experimentally derived ENDOR parameters into structures for species bound to FeMo-co during N2 reduction. The first report of a paramagnetic Fe-S compound that binds reduced forms of N2 involved Fe complexes stabilized by a bulky beta-diketiminate ligand (Vela, J.; Stoian, S.; Flaschenriem, C. J.; Münck, E.; Holland, P. L. J. Am. Chem. Soc. 2004, 126, 4522-4523). Treatment of a sulfidodiiron(II) complex with phenylhydrazine gave an isolable mixed-valence FeII-Fe(III) complex with a bridging phenylhydrazido (PhNNH2) ligand, and this species now has been characterized by ENDOR spectroscopy. Using both 15N, 2H labeled and unlabeled forms of the hydrazido ligand, the hyperfine and quadrupole parameters of the -N-NH2 moiety have been derived by a procedure that incorporates the (near-) mirror symmetry of the complex and involves a strategy which combines experiment with semiempirical and DFT computations. The results support the use of DFT computations in identifying nitrogenous species bound to FeMo-co of nitrogenase turnover intermediates and indicate that 14N quadrupole parameters from nitrogenase intermediates will provide a strong indication of the nature of the bound nitrogenous species. Comparison of the large 14N hyperfine couplings measured here with that of a hydrazine-derived species bound to FeMo-co of a trapped nitrogenase intermediate suggests that the ion(s) are not high spin and/or that the spin coupling coefficients of the coordinating cofactor iron ion(s) in the intermediate are exceptionally small.     

Theoretical investigation on the electronic and geometric structure of GaN2+ and GaN4+

The electronic and geometric structures of gallium dinitride cation, GaN2+ and gallium tetranitride cation, GaN4+ were systematically studied by employing density functional theory (DFT-B3LYP) and perturbation theory (MP2, MP4) in conjunction with large basis sets, (aug-)cc-pVxZ, x = T, Q. A total of 7 structures for GaN2+ and 24 for GaN4+ were identified, corresponding to minima, transition states, and saddle points. We report geometries and dissociation energies for all the above structures as well as potential energy profiles, potential energy surfaces, and bonding mechanisms for some low-lying electronic states. The calculated dissociation energy (De) of the ground state of GaN2+, X1Sigma+, is 5.6 kcal/mol with respect to Ga+(1S) + N2(X1Sigmag+) and that of the excited state, ?3Pi, is 24.8 kcal/mol with respect to Ga+(3P) + N2(X1Sigmag+). The ground state and the first excited minimum of GaN4+ are of 1A1(C2v) and 3B1(C2v) symmetry with corresponding De of 11.0 and 43.7 kcal/mol with respect to Ga+(1S) + 2N2(X1Sigmag+) for X1A1 and Ga+(3P) + 2N2(X1Sigmag+) for 3B1.     

NFeN结构和性质及反应势能曲线

用B3LYP/6-311+G(d)方法对化合物NFeN弯曲型和直线型的不同自旋多重度多个电子态的几何结构、电子结构、能量和振动光谱进行了计算研究. 结果表明, 单重态中Fe―N键长普遍比三重态和五重态中的短, 在155 pm左右; NFeN三重态电子结构最丰富, 自然键轨道和Mulliken布居显示Fe―N键具有部分离子键特征; 两种结构的所有稳定态中能量最低的是1⁵A₂态, 能量相近的有1³B₁、1³A₂、1³B₂和1¹A₁态, 直线型中能量最低是³Δ_g态; 相对于分子基态反应物Fe(a⁵D)+N₂(X¹Σ_g⁺)所有电子态的能量都偏高, 该反应在热力学上是不利的, 但是对于原子态反应物Fe(a⁵D)+2N(⁴S)则是放热反应; 计算振动频率和强度与实验较吻合的是1³B₁态; 复合物FeN₂与化合物NFeN结构差异明显; Fe原子直接插入N₂分子的势能曲线表明该反应能垒很高, 在动力学上也是不利的.