Next Article in Journal
Effect of Fiber Content on the Mechanical Properties of Engineered Cementitious Composites with Recycled Fine Aggregate from Clay Brick
Next Article in Special Issue
Spectroscopic, Thermal, Microbiological, and Antioxidant Study of Alkali Metal 2-Hydroxyphenylacetates
Previous Article in Journal
Breakage Strength of Wood Sawdust Pellets: Measurements and Modelling
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Materials
Volume 14
Issue 12
10.3390/ma14123274
materials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Synthesis, Spectroscopic, and Theoretical Study of Copper and Cobalt Complexes with Dacarbazine

by
Grzegorz Świderski
1,*,
Ryszard Łaźny
2,
Michał Sienkiewicz
2,
Monika Kalinowska
1,
Renata Świsłocka
1,
Ali Osman Acar
3,
Aleksandra Golonko
4,
Marzena Matejczyk
1 and
Włodzimierz Lewandowski
1,*
1
Department of Chemistry, Biology and Biotechnology, Bialystok University of Technology, Wiejska Street 45E, 15-351 Bialystok, Poland
2
Faculty of Chemistry, University of Bialystok, Ciolkowskiego Street 1K, 15-245 Bialystok, Poland
3
Material Science and Nanotechnology Engineering, TOBB University of Economics and Technology, 06-560 Ankara, Turkey
4
Institute of Agricultural and Food Biotechnology, Department of Microbiology, Rakowiecka 36, 02-532 Warsaw, Poland
*
Authors to whom correspondence should be addressed.
Materials 2021, 14(12), 3274; https://doi.org/10.3390/ma14123274
Submission received: 16 April 2021 / Revised: 7 June 2021 / Accepted: 10 June 2021 / Published: 13 June 2021
(This article belongs to the Special Issue Comparison of Structural, Spectroscopic, Theoretical and Thermal Properties of Metal Complexes and Their Application)
Browse Figures
Versions Notes

Abstract

:
Dacarbazine (DAC) 5-(3,3-dimethyl-1-triazenyl)imidazole-4-carboxamide is an imidazole-carboxamide derivative that is structurally related to purines. DAC belongs to the triazene compounds, which are a group of alkylating agents with antitumor and mutagenic properties. DAC is a non-cell cycle specific drug, active in all phases of the cellular cycle. In the frame of this work the 3d metal complexes (cobalt and copper) with dacarbazine were synthesized. Their spectroscopic properties by the use of FT-IR, FT-Raman, and 1HNMR were studied. The structures of dacarbazine and its complexes with copper(II) and cobalt(II) were calculated using DFT methods. The effect of metals on the electronic charge distribution of dacarbazine was discussed on the basis of calculated NBO atomic charges. The reactivity of metal complexes in relation to ligand alone was estimated on the basis of calculated energy of HOMO and LUMO orbitals. The aromaticity of the imidazole ring in dacarbazine and the complexes were compared (on the basis of calculated geometric indices of aromaticity). Thermal stability of the investigated 3d-metal complexes with dacarbazine and the products of their thermal decomposition were analyzed.

1. Introduction

Dacarbazine (DAC) [5-(3,3-Dimethyl-1-triazenyl) imidazole-4-carboxamide] belongs to the alkylating agents of triazene group use in anticancer therapy [1,2,3,4]. DAC is an imidazole-carboxamide derivative that is structurally related to purines. DAC is a non-cell cycle specific drug, active in all phases of the cellular cycle. This drug is used for chemotherapy in different types of cancer malignant melanoma, soft-tissue sarcoma, osteogenic sarcoma, neuroblastomas, and Hodgkin’s disease [5,6]. The detailed molecular mechanism of DAC activity is unclear. Research to date indicates three possible mechanisms of action of this drug: (a) inhibition of DNA synthesis by acting as a purine analog, (b) action as alkylating agent, and (c) interaction with SH groups [7]. Hayward et al. [8], proposed another mechanism of DAC activity based on antimetabolite-like activity. In this mechanism, DAC plays the role of an inhibitor of purinic base incorporation into DNA during DNA synthesis. However, this mechanism does not seem to play a key role in the antitumor activity of DAC [8].
Based on experimental data it seems that alkylation of nucleic acid is the primary mode of antitumor DAC action [7]. DAC as an alkylator donates methyl groups to many sites in DNA. A very common place of attack of DAC is the O(6)-position of guanine resulting in O(6)-methylguanine (O(6)-MeG) which frequently mispairs with thymine during DNA duplication. Each cell is equipped with repair systems, the mismatch repair (MMR) system recognizes O(6)-MeG:T mismatches which lead to cell cycle arrest and cell death. If MMR is defective, the cell will divide and GC → AT transition mutations occur (Figure 1) [5,6]. Generally, DAC causes DNA damage, cell cycle arrest and apoptosis [9]. Hayward et al. [8] proposed another mechanism of DAC activity based on antimetabolite-like activity. In this mechanism, DAC plays the role of an inhibitor of purinic base incorporation into DNA during DNA synthesis. However, this mechanism does not seem to play a key role in the antitumor activity of DAC [8]. DAC is a prodrug that is activated by N-demethylation in liver microsomes with participation of the liver’s microsomal enzymes (CYP1A1, CYP1A2, and CYP2E1) that hydroxilate and N-demethylate the DAC [7].
Studies by earlier authors have shown that DAC, in combination with other drugs is used for treating renal adenocarcinoma, soft-tissue sarcoma and malignant lymphomas. Al-Qatati and Aliwaini [9] in their work, presented that combined pitavastatin with DAC treatment induces cell death through intrinsic apoptosis in melanoma cells. Naserian et al. [10] investigated the cytotoxicity and apoptosis of the metformin alone and in combination with DAC in Raji and Ramos lymphoma cell lines. Metformin showed synergistic cytotoxic effects in combination with DAC, reduced cell viability, and increased apoptosis in Raji and Ramos lymphoma cells in comparison with the use of each drug alone. Additionally, Finotello et al. [11] reported that results of in vitro and in vivo studies suggest that dacarbazine acts synergistically with anthracyclines and has a moderate effect on the treatment of high-grade sarcomas in humans.
The greatest obstacle in the treatment of patients with advanced stages of melanoma and other cancers is the unique tumor resistance, both primary and secondary, to all known compounds including dacarbazine. The complexation of metal ions with anticancer drugs creates the possibility of an increase in the bioavailability and activity of drugs [12,13]. Metal complexes play a very important role in modern cancer therapy. The complexation of metal ions with biologically active ligands can change their therapeutic activity. Metal ions affect the distribution of the electronic system of ligands which in turn changes the reactivity and biological activity of the ligands [14,15]. Complexes with some of the transition metals, e.g., Au, Ag, Co, Cu, Ni, Fe Pt(II) and (IV), Pd(II), and Ru(III) show high anticancer activity [16].
DAC possesses a triazene group connected to the imidazole ring which is attached to the carboxamide group. In the structure of the DAC, they are atoms that may be donors of electron pairs, making complexation of the ligand with metal ions possible. Different types of metal-ligand complexes can be formed when the molecule of dacarbazine coordinates metal through seven potential donor atoms.
Several studies have shown that complexation of dacarbazine with metal ions may change the reactivity and stability of the molecule, thereby improving its therapeutic properties. In vitro and in vivo studies have shown that the iron(III) complex of dacarbazine may exhibit an increased antitumor effect compared to the ligand alone [17,18]. Temerk and Ibrahim [19] in their work investigated the influence of DAC complexation with Cu(II) on the drug intercalation to DNA. The results of that experiment showed that the interaction of DAC with dsDNA in the presence of Cu(II) led to a much stronger intercalation, compared to DAC alone [19].
Our studies showed that there is a dependency between the electronic charge distribution in metal complex and salt molecules and the location of these metals in the periodic table [20]. Alkali metal cations and some heavy metal cations (Ag(I), Pb(II), Hg(I), Hg(II)) perturb the electronic system of ligands (benzoic, salicylic, pyridinecarboxylic acids and others), whereas 3d and 4f metal cations stabilize it [21,22,23,24,25,26]. This conclusion allows us to foresee the changes in the electronic structure under the influence of metal cations and estimation of the physicochemical and biological properties of complexes. Thermal analysis allows us to determine the thermal stability of the investigated metal complexes and to determine the residues of this process. We have repeatedly tested the thermal stability of ligands of high biological importance and their complexes with metals. We found that metal oxides are the final product of the decomposition of metal complexes under the oxygen conditions of the process, with the use of a sufficiently high temperature [27,28,29,30,31]. This may be of practical importance when thermal disposal of drugs based on metal complexes with particularly toxic ligands, such as dacarbazine, is used.
In the frame of this work the 3d metal complexes with dacarbazine were synthesized. The aim of the work was to estimate the effect of metal ions (copper(II) and cobalt(II)) on the reactivity and electronic stability of dacarbazine. For this purpose, the spectroscopic properties of compounds were studied by the use of FT-IR, FT-Raman, and 1HNMR. Thermal stability of the investigated 3d-metal complexes (cobalt(II) and copper(II)) with dacarbazine and the products of their thermal decomposition were analyzed.
Synthesis of complexes with metals such as copper and cobalt is a noticeable trend in many fields of science, such as pharmaceuticals and technologies related to catalysis and sensors. Padnya et al., with high efficiency, synthesized a series of catechol-containing Schiff bases, tetrasubstituted on the lower rim thiacalix [32] arene derivatives in three stereoisomeric forms, cone, partial cone and 1,3-alternate, which can be used as a component of antimicrobial agents, catalysts, chemical, and electrochemical sensors [32]. On the other hand, the copper complexes described by Colombo et al. can be applied as eco-sustainable redox mediators in sye-sensitized solar cells [33].
In the case of compounds with potential biological applications, it should be noted that some metal complexes with many oxidation states, e.g., vanadium, undergo many chemical transformations under physiological conditions. So the study of the mechanism of action should assume that the complexes do not exist in cells/tissues in its original form. The authors of the review suggest the necessity of speciation of the vanadium species formed in each biological compartment and defining the defined reaction conditions under which the complex compound acts [34]. This is especially important when assessing the interaction of such complexes with biological macromolecules, e.g., proteins or DNA. Currently, there are many inorganic coordination compounds used as antibiotics (bleomycin), in the treatment of neoplastic diseases (cisplatin) or anti-arthritis drugs (auranofin) [35]. Metal complexes also offer hope for a more effective treatment of solid tumors in which there is a hypoxic state with a more acidic and reducing environment. Studies on compounds containing Pt4+, Ru3+, and Co3+ have shown the potential for selective action in cancer cells. Cobalt, unlike ruthenium or platinum, is non-toxic and, when reduced from Co3+ to Co2+ under hypoxic conditions, the complex becomes labile, which enables drug dissociation. In experimental research, de Souza et al. proposed cobalt(III) L-Phe complexes as a model for the hypoxia-activated delivery of melphalan [36].

