Exploring the interstitial atom in the FeMo cofactor of nitrogenase: insights from QM and QM/MM calculations

Density functional theory and combined quantum mechanics and molecular mechanics (QM/MM) calculations have been used to explore structural features of the FeMo cofactor with an interstitial atom X (X = N, C, or O) and its interactions with CO and N 2. Predicted frequencies of the metal-bound CO, QM/MM-optimized geometries, and calculated redox potentials of the FeMo cofactor with different central ligands show that the oxygen atom is the candidate for the interstitial atom. Calculations on the interactions of the FeMo cofactor with CO and N 2 reveal that there is a remarkable dependence of the binding energy on the binding site and the interstitial atom. Generally, the Fe2 site of the FeMo cofactor has stronger interactions with CO and N 2 than Fe6, and both the Fe2 and Fe6 sites in the N-centered and O-centered clusters of the FeMo cofactor can effectively bind N 2 while the coordination of N 2 to the Fe6 site of the C-centered active cluster is unfavorable energetically. Present results indicate that the protein environment is important for computational characterization of the structure of the FeMo cofactor and properties of the metal-bound CO and N 2 are sensitive to the interstitial atom.     

Chemical activity of the nitrogenase FeMo cofactor with a central nitrogen ligand: density functional study

We investigate the chemical consequences of a central ligand in the nitrogenase FeMo cofactor using density functional calculations. Several studies have shown that the central ligand most probably is a nitrogen atom, but the consequences for the chemical reactivity of the cofactor are unknown. We investigate several possible routes for insertion of the central nitrogen ligand and conclude that all routes involve barriers and intermediate states, which are inaccessible at ambient conditions. On this basis we suggest that the central nitrogen ligand is present at all times during the reaction. Furthermore, we investigate how the FeMoco with the central ligand can interact with N(2) and reduce it.     

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.     

Enzymatic and catalytic reduction of dinitrogen to ammonia: Density functional theory characterization of alternative molybdenum active sites

We used density functional calculations to model dinitrogen reduction by a FeMo cofactor containing a central nitrogen atom and by a Mo‐based catalyst. Plausible intermediates, reaction pathways, and relative energetics in the enzymatic and catalytic reduction of N₂ to ammonia at a single Mo center are explored. Calculations indicate that the binding of N₂ to the Mo atom and the subsequent multiple proton–electron transfer to dinitrogen and its protonated species involved in the conversion of N₂ are feasible energetically. In the reduction of N₂ the Mo atom experiences a cycled oxidation state from Mo(IV) to Mo(VI) by nitrogenase and from Mo(III) to Mo(VI) by the molybdenum catalyst, respectively, tuning the gradual reduction of N₂. Such a wide range of oxidation states exhibited by the Mo center is crucial for the gradual reduction process via successive proton–electron transfer. Present results suggest that the Mo atom in the N‐centered FeMo cofactor is a likely alternative active site for dinitrogen binding and reduction under mild conditions once there is an empty site available at the Mo site. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

固氮酶的固氮机理和其人工模拟问题的探讨

固氮酶将N₂ 还原为NH₃ 的过程是自然界实现氮循环的重要环节。本文着重对固氮酶的固氮机理和其活性中心FeMo 辅基的人工模拟合成进行探讨, 其中包括FeMo蛋白中的质子和电子的传递, FeMo 辅基对N₂ 的活化方式,Mo 原子的作用, 固氮活性的测试。最后还就固氮酶的活性中心FeMo 辅基的人工模拟合成进行了探讨。     

The Interstitial Carbon of the Nitrogenase FeMo Cofactor is Far Better Stabilized than Previously Assumed

The first quantum-mechanical calculations of all relevant potential constants in both the iron-molybdenum cofactor and the iron-vanadium cofactor of nitrogenase suggest that the carbide is bound to the center of the enzyme much more strongly than hitherto assumed. Previous studies seemed to indicate a dummy function of the interstitial carbon, with a weak force constant (ca. 0.32 N cm⁻¹). Our new investigations confirm a different picture: the central carbon atom binds the iron-sulfur cluster through six covalent C−Fe bonds. With a potential constant of more than 1.3 N cm⁻¹, the interstitial carbon also appears to be dynamically persistent. According to our investigations, the values for the elasticity within the iron-sulfur cluster have to be corrected too. These new details on the mechano-chemical properties of the FeMo cofactor will be important for elucidating the catalytic cycle of nitrogen fixation. By implementing our new algorithm in the freely available COMPLIANCE program, the dependence on the coordinates during the calculation of Hesse matrices is eliminated completely.     

