Density functional theory studies of model complexes for molybdenum-dependent nitrate reductase active sites

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.

Microbial Denitrification: Active Site and Reaction Path Models Predict New Isotopic Fingerprints

The study of isotopic fingerprints in nitrate (δ15N, δ18O, Δ17O) has enabled pivotal insights into the global nitrogen cycle and revealed new knowledge gaps. Measuring populations of isotopic homologs of intact NO3 - ions (isotopologues) shows promise to advance the understanding of nitrogen cycling processes; however, we need new theory and predictions to guide laboratory experiments and field studies. We investigated the hypothesis that the isotopic composition of the residual nitrate pool is controlled by the N-O bond-breaking step in Nar dissimilatory nitrate reductase using molecular models of the enzyme active sites and associated kinetic isotope effects (KIEs). We integrated the molecular model results into reaction path models representing the reduction of nitrate under either closed-system or steady-state conditions. The predicted intrinsic KIE (15ε and 18ε) of the Nar active site matches observed fractionations in both culture and environmental studies. This is what would be expected if the isotopic composition of marine nitrate were controlled by dissimilatory nitrate reduction by Nar. For a closed system, the molecular models predict a pronounced negative 15N-18O clumping anomaly in residual nitrate. This signal could encode information about the amount of nitrate consumed in a closed system and thus constrain initial nitrate concentration and its isotopic composition. Similar clumped isotope anomalies can potentially be used to distinguish whether a system is open or closed to new nitrate addition. These mechanistic predictions can be tested and refined in combination with emerging ESI-Orbitrap measurements.

Addressing Ligand-Based Redox in Molybdenum-Dependent Methionine Sulfoxide Reductase

A combination of pulsed EPR, CW EPR, and X-ray absorption spectroscopies has been employed to probe the geometric and electronic structure of the E. coli periplasmic molybdenum-dependent methionine sulfoxide reductase (MsrP). ¹⁷O and ¹H pulsed EPR spectra show that the as-isolated Mo(V) enzyme form does not possess an exchangeable H₂O/OH^- ligand bound to Mo as found in the sulfite oxidizing enzymes of the same family. The nature of the unusual CW EPR spectrum has been re-evaluated in light of new data on the MsrP-N45R variant and related small-molecule analogues of the active site. These data point to a novel "thiol-blocked" [(PDT)Mo^VO(S_Cys)(thiolate)]^- structure, which is supported by new EXAFS data. We discuss these new results in the context of ligand-based and metal-based redox chemistry in the enzymatic oxygen atom transfer reaction.

Theoretically exploring the key role of the Lys412 residue in the conversion of N2O to N2by nitrous oxide reductase from Achromobacter cycloclastes

The active Cu_Zcluster of NOR provides strong back-donation to coordinated N₂O and activates the O atom of the N₂O group facilitating H-bonding and protonationviathe Lys412 residue.

Why is the molybdenum-substituted tungsten-dependent formaldehyde ferredoxin oxidoreductase not active? A quantum chemical study

Formaldehyde ferredoxin oxidoreductase is a tungsten-dependent enzyme that catalyzes the oxidative degradation of formaldehyde to formic acid. The molybdenum ion can be incorporated into the active site to displace the tungsten ion, but is without activity. Density functional calculations have been employed to understand the incapacitation of the enzyme caused by molybdenum substitution. The calculations show that the enzyme with molybdenum (Mo-FOR) has higher redox potential than that with tungsten, which makes the formation of the Mo(VI)=O complex endothermic by 14 kcal/mol. Following our previously suggested mechanism for this enzyme, the formaldehyde substrate oxidation was also investigated for Mo-FOR using the same quantum-mechanics-only model, except for the displacement of tungsten by molybdenum. The calculations demonstrate that formaldehyde oxidation occurs via a sequential two-step mechanism. Similarly to the tungsten-catalyzed reaction, the Mo(VI)=O species performs the nucleophilic attack on the formaldehyde carbon, followed by proton transfer in concert with two-electron reduction of the metal center. The first step is rate-limiting, with a total barrier of 28.2 kcal/mol. The higher barrier is mainly due to the large energy penalty for the formation of the Mo(VI)=O species.