2. Materials and Methods

2.1. Synthesis

Complexes dacarbazine-CoCl2 (DAC)2*1.75CH3OH and CuCl2(DAC)2*1.5CH3OH were obtained by the same method of synthesis by reaction dacarbazine of metal chlorides in a methanolic medium. The weight amount of dacarbazine (0.2 mol) was dissolved in cold methanol and it was added to 5 mL of methanolic solution of metal chloride (0.1 mol) (CuCl2, CoCl2). Then the mixture was shaken for 2 h in the water bath in room temperature. Next, 50 mL of diethyl ether was added and it was left for 2 h. The precipitate was washed three times with 50 mL of diethyl ether. The obtained complexes were dried in a vacuum for 72 h.
Elemental analysis: Anal. Calc. for CoCl2(DAC)2*1.75CH3OH: C, 29.62; H, 4.75; N, 30.72. Found: C, 29.35; H, 5.03; N, 30.318. Calc. for CuCl2(DAC)2*1.5CH3OH: C, 29.98; H, 4.91; N, 30.53. Found: C, 30.13; H, 5.30; N, 30.82.

2.2. Methods

The FT-IR spectra were recorded with an Alfa (Bruker) spectrometer (Billerica, MA, USA) within the range of 400–4000 cm−1. Samples in the solid state were measured in KBr matrix pellets and ATR technique. FT-Raman spectra of solid samples were recorded in the range of 100–4000 cm−1 with a MultiRam (Bruker) spectrometer (Billerica, MA, USA). The resolution of the spectrometer was 1 cm−1. The 1H spectra of DMSO solution of studied compounds were recorded with a Bruker Avance II 400 MHz (Billerica, MA, USA) unit at room temperature (sample content 5 mg/1 mL of solvent). TMS was used as an internal reference. Elemental analysis for the weight percentages of carbon and hydrogen was done with Perkin-Elmer 240 (Waltham, MA, USA) equipment. To calculate optimized geometrical structures, NMR and IR spectra, NBO analysis, HOMO/LUMO orbitals of dacarbazine, and complexes with copper and cobalt were used with the density functional (DFT) hybrid method B3LYP with nonlocal correlation provided by the Lee–Young–Parr expression. All calculations were carried out with 6-311++G(d,p) basis set and performed with GAUSSIAN 09 [37] packed.
The aromaticity indices (HOMA, GEO, EN, I6) were calculated for geometric structures (theoretical and calculated) of dacarbazine and its complexes. The HOMA index (harmonic oscillator model of aromaticity) differs from all other geometry-based ones by assuming another reference bond length. In this model, instead of the mean bond length a concept of the optimal bond length is applied [38]:
H O M A = 1 − [ α R o p t − R a v 2 + α n ∑ R a v − R i 2 ] = 1 − E N − G E O
Within the confines of the HOMA model, it is possible to obtain two components which describe different contributions to decrease in aromaticity, i.e., (a) due to bond elongation (the EN component), and (b) due to bond length alternation (the GEO component). The value of HOMA index is equal to 1 for the entire aromatic system; HOMA = 0 when the structure is non-aromatic and HOMA < 0 for an anti-aromatic ring.
The value of the Bird’s aromaticity index (I5, I6) describes the equation [39]:
I = 100 1 − V V k
where: Vk is for the five-membered rings 35 and the six-membered 33.3, and V is calculated from the equation:
V = 100 n a v [ ∑ r = 1 n n r − n a v 2 / n ] 1 / 2
where: nav—average binding order, n—bond order based on bond length: n = (a/R) − b, a and b—parameters depending on the type of atoms in the bond.
NBO analysis was performed for the optimized structures to determine the electronic charge distribution [40]. Calculations were made using the B3LYP/6-311++G(d, p) method.
Thermal analyses of the prepared complexes were performed by the thermogravimetric (TG) methods using the Perkin-Elmer analyzer in a dynamic air atmosphere. Therefore, 6.0–8.0 mg samples were heated in the range of 50–890 °C in the ceramic crucibles using a heating rate of 10 °C min−1. The products of dehydration and decomposition processes were determined from the TG curves.
The crystal structure of the DNA dodecamer was obtained from Protein Database (PDB ID: 1BDNA) and used for molecular modeling studies carried out in AutoDock 4.2 MGL Tools package supplied with AutoGrid4.0 and AutoDock4.0 [41]. In order to analyze the obtained results and compare the number of hydrogen bonds formed, the VMD (Visual Molecular Dynamics) program was used [42].
AutoGrid4.0 was used to compute grid maps using a grid box. The molecular docking was carried out by setting the grid box size to cover the predicted binding sites, using 52, 56, 118  x, y, z points with a grid spacing of 0.375 Å. The grid center was set to 14.72, 21.006, 8.801 x, y, and z dimensions, respectively. The Lamarckian genetic algorithm (LGA) was selected to generate the best ligand conformers. We successfully performed 1000 docking runs with AutoDock4.0 for each ligand (DAC1, DAC2) and Co(II)-DAC, Cu(II)-DAC complexes. As standard values for hydrogen bond formation assumed distance 2.9 Å and cutoff angle 60 degrees (120 to 180 degrees). To visualize results BIOVIA Discovery Studio and USCF Chimera 1.10.2 were used [43,44].

3. Results and Discussion

3.1. IR and Raman Spectra