Modeling the nitrogenase FeMo cofactor with high‐spin Fe8S9X+ (XN,C) clusters. Is the first step for N2 reduction to NH3 a concerted dihydrogen transfer?

A high-spin Fe(8)S(9)X(+) (X=N, C) cluster is used to model the reduction of molecular nitrogen to ammonia by the nitrogenase FeMo cofactor at the B3LYP/6-311G(d,p)/ECP(Fe,SDD) level of theory. A total of seventy-three structures were optimized (including three transition state optimizations) to explore the structure and energetic of N(2), C(2)H(2), and CO coordination to the Fe(8)S(9)X(+) cluster. After three protonation-reduction (PR) steps (modeled by addition of hydrogen atoms), N(2), C(2)H(2), and CO are predicted to bind to a Fe atom in the exo (cage does not open) position with binding energies of 7.6, 14.7, and 11.7 kcal/mol. With additional PR steps the coordination number of the core nitrogen atom is reduced from six to five and the bridging thiol group becomes a terminal SH(2) group. The fifth and sixth PR steps occur on the core nitrogen and the open Fe site. Coordination of N(2) is enhanced after six PR steps to give an intermediate ideally suited for a concerted dihydrogen transfer from the Fe and core nitrogen atoms to the coordinated N(2). The identity of the central atom (nitrogen or carbon) has only a minor effect on the reaction steps.     

Structural,spectroscopic, and redox consequences of a central ligand in the FeMoco of nitrogenase: a density functional theoretical study

Broken symmetry density functional and electrostatics calculations have been used to shed light on which of three proposed atoms, C, N, or O, is most likely to be present in the center of the FeMoco, the active site of nitrogenase. At the Mo(4+)4Fe(2+)3Fe(3+) oxidation level, a central N(3-) anion results in (1) calculated Fe-N bond distances that are in very good agreement with the recent high-resolution X-ray data of Einsle et al.; (2) a calculated redox potential of 0.19 eV versus the standard hydrogen electrode (SHE) for FeMoco(oxidized) + e(-) --> FeMoco(resting), in good agreement with the measured value of -0.042 V in Azotobacter vinelandii; and (3) average M?ssbauer isomer shift values (IS(av) = 0.48 mm s(-1)) compatible with experiment (IS(av) = 0.40 mm s(-1)). At the more reduced Mo(4+)6Fe(2+)1Fe(3+) level, the calculated geometry around a central N(3-) anion still correlates well with the X-ray data, but the average M?ssbauer isomer shift value (IS(av) = 0.54 mm s(-1)) and the redox potential of -2.21 eV show a much poorer agreement with experiment. These calculated structural, spectroscopic, and redox data indicate the most likely iron oxidation state for the resting FeMoco of nitrogenase to be 4Fe(2+)3Fe(3+). At this favored oxidation state, oxygen or carbon coordination leads to (1) Fe-O distances in poor agreement and Fe-C distances in good agreement with experiment and (2) calculated redox potentials of +0.97 eV for O(2-) and -1.31 eV for C(4-). The calculated structural parameters and/or redox data suggest either O(2-) or C(4-) is unlikely as a central anion.     

Intermediates trapped during nitrogenase reduction of N triple bond N, CH3-N=NH, and H2N-NH2

A high population intermediate has been trapped on the nitrogenase active site FeMo cofactor during reduction of N2. In addition, intermediates have been trapped during reduction of CH3-N=NH by the alpha-195Gln variant and during reduction of H2N-NH2 by the alpha-70Ala/alpha-195Gln variant. Each of these trapped states shows an EPR signal arising from an S = 1/2 state of the FeMo cofactor. 15N ENDOR shows that each intermediate has a nitrogenous species bound to the FeMo cofactor, with a single type of N seen for each bound intermediate. The g tensors are unique to each intermediate, g(e) = [2.084, 1.993, 1.969], g(m) = [2.083, 2.021, 1.993], g(l) = [2.082, 2.015, 1.987], as are the 15N hyperfine couplings at g1, which suggests that three distinct stages of NN reduction may have been trapped. The 1H ENDOR spectra show that the N2 intermediate is at a distinct and earlier stage of reduction from the other two, so at least two stages of NN reduction have been trapped. Some possible structures of the hydrazine intermediate are presented.     