The prokaryotic Mo/W-bisPGD enzymes family: a catalytic workhorse in bioenergetic

Over the past two decades, prominent importance of molybdenum-containing enzymes in prokaryotes has been put forward by studies originating from different fields. Proteomic or bioinformatic studies underpinned that the list of molybdenum-containing enzymes is far from being complete with to date, more than fifty different enzymes involved in the biogeochemical nitrogen, carbon and sulfur cycles. In particular, the vast majority of prokaryotic molybdenum-containing enzymes belong to the so-called dimethylsulfoxide reductase family. Despite its extraordinary diversity, this family is characterized by the presence of a Mo/W-bis(pyranopterin guanosine dinucleotide) cofactor at the active site. This review highlights what has been learned about the properties of the catalytic site, the modular variation of the structural organization of these enzymes, and their interplay with the isoprenoid quinones. In the last part, this review provides an integrated view of how these enzymes contribute to the bioenergetics of prokaryotes. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

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

A density functional theory study of the influence of the various functional groups of the molybdopterin ligand on electronic and geometric properties of active-site models for the molybdenum and tungsten cofactors has been undertaken. We used analogous molybdenum and tungsten complexes with increasingly accurate representation of the molybdopterin ligands and compared bond lengths, angles, charge distribution, composition of the binding orbitals, as well as the redox potentials in relation to each other. On the basis of our findings, we suggest using ligand systems including the pyrane and the pyrazine rings, besides the dithiolene function, to obtain sufficiently reliable computational, but also synthetic, models for the molybdenum and tungsten cofactors, whereas the second ring of the pterin might be neglected for efficiency reasons.

Periplasmic nitrate reductase revisited: a sulfur atom completes the sixth coordination of the catalytic molybdenum

Nitrate reductase from Desulfovibrio desulfuricans ATCC 27774 (DdNapA) is a monomeric protein of 80 kDa harboring a bis(molybdopterin guanine dinucleotide) active site and a [4Fe-4S] cluster. Previous electron paramagnetic resonance (EPR) studies in both catalytic and inhibiting conditions showed that the molybdenum center has high coordination flexibility when reacted with reducing agents, substrates or inhibitors. As-prepared DdNapA samples, as well as those reacted with substrates and inhibitors, were crystallized and the corresponding structures were solved at resolutions ranging from 1.99 to 2.45 A. The good quality of the diffraction data allowed us to perform a detailed structural study of the active site and, on that basis, the sixth molybdenum ligand, originally proposed to be an OH/OH(2) ligand, was assigned as a sulfur atom after refinement and analysis of the B factors of all the structures. This unexpected result was confirmed by a single-wavelength anomalous diffraction experiment below the iron edge (lambda = 1.77 A) of the as-purified enzyme. Furthermore, for six of the seven datasets, the S-S distance between the sulfur ligand and the Sgamma atom of the molybdenum ligand Cys(A140) was substantially shorter than the van der Waals contact distance and varies between 2.2 and 2.85 A, indicating a partial disulfide bond. Preliminary EPR studies under catalytic conditions showed an EPR signal designated as a turnover signal (g values 1.999, 1.990, 1.982) showing hyperfine structure originating from a nucleus of unknown nature. Spectropotentiometric studies show that reduced methyl viologen, the electron donor used in the catalytic reaction, does not interact directly with the redox cofactors. The turnover signal can be obtained only in the presence of the reaction substrates. With use of the optimized conditions determined by spectropotentiometric titration, the turnover signal was developed with (15)N-labeled nitrate and in D(2)O-exchanged DdNapA samples. These studies indicate that this signal is not associated with a Mo(V)-nitrate adduct and that the hyperfine structure originates from two equivalent solvent-exchangeable protons. The new coordination sphere of molybdenum proposed on the basis of our studies led us to revise the currently accepted reaction mechanism for periplasmic nitrate reductases. Proposals for a new mechanism are discussed taking into account a molybdenum and ligand-based redox chemistry, rather than the currently accepted redox chemistry based solely on the molybdenum atom.