The IR and Raman spectra of dacarbazine and IR spectra of cobalt(II) and copper(II) complexes of dacarbazine are shown in Figure 2.
The wavenumbers and intensities of selected bands from the experimental and theoretical IR and Raman spectra of ligand and experimental IR spectra of complexes are gathered in Table 1. The assignment of bands was done on the basis of theoretical wavenumbers obtained at the B3LYP/6-311++G** level as well as literature data [45]. The theoretical wavenumbers were scaled to reproduce adequately the experimental wavenumbers (f = 0.967). Some characteristic bands occurred in the experimental spectra of dacarbazine, i.e., bands assigned to the stretching vibrations of the carbonyl group νC=O at 1609 cm−1 (IR) and 1604 cm−1 (Raman). The theoretical wavenumbers of the νC=O band were 1691 cm−1 (for the DAC1 conformer) and 1683 cm−1 (in the case of the DAC2 conformer) (Figure 3) (the values scaled by the factor f = 0.967). In the IR spectra of Co(II) and Cu(II) complexes of dacarbazine the νC=O band was of similar wavenumbers and reduced intensities compared with the spectra of ligand (i.e., 1605 cm−1—Co(II) complex and 1609 cm−1—Cu(II) complex). This suggested that the carbonyl group did participate in metal ion coordination. In the spectra of dacarbazine and its metal complexes, the characteristic bands assigned to the vibrations of NH2 (of amide group) occurred. There were bands of the stretching asymmetric vibrations νasNH2 located at 3383 cm−1 (IR) and 3371 cm−1 (Raman) as well as stretching symmetric vibration νsNH2 near 3269 cm−1 (IR). In the spectra of complexes, the bands assigned to the stretching vibrations νasNH2 were significantly shifted whereas bands of νasNH2 disappeared compared with the appropriate bands in the spectra of ligand. Significant changes were also observed in the case of deforming out-of-plane vibrations of NH2 group (ρNH2), i.e., 1436 cm−1 in the spectrum of ligand, 1440 cm−1 for Cu and Co complexes as well as 542 cm−1 for ligand and 565 cm−1 for the Co complex and 568 cm−1 for the Cu complex. The metal complexation affected the location of the band derived from the vibrations of the C-NH2 group as well. Namely, the band of deforming out-of-plane vibrations ρC-NH2 was located at 1476 cm−1 and underwent movement to the 1487 cm−1 and 1488 cm−1 in the cobalt and copper complexes of dacarbazine, respectively. The characteristic changes in the wavenumbers of the bands assigned to the vibrations of N=N–N and CH3 groups were observed.
The band assigned to the deformations of the triazene group αNNN was located at 630 cm−1 in the spectra of ligand, and then was shifted toward higher wavenumbers in the spectra of complexes (649 cm−1—Cu complex and 646 cm−1—Co complex). The symmetric stretching vibrations of the methyl group νsCH3, which were present in the spectra of ligand at 3147, 2946, 2753, and 2612 cm−1 disappeared in the spectra of metal complexes. Whereas band of asymmetric stretching vibrations of the methyl group νsCH3 were significantly shifted in the spectra of complexes compared with the spectra of ligand, i.e., 2905 cm−1 in the IR spectra of dacarbazine and 2925 cm−1 and 2923 cm−1 in the IR spectra of Co and Cu complexes, respectively. The band of the deforming in-plane vibrations of the methyl group of dacarbazine disappeared or decreased/increased in their wavenumbers in the spectra of Cu and Co complexes (Table 1).
Moreover, the coordination of copper(II) and cobalt(II) by dacarbazine caused changes in the spectra of complexes in the region of bands assigned to the stretching (νring) and deforming (ringdef) vibrations of the ring (compared with the appropriate bands in the spectra of dacarbazine). Namely, νring bands located at 1381, 1270, 1231 cm−1, and ringdef bands near 882 and 450 cm−1 disappeared in the spectra of metal complexes. Whereas νring bands situated in the spectra of ligand at 1344 and 1402 cm−1 were slightly shifted to the higher wavenumbers in the spectra of complexes. The band of stretching vibrations of the CH bond of the aromatic ring of dacarbazine located at 3174 cm−1 was shifted toward higher wavenumbers in the spectra of complexes (i.e., 3188 cm−1—Co(II) complex and 1386 cm−1—Cu(II) complex). The deforming vibrations of the CH bond of the ring slightly moved to the higher wavenumbers in the spectra of complexes compared with the spectra of ligand.
The difference between the location of the appropriate bands in the spectra of dacarbazine and copper(II) and cobalt(II) complexes of dacarbazine might be caused by the coordination of metal ions by the oxygen atom of carbonyl group and nitrogen atom of imidazole ring.

3.2. NMR Spectra

In the experimental 1HNMR spectra of dacarbazine, four groups signals from protons were observed (Table 2 and Figure S1). The values of the chemical shift for two protons of the amide group have been values 7.41 and 7.29 ppm. This signal possessed very low intensity in the spectra of the cobalt complex, and was shifted to 6.20, 5.40, giving extended signals, due to changes in the electronic density in the amide group, due to the attachment of a metal to this group. In the spectra of ligand, the signals at 3.13–3.50 ppm were assigned to the protons of the methyl group.
Complexation of dacarbazine by the cobalt caused insignificant displacement of this signal. On the basis of the value of signals deriving from the aromatic protons (imidazole ring), the effect of metal ions on the electronic charge distribution of ligand might be discussed. The movement of the chemical shifts of aromatic protons toward higher wavenumbers showed an increase in the aromaticity of the ring as well as the stabilization of the aromatic system [20]. In the case of cobalt(II) complexes of dacarbazine, an increase in the value of the chemical shift of protons of CH bond of the aromatic ring was observed (Figure S1). Proton NMR spectrum for the copper complex could not be recorded. The theoretically calculated chemical shifts in the proton spectrum of the copper complex show similar trends in proton shifts relative to the ligand as in the cobalt complex. Cobalt and copper possess a stabilizing effect on the electronic system of ligand. The same effect was observed on the basis of the IR spectra of Co and Cu complexes of dacarbazine (an increase in the values of the wavenumbers and intensities of the aromatic ring). Our long-term study on the influence of metals on the electronic system of ligand (aromatic and heteroaromatic system of five- and six-membered carboxylic acids) showed that 3d transition metals (including cobalt and copper) stabilize the electronic system of ligands [20,21,22,23,24]. In the spectra of such complexes (a) an increase in the values of signals of the aromatic protons (in 1HNMR spectra), and (b) an increase in wavenumbers and intensity of bands assigned to the aromatic system (in the IR spectra) compared with the appropriate signals/bands in the spectra of ligand were observed.

3.3. Structure, Aromaticity, and NBO Analysis of Dacarbazine and Copper and Cobalt Complexes

The dacarbazine molecule has a triazene moiety linked to an imidazole ring. IR, NMR studies and literature X-ray data on the structure of dacarbazine complexes [46] showed that metals such as copper and cobalt coordinate the ligand (dacarbazine) through the carbonyl group and the nitrogen atom from the imidazole ring. A similar type of metal-ligand coordination has been observed for 3d-transition metal complexes with imidazole carboxylic acid [38]. On the basis of spectroscopic studies (IR, Raman) and theoretical calculations, it was found that 3d metals affect the electronic charge distribution of the aromatic ring of ligand. With the increase in the electronic stability of the ligand, the reactivity of the compound changes, which is important in the case of biologically active ligands. Theoretical structures of dacarbazine complexes with cobalt and copper were calculated. The calculated bond lengths and the values of the selected angles were compared with these ones obtained for experimentally solved structures (X-ray) available in the literature. The geometric indices of aromaticity (HOMA, GEO, EN, and I5) and the electronic charge distribution by the NBO method were also calculated. Figure 3 shows the experimentally determined structure [46] and the calculated two dacarbazine conformers (at B3LYP method). Figure 4 shows the experimental structure of the complex of copper with dacarbazine (A) [47] and the theoretically modeled structures of the complexes with copper (B) and cobalt (C).
Both theoretically modeled dacarbazine conformational structures have similar energy values. According to the calculated geometric aromaticity indices, the imidazole ring of the DAC1 structure is more stable than the DAC2 structure (Table 3). The system is stabilized by an intramolecular hydrogen bond between the hydrogen atom of the amino group and the nitrogen atom of the triazene group. The length of this bond is 2.154 Å. In the DAC2 structure, the amide moiety is rotated 180° along the C1–C2 carbons as compared to the DAC1 structure. The aromatic system of imidazole in this structure is less stable, as evidenced by the values of the aromaticity indices. The structure of the DAC1 monomer is the same as that of experimentally determined dacarbazine. In the structure determined by the X-ray diffraction method [42], there is an intramolecular hydrogen bond that stabilizes the aromatic system between the hydrogen atom of the amino group and the triazene nitrogen atom with a length of 2.294 Å, and an intermolecular hydrogen bond with a length of 1.945 Å. The values of geometric aromaticity indices calculated for the experimental structure and the theoretically modeled DAC1 are at a similar level (Table 3). The differences in the lengths of the corresponding bonds in the experimentally determined dacarbazine structure and the theoretically modeled DAC1 monomer range from 0.001 to 0.026 Å.
The complexation of dacarbazine with 3d transition metals (copper and cobalt) influences the stability of the ligand’s electronic system, which is observed by changing the aromaticity of the system expressed in the values of aromaticity indices, changes in the distribution of electronic charge, and changes in the bond lengths in the studied systems. The comparison of the structure of the dacarbazine complex with copper [47] to dacarbazine [46] shows that the aromatic system of the imidazole ring is stabilized by complexation (Table 3). It was observed both when comparing the ligand structure to the complex for the experimental and theoretically modeled structures (the effect of the metal is more noticeable when comparing the experimental structures than the theoretically modeled ones). Under the influence of copper, there is a significant change in the length of the imidazole ring bonds, which results in an increase in the HOMA and I5 indices. Complexation of metals by ligand slightly changes the bond length of the triazene moiety (Table 3). The complexation of cobalt with dacarbazine influences the change of aromaticity of the imidazole ligand ring to a lesser extent than complexation of copper (when comparing theoretical structures). A slight increase in the value of the I5 index and a slight decrease in the value of the HOMA index were observed. The theoretical structures of copper and cobalt complexes have similar bond lengths and angles.
It was observed that the complexation of dacarbazine with copper and cobalt slightly changed the bond length of the N4-N5-N6 triazene moiety and the length of the N6-C5 and N6-C6 bonds (nitrogen N6 with methyl groups). The analysis of the NBO electronic charge distribution (Table 4) showed that the values of the electronic charges on the nitrogen atoms of the triazene group and the carbons of the methyl groups are similar in the complexes and the ligand.
On the other hand, the influence of metals on the distribution of electronic charge on the imidazole ring of dacarbazine was observed. The electronic density around the carbon atom labeled C2 and the nitrogen atom N3 increased in copper and cobalt complexes compared to the ligand, while around the C3, C4, and N2 atoms the electronic density decreased (Table 4). The effect of copper on the electronic charge distribution in the imidazole ring of the ligand is greater than that of cobalt. The total electronic density of the C-ring of the copper complex is greater than that of the cobalt complex. Studies of 4-imidazole carboxylic acid complexes with 3d-transition metals (including copper and cobalt) also showed that copper and cobalt had a stabilizing effect on the aromatic system of the imidazole ring [48].