Ammonia production at the FeMo cofactor of nitrogenase: results from density functional theory

Biological nitrogen fixation has been investigated beginning with the monoprotonated dinitrogen bound to the FeMo cofactor of nitrogenase up to the formation of the two ammonia molecules. The energy differences of the relevant intermediates, the reaction barriers, and potentially relevant side branches are presented. During the catalytic conversion, nitrogen bridges two Fe atoms of the central cage, replacing a sulfur bridge present before dinitrogen binds to the cofactor. A transformation from cis- to trans-diazene has been found. The strongly exothermic cleavage of the dinitrogen bond takes place, while the Fe atoms are bridged by a single nitrogen atom. The dissociation of the second ammonia from the cofactor is facilitated by the closing of the sulfur bridge following an intramolecular proton transfer. This closes the catalytic cycle.     

Model for acetylene reduction by nitrogenase derived from density functional theory

The catalytic cycle of acetylene reduction at the FeMo cofactor of nitrogenase has been investigated on the basis of density functional theory. C2H2 binds to the same site as N2, but it binds to a less reduced state of the cofactor. In a manner similar to that of N2 binding, one of the sulfur bridges opens during acetylene binding. The model explains the strong noncompetitive inhibition of N2 reduction by C2H2 and the weak competitive inhibition of C2H2 reduction by N2. Our proposed mechanism is consistent with experimentally observed stereoselectivity and the ability of C2H2 to suppress H2 production by nitrogenase.     

Studies on the hydrogenation steps of the nitrogen molecule at the Azotobacter vinelandii nitrogenase site

We follow the initial activation of the nitrogen molecule at the FeMo cofactor of nitrogenase and subsequently model the hydrogenation of N₂ up to the fourth protonation step using the intermediate neglect of differential overlap quantum-chemical model. The results obtained favor a reaction mechanism going through hydrazido intermediates on the 4-Fe surfaces, externally to the FeMo cofactor. Calculations using density functional theory on smaller model systems also support the suggested mechanism over other possible schemes that involve early release of the first molecule of ammonia as a product of the enzymatic reaction. We also demonstrate that dielectric stabilization due to the protein around the cofactor could lower markedly the barrier for the product release as an ammonium ion. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 70: 1159–1168, 1998     

The Nitrogenase FeMo‐Cofactor Precursor Formed by NifB Protein: A Diamagnetic Cluster Containing Eight Iron Atoms

The biological activation of N₂ occurs at the FeMo‐cofactor, a 7Fe–9S–Mo–C–homocitrate cluster. FeMo‐cofactor formation involves assembly of a Fe_6–8–S_X–C core precursor, NifB‐co, which occurs on the NifB protein. Characterization of NifB‐co in NifB is complicated by the dynamic nature of the assembly process and the presence of a permanent [4Fe–4S] cluster associated with the radical SAM chemistry for generating the central carbide. We have used the physiological carrier protein, NifX, which has been proposed to bind NifB‐co and deliver it to the NifEN protein, upon which FeMo‐cofactor assembly is ultimately completed. Preparation of NifX in a fully NifB‐co‐loaded form provided an opportunity for Mössbauer analysis of NifB‐co. The results indicate that NifB‐co is a diamagnetic (S=0) 8‐Fe cluster, containing two spectroscopically distinct Fe sites that appear in a 3:1 ratio. DFT analysis of the ⁵⁷Fe electric hyperfine interactions deduced from the Mössbauer analysis suggests that NifB‐co is either a 4Fe²⁺–4Fe³⁺ or 6Fe²⁺–2Fe³⁺ cluster having valence‐delocalized states.     

Role of the pyridine nitrogen in pyridoxal 5'-phosphate catalysis: activity of three classes of PLP enzymes reconstituted with deazapyridoxal 5'-phosphate

