Article pubs.acs.org/IC
Relation between the Catalytic Efficiency of the Synthetic Analogues of Catechol Oxidase with Their Electrochemical Property in the Free State and Substrate-Bound State Prateeti Chakraborty,† Jaydeep Adhikary,† Bipinbihari Ghosh,∥ Ria Sanyal,† Shyamal Kumar Chattopadhyay,*,∥ Antonio Bauzá,‡ Antonio Frontera,*,‡ Ennio Zangrando,*,§ and Debasis Das*,† †
Department of Chemistry, University of Calcutta, 92 A. P. C. Road, Kolkata-700 009, India Departament de Química, Universitat de les Illes Balears, Crta. De Valldemossa km 7.5, 07122 Palma, Baleares, Spain § Dipartimento di Scienze Chimiche, University of Trieste, Via L. Giorgieri 1, 34127 Trieste, Italy ∥ Department of Chemistry, Bengal Engineering and Science University, Howrah, India ‡
S Supporting Information *
ABSTRACT: A library of 15 dicopper complexes as synthetic analogues of catechol oxidase has been synthesized with the aim to determine the relationship between the electrochemical behavior of the dicopper(II) species in the absence as well as in the presence of 3,5-di-tert-butylcatechol (3,5-DTBC) as model substrate and the catalytic activity, kcat, in DMSO medium. The complexes have been characterized by routine physicochemical techniques as well as by X-ray single-crystal structure analysis in some cases. Fifteen “end-off” compartmental ligands have been designed as 1 + 2 Schiff-base condensation product of 2,6-diformyl-4-R-phenol (R = Me, t Bu, and Cl) and five different amines, N-(2-aminoethyl)piperazine, N-(2-aminoethyl)pyrrolidine, N-(2-aminoethyl)morpholine, N-(3-aminopropyl)morpholine, and N-(2-aminoethyl)piperidine. Interestingly, in case of the combination of 2,6diformyl-4-methylphenol and N-(2-aminoethyl)morpholine/N-(3-aminopropyl)morpholine/N-(2-aminoethyl)piperidine 1 + 1 condensation becomes the reality and the ligands are denoted as L21−3. On reaction of copper(II) nitrate with L21−3 in situ complexes 3, 12, and 13 are formed having general formula Cu2(L21−3)2(NO3)2. The remaining 12 ligands obtained as 1 + 2 condensation products are denoted as L11−12, which produce complexes having general formula Cu2(L11−12)(NO3)2. Catecholase activity of all 15 complexes has been investigated in DMSO medium using 3,5-DTBC as model substrate. Treatment on the basis of Michaelis−Menten model has been applied for kinetic study, and thereby turnover number, kcat, values have been evaluated. Cyclic voltametric (CV) and differential pulse voltametric (DPV) studies of the complexes in the presence as well as in the absence of 3,5-DTBC have been thoroughly investigated in DMSO medium. From those studies it is evident that oxidation of 3,5-DTBC catalyzed by dicopper(II) complexes proceed via two steps: first, semibenzoquinone followed by benzoquinone with concomitant reduction of CuII to CuI. Our study reveals that apparently there is nearly no linear relationship between kcat and E° values of the complexes. However, a detailed density functional theory (DFT) calculation sheds light on this subject. A very good correlation prevails in terms of the energetics associated with the CuII to CuI reduction process and kcat values, as revealed from the combined theoretical and experimental approach.
■
INTRODUCTION
conversion is of great importance in medical diagnosis for determination of the hormonally active catecholamines adrenaline, noradrenaline, and dopa.4 Studies on the synthetic analogues of catechol oxidase have become an emerging field of modern research in bioinorganic chemistry to understand the in-depth functional mechanisms of catechol oxidase and elucidate different internal and external
Catechol oxidase (CO), the less well-known member of the type-3 copper proteins, found in plant tissues and crustaceans, catalyzes the oxidation of a wide range of o-diphenols (catechols) to the corresponding o-quinones, a process known as catecholase activity (Scheme 1).1−3 The o-quinones are highly reactive compounds and readily undergo an autopolymerization leading to formation of melanin, a brown polyphenolic pigment, thought to protect damaged tissues against pathogens or insects.1 The catechol to quinone © 2014 American Chemical Society
Received: March 5, 2014 Published: July 29, 2014 8257
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
Scheme 1. Reaction Pathway of the Oxygenation and Oxidation Catalyzed by TYR and CO
Scheme 2. Structure and Labels of Synthesized Compartmental Ligands
factors which may influence the activity.1−3,5−18 It is now well documented that among the external factors, the most important contributors are the pH of the medium, the nature of the model substrate, and the nature of the solvent. As far as internal factors are concerned the dominating effects come from the Cu−Cu distance, the flexibility of the coordinating ligands, and the coordination environment of the copper centers. Recently, it was documented that the presence of a positive charge center close to the metal center(s) may enhance the activity. Oxidation of catechol to quinone is a two-electron oxidation process. Most of the researchers working in elucidating the mechanistic pathway involved in catecholase activity believe that metal-centered redox participation, i.e., reduction of CuII to CuI, is the key step responsible for oxidation of catechol, although a radical pathway is now emerging to be a reasonable alternative, at least for nickel and zinc systems.19 However, since redox participation of two copper centers appears to be crucial in catecholase activity, the electrochemical behavior of the synthetic analogues of CO could clarify the mechanism of their catecholase activity. Surprisingly, only a few scattered reports are in the literature5,9,20−24 where this particular issue has been pointed out. Krebs et al., in one of their classical papers,1 after thorough electrochemical analysis came to the conclusion that the poorly defined redox chemistry of their Mannich base dicopper complexes could not allow them to establish any relationship between the electrochemical properties of the complexes and the oxidation potential of the model substrate, 3,5-di-tertbutylcatechol (3,5-DTBC). However, Neves et al.10 observed a very good correlation between (ΔE)1,2 (= E(red)1 − E(red)2) and kinetic parameter k2 (second-order rate constant = kcat/ KM) in their dinuclear Cu(II) complexes, although they failed to find any correlation between kcat/KM and the individual E(red)1 and E(red)2 of the complexes. In order to get a better understanding on this intriguing issue we planned to investigate the electrochemical property and catecholase activity of as many as 15 closely related dicopper(II) complexes in DMSO medium extensively. To fulfill our purpose we designed 15 “end-off” compartmental ligands as 1 + 2 the Schiff-base condensation product of 2,6diformyl-4-R-phenol (R= Me, tBu and Cl) and five different amines, N-(2-aminoethyl)piperazine, N-(2-aminoethyl)pyrrolidine, N-(2-aminoethyl)morpholine, N-(3-aminopropyl)morpholine, and N-(2-aminoethyl)piperidine (Scheme 2). In three cases, instead of 1 + 2 condensation, 1 + 1 condensation becomes a reality. In those cases the diformyl species is 2,6diformyl-4-methylphenol and amines are N-(2-aminoethyl)piperidine. N-(2-Aminoethyl)morpholine and N-(3aminopropyl)morpholine and the ligands are denoted as L21−3. On reaction of copper(II) nitrate with L21−3 in situ complexes 3, 12, and 13 are generated having general formula Cu2(L21−3)2(NO3)2. The remaining 12 ligands obtained as 1 + 2 condensation products are denoted as L11−12 which produce complexes 1, 2, 4−11, 14, and 15 having general formula Cu2(L11−12)(NO3)2. Catecholase activity of all 15 complexes
has been investigated in DMSO medium using 3,5-DTBC as model substrate. Cyclic voltametric (CV) and differential pulse voltametric (DPV) studies of the complexes in the presence as well as in the absence of 3,5-DTBC have been thoroughly investigated in DMSO medium. A detailed DFT calculation has been performed to make a breakthrough on this subject. Although no linear relationship has apparently been established between E° and kcat values, a very good correlation prevails in terms of the energetics associated with the electrochemical behavior of the catechol-bound dicopper(II) system involved in the catalytic process and kcat values, as revealed from the combined theoretical and experimental approach. All those interesting findings are vividly portrayed in this article.
■
EXPERIMENTAL SECTION
Physical Methods and Materials. Elemental analyses (carbon, hydrogen, and nitrogen) were performed using a Perkin−Elmer 240C elemental analyzer, and copper content was estimated gravimetrically. Infrared spectra were recorded on KBr disks (400−4000 cm−1) with a Perkin−Elmer RXI FTIR spectrophotometer. Electronic spectra (200−800 nm) were measured at room temperature on a Shimadzu UV-3101PC using dry acetonitrile/dry methanol/nujol as medium. Cyclic voltammetric and DPV measurements were performed using a CH1106A potentiostat with glassy carbon (GC) as the working electrode, Pt wire as the counter electrode, and Ag,AgCl/sat KCl as the reference electrode. All solutions were purged with dinitrogen prior to measurements. High-purity N-(2-aminoethyl)piperazine, N-(2-aminoethyl)pyrrolidine, N-(2-aminoethyl)morpholine, N-(3-aminopropyl)morpholine, and N-(2-aminoethyl)piperidine were purchased from commercial sources (Fluka, Lancaster Chemical Co. Inc., Aldrich) and used as received. The 2,6-diformyl-4-R-phenols (R = methyl, tert-butyl, chloro) were prepared according to the literature method.25 Solvents were dried according to standard procedures and distilled prior to use. All other chemicals used were of AR grade. Synthesis of the Complexes. The following general template synthetic route was adopted for preparing all complexes. A methanolic solution of copper nitrate was added to the ligand solution formed in situ via condensation of 2,6-diformyl-4-R-phenol (where R= methyl, tert-butyl, and chloro) with the corresponding amines maintaining the same molar ratio. The preparation, composition, and other physicochemical characteristics of all complexes using the template technique are given below. Syntheses. a. [Cu2(L11)(OH)(H2O)(NO3)2], 1. To a methanolic solution (5 mL) of N-(2-aminoethyl)piperidine (0.256 g, 2 mmol) a methanolic solution (10 mL) of 4-chloro-2,6-diformylphenol (0.180 g, 1 mmol) was added in a dropwise manner, and the mixture was refluxed 30 min. Then a methanolic solution (15 mL) of copper nitrate (0.604 g, 2.5 mmol) was added to it, and the resulting mixture was allowed to reflux for 2 h. The deep green solution was filtered and 8258
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
Table 1. Crystallographic Data and Details of Refinements for Complexes 1−3 and 8−10a empirical formula M cryst syst space group a/Å b/Å c/Å α/deg β/deg γ/deg V/Å3 Z Dcalcd/g cm−3 μ/mm−1 F(000) total data unique data Rint reflns I > 2σ(I) params goodness-of-fit R1 wR2 (I > 2σ(I))b residuals/e Å−3
1·H2O
2·H2O
3·2H2O
8
9·H2O
10·2H2O
C22H35ClCu2N6O9 690.09 triclinic P1̅ 9.531(6) 12.534(6) 13.472(7) 100.526(5) 103.720(5) 106.472(4) 1444.3(13) 2 1.587 1.623 712 7223 4207 0.0510 2450 388 1.094 0.0757 0.1462 0.767, −0.710
C26H44Cu2N6O9 711.75 monoclinic P21/n 15.907(3) 11.836(2) 17.008(3)
C32H50Cu2N8O20 993.88 triclinic P1̅ 10.059(3) 10.470(3) 11.463(4) 77.391(3) 89.205(3) 66.130(3) 1073.6(6) 1 1.537 1.077 516 7278 3993 0.0284 2599 300 1.039 0.0756 0.2114 1.367, −0.525
C20H33ClCu2N6O10 680.05 monoclinic P21/m 8.1106(3) 12.3324(5) 13.7055(5)
C48H40Cu4N12O18 1367.40 tetragonal P42/n 24.820(3) 10.501(2)
C20H39ClCu2N10O17 854.14 orthorhombic Pna21 8.7183(5) 40.122(3) 9.5437(6)
6469.0(18) 4 1.404 1.369 2848 24566 4543 0.0667 2455 427 1.167 0.0703 0.2241 0.744, −0.355
3338.3(4) 4 1.699 1.443 1760 33586 5005 0.1276 3588 472 0.951 0.0488 0.0967 0.492, −0.271
94.178(2) 3193.5(10) 4 1.480 1.390 1488 21837 5720 0.0348 4446 403 1.051 0.0519 0.1594 1.249, −0.401
102.507(2) 1338.34(9) 2 1.688 1.753 684 24430 7253 0.0358 5172 218 1.015 0.0445 0.1126 0.685, −0. 964
a
CCDC 988842−988847 contain the supplementary crystallographic data for this paper. These data can be obtained free charge from the Cambridge Crystallographic Data Centre via http://www.ccdc.cam.ac.uk/5datarequest/cif. bR1 = Σ|Fo| − |Fc||/Σ|Fo|, wR2 = [Σw (Fo2− Fc2)2/Σw (Fo2)2]1/2. Found: C, 40.24; H, 5.66; N, 11.69. IR (KBr): ν(CN) 1636 cm−1; ν(skeletal vibration) 1553 cm−1; ν(H2O) 3435 cm−1; ν(NO3−) 1384 cm−1. f. [Cu2(L15)(OH)(H2O)(NO3)2], 6. 6 was prepared adopting the same procedure as for 4 using N-(2-aminopropyl) morpholine (0.144g, 1 mmol) instead of N-(2-aminoethyl) morpholine. A deep green colored crystalline solid mass was obtained from the filtrate after a few days (yield 80%). Anal. Calcd for C22H35 ClCu2N6O11: C, 36.64; H, 4.85; N, 11.65. Found: C, 36.69; H, 4.69; N, 11.69. IR (KBr): ν(CN) 1643 cm−1; ν(skeletal vibration) 1550 cm−1; ν(H2O) 3413 cm−1; ν(NO3−) 1384 cm−1. g. [Cu2(L16)(OH)(H2O)(NO3)2], 7. 7 was prepared following the same procedure as for 5 using N-(2-aminopropyl)morpholine (0.144g, 1 mmol) instead of N-(2-aminoethyl)morpholine. Deep green colored single crystals suitable for X-ray data collection were obtained from the filtrate after a week on keeping the filtrate in a CaCl2 desiccator (yield 80%). Anal. Calcd for C26H44Cu2N6O11: C, 42.04; H, 4.71; N, 11.32. Found: C, 42.14; H, 4.65; N, 11.39. IR (KBr): ν(CN) 1642 cm−1; ν(skeletal vibration) 1551 cm−1; ν(H2O) 3434 cm−1; ν(NO3−) 1384 cm−1. h. [Cu2(L17)(OH)(H2O)2(NO3)2], 8. 8 was prepared following the same procedure as for 1 using N-(2-aminoethyl)pyrrolidine (0.114g, 1 mmol) instead of N-(2-aminoethyl)piperidine. Crystals suitable for Xray diffraction were obtained from the filtrate after 2 days on keeping the filtrate at room temperature (yield 80%). Anal. Calcd for C20H33ClCu2N6O10: C, 35.21; H, 4.86; N, 12.4. Found: C, 35.29; H, 4.76; N, (12.35. IR (KBr): ν(CN) 1650 cm−1; ν(skeletal vibration) 1551 cm−1; ν(H2O) 3435 cm−1; ν(NO3−) 1384 cm−1. I. [Cu2(L19)(OH)(H2O)(NO3)2], 9. 9 was prepared following the same procedure as for 2 using N-(2-aminoethyl)pyrrolidine (0.114g, 1 mmol) instead of N-(2-aminoethyl)piperidine. Crystals suitable for Xray diffraction were obtained from the filtrate after several days on keeping the filtrate in a CaCl2 desiccator. A deep green colored solid mass was obtained from the filtrate after a few days on keeping the filtrate at room temperature (yield 80%). Anal. Calcd for C24H39Cu2N6O9: C, 42.18; H, 5.71; N, 12.30. Found: C, 42.14; H,
kept in a CaCl2 desiccator. Deep green crystals suitable for X-ray data collection were obtained from the filtrate after a few days (yield 75%). Anal. Calcd for C23H33ClCu2N5O9: C, 31.90; H, 3.81; N, 8.08. Found: C, 31.87; H, 3.86; N, 7.99. IR (KBr): ν(CN) 1651 cm−1; ν(skeletal vibration) 1552 cm−1; ν(H2O) 3413 cm−1; ν(NO3−) 1384 cm−1. b. [Cu2(L12)(OH)(H2O)(NO3)2], 2. 2 was prepared by adopting the same procedure as in the case of 1 using 4-tert-butyl-2,6diformylphenol (0.206 g, 1 mmol) instead of 4-chloro-2,6diformylphenol. Deep green colored single crystals suitable for X-ray data collection were obtained from the filtrate after a week (yield 65%). Anal. Calcd for C26H44Cu2N6O9: C, 43.83; H, 6.18; N, 11.98. Found: C, 43.75; H, 6.10; N, 11. 9. IR (KBr): ν(CN) 1639 cm−1; ν(skeletal vibration) 1550 cm−1; ν(H2O) 3435 cm−1; ν(NO3−) 1384 cm−1. c. [Cu(L21)(H2O)(NO3)]2, 3. 3 was prepared following the same procedure as for 1 using 4-methyl-2,6-diformylphenol (0.184 g, 1 mmol) instead of 4-chloro-2,6-diformylphenol. Single crystals suitable for diffraction were obtained from the filtrate after 2 days (yield 80%). Anal. Calcd for C32H50Cu2N8O20: C, 38.62; H, 5.03; N, 11.40. Found: C, 38.63; H, 5.09; N, 11.26. IR (KBr): ν(CO) 1680 cm−1; ν(CN) 1653 cm−1; ν(skeletal vibration) 1552 cm−1; ν(H2O) 3429 cm−1; ν(NO3−) 1383 cm−1. d. [Cu2(L13)(OH)(H2O)(NO3)2], 4. 4 was prepared following the same procedure as for 1 just using N-(2-aminoethyl) morpholine (0.130 g,1 mmol) in place of N-(2-aminoethyl)piperidine. A deep green colored solid mass was obtained from the filtrate after a few days on keeping the filtrate at room temperature (yield 80%.) Anal. Calcd for C20H31ClCu2N6O11: C, 34.65; H, 4.47; N, 12.13. Found: C, 34.54; H, 4.58; N, 12.19. IR (KBr): ν(CN) 1650 cm−1; ν(skeletal vibration) 1552 cm−1; ν(H2O) 3436 cm−1; ν(NO3−) 1384 cm−1. e. [Cu2(L14)(OH)(H2O)(NO3)2], 5. 5 was prepared adopting the same procedure as for 2 by replacing N-(2-aminoethyl) piperidine with N(2-aminoethyl) morpholine (0.130 g, 1 mmol). Deep green colored single crystals suitable for X-ray data collection were obtained from the filtrate after a week on keeping the filtrate in a CaCl2 desiccator (yield 80%). Anal. Calcd for C24H40Cu2N6O11: C, 40.33; H, 5.60; N, 11.76. 8259
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
5.68; N, 12.39. IR (KBr): ν(CN) 1644 cm−1; ν(skeletal vibration) 1552 cm−1; ν(H2O) 3422 cm−1; ν(NO3−) 1371 cm−1. j. [Cu2(H2L110)(OH)(H2O)(NO3)4], 10. 10 was prepared following the same procedure as 4 using N-(2-aminoethyl) piperazine (0.129g, 1 mmol) instead of N-(2-aminoethyl)morpholine. Single crystals suitable for diffraction were obtained after a week from the filtrate kept in a CaCl 2 desiccator (yield 80%). Anal. Calcd for C20H39ClCu2N10O17: C, 28.12; H, 4.6; N, 16.4. Found: C, 28.09; H, 4.56; N, 16.39. IR (KBr): ν(CN) 1647 cm−1; ν(skeletal vibration) 1551 cm−1; ν(H2O) 3419 cm−1; ν(NO3−) 1368 cm−1. k. [Cu2(L112)(OH)(H2O)(NO3)2], 11. 11 was prepared following the same procedure like 8 using 4-tert-butyl-2,6-diformylphenol (0.206 g, 1 mmol) instead of 4-chloro-2,6-diformylphenol. Deep green colored crystals were obtained from the filtrate after a few days on keeping the filtrate at room temperature (yield 80%). Anal. Calcd for C24H48Cu2N10O17: C, 32.95; H, 5.49; N, 16.01. Found: C, 32.88; H, 5.6; N, 15.89. IR (KBr): ν(CN) 1645 cm−1; ν(skeletal vibration) 1549 cm−1; ν(H2O) 3413 cm−1; ν(NO3−) 1384 cm−1. Complexes 12−15 of composition [Cu(HL2 2 )(H 2 O)(NO3)]2(NO3)2·2H2O, [Cu(L23)(H2O)(NO3)]2, [Cu2(L18)(OH)(H2O)2](NO 3)2, and [Cu2(H2L111)(OH)(H2O)(NO3)](NO3)3· 2H2O were prepared by following the method previously reported by us.18 X-ray Data Collection and Crystal Structure Determinations. Data collection of the structures reported was carried out on a Bruker Smart CCD diffractometer equipped with graphite-monochromated Mo Kα radiation (λ = 0.71073 Å) at room temperature. Cell refinement, indexing, and scaling of the data set were carried out using Bruker Smart Apex and Bruker Saint packages.26 Structures were solved by direct methods and subsequent Fourier analyses27 and refined by the full-matrix least-squares method based on F2 with all observed reflections.27 In 1 a nitrate anion is disordered over two positions at half occupancy with an oxygen atom (O8 and O11) sharing the same site with a water molecule (O1w and O2w, respectively, with occupancy = 0.5). In 2 two residuals in the difference Fourier map were treated as lattice water molecules at half occupancy (H atoms not located). In 3 the piperidine ring and a nitrate anion were found disordered over two positions (occupancies of 0.57/0.43 and 0.76/0.24, respectively); thus, the piperidine carbon atoms were isotropically refined. In 8 the uncoordinated nitrate anion was found disordered over two positions: one having all atoms and the second only the nitrogen and one oxygen sitting in a crystallographic symmetry plane. The crystal of 9 was low diffracting due to the disorder over two positions observed in the tert-butyl group and in the methylene groups attached at N3 and N4 (occupancies of 0.64/0.36, 0.63/37, and 0.64/0.36). Moreover, beside nitrate N3, located on a center of symmetry, to nitrate N7 an occupancy of 0.50 was assigned for the sake of neutrality. Non-hydrogen atoms were anisotropically refined except the tert-butyl in complex 9. All calculations were performed using WinGX System, Ver 1.80.05.28 Pertinent crystallographic data and refinement details are summarized in Table 1. Theoretical Methods. Energies of all complexes included in this study were computed at the BP86-D3/def2-TZVPD level of theory using the optimized geometries within the program TURBOMOLE version 6.4.29 Solvent effects were evaluated using the COSMO (Conductor like Screening Model) solvation model.30 For the calculations we used the BP86 functional with the latest available correction for dispersion (D3).
ray analysis were obtained. FT-IR spectral study reveals that all complexes exhibit bands due to the CN stretch in the range 1630−1648 cm−1 and skeletal vibration in the range 1545− 1555 cm−1. The presence of weakly coordinated NO3− anion in all complexes is indicated by the broad band centered in the range 1365−1385 cm−1.31 Electronic spectra of all complexes studied in DMSO medium display very similar absorption bands in the range 370−415 and 640−650 nm. The observed lower energy band may be assigned due to the d−d transition, and the corresponding strong higher energy single band (between 370 and 415 nm) is due to combination of both phenoxido−Cu(II) and hydroxido−Cu(II) LMCT bands. It is well known that for a d9 system the electronic transition 2Eg → 2 T2g is expected to take place at around 800 nm for octahedral coordination geometry, and this band undergoes a significant blue shift for octahedral distortions to square-pyramidal and square-planar structures.32 In all complexes the d−d transition positions are in agreement with a square-pyramidal geometry around the metal centers. Description of Crystal Structures. Structural analysis of complexes 1, 2, 8, 9, and 10 reveals for all a [Cu2(L1)(OH)] unit with metals chelated by the pentacoordinating ligand L1 and bridged by a hydroxido oxygen. In each case the copper atoms complete their square planar coordination sphere through nitrate anions or aqua ligands. ORTEP drawings of complexes 1, 8, and 9 are shown in Figures 1−3; those of 2 and 10 are given in Figures S12 and S13, Supporting Information.
Figure 1. ORTEP drawing of the complex cation of complex 1.
Complexes 1 and 2 show a close comparable metal coordination environment and differ only for the R group at the Schiff base (Cl and tBu, respectively). Copper(II) atoms are chelated by the phenoxido oxygen, the imino and amino piperidine nitrogens of the ligand, and by a bridging hydroxide, completing the coordination sphere by a nitrate anion that acts as a bridging species toward the metals. In both complexes the basal Cu−O distances are comparable in lengths ranging from 1.897(3) to 2.006(7) Å (Table 2), showing a shorter value for the hydroxido oxygen. On the other hand, the Cu−N(imino) bond distances (range 1.919(4)−1.959(9) Å) are shorter than those involving the amino nitrogen of the heterocycle ring, falling between 2.036(4) and 2.097(8) Å) Å. The Cu− O(nitrate) distance at the apical position of the square planar
■
RESULTS AND DISCUSSION Four structurally characterized copper(II) complexes (12−15) were prepared by adopting the template synthesis technique and identified by their reported physicochemical properties.18 The 11 new complexes, namely, complexes 1−11, were synthesized by treating a methanolic solution of copper(II) nitrate trihydrate with the Schiff base formed in situ between 2,6-diformyl-4-R-phenol (where R= tert-butyl, chloro, methyl) and the diamines. In some cases single crystals suitable for X8260
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
Table 2. Coordination Bond Lengths (Angstroms) and Angles (degrees) for Complexes 1 and 2 Cu(1)−O(1) Cu(1)−O(2) Cu(1)−N(1) Cu(1)−N(3) Cu(1)−O(3) Cu(2)−O(1) Cu(2)−O(2) Cu(2)−N(2) Cu(2)−N(4) Cu(2)−O(4) Cu(1)−Cu(2) O(2)−Cu(1)−N(1) O(2)−Cu(1)−O(1) N(1)−Cu(1)−O(1) O(2)−Cu(1)−N(3) N(1)−Cu(1)−N(3) O(1)−Cu(1)−N(3) O(2)−Cu(1)−O(3) N(1)−Cu(1)−O(3) O(1)−Cu(1)−O(3) N(3)−Cu(1)−O(3) O(2)−Cu(2)−N(2) O(2)−Cu(2)−O(1) N(2)−Cu(2)−O(1) O(2)−Cu(2)−N(4) N(2)−Cu(2)−N(4) O(1)−Cu(2)−N(4) O(2)−Cu(2)−O(4) N(2)−Cu(2)−O(4) O(1)−Cu(2)−O(4) N(4)−Cu(2)−O(4)
Figure 2. ORTEP drawing of the complex cation of complex 8 located on the crystallographic symmetry plane.
1
2
2.006(7) 1.925(6) 1.951(8) 2.097(8) 2.360(8) 1.980(6) 1.937(7) 1.959(9) 2.075(7) 2.437(9) 2.944(2) 167.6(3) 80.7(3) 90.0(3) 103.9(3) 85.4(4) 175.4(3) 97.2(3) 90.4(3) 87.4(3) 92.9(3) 170.6(3) 81.0(3) 90.6(3) 102.2(3) 85.9(3) 175.0(3) 93.3(3) 89.8(3) 83.1(3) 100.4(3)
1.961(3) 1.897(3) 1.928(4) 2.037(4) 2.434(3) 1.970(3) 1.896(3) 1.919(4) 2.036(4) 2.413(3) 2.9049(8) 165.85(16) 79.69(12) 90.08(14) 103.25(13) 86.17(15) 174.41(13) 94.20(14) 95.40(16) 89.12(12) 95.37(13) 165.44(16) 79.52(12) 89.86(14) 103.40(13) 86.54(15) 175.03(13) 95.54(13) 94.28(15) 89.46(12) 94.22(13)
behave as H-bond donors toward nitrate oxygens and a lattice water molecule (N···O distance of 2.99−3.03 Å). The intermetallic distance in complex 9 is 2.934(2) Å, while in 10 it is 2.9596(11) Å, which is the longest among those here reported. Complex 8 is located on a crystallographic symmetry plane passing in between the metals and referring the two halves of the phenolato ligand. The crystal structure shows that the equatorial plane of copper atoms is realized by doubly bridging phenoxido and hydroxido oxygen atoms and by the imine- and amine-nitrogen donors from L17 and from the pyrrolidine ring. The coordination sphere of the crystallographically independent copper ion is completed at axial sites by a water molecule (Cu−O1w = 2.5495(16) Å) and by a nitrate oxygen (Cu−O3′) at 2.7669(18) Å, indicating a weaker interaction and providing a pseudo-octahedral geometry as a consequence of the Jahn− Teller effect. The basal Cu−O distances are 1.9063(10) and 1.9726(9) Å (Table 4), the shorter value being relative to the hydroxido oxygen. On the other hand, the Cu−N(imino) bond distance (1.9393(13) Å) is slightly shorter than that involving the amino nitrogen of the heterocycle ring 2.0278(13) Å. The piperazine moieties attain a chair conformation. The bridging oxygen atoms locate the metals separated by 2.9450(4) Å. Complexes 12 and 13 built with the methyl-phenolato moiety, already reported by us,18 will not be described as they are isomorphous and isostructural with chloro derivative complexes 10 and 8, respectively. X-ray structural data of all
Figure 3. ORTEP drawing of the complex cation of complex 9.
coordination geometry is sensibly longer falling in a range from 2.360(8) to 2.437(9) Å. The Cu(1)−O(1)−Cu(2) bond angle is ca. 95.3°, whereas the Cu(1)−O(2)−Cu(2) one, subtended by the hydroxido group, is slightly larger with an average value of 99.7°, leading to a metal−metal separation of 2.944(2) and 2.9049(8) Å in 1 and 2, respectively. The metal coordination environment of complexes 9 and 10 comprises the phenoxido oxygen and the imine- and aminenitrogen donors, the bridging hydroxo group, and the square pyramidal geometry is completed by a monocoordinated nitrate and an aqua ligand on opposite side of the Schiff base ligand plane. Bond distances follow the trend observed for complexes 1 and 2 and typical for these types of dicompartmental copper species. A slight difference is detected between the Cu− N(pyrrolidine) bond lengths in complex 9 of 2.043(12) and 2.007(12) Å with respect to the Cu−N(piperazine) ones of 2.113(6) and 2.083(6) Å in 10, Table 3. In complex 9 the nitrate N5 is located on a crystallographic inversion center, and it appears disordered, connecting two symmetry related complexes (Figure S14, Supporting Information). In complex 10 the piperazine moieties, which attain a chair conformation, are protonated at the amine nitrogens leading to a tricationic complex species. This feature is confirmed by the number of nitrate counteranions and the NH2 groups which 8261
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
independent part is shown in Figure 4, and a selection of bond lengths and angles is given in Table 5. The metal
Table 3. Coordination Bond Lengths (Angstroms) and Angles (degrees) for Complexes 9 and 10 Cu(1)−O(1) Cu(1)−O(2) Cu(1)−N(1) Cu(1)−N(3) Cu(1)−O(1w) Cu(2)−O(1) Cu(2)−O(2) Cu(2)−N(2) Cu(2)−N(4) Cu(2)−O(4) Cu(1)−Cu(2) O(2)−Cu(1)−N(1) O(2)−Cu(1)−O(1) N(1)−Cu(1)−O(1) O(2)−Cu(1)−N(3) N(1)−Cu(1)−N(3) O(1)−Cu(1)−N(3) O(1)−Cu(1)−O(1w) O(2)−Cu(1)−O(1w) N(1)−Cu(1)−O(1w) N(3)−Cu(1)−O(1w) N(2)−Cu(2)−O(2) N(2)−Cu(2)−O(1) O(2)−Cu(2)−O(1) N(2)−Cu(2)−N(4) O(2)−Cu(2)−N(4) O(1)−Cu(2)−N(4) N(2)−Cu(2)−O(4) O(2)−Cu(2)−O(4) O(1)−Cu(2)−O(4) N(4)−Cu(2)−O(4)
9
10
1.977(9) 1.933(10) 1.932(11) 2.043(12) 2.393(14) 1.973(9) 1.918(10) 1.934(11) 2.007(12) 2.38(3) 2.934(2) 170.8(5) 80.6(4) 90.7(5) 102.8(5) 85.2(5) 168.9(5) 96.6(4) 92.4(4) 91.6(5) 93.9(5) 168.0(5) 89.7(4) 81.1(4) 85.7(5) 103.2(5) 174.8(4) 108.7(6) 77.2(6) 79.5(8) 104.3(8)
1.979(4) 1.929(5) 1.931(6) 2.113(6) 2.305(6) 1.972(5) 1.919(5) 1.919(6) 2.083(6) 2.398(6) 2.9596(11) 170.0(2) 80.4(2) 90.7(2) 103.3(2) 84.7(3) 169.5(2) 90.9(2) 89.7(2) 95.0(2) 98.9(2) 169.1(2) 90.3(2) 80.80(18) 85.9(2) 102.4(2) 173.2(2) 92.5(2) 93.5(2) 88.4(2) 97.4(2)
Figure 4. ORTEP drawing of the complex cation of complex 3 located on a crystallographic inversion center.
Table 5. Coordination Bond Lengths (Angstroms) and Angles (degrees) for Complex 3a Cu−O(1) Cu−O(1′) Cu−O(2′) Cu−N(1) O(1)−Cu−O(2′) O(1)−Cu−O(1′) O(2′)−Cu−O(1′) O(1)−Cu−N(1) O(2′)−Cu−N(1) O(1′)−Cu−N(1) O(1)−Cu−O(1w) O(2′)−Cu−O(1w) a
1.940(4) 1.951(4) 1.944(4) 1.957(5) 170.23(19) 78.48(18) 92.30(19) 93.28(18) 95.5(2) 169.84(19) 94.7(3) 89.1(3)
Cu−O(1w) Cu−O(21) Cu−Cu#1 O(1′)−Cu−O(1w) N(1)−Cu−O(1w) O(1)−Cu−O(21) O(1w)−Cu−O(21) O(21)−Cu−N(1) O(21)−Cu−O(1′) O(21)−Cu−O(2′)
2.231(6) 2.582(9) 3.0137(14) 93.9(3) 92.6(3) 84.7(2) 171.5(3) 78.9(3) 94.3(3) 92.9(3)
Primed atoms at #1 −x + 1, −y, −z + 1.
Table 4. Coordination Bond Lengths (Angstroms) and Angles (degrees) for Complex 8a Cu−O(1) Cu−O(2) Cu−N(1) O(2)−Cu−N(1) O(2)−Cu−O(1) N(1)−Cu−O(1) O(2)−Cu−N(3) N(1)−Cu−N(3) a
1.9726(9) 1.9063(10) 1.9393(13) 170.89(5) 80.81(5) 90.84(5) 102.27(6) 85.96(6)
Cu−N(3) Cu−O(1w) Cu−O(3)′ Cu−Cu′ O(1w)−Cu−O(1) O(1w)−Cu−O(2) O(1w)−Cu−N(1) O(1w)−Cu−N(3) O(1w)−Cu−O(3)′ O(1)−Cu−N(3)
2.0278(13) 2.5495(16) 2.7669(18) 2.9450(4) 89.80(6) 95.20(6) 88.37(6) 91.77(5) 166.46(6) 176.40(5)
possesses an octahedral coordination sphere with two phenoxido-bridged oxygens, an imine nitrogen donor, and a carbonyl oxygen located in the equatorial plane and an aqua and a nitrate oxygen occupying the axial positions at longer distances. The Cu−N and Cu−O bond distances are comparable to those reported for the complexes containing the symmetric ligand L1 and one ligand ranging from 1.940(4) to 1.957(4) Å. The axial distances are significantly elongated for the Jahn−Teller effect, the Cu−O(1w) bond length, 2.231(6) Å, is significantly shorter than that involving the nitrate oxygen (Cu−O(4) 2.582(9) Å). The bond angle Cu−O1−Cu′ (101.52(18)°) leads to an intermetallic separation of 3.0137(14) Å, which is ca. 0.1 Å longer with respect to the values measured in the other complexes here reported. The piperidine moieties, in the usual chair conformation, are far apart from the metal centers with the alkyl chain in an anti periplanar conformation (torsion angle N1−C−C−N2 = 168.5(6)°). From the above structural description of the complexes it is evident that in a few cases the X-ray data suffers from disorder. Our all efforts to get high-quality crystals on those cases failed, and we were unsuccessful to resolve the issue. However, despite the relatively poor figures of merit in some cases, which makes the absolute values of the metrical parameters for them somewhat inaccurate, the relative positions of the atoms within
Primed atoms at x, −y + 1/2, z.
these complexes do not indicate any significant effect of the group in the para position of the phenolato ligand (Cl, Me, and t-Bu) on the coordination distances and on the intermetallic separation, and the slight differences observed may be attributed to a packing effect and to a lesser extent to H bonds realized in the crystal structure. Different from previous complexes, the crystal structure of complex 3 reveals that it comprises a dinuclear complex cation located on a crystallographic inversion center, beside nitrate anions, and lattice water molecules. The complex is formed by two unsymmetrical tridentate ligands (L2) chelating the metal ions and forming a phenoxido-bridged Cu(II) dimer. An ORTEP view of 3 with atom-labeling scheme of the 8262
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
Table 6. Electrochemical Data of Complexes 1−15 in DMSO Medium DPV data Epc/V
a
complex
CV dataa E0/V or Epc/V
Cu(II)/Cu(I)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
0.08 (178), −0.56 (Epc) 0.04(175), −0.51(212) 0.33(210), −0.56 (Epc), −0.81 (Epc) 0.06(230), −0.41(140) 0.06(98) 0(320), −0.45 (Epc), −1.26 (Epc) 0.14(100)b, −0.31 (Epc), −0.63(240) 0.07(165), −0.47 (230) 0.11(106), −0.67 (Epc), −0.95 (Epc) 0.09(183), −0.49 0.05(140), −0.49 (Epc) 0.04 (230), −0.42(180) −0.06(320), −0.55 (Epc) −0.75 (Epc) 0.06(230), −0.45
0.06 0.02 0.30, 0.15 0.08 0.22,b 0.08 0.28,b 0.10, −0.06 0.19,b −0.16 0.08, −0.06 0.13 0.10 0.07 0.04 0.23b, 0.06, −0.11 0.07
Cu(I)/Cu(0) −0.42 −0.44, −0.44, −0.37 −0.40, −0.36 −0.63 −0.46 −0.48, −0.33 −0.33, −0.41 −0.45 −0.50 −0.40
−0.48 −0.70 −0.46
−0.83 −0.42
Error range: ± 0.005 V. bPeaks are probably due to some unidentified degradation product.
Figure 5. CV and DPV spectrum of complex 10 at the GC electrode at 100 mV s−1 scan rate as representative.
value (ΔEp = Epa − Epc), or in few cases as an irreversible reduction redox process with ipc/ipa ≪ 1, which is typical of Cu(II)/Cu(I) reduction in many of its complexes.10,11,20 It has been argued that though the formal potential (E°′) is equal to the mid peak potential (Epa + Epc)/2 determined by cyclic voltammetry, and this relationship is strictly valid for ideal reversible couples; for quasi-reversible systems the mid peak potential may also be used as the formal potential. This is because for quasi-reversible systems the amount of shift of the cathodic peak toward more negative potential values is compensated by an almost equal shift of the reoxidation peak toward more positive values.33 Thus, in our case also we assumed the mid peak potential in CV experiments is very close to the formal potential (E°′). In fact, from Table 6 it can be seen that the DPV peak potentials are very close to the E°′ value obtained from CV experiments. As in our case DPV peak potentials can be more accurately measured than that of CV peaks, we have taken the DPV peak potentials for our attempted correlation study with the kcat values. Where the two Cu(II)/Cu(I) potentials were observed in DPV we have taken an average value of the two potentials for correlation purpose. For each complex it is expected that the Cu(II)/Cu(I) reduction potentials should follow the order of Hammett σp values of the substituent on the diformyl phenyl ring; a more positive σp value (meaning more electron withdrawing character) should lead to more positive Cu(II)/Cu(I) potential. σp values for Cl, tBu, and Me groups are 0.227, −0.151, and
these complexes are sufficiently accurate to describe the shape of these molecules. Electrochemistry of the Complexes. The redox behavior of the complexes was explored using cyclic voltammetry (CV) and differential pulse voltammetry (DPV), and the data is tabulated in Table 6. The CV and DPV pictures of complex 10 are presented in Figure 5 as a representative one, and some other representative voltammograms are given in the Supporting Information (Figures S14−S17). The Cu(II)/ Cu(I) couple appears near 0.0 V with respect to the Ag/ AgCl electrode, and in most cases in the CV experiments the two Cu(II)/Cu(I) couples are not resolved, though in DPV these two couples are sometimes resolved. It has been argued that DPV is a very good technique for resolving redox responses having small differences in peak potentials, provided the two peaks differ in their formal potential by more than 180 mV.33 The fact that in our case we were unable to resolve the two Cu(II)/Cu(I) couples for most of the compounds by CV or DPV experiments suggests that the formal potentials for the two successive metal reductions differ by less than 180 mV, and hence, the mixed-valent CuICuII species will be unstable. Thus, during the catalytic cycle involving two-electron oxidation of catechol a one-electron redox intermediate is disfavored. The reductive process at around −0.4 V is assigned to Cu(I)/Cu(0) reduction as after cycling through the potential a strong stripping current is observed in the reverse (anodic) scan in CV experiments.11 In CV experiments the Cu(II)/Cu(I) couple(s) appears either as quasi-reversible, with a relatively high ΔEp 8263
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
Figure 6. Changes observed in UV−vis spectra of complexes (a) 2 and (b) 14 (conc. 1 × 10−4 M) upon addition of 100-fold 3,5-DTBC (1 × 10−2 M) in DMSO medium.
−0.170, respectively,34 and the Cu(II)/Cu(I) potential is expected to be most positive for the chloro derivatives. On the other hand, the σp values for tBu and Me are very comparable (a very small difference of ca. 0.02), so it is expected that these groups affect the potentials at a comparable degree and any difference between them will be largely determined by steric factors. This expected trend was found to be confirmed by complexes with six ligands L15, L16, L23, L110, L112, and H2L111, namely, complexes 6, 7, 13, 10, 11, and 15, respectively. In the case of complexes 1, 2, and 3 derived from three different ligands, namely, L11, L12, and L21, the Clsubstituted compound (complex 1) has an expectedly higher E° value than the tBu-substituted compound 2; however, the Mesubstituted compound 3 has much higher E° value than even the Cl-substituted compound 1, and this may be due to the different structure of the Me compound (complex 3). The electrochemical behavior of complexes 4, 5, and 12 is anomalous. Complex 5 having a tBu substituent gives the most positive E° value, whereas E° values of complex 4 (having a Cl substituent) and complex 12 (having a Me substituent) are very comparable. Among complex 8, 9, and 14 again the tBu substituent (complex 9) gives a higher E° value than the chloro species (complex 8); for the Me-substituent complex 12 of this series the Cu(II)/Cu(I) couple could not be detected. It is expected that complexes 6, 7, and 13 having a propyl group in the ligand backbone should show more facile Cu(II)/ Cu(I) reduction (more positive E° value) compared to the analogous complexes 4, 5, and 12, respectively, having an ethyl group in the ligand backbone, because of their greater steric flexibility, which allows structural rearrangement on electron transfer. This is indeed the case for the chloro (complexes 4 and 6) and methyl (complexes 13 and 12) series of compounds, but this trend is not observed for tBu series, where compound 5 has an abnormally high positive value of the Cu(II)/Cu(I) potential compared to the other two compounds of that series. Thus, the results discussed above indicate that the R group at the para position affects the electrochemical property of the complexes, although in some cases the relationship is not unambiguously interpreted. It is worth noting that when the binuclear Cu(II) complexes were reacted with 3,5-DTBC, two new peaks at 0.3−0.4 and 0.1−0.2 V were detected in DPV experiments. It has been reported that CuII(bpy)(DTBC) and CuII(phen)(DTBC) complexes undergo a two-electron oxidation leading to
DTBQ at around 0.05 V vs SCE and around 0.1 V vs Ag/ AgCl. In the same potential region CuII(en)(DTBC) undergoes two successive one-electron oxidations.35 On the basis of the above observations we assign the new peaks at 0.1−0.2 and 0.3−0.4 V to Cu(II)-bound catechol to semiquinone and semiquinone to quinone oxidations (Figure S19, Supporting Information). Catecholase Activity. Catechol oxidase, a type-3 copper protein, can bind oxygen reversibly at room temperature, and so it can be used to oxidize phenols to the respective obenzoquinones. As a model of the enzyme we have taken 15 dinuclear complexes of Cu(II) and studied their efficiency toward oxidation of 3,5-di-tert-butylcatechol (3,5-DTBC) to 3,5-di-tert-butylbenzoquinone (3,5-DTBQ). Interestingly, all of them display significant catalytic activity toward oxidation of 3,5-DTBC to 3,5-DTBQ in DMSO medium. Before proceeding into a detailed kinetic study to evaluate the ability of our complexes to mimic the active site of catechol oxidase, 1 × 10−4 mol dm−3 solutions of complexes 1−15 were treated with 1 × 10−2 mol dm−3 (100 equiv) of 3,5-DTBC under aerobic condition. The course of the reaction was followed by UV−vis spectroscopy (Figures S20−S32, Supporting Information). The time-dependent UV−vis spectral scan was performed in pure DMSO medium. All complexes behave similarly, showing a smooth conversion of 3,5-DTBC to 3,5-DTBQ. Figure 6a and 6b shows the spectral change for 2 (as representative of tertbutyl complexes) and 14 (as representative of methyl complexes) upon addition of 100-fold 3,5-DTBC (1 × 10−2 M) observed at an interval of 5 min in DMSO medium. The kinetics of the 3,5-DTBC oxidation was determined by monitoring the increase of the product 3,5-DTBQ. The experimental conditions were the same as we reported earlier using ε of 3,5-DTBQ in DMSO = 2100 M−1 cm−1.18 All complexes showed saturation kinetics, and a treatment based on the Michaelis−Menten model seemed to be appropriate. The binding constant (KM), maximum velocity (Vmax), and rate constant for dissociation of substrates (i.e., turnover number, kcat) were calculated for all complexes using the Lineweaver− Burk graph of 1/V vs 1/[S] (Figures S33−S62, Supporting Information) using the eq 1/V = {KM/Vmax}{1/[S]} + 1/Vmax. The calculated values of the kinetic parameters Vmax ranges from 7.483 × 10−7 to 6.858 × 10−6 (standard deviation ranges ±0.0001 × 10−6 to ±0.133 × 10−6), KM ranges from 1.918 × 10−2 to 7.898 × 10−4 (standard deviation ranges ±0.02 × 10−3 to ±0.404 × 10−3), and kcat ranges from 0.501 × 102 to 2.468 × 8264
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
102 (standard deviation ranges from 0.00009 × 102 to 0.09 × 102). Data from Table 7 indicate that the group at the para position has a significant high impact on the catecholase
that the poorly defined redox chemistry of that class of complexes would never allow one to establish such a correlation. In 2002, Belle and co-workers tried to evaluate the effect of the para substituent groups of a ligand system similar to that of Krebs on the catecholase activity and on the electrochemical property of their dicopper(II) complexes. They observed that the catalytic properties of their complexes depend on the electronic effects reflected by their redox potential values, even though the electronic effects did not control solely the reactivity. Here again we are trying to unfold the effect of the substituent groups in the para position of the bridging phenolate group on catecholase and on the electronic property of dicopper(II) complexes. Belle et al. reported that the strongly electron-withdrawing trifluoromethyl group completely inhibits the activity, the poorly electron-withdrawing fluorine atom inhibits to a moderate extent, and the strongly electron-donating methoxyl group increases notably the activity. On the contrary, here we observed that the activity remarkably increases with the electron-withdrawing chloro group whereas it decreases with the electron-donating tert-butyl group considering the corresponding complex with a methyl at the para position as a reference especially for nine complexes, namely, complexes 4, 5, and 12, complexes 6, 7, and 13, and complexes 10, 11, and 15. Among the several possibilities, (i) enzyme−substrate adduct formation and (ii) reduction of Cu(II) to Cu(I) are the possible rate-determining steps. If the former is the key step in catalysis the effect of electrochemical property changes on changing the para substituent groups should not be prominent on catecholase activity. Belle and coworkers explained their observed effect on considering enzyme−substrate adduct formation as the key step in their catalysis. On the other hand, our observation clearly suggests that possibility (ii) should be the key step in our catalysis, especially with complexes 6, 7, and 13, omplexes 10, 11, and
Table 7. kcat Values for the 15 Dinuclear Complexes complex
kcat × 102 (in h−1)
complex
kcat × 102 (in h−1)
1 2 3 4 5 6 7 8
0.737 0.505 0.846 1.272 0.496 2.468 0.507 0.269
9 10 11 12 13 14 15
0.518 0.683 0.501 0.907 0.546 1.178 0.542
activity, and the kcat values for complexes containing L13, L14, and HL22 ligands (namely, complexes 4, 5, and 12), L15, L16, and L23 ligands (namely, complexes 6, 7, and 13) and L110, L112, and H2L111 ligands (namely, complexes 10, 11, and 15) follow the influence of R groups in the order chloro > methyl > t-butyl, in accord with the electron-withdrawing property, i.e., complexes with the most electron-withdrawing chloro group show the highest activity and vice versa. However, a similar trend is not observed in other cases. Correlation between kcat and Electrochemical Property of the Complexes. In 1997, Krebs et al., one of the pioneers of catechol oxidase model study, tried to establish a correlation between catecholase activity and electrochemical property of a series of dicopper(II) complexes of phenol-based Mannich ligands. However, they failed to find any clear relationship between the electrochemical properties and the catecholase activity of the complexes. In addition, they stated
Scheme 3. Catalytic Cycle of Catechol Oxidation by a Generic Dicopper Complex Used in This Studya
a
Two molecules of catechol are oxidized coupled with the reduction of molecular oxygen to water. 8265
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
Figure 7. BP86/def2-TZVP-optimized structures of A−E intermediates of the catalytic cycle. Distances in Angstroms. Hydrogen atoms have been omitted for clarity apart from the OH group of the coordinated catechol to differentiate it from the anionic oxygen atom.
of catechol and one copper metal center, and we indicated this interaction as a dashed line (see B complex in Scheme 3). We focused our study in the first part of the mechanism, i.e., from the met state (A) to the deoxy state (C) and subsequent oxidation to D. We studied a possible correlation between the energy associated with each step and the experimental kcat values. We optimized compounds A, B, and C for the Cl and Me series of complexes in order to see differences between an electron-withdrawing and electron-donating substituent. Since formation of the catechol complex (B) is energetically favorable for all complexes it is expected to be a fast step of the mechanism. In contrast, conversion of compound B to C is unfavorable energetically because both Cu ions are reduced and the coordinated ligands are released (oxidized catechol and a water molecule). Experimental observations (vide supra) suggest that oxidation of 3,5-DTBC is produced via two steps, semibenzoquinone followed by quinone formation with concomitant reduction of Cu(II) to Cu(I). Therefore, conversion of B to C is likely formed by two consecutive one-electron transfer reactions with formation of an intermediate denoted as B′ in Scheme 3. We computed the energetic cost of this process for two series of complexes (R = Cl and Me) in order to investigate if there is a correlation between the energetic cost of this step and the catalytic constant. We have not performed the study for the R = tBu series because four out of five complexes present very similar kcat ≈ 50 h−1, and consequently, it was impracticable to correlate the data. In Figure 7 we represent the optimized geometries for L3 (R = Cl) complexes A−E, as an exemplifying model system. It can be observed that the monocoordinated catechol molecule forms the aforementioned ancillary interaction with the other copper metal center (Cu−OH distance 2.61 Å). This interaction is important to fix the final geometry of the catechol in the
15, and complexes 8, 9, and 14 series, and therefore, we are supposed to get a better correlation between electrochemical property and catecholase activity in those cases. On comparing Tables 6 and 7 another interesting point may be stated that the E° value is proportional to the kcat value for certain ligand series complexes, e.g., complexes 1−3 and complexes 8 and 9. However, in other cases we cannot see any direct relation between the electrochemical property and the catecholase activity of the complexes, and the reason is not at all clear to us at this moment. However, Neves et al. in one of their classical works tried to establish a relation between the kinetic parameters and the electrochemical properties. They observed a very good relation between kinetic parameter k2 (secondorder rate constant = kcat/KM) and (ΔE)1,2 (=E(red)1 − E(red)2), although no good correlation was noticed between kcat and E°. Theoretical Study. From the above discussion it is clear that we do not get a very linear relationship between E° and kcat all the time. As stated earlier, redox participation of two copper centers appears to be crucial in catecholase activity. Thus, we do believe that the kcat value should be related to a certain factor associated with electrochemical change of the catalyst in the catalytic cycle, and that very idea prompted us to perform a detailed DFT theoretical study especially on the reaction pathway involved in this catalysis. It should be mentioned that the E° values are computed in the original complexes, and the reduction of Cu(II) to Cu(I) in the catalytic cycle occurs in a species where the catechol is coordinated to the complex; consequently, the E° value of this species may differ considerably from that of the uncoordinated complex. To perform the theoretical study we used the catalytic cycle proposed by Krebs and co-workers36 and redrawn by Reedijk and co-workers.37 We observed in the optimized catechol complexes an ancillary interaction between the phenolic group 8266
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
complex and likely facilitates the electron transfer from the catechol to the Cu(II) ions. It is interesting to note the optimized geometry of compound D, where the copper ions have been reoxidized, the O−O distance is 1.42 Å, and all Cu− O distances are very similar, ranging from 2.11 to 2.13 Å. We also studied the subsequent complexation of catechol to form complex E. The global transformation of C to E that involves incorporation of molecular oxygen, oxidation of the Cu(I) ions, and complexation of catechol is energetically very favorable (−68.9 kcal/mol) for the system shown in Figure 8.
Figure 9. Spin density (left) and SOMOs of intermediate B′.
onto the atoms of the ligands directly bonded to the metal centers and the catechol ring. The SOMO shows the contribution of the π system of the semiquinonate moiety to this orbital with partial contribution of the d atomic orbitals of the metal centers. In contrast, the catechol does not participate in the SOMO-1, which has the contribution of the atomic orbitals of the metal centers and the atoms directly bonded to them. In Table 8 we summarize the energetic cost on going from intermediate B to C for the R = Cl and R = Me series and the Table 8. Reaction Energy Necessary for Conversion of B to C
Figure 8. Partial views of the X-ray structures of GUGPUT, NAPGEQ, QOSRAS, and the optimized intermediate D for comparison purposes. Hydrogen atoms have been omitted for clarity. Distances in Angstroms.
complex
ΔE (kcal/mol)
ΔErel (kcal/mol)
kcat × 102 (h−1)
1 3 4 6 8 10 12 13 14 15
38.2 42.3 36.1 26.9 42.1 39.5 39.2 40.4 37.6 46.7
11.3 4.6 9.2 0.0 15.1 12.6 1.6 2.8 0.0 9.0
0.7372 0.8451 1.272 2.45 0.2690 0.6827 0.9090 0.754 1.1681 0.5421
experimental catalytic constant values for comparison purposes. For each series we also included the relative energies with respect to the most favorable conversion. It should be mentioned that the energies have been computed taking into account solvation effects by means of a continuum model. As aforementioned, this process is unfavorable due to the loss of two coordinating ligands in this mechanistic step. We are interested in the variation of the energy required to oxidize catechol and concomitant reduction of Cu(II) to Cu(I) depending on the substituents of the ligands. We represented in Figure 10 the relative energy values versus the experimental kcat values resulting in interesting correlations. As it can be observed we found a very good correlation between the relative energy and kcat (R = 0.995) for the R = Cl series that strongly support the hypothesis that the key mechanistic step is the reduction of Cu(II) to Cu(I). It is worth mentioning that such correlation is not obtained using the experimental E° reduction potential values. This likely indicates that the coordination of the catechol to the complex influences the reduction ability of the copper ions due to a combination of the coordination itself and steric effects. The same representation for the R = Me series also shows a rather good correlation between the relative energy and the experimental kcat values (R = 0.903). If both series are gathered in the same representation the regression coefficient is also good (R = 0.912), which further confirms that there is a clear correlation between the kcat values and the Cu(II) to Cu(I) reduction, which does not necessarily
It should be mentioned that we studied theoretically and experimentally other orientations and conformations for the intermediate D. In Figure S63, Supporting Information, we show three different orientations that we studied in combination to searches in the Cambridge Structural Database (CSD) to give experimental support to the DFT calculations. In Figure 8 we show several X-ray structures that exhibit the coordination mode of O2 proposed in the intermediate D. In particular, (μ2-(1,2-bis(2-(bis(6-methyl-2-pyridyl)methyl)-6pyridyl)ethane-N,N′,N″))-(μ2-peroxo)-dicopper (GUGPUT)38 and (μ2-peroxo)-bis(1,4,7-tris(isopropyl)-1,4,7-triazacyclodecane)-dicopper (NAPGEQ)39 structures do not present additional oxo-bridging ligands in the structure. Conversely, it is remarkable that the geometry of the (μ2-benzoato)-(μ2peroxo)-bis(α-isosparteine)-dicopper (QOSRAS)40 structure is similar to the optimized geometry of intermediate D (see Figure 8, bottom right) concerning bond distances and angles, keeping in mind that the two copper centers are bridged by a carboxylate instead of phenolate (the Cu−O2 unit is less butterfly distorted). We also analyzed the spin density and single-occupied molecular orbitals (SOMOs) of intermediate B′ that are shown in Figure 9. In this electronic configuration one electron has been transferred from the catechol to one metal center. The spin density plot simply reflects the sum of both SOMOs, and it can be clearly observed that some spin density is delocalized 8267
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
Figure 10. Regression plots of kcat versus the relative energy cost for conversion of B to C (see Scheme 3).
■
ACKNOWLEDGMENTS The authors thank CSIR, New Delhi [01(2464)/11/EMR-II dt 16-05-11 to DD], for financial support and the University of Calcutta for providing the facility of the single-crystal X-ray diffractometer from the DST FIST program. A.B. and A.F. acknowledge the Spanish Ministerio de Economiá y Competitividad (MINECO) (Projects, and CONSOLIDER INGENIO 2010, CSD2010-00065, FEDER funds) and the Direcció General de Recerca i Innovació del Govern Balear (project 23/2011, FEDER funds) for financial support. We thank the CTI (UIB) for computational facilities. This paper is dedicated to the memory of Pali Maiti, whom we lost very untimely and unfortunately.
correspond to the electrochemical properties of the free complexes.
■
CONCLUSIONS
In order to explore the deep insight on the correlation between the electrochemical behavior and the catecholase activity of synthetic analogues of catechol oxidase, a nearly unanswered corner of catecholase activity, we designed and synthesized 15 closely related dicopper(II) complexes of phenol-based compartmental ligands. Eleven of them are newly introduced, and four were reported earlier. The 11 new complexes have been characterized by routine physicochemical techniques, and complexes 1−3 and 8−10 have further been structurally characterized by X-ray single-crystal structural analyses. Electrochemical and catecholase activities of all complexes have been thoroughly investigated in DMSO medium. CV and DPV analyses reveal that oxidation of 3,5-DTBC catalyzed by dicopper(II) complexes proceed via two steps: first semibenzoquinone followed by benzoquinone with concomitant reduction of CuII to CuI. Our critical analysis reveals that apparently there is nearly no linear relationship between kcat and E° values of the complexes. That very observation prompted us to calculate the energetics associated with each step of the catalysis, and interestingly, we were successful to establish for the first time quite a good linear relationship between the energetic cost for the catechol-bound CuII to CuI reduction and kcat values.
■
■
(1) (a) Eicken, C.; Zippel, F.; Buldt-Karentzopoulos, K.; Krebs, B. FEBS Lett. 1998, 436, 293. (b) Than, R.; Feldmann, A. A.; Krebs, B. Coord. Chem. Rev. 1999, 182, 211. (c) Gentschev, P.; M?ller, N.; Krebs, B. Inorg. Chim. Acta 2000, 300, 442. (d) Merkel, M.; Mçller, N.; Piacenza, M.; Grimme, S.; Rompel, A.; Krebs, B. Chem.Eur. J. 2005, 11, 1201. (2) Monzani, E.; Quinti, L.; Perotti, A.; Casella, L.; Gullotti, M.; Randaccio, L.; Geremia, S.; Nardin, G.; Faleschini, P.; Tabbi, G. Inorg. Chem. 1998, 37, 553. (3) Manzur, J.; Garcia, A. M.; Vega, A.; Spodine, E. Polyhedron 1999, 18, 2399. (4) (a) Torelli, S.; Belle, C.; Gautier-Luneau, I.; Pierre, J. L.; SaintAman, E.; Latour, J. M.; Le Pape, L.; Luneau, D. Inorg. Chem. 2000, 39, 3526. (b) Belle, C.; Selmeczi, K.; Torelli, S.; Pierre, J.-L. C. R. Chim. 2007, 10, 271. (5) Gupta, M.; Mathur, P.; Butcher, R. J. Inorg. Chem. 2001, 40, 878. (6) Murthy, N. N.; Mahroof-Tahir, M.; Karlin, K. D. Inorg. Chem. 2001, 40, 628. (7) Kao, C. H.; Wei, H. H.; Liu, Y. H.; Lee, G. H.; Wang, Y.; Lee, C. J. J. Inorg. Biochem. 2001, 84, 171. (8) Fernandes, C.; Neves, A.; Bortoluzzi, A. J.; Mangrich, A. S.; Rentschler, E.; Szpoganicz, B.; Schwingel, E. Inorg. Chim. Acta 2001, 320, 12. (9) Mukherjee, J.; Mukherjee, R. Inorg. Chim. Acta 2002, 337, 429. (10) Neves, A.; Rossi, L. M.; Bortoluzzi, A. J.; Szpoganicz, B.; Wiezbicki, C.; Schwingel, E.; Haase, W.; Ostrovsky, S. Inorg. Chem. 2002, 41, 1788. (11) (a) Ackermann, J.; Meyer, F.; Kaifer, E.; Pritzkow, H. Chem. Eur. J. 2002, 8, 247. (b) Ackermann, J.; Buchler, S.; Meyer, F. C. R. Chim. 2007, 10, 421. (12) Seneque, O.; Campion, M.; Douziech, B.; Giorgi, M.; Riviere, E.; Journaux, Y.; Mest, Y. L.; Reinaud, O. Eur. J. Inorg. Chem. 2002, 2007. (13) (a) Koval, I. A.; Pursche, D.; Stassen, A. F.; Gamez, P.; Krebs, B.; Reedjik, J. Eur. J. Inorg. Chem. 2003, 1669. (b) Koval, I. A.; Huisman, M.; Stassen, A. F.; Gamez, P.; Lutz, M.; Spek, A. L.; Reedijk, J. Eur. J. Inorg. Chem. 2004, 591.
ASSOCIATED CONTENT
S Supporting Information *
X-ray crystallographic data in CIF format of complexes 1−3 and 8−10, FT-IR spectra, wavelength scans for catecholase activity, kinetic plots, ORTEP views of complexes 2 and 10, cyclic voltammograms, and differential pulse voltammograms. This material is available free of charge via the Internet at http://pubs.acs.org.
■
AUTHOR INFORMATION
Corresponding Authors
*E-mail: *E-mail: *E-mail: *E-mail:
REFERENCES
[email protected].
[email protected].
[email protected].
[email protected] .
Notes
The authors declare no competing financial interest. 8268
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269
Inorganic Chemistry
Article
(14) Sreenivasulu, B.; Vetrichelvan, M.; Zhao, F.; Gao, S.; Vittal, J. J. Eur. J. Inorg. Chem. 2005, 4635. (15) Wegner, R.; Gottschaldt, M.; Wolfgang, P.; Jager, E.-G.; Klemm, D. J. Mol. Catal. A 2003, 201, 93. (16) Thirumavalavan, M.; Akilan, P.; Kandaswamy, M.; Chinnakali, K.; Kumar, G. S.; Fun, H. K. Inorg. Chem. 2003, 42, 3308. (17) Kaizer, J.; Csonka, R.; Speir, G.; Giorgi, M.; Reglier, M. J. Mol. Catal. A 2005, 235, 81. (18) Banu, K. S.; Chattopadhyay, T.; Banerjee, A.; Bhattacharya, S.; Suresh, E.; Nethaji, M.; Zangrando, E.; Das, D. Inorg. Chem. 2008, 47, 7083. (19) (a) Guha, A.; Chattopadhyay, T.; Paul, N. D.; Mukherjee, M.; Goswami, S.; Mondal, T. K.; Zangrando, E.; Das, D. Inorg. Chem. 2012, 51, 8750. (b) Adhikary, J.; Chakraborty, P.; Das, S.; Chattopadhyay, T.; Bauza, A.; Chattopadhyay, S. K.; Ghosh, B.; Mautner, F. A.; Frontera, A.; Das, D. Inorg. Chem. 2013, 52, 13442. (c) Ghosh, T.; Adhikary, J.; Chakraborty, P.; Sukul, P. K.; Jana, M. S.; Mondal, T. K.; Zangrando, E.; Das, D. Dalton Trans. 2014, 43, 841. (20) Reim, J.; Krebs, B. J. Chem. Soc., Dalton Trans. 1997, 3793. (21) Torelli, S.; Belle, C.; Hamman, S.; Pierre, J.-L. Inorg. Chem. 2002, 41, 3983. (22) Meiwes, D.; Ross, B.; Kiesshauer, M.; Cammann, K.; Witzel, H.; Knoll, M.; Borchardt, M.; Sandermaier, C. Lab. Med. 1992, 15, 24. (23) Sreenivasulu, B.; Zhao, F.; Gao, S.; Vittal, J. J. Eur. J. Inorg. Chem. 2006, 2656. (24) Saha, S.; Koner, S.; Tuchagues, J.-P.; Boudalis, A. K.; Okamoto, K.-I.; Banerjee, S.; Mal, D. Inorg. Chem. 2005, 44, 6379. (25) Gagne, R. R.; Spiro, C. L.; Smith, T. J.; Hamann, C. A.; Thies, W. R.; Shiemeke, A. K. J. Am. Chem. Soc. 1981, 103, 4073. (26) SMART, SAINT, Software Reference Manual; Bruker AXS Inc.: Madison, WI, 2000. (27) Sheldrick, G. M. Acta Crystallogr. 2008, 64, 112. (28) Farrugia, L. J. J. Appl. Crystallogr. 2012, 45, 849. (29) Ahlrichs, R.; Bar, M.; Hacer, M.; Horn, H.; Komel, C. Chem.Phys.Lett. 1989, 162, 165. (30) Klamt, A.; Schuurmann, G. J. Chem. Soc., Perkin Trans. 2 1993, 799. (31) Nakamoto, K. Infrared and Raman Spectra of Inorganic and Coordination Compounds, 3rd ed.; Wiley: New York, 1978. (32) Lever, A. B. P. Inorganic Electronic Spectroscopy; Elsevier: Amsterdam, The Netherlands, 1984; p 553. (33) (a) Zanello, Inorganic Electrochemistry: Theory, Practice and Applications; The Royal society of Chemistry, U. K., 2003; pp110− 113. (b) Richardson, D. E.; Taube, H. Inorg. Chem. 1981, 20, 1278. (34) McDaniel, D. H.; Brown, H. C. J. Org. Chem. 1958, 23, 420. (35) Harmalker, S.; Jones, S. E.; Sawyer, D. T. Inorg. Chem. 1983, 22, 2790. (36) Gerdemann, C.; Eicken, C.; Krebs, B. Acc. Chem. Res. 2002, 35, 183. (37) Koval, I. A.; Belle, C.; Selmeczi, K.; Philouze, C.; Saint-Aman, E.; Schuitema, A. M.; Gamez, P.; Pierre, J. L.; Reedijk, J. J. Biol. Inorg. Chem. 2005, 10, 739. (38) Kodera, M.; Katayama, K.; Tachi, Y.; Kano, K.; Hirota, S.; Fujinami, S.; Suzuki, M. J. Am. Chem. Soc. 1999, 121, 11006. (39) Lam, B. M. T.; Halfen, J. A.; Young Junior, V. G.; Hagadorn, J. R.; Holland, P. L.; Lledos, A.; Cucurull-Sanchez, L.; Novoa, J. J.; Alvarez, S.; Tolman, W. B. Inorg.Chem. 2000, 39, 4059. (40) Funahashi, Y.; Nishikawa, T.; Wasada-Tsutsui, Y.; Kajita, Y.; Yamaguchi, S.; Arii, H.; Ozawa, T.; Jitsukawa, K.; Tosha, T.; Hirota, S.; Kitagawa, T.; Masuda, H. J. Am. Chem. Soc. 2008, 130, 16444.
8269
dx.doi.org/10.1021/ic5005177 | Inorg. Chem. 2014, 53, 8257−8269