Carbon Monoxide Binding to the Iron-Molybdenum Cofactor of Nitrogenase: a Detailed Quantum Mechanics/Molecular Mechanics Investigation

Highly Activated Terminal Carbon Monoxide Ligand in an Iron-Sulfur Cluster Model of FeMco with Intermediate Local Spin State at Fe

Nitrogenases, the enzymes that convert N2 to NH3, also catalyze the reductive coupling of CO to yield hydrocarbons. CO-coordinated species of nitrogenase clusters have been isolated and used to infer mechanistic information. However, synthetic FeS clusters displaying CO ligands remain rare, which limits benchmarking. Starting from a synthetic cluster that models a cubane portion of the FeMo cofactor (FeMoco), including a bridging carbyne ligand, we report a heterometallic tungsten-iron-sulfur cluster with a single terminal CO coordination in two oxidation states with a high level of CO activation (νCO = 1851 and 1751 cm-1). The local Fe coordination environment (2S, 1C, 1CO) is identical to that in the protein making this system a suitable benchmark. Computational studies find an unusual intermediate spin electronic configuration at the Fe sites promoted by the presence the carbyne ligand. This electronic feature is partly responsible for the high degree of CO activation in the reduced cluster.

Mechanism of Nitrogen Reduction to Ammonia in a Diiron Model of Nitrogenase

