Computational studies of EPR parameters for paramagnetic molybdenum complexes. II. Larger MoV systems relevant to molybdenum enzymes

Exploring Ligand Effects on Structure, Bonding, and Photolytic Hydride Transfer of Cationic Gold(I) Bridging Hydride Complexes of Molybdocene and Tungstenocene

A diverse family of heterobimetallic bridging hydride adducts of the type [LAu(μ-H)2MCp2][X] (L = 1,3-bis(2,6-diisopropylphenyl)imidazole-2-ylidene, IPr; 1,3-bis(1-adamantyl)imidazole-2-ylidene, IAd; 1,3-bis(2,6-di-iso-propylphenyl)-5,5-dimethyl-4,6-diketopyrimidinyl-2-ylidene, DippDAC; triphenylphosphine, PPh3; 2-di-tert-butylphosphino-2',4',6'-triisopropylbiphenyl, tBuXPhos; X = SbF6-, BF4- or TfO-) was synthesized by reacting group VI metallocene dihydrides Cp2MH2 (Cp = cyclopentadienyl anion; M = Mo, W) with cationic gold(I) complexes [LAu(NCMe)][X]. Trimetallic [L'Au2(μ-H)2WCp2][X]2 and tetrametallic [L'Au2{(μ-H)2WCp2}2] [X]2 complexes (L' = rac-2,2'-bis(diphenylphosphino)-1,1'-binaphthalene or bis(diphenylphosphinomethane)) were obtained by reacting digold [L'{Au(NCMe)}2][X]2 with Cp2WH2 in a 1:1 and a 1:2 stoichiometry. Accessing such a broad structural diversity allowed us to pinpoint roles played by the ancillary ligands and group VI metals on the bonding properties of this family of bridging hydrides. In particular, a clear effect of the ligand on the interaction energy and electronic structure was observed, with important implications on photolytic reactivity. UV or visible light irradiation, indeed, leads to the selective cleavage of the heterobimetallic Au(μ-H)2M arrangement and formation of molecular gold hydrides. The photolysis was found to be chromoselective (wavelength-dependent), which can be ascribed to different charge redistributions upon excitation to the first (Kasha's reactivity) and higher (anti-Kasha's reactivity) excited states.

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.

Decoding the Ambiguous Electron Paramagnetic Resonance Signals in the Lytic Polysaccharide Monooxygenase from Photorhabdus luminescens

Understanding the structure and function of lytic polysaccharide monooxygenases (LPMOs), copper enzymes that degrade recalcitrant polysaccharides, requires the reliable atomistic interpretation of electron paramagnetic resonance (EPR) data on the Cu(II) active site. Among various LPMO families, the chitin-active PlAA10 shows an intriguing phenomenology with distinct EPR signals, a major rhombic and a minor axial signal. Here, we combine experimental and computational investigations to uncover the structural identity of these signals. X-band EPR spectra recorded at different pH values demonstrate pH-dependent population inversion: the major rhombic signal at pH 6.5 becomes minor at pH 8.5, where the axial signal dominates. This suggests that a protonation change is involved in the interconversion. Precise structural interpretations are pursued with quantum chemical calculations. Given that accurate calculations of Cu g-tensors remain challenging for quantum chemistry, we first address this problem via a thorough calibration study. This enables us to define a density functional that achieves accurate and reliable prediction of g-tensors, giving confidence in our evaluation of PlAA10 LPMO models. Large models were considered that include all parts of the protein matrix surrounding the Cu site, along with the characteristic second-sphere features of PlAA10. The results uniquely identify the rhombic signal with a five-coordinate Cu ion bearing two water molecules in addition to three N-donor ligands. The axial signal is attributed to a four-coordinate Cu ion where only one of the waters remains bound, as hydroxy. Alternatives that involve decoordination of the histidine brace amino group are unlikely based on energetics and spectroscopy. These results provide a reliable spectroscopy-consistent view on the plasticity of the resting state in PlAA10 LPMO as a foundation for further elucidating structure-property relationships and the formation of catalytically competent species. Our strategy is generally applicable to the study of EPR parameters of mononuclear copper-containing metalloenzymes.

EPR Spectroscopy of Cu(II) Complexes: Prediction of g-Tensors Using Double-Hybrid Density Functional Theory

Computational electron paramagnetic resonance (EPR) spectroscopy is an important field of applied quantum chemistry that contributes greatly to connecting spectroscopic observations with the fundamental description of electronic structure for open-shell molecules. However, not all EPR parameters can be predicted accurately and reliably for all chemical systems. Among transition metal ions, Cu(II) centers in inorganic chemistry and biology, and their associated EPR properties such as hyperfine coupling and g-tensors, pose exceptional difficulties for all levels of quantum chemistry. In the present work, we approach the problem of Cu(II) g-tensor calculations using double-hybrid density functional theory (DHDFT). Using a reference set of 18 structurally and spectroscopically characterized Cu(II) complexes, we evaluate a wide range of modern double-hybrid density functionals (DHDFs) that have not been applied previously to this problem. Our results suggest that the current generation of DHDFs consistently and systematically outperform other computational approaches. The B2GP-PLYP and PBE0-DH functionals are singled out as the best DHDFs on average for the prediction of Cu(II) g-tensors. The performance of the different functionals is discussed and suggestions are made for practical applications and future methodological developments.

Estimating the accuracy of calculated electron paramagnetic resonance hyperfine couplings for a lytic polysaccharide monooxygenase

Lytic polysaccharide monooxygenases (LPMOs) are enzymes that bind polysaccharides followed by an (oxidative) disruption of the polysaccharide surface, thereby boosting depolymerization. The binding process between the LPMO catalytic domain and polysaccharide is key to the mechanism and establishing structure-function relationships for this binding is therefore crucial. The hyperfine coupling constants (HFCs) from EPR spectroscopy have proven useful for this purpose. Unfortunately, EPR does not provide direct structural data and therefore the experimental EPR parameters have to be supported with parameters calculated with density functional theory. Yet, calculated HFCs are extremely sensitive to the employed computational setup. Using the LPMO Ls(AA9)A catalytic domain, we here quantify the importance of several choices in the computational setup, ranging from the use of specialized basis, the underlying structures, and the employed exchange-correlation functional. We show that specialized basis sets are an absolute necessity, and also that care has to be taken in the optimization of the underlying structure: only by allowing large parts of the protein around the active site to structurally relax could we obtain results that uniformly reproduced experimental trends. We compare our results to previously published X-ray structures and experimental HFCs for Ls(AA9)A as well as to recent experimental/theoretical results for another (AA10) family of LPMOs.

A Straightforward Access to Stable, 16 Valence-electron Phosphine-Stabilized Fe0 Olefin Complexes and their Reactivity

The use of the dialkene divinyltetramethyldisiloxane (dvtms) allows easy access to the reactive 16 valence-electron complexes [Fe⁰(L-L)(dvtms)], (L-L) = dppe (1,2-bis(diphenylphosphino)ethane), (1), dppp (1,2-bis(diisopropylphosphino)propane), (2), pyNMeP(ⁱPr)₂ (N-(diisopropylphosphino)-N-methylpyridin-2-amine), (4), dipe (1,2-bis(diisopropylphosphino)ethane), (5), and [Fe⁰(L)₂(dvtms)], L = PMe₃, (3), by a mild reductive route using AlEt₂(OEt) as reducing agent. In contrast, by the same methodology, the 18 valence-electron complexes [Fe⁰(L-L)₂(ethylene)], (L-L) = dppm (1,2-bis(diphenylphosphino)methane), 6, (L-L) = dppa (1,2-bis(diphenylphosphino)amine) 7 or (L-L)=dppe, 8, were obtained, which do not contain dvtms. In addition, a combined DFT and solid-state paramagnetic NMR methodology is introduced for the structure determination of 5. A comparative study of the reactivity of 1,2,4-6 and 8 with 3-hexyne highlights emerging mechanistic implications for C-C coupling reactions using these complexes as catalysts.

d(1) Oxosulfido-Mo(V) Compounds: First Isolation and Unambiguous Characterization of an Extended Series

Reaction of Tp(iPr)Mo(VI)OS(OAr) with cobaltocene in toluene results in the precipitation of brown, microcrystalline oxosulfido-Mo(V) compounds, [CoCp2][Tp(iPr)Mo(V)OS(OAr)] (Cp(-) = η(5)-C5H5(-), Tp(iPr)(-) = hydrotris(3-isopropylpyrazol-1-yl)borate, OAr(-) = phenolate or 2-(s)Bu, 2-(t)Bu, 3-(t)Bu, 4-(s)Bu, 4-Ph, 3,5-(s)Bu2, 2-CO2Me, 2-CO2Et or 2-CO2Ph derivative thereof). The compounds are air- and water-sensitive and display ν(Mo═O) and ν(Mo[Formula: see text]S) IR absorption bands at ca. 890 and 435 cm(-1), respectively, 20-40 cm(-1) lower in energy than the corresponding bands in Tp(iPr)MoOS(OAr). They are electrochemically active and exhibit three reversible cyclovoltammetric waves (E(Mo(VI)/Mo(V)) = -0.40 to -0.66 V, E([CoCp2](+)/CoCp2) = -0.94 V and E(CoCp2/[CoCp2](-)) = -1.88 V vs SCE). Structural characterization of [CoCp2][Tp(iPr)MoOS(OC6H4CO2Et-2)]·2CH2Cl2 revealed a distorted octahedral Mo(V) anion with Mo═O and Mo[Formula: see text]S distances of 1.761(5) and 2.215(2) Å, respectively, longer than corresponding distances in related Tp(iPr)MoOS(OAr) compounds. The observation of strong S(1s) → (S(3p) + Mo(4d)) S K-preedge transitions indicative of a d(1) sulfido-Mo(V) moiety and the presence of short Mo═O (ca. 1.72 Å) and Mo[Formula: see text]S (ca. 2.25 Å) backscattering contributions in the Mo K-edge EXAFS further support the oxosulfido-Mo(V) formulation. The compounds are EPR-active, exhibiting highly anisotropic (Δg 0.124-0.150), rhombic, frozen-glass spectra with g1 close to the value observed for the free electron (ge = 2.0023). Spectroscopic studies are consistent with the presence of a highly covalent Mo[Formula: see text]S π* singly occupied molecular orbital. The compounds are highly reactive, with reactions localized at the terminal sulfido ligand. For example, the compounds react with cyanide and PPh3 to produce thiocyanate and SPPh3, respectively, and various (depending on solvent) oxo-Mo(V) species. Reactions with copper reagents also generally lead to desulfurization and the formation of oxo-Mo(V) or -Mo(IV) complexes.

The world as viewed by and with unpaired electrons

Recent advances in electron paramagnetic resonance (EPR) include capabilities for applications to areas as diverse as archeology, beer shelf life, biological structure, dosimetry, in vivo imaging, molecular magnets, and quantum computing. Enabling technologies include multifrequency continuous wave, pulsed, and rapid scan EPR. Interpretation is enhanced by increasingly powerful computational models.

Structural studies of the molybdenum center of mitochondrial amidoxime reducing component (mARC) by pulsed EPR spectroscopy and 17O-labeling

Mitochondrial amidoxime reducing components (mARC-1 and mARC-2) represent a novel group of Mo-containing enzymes in eukaryotes. These proteins form the catalytic part of a three-component enzyme complex known to be responsible for the reductive activation of several N-hydroxylated prodrugs. No X-ray crystal structures are available for these enzymes as yet. A previous biochemical investigation [Wahl, B., et al. (2010) J. Biol. Chem., 285, 37847-37859 ] has revealed that two of the Mo coordination positions are occupied by sulfur atoms from a pyranopterindithiolate (molybdopterin, MPT) cofactor. In this work, we have used continuous wave and pulsed electron paramagnetic resonance (EPR) spectroscopy and density functional theoretical (DFT) calculations to determine the nature of remaining ligands in the Mo(V) state of the active site of mARC-2. Experiments with samples in D(2)O have identified the exchangeable equatorial ligand as a hydroxyl group. Experiments on samples in H(2)(17)O-enriched buffer have shown the presence of a slowly exchangeable axial oxo ligand. Comparison of the experimental (1)H and (17)O hyperfine interactions with those calculated using DFT has shown that the remaining nonexchangeable equatorial ligand is, most likely, protein-derived and that the possibility of an equatorial oxo ligand can be excluded.

EPR parameters of radical ions and polarizability effect

The literature data on substituent influence on the g factors, hyperfine coupling (hfc) constants (a) of the EPR signals and the spin densities (ρ) have been analyzed for 25 series of the organic radical cations and radical anions as well as of the transition metal complexes. The g, a, and ρ values were first established to depend not only on the inductive and resonance effects but also on the polarizability of substituents which can be characterized by the σ(α) constants. The polarizability effect is caused by the partial charge on the radical center. This effect consists in an electrostatic attraction between the charge and the dipole moments induced by this charge in the substituents. The polarizability contribution ranges up to 92%.

Comparative calculation of EPR spectral parameters in [Mo(V)OX4]-, [Mo(V)OX5]2-, and [Mo(V)OX4(H2O)]- complexes

The EPR spectral parameters, i.e. g-tensors and molybdenum hyperfine couplings, for several d(1) systems of the general formula [Mo(V)EX(4)](n-), [Mo(V)OX(5)](2-), and [Mo(V)OX(4)(H(2)O)](-) (E = O, N; X = F, Cl, Br; n = 1 or 2) were calculated using Density Functional Theory. The influence of basis sets, their contraction scheme, the type of exchange-correlation functional, the amount of Hartree-Fock exchange, molecular geometry, and relativistic effects on the calculated EPR spectra parameters have been discussed. The g-tensors and molybdenum hyperfine coupling parameters were calculated using a relativistic Hamiltonian coupled with several LDA, GGA, and 'hybrid' exchange-correlation functionals and uncontracted full-electron DGauss DZVP basis sets. The calculated EPR parameters are found to be sensitive to the Mo=E distance and E=Mo-Cl angle, and thus the choice of starting molecular geometry should be considered as an important factor in predicting the g-tensors and hyperfine coupling constants in oxo-molybdenum compounds. In the present case, the GGA exchange-correlation functionals provide a better agreement between the theory and the experiment.

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.