Pyridoxal 5'-phosphate (PLP; vitamin B(6))-catalyzed reactions have been well studied, both on enzymes and in solution, due to the variety of important reactions this cofactor catalyzes in nitrogen metabolism. Three functional groups are central to PLP catalysis: the C4' aldehyde, the O3' phenol, and the N1 pyridine nitrogen. In the literature, the pyridine nitrogen has traditionally been assumed to be protonated in enzyme active sites, with the protonated pyridine ring providing resonance stabilization of carbanionic intermediates. This assumption is certainly correct for some PLP enzymes, but the structures of other active sites are incompatible with protonation of N1, and, consequently, these enzymes are expected to use PLP in the N1-unprotonated form. For example, aspartate aminotransferase protonates the pyridine nitrogen for catalysis of transamination, while both alanine racemase and O-acetylserine sulfhydrylase are expected to maintain N1 in the unprotonated, formally neutral state for catalysis of racemization and β-elimination. Herein, kinetic results for these three enzymes reconstituted with 1-deazapyridoxal 5'-phosphate, an isosteric analogue of PLP lacking the pyridine nitrogen, are compared to those for the PLP enzyme forms. They demonstrate that the pyridine nitrogen is vital to the 1,3-prototropic shift central to transamination, but not to reactions catalyzed by alanine racemase or O-acetylserine sulfhydrylase. Not all PLP enzymes require the electrophilicity of a protonated pyridine ring to enable formation of carbanionic intermediates. It is proposed that modulation of cofactor electrophilicity plays a central role in controlling reaction specificity in PLP enzymes.     

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.     

Double sandwich polyoxometalate and its Fe(III) substituted derivative, [As2Fe5Mo21O82]17– and [As2Fe6Mo20O80(H2O)2]16–

A double sandwich polyoxometalate and its Fe(III) substituted derivative, [As(2)Fe(5)Mo(21)O(82)](17-) (1) and [As(2)Fe(6)Mo(20)O(80)(H(2)O)(2)](16-) (2), were synthesized and characterized by single-crystal X-ray diffraction, infrared spectroscopy, fluorescent spectroscopy, UV spectra, thermogravimetry-differential scanning calorimetry analyses, electrospray ionization mass spectrometry, and magnetism measurements. The polyoxoanion is composed of a central fragment FeMo(7)O(28) for 1 (Fe(2)Mo(6)O(26)(H(2)O)(2) for 2) and two external AsMo(7)O(27) fragments linked together by two distinct edge-sharing dimeric clusters Fe(2)O(10) to lead to a C(2v) molecular symmetry. The central FeMo(7)O(28) fragment and external AsMo(7)O(27) fragment have a similar structure, and both of them can be viewed as a monocapped hexavacant α-Keggin subunit with a central FeO(4) group or a central AsO(3) group. Both of the polyoxoanions contain a oxo-bridged Fe(III)(5) magnetic core with the angles of Fe-O-Fe in the range of 96.4(4)-125.7(5)°, and magnetism measurements show an overall ferromagnetic interactions among the five-nuclearity cluster Fe(5) with the spin ground state S = 15/2.     

Die Reaktionen der Benzoyl-und Salicoylhydrazone des Vanillins,Furfurols und Zimtaldehyds mit zweiwertigen Ionen der Übergangsmetalle

The reaction of benzoyl and salicoylhydrazones of vanilline, furfural, and cinnamaldehyde with some divalent metal ions of the first transition series is investigated. The structure of the ligand in the solid chelates is studied by IR-spectrophotometry which showed that the ligand is coordinated to the central metal ion through the oxygen atom of the C=O and nitrogen atom of the C=N group. Spectrophotometric, and conductometric studies were carried out to investigate the stoichiometry of the complexes under consideration. The shifts in the C=O and C=N IR-bands are utilised for the determination of the coordination bond length.

---

---

Mit 6 Abbildungen     

Steric effects on uranyl complexation: synthetic, structural, and theoretical studies of carbamoyl pyrazole compounds of the uranyl(VI) ion

