Which functional groups of the molybdopterin ligand should be considered when modeling the active sites of the molybdenum and tungsten cofactors? A density functional theory study

Delving into the catalytic mechanism of molybdenum cofactors: a novel coupled cluster study

In this work, we use modern electronic structure methods to model the catalytic mechanism of different variants of the molybdenum cofactor (Moco). We investigate the dependence of various Moco model systems on structural relaxation and the importance of environmental effects for five critical points along the reaction coordinate with the DMSO and NO3- substrates. Furthermore, we scrutinize the performance of various coupled-cluster approaches for modeling the relative energies along the investigated reaction paths, focusing on several pair coupled cluster doubles (pCCD) flavors and conventional coupled cluster approximations. Moreover, we elucidate the Mo-O bond formation using orbital-based quantum information measures, which highlight the flow of σM-O bond formation and σN/S-O bond breaking. Our study shows that pCCD-based models are a viable alternative to conventional methods and offer us unique insights into the bonding situation along a reaction coordinate. Finally, this work highlights the importance of environmental effects or changes in the core and, consequently, in the model itself to elucidate the change in activity of different Moco variants.

N2 Reduction versus H2 Evolution in a Molybdenum- or Tungsten-Based Small-Molecule Model System of Nitrogenase

Molybdenum dinitrogen complexes have played a major role as catalytic model systems of nitrogenase. In comparison, analogous tungsten complexes have in most cases found to be catalytically inactive. Herein, a tungsten complex was shown to be supported by a pentadentate tetrapodal (pentaPod) phosphine ligand, under conditions of N₂ fixation, primarily catalyzes the hydrogen evolution reaction (HER), in contrast to its Mo analogue, which catalytically mediates the nitrogen-reduction reaction (N₂ RR). DFT calculations were employed to evaluate possible mechanisms and identify the most likely pathways of N₂ RR and HER activities exhibited by Mo- and W-pentaPod complexes. Two mechanisms for N₂ RR by PCET are considered, starting from neutral (M(0) cycle) and cationic (M(I) cycle) dinitrogen complexes (M=Mo, W). The latter was found to be energetically more favorable. For HER three scenarios are treated; that is, through bimolecular reactions of early M-N_x H_y intermediates, pure hydride intermediates or mixed M(H)(N_x H_y ) species.

Formate Dehydrogenase Mimics as Catalysts for Carbon Dioxide Reduction

Formate dehydrogenases (FDH) reversibly catalyze the interconversion of CO₂ to formate. They belong to the family of molybdenum and tungsten-dependent oxidoreductases. For several decades, scientists have been synthesizing structural and functional model complexes inspired by these enzymes. These studies not only allow for finding certain efficient catalysts but also in some cases to better understand the functioning of the enzymes. However, FDH models for catalytic CO₂ reduction are less studied compared to the oxygen atom transfer (OAT) reaction. Herein, we present recent results of structural and functional models of FDH.

Carbon Dioxide Reduction: A Bioinspired Catalysis Approach

While developed in a number of directions, bioinspired catalysis has been explored only very recently for CO₂ reduction, a challenging reaction of prime importance in the context of the energetic transition to be built up. This approach is particularly relevant because nature teaches us that CO₂ reduction is possible, with low overpotentials, high rates, and large selectivity, and gives us unique clues to design and discover new interesting molecular catalysts. Indeed, on the basis of our relatively advanced understanding of the structures and mechanisms of the active sites of fascinating metalloenzymes such as formate dehydrogenases (FDHs) and CO dehydrogenases (CODHs), it is possible to design original, active, selective, and stable molecular catalysts using the bioinspired approach. These metalloenzymes use fascinating metal centers: in FDHs, a Mo(W) mononuclear ion is coordinated by four sulfur atoms provided by a specific organic ligand, molybdopterin (MPT), containing a pyranopterin heterocycle (composed of a pyran ring fused with a pterin unit) and two sulfhydryl groups for metal chelation; in CODHs, catalytic activity depends on either a unique nickel-iron-sulfur cluster or a dinuclear Mo-Cu complex in which the Mo ion is chelated by an MPT ligand. As a consequence, the novel class of catalysts, designed by bioinspiration, consists of mononuclear Mo, W, and Ni and as well as dinuclear Mo-Cu and Ni-Fe complexes in which the metal ions are coordinated by sulfur ligands, more specifically, dithiolene chelates mimicking the natural MPT cofactor. In general, their activity is evaluated in electrochemical systems (cyclic voltammetry and bulk electrolysis) or in photochemical systems (in the presence of a photosensitizer and a sacrificial electron donor) in solution. This research is multidisciplinary because it implies detailed biochemical, functional, and structural characterization of the inspiring enzymes together with synthetic organic and organometallic chemistry and molecular catalysis studies. The most important achievements in this direction, starting from the first report of a catalytically active biomimetic bis-dithiolene-Mo complex in 2015, are discussed in this Account, highlighting the challenging issues associated with synthesis of such sophisticated ligands and molecular catalysts as well as the complexity of reaction mechanisms. While the very first active biomimetic catalysts require further improvement, in terms of performance, they set the stage in which molecular chemistry and enzymology can synergistically cooperate for a better understanding of why nature has selected these sites and for developing highly active catalysts.

