key: cord-0843085-6flskoob
authors: Ghasemi, Liana; Behzad, Mahdi; Khaleghian, Ali; Abbasi, Alireza; Abedi, Anita
title: Synthesis and characterization of two new mixed‐ligand Cu(II) complexes of a tridentate NN'O type Schiff base ligand and N‐donor heterocyclic co‐ligands: In vitro anticancer assay, DNA/human leukemia/COVID‐19 molecular docking studies, and pharmacophore modeling
date: 2022-02-24
journal: Appl Organomet Chem
DOI: 10.1002/aoc.6639
sha: f2f6a752c2779636bb78c310218b33a45651915f
doc_id: 843085
cord_uid: 6flskoob

Two new mixed‐ligand complexes with general formula [Cu(SB)(L')]ClO(4) (1 and 2) were synthesized and characterized by different spectroscopic and analytical techniques including Fourier transform infrared (FT‐IR) and UV–Vis spectroscopy and elemental analyses. The SB ligand is an unsymmetrical tridentate NN'O type Schiff base ligand that was derived from the condensation of 1,2‐ethylenediamine and 5‐bromo‐2‐hydroxy‐3‐nitrobenzaldehyde. The L' ligand is pyridine in (1) and 2,2′‐dimethyl‐4,4′‐bithiazole (BTZ) in (2). Crystal structure of (2) was also obtained. The two complexes were used as anticancer agents against leukemia cancer cell line HL‐60 and showed considerable anticancer activity. The anticancer activity of these complexes was comparable with the standard drug 5‐fluorouracil (5‐FU). Molecular docking and pharmacophore studies were also performed on DNA (PDB:1BNA) and leukemia inhibitor factor (LIF) (PDB:1EMR) to further investigate the anticancer and anti‐COVID activity of these complexes. The molecular docking results against DNA revealed that (1) preferentially binds to the major groove of DNA receptor whereas (2) binds to the minor groove. Complex (2) performed better with 1EMR. The experimental and theoretical results showed good correlation. Molecular docking and pharmacophore studies were also applied to study the interactions between the synthesized complexes and SARS‐CoV‐2 virus receptor protein (PDB ID:6LU7). The results revealed that complex (2) had better interaction than (1), the free ligands (SB and BTZ), and the standard drug favipiravir.

Schiff base ligands and their transition metal complexes have played a great role in the design and development of metallodrugs. [1] [2] [3] [4] [5] [6] These ligands and complexes are, to some extent, easily synthesized and their structural and electronic features could be smoothly tuned by the presence of various substituents. From a different point of view, mixed-ligand strategy also provides a new route to the synthesis of unique complexes with tuned stereo-electronic characteristics. In this regard, the synthesis of mixed-ligand Schiff base complexes could provide a distinctive opportunity for the study of a variety of complexes with special applications.

Thiazoles are heterocyclic compounds that contain both sulfur and nitrogen in the same ring. Thiazoles are bioactive materials that are found in a wide range of drugs and natural molecules. [7] [8] [9] [10] Thiamine is perhaps the most well-known compound containing thiazole ring, [11] but other important compounds such as dasatinib (used for treatment of leukemia), [12] bleomycin and its analogues (used for treatment of different cancers), [13] sulfathiazole and sulfamethizole (used as antibiotics), [14] and pramipexole (antidepressant) [15] are also known.

We have previously synthesized a series of mixedligand Schiff base complexes [16] [17] [18] and have studied their antibacterial properties. Those complexes had the general formula [Cu(SB)(L')]ClO 4 , in which the SBs were unsymmetrical tridentate NN'O type Schiff base ligands and the L's were either pyridine or bipyridine/ phenanthroline derivatives. The interesting features of the Schiff base ligands, along with the bithiazole ligands and their complexes, have motivated us to synthesize a series of new complexes that contain both of these ligands. Hence, herein, we report the synthesis and characterization of two new mixed-ligand Schiff base complexes with the same general formula, in which the SB is the Schiff base ligand derived from the condensation of 5-bromo 2-hydroxy 3-nitrobenzaldehyde with ethylenediamine and the co-ligand L' is pyridine in (1) or 2,2 0 -dimethyl-4,4 0 -bithiazole (which will be abbreviated as Btz) in (2) (Scheme 1). The crystal structure of (2) was also obtained. The two complexes were studied as anticancer agent against HL-60 leukemia cell line and showed considerable activity that is comparable with the standard drug 5-fluorouracil (5-FU). Molecular docking and pharmacophore studies were also performed to further investigate the anticancer properties of the studied complexes and their corresponding free ligands. The theoretical and experimental results show good correlation.

From a different point, designing new drugs for the COVID-19 pandemic and the investigation of their inhibitory action could effectively accelerate the drug discovery. [19] SARS-CoV-2, also known as 2019-nCoV, is a novel coronavirus that originated in Wuhan, China, and swept the globe in late 2019. [20] [21] [22] COVID-19 was the disease caused by the SARS-CoV-2 virus, and it was similar to diseases caused by past coronavirus epidemics, such as SARS-CoV [23] and MERS. [24] It presents as a dry cough, fever, exhaustion, and a loss of smell or taste in the respiratory system. In this regard, herein, we report the use of molecular docking and pharmacophore studies to investigate the interaction of the synthesized complexes and the corresponding free ligands with SARS-CoV-2 virus receptor protein (PDB ID:6LU7). [25] [26] [27] [28] The results revealed that complex (2) had better interaction than (1), the free ligands (SB and BTZ), and the standard drug favipiravir.

2 | EXPERIMENTAL 2.1 | Materials and instrumentation 1,2-Ethylenediamine, 5-bromo-2-hydroxybenzaldehyde, and the solvents were purchased from Merck Millipore. Cu(ClO 4 ) 2 Á6H 2 O was purchased from Sigma-Aldrich. These chemicals were used as received. 5-Bromo-2-hydroxy-3-nitrobenzaldehyde was synthesized as described elsewhere. [29, 30] 2,2 0 -Dimethyl-4,4 0 -bithiazole was also synthesized following a literature procedure. [31] Fourier transform infrared (FT-IR) spectra were obtained as KBr pellets by a Bruker FT-IR instrument. UV-Vis spectra were recorded by a Shimadzu UV-1650 PC spectrophotometer. Elemental analyses were performed on a Perkin-Elmer 2400II CHNS-O analyzer.

2.1.1 | X-ray crystallography SCXRD data were collected on a MAR345 dtb diffractometer equipped with an image plate detector using Mo-Kα X-ray radiation. The structure was solved by direct methods using SHELXS-97 and refined using the full-matrix least-squares method on F 2 , SHELXL-2016. All non-hydrogen atoms were refined anisotropically. The structure was deposited in Cambridge Crystallographic Data Center with deposition number CCDC 2089608.

The human leukemia cancer cells (HL-60) were maintained in RPMI 1640 media supplemented with 10% of fetal bovine serum and a pen strep antibiotic (100 μg ml À1 of penicillin and 100 μg ml À1 of streptomycin), and this test was performed at 37 C with a humid atmosphere in a CO 2 incubator. When cancer cells reached 80% confluence, they could be used to treat cancer cells. To perform these experiments, we first prepared stock solutions of the complexes that were diluted in RPMI 1640. [32, 33] For cytotoxicity assessments, cancer cells were incubated at 37 C and were plated in 96-well plates overnight. Then, 200 μl of the mixture of compounds and culture medium with different concentrations (100, 50, 25, 12.5, and 6.25 μM) were plated in 96-well plates to investigate inhibition of cell growth for 48 h. These experiments were replicated three times. After the end of treatments and complication incubation, 10 μl of MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazoliumbromide) (5 mg ml À1 ) solution in phosphate buffered saline (PBS) (5 mg ml À1 ) was added to each well. Then, the 96-well plates were covered with aluminum foil and incubated for 3 h at 37 C. After draining the culture medium, 100 μl of dimethyl sulfoxide (DMSO) was added to 96-well plates. Absorbance values were determined at 570 nm using an enzymelinked immunosorbent assay (ELISA) reader. Cell viability is expressed using control wells as 100% viability and their absorption values are expressed. The IC 50 values were evaluated for the sensitivity of treated compounds. [34, 35] 2.1.3 | Docking studies Molecular docking simulation was used for further studying the anticancer activity of the complexes as well as their possible anti-COVID-19 activity. The crystal structures of 1BNA (BDNA Dodecamer: righthanded double-stranded B helix as rigid molecule with the sequence [5 0 -D(CpGpCpGpApApTpTpCpGpCpG)-3 0 ]), human leukemia inhibitor factor (LIF) as receptor (PDB ID:1emr as rigid molecule), and SARS-CoV-2 virus receptor (PDB ID:6LU7 as rigid molecule) protein were downloaded from PubChem. The complex (1) and the free ligands (SB) and (BTZ) were optimized via standard 6-311G** basis sets that were used for C, H, N, O, and Br atoms whereas the LANL2DZ basis set along with the effective core potential (ECP) functions was employed for Cu. Crystallographic data for complex (2) were taken as a CIF file and converted to the PDB format using Mercury software. The molecular docking simulation and calculations were performed by AutoDock 4.2 (The Scripps Research Institute, La Jolla, CA, USA) with AutoGrid 4, meaning all non-ring torsions were maintained. Primarily, the heteroatoms including water molecules around the duplex were removed using AutoDock tools; Kollman united atom type charges, polar hydrogen atoms, and Gasteiger partial charges were then added to the receptor molecule. All the docking simulations were defined by the use of a grid box with 92 Â 124 Â 126 Å points with a grid-point spacing of 0.375 Å for BNA, a grid box with a size of 122 Â 98 Â 126 Å and a grid-point spacing of 0.375 Å for LIF, and a grid box with a size of 100 Â 86 Â 126 Å and a grid-point spacing of 0.758 Å for 6LU7. In this study, we used a Lamarckian genetic algorithm method. The number of evaluations and the number of genetic algorithm runs were set to 200. We examined the structures and selected the structures with the lowest energy from similar structures. After that, the interactions of BNA, 1EMR, and 6LU7 and their binding modes with compounds were analyzed using the AutoDock program, UCSF Chimera 1.5.1 software, Discovery Studio 3.0 from Accelrys, and DS Visualizer. [36] [37] [38] [39] Morever, docking studies were performed on all of the complexes and their corresponding ligands for three times and the average ± SD are reported.

To explain the interactions between the synthesized complexes and standard drugs with 1EMR, DNA, and COVID-19, a pharmacophore modeling was conducted. Controlling DNA and protein are key components for successful treatment, particularly in cancer cells. Pharmit link (http://pharmit.csb.edu) created this interaction profile. Moreover, the interaction bonds were categorized as H-bond acceptor (H-Acc.), H-bond donor (H-Don.), and hydrophobic (Hyd). [40] 2.2 | Synthesis of the complexes

Synthesis of unsymmetrical Schiff base ligands with symmetrical aliphatic diamines such as ethylenediamine, propanediamine, and other symmetrical aliphatic diamines is extremely difficult. Attempts to synthesize such unsymmetrical so-called half-units usually result in symmetrical tetradentate Schiff base ligands, making the synthesis of their corresponding complexes even more difficult. Many researchers have attempted to discover new synthetic methods for such unsymmetrical Schiff base complexes. We were unable to synthesize free unsymmetrical Schiff base ligands from ethylenediamine, but we were able to apply a template technique that works well for the synthesis of Cu(II) complexes containing the required unsymmetrical NN'O type ligand. The related co-ligand, pyridine, could be easily replaced with chelating co-ligands, allowing for the creation of a wide range of Cu(II) complexes.

In a typical experiment, 1.23 g of 5-bromo-2-hydroxy-3-nitrobenzaldehyde (5.00 mmol) was dissolved in 30 ml of methanol in a round-bottom flask. While this solution was being continuously stirred at room temperature, 5 ml of an aqueous solution of Cu(ClO 4 ) 2 Á6H 2 O (5.00 mmol, 1.85 g) was slowly added, followed by the addition of 0.80 g (10.00 mmol) of pyridine. The reaction mixture was stirred for 1 h. A total of 0.30 g of 1,2-ethylenediamine (5.00 mmol) dissolved in 5 ml of methanol was then added dropwise to the reaction mixture. This reaction mixture was further stirred for 3 h without heating and the reaction progress was followed by thin-layer chromatography (TLC). The green powder precipitate was collected by filtration and washed with diethyl ether and air dried. The yield was 90%. The powder was recrystallized from methanol and needle shaped green pure polycrystals were collected by filtration and air dried. Selected infrared (IR) ( (2) A total of 0.14 g (0.33 mmol) of (1) was dissolved in 20 ml of methanol in a round-bottom flask while being magnetically stirred. Afterwards, 20 ml of a methanolic solution of 2,2 0 -dimethyl-4,4 0 -bithiazole (0.66 mmol, 0.12 g) was slowly added to the reaction mixture. The obtained reaction mixture was then heated at 50 C. The reaction progress was checked by TLC, and after 3 h, the dark green powder precipitate was collected by filtration and was washed with 10 ml of diethyl ether and air dried. This solid was recrystallized from methanol/acetonitrile 1:1 solution. 

The synthesis procedure is shown in Scheme 1. The synthesis of unsymmetrical Schiff base ligands from aliphatic diamines is extremely difficult and usually results in symmetrical ligands even when the aldehyde:amine ratio is 1:1 at low temperatures. Hence, we used a modified template method to synthesize the complexes. Complex (1) was prepared by the reaction between ethylenediamine, 5-bromo-2-hydroxy-3-nitrobenzaldehyde, Cu(ClO 4 ) 2 Á6H 2 O, and pyridine in methanol at room temperature. Complex (2) was readily prepared by ligand substitution of the monodentate pyridine co-ligand by bidentate N-donor heterocyclic co-ligand (Btz) in methanol at room temperature.

Dark green single crystals were obtained by slow evaporation of methanol/acetonitrile (1:1) solution after several days. Table 1 summarizes the crystallographic data for (2) . Selected bond lengths and bond angles around the central metal ion are also collected in Table 2 . Figure 1 shows the labeled molecular structure of complex (2) with atom numbering scheme. In the structure of the cationic part of the complexes, an NN'O type unsymmetrical Schiff base ligand exists along with a 2,2 0 -dimethyl-4,4 0 -bithiazole co-ligand. The Schiff base ligand is the monocondensed form of the reaction between 1,2-ethylenediamine and 5-bromo-2-hydroxy-3-nitrobenzaldehyde. In this complex, one of the NH 2 groups of ethylenediamine moiety is directly coordinated to the metal center, and the other NH 2 group is condensed to 5-bromo-2-hydroxy-3-nitrobenzaldehyde, forming the unsymmetrical Schiff base ligand. The Schiff base ligand is deprotonated from the O atom of the salicylaldehyde moiety and has a (À1) charge, Hence, the Cu 2+ central metal ion's charge is balanced by the presence of one perchlorate anion and the overall charge of the complex is neutral. The central metal ion in (2) is penta-coordinated and the geometry around the Cu(II) ion is slightly distorted square-based pyramid (SBP). The trigonality index (τ) that was calculated from τ ¼ β À α ð Þ=60 is equal to 0.197. [41] In this equation, β is the greatest angel and α is the second greatest angel. The τ values range from 0 to 1, which corresponds to ideal SBP and trigonal bipyramidal (TBP) structures, respectively. Three of the basal positions in the SBP structure are occupied by the NN'O coordinating atoms of the Schiff base ligand, including N atom from the NH 2 group, the N 0 from the azomethine group, and the O atom from the phenoxide group. The fourth basal position is occupied by the N atom from the bithiazole ligand. The apical position is also occupied by the other N atom of the bithiazole ligand.

As could be seen from Table 2 , the N1 Cu1 N3 and O1 Cu1 N4 angles are 175.01 and 163.21 , respectively, which shows deviation from ideal SBP structure. The Cu(II) ion is also located about 0.13 Å above the plane defined by the four basal coordinating atoms, which further confirms the mentioned slight deviation from SBP structure. These data are similar to previously published complexes with the same scaffold. [16] [17] [18] Besides, there are two moderate and two weak hydrogen bonds (H. bonds) (N4-H4A … O5, N4-H4B … O4, and C11-H11 … O4, C15-H15A … O6 with distances D … A of 3.198, 3.220, and 3.496, 3.411 Å, respectively), together with one relatively weak halogen bond of Br1 … O7, 3.263 Å, between perchlorate ions and complexes, stabilizing the 3D crystal structure.

As shown in Scheme 1, the three most specific functional groups in the complexes are azomethine (C═N), NH 2 , and perchlorate anion. The FT-IR spectrum of the complexes showed signals related to these groups. The sharp bands at 1629 and 1635 cm À1 for complexes (1) and (2) were assigned to the (C═N) stretching vibration. The two bands at around 3300 and 3200 cm À1 correspond to the asymmetric and symmetric stretching vibrations of the NH 2 group, respectively. A sharp band at around 1100 cm À1 was also assigned to the stretching vibration of the uncoordinated perchlorate anion. [42] Other signals such as δ(CCH) of the aromatic rings, aliphatic (CCH) (14) and δ(C═C) were observed at around 3100, 3000, and 1600-1400 cm À1 . In the UV-Vis spectra of these complexes, π ! π* transitions of the phenyl rings and azomethine groups are observed as two intense bands at about 231 and 264 nm. The other signals above 392 nm could be assigned as the ligand to metal charge transfer (LMCT) transition. The very weak signals above 600 nm, which could only be observed at higher concentrations, were assigned to the d ! d transition. [43] [44] [45] [46] [47] [48] 3.4 | Biological studies

The HL-60 cell line is a human leukemia cell line that has been used to study blood cell development and physiology in the laboratory. The early detection and treatment of leukemia cancer with appropriate medications such as Schiff base compounds could help reduce the number of leukemia cancer victims. To investigate cytotoxic effects, the compounds were incubated for 48 h against HL-60 cancer cell line. The charts of dosedependent effects of complexes on the viability of each HL-60 cell line are presented in Figure 2 . A logical relationship was found between IC 50 and concentration values. In this figure, the cell viability decreased with the increasing concentration of the complexes. These complexes showed good cytotoxic and dose-dependent activity against HL-60 cancer cells. Although at higher concentrations, the performance of the two complexes was almost identical, it could be seen that complex (2), which is penta-coordinated with BTZ ligand, can slightly better destroy cancer cells at lower concentrations. The standard drug was 5-FU. As you can see in Table 3 , the inhibitory concentration of complexes against the HL-60 cancer cell line and the IC 50 value for the complexes were almost similar. The results are also comparable with the standard drug 5-FU. [49, 50] F I G U R E 1 Molecular structure of (2) with atom numbering scheme. Thermal ellipsoids are drawn at 50% probability level. [51, 52] Computer-based simulations help us predict the pharmacological action of drugs at the molecular level. [26, 53, 54] The lowest energy conformations were rated based on docking calculations of the lowest free binding energy. The docking results are shown in Figure 3 and are listed in Table 4 . The minimum energy of each compound was applied inside the DNA double helix. The standard anticancer chemotherapeutic drugs including cyclophosphamide, 5-FU, and doxorubicin were also used. The calculated binding free energies (ΔG binding ) for the complexes (1 and 2), ligands (SB and BTZ), and drugs to DNA (cyclophosphamide, 5-FU, and doxorubicin) were À10.10 ± 0.01, À8.39 ± 0.01, À7.95 ± 0.14, À6.37 ± 0.01, À4.66 ± 0.18, À10.90 ± 0.06, and À7.00 ± 0.01 kcal mol À1 , respectively, and the binding of the compounds to the DNA was spontaneous. Comparison of the ΔG binding reveals the higher tendency of (1) to bind to DNA among the studied drugs. [55] The final energy value was calculated from vdW + H-bond + desolve energy (kcal mol À1 ), which is related to the formation of H. bonds between the compounds and the receptor. As could be seen from Figure 3 , complex (1) and BTZ ligand are located in the major groove of DNA whereas complex (2) and SB ligand are located in the minor groove. The studied standard drugs are also in the major groove and are more compatible with complex (1) in binding with DNA. This could be the reason for slightly better activity of complex (1) than (2). The hydrogen-bonding interactions of the compounds and binding sites are illustrated in Figure 4 . There are four H. bonding interactions between the NH 2 group of complex (1) The drug-receptor residue interactions between the studied complexes, ligands, and standard drugs, with DNA, are summarized in Table 5 . The binding of the compounds to the DNA was spontaneous. According to data from density functional theory, positive electron density on the surface of complexes facilitates the metal drug's access to DNA phosphate groups. [56] The DNA docking data for (1) and (2) indicated that both complexes had a high affinity for binding to the target DNA duplex. In consideration of the interaction strength, however, complex (1) had higher affinity (À10.10 ± 0.01 kcal mol À1 ) than complex (2) (À8.39 ± 0.01 kcal mol À1 ). Furthermore, (1) interacted with the DNA binding sites via a H. bond and π-alkyl hydrophobic interactions. On the other hand, (2) was linked to DNA via π-sulfur and H. bond interactions. The ligand binding energies were also À7.95 ± 0.14 and À6.37 ± 0.01 for the SB and BTZ, respectively. Compared with copper complexes reported in literature, we found that the presence of withdrawing substitutions such as bromine and nitro groups shows better binding energy, compared with the electron-donating groups. [36, 57] 3.4.3 | Molecular docking with human LIF protein as receptor (PDB ID:1emr) LIF was first purified and cloned because of its capacity to induce the differentiation of the M1 murine myeloid leukemia cell line. [58] [59] [60] [61] [62] In this study, the docking calculations were performed to describe the mode of binding to LIF. The values of docking energies were À7.20 ± 0.04, À8.43 ± 0.01, À7.07 ± 0.20, À5.52 ± 0.01, À7.39 ± 0.08, and À4.60 ± 0.01 kcal mol À1 for (1), (2), (SB), (BTZ), lenalidomide, and 5-FU, respectively ( Figure 5 and Table 6 ). Comparing the relative binding energies (ΔG binding ) proved that complex (2) was more effective than (1) 

Estimated free binding energy [a] À10. 10 In ligand (BTZ), two H. bondings between S atom to THR6 and ASN97 were found at 2.85 and 2.73 Å. One π-Don. H. bond between a ring of the ligand to SER94 was found at 2.76 Å. One π-σ-Hyd. bond between a ring of the ligand to VAL65 was found at 3.79 Å. Four alkyl Hyd. bondings between C of the ligand and ILE98, LEU166, VAL65, and LEU162 were found at 4.03, 5.14, 3.93, and 4.95 Å. Four π-alkyl Hyd. bondings between a ring of the ligand to PRO149 (two), LEU166, and LEU162 were found at 5.06, 4.93, 5.28, and 4.94 Å, respectively. Moreover, further pi-alkyl bonds were observed for (1), cyclophosphamide, and doxorubicin. 5-FU showed π-hydrophobic bonds and (2) showed pi-sulfur bonds. Molecular docking studies by BDNA and 1emr for (1) and (2) proved that the complexes could be potent anticancer agents. The higher tendency of (1) to interact with BDNA and (2) with 1mer may clarify the similar anticancer activity of these complexes but with a different mechanism of actions. The drug-receptor residue interactions between complexes and standard drugs, with 1EMR, are shown in Table 7 . As you can see in Table 7 , H. bonds and alkyl hydrophobic interactions result in better binding of complexes to proteins. Complex (2) has more protein interactions and also has more negative binding energy. [63] 3.4.4 | Molecular docking with SARS-CoV-2 virus (COVID-19 disease) as receptor (PDB ID:6LU7) protein Molecular docking was used to investigate how the complexes interacted with 6LU7 as the receptor protein. The values of docking energies for complexes (1), (2), (SB), (BTZ), and favipiravir were À7.58 ± 0.04, À8.28 ± 0.07, À6.25 ± 0.22, À4.98 ± 0.01, and À3.77 ± 0.06 kcal mol À1 , respectively (Table 8 ). Comparing the relative binding energies (ΔG binding ) proved that complex (2) was more effective than (1) and the standard drug favipiravir. The order of the binding energies was (2) > (1) > (SB) > (BTZ). The drug-receptor residue interactions between the complexes and favipiravir with 6LU7 are shown in Table 9 . In complex (1), three H. bondings between one of the H atoms of the NH 2 group as well as the H15 and H16 to ASN142, GLU166 were found at 1.91, 2.12, and 2.06 Å. C-H. bonding between C pyridine and ASN142 was also found at 3.24 Å. Alkyl hydrophobic bonding between phenyl ring from aldehyde moiety and MET165 was observed at 5.24 Å and another alkyl hydrophobic bonding to HIS163 was also found at 5.29 Å. Br to ASP 187 interaction was found at 3.53 Å. In complex (2) 1.77 Å. One π-lone pair bonding between a ring of ligand and TYR101 was found at 2.99 Å. One alkyl Hyd. and π-alkyl Hyd. to LYS102 were found at 4.24 and 4.75 Å. In (BTZ) ligand, C-H. bond between the N atom to VAL18 was found at 3.78 Å. One π-donor-H. bond between a ring of ligand to GLY71 was found at 2.37 Å. One π-lone pair bond to MET17 was found. Two alkyl Hyd. bonds between C atom to ALA70 and VAL18 were found at 3.49 and 5.40 Å. Two π-alkyl Hyd. bonds to TRP31 and ALA70 were found at 4.71 and 4.77 Å. The graphical representation of the molecular docking, that is, drugs-receptor H. bond interactions, receptor side surface H. bond interactions, and 2D-ligand-receptor H. bond interactions, was shown for all of the studied materials [64, 65] (Figure 6a-c) . Compared with the other copper complexes reported in the literature, [66] these complexes showed more negative binding energies, which could be attributed to the electron-withdrawing groups of the aldehyde and BTZ moieties.

F I G U R E 7 Interaction profile according to pharmacophore study of the complexes (1), (2), (SB), and (BTZ)

Using the AutoDock 4.2 program, pharmacophore modeling for the Schiff base drugs was created. To find pharmacophore modeling, we used the pharmit website (http://pharmit.csb.edu). (1) and (2) were saved as PDB files and then imported into the pharmit as a ligand to be docked with DNA, 6LU7, and 1EMR. All of the studied chemicals interacted well with DNA, 6LU7, and 1EMR protein pockets as shown in Figure 7 . The number of bonds counted with 1EMR, 6LU7, and DNA appeared as follows [67] : From pharmacophore results, it was found that the SB ligand has high number of interactions with DNA and 1EMR. These results suggest that the SB ligands and their corresponding complexes may play an important rule in the drug design for anticancer purposes. The results for 6LU7 suggest that the complexes acted better than the ligands and the presence of the SB ligands are again important. Table 10 compares the binding energies of other comparable complexes to 6LU7 to demonstrate the merit of this work. Table 10 shows that our results are similar to, and in some cases superior than, those previously published.

Two new mixed-ligand complexes containing an NN'O type unsymmetrical tridentate Schiff base ligand as the main ligand and pyridine (1) or 2,2 0 -dimethyl-4,4 0bithiazole (BTZ) co-ligands (2) were synthesized and characterized. These two complexes were used as anticancer agents against leukemia cancer cell line HL-60. Molecular docking and pharmacophore studies were performed on BDNA and LIF to further investigate the anticancer activities of the complexes. The same computational studies were performed on SARS-CoV-2 virus receptor protein (PDB ID:6LU7).

From molecular docking studies, the obtained data reveal the following remarks:

1. Complex (1) showed better docking score with DNA (À10.10 ± 0.005). On the other hand, complex (2) showed better results with 1EMR (À8.43 ± 0.005 kcal mol À1 ) and with COVID-19 (À8.28 ± 0.070 kcal mol À1 ). The order of the binding best energy for (1) is as follows: DNA > 6LU7 > 1EMR, and for (2) À6. 18 69 interacted with protein pockets by predominantly Hdon. and H-acc. interactions. 3. The C═N groups and the BTZ both contributed to stable docking positions. 4. Finally, this in silico approach confirmed the activity of Cu(II) complexes in controlling human cancer cells, which agreed with in vitro results. 5. The anticancer activity of these complexes was comparable with the standard drug 5-FU. 6. The results of pharmacophore are consistent with the results of molecular docking and experimental results, and complex (2) showed the best performance.

Press, 1250

Polyhedron

Mater. Today: Proc

SUPPORTING INFORMATION Additional supporting information may be found in the online version of the article at the publisher's website. How to cite this article

ORCID Liana Ghasemi https://orcid.org/0000-0002-4175-3661 Mahdi Behzad https://orcid.org/0000-0001-5957-0778 Ali Khaleghian https://orcid.org/0000-0001-8520-6654 Alireza Abbasi https://orcid.org/0000-0002-0331-0202 Anita Abedi https://orcid.org/0000-0001-7562-6241