A comparison of the thermodynamics of O[bond]O cleavage for dicopper complexes in enzymes and synthetic systems

Hydroxypyridinone-Based Metal Chelators towards Ecotoxicity: Remediation and Biological Mechanisms

Hydroxypyridinones (HPs) are recognized as excellent chemical tools for engineering a diversity of metal chelating agents, with high affinity for hard metal ions, exhibiting a broad range of activities and applications, namely in medical, biological and environmental contexts. They are easily made and functionalizable towards the tuning of their pharmacokinetic properties or the improving of their metal complex thermodynamic stabilities. In this review, an analysis of the recently published works on hydroxypyridinone-based ligands, that have been mostly addressed for environmental applications, namely for remediation of hard metal ion ecotoxicity in living beings and other biological matrices is carried out. In particular, herein the most recent developments in the design of new chelating systems, from bidentate mono-HP to polydentate multi-HP derivatives, with a structural diversity of soluble or solid-supported backbones are outlined. Along with the ligand design, an analysis of the relationship between their structures and activities is presented and discussed, namely associated with the metal affinity and the thermodynamic stability of the corresponding metal complexes.

Mechanistic Insights into the ortho-Defluorination-Hydroxylation of 2-Halophenolates Promoted by a Bis(μ-oxo)dicopper(III) Complex