New bifunctional pyrazole based ligands of the type [C(3)HR(2)N(2)CONR'] (where R = H or CH(3); R' = CH(3), C(2)H(5), or (i)C(3)H(7)) were prepared and characterized. The coordination chemistry of these ligands with uranyl nitrate and uranyl bis(dibenzoyl methanate) was studied with infrared (IR), (1)H NMR, electrospray-mass spectrometry (ES-MS), elemental analysis, and single crystal X-ray diffraction methods. The structure of compound [UO(2)(NO(3))(2)(C(3)H(3)N(2)CON{C(2)H(5)}(2))] (2) shows that the uranium(VI) ion is surrounded by one nitrogen atom and seven oxygen atoms in a hexagonal bipyramidal geometry with the ligand acting as a bidentate chelating ligand and bonds through both the carbamoyl oxygen and pyrazolyl nitrogen atoms. In the structure of [UO(2)(NO(3))(2)(H(2)O)(2)(C(5)H(7)N(2)CON {C(2)H(5)}(2))(2)], (5) the pyrazole ligand acts as a second sphere ligand and hydrogen bonds to the water molecules through carbamoyl oxygen and pyrazolyl nitrogen atoms. The structure of [UO(2)(DBM)(2)C(3)H(3)N(2)CON{C(2)H(5)}(2)] (8) (where DBM = C(6)H(5)COCHCOC(6)H(5)) shows that the pyrazole ligand acts as a monodentate ligand and bonds through the carbamoyl oxygen to the uranyl group. The ES-MS spectra of 2 and 8 show that the ligand is similarly bonded to the metal ion in solution. Ab initio quantum chemical studies show that the steric effect plays the key role in complexation behavior.     

Theoretical Study on the Dissociation of Ligands in the Rhodium and Iridium Complexes Containing 1,1,1,5,5,5-Hexafluoroacetylacetonato

Functionalization of the inert C? H bonds of unsaturated molecules by transition metal complex is an important means to form new C? C bonds. The functionalization is usually initiated by the ligand dissociation of a complex. In this paper we employ both ab initio and density functional methods to explore the influence of central metals, conformation, solvent and protonation on the ligand dissociation of the (hfac‐O,O)₂M(L)(py) complexes [M=Rh(III) or Ir(III), hfac‐O,O=k²‐O,O‐1,1,1,5,5,5‐hexafluoroacetylacetonato, L=CH₃, CH₃CO₂, (CH₃CO)₂CH, CH₃O or OH, py=pyridine]. We demonstrate that ligand pyridine dissociates more easily than the "L" ligands under study in aprotic solvent and gas phase and the dissociation of pyridine is more facile in the trans‐conformation than in the cis‐isomer. These phenomena are rationalized based on electronic structure and molecular orbital interactions. We show that solvation only slightly stabilizes the complexes and does not change the ligand dissociation ordering. In particular, we show that pyridine is no longer the labile ligand in protic media. Instead, the oxygen‐containing ligands (apart from those like hfac that form a cyclic structure with the central metal) that coordinate to the central metal via oxygen atom become the labile ones. Finally our calculations indicate that hfac is a stable ligand, even in protic media.     

Controlled Incorporation of Nitrides into W‐Fe‐S Clusters

Incorporation of monatomic 2p ligands into the core of iron–sulfur clusters has been researched since the discovery of interstitial carbide in the FeMo cofactor of Mo‐dependent nitrogenase, but has proven to be a synthetic challenge. Herein, two distinct synthetic pathways are rationalized to install nitride ligands into targeted positions of W‐Fe‐S clusters, generating unprecedented nitride‐ligated iron–sulfur clusters, namely [(Tp*)₂W₂Fe₆(μ₄‐N)₂S₆L₄]^2? (Tp*=tris(3,5‐dimethyl‐1‐pyrazolyl)hydroborate(1?), L=Cl^? or Br^?). ⁵⁷Fe Mössbauer study discloses metal oxidation states of W^IV₂Fe^II₄Fe^III₂ with localized electron distribution, which is analogous to the mid‐valent iron centres of FeMo cofactor at resting state. Good agreement of Mössbauer data with the empirical linear relationship for Fe–S clusters indicates similar ligand behaviour of nitride and sulfide in such clusters, providing useful reference for reduced nitrogen in a nitrogenase‐like environment.