Oxygen Atom Transfer Reactivity of Molybdenum(VI) Complexes Employing Pyrimidine- and Pyridine-2-thiolate Ligands

Four dioxidomolybdenum(VI) complexes of the general structure [MoO₂L₂] employing the S,N-bidentate ligands pyrimidine-2-thiolate (PymS, 1), pyridine-2-thiolate (PyS, 2), 4-methylpyridine-2-thiolate (4-MePyS, 3) and 6-methylpyridine-2-thiolate (6-MePyS, 4) were synthesized and characterized by spectroscopic means and single-crystal X-ray diffraction analysis (2-4). Complexes 1-4 were reacted with PPh₃ and PMe₃, respectively, to investigate their oxygen atom transfer (OAT) reactivity and catalytic applicability. Reduction with PPh₃ leads to symmetric molybdenum(V) dimers of the general structure [Mo₂O₃L₄] (6-9). Kinetic studies showed that the OAT from [MoO₂L₂] to PPh₃ is 5 times faster for the PymS system than for the PyS and 4-MePyS systems. The reaction of complexes 1-3 with PMe₃ gives stable molybdenum(IV) complexes of the structure [MoOL₂(PMe₃)₂] (10-12), while reduction of [MoO₂(6-MePyS)₂] (4) yields [MoO(6-MePyS)₂(PMe₃)] (13) with only one PMe₃ coordinated to the metal center. The activity of complexes 1-4 in catalytic OAT reactions involving Me₂SO and Ph₂SO as oxygen donors and PPh₃ as an oxygen acceptor has been investigated to assess the influence of the varied ligand frameworks on the OAT reaction rates. It was found that [MoO₂(PymS)₂] (1) and [MoO₂(6-MePyS)₂] (4) are similarly efficient catalysts, while complexes 2 and 3 are only moderately active. In the catalytic oxidation of PMe₃ with Me₂SO, complex 4 is the only efficient catalyst. Complexes 1-4 were also found to catalytically reduce NO₃^- with PPh₃, although their reactivity is inhibited by further reduced species such as NO, as exemplified by the formation of the nitrosyl complex [Mo(NO)(PymS)₃] (14), which was identified by single-crystal X-ray diffraction analysis. Computed ΔG^⧧ values for the very first step of the OAT were found to be lower for complexes 1 and 4 than for 2 and 3, explaining the difference in catalytic reactivity between the two pairs and revealing the requirement for an electron-deficient ligand system.

Computational Investigation into the Ni(SeNHC2(CN)2)2 and Ni(SNHC2(CN)2)2 Complexes as Potential Catalysts for Hydrogen Production

To reduce our carbon footprint, we must look at alternative non-carbon-containing fuels to prevent continued global climate change. One environmentally friendly alternative fuel is molecular hydrogen. Herein the Ni(SeNHC₂(CN)₂)₂ complex was studied using DFT to determine the thermodynamics associated with the electrocatalytic formation of H₂(g). From the calculated thermodynamics, it appears that the Ni(SeNHC₂(CN)₂)₂ complex is predicted to catalyze the production of H₂ gas under mildly reducing conditions relative to the SHE. Notably, the thermodynamics are better than the values calculated for the analogous Ni(SNHC₂(CN)₂)₂ complex which has been shown experimentally to catalyze the formation of H₂ gas in aqueous solution. Regarding possible kinetic reactivity, the HOMO-LUMO gap energies were calculated. From the gap energies, it is expected that the Se-containing compounds would be more reactive to electron transfer in the third reduction step, meaning therefore that a smaller overpotential would be needed to drive the reduction of Red2-H₂ relative to ^SRed2-H₂ in agreement with past experimental work. Thus, the use of Se in such compounds may offer a means to improve the catalysts for H₂ production.

A computational investigation into the catalytic activity of a diselenolene sulfite oxidase biomimetic complex

