Iron-sulfur bond covalency from electronic structure calculations for classical iron-sulfur clusters

Computational Study on the Influence of Mo/V Centers on the Electronic Structure and Hydrazine Reduction Capability of [MFe3S4]3+/2+ Complexes

[MFe3S4] cubanes have for some time been of interest for their ability to mimic the electronic and geometric structure of the active site of nitrogenase, the enzyme responsible for fixing N2 to NH3. Nitrogenase naturally occurs in three forms, with the major difference being that the metal ion present in the cofactor active site is either molybdenum (FeMoco), vanadium (FeVco), or iron. The molybdenum and vanadium versions of these cofactors are more closely studied, owing to their larger abundance and rate of catalysis. In this study, we compare free energy profiles and electronic properties of the Mo/V cubanes at various stages during the reduction of N2H4 to NH3. Our findings highlight the differences in how the complexes facilitate the reaction, in particular, vanadium's comparatively weaker ability to interact with the Fe/S network and stabilize reducing electrons prior to N-N bond cleavage, which may have implications when considering the lower efficiency of the vanadium-dependent nitrogenase.

Evaluation of the Antibacterial Properties of Iron Oxide, Polyethylene Glycol, and Gentamicin Conjugated Nanoparticles against Some Multidrug-Resistant Bacteria

Antibacterial resistance is observed as a public health issue around the world. Every day, new resistance mechanisms appear and spread over the world. For that reason, it is imperative to improve the treatment schemes that have been developed to treat infections caused by wound infections, for instance, Staphylococcus epidermidis (S. epidermidis), Proteus mirabilis (P. mirabilis), and Acinetobacter baumannii (A. baumannii). In this case, we proposed a method that involves mixing the Gentamicin (Gen) with iron oxide nanoparticles (Fe3O4 NPs) and a polymer (polyethylene glycol (PEG)) with Fe3O4 NPs. X-ray diffraction (XRD), Fourier-transform infrared spectroscopy (FTIR), energy dispersive X-ray (EDX), scanning electron microscope (SEM), and transmission electron microscope (TEM) were used to characterize Fe3O4 NPs. Zeta potential and dynamic light scattering (DLS) were also assessed. The antibacterial activity of Fe3O4 NPs, Fe3O4 NPs+PEG, Fe3O4 NPs+Gen, and Fe3O4 NPs+PEG+Gen composites was investigated. The results showed a significant improvement in the antibacterial activity of nanoparticles against bacterial isolates, especially for the Fe3O4 NPs+PEG+Gen as the diameter of the inhibition zone reached 26.33 ± 0.57 mm for A. baumannii, 25.66 ± 0.57 mm for P. mirabilis, and 23.66 ± 0.57 mm for S. epidermidis. The Fe3O4 NPs, Fe3O4 NPs+PEG, Fe3O4+Gen, and Fe3O4+PEG+Gen also showed effectiveness against the biofilm produced by these isolated bacteria. The minimum inhibitory concentration (MIC) of Fe3O4 NPs for S. epidermidis was 25 µg mL−1 and for P. mirabilis and A. baumannii was 50 µg mL−1. The findings suggest that the prepared nanoparticles could be potential therapeutic options for treating wound infections caused by S. epidermidis, P. mirabilis, and A. baumannii.

Analysis of the Geometric and Electronic Structure of Spin-Coupled Iron-Sulfur Dimers with Broken-Symmetry DFT: Implications for FeMoco

The open-shell electronic structure of iron–sulfur clusters presents considerable challenges to quantum chemistry, with the complex iron–molybdenum cofactor (FeMoco) of nitrogenase representing perhaps the ultimate challenge for either wavefunction or density functional theory. While broken-symmetry density functional theory has seen some success in describing the electronic structure of such cofactors, there is a large exchange–correlation functional dependence in calculations that is not fully understood. In this work, we present a geometric benchmarking test set, FeMoD11, of synthetic spin-coupled Fe–Fe and Mo–Fe dimers, with relevance to the molecular and electronic structure of the Mo-nitrogenase FeMo cofactor. The reference data consists of high-resolution crystal structures of metal dimer compounds in different oxidation states. Multiple density functionals are tested on their ability to reproduce the local geometry, specifically the Fe–Fe/Mo–Fe distance, for both antiferromagnetically coupled and ferromagnetically coupled dimers via the broken-symmetry approach. The metal–metal distance is revealed not only to be highly sensitive to the amount of exact exchange in the functional but also to the specific exchange and correlation functionals. For the antiferromagnetically coupled dimers, the calculated metal–metal distance correlates well with the covalency of the bridging metal–ligand bonds, as revealed via the corresponding orbital analysis, Hirshfeld S/Fe charges, and Fe–S Mayer bond order. Superexchange via bridging ligands is expected to be the dominant interaction in these dimers, and our results suggest that functionals that predict accurate Fe–Fe and Mo–Fe distances describe the overall metal–ligand covalency more accurately and in turn the superexchange of these systems. The best performing density functionals of the 16 tested for the FeMoD11 test set are revealed to be either the nonhybrid functionals r²SCAN and B97-D3 or hybrid functionals with 10–15% exact exchange: TPSSh and B3LYP*. These same four functionals are furthermore found to reproduce the high-resolution X-ray structure of FeMoco well according to quantum mechanics/molecular mechanics (QM/MM) calculations. Almost all nonhybrid functionals systematically underestimate Fe–Fe and Mo–Fe distances (with r²SCAN and B97-D3 being the sole exceptions), while hybrid functionals with >15% exact exchange (including range-separated hybrid functionals) overestimate them. The results overall suggest r²SCAN, B97-D3, TPSSh, and B3LYP* as accurate density functionals for describing the electronic structure of iron–sulfur clusters in general, including the complex FeMoco cluster of nitrogenase.

Effect of Oriented Electric Fields on Biologically Relevant Iron-Sulfur Clusters: Tuning Redox Reactivity for Catalysis

Enzyme-based iron-sulfur clusters, exemplified in families such as hydrogenases, nitrogenases, and radical S-adenosylmethionine enzymes, feature in many essential biological processes. The functionality of biological iron-sulfur clusters extends beyond simple electron transfer, relying primarily on the redox activity of the clusters, with a remarkable diversity for different enzymes. The active-site structure and the electrostatic environment in which the cluster resides direct this redox reactivity. Oriented electric fields in enzymatic active sites can be significantly strong, and understanding the extent of their effect on iron-sulfur cluster reactivity can inform first steps toward rationally engineering their reactivity. An extensive systematic density functional theory-based screening approach using OPBE/TZP has afforded a simple electric field-effect representation. The results demonstrate that the orientation of an external electric field of strength 28.8 MV cm^-1 at the center of the cluster can have a significant effect on its relative stability in the order of 35 kJ mol^-1. This shows clear implications for the reactivity of iron-sulfur clusters in enzymes. The results also demonstrate that the orientation of the electric field can alter the most stable broken-symmetry state, which further has implications on the directionality of initiated electron-transfer reactions. These insights open the path for manipulating the enzymatic redox reactivity of iron-sulfur cluster-containing enzymes by rationally engineering oriented electric fields within the enzymes.

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

Carbon monoxide (CO) is a well-known inhibitor of nitrogenase activity. Under turnover conditions, CO binds to FeMoco, the active site of Mo nitrogenase. Time-resolved IR measurements suggest an initial terminal CO at 1904 cm^-1 that converts to a bridging CO at 1715 cm^-1, and an X-ray structure shows that CO can displace one of the bridging belt sulfides of FeMoco. However, the CO-binding redox state(s) of FeMoco (E_n) and the role of the protein environment in stabilizing specific CO-bound intermediates remain elusive. In this work, we carry out an in-depth analysis of the CO-FeMoco interaction based on quantum chemical calculations addressing different aspects of the electronic structure. (1) The local electronic structure of the Fe-CO bond is studied through diamagnetically substituted FeMoco. (2) A cluster model of FeMoco within a polarizable continuum illustrates how CO binding may affect the spin-coupling between the metal centers. (3) A QM/MM model incorporates the explicit influence of the amino acid residues surrounding FeMoco in the MoFe protein. The QM/MM model predicts both a terminal and a bridging CO in the E₁ redox state. The scaled calculated CO frequencies (1922 and 1716 cm^-1, respectively) are in good agreement with the experimentally observed IR bands supporting CO binding to the E₁ state. Alternatively, an E₂ state QM/MM model, which has the same atomic structure as the CO-bound X-ray structure, features a semi-bridging CO with a scaled calculated frequency (1718 cm^-1) similar to the bridging CO in the E₁ model.

The two redox states of the human NEET proteins' [2Fe-2S] clusters

The NEET proteins constitute a unique class of [2Fe–2S] proteins. The metal ions bind to three cysteines and one histidine. The proteins’ clusters exist in two redox states; the oxidized protein (containing two Fe^III ions) can transfer the cluster to apo-acceptor protein(s), while the reduced form (containing one ferrous ion) remains bound to the protein frame. Here, we perform in silico and in vitro studies on human NEET proteins in both reduced and oxidized forms. Quantum chemical calculations on all available human NEET proteins structures suggest that reducing the cluster weakens the Fe–N^His and Fe–S^Cys bonds, similar to what is seen in other Fe–S proteins (e.g., ferredoxin and Rieske protein). We further show that the extra electron in the [2Fe–2S]⁺ clusters of one of the NEET proteins (mNT) is localized on the His-bound iron ion, consistently with our previous spectroscopic studies. Kinetic measurements demonstrate that the mNT [2Fe–2S]⁺ is released only by an increase in temperature. Thus, the reduced state of human NEET proteins [2Fe–2S] cluster is kinetically inert. This previously unrecognized kinetic inertness of the reduced state, along with the reactivity of the oxidized state, is unique across all [2Fe–2S] proteins. Finally, using a coevolutionary analysis, along with molecular dynamics simulations, we provide insight on the observed allostery between the loop L2 and the cluster region. Specifically, we show that W75, R76, K78, K79, F82 and G85 in the latter region share similar allosteric characteristics in both redox states.

Nutrition and sulfur

Sulfur is unusual in that it is a mineral that may be taken into the body in both inorganic and organic combinations. It has been available within the environment throughout the development of lifeforms and as such has become integrated into virtually every aspect of biochemical function. It is essential for the nature and maintenance of structure, assists in communication within the organism, is vital as a catalytic assistant in intermediary metabolism and the mechanism of energy flow as well as being involved in internal defense against potentially damaging reactive species and invading foreign chemicals. Recent studies have suggested extended roles for sulfur-containing molecules within living systems. As such, questions have been raised as to whether or not humans are receiving sufficient sulfur within their diet. Sulfur appears to have been the "poor relation" with regards to mineral nutrition. This may be because of difficulties encountered over its multifarious functions, the many chemical guises in which it may be ingested and its complex biochemical interconversions once taken into the body. No established daily requirements have been determined, unlike many minerals, although suggestions have been proposed. Owing to its widespread distribution within dietary components its intake has almost been taken for granted. In the majority of individuals partaking of a balanced diet the supply is deemed adequate, but those opting for specialized or restrictive diets may experience occasional and low-level shortages. In these instances, the careful use of sulfur supplements may be of benefit.

From Synthetic to Biological Fe4 S4 Complexes: Redox Properties Correlated to Function of Radical S-Adenosylmethionine Enzymes

By employing the computational protocol for calculation of reduction potentials of the Fe₄ S₄ -containing species validated using a representative series of well-defined synthetic complexes, we focused on redox properties of two prototypical radical SAM enzymes to reveal how they transform SAM into the reactive 5'-deoxyadenosyl radical, and how they tune this radical for its proper biological function. We found the reduction potential of SAM is indeed elevated by 0.3-0.4 V upon coordination to Fe₄ S₄ , which was previously speculated in the literature. This makes a generation of 5'-deoxyadenosyl radical from SAM less endergonic (by ca. 7-9 kcal mol^-1 ) and hence more feasible in both enzymes as compared to the identical process in water. Furthermore, our calculations indicate that the enzyme-bound 5'-deoxyadenosyl radical has a significantly lower reduction potential than in referential aqueous solution, which may help the enzymes to suppress potential side redox reactions and simultaneously elevate its proton-philic character, which may, in turn, promote the radical hydrogen-atom abstraction ability.

The Spectroscopy of Nitrogenases

