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.

Intermediate Detection in the Casiopeina-Cysteine Interaction Ending in the Disulfide Bond Formation and Copper Reduction

A strategy to improve the cancer therapies involves agents that cause the depletion of the endogenous antioxidant glutathione (GSH), increasing its efflux out of cells and inducing apoptosis in tumoral cells due to the presence of reactive oxygen species. It has been shown that Casiopeina copper complexes caused a dramatic intracellular GSH drop, forming disulfide bonds and reducing Cu^II to Cu^I. Herein, through the determination of the [Cu^II]-SH bond before reduction, we present evidence of the adduct between cysteine and one Casiopeina as an intermediate in the cystine formation and as a model to understand the anticancer activity of copper complexes. Evidence of such an intermediate has never been presented before.

Comparison of Density Functional and Correlated Wave Function Methods for the Prediction of Cu(II) Hyperfine Coupling Constants

The reliable prediction of Cu(II) hyperfine coupling constants remains a challenge for quantum chemistry. Until recently only density functional theory (DFT) could target this property for systems of realistic size. However, wave function based methods become increasingly applicable. In the present work, we define a large set of Cu(II) complexes with experimentally known hyperfine coupling constants and use it to investigate the performance of modern quantum chemical methods for the prediction of this challenging spectroscopic parameter. DFT methods are evaluated against orbital‐optimized second‐order Møller‐Plesset (OO‐MP2) theory and coupled cluster calculations including singles and doubles excitations, driven by the domain‐based local pair natural orbital approach (DLPNO‐CCSD). Special attention is paid to the definition of a basis set that converges adequately toward the basis set limit for the given property for all methods considered in this study, and a specifically optimized basis set is proposed for this purpose. The results suggest that wave function based methods can supplant but do not outcompete DFT for the calculation of Cu(II) hyperfine coupling constants. Mainstream hybrid functionals such as B3PW91 remain on average the best choice.

Unveiling the Role of Hydrogen Bonding and g-Tensor in the Interaction of Ru-Bis-DMSO with Amino Acid Residue and Human Serum Albumin

Density functional theory calculations have been carried out to observe the role of hydrogen bonding in hydrolysis and the coordination mechanism of three amino acid residues (histidine, cysteine, and alanine) with Ru-bis-DMSO complex via which the complex tends to interact with the HSA protein receptor. The interaction mechanism shows that ruthenium complexes prefer to bind protein receptor through cysteine and histidine residues rather than through alanine, which has been confirmed by DFT evaluated H-bonding and g-tensor analysis. The number of H-bonds plays a major role in stabilizing the intermediates and transition states involved in the Ru-bis-DMSO and amino acid residue interactions. Our theoretical g-tensor values are in good agreement with the available experimental results. Further QM/MM calculation on the Ru-bis-DMSO-HSA adducts reveals that the adduct is more stable when Ru gets coordinated with histidine imidazole rather than cysteine. These investigations helped us in understanding the type of amino acid residue responsible for binding the metal complex Ru-bis-DMSO with the carrier protein HSA.

Correlation between the C-C Cross-Coupling Activity and C-to-Ni Charge Transfer Transition of High-Valent Ni Complexes

High-valent Ni complexes have proven to be good platforms for diverse cross-coupling reactions that are otherwise difficult to be achieved with conventional low-valent catalysts. However, their reductive elimination (RE) activities are still significantly variable by up to 5 orders of magnitude, depending on the supporting ligand and oxidation state of the Ni center. To elucidate frontier molecular orbitals (FMOs) that determine the RE activity of the Ni center, the electronic structures of cycloneophyl (CH₂C(CH₃)₂-o-C₆H₄) Ni^III and Ni^IV complexes have been characterized by utilizing various transition metal-based spectroscopic techniques such as electronic absorption, magnetic circular dichroism, electron paramagnetic resonance, resonance Raman, and X-ray absorption spectroscopies. In combination with density functional theory computations, the spectroscopic analyses have shown that the energies of the C-to-Ni charge-transfer (CT) electronic transitions are strongly correlated to the rates of C-C bond-forming RE reaction. This correlation suggests that the kinetic barrier of the RE reaction is determined by energy cost for internal CT (ICT) from the coordinated carbon moiety to the Ni center, and that FMOs involved in the RE reaction and the C-to-Ni CT electronic transitions are essentially identical. This FMO determination has led us to discover that photoexcitation to the C-to-Ni CT excited states accelerates the C-C cross-coupling reaction by up to 10⁵ times, as the CT electronic transition can substitute for the rate-determining ICT step of the RE reaction at the ground electronic state.

DFT Protocol for EPR Prediction of Paramagnetic Cu(II) Complexes and Application to Protein Binding Sites

With the aim to provide a general protocol to interpret electron paramagnetic resonance (EPR) spectra of paramagnetic copper(II) coordination compounds, density functional theory (DFT) calculations of spin Hamiltonian parameters g and A for fourteen Cu(II) complexes with different charges, donor sets, and geometry were carried out using ORCA software. The performance of eleven functionals was tested, and on the basis of the mean absolute percent deviation (MAPD) and standard deviation (SD), the ranking of the functionals for Az is: B3LYP > B3PW91 ~ B3P86 > PBE0 > CAM-B3LYP > TPSSh > BH and HLYP > B2PLYP > MPW1PW91 > ω-B97x-D >> M06; and for gz is: PBE0 > BH and HLYP > B2PLYP > ω-B97x-D > B3PW91~B3LYP~B3P86 > CAM-B3LYP > TPSSh~MPW1PW91 >> M06. With B3LYP the MAPD with respect to A z exp t l is 8.6% with a SD of 4.2%, while with PBE0 the MAPD with respect to g z exp t l is 2.9% with a SD of 1.1%. The results of the validation confirm the fundamental role of the second order spin-orbit contribution to Az. The computational procedure was applied to predict the values of gz and Az of the adducts formed by Cu(II) with albumin and two fragments of prion protein, 106–126 and 180–193.

Copper(II)-bis-histidine coordination structure in a fibrillar amyloid β-peptide fragment and model complexes revealed by electron spin echo envelope modulation spectroscopy

Truncated and mutated amyloid-β (Aβ) peptides are models for systematic study-in homogeneous preparations-of the molecular origins of metal ion effects on Aβ aggregation rates, types of aggregate structures formed, and cytotoxicity. The 3D geometry of bis-histidine imidazole coordination of Cu(II) in fibrils of the nonapetide acetyl-Aβ(13-21)H14A has been determined by powder (14) N electron spin echo envelope modulation (ESEEM) spectroscopy. The method of simulation of the anisotropic combination modulation is described and benchmarked for a Cu(II) -bis-cis-imidazole complex of known structure. The revealed bis-cis coordination mode, and the mutual orientation of the imidazole rings, for Cu(II) in Ac-Aβ(13-21)H14A fibrils are consistent with the proposed β-sheet structural model and pairwise peptide interaction with Cu(II) , with an alternating [-metal-vacancy-]n pattern, along the N-terminal edge. Metal coordination does not significantly distort the intra-β-strand peptide interactions, which provides a possible explanation for the acceleration of Ac-Aβ(13-21)H14A fibrillization by Cu(II) , through stabilization of the associated state and low-reorganization integration of β-strand peptide pair precursors.

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.

Application of DFT methods to the study of the coordination environment of the VO2+ ion in V proteins

Density functional theory (DFT) methods were used to simulate the environment of vanadium in several V proteins, such as vanadyl-substituted carboxypeptidase (sites A and B), vanadyl-substituted chloroplast F(1)-ATPase (CF(1); site 3), the reduced inactive form of vanadium bromoperoxidase (VBrPO; low- and high-pH sites), and vanadyl-substituted imidazole glycerol phosphate dehydratase (IGPD; sites α, β, and γ). Structural, electron paramagnetic resonance, and electron spin echo envelope modulation parameters were calculated and compared with the experimental values. All the simulations were performed in water within the framework of the polarizable continuum model. The angular dependence of [Formula: see text] and [Formula: see text] on the dihedral angle θ between the V=O and N-C bonds and on the angle φ between the V=O and V-N bonds, where N is the coordinated aromatic nitrogen atom, was also found. From the results it emerges that it is possible to model the active site of a vanadium protein through DFT methods and determine its structure through the comparison between the calculated and experimental spectroscopic parameters. The calculations confirm that the donor sets of sites B and A of vanadyl-substituted carboxypeptidase are [[Formula: see text], H(2)O, H(2)O, H(2)O] and [N(His)(||), N(His)(⊥), [Formula: see text], H(2)O], and that the donor set of site 3 of CF(1)-ATPase is [[Formula: see text], OH(Thr), H(2)O, H(2)O, [Formula: see text]]. For VBrPO, the coordination modes [N(His)(||), N(His)(∠), OH(Ser), H(2)O, H(2)O(ax)] for the low-pH site and [N(His)(||), N(His)(∠), OH(Ser), OH(-), H(2)O(ax)] or [N(His)(||), N(His)(∠), [Formula: see text], H(2)O] for the high-pH site, with an imidazole ring of histidine strongly displaced from the equatorial plane, can be proposed. Finally, for sites α, β, and γ of IGPD, the subsequent deprotonation of one, two, and three imidazole rings of histidine and the participation of a carboxylate group of a glutamate residue ([N(His)(||), [Formula: see text], H(2)O, H(2)O], [N(His)(||), N(His)(||), [Formula: see text], H(2)O], and [N(His)(||), N(His)(||), [Formula: see text], OH(-), [Formula: see text]], respectively) seems to be the most plausible hypothesis.