3.4. HOMO and LUMO Orbitals

The energy values of the HOMO and LUMO orbitals for dacarbazine as well as copper and cobalt complexes were calculated. The shapes of the orbitals are shown in Figure 5 and Figure 6, while in Table 5 the values of the energy of the orbitals and other electronic parameters are gathered. The highest-occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) play an important role in predicting charge transfer in the molecule, chemical reactivity, bioactivity, and compound stability [46]. The shapes of the LUMO and HOMO molecular orbitals of the dacarbazine complexes with copper and cobalt are shown in Figure 6. The energy values of the HOMO and LUMO orbitals and the value of the energy difference between the HOMO-LUMO levels for the studied compounds are presented in Table 5. The value of the HOMO-LUMO energy difference (Energy gap) for the copper complex is 0.3761 eV, and for the cobalt complex it is 1.5894 eV and the values are lower than those obtained for dacarbazine (3.057 eV for DAC1). This indicates a decrease in kinetic stability and an increase in the reactivity of the complexes toward the ligand. Other general reactivity descriptors [49,50] such as ionization potential (I), electron affinity (A), electronegativity (χ), chemical hardness (η), softness (e), and electrophilicity index (ω) calculated on the basis of the energy of HOMO and LUMO are summarized in Table 5. The data show that the chemical hardness of the complexes is lower than that of the ligand. The electrophilicity index (ω) provides information not only about the reactivity but also about the toxicity of the molecule. This indicator is related to the stabilization of energy when the system acquires an additional electrostatic charge from the environment and quantifies the global electrophilic force of the molecule [51]. The electrophilicity index for the copper complex is much higher than that of the free ligand (four times higher), while for the cobalt complex it is twice as high as obtained for the ligand.

3.5. Docking Studies

Comparing the average number of hydrogen bonds formed with the B-DNA structure, it was observed that the number of H-Bonds does not change significantly in the case of decarbazine metal complexes compared to ligands.
The obtained results indicate that the significant decrease (~4 kcal/M) in the free energy of binding the complexes with B-DNA is not due to formed hydrogen bonds (Figure 7). We conclude that formation of the complex causes favorable changes to occur in the distribution of the electronic charge. These changes in the electron density distribution affect the affinity to biological macromolecules like DNA. The tested compounds had an entropy value of S = 13.75 kcal/mol and an internal energy with values Dac1 (−4.37 kcal/mol), Dac2 (−4.26 kcal/mol), Cu-DAC (−8.13 kcal/mol), and Co-DAC (−8.19 kcal/mol), respectively (Table 6). The most energetically advantageous complexes with DNA formed two hydrogen bonds with thymine (Figure 8). In the case of constant inhibition, the value for metal complexes decreased by more than 1000-fold, and the free-binding energy decreased by about 4 kcal/mol.
In Ahmad and Ahmad’s studies, the binding constant (Kb) was 7.89 × 104 experimental isothermal titration calorimetry of ct-DNA-DAC, and Kb = 6.99 × 104 in docking studies for B-DNA-DAC [52]. Dacarbazine–ctDNA binding energy obtained from ITC experiments was –5.49 kcal/mol and it is comparable with minimized free energy of DAC–DNA obtained in docking studies (−5.35 kcal/mol). The most favorable interactions were found in the A-T residue of the minor DNA groove and hydrogen bonding and van der Waals interactions likely play a major role in these interactions. Previous research described by Wang et al. [53] obtained by atomic force microscopy and nanostructural image of DAC–DNA complexes indicate interaction at some specific site of the DNA sequence without intercalation. Additionally, DAC shows poor binding affinity to ssDNA [53]. Based on cyclic voltammetry, differential pulse voltammetry studies, and using spectroscopic methods performed by Radi et al., it has been concluded that decarbazine binds to double-stranded DNA by the combined effect of intercalation and electrostatic interactions [54].
Temerk and Ibrahim [19] describe that both DAC and DAC-Cu(II) complexes intercalate with base stacking of dsDNA, independent of ionic strength. In contrast, the interaction of a DAC with ssDNA, which is negatively charged, is electrostatic attractions. However, in the presence of Cu(II), these interactions are much stronger. Values of standard Gibbs free energy (−ΔG◦) in 278 K for DAC-ssDNA is 28.74 kJ/M (6.87 kcal/M) and for DAC-Cu(II)-ssDNA is 29.98 kJ/M (7.17 kcal/M) with enthalpy ΔH◦ = 21.98 kJ/M and 22.29 kJ/M, respectively. For comparison −ΔG◦ for DAC-dsDNA is 25.75 kJ/M and 27.55 kJ/mol for DAC-Cu-dsDNA [19].
The presented results indicate the formation of beneficial interactions and changes in the distribution of the electronic charge in decarbazine complexes with copper and cobalt, which in turn may enhance the interaction with macromolecules such as DNA.

3.6. Thermogravimetric Study

Thermogravimetric curves of dacarbazine and its metal complexes with copper(II) and cobalt(II) are shown in Figure 9. Dacarbazine underwent immediate decomposition above 205 °C. Rapid mass loss observed on the TG curve was accompanied by a narrow signal on the DTG curve. In the first stage of thermal decomposition of dacarbazine (in the range of 205–210 °C), a 50% loss in weight occurred. The mid weight loss was due to degradation of the triazene group and thermal dissociation of the amide group. The second stage of thermal decomposition of decarbazine was less rapid (above 201 °C)—the imidazole ring and C=O group totally decomposed, and the total mass loss was at 675 °C (Table 7).
In the first stage of thermal decomposition of copper(II) and cobalt(II) complexes of dacarbazine the process of methanol desolvation occurred. The cobalt complex totally lost methanol (1.75 mole of CH3OH) in the range of 60–120 °C, and the copper complex (1.5 mole of CH3OH) in the range of 60–110 °C. The thermal decomposition of cobalt compound started at 175 °C (at a much lower temperature than decarbazine). The first stage of decomposition was as rapid as it was in the case of ligand (with a 33% loss in weight). The second stage of decomposition ran gently until achievement of 540 °C, and then the rapid mass loss occurred. The TG curve possessed a similar shape as the TG curve of decarbazine. At 630 °C the cobalt complex totally decomposed to CoO.
At the beginning of the thermal decomposition of copper(II) complex of dacarbazine the methanol desolvation occurred. Next the thermal decomposition of the complex started at about 160 °C (a slightly lower temperature than in the case of dacarbazine). During the thermal decomposition of copper complex, no stable intermediate products were observed (similarly as in the case of the cobalt complex). Above 580 °C the final product of decomposition CuO occurred. Comparing the curves of thermal decomposition of studied compounds it appeared that complexes of dacarbazine were less stable than ligand alone. The thermal decomposition of cobalt and copper complexes started at temperatures 35 and 50 degrees lower than ligand, respectively. The cobalt complex of dacarbazine was more thermal stable than the copper complex.