Nitrogenases are responsible for biological nitrogen fixation, a crucial step in the biogeochemical nitrogen cycle. These enzymes utilize a two-component protein system and a series of iron-sulfur clusters to perform this reaction, culminating at the FeMco active site (M = Mo, V, Fe), which is capable of binding and reducing N2 to 2NH3. In this review, we summarize how different spectroscopic approaches have shed light on various aspects of these enzymes, including their structure, mechanism, alternative reactivity, and maturation. Synthetic model chemistry and theory have also played significant roles in developing our present understanding of these systems and are discussed in the context of their contributions to interpreting the nature of nitrogenases. Despite years of significant progress, there is still much to be learned from these enzymes through spectroscopic means, and we highlight where further spectroscopic investigations are needed.

Potential Functions and Thermodynamic Properties of UC, UN, and UH

Potential energy surface scanning for UC, UN, and UH was performed by configuration interaction (CI), coupled cluster singles and doubles (CCSD) excitation, quadratic configuration interaction (QCISD (T)), and density functional theory PBE1 (DFT-PBE1) methods in coupling with the ECP80MWB_AVQZ + 2f basis set for uranium and 6 − 311 + G∗for carbon, hydrogen, and nitrogen. The dissociation energies of UC, UN, and UH are 5.7960, 4.5077, and 2.6999 eV at the QCISD (T) levels, respectively. The calculated energy was fitted to the potential functions of Morse, Lennard-Jones, and Rydberg by using the least square method. The anharmonicity constant of UC is 0.0047160. The anharmonic frequency of UC is 780.27 cm−1which was obtained based on the PBE1 results. For UN, the anharmonicity constant is 0.0049827. The anharmonic frequency is 812.65 cm−1which was obtained through the PBE1 results. For UH, the anharmonicity constant is 0.017300. The anharmonic frequency obtained via the QCISD (T) results is 1449.8 cm−1. The heat capacity and entropy in different temperatures were calculated using anharmonic frequencies. These properties are in good accordance with the direct DFT-UPBE1 results (for UC and UN) and QCISD (T) results (for UH). The relationship of entropy with temperature was established.

The influences of carbon donor ligands on biomimetic multi-iron complexes for N2 reduction

The active site clusters of nitrogenase enzymes possess the only examples of carbides in biology. These are the only biological FeS clusters that are capable of reducing N₂ to NH₄⁺, implicating the central carbon and its interaction with Fe as important in the mechanism of N₂ reduction. This biological question motivates study of the influence of carbon donors on the electronic structure and reactivity of unsaturated, high-spin iron centers. Here, we present functional and structural models that test the impacts of carbon donors and sulfide donors in simpler iron compounds. We report the first example of a diiron complex that is bridged by an alkylidene and a sulfide, which serves as a high-fidelity structural and spectroscopic model of a two-iron portion of the active-site cluster (FeMoco) in the resting state of Mo-nitrogenase. The model complexes have antiferromagnetically coupled pairs of high-spin iron centers, and sulfur K-edge X-ray absorption spectroscopy shows comparable covalency of the sulfide for C and S bridged species. The sulfur-bridged compound does not interact with N₂ even upon reduction, but upon removal of the sulfide it becomes capable of reducing N₂ to NH₄⁺ with the addition of protons and electrons. This provides synthetic support for sulfide extrusion in the activation of nitrogenase cofactors.

High-spin diiron alkylidenes give insight into the electronic structure and functional relevance of carbon in the FeMoco active site of nitrogenase.

Evaluating Density Functionals by Examining Molecular Structures, Chemical Bonding, and Relative Energies of Mononuclear Ru-Cl-H-PR3 Isomers

In order to define a robust level of theory using density functionals for investigating the reactivity of ruthenium complexes, we used benchmark wave function theory, with saturated basis sets to validate generalized gradient approximation (GGA), meta-GGA, and hyper-GGA functionals in the presence and absence of empirical dispersion and range-separated corrections. We first selected potentially suitable functionals that gave accurate predictions of the relative energetics of coordination isomers. These functionals were further evaluated for the chemical accuracy of their predicted geometric and electronic structures. For the latter, both the ionic and covalent interactions were considered. The reference level of theory for comparison was coupled-cluster perturbation theory using full treatment of singles and doubles (CCSD) with a saturated triple-ζ quality basis set (TZVP) and corresponding small-core, effective core potentials for ruthenium. Several population analysis methods were evaluated to predict the ionic interactions. We found that the atomic charges obtained from fitting the electrostatic potential provided the most reasonable estimates for the ruthenium complexes. The covalent interactions were quantified by considering the atomic compositions of Ru 4d _{x²- y²}- and 4d _z²-based frontier unoccupied orbitals. Comparison of more than two dozen functionals with reference data from high-level wave function calculations revealed trends that allowed for the formulation of an optimal hybrid density functional: PBE exchange and correlation functionals with 50% HF exchange component. This level of theory was found to reproduce the experimental structure of Ru^(II) complexes. These complexes were used to investigate chemical speciation in a simplified model for an ionic liquid environment.

Magnetically Separable MoS₂/Fe₃O₄/nZVI Nanocomposites for the Treatment of Wastewater Containing Cr(VI) and 4-Chlorophenol

With a large specific surface area, high reactivity, and excellent adsorption properties, nano zerovalent iron (nZVI) can degrade a wide variety of contaminants in wastewater. However, aggregation, oxidation, and separation issues greatly impede its wide application. In this study, MoS₂/Fe₃O₄/nZVI nanocomposites were successfully synthesized by a facile step-by-step approach to overcome these problems. MoS₂ nanosheets (MNs) acted as an efficient support for nZVI and enriched the organic pollutants nearby, leading to an enhanced removal efficiency. Fe₃O₄ nanoparticles (NPs) could not only suppress the agglomeration and restacking of MNs, but also facilitate easy separation and recovery of the nanocomposites. The synergistic effect between MNs and Fe₃O₄ NPs effectively enhanced the reactivity and efficiency of nZVI. In the system, Cr(VI) was reduced to Cr(III) by nZVI in the nanocomposites, and Fe2+ produced in the process was combined with H₂O₂ to further remove 4-Chlorophenol (4-CP) through a Fenton reaction. Furthermore, the nanocomposites could be easily separated from wastewater by a magnet and be reused for at least five consecutive runs, revealing good reusability. The results demonstrate that the novel nanocomposites are highly efficient and promising for the simultaneous removal of Cr(VI) and 4-CP in wastewater.

Spin-Polarization-Induced Preedge Transitions in the Sulfur K-Edge XAS Spectra of Open-Shell Transition-Metal Sulfates: Spectroscopic Validation of σ-Bond Electron Transfer

Sulfur K-edge X-ray absorption spectroscopy (XAS) spectra of the monodentate sulfate complexes [MII(itao)(SO4)(H2O)0,1] (M = Co, Ni, Cu) and [Cu(Me6tren)(SO4)] exhibit well-defined preedge transitions at 2479.4, 2479.9, 2478.4, and 2477.7 eV, respectively, despite having no direct metal-sulfur bond, while the XAS preedge of [Zn(itao)(SO4)] is featureless. The sulfur K-edge XAS of [Cu(itao)(SO4)] but not of [Cu(Me6tren)(SO4)] uniquely exhibits a weak transition at 2472.1 eV, an extraordinary 8.7 eV below the first inflection of the rising K-edge. Preedge transitions also appear in the sulfur K-edge XAS of crystalline [MII(SO4)(H2O)] (M = Fe, Co, Ni, and Cu, but not Zn) and in sulfates of higher-valent early transition metals. Ground-state density functional theory (DFT) and time-dependent DFT (TDDFT) calculations show that charge transfer from coordinated sulfate to paramagnetic late transition metals produces spin polarization that differentially mixes the spin-up (α) and spin-down (β) spin orbitals of the sulfate ligand, inducing negative spin density at the sulfate sulfur. Ground-state DFT calculations show that sulfur 3p character then mixes into metal 4s and 4p valence orbitals and various combinations of ligand antibonding orbitals, producing measurable sulfur XAS transitions. TDDFT calculations confirm the presence of XAS preedge features 0.5-2 eV below the rising sulfur K-edge energy. The 2472.1 eV feature arises when orbitals at lower energy than the frontier occupied orbitals with S 3p character mix with the copper(II) electron hole. Transmission of spin polarization and thus of radical character through several bonds between the sulfur and electron hole provides a new mechanism for the counterintuitive appearance of preedge transitions in the XAS spectra of transition-metal oxoanion ligands in the absence of any direct metal-absorber bond. The 2472.1 eV transition is evidence for further radicalization from copper(II), which extends across a hydrogen-bond bridge between sulfate and the itao ligand and involves orbitals at energies below the frontier set. This electronic structure feature provides a direct spectroscopic confirmation of the through-bond electron-transfer mechanism of redox-active metalloproteins.

Charge-density analysis of an iron-sulfur protein at an ultra-high resolution of 0.48 Å

The fine structures of proteins, such as the positions of hydrogen atoms, distributions of valence electrons and orientations of bound waters, are critical factors for determining the dynamic and chemical properties of proteins. Such information cannot be obtained by conventional protein X-ray analyses at 3.0-1.5 Å resolution, in which amino acids are fitted into atomically unresolved electron-density maps and refinement calculations are performed under strong restraints. Therefore, we usually supplement the information on hydrogen atoms and valence electrons in proteins with pre-existing common knowledge obtained by chemistry in small molecules. However, even now, computational calculation of such information with quantum chemistry also tends to be difficult, especially for polynuclear metalloproteins. Here we report a charge-density analysis of the high-potential iron-sulfur protein from the thermophilic purple bacterium Thermochromatium tepidum using X-ray data at an ultra-high resolution of 0.48 Å. Residual electron densities in the conventional refinement are assigned as valence electrons in the multipolar refinement. Iron 3d and sulfur 3p electron densities of the Fe4S4 cluster are visualized around the atoms. Such information provides the most detailed view of the valence electrons of the metal complex in the protein. The asymmetry of the iron-sulfur cluster and the protein environment suggests the structural basis of charge storing on electron transfer. Our charge-density analysis reveals many fine features around the metal complex for the first time, and will enable further theoretical and experimental studies of metalloproteins.

A New Nitrogenase Mechanism Using a CFe8S9 Model: Does H2 Elimination Activate the Complex to N2 Addition to the Central Carbon Atom?

A truncated model of the FeMo cofactor is used to explore a new mechanism for the conversion of N2 to NH3 by the nitrogenase enzyme. After four initial protonation/reduction steps, the H4CFe8S9 cluster has two hydrogen atoms attached to sulfur, one hydrogen bridging two iron centers and one hydrogen bonded to carbon. The loss of the CH and FeHFe hydrogens as molecular hydrogen activates the cluster to addition of N2 to the carbon center. This unique step takes place at a nearly planar four-coordinate carbon center and leads to an intermediate with a significantly weakened N-N bond. A hydrogen attached to a sulfur atom is then transferred to the distal nitrogen atom. Additional prontonation/reduction steps are modeled by adding a hydrogen atom to sulfur and locating the transition states for transfer to nitrogen. The first NH3 is lost in a thermal neutral step, while the second step is endothermic. The loss of H2 activates the complex by reducing the barrier for N2 addition by 3.5 kcal/mol. Since this is the most difficult step in the mechanism, reducing the barrier for this step justifies the "extra expense" of H2 production.

Mössbauer spectroscopy of Fe/S proteins

Iron-sulfur (Fe/S) clusters are structurally and functionally diverse cofactors that are found in all domains of life. (57)Fe Mössbauer spectroscopy is a technique that provides information about the chemical nature of all chemically distinct Fe species contained in a sample, such as Fe oxidation and spin state, nuclearity of a cluster with more than one metal ion, electron spin ground state of the cluster, and delocalization properties in mixed-valent clusters. Moreover, the technique allows for quantitation of all Fe species, when it is used in conjunction with electron paramagnetic resonance (EPR) spectroscopy and analytical methods. (57)Fe-Mössbauer spectroscopy played a pivotal role in unraveling the electronic structures of the "well-established" [2Fe-2S](2+/+), [3Fe-4S](1+/0), and [4Fe-4S](3+/2+/1+/0) clusters and -more-recently- was used to characterize novel Fe/S clustsers, including the [4Fe-3S] cluster of the O2-tolerant hydrogenase from Aquifex aeolicus and the 3Fe-cluster intermediate observed during the reaction of lipoyl synthase, a member of the radical SAM enzyme superfamily.