Molybdenum is the only 4d metal found in almost all life. One such molybdenum-containing enzyme is sulfite oxidase, which also contains the dithiolene-molybdopterin ligand. Sulfite oxidase is essential in the degradation of sulfur-containing compounds such as cysteine and methionine. Past work has shown parallels in the chemistry of dithiolene–metal and diselenolene–metal complexes. Thus, in this present work, the oxygen atom transfer mechanism for a diselenolene sulfite oxidase biomimetic complex was investigated using computational tools, the results of which were compared to the analogous dithiolene biomimetic complex. From the results obtained, the molybdenum-diselenolene sulfite oxidase biomimetic complex is able to catalyse the oxygen atom transfer and does so with a marginally lower value of ΔrG‡ than that for the analogous dithiolene complex. In particular, it was found that on average, the diselenolene complex had an activation energy 1.2 kJ mol–1 lower in energy than the analogous dithiolene complex. However, the calculated value of ΔrG suggests that the oxidation of sulfite is more favourable for the dithiolene complex where the average difference in reaction aqueous Gibbs reaction energy was –9.4 kJ mol–1 relative to the diselenolene complex. It is noted that with the use of D3 and D3BJ corrections in combination with the B3LYP functional, the barrier for oxygen atom transfer is lowered by more than 30.0 kJ mol–1 for both the diselenolene and dithiolene complexes. Such results suggest that to study such oxo-transfer reactions, the proper treatment of dispersion interaction is necessary.

A Bioinspired Molybdenum Complex as a Catalyst for the Photo- and Electroreduction of Protons

A molybdenum-dithiolene-oxo complex was prepared as a model of some active sites of Mo/W-dependent enzymes. The ligand, a quinoxaline-pyran-fused dithiolene, mimics molybdopterin present in these active sites. For the first time, this type of complex was shown to be active as a catalyst for the photoreduction of protons with excellent turnover numbers (500) and good stability in aqueous/organic media and for the electroreduction of protons in acetonitrile with remarkable rate constants (1030 s(-1) at -1.3 V versus Ag/AgCl). DFT calculations provided insight into the catalytic cycle of the reaction, suggesting that the oxo ligand plays a key role in proton exchange. These results provide a basis to optimize this new class of H2 -evolving catalysts.

Pulsed electron paramagnetic resonance spectroscopy of (33)S-labeled molybdenum cofactor in catalytically active bioengineered sulfite oxidase

Molybdenum enzymes contain at least one pyranopterin dithiolate (molybdopterin, MPT) moiety that coordinates Mo through two dithiolate (dithiolene) sulfur atoms. For sulfite oxidase (SO), hyperfine interactions (hfi) and nuclear quadrupole interactions (nqi) of magnetic nuclei (I ≠ 0) near the Mo(V) (d(1)) center have been measured using high-resolution pulsed electron paramagnetic resonance (EPR) methods and interpreted with the help of density functional theory (DFT) calculations. These have provided important insights about the active site structure and the reaction mechanism of the enzyme. However, it has not been possible to use EPR to probe the dithiolene sulfurs directly since naturally abundant (32)S has no nuclear spin (I = 0). Here we describe direct incorporation of (33)S (I = 3/2), the only stable magnetic sulfur isotope, into MPT using controlled in vitro synthesis with purified proteins. The electron spin echo envelope modulation (ESEEM) spectra from (33)S-labeled MPT in this catalytically active SO variant are dominated by the "interdoublet" transition arising from the strong nuclear quadrupole interaction, as also occurs for the (33)S-labeled exchangeable equatorial sulfite ligand [ Klein, E. L., et al. Inorg. Chem. 2012 , 51 , 1408 - 1418 ]. The estimated experimental hfi and nqi parameters for (33)S (aiso = 3 MHz and e(2)Qq/h = 25 MHz) are in good agreement with those predicted by DFT. In addition, the DFT calculations show that the two (33)S atoms are indistinguishable by EPR and reveal a strong intermixing between their out-of-plane pz orbitals and the dxy orbital of Mo(V).

Pterin chemistry and its relationship to the molybdenum cofactor

The molybdenum cofactor is composed of a molybdenum coordinated by one or two rather complicated ligands known as either molybdopterin or pyranopterin. Pterin is one of a large family of bicyclic N-heterocycles called pteridines. Such molecules are widely found in Nature, having various forms to perform a variety of biological functions. This article describes the basic nomenclature of pterin, their biological roles, structure, chemical synthesis and redox reactivity. In addition, the biosynthesis of pterins and current models of the molybdenum cofactor are discussed.