C-F bonds are one of the most inert functionalities. Nevertheless, some [Cu₂O₂]²⁺ species are able to defluorinate-hydroxylate ortho-fluorophenolates in a chemoselective manner over other ortho-halophenolates. Albeit it is known that such reactivity is promoted by an electrophilic attack of a [Cu₂O₂]²⁺ core over the arene ring, the crucial details of the mechanism that explain the chemo- and regioselectivity of the reaction remain unknown, and it has not being determined either if Cu^II₂(η²:η²-O₂) or Cu^III₂(μ-O)₂ species are responsible for the initial attack on the arene. Herein, we present a combined theoretical and experimental mechanistic study to unravel the origin of the chemoselectivity of the ortho-defluorination-hydroxylation of 2-halophenolates by the [Cu₂(O)₂(DBED)₂]²⁺ complex (DBED = N,N'-di-tert-butylethylenediamine). Our results show that the equilibria between (side-on)peroxo (P) and bis(μ-oxo) (O) isomers plays a key role in the mechanism, with the latter being the reactive species. Furthermore, on the basis of quantum-mechanical calculations, we were able to rationalize the chemoselective preference of the [Cu₂(O)₂(DBED)₂]²⁺ catalyst for the C-F activation over C-Cl and C-H activations.

Theoretical rationalization for the equilibrium between (μ-η2:η2-peroxido)CuIICuII and bis(μ-oxido)CuIIICuIII complexes: perturbational effects from ligand frameworks

DFT calculations are carried out to investigate the geometric effects of the supporting ligands in the relative energies of the (μ-η2:η2-peroxido)CuIICuII complex 1 and the bis(μ-oxido)CuIIICuIII complex 2. The N3-tridentate ligand bearing acyclic propane diamine framework La preferentially provided 1, whereas the N3-tridentate ligand with cyclic diamine framework such as 1,4-diazacycloheptane Lb gave 2 after the oxygenation of the corresponding CuI complexes as reported previously [S. Itoh, et al., Inorg. Chem., 2014, 53, 8786-8794]. Calculations at the B3LYP*-D3 level of theory can reasonably explain the experimental results in relative energies, structures and harmonic frequencies of 1 and 2. Perturbational effects of the diamine chelates of La and Lb especially on the equilibrium of 1 and 2 are investigated in detail. In the range from 2.30 Å to 3.40 Å of the N-N distance in the diamine moiety, 1 is more stable than 2 by 8.4 kcal mol-1 at the distance of 3.40 Å. Calculated potential energies indicate that the decrease in the N-N distance is associated with a decrease in energy of 2, leading that 2 can be most stabilized at the N-N distance of 2.60 Å. Furthermore, molecular orbitals analyses are performed to explain that the energy gaps between the σ* orbital of the O-O bond and the dx2-y2 orbitals of the CuII ions of 1 get small as the diamine moiety is shrunk, leading to facilitate the O-O bond cleavage from 1 to 2.

Mechanistic analysis of water oxidation catalyzed by mononuclear copper in aqueous bicarbonate solutions

We characterize a mechanism for a monomeric copper catalyst reported to oxidize water in bicarbonate solution when subject to sufficiently high external potentials at near neutral pH values.

Electrostatic effects on proton coupled electron transfer in oxomanganese complexes inspired by the oxygen-evolving complex of photosystem II

The influence of electrostatic interactions on the free energy of proton coupled electron transfer in biomimetic oxomanganese complexes inspired by the oxygen-evolving complex (OEC) of photosystem II (PSII) are investigated. The reported study introduces an enhanced multiconformer continuum electrostatics (MCCE) model, parametrized at the density functional theory (DFT) level with a classical valence model for the oxomanganese core. The calculated pKa's and oxidation midpoint potentials (E(m)'s) match experimental values for eight complexes, indicating that purely electrostatic contributions account for most of the observed couplings between deprotonation and oxidation state transitions. We focus on pKa's of terminal water ligands in [Mn(II/III)(H2O)6](2+/3+) (1), [Mn(III)(P)(H2O)2](3-) (2, P = 5,10,15,20-tetrakis(2,6-dichloro-3-sulfonatophenyl)porphyrinato), [Mn2(IV,IV)(μ-O)2(terpy)2(H2O)2](4+) (3, terpy = 2,2':6',2″-terpyridine), and [Mn3(IV,IV,IV)(μ-O)4(phen)4(H2O)2](4+) (4, phen = 1,10-phenanthroline) and the pKa's of μ-oxo bridges and Mn E(m)'s in [Mn2(μ-O)2(bpy)4] (5, bpy = 2,2'-bipyridyl), [Mn2(μ-O)2(salpn)2] (6, salpn = N,N'-bis(salicylidene)-1,3-propanediamine), [Mn2(μ-O)2(3,5-di(Cl)-salpn)2] (7), and [Mn2(μ-O)2(3,5-di(NO2)-salpn)2] (8). The analysis of complexes 6-8 highlights the strong coupling between electron and proton transfers, with any Mn oxidation lowering the pKa of an oxo bridge by 10.5 ± 0.9 pH units. The model also accounts for changes in the E(m)'s by ligand substituents, such as found in complexes 6-8, due to the electron withdrawing Cl (7) and NO2 (8). The reported study provides the foundation for analysis of electrostatic effects in other oxomanganese complexes and metalloenzymes, where proton coupled electron transfer plays a fundamental role in redox-leveling mechanisms.

Complete σ* intramolecular aromatic hydroxylation mechanism through O2 activation by a Schiff base macrocyclic dicopper(I) complex

In this work we analyze the whole molecular mechanism for intramolecular aromatic hydroxylation through O2 activation by a Schiff hexaazamacrocyclic dicopper(I) complex, [Cu(I) 2(bsH2m)](2+). Assisted by DFT calculations, we unravel the reaction pathway for the overall intramolecular aromatic hydroxylation, i.e., from the initial O2 reaction with the dicopper(I) species to first form a Cu(I)Cu(II)-superoxo species, the subsequent reaction with the second Cu(I) center to form a μ-η(2):η(2)-peroxo-Cu(II) 2 intermediate, the concerted peroxide O-O bond cleavage and C-O bond formation, followed finally by a proton transfer to an alpha aromatic carbon that immediately yields the product [Cu(II) 2(bsH2m-O)(μ-OH)](2+).

Dioxygen reactivity of new bispidine-copper complexes

The reactivity of copper complexes of three different second-generation bispidine-based ligands (bispidine = 3,7-diazabicyclo[3.3.1]nonane; mono- and bis-tetradentate; exclusively tertiary amine donors) with dioxygen [(reversible) binding of dioxygen by copper(I)] is reported. The UV-vis, electrospray ionization mass spectrometry, electron paramagnetic resonance, and vibrational spectra (resonance Raman) of the dioxygen adducts indicate that, depending on the ligand and reaction conditions, several different species (mono- and dinuclear, superoxo, peroxo, and hydroperoxo), partially in equilibrium with each other, are formed. Minor changes in the ligand structure and/or experimental conditions (solvent, temperature, relative concentrations) allow switching between the different forms. With one of the ligands, an end-on peroxodicopper(II) complex and a mononuclear hydroperoxocopper(II) complex could be characterized. With another ligand, reversible dioxygen binding was observed, leading to a metastable superoxocopper(II) complex. The amount of dioxygen involved in the reversible binding to Cu(I) was determined quantitatively. The mechanism of dioxygen binding as well as the preference of each of the three ligands for a particular dioxygen adduct is discussed on the basis of a computational (density functional theory) analysis.

Correlated wavefunction methods in bioinorganic chemistry

In this commentary the challenges faced in the application of wavefunction-based ab initio methods to (open-shell) transition metal complexes of (bio)inorganic interest are briefly touched on. Both single-reference and multireference methods are covered. It is stressed that the generation and nature of the reference wavefunction is a subject of major importance. How erroneous results can be easily obtained even with coupled-cluster theory is illustrated through the example of the septet-quintet separation in iron(IV)-oxo complexes. Second, the interplay between relativistic and correlation effects is important. This is demonstrated with coupled-cluster calculations on models for dinuclear copper active sites, where relativity has a major influence on the relative stabilities of the bis(μ-oxo) and side-on peroxo species.

Mechanistic insight into the alcohol oxidation mediated by an efficient green [CuBr(2)(2,2'-bipy)]-TEMPO catalyst by density functional method

Density functional theory (DFT) calculations have been performed to investigate the alcohol oxidation to acetaldehyde catalyzed by [CuBr(2)(2,2'-bipy)]-TEMPO (TEMPO stands for 2,2,6,6-tetramethylpiperidinyloxy; bipy stands for bipyridine). The total charge for the studied catalytic system is +1. The catalytic cycle consists of two parts, namely, alcohol oxidation and TEMPO regeneration. In alcohol oxidation, the reaction follows the Sheldon's mechanism for the proposed two mechanisms, i.e., Semmelhack's mechanism and Sheldon's mechanism. The water participation plays minor role in the H atom abstraction step. In TEMPO regeneration, the proposed three paths are competitive in energy. By comparing with experimental observation, it is found that the path, in which alcohol provides the proton to TEMPO(-) to produce TEMPOH followed by the oxidation of TEMPOH directly to TEMPO by O(2), is favored. In TEMPO regeneration, CH(3)CN acts as the ligand to stabilize the Cu(I) species during the catalytic cycle.

Probing the mechanism of O2 activation by a copper(I) biomimetic complex of a C-H hydroxylating copper monooxygenase

In this paper, we report, for the first time, a plausible full reaction pathway for the activation of O(2) by a tetraazamacrocyclic monocopper(I) complex and for the subsequent intramolecular alkylic hydroxylation to yield the alkoxide product. This theoretical insight offers remarkable support to the fundamental hypothesis in the field that a hydroperoxo complex of the type Cu(II)OOH intermediate is the key intermediate in this class of reactions. Overall, we give insight into an intramolecular alkylic C-H bond activation due to the O(2) binding to copper(I) with an end-on eta(1)-O(2) ligation. The loss of a water molecule involves the final substrate oxygenation. The complex we consider is a biomimetic of several systems of biological relevance, such as amine oxidases, peptidylglycine-alpha-hydroxylating monooxygenase, and dopamine-beta monooxygenases.

Theoretical study of the hydroxylation of phenolates by the Cu(2)O (2)(N,N'-dimethylethylenediamine) (2) (2+) complex

Tyrosinase catalyzes the ortho hydroxylation of monophenols and the subsequent oxidation of the diphenolic products to the resulting quinones. In efforts to create biomimetic copper complexes that can oxidize C-H bonds, Stack and coworkers recently reported a synthetic mu-eta(2):eta(2)-peroxodicopper(II)(DBED)(2) complex (DBED is N,N'-di-tert-butylethylenediamine), which rapidly hydroxylates phenolates. A reactive intermediate consistent with a bis-mu-oxo-dicopper(III)-phenolate complex, with the O-O bond fully cleaved, is observed experimentally. Overall, the evidence for sequential O-O bond cleavage and C-O bond formation in this synthetic complex suggests an alternative mechanism to the concerted or late-stage O-O bond scission generally accepted for the phenol hydroxylation reaction performed by tyrosinase. In this work, the reaction mechanism of this peroxodicopper(II) complex was studied with hybrid density functional methods by replacing DBED in the mu-eta(2):eta(2)-peroxodicopper(II)(DBED)(2) complex by N,N'-dimethylethylenediamine ligands to reduce the computational costs. The reaction mechanism obtained is compared with the existing proposals for the catalytic ortho hydroxylation of monophenol and the subsequent oxidation of the diphenolic product to the resulting quinone with the aim of gaining some understanding about the copper-promoted oxidation processes mediated by 2:1 Cu(I)O(2)-derived species.

The restricted active space followed by second-order perturbation theory method: theory and application to the study of CuO2 and Cu2O2 systems

A multireference second-order perturbation theory using a restricted active space self-consistent field wave function as reference (RASPT2/RASSCF) is described. This model is particularly effective for cases where a chemical system requires a balanced orbital active space that is too large to be addressed by the complete active space self-consistent field model with or without second-order perturbation theory (CASPT2 or CASSCF, respectively). Rather than permitting all possible electronic configurations of the electrons in the active space to appear in the reference wave function, certain orbitals are sequestered into two subspaces that permit a maximum number of occupations or holes, respectively, in any given configuration, thereby reducing the total number of possible configurations. Subsequent second-order perturbation theory captures additional dynamical correlation effects. Applications of the theory to the electronic structure of complexes involved in the activation of molecular oxygen by mono- and binuclear copper complexes are presented. In the mononuclear case, RASPT2 and CASPT2 provide very similar results. In the binuclear cases, however, only RASPT2 proves quantitatively useful, owing to the very large size of the necessary active space.

Validation of density functional modeling protocols on experimental bis(μ-oxo)/μ-η2:η2-peroxo dicopper equilibria

The bis(mu-oxo)/mu-eta(2):eta(2)-peroxo equilibria for seven supported Cu(2)O(2) cores were studied with different hybrid and nonhybrid density functional theory models, namely, BLYP, mPWPW, TPSS, TPSSh, B3LYP, mPW1PW, and MPW1K. Supporting ligands 3,3'-iminobis(N,N-dimethylpropylamine), N,N,N',N',N''-pentamethyldipropylenetriamine, N-[2-(pyridin-2-yl)ethyl]-N,N,N'-trimethylpropane-1,3-diamine, bis[2-(2-pyridin-2-yl)ethyl]methylamine, bis[2-(4-methoxy-2-pyridin-2-yl)ethyl]methylamine, bis[2-(4-N,N-dimethylamino-2-pyridin-2-yl)ethyl]methylamine, and 1,4,7-triisopropyl-1,4,7-triazacyclononane were chosen on the basis of the availability of experimental data for comparison. Density functionals were examined with respect to their ability accurately to reproduce experimental properties, including, in particular, geometries and relative energies for the bis(mu-oxo) and side-on peroxo forms. While geometries from both hybrid and nonhybrid functionals were in good agreement with experiment, the incorporation of Hartree-Fock (HF) exchange in hybrid density functionals was found to have a large, degrading effect on predicted relative isomer energies. Specifically, hybrid functionals predicted the mu-eta(2):eta(2)-peroxo isomer to be too stable by roughly 5-10 kcal mol(-1) for each 10% of HF exchange incorporated into the model. Continuum solvation calculations predict electrostatic effects to favor bis(mu-oxo) isomers by 1-4 kcal mol(-1) depending on ligand size, with larger ligands having smaller differential solvation effects. Analysis of computed molecular partition functions suggests that nonzero measured entropies of isomerization are likely to be primarily associated with interactions between molecular solutes and their first solvation shell.

Theoretical Characterization of End-On and Side-On Peroxide Coordination in Ligated Cu2O2 Models

The relative energetics of mu-eta1:eta1 (trans end-on) and mu-eta2:eta2 (side-on) peroxo isomers of Cu2O2 fragments supported by 0, 2, 4, and 6 ammonia ligands have been computed with various density functional, coupled-cluster, and multiconfigurational protocols. There is substantial disagreement between the different levels for most cases, although completely renormalized coupled-cluster methods appear to offer the most reliable predictions. The significant biradical character of the end-on peroxo isomer proves problematic for the density functionals, while the demands on active space size and the need to account for interactions between different states in second-order perturbation theory prove challenging for the multireference treatments. In the latter case, it proved impossible to achieve any convincing convergence.

The performance of hybrid DFT for mechanisms involving transition metal complexes in enzymes

The accuracy of density functional theory with the B3LYP functional is reviewed for systems of relevance to transition-metal-containing enzymes. Calculated energies are commonly within 3-5 kcal/mol of the correct values; however, some exceptions have appeared in the literature and are discussed here. For example, the binding of NO and that of O(2) to metal centers have for some time been known to be underestimated. Most barriers for chemical reactions are overestimated except those involving hydrogen (or proton) transfer, which instead tend to be underestimated. A minor general improvement of the accuracy can probably be obtained by slightly reducing the amount of exact exchange in the B3LYP functional.

Theoretical studies of enzyme mechanisms involving high-valent iron intermediates

Recent theoretical contributions to the elucidation of mechanisms for iron containing enzymes are reviewed. The method used in most of these studies is hybrid density functional theory with the B3LYP functional. Three classes of enzymes are considered, the mononuclear non-heme enzymes, enzymes containing iron dimers, and heme-containing enzymes. Mechanisms for both dioxygen and substrate activations are discussed. The reactions usually go through two half-cycles, where a high-valent intermediate Fe(IV)O species is created in the first half-cycle, and the substrate reactions involving this intermediate occur in the second half-cycle. Similarities between the three classes of enzymes dominate, but significant differences also exist.

Effects of thioether substituents on the O2 reactivity of beta-diketiminate-Cu(I) complexes: probing the role of the methionine ligand in copper monooxygenases

The activation of dioxygen by dopamine beta-monooxygenase (DbetaM) and peptidylglycine alpha-hydroxylating monooxygenase (PHM) is postulated to occur at a copper site ligated by two histidine imidazoles and a methionine thioether, which is unusual because such thioether ligation is not present in other O2-activating copper proteins. To assess the possible role of the thioether ligand in O2 activation by DbetaM and PHM, two new ligands comprising beta-diketiminates with thioether substituents were synthesized and Cu(I) and Cu(II) complexes were isolated. The Cu(II) compounds are monomeric and exhibit intramolecular thioether coordination. While the Cu(I) complexes exhibit a multinuclear topology in the solid state, variable-temperature 1H NMR studies implicate equilibria in solution, possibly including monomers with intramolecular thioether coordination that are structurally defined by DFT calculations. Low-temperature oxygenation of solutions of the Cu(I) complexes generates stable 1:1 Cu/O2 adducts, which on the basis of combined experimental and theoretical studies adopt side-on "eta(2)" structures with negligible Cu-thioether bonding and significant peroxo-Cu(III) character. In contrast to previously reported findings with related ligands lacking the thioether group, however (cf., Aboelella; et al. J. Am. Chem. Soc. 2004, 126, 16896), purging the solutions of the thioether-containing adducts with argon results in conversion to bis(mu-oxo)dicopper(III) species. A role for the thioether in promoting loss of O2 from the 1:1 Cu/O2 adduct and facilitating trapping of the resulting Cu(I) complex to yield the bis(mu-oxo) species is proposed, and the possible relevance of this role to that of the methionine in the active sites of DbetaM and PHM is discussed.

Theoretical Models on the Cu2O2 Torture Track: Mechanistic Implications for Oxytyrosinase and Small-Molecule Analogues

Accurately describing the relative energetics of alternative bis(mu-oxo) and mu-eta2:eta2 peroxo isomers of Cu2O2 cores supported by 0, 2, 4, and 6 ammonia ligands is remarkably challenging for a wide variety of theoretical models, primarily owing to the difficulty of maintaining a balanced description of rapidly changing dynamical and nondynamical electron correlation effects and a varying degree of biradical character along the isomerization coordinate. The completely renormalized coupled-cluster level of theory including triple excitations and extremely efficient pure density functional levels of theory quantitatively agree with one another and also agree qualitatively with experimental results for Cu2O2 cores supported by analogous but larger ligands. Standard coupled-cluster methods, such as CCSD(T), are in most cases considerably less accurate and exhibit poor convergence in predicted relative energies. Hybrid density functionals significantly underestimate the stability of the bis(mu-oxo) form, with the magnitude of the error being directly proportional to the percentage Hartree-Fock exchange in the functional. Single-root CASPT2 multireference second-order perturbation theory, by contrast, significantly overestimates the stability of bis(mu-oxo) isomers. Implications of these results for modeling the mechanism of C-H bond activation by supported Cu2O2 cores, like that found in the active site of oxytyrosinase, are discussed.

4,4'-Dihydroxybiphenyl as a new potent tyrosinase inhibitor

The color of mammalian skin is determined by many factors, for which visible ones are the degree and distribution of melanin pigmentation. Because tyrosinase, (polyphenol oxidase) is the key enzyme for melanin biosynthesis, the use of various tyrosinase inhibitors is a common practice for whitening purpose in cosmetics. In the present study, the inhibition of tyrosinase by 4,4'-dihydroxybiphenyl (44'-BP) was investigated. In addition to tyrosinase inhibiting activity, melanin biosynthesis was assessed in B16F10 melanoma cells (B16 cells). The results showed that 44'-BP exhibits a strong anti-tyrosinase activity with IC50=1.91 microM. The kinetic analysis of tyrosinase inhibition revealed that 44'-BP acts a competitive inhibitor (Ki=4.0 x 10(-4) M at 2.5 microM and Ki =21 x 10(-5) M at 5 microM). Furthermore, data on melanin biosynthesis indicated that the amount of melanin was clearly suppressed by 44'-BP.

The catalytic cycle of catechol oxidase

Hybrid density functional theory with the B3LYP functional has been used to investigate the catalytic mechanism of catechol oxidase. Catechol oxidase belongs to a class of enzymes that has a copper dimer with histidine ligands at the active site. Another member of this class is tyrosinase, which has been studied by similar methods previously. An important advantage for the present study compared to the one for tyrosinase is that X-ray crystal structures exist for catechol oxidase. The most critical step in the mechanism for catechol oxidase is where the peroxide O-O bond is cleaved. In the suggested mechanism this cleavage occurs in concert with a proton transfer from the substrate. Shortly after the transition state is passed there is another proton transfer from the substrate, which completes the formation of a water molecule. An important feature of the mechanism, like the one for tyrosinase, is that no proton transfers to or from residues outside the metal complex are needed. The calculated energetics is in reasonable agreement with experiments. Comparisons are made to other similar enzymes studied previously.