Nitrogenase is a fascinating enzyme in biology that reduces dinitrogen from air to ammonia through stepwise reduction and protonation. Despite it being studied in detail by experimental and computational groups, there are still many unknown factors in the catalytic cycle of nitrogenase, especially related to the addition of protons and electrons and their order. A recent biomimetic study characterized a potential dinitrogen-bridged diiron cluster as a synthetic model of nitrogenase. Using strong acid and reductants, the dinitrogen was converted into ammonia molecules, but details of the mechanism remains unknown. In particular, it was unclear from the experimental studies whether the proton and electron transfer steps are sequential or alternating. Moreover, the work failed to establish what the function of the diiron core is and whether it split into mononuclear iron fragments during the reaction. To understand the structure and reactivity of the biomimetic dinitrogen-bridged diiron complex [(P2P'PhFeH)2(μ-N2)] with triphenylphosphine ligands, we performed a density functional theory study. Our computational methods were validated against experimental crystal structure coordinates, Mössbauer parameters, and vibrational frequencies and show excellent agreement. Subsequently, we investigated the alternating and consecutive addition of electrons and protons to the system. The calculations identify a number of possible reaction channels, namely, same-site protonation, alternating protonation, and complex dissociation into mononuclear iron centers. The calculations show that the overall mechanism is not a pure sequential set of electron and proton transfers but a mixture of alternating and consecutive steps. In particular, the first reaction steps will start with double proton transfer followed by an electron transfer, while thereafter, there is another proton transfer and a second electron transfer to give a complex whereby ammonia can split off with a low energetic barrier. The second channel starts with alternating protonation of the two nitrogen atoms, whereafter the initial double proton transfer, electrons and protons are added sequentially to form a hydrazine-bound complex. The latter split off ammonia spontaneously after further protonation. The various reaction channels are analyzed with valence bond and orbital diagrams. We anticipate the nitrogenase enzyme to operate with mixed alternating and consecutive protonation and electron transfer steps.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

The E3 state of FeMoco: one hydride, two hydrides or dihydrogen?

Hydrides are present in the reduced states of the iron-molybdenum cofactor (FeMoco) of Mo nitrogenase and are believed to play a key mechanistic role in the dinitrogen reduction reaction catalyzed by the enzyme. Two hydrides are present in the E4 state according to 1H ENDOR and there is likely a single hydride in the E2 redox state. The 2-hydride E4 state has been experimentally observed to bind N2 and it has been speculated that E3 may bind N2 as well. However, the E3 state has not been directly observed and very little is known about its molecular and electronic structure or reactivity. In recent computational studies, we have explored the energy surfaces of the E2 and E4 by QM/MM modelling, and found that the most stable hydride isomers contain bridging or partially bridging hydrides with an open protonated belt sulfide-bridge. In this work we systematically explore the energy surface of the E3 redox state, comparing single hydride and two-hydride isomers with varying coordination and bridging vs. terminal sulfhydryl groups. We also include a model featuring a triply protonated carbide. The results are only mildly dependent on the QM-region size. The three most stable E3 isomers at the r2SCAN level of theory have in common: an open belt sulfide-bridge (terminal sulfhydryl group on Fe6) and either 2 bridging hydrides (between Fe2 and Fe6), 1 bridging-1-terminal hydride (around Fe2 and Fe6) or a dihydrogen ligand bound at the Fe2 site. Analyzing the functional dependency of the results, we find that functionals previously found to predict accurate structures of spin-coupled Fe/Mo dimers and FeMoco (TPSSh, B97-D3, r2SCAN, and B3LYP*) are in generally good agreement about the stability of these 3 E3 isomers. However, B3LYP*, similar to its parent B3LYP method, predicts a triply protonated carbide isomer as the most stable isomer, an unlikely scenario in view of the lack of experimental evidence for carbide protonation occurring in reduced FeMoco states. Distinguishing further between the 3 hydride isomers is difficult and this flexible coordination nature of hydrides suggests that multiple hydride isomers could be present during experimental conditions. N2 binding was explored and resulted in geometries with 2 bridging hydrides and N2 bound to either Fe2 or Fe6 with a local low-spin state on the Fe. N2 binding is predicted to be mildly endothermic, similar to the E2 state, and it seems unlikely that the E3 state is capable of binding N2.

Influence of the ligand-field on EPR parameters of cis- and trans-isomers in MoV systems relevant to molybdenum enzymes: Experimental and density functional theory study

The electron paramagnetic resonance (EPR) investigation of mononuclear cis- and trans-(L1O)MoOCl2 complexes [L1OH = bis(3,5-dimethylpyrazolyl)-3-tert-butyl-2-hydroxy-5-methylphenyl)methane] reveals a significant difference in their spin Hamiltonian parameters which reflect different equatorial and axial ligand fields created by the heteroscorpionate donor atoms. Density functional theory (DFT) was used to calculate the values of principal components and relative orientations of the g and A tensors, and the molecular framework in four pairs of isomeric mononuclear oxo‑molybdenum(V) complexes (cis- and trans-(L1O)MoOCl2, cis,cis- and cis,trans-(L-N2S2)MoOCl [L-N2S2H2 = N,N'-dimethyl-N,N'-bis(mercaptophenyl)ethylenediamine], cis,cis- and cis,trans-(L-N2S2)MoO(SCN), and cis- and trans-[(dt)2MoO(OMe)]2- [dtH2 = 2,3-dimercapto-2-butene]). Scalar relativistic DFT calculations were conducted using three different exchange-correlation functionals. It was found that the use of hybrid exchange-correlation functional with 25% of the Hartree-Fock exchange leads to the best quantitative agreement between theory and experiment. A simplified ligand-field approach was used to analyze the influence of the ligand fields in all cis- and trans-isomers on energies and contributions of molybdenum d-orbital manifold to g and A tensors and relative orientations. Specifically, contributions that originated from the spin-orbit coupling of the dxz, dyz, and dx2-y2 orbitals into the ground state have been discussed. The new findings are discussed in the context of the experimental data of mononuclear molybdoenzyme, DMSO reductase.

The HD Reaction of Nitrogenase: a Detailed Mechanism

Nitrogenase is the enzyme that converts N₂ to NH₃ under ambient conditions. The chemical mechanism of this catalysis at the active site FeMo-co [Fe₇ S₉ CMo(homocitrate)] is unknown. An obligatory co-product is H₂ , while exogenous H₂ is a competitive inhibitor. Isotopic substitution using exogenous D₂ revealed the N₂ -dependent reaction D₂ +2H⁺ +2e^- →2HD (the 'HD reaction'), together with a collection of additional experimental characteristics and requirements. This paper describes a detailed mechanism for the HD reaction, developed and elaborated using density functional simulations with a 486-atom model of the active site and surrounding protein. First D₂ binds at one Fe atom (endo-Fe6 coordination position), where it is flanked by H-Fe6 (exo position) and H-Fe2 (endo position). Then there is synchronous transfer of these two H atoms to bound D₂ , forming one HD bound to Fe2 and a second HD bound to Fe6. These two HD dissociate sequentially. The final phase is recovery of the two flanking H atoms. These H atoms are generated, sequentially, by translocation of a proton from the protein surface to S3B of FeMo-co and combination with introduced electrons. The first H atom migrates from S3B to exo-Fe6 and the second from S3B to endo-Fe2. Reaction energies and kinetic barriers are reported for all steps. This mechanism accounts for the experimental data: (a) stoichiometry; (b) the N₂ -dependence results from promotional N₂ bound at exo-Fe2; (c) different N₂ binding K_m for the HD reaction and the NH₃ formation reaction results from involvement of two different sites; (d) inhibition by CO; (e) the non-occurrence of 2HD→H₂ +D₂ results from the synchronicity of the two transfers of H to D₂ ; (f) inhibition of HD production at high pN₂ is by competitive binding of N₂ at endo-Fe6; (g) the non-leakage of D to solvent follows from the hydrophobic environment and irreversibility of proton introduction.

Understanding the tethered unhooking and rehooking of S2B in the reaction domain of FeMo-co, the active site of nitrogenase

The active site of the nitrogen fixing enzyme nitrogenase is an Fe₇MoS₉C cluster, and investigations of the enigmatic chemical mechanism of the enzyme have focussed on a pair of Fe atoms, Fe2 and Fe6, and the S2B atom that bridges them. There are three proposals for the status of the Fe2-S2B-Fe6 bridge during the catalytic cycle: one that it remains intact, another that it is completely labile and absent during catalysis, and a third that S2B is hemilabile, unhooking one of its bonds to Fe2 or Fe6. This report examines the tethered unhooking of S2B and factors that affect it, using DFT calculations of 50 geometric/electronic possibilities with a 485 atom model including all relevant parts of surrounding protein. The outcomes are: (a) unhooking the S2B-Fe2 bond is feasible and favourable, but alternative unhooking of the S2B-Fe6 bond is unlikely for steric reasons, (b) energy differences between hooked and unhooked isomers are generally <10 kcal mol^-1, usually with unhooked more stable, (c) ligation at the exo-Fe6 position inhibits unhooking, (d) unhooking of hydrogenated S2B is more favourable than that of bare S2B, (e) hydrogen bonding from the NεH function of His195 to S2B occurs in hooked and unhooked forms, and possibly stabilises unhooking, (f) unhooking is reversible with kinetic barriers ranging 10-13 kcal mol^-1. The conclusion is that energetically accessible reversible unhooking of S2B or S2BH, as an intrinsic property of FeMo-co, needs to be considered in the formulation of mechanisms for the reactions of nitrogenase.

Nitrogenase Chemistry at 10 Kelvin─Phototautomerization and Recombination of CO-Inhibited α-H195Q Enzyme

CO-bound forms of nitrogenase are N₂-reduction inhibited and likely intermediates in Fischer-Tropsch chemistry. Visible-light photolysis at 7 K was used to interrogate all three known CO-related EPR-active forms as exhibited by the α-H195Q variant of Azotobacter vinelandii nitrogenase MoFe protein. The hi(5)-CO EPR signal converted to the hi-CO EPR signal, which reverted at 10 K. FT-IR monitoring revealed an exquisitely light-sensitive "Hi-2" species with bands at 1932 and 1866 cm^-1 that yielded "Hi-1" with bands at 1969 and 1692 cm^-1, which reverted to "Hi-2". The similarities of photochemical behavior and recombination kinetics showed, for the first time, that hi-CO EPR and "Hi-1" IR signals arise from one chemical species. hi(5)-CO EPR and "Hi-2" IR signals are from a second species, and lo-CO EPR and "Lo-2" IR signals, formed after prolonged illumination, are from a third species. Comparing FT-IR data with CO-inhibited MoFe-protein crystal structures allowed assignment of CO-bonding geometries in these species.

Pulsed Multifrequency Electron Paramagnetic Resonance Spectroscopy Reveals Key Branch Points for One- vs Two-Electron Reactivity in Mn/Fe Proteins

Traditionally, the ferritin-like superfamily of proteins was thought to exclusively use a diiron active site in catalyzing a diverse array of oxygen-dependent reactions. In recent years, novel redox-active cofactors featuring heterobimetallic Mn/Fe active sites have been discovered in both the radical-generating R2 subunit of class Ic (R2c) ribonucleotide reductases (RNRs) and the related R2-like ligand-binding oxidases (R2lox). However, the protein-specific factors that differentiate the radical reactivity of R2c from the C-H activation reactions of R2lox remain unknown. In this work, multifrequency pulsed electron paramagnetic resonance (EPR) spectroscopy and ligand hyperfine techniques in conjunction with broken-symmetry density functional theory calculations are used to characterize the molecular and electronic structures of two EPR-active intermediates trapped during aerobic assembly of the R2lox Mn/Fe cofactor. A Mn^III(μ-O)(μ-OH)Fe^III species is identified as the first EPR-active species and represents a common state between the two classes of redox-active Mn/Fe proteins. The species downstream from the Mn^III(μ-O)(μ-OH)Fe^III state exhibits unique EPR properties, including unprecedented spectral breadth and isotope-dependent g-tensors, which are attributed to a weakly coupled, hydrogen-bonded Mn^III(μ-OH)Fe^III species. This final intermediate precedes formation of the Mn^III/Fe^III resting state and is suggested to be relevant to understanding the endogenous reactivity of R2lox.

Evidence for Low-Valent Electronic Configurations in Iron-Sulfur Clusters

Although biological iron-sulfur (Fe-S) clusters perform some of the most difficult redox reactions in nature, they are thought to be composed exclusively of Fe²⁺ and Fe³⁺ ions, as well as mixed-valent pairs with average oxidation states of Fe^2.5+. We herein show that Fe-S clusters formally composed of these valences can access a wider range of electronic configurations─in particular, those featuring low-valent Fe¹⁺ centers. We demonstrate that CO binding to a synthetic [Fe₄S₄]⁰ cluster supported by N-heterocyclic carbene ligands induces the generation of Fe¹⁺ centers via intracluster electron transfer, wherein a neighboring pair of Fe²⁺ sites reduces the CO-bound site to a low-valent Fe¹⁺ state. Similarly, CO binding to an [Fe₄S₄]⁺ cluster induces electron delocalization with a neighboring Fe site to form a mixed-valent Fe^1.5+Fe^2.5+ pair in which the CO-bound site adopts partial low-valent character. These low-valent configurations engender remarkable C-O bond activation without having to traverse highly negative and physiologically inaccessible [Fe₄S₄]⁰/[Fe₄S₄]^- redox couples.