4. Conclusions

The spectroscopic IR and NMR studies revealed that dacarbazine coordinates metal cations (cobalt(II) and copper(II)) through the oxygen atom of carbonyl group and nitrogen atoms from the imidazole ring. It was stated on the basis of significant shifts of bands assigned to the vibrations of the aforementioned groups from the infrared spectra and chemical shifts of protons from the 1HNMR spectra (changes in the electronic charge density under the influence of complexation). The thermal analysis showed that cobalt(II) and copper(II) complexes of dacarbazine were less stable than ligand. The cobalt complex of dacarbazine was more thermal stable than the copper complex.
Based on the literature data (X-ray diffraction) and spectroscopic data, the theoretical structures of the complexes were modeled. Experimentally determined and theoretically calculated structures were compared. For the modeled structures, the distribution of electronic charge, aromaticity, and energies of the HOMO and LUMO orbitals were calculated. This allowed the assessment of the effect of complexation with copper and cobalt on the stabilization of the aromatic system and the reactivity of the investigated complexes. Complexation with copper and cobalt changes the electronic charge distribution of dacarbazine. The aromatic system of dacarbazine is stabilized (aromaticity indices increase) due to metal complexation. The values of the calculated energy of the HOMO and LUMO orbitals indicate an increase in the reactivity of the dacarbazine molecule after complexation. This has a significant impact on the biological activity of these substances. The electron charge distribution of a ligand provides information about the stability of a given compound as well as its reactivity and biological activity. By complexing ligands of high biological importance with metals, we influence the electronic stability of a given ligand and, at the same time, we can influence its biological activity. Additionally, by examining the thermal properties of given chemical compounds, we can potentially check their susceptibility to thermal utilization and test the decomposition products.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ma14123274/s1, Figure S1. 1HNMR spectra for dacarbazine and dacarbazine complex with cobalt (II).

Author Contributions

Conceptualization, G.Ś., M.K., and R.Ł., methodology, G.Ś. and M.K., validation, W.L. and R.Ł.; formal analysis, A.O.A., M.M., M.K., G.Ś., and M.S., investigation, W.L. and R.Ł.; writing—original draft preparation, G.Ś., M.K., R.Ś., A.G., and M.M., visualization, G.Ś., A.G., and A.O.A.; supervision, G.Ś.; project administration, W.L. All authors have read and agreed to the published version of the manuscript.

Funding

Studies have been carried out in the framework of the work no. WZ/WB-IIS/5/2020 in Bialystok University of Technology and financed from the funds for Science, Ministry of Science and Higher Education of Poland.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

No potential conflict of interest are reported by the authors.

References

  1. Serrone, L.; Zeuli, M.; Sega, F.M.; Cognetti, F. Dacarbazine-based chemotherapy for metastatic melanoma: Thirty-year experience overview. J. Exp. Clin. Cancer Res. 2000, 19, 21–34. [Google Scholar]
  2. Reid, J.M.; Kuffel, M.J.; Miller, J.K.; Rios, R.; Ammes, M.M. Metabolic activation of dacarbazine by human cytochromes p450: The role of CYP1A1, CYP1A2e CYP2E1. Clin. Cancer Res. 1999, 5, 2192–2197. [Google Scholar]
  3. Nussbaumera, S.; Bonnabry, P.; Veuthey, J.L.; Fleury-Souverain, S. Analysis of anticancer drugs: A review. Talanta 2011, 85, 2265–2289. [Google Scholar] [CrossRef] [PubMed]
  4. Moody, C.L.; Wheelhouse, R.T. The Medicinal Chemistry of Imidazotetrazine Prodrugs. Pharmaceuticals 2014, 7, 797–838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Bonmassar, L.; Marchesi, F.; Pascale, E.; Franzese, O.; Margison, G.P.; Bianchi, A.; D’Atri, S.; Bernardini, S.; Lattuada, D.; Bonmassar, E.; et al. Triazene compounds in the treatment of acute myeloid leukemia: A short review and a case report. Curr. Med. Chem. 2013, 20, 2389–2401. [Google Scholar] [CrossRef] [PubMed]
  6. Marchesi, F.; Turriziani, M.; Tortorelli, G.; Avvisati, G.; Torino, F.; Vecchis, L.D. Triazene compounds: Mechanism of action and related DNA repair systems. Pharmacol. Res. 2007, 56, 275–287. [Google Scholar] [CrossRef]
  7. Pourahmad, J.; Amirmostofian, M.; Kobarfard, F.; Shahraki, J. Biological reactive intermediates that mediate dacarbazine cytotoxicity. Cancer Chemother. Pharmacol. 2009, 65, 89–96. [Google Scholar] [CrossRef]
  8. Hayward, I.P.; Parson, P.G. Epigenetic effects of the methylating agent 5-(3-dimethyl-1-triazeno)-imidazole-4-carboxamide in human melanoma cells. Aust. J. Exp. Biol. Med. Sci. 1984, 62, 597–606. [Google Scholar] [CrossRef]
  9. Al-Qatati, A.; Aliwaini, S. Combined pitavastatin and dacarbazine treatment activates apoptosis and autophagy resulting in synergistic cytotoxicity in melanoma cells. Oncol. Lett. 2017, 14, 7993–7999. [Google Scholar] [CrossRef] [Green Version]
  10. Naserian, M.; Ramazani, E.; Iranshahi, M.; Tayarani-Najaran, Z. The Role of SAPK/JNK pathway in the synergistic effects of metformin and dacarbazine on apoptosis in Raji and Ramos lymphoma cells. Curr. Mol. Pharmacol. 2018, 11, 336–342. [Google Scholar] [CrossRef]
  11. Finotello, R.; Stefanello, D.; Zini, E.; Marconato, L. Comparison of doxorubicin–cyclophosphamide with doxorubicin–dacarbazine for the adjuvant treatment of canine hemangiosarcoma. Vet. Comp. Oncol. 2015, 15, 25–35. [Google Scholar] [CrossRef] [Green Version]
  12. Song, M.; Zhang, R.; Wang, X. Nano-titanium dioxide enhanced biosensing of the interaction of dacarbazine with DNA and DNA bases. Mater. Lett. 2006, 60, 2143–2147. [Google Scholar] [CrossRef]
  13. Shen, Q.; Wang, X.; Fu, D. The amplification effect of functionalized gold nanopar-ticles on the binding of anticancer drug dacarbazine to DNA and DNA bases. Appl. Surf. Sci. 2008, 255, 577–580. [Google Scholar] [CrossRef]
  14. Matejczyk, M.; Świsłocka, R.; Golonko, A.; Lewandowski, W.; Hawrylik, E. Cytotoxic, genotoxic and antimicrobial activity of caffeic and rosmarinic acids and their lithium, sodium and potassium salts as potential anticancer compounds. Adv. Med. Sci. 2008, 63, 14–21. [Google Scholar] [CrossRef]
  15. Jabłońska-Trypuć, A.; Świderski, G.; Krętowski, R.; Lewandowski, W. Newly synthesized doxorubicin complexes with selected metals—Synthesis, structure and anti-breast cancer activity. Molecules 2017, 22, 1106. [Google Scholar] [CrossRef]
  16. Trynda-Lemiesz, L.; Śliwińska-Hill, U. Kompleksy metali w terapii nowotworowej. Teraźniejszość i przyszłość. J. Oncol. 2011, 61, 465–474. [Google Scholar]
  17. Kumari, T.; Shukla, J.; Joshin, S. Study on the complex formation and anticancer effect of complex, zinc(II)-dacarbazine. Int. J. Chem. Sci. 2011, 9, 1751–1762. [Google Scholar]
  18. Shukla, J.; Pitre, K.S. Role of bio-metal Fe(III) in anticancer effect of dacarbazine. Indian J. Physiol. Pharmacol. 1998, 42, 223–230. [Google Scholar]
  19. Temerk, Y.; Ibrahim, H.I. Binding mode and thermodynamic studies on the interaction of the anticancer drug dacarbazine and dacarbazine–Cu(II) complex with single and double stranded DNA. J. Pharm. Biomed. Anal. 2014, 95, 26–33. [Google Scholar] [CrossRef]
  20. Lewandowski, W.; Kalinowska, M.; Lewandowska, H. The influence of metals on the electronic system of biologically important ligands, Spectroscopic study of benzoates, salicylates, nicotinates and isoorates. Review. J. Inorg. Biochem. 2005, 99, 1407–1423. [Google Scholar] [CrossRef]
  21. Koczoń, P.; Hrynaszkiewicz, T.; Świsłocka, R.; Samsonowicz, M.; Lewandowski, W. Spectroscopic (Raman, FT-IR and NMR) study of alkaline metal nicotinates and isonicotinates. Vib. Spectrosc. 2003, 33, 215–222. [Google Scholar] [CrossRef]
  22. Lewandowska, M.; Janowski, A.; Lewandowski, W. Spectroscopic Investigations on Lanthanide Complexes with Salicylic Acid. Can. J. Spectr. 1984, 29, 87–92. [Google Scholar]
  23. Lewandowski, W.; Kalinowska, M.; Lewandowska, H. The influence of halogens on the electronic system of biologically important ligands. Spectroscopic study of halogenobenzoic acids, halogenobenzoates and 5-halogenouracils. Review. Inorg. Chim. Acta 2005, 358, 2155–2166. [Google Scholar] [CrossRef]
  24. Koczoń, P.; Piekut, J.; Borawska, M.; Lewandowski, W. Vibrational structure and antimicrobial activity of selected isonicotinates, potassium picolinate and nicotinate. J. Mol. Struct. 2003, 651–653, 651–656. [Google Scholar] [CrossRef]
  25. Kalinowska, M.; Borawska, M.; Świsłocka, R.; Piekut, J.; Lewandowski, W. Spectroscopic (IR, Raman, UV, 1H and 13C NMR) and microbiological studies of Fe(III), Ni(II), Cu(II), Zn(II) and Ag(I) picolinates. J. Mol. Struct. 2007, 834–836, 419–425. [Google Scholar] [CrossRef]
  26. Świderski, G.; Świsłocka, R.; Łyszczek, R.; Wojtulewski, S.; Samsonowicz, M.; Lewandowski, W. Thermal, spectroscopic, X-ray and theoretical studies of metal complexes (sodium, manganese, copper, nickel, cobalt and zinc) with pyrimidine-5-carboxylic and pyrimidine-2-carboxylic acids. J. Therm. Anal. Calorim. 2019, 138, 2813–2837. [Google Scholar] [CrossRef] [Green Version]
  27. Świderski, G.; Wojtulewski, S.; Kalinowska, M.; Świsłocka, R.; Wilczewska, A.Z.; Pietryczuk, A.; Cudowski, A.; Lewandowski, W. The influence of selected transition metal ions on the structure, thermal and microbiological properties of pyrazine-2-carboxylic acid. Polyhedron 2020, 175, 114173. [Google Scholar] [CrossRef]
  28. Świderski, G.; Łyszczek, R.; Wojtulewski, S.; Kalinowska, M.; Świsłocka, R.; Lewandowski, W. Comparison of structural, spectroscopic, theoretical and thermal properties of metal complexes (Zn, Mn (II), Cu (II), Ni (II) and Co (II)) of pyridazine-3-carboxylic acid and pyridazine-4-carboxylic acids. Inorg. Chim. Acta 2020, 512, 119865. [Google Scholar] [CrossRef]
  29. Świderski, G.; Kalinowska, M.; Wilczewska, A.Z.; Malejko, J.; Lewandowski, W. Lanthanide complexes withpyridinecarboxylic acids–Spectroscopic and thermal studies. Polyhedron 2018, 150, 97–109. [Google Scholar] [CrossRef]
  30. Świderski, G.; Kalinowska, M.; Rusinek, I.; Samsonowicz, M.; Rzączyńska, Z.; Lewandowski, W. Spectroscopic (IR, Raman) and thermogravimetric studies of 3d-metal cinchomeronates and dinicotinates. J. Therm. Anal. Calorim. 2016, 126, 1521–1532. [Google Scholar] [CrossRef] [Green Version]
  31. Świderski, G.; Lewandowska, H.; Świsłocka, R.; Wojtulewski, S.; Siergiejczyk, L.; Wilczewska, A.Z.; Misztalewska, I. Spectroscopic (IR, Raman, NMR), thermal and theoretical (DFT) study of alkali metal dipicolinates (2,6) and quinolinates (2,3). Arab. J. Chem. 2019, 12, 4414–4426. [Google Scholar] [CrossRef] [Green Version]
  32. Padnya, P.; Shibaeva, K.; Arsenyev, M.; Baryshnikova, S.; Terenteva, O.; Shiabiev, I.; Khannanov, A.; Boldyrev, A.; Gerasimov, A.; Grishaev, D.; et al. Catechol-Containing Schiff Bases on Thiacalixarene: Synthesis, Copper (II) Recognition, and Formation of Organic-Inorganic Copper-Based Materials. Molecules 2021, 26, 2334. [Google Scholar] [CrossRef]
  33. Colombo, A.; Dragonetti, C.; Roberto, D.; Fagnani, F. Copper complexes as alternative redox mediators in dye-sensitized solar cells. Molecules 2021, 26, 194. [Google Scholar] [CrossRef]
  34. Pessoa, J.C.; Correia, I. Misinterpretations in Evaluating Interactions of Vanadium Complexes with Proteins and Other Biological Targets. Inorganics 2021, 9, 17. [Google Scholar] [CrossRef]
  35. Soldatović, T. Mechanism of Interactions of Zinc(II) and Copper(II) Complexes with Small Biomolecules. In Basic Concepts Viewed from Frontier in Inorganic Coordination Chemistry; IntechOpen: London, UK, 2018. [Google Scholar]
  36. De Souza, I.C.A.; De Souza Santana, S.; Gómez, J.G.; Guedes, G.P.; Madureira, J.; De Ornelas Quintal, S.M.; Lanznaster, M. Investigation of cobalt(iii)-phenylalanine complexes for hypoxia-activated drug delivery. Dalton Trans. 2020, 49, 16425–16439. [Google Scholar] [CrossRef]
  37. Gaussian; Version 9, Revision A.02; Gaussian Inc.: Wallingford, CT, USA, 2016.
  38. Krygowski, T.M.; Cyrański, M. Separation of the energetic and geometric contributions to the aromaticity. Part IV. A general model for thep-electron systems. Tetrahedron 1996, 52, 10255–10264. [Google Scholar] [CrossRef]
  39. Bird, C. A new aromaticity index and its application to fivemembered ring heterocycles. Tetrahedron 1985, 41, 1409–1414. [Google Scholar] [CrossRef]
  40. Weinhold, F.; Landis, C.R. Natural bond orbitals and extensions of localized bonding concepts. Chem. Educ. Res. Pract. 2001, 2, 91–104. [Google Scholar] [CrossRef]
  41. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  42. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. Software news and updates AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef] [Green Version]
  43. Dassault Systèmes. Biovia Discovery Studio Modeling Environment; Dassault Systèmes Biovia: San Diego, CA, USA, 2016; Available online: https://www.3ds.com/products-services/biovia/products/molecular-modeling-simulation/biovia-discovery-studio/ (accessed on 13 June 2021).
  44. Pettersen, E.F.; Goodard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Gunasekaran, S.; Kumaresan, S.; Arunbalaji, R.; Anand, G.; Srinivasan, S. Density functional theory study of vibrational spectra, and assignment of fundamental modes of dacarbazine. J. Chem. Sci. 2008, 120, 315–324. [Google Scholar] [CrossRef]
  46. Freeman, H.C.; Hutchinson, N.D. The crystal structure of the anti-tumor agent 5-(3,3-dimethyl-1-triazenyl) imidazole-4-carboxamide. Acta Crystall. Sec. B 1979, 35, 2051–2054. [Google Scholar] [CrossRef]
  47. Freeman, H.C.; Hutchinson, N.D. The crystal structures of two copper(II) complexes of the antitumor agent 5-(3,3-dimethyl-1-triazenyl) imidazole-4-carboxamide. Acta Crystall. Sec. B 1979, 35, 2045–2050. [Google Scholar] [CrossRef]
  48. Świderski, G.; Wilczewska, A.Z.; Świsłocka, R.; Kalinowska, M.; Lewandowski, W. Spectroscopic (IR, Raman, UV–Vis) study and thermal analysis of 3d-metal complexes with 4-imidazolecarboxylic acid. J. Therm. Anal. Calorim. 2018, 134, 513–525. [Google Scholar] [CrossRef] [Green Version]
  49. Fukui, K. Role of frontier orbitals in chemical reactions. Science 1982, 218, 747–754. [Google Scholar] [CrossRef] [Green Version]
  50. Parr, R.G.; Szentpály, L.V.; Liu, S. Electrophilicity Index. J. Am. Chem. Soc. 1999, 121, 1922–1924. [Google Scholar] [CrossRef]
  51. Samsonowicz, M. Molecular structure of phenyl- and phenoxyacetic acids—Spectroscopic and theoretical study. Spectrochim. Acta A 2014, 118, 1386–1425. [Google Scholar] [CrossRef]
  52. Ahmad, I.; Ahmad, M. Dacarbazine as a minor groove binder of DNA: Spectroscopic, biophysical and molecular docking studies. Int. J. Biol. Macromol. 2015, 79, 193–200. [Google Scholar] [CrossRef]
  53. Wang, X.; Li, Y.; Gong, S.; Fu, D. A spectroscopic study on the DNA binding behavior of the anticancer drug dacarbazine. Spectrosc. Lett. 2002, 35, 751–756. [Google Scholar] [CrossRef]
  54. Radi, A.E.; Eissa, A.; Nassef, H.M. Voltammetric and spectroscopic studies on the binding of the antitumor drug dacarbazine with DNA. J. Electroanal. Chem. 2014, 717–718, 24–28. [Google Scholar] [CrossRef]
Materials 14 03274 g001 550
Figure 1. Mechanism of action of dacarbazine.
Figure 1. Mechanism of action of dacarbazine.
Materials 14 03274 g001
Materials 14 03274 g002 550
Figure 2. IR spectra of dacarbazine (blue line), cobalt complex (red line), copper complex (green line), and Raman spectra of dacarbazine (black line).
Figure 2. IR spectra of dacarbazine (blue line), cobalt complex (red line), copper complex (green line), and Raman spectra of dacarbazine (black line).
Materials 14 03274 g002
Materials 14 03274 g003 550
Figure 3. Crystal and theoretical structure (B3LYP/6-311++g(d,p)) of dacarbazine.
Figure 3. Crystal and theoretical structure (B3LYP/6-311++g(d,p)) of dacarbazine.
Materials 14 03274 g003
Materials 14 03274 g004 550
Figure 4. Crystal (A), theoretical structure (B3LYP/6-311++g(d,p)) of cobalt (B) and copper (C) complex with dacarbazine.
Figure 4. Crystal (A), theoretical structure (B3LYP/6-311++g(d,p)) of cobalt (B) and copper (C) complex with dacarbazine.
Materials 14 03274 g004
Materials 14 03274 g005 550
Figure 5. HOMO and LUMO orbitals of dacarbazine.
Figure 5. HOMO and LUMO orbitals of dacarbazine.
Materials 14 03274 g005
Materials 14 03274 g006 550
Figure 6. HOMO and LUMO orbitals of copper and cobalt complexes with dacarbazine.
Figure 6. HOMO and LUMO orbitals of copper and cobalt complexes with dacarbazine.
Materials 14 03274 g006
Materials 14 03274 g007 550
Figure 7. (Upper) Representation of the number of conformations formed by both ligand conformers (DAC1, DAC2) and metal complexes Co(II)-DAC and Cu(II)-DAC. (Lower) Average number of hydrogen bonds formed in complexes and ligands and free binding energy.
Figure 7. (Upper) Representation of the number of conformations formed by both ligand conformers (DAC1, DAC2) and metal complexes Co(II)-DAC and Cu(II)-DAC. (Lower) Average number of hydrogen bonds formed in complexes and ligands and free binding energy.
Materials 14 03274 g007
Materials 14 03274 g008a 550Materials 14 03274 g008b 550Materials 14 03274 g008c 550
Figure 8. Best docking poses of decarbazine and complexes with B-DNA structure (Upper) in 2D with hydrogen bonds (green), Van der Waals (light green), and Pi-Alkyl (pink) interactions. (Lower) 3D visualization of location in the minor groove of DNA.
Figure 8. Best docking poses of decarbazine and complexes with B-DNA structure (Upper) in 2D with hydrogen bonds (green), Van der Waals (light green), and Pi-Alkyl (pink) interactions. (Lower) 3D visualization of location in the minor groove of DNA.
Materials 14 03274 g008aMaterials 14 03274 g008bMaterials 14 03274 g008c
Materials 14 03274 g009 550
Figure 9. TG/DTG curves of dacarbazine (A), cobalt (B), and copper (C) complexes.
Figure 9. TG/DTG curves of dacarbazine (A), cobalt (B), and copper (C) complexes.
Materials 14 03274 g009
Table
Table 1. The wavenumbers from the IR (theoretical and experimental) and Raman spectra of dacarbazine and IR spectra of dacarbazine complexes with cobalt and copper.
Table 1. The wavenumbers from the IR (theoretical and experimental) and Raman spectra of dacarbazine and IR spectra of dacarbazine complexes with cobalt and copper.
Calculated IR (B3LYP 6-311G(d,p)Experimental Spectra (cm−1)
DAC 1 DAC 2 DAC Cobalt ComplexCopper Complex
Unscaled Scaled = 0.967 IntensityUnscaledScaled = 0.967IntensityIR RamanIRIRAssignments
37033581105.503728360572.843383 s3371 w3406 s3406 sνasNH2
3647352758.233560344365.93--3322 s-νNH
3549343260.673585346753.443269 m---νsNH2
323831312.26324031332.093174 m-3188 s3186 sνCHar
314430407.763143303911.283147 m3141 vw--νsCH3
313930358.76313930355.322964 m---ρsCH
3085298313.513097299515.072905 m 2925 m2923 mνasCH3
3078297626.133074297328.842793 m--2796 wνasCH3
3014291547.623021292140.932753 m---νsCH3
3003290478.482997289894.762612 m---νsCH3
17491691494.0717401683326.181609 vs1604 w1605 s1609 sνC=O
------1656 vs-1637 vs1639 vsνCH3
16201567141.6916071554256.98----ρC-NH2
------1561 w1559 w1561 m1566 mρsCH3
156915173.591586153455.591510 m 1525 s1528 sνCN
15271477197.9015251475155.811476 s1488 vw1487 vs1488 vsρC-NH2
1483143427.901481143225.741436 m1449 vs1440 m1440 mρNH2
1451140375.221453140532.351402 s-1405 m1405 mνCC, νring, ρNH2
143813915.80143813919.68----νCC, αN-CH2, βNCH
1431138447.981430138374.131381 s1388 w--βasCH3, νring
1389134321.231388134224.581344 s-1352 vs1351 vsβsCH3, νring
13631318114.3913671322408.841304 m-1306 s1307 sβsCH3, νring
1331128736.741328128440.871270 w1289 vw--βsCH3, νring, βNCH
127012282.301260121816.531259 w1258 vw1253 vs1255 mνCN, βCH
1231119067.151233119220.711231 m1204 vw--νring, βCH
1153111525.501152111439.121183 w---βCH + νNN
1136109926.51113710999.081072 s1091 vw1091 s1091 sνNN, βCH
1114107724.06111010737.62----βCH + βCNC
10911055202.221093105717.911049 w1012 m--ρCH2, βCH
1084104812.7210751040117.70986 w---ρC-NH2, βCNC
1062102714.221061102618.40962 w-969 m980 wρCH2, βCNC
9519204.3295592315.06904 w913 vw902 w902 wρCH3, βCH
91588514.8691288216.25896 w---ρCH3, νring
8177909.748127855.02882 w---Δring, νring(CN)
7997732.217997732.40868 w873 vw--ρCH3, νring
79777110.9579877214.79796 w---βNCH, βNCC
67265013.4468366014.27690 w691 w--βNH, βCCO
6686460.826796570.01653 vw---β(CCN + NH)
63261125.6260258243.51630 m-646 m649 mαNNN, αCNC
5865670.505965760.65542 m-565 m568 mρNH2, αNNH, αCCO
55253473.21574555104.33450 w445 vw--ringdef
55053239.8256854910.64----ringdef
44743210.0144943411.86416 vw---αC-NH2
* Abbreviation: band intensity: s—strong, m—medium, w—weak, vw—very weak, type of vibrations: ν—stretching, ρ—bending out-of-plane β—bending in-plane, α—deformation in-plane, Δ—deformation out-of-plane, ringdef—deformation of ring, s—symmetric, as—asymmetric.
Table
Table 2. The chemical shifts (ppm) from the 1HNMR theoretical (GIAO/B3LYP/6-311++G**) and experimental spectra of dacarbazine and experimental spectra of cobalt and copper complexes with dacarbazine.
Table 2. The chemical shifts (ppm) from the 1HNMR theoretical (GIAO/B3LYP/6-311++G**) and experimental spectra of dacarbazine and experimental spectra of cobalt and copper complexes with dacarbazine.
DAC 1DAC 2Dacarbazine Complexes
CobaltCopper
H (amide group)Exp.7.29, 7.41-6.20, 5.40-
Theoret.6.86, 4.576.87, 4.406.388.29, 5.42
H (CHring)Exp.7.83-8.04-
Theoret.7.206.988.907.74
H (NHring)Exp.11.52---
Theoret.8.568.729.438.81
H (CH3-triazene group)Exp.3.13, 3.50-3.49-
Theoret.3.67, 3.48, 3.36, 3.19, 2.29, 2.594.73, 4.37, 3.91, 3.88, 3.69, 3.173.78, 3.74, 3.71, 3.44, 3.34, 3.143.62, 3.61, 3.61, 3.24, 3.23, 3.03
Table
Table 3. Geometrical aromaticity parameters and selected bond lengths and angle values obtained for dacarbazine and copper(II) and cobalt(II) complexes of dacarbazine.
Table 3. Geometrical aromaticity parameters and selected bond lengths and angle values obtained for dacarbazine and copper(II) and cobalt(II) complexes of dacarbazine.
DacarbazineCopper ComplexCobalt Complex
Exp [46]DAC1DAC2Exp [47]CalcCalc
Energy-−638.54−638.54-−3199.60−2941.88
Dipole m-10.74985.5647-21.841122.5754
Aromaticity indices
HOMA0.9050.8930.8800.9250.8970.890
GEO0.0600.0690.0790.0570.0650.067
EN0.0350.0390.0410.0170.0390.043
I567.9767.3965.3570.4468.4267.68
Bond lengths [A]
C1-O11.2301.2181.2221.3621.2361.247
C1-N11.3381.3701.3661.3091.3451.342
C1-C21.4701.4921.4841.4601.4771.465
C2-C31.3791.3921.3881.3871.3931.392
C3-N21.3751.3771.3811.3701.3831.381
N2-C41.3521.3701.3671.3451.3561.359
C4-N31.3331.3071.3091.3111.3111.312
N3-C21.3871.3761.3831.3791.3791.384
C3-N41.3831.3841.3801.3871.3821.380
N4-N51.2851.2701.2791.2861.2731.274
N5-N61.3041.3251.3261.2991.3141.313
N6-C51.4511.4561.4581.4501.4551.455
N6-C61.4491.4531.4511.4501.4581.458
Angles (°)
N2-C3-N4127.50125.36125.07116.37115.48116.15
C3-N4-N5112.80113.71113.23112.17115.88115.46
N4-N5-N6113.42115.20 114.26114.66114.71114.64
N5-N6-C6116.55115.90116.01122.46121.71121.70
C5-N6-C6120.96120.43120.46120.31120.82120.76
C1-C2-N3120.18120.93121.18112.88113.31112.19
O1-C1-N1123.28122.53123.17121.95122.79122.19
Table
Table 4. NBO atomic charges calculated for dacarbazine and copper(II) and cobalt(II) complexes of dacarbazine.
Table 4. NBO atomic charges calculated for dacarbazine and copper(II) and cobalt(II) complexes of dacarbazine.
DAC1DAC2Copper ComplexCobalt Complex
O1−0.605−0.805−0.607−0.582
N1−0.820−0.627−0.775−0.764
C10.6350.6350.6470.648
C20.0210.0130.008−0.002
C30.2830.3070.3140.321
N2−0.555−0.552−0.514−0.513
C40.2190.2200.2790.276
N3−0.444−0.509−0.518−0.470
qring−0.476−0.521−0.431−0.388
N4−0.346−0.296−0.343−0.343
N5−0.034−0.047−0.032−0.031
N6−0.248−0.246−0.2150.213
C5−0.390−0.392−0.390−0.390
C6−0.349−0.349−0.349−0.349
Table
Table 5. The values of energy of HOMO and LUMO orbitals and other electronic parameters calculated for dacarbazine and copper(II) and cobalt(II) complexes of dacarbazine.
Table 5. The values of energy of HOMO and LUMO orbitals and other electronic parameters calculated for dacarbazine and copper(II) and cobalt(II) complexes of dacarbazine.
DAC1DAC2Copper ComplexCobalt Complex
HOMO−7.8342−7.8391−4.7827−7.1092
LUMO−4.7772−4.7851−4.4066−5.5198
Energy gap3.0573.0540.37611.5894
Ionization potential7.83427.83914.78277.1092
Electron affinity4.77724.78514.40665.5198
Electronegativity6.30576.31214.594656.3145
Chemical potential−6.3057−6.3121−4.59465−6.3145
Chemical hardness1.52851.5270.188050.7947
Chemical softness0.3271180.3274392.6588670.629168
Electrophilicity index13.0068213.0460456.1308425.08677
Table
Table 6. Molecular docking analysis results.
Table 6. Molecular docking analysis results.
Free Energy of Binding (kcal/mol)Inhibition Constant (µM)Intermolecular EnergyTorsional EnergyUnbound Extended EnergyInteracting Residues
(H-Bonds) and Distance [Å]		Reference RMSD
DAC 1 −5.17163 −5.990.82−0.38N1-THY B:19 (1.92)
N2-THY A:7 (2.11)		25.45
DAC 2−5.13173.91−5.950.82−0.37N1-THY B:19 (1.96)
N2:THY A:7 (2.11)		25.43
Cu-DAC−9.310.15106 −10.40+1.10−0.63N1-THY:A8 (2.19)
N2-THY A:7 (1.86)		25.26
Co-DAC−9.330.14391−10.43+1.10−0.56N2-THY A:7 (1.86)
N1-THY B:19 (1.87)		26.06
Table
Table 7. The results of thermal decomposition of dacarbazine and their complexes with cobalt and copper.
Table 7. The results of thermal decomposition of dacarbazine and their complexes with cobalt and copper.
Compound Range of DecompositionWeight Loss (%)Products of Decomposition
Calcd.Found
Dacarbazine205–21047.7048.50Imidazole +C*
210–660100100-
CoCl2(Dac)2·75CH3OH60–120
175–630		10.19
86.38		10.34
85.61		CoCl2(Dac)2
CoO
CuCl2(Dac)2·1.5CH3OH60–110
160–580		8.78
85.40		8.39
87.16		CuCl2(Dac)2
CuO
C*—carbon residue.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Świderski, G.; Łaźny, R.; Sienkiewicz, M.; Kalinowska, M.; Świsłocka, R.; Acar, A.O.; Golonko, A.; Matejczyk, M.; Lewandowski, W. Synthesis, Spectroscopic, and Theoretical Study of Copper and Cobalt Complexes with Dacarbazine. Materials 2021, 14, 3274. https://doi.org/10.3390/ma14123274

AMA Style

Świderski G, Łaźny R, Sienkiewicz M, Kalinowska M, Świsłocka R, Acar AO, Golonko A, Matejczyk M, Lewandowski W. Synthesis, Spectroscopic, and Theoretical Study of Copper and Cobalt Complexes with Dacarbazine. Materials. 2021; 14(12):3274. https://doi.org/10.3390/ma14123274

Chicago/Turabian Style

Świderski, Grzegorz, Ryszard Łaźny, Michał Sienkiewicz, Monika Kalinowska, Renata Świsłocka, Ali Osman Acar, Aleksandra Golonko, Marzena Matejczyk, and Włodzimierz Lewandowski. 2021. "Synthesis, Spectroscopic, and Theoretical Study of Copper and Cobalt Complexes with Dacarbazine" Materials 14, no. 12: 3274. https://doi.org/10.3390/ma14123274

APA Style

Świderski, G., Łaźny, R., Sienkiewicz, M., Kalinowska, M., Świsłocka, R., Acar, A. O., Golonko, A., Matejczyk, M., & Lewandowski, W. (2021). Synthesis, Spectroscopic, and Theoretical Study of Copper and Cobalt Complexes with Dacarbazine. Materials, 14(12), 3274. https://doi.org/10.3390/ma14123274

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop