Revisiting the catalytic mechanism of Mo-Cu carbon monoxide dehydrogenase using QM/MM and DFT calculations

Molybdenum-Copper Antagonism In Metalloenzymes And Anti-Copper Therapy

The connection between 3d (Cu) and 4d (Mo) via the "Mo-S-Cu" unit is called Mo-Cu antagonism. Biology offers case studies of such interactions in metalloproteins such as Mo/Cu-CO Dehydrogenases (Mo/Cu-CODH), and Mo/Cu Orange Protein (Mo/Cu-ORP). The CODH significantly maintains the CO level in the atmosphere below the toxic level by converting it to non-toxic CO2 for respiring organisms. Several models were synthesized to understand the structure-function relationship of these native enzymes. However, this interaction was first observed in ruminants, and they convert molybdate (MoO4 2- ) into tetrathiomolybdate (MoS4 2- ; TTM), reacting with cellular Cu to yield biological unavailable Mo/S/Cu cluster, then developing Cu-deficiency diseases. These findings inspire the use of TTM as a Cu-sequester drug, especially for treating Cu-dependent human diseases such as Wilson diseases (WD) and cancer. It is well known that a balanced Cu homeostasis is essential for a wide range of biological processes, but negative consequence leads to cell toxicity. Therefore, this review aims to connect the Mo-Cu antagonism in metalloproteins and anti-copper therapy.

Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set

We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.

Studies Relevant to the Functional Model of Mo-Cu CODH: In Situ Reactions of Cu(I)-L Complexes with Mo(VI) and Synthesis of Stable Structurally Characterized Heterotetranuclear MoVI2CuI2 Complex

In this study, we report the synthesis, characterization, and reactions of Cu(I) complexes of the general form Cu(L)(LigH₂) (LigH₂ = xanthene-based heterodinucleating ligand (E)-3-(((5-(bis(pyridin-2-ylmethyl)amino)-2,7-di-tert-butyl-9,9-dimethyl-9H-xanthen-4-yl)imino)methyl)benzene-1,2-diol); L = PMe₃, PPh₃, CN(2,6-Me₂C₆H₃)). New complexes [Cu(PMe₃)(LigH₂)] and [CuCN(2,6-Me₂C₆H₃)(LigH₂)] were synthesized by treating [Cu(LigH₂)](PF₆) with trimethylphosphine and 2,6-dimethylphenyl isocyanide, respectively. These complexes were characterized by multinuclear NMR spectroscopy, IR spectroscopy, high-resolution mass spectrometry (HRMS), and X-ray crystallography. In contrast, attempted reactions of [Cu(LigH₂)](PF₆) with cyanide or styrene failed to produce isolable crystalline products. Next, the reactivity of these and previously synthesized Cu(I) phosphine and isocyanide complexes with molybdate was interrogated. IR (for isocyanide) and ³¹P NMR (for PPh₃/PMe₃) spectroscopy demonstrates the lack of oxidation reactivity. We also describe herein the first example of a structurally characterized multinuclear complex combining both Mo(VI) and Cu(I) metal ions within the same system. The heterobimetallic tetranuclear complex [Cu₂Mo₂O₄(μ₂-O)(Lig)₂]·HOSiPh₃ was obtained by the reaction of the silylated Mo(VI) precursor (Et₄N)(MoO₃(OSiPh₃)) with LigH₂, followed by the addition of [Cu(NCMe)₄](PF₆). This complex was characterized by NMR spectroscopy, high-resolution mass spectrometry, and X-ray crystallography.

Resonance Raman spectroscopy of pyranopterin molybdenum enzymes

Resonance Raman spectroscopy (rR) is a powerful spectroscopic probe that is widely used for studying the geometric and electronic structure of metalloproteins. In this focused review, we detail how resonance Raman spectroscopy has contributed to a greater understanding of electronic structure, geometric structure, and the reaction mechanisms of pyranopterin molybdenum enzymes. The review focuses on the enzymes sulfite oxidase (SO), dimethyl sulfoxide reductase (DMSOR), xanthine oxidase (XO), and carbon monoxide dehydrogenase. Specifically, we highlight how Mo-Ooxo, Mo-Ssulfido, Mo-Sdithiolene, and dithiolene CC vibrational modes, isotope and heavy atom perturbations, resonance enhancement, and associated Raman studies of small molecule analogs have provided detailed insight into the nature of these metalloenzyme active sites.

Microbial oxidation of atmospheric trace gases

The atmosphere has recently been recognized as a major source of energy sustaining life. Diverse aerobic bacteria oxidize the three most abundant reduced trace gases in the atmosphere, namely hydrogen (H₂), carbon monoxide (CO) and methane (CH₄). This Review describes the taxonomic distribution, physiological role and biochemical basis of microbial oxidation of these atmospheric trace gases, as well as the ecological, environmental, medical and astrobiological importance of this process. Most soil bacteria and some archaea can survive by using atmospheric H₂ and CO as alternative energy sources, as illustrated through genetic studies on Mycobacterium cells and Streptomyces spores. Certain specialist bacteria can also grow on air alone, as confirmed by the landmark characterization of Methylocapsa gorgona, which grows by simultaneously consuming atmospheric CH₄, H₂ and CO. Bacteria use high-affinity lineages of metalloenzymes, namely hydrogenases, CO dehydrogenases and methane monooxygenases, to utilize atmospheric trace gases for aerobic respiration and carbon fixation. More broadly, trace gas oxidizers enhance the biodiversity and resilience of soil and marine ecosystems, drive primary productivity in extreme environments such as Antarctic desert soils and perform critical regulatory services by mitigating anthropogenic emissions of greenhouse gases and toxic pollutants.

Can water act as a nucleophile in CO oxidation catalysed by Mo/Cu CO-dehydrogenase? Answers from theory

The aerobic CO dehydrogenase from Oligotropha carboxidovorans is an environmentally crucial bacterial enzyme for maintenance of subtoxic concentration of CO in the lower atmosphere, as it allows for the oxidation of CO to CO₂ which takes place at its Mo−Cu heterobimetallic active site. Despite extensive experimental and theoretical efforts, significant uncertainties still concern the reaction mechanism for the CO oxidation. In this work, we used the hybrid quantum mechanical/molecular mechanical approach to evaluate whether a water molecule present in the active site might act as a nucleophile upon formation of the new C−O bond, a hypothesis recently suggested in the literature. Our study shows that activation of H₂O can be favoured by the presence of the Mo=O_eq group. However, overall our results suggest that mechanisms other than the nucleophilic attack by Mo=O_eq to the activated carbon of the CO substrate are not likely to constitute reactive channels for the oxidation of CO by the enzyme.

Synthesis and Cu(I)/Mo(VI) Reactivity of a Bifunctional Heterodinucleating Ligand on a Xanthene Platform

In an effort to probe the feasibility of a model of Mo-Cu CODH (CODH = carbon monoxide dehydrogenase) lacking a bridging sulfido group, the new heterodinucleating ligand LH₂ was designed and its Cu(I)/Mo(VI) reactivity was investigated. LH₂ ((E)-3-(((5-(bis(pyridin-2-ylmethyl)amino)-2,7-di-tert-butyl-9,9-dimethyl-9H-xanthen-4-yl)imino)methyl)benzene-1,2-diol) features two different chelating positions bridged by a xanthene linker: bis(pyridyl)amine for Cu(I) and catecholate for Mo(VI). LH₂ was synthesized via the initial protection of one of the amine positions, followed by two consecutive alkylations of the second position, deprotection, and condensation to attach the catechol functionality. LH₂ was found to exhibit dynamic cooperativity between two reactive sites mediated by H-bonding of the catechol protons. In the free ligand, catechol protons exhibit H-bonding with imine (intramolecular) and with pyridine (intermolecular in the solid state). The reaction of LH₂ with [Cu(NCMe)₄]⁺ led to the tetradentate coordination of Cu(I) via all nitrogen donors of the ligand, including the imine. Cu(I) complexes were characterized by multinuclear NMR spectroscopy, high-resolution mass spectrometry (HRMS), X-ray crystallography, and DFT calculations. Cu(I) coordination to the imine disrupted H-bonding and caused rotation away from the catechol arm. The reaction of the Cu(I) complex [Cu(LH₂)]⁺ with a variety of monodentate ligands X (PPh₃, Cl^-, SCN^-, CN^-) released the metal from coordination to the imine, thereby restoring imine H-bonding with the catechol proton. The second catechol proton engages in H-bonding with Cu-X (X = Cl, CN, SCN), which can be intermolecular (XRD) or intramolecular (DFT). The reaction of LH₂ with molybdate [MoO₄]^2- led to incorporation of [Mo^VIO₃] at the catecholate position, producing [MoO₃(L)]^2-. Similarly, the reaction of [Cu(LH₂)]⁺ with [MoO₄]^2- formed the heterodinuclear complex [CuMoO₃(L)]^-. Both complexes were characterized by multinuclear NMR, UV-vis, and HRMS. HRMS in both cases confirmed the constitution of the complexes, containing molecular ions with the expected isotopic distribution.

A W/Cu Synthetic Model for the Mo/Cu Cofactor of Aerobic CODH Indicates That Biochemical CO Oxidation Requires a Frustrated Lewis Acid/Base Pair

Constructing synthetic models of the Mo/Cu active site of aerobic carbon monoxide dehydrogenase (CODH) has been a long-standing synthetic challenge thought to be crucial for understanding how atmospheric concentrations of CO and CO₂ are regulated in the global carbon cycle by chemolithoautotrophic bacteria and archaea. Here we report a W/Cu complex that is among the closest synthetic mimics constructed to date, enabled by a silyl protection/deprotection strategy that provided access to a kinetically stabilized complex with mixed O^2-/S^2- ligation between (bdt)(O)W^VI and Cu^I(NHC) (bdt = benzene dithiolate, NHC = N-heterocyclic carbene) sites. Differences between the inorganic core's structural and electronic features outside the protein environment relative to the native CODH cofactor point to a biochemical CO oxidation mechanism that requires a strained active site geometry, with Lewis acid/base frustration enforced by the protein secondary structure. This new mechanistic insight has the potential to inform synthetic design strategies for multimetallic energy storage catalysts.

Energetics for the Mechanism of Nickel-Containing Carbon Monoxide Dehydrogenase

Nickel-containing carbon monoxide (CO) dehydrogenase is an enzyme that catalyzes the important reversible carbon dioxide reduction. Several high-resolution structures have been determined at various stages of the reduction, which can be used as good starting points for the present computational study. The cluster model is used in combination with a systematic application of the density functional theory as recently described. The results are in very good agreement with experimental evidence. There are a few important results. To explain why the X-ray structure for the reduced C_red1 state has an empty site on nickel, it is here suggested that the cluster has been over-reduced by X-rays and is therefore not the desired reduced state, which instead contains a bound CO on nickel. After an additional reduction, a hydride bound to nickel is suggested to play a role. In order to obtain energetics in agreement with experiments, it is concluded that one sulfide bridge in the Ni-Fe cluster should be protonated. The best test of the accuracy obtained is to compare the computed rate for reduction using -0.6 V with that for oxidation using -0.3 V, where good agreement was obtained. Obtaining a mechanism that is easily reversible is another demanding aspect of the modeling. Nickel oscillates between nickel(II) and nickel(I), while nickel(0) never comes in.

The Challenging in silico Description of Carbon Monoxide Oxidation as Catalyzed by Molybdenum-Copper CO Dehydrogenase

Carbon monoxide (CO) is a highly toxic gas to many living organisms. However, some microorganisms are able to use this molecule as the sole source of carbon and energy. Soil bacteria such as the aerobic Oligotropha carboxidovorans are responsible for the annual removal of about 2x108 tons of CO from the atmosphere. Detoxification through oxidation of CO to CO2 is enabled by the MoCu-dependent CO-dehydrogenase enzyme (MoCu-CODH) which-differently from other enzyme classes with similar function-retains its catalytic activity in the presence of atmospheric O2. In the last few years, targeted advancements have been described in the field of bioengineering and biomimetics, which is functional for future technological exploitation of the catalytic properties of MoCu-CODH and for the reproduction of its reactivity in synthetic complexes. Notably, a growing interest for the quantum chemical investigation of this enzyme has recently also emerged. This mini-review compiles the current knowledge of the MoCu-CODH catalytic cycle, with a specific focus on the outcomes of theoretical studies on this enzyme class. Rather controversial aspects from different theoretical studies will be highlighted, thus illustrating the challenges posed by this system as far as the application of density functional theory and hybrid quantum-classical methods are concerned.