Synthesis, Structure, and Biologic Activity of Some Copper, Nickel, Cobalt, and Zinc Complexes with 2-Formylpyridine N4-Allylthiosemicarbazone

A series of zinc(II) ([Zn(H2O)(L)Cl] (1)), copper (II) ([Cu(L)Cl] (2), [Cu(L)Br] (3), [Cu2(L)2(CH3COO)2]·4H2O (4)), nickel(II) ([Ni(HL)2]Cl2·H2O (5)), and cobalt(III) ([Co(L)2]Cl (6)) complexes were obtained with 2-formylpyridine N4-allylthiosemicarbazone (HL). In addition another two thiosemicarbazones (3-formylpyridine N4-allylthiosemicarbazone (HLa) and 4-formylpyridine N4-allylthiosemicarbazone (HLb)) have been obtained. The synthesized thiosemicarbazones have been studied using 1H and 13C NMR spectroscopy, IR spectroscopy, and X-ray diffraction analysis. The composition and structure of complexes were studied using elemental analysis, IR and UV-Vis spectroscopies, molar conductivity, and magnetic susceptibility measurements. Single crystal X-ray diffraction analysis elucidated the structure of thiosemicarbazones HL, HLa, and HLb, as well as complexes 4 and 5. The antiproliferative properties of these compounds toward a series of cancer cell lines (HL-60, HeLa, BxPC-3, RD) and a normal cell line (MDCK) have been investigated. The nickel complex shows high selectivity (SI > 1000) toward HL-60 cell line and is the least toxic. The zinc complex shows the highest selectivity toward RD cell line (SI = 640). The copper complexes (2–4) are the most active molecular inhibitors of proliferation of cancer cells, but exhibit not such a high selectivity and are significantly more toxic. Zinc and copper complexes manifest high antibacterial activity. It was found that calculated at B3LYP level of theory different reactivity descriptors of studied compounds strongly correlate with their biological activity.


Introduction
e steady increase in the number of diseases poses a challenge for modern chemistry in the search for new biologically active substances that will allow to treat these diseases and enhance the duration and quality of life [1][2][3][4][5][6].
e number of cancers continues to rise [7]. erefore, a great attention is paid to the new anticancer compounds [8,9]. Because the main disadvantage of substances, which are used for cancer treatment, is their general toxicity for the human organism, special attention should be paid to the study of selectivity of antiproliferative activity and the toxicity of these substances. e expansion of the arsenal of such substances will allow in the future to personally select the active substance for the treatment of this disease, which should increase the effectiveness of treatment [10].
Copper complexes are the most studied among the transition metal complexes with thiosemicarbazones from the biological point of view. In many cases, the coordination of thiosemicarbazones to the copper(II) ion leads to the greatest enhancement of anticancer and antimicrobial activity [31][32][33][34]. Despite the fact that the selectivity of the described copper complexes in many cases exceeds the selectivity of anticancer drugs such as doxorubicin and cisplatin, in most cases the selectivity index does not exceed 10 [35]. Among N 4 -substituted thiosemicarbazones more attention is paid to N 4 -alkyl, N 4 -arylthiosemicarbazones. N 4allylthiosemicarbazones are described less frequently. However, N 4 -allylthiosemicarbazones exhibit high biological activity and are comparatively more soluble in water. Based on this, it is important to continue the systematic study of thiosemicarbazones and biometal complexes with them in order to identify new substances with high anticancer, antimicrobial, and antifungal properties, which at the same time will be free from such a significant drawback as high toxicity.
In order to find such substances, the present study describes the chemical synthesis and characterization of Zn(II), Cu(II), Ni(II), and Co(III) complexes with 2-formylpyridine N 4 -allylthiosemicarbazone (Scheme 1), in vitro biological evaluation, and in vivo toxicity of these substances.

Synthesis of the Complexes 1-6.
e direct interaction of biometal salts with the ethanolic solutions of 3-formylpyridine N 4 -allylthiosemicarbazone was used to obtain complexes 1-6.

Synthesis of the Complex [Zn(H 2 O)(L)Cl] (1).
2-Formylpyridine N 4 -allylthiosemicarbazone (HL) (0.220 g, 1 mmol) was dissolved in 25 mL of ethanol upon heating (55°C). After that, a sample of solid zinc chloride (ZnCl 2 ) (0.136 g, 1 mmol) was added to the hot solution. e resultant mixture was stirred for 40 min at 55°C. e yellow precipitate of complex 1 was obtained upon cooling to room temperature. It was filtered out, washed with cold ethanol, and dried in vacuo. Yellow

Synthesis of the Complex [Cu(L)Cl] (2).
e complex 2 has been obtained as described by Zeglis B.M. et al. [39]. Solid copper(II) chloride dihydrate (0.170 g; 1 mmol) was added to a hot ethanolic solution of 2-formylpyridine N 4allylthiosemicarbazone (HL) (0.220 g; 1 mmol). e resultant green solution was stirred and heated (55°C) for 40 min. Green precipitate was formed upon cooling to room temperature. It was filtered out, washed with cold ethanol, and dried in vacuo.

Physical Measurements.
e chemical elemental analysis, melting points of the thiosemicarbazones, molar conductivity values, and magnetic susceptibilities were Bioinorganic Chemistry and Applications determined by specific methods [40]. IR spectra of the studied organic and complexes were recorded on a Bruker ALPHA FTIR spectrophotometer. UV-Vis spectra were measured on Agilent Cary 300 UV-Vis spectrophotometer in 0.1 mM DMSO solutions. e NMR ( 1 H and 13 C) spectra of the studied thiosemicarbazones (HL, HL a , HL b ) were recorded on a Bruker DRX-400 spectrometer (Billerica, MA, USA) DMSO-d 6 and acetone-d 6 solutions.
EPR data were recorded using an Elexsys E 500 Bruker spectrometer, operating at a microwave frequency of approximately 9.47 GHz. Spectra were recorded using a microwave power of 10 mW across a sweep width of 200 mT (centered at 320 mT) with modulation amplitude of 0.4 mT. Experiments were carried out at 110 K using a liquid nitrogen cryostat. About 10 μmol of the corresponding compound were dissolved in 1 mL of DMSO to obtain a 10-mM stock solution. e 0.1, 0.2, and 0.5 mM solutions were prepared by dilution of the corresponding stock solution with DMSO. For low concentrations, several scans (4 to 6) were averaged.

X-Ray Crystallography.
e single crystal X-ray analysis of three thiosemicarbazones (HL, HL a and HL b ) and complexes 4 and 5 were carried out on a Xcalibur E CCD diffractometer equipped with a CCD area detector and a graphite monochromator, MoKα radiation, at room temperature. Data collection and reduction, and unit cell determination were done by CrysAlis PRO CCD (Oxford Diffraction); the SHELXS97 and SHELXL2014 program packages [41,42] were used to solve and refine the structures. e non-hydrogen atoms were treated anisotropically (fullmatrix least squares method on F 2 ). e hydrogen atoms were placed in calculated positions and were treated using riding model approximations with Uiso(H) � 1.2Ueq(C), while the oxygen bounded H-atoms were found from differential Fourier maps at an intermediate stage of the structure refinement. ese hydrogen atoms were refined with the isotropic displacement parameter Uiso(H) � 1.5Ueq(O). e X-ray data and the details of the refinement of HL, HL a , HL b , 4 and 5 are summarized in Table 1, and the selected bond lengths are given in Table 2. e selected bond lengths for HL a , HL b , and hydrogen bond parameters are given in Tables S1, S2. e geometric parameters were calculated by PLATON program [43] and Mercury software [44] was used for the visualization of structures. e hydrogen atoms that were not involved in the hydrogen bonding were omitted from the generation of the packing diagrams.

Biological Studies.
e antiproliferative activity of the studied substances was determined toward a series of cancer cell lines (HL-60 (ATCC CCL-240), BxPC-3 (ATCC CRL-1687), HeLa (ATCC CCL-2), RD (ATCC CCL-136)), and a normal cell line (MDCK (ATCC CCL-34)). ese cell lines were cultured in accordance with the previously described method [45]. Determination of the number of viable cells after cultivation in the presence of various concentrations of the studied substances was carried out by the MTS and Resazurin assays according to the previously described method [45]. e 10-mM stock solutions of the studied compounds as well as the reference compound (doxorubicin) were prepared by dissolving 10 μmol of the corresponding compound in 1 mL of DMSO. e 0.1, 1, 10, and 100 μM solutions were prepared by dilution of the stock solutions with the corresponding medium and used for the study of antiproliferative activity.
Antibacterial and antifungal properties of the studied substances were determined on a series of standard strains of Gram-positive microorganisms (Staphylococcus aureus (ATCC 25923)), Gram-negative microorganisms (Escherichia coli (ATCC 25922)), and fungi (Candida albicans (ATCC 10231)) using the method of serial dilutions in liquid broth. Stock solutions of the studied compounds as well as the reference compounds (furacillinum and nystatin) were prepared by dissolving 10 mg of the corresponding compound in 1 mL of DMSO. After that, the starting solutions (1 mg/mL) were prepared by tenfold dilution of the stock solution with nutrition medium (2% of peptonate bullion). e next dilutions were prepared by serial double dilution with the nutrition medium.
e ABTS •+ method [46] was applied for the determination of antiradical activity of the studied substances. Standard solutions of ABTS •+ radical cations and tested substances as well as the conditions of spectrophotometric measurements and calculations of inhibition were made as described [6]. e 10-mM stock solutions of the studied compounds as well as the reference compound (trolox) were prepared by dissolving 10 μmol of the corresponding compound in 1 mL of DMSO. e 1, 10, 100, and 1000 μM solutions were prepared by the dilution of stock solutions with DMSO. After that, 20 μL of each tested compound solutions was transferred in a 96-well microtiter plate, and 180 μL of ABTS •+ working solution was added to receive solutions with final concentrations of tested compounds 0.1, 1, 10, and 100, respectively.

Acute Toxicity Assay against Daphnia magna.
e general toxicity of the tested compounds was evaluated using Daphnia magna (Straus, 1820). Daphnia magna was originated from a culture maintained parthenogenetically at the Institute of Zoology (Laboratory of Systematics and Molecular Phylogeny). e test design was based on ISO 6341 : 2012. e testorganisms Daphnia magna Straus were fed with Chlorella vulgaris.
ese unicellular algae were grown using aseptic technology to exclude contamination of the culture by bacteria, algae or protozoa. Chlorella vulgaris were cultivated in Prat's growth medium containing KNO 3 in H 2 O distilled (adjusted the pH to 7.0, autoclaved and stored at 5°C).
Juveniles were selected according to their size and kept in fresh medium for 24 h. D. magna was cultured in Costar ® 24-well culture clear sterile multiple well plates covered by a lid to prevent the possibility of contamination and evaporation but at the same time to allow gaseous exchange between air and culture medium. Each well contained 10 daphnids in 1000 μL final volume of each dilution of the tested compounds. e bioassay was then repeated at the concentrations ranging from 0.1 to 100 μM (0.1, 1, 10, and 100 μM) in order to determine LC 50 for each compound, including doxorubicin (the positive control). Aqueous straw infusion growth media was used to dilute the stock solutions to the required concentrations. e final test solutions contained up to 0.1% DMSO and had a final volume of 1 mL. A 0.1% solution of DMSO in aerated medium (pH∼7.5 ± 0.2; O 2 ≥ 6.0 mg/L) was used as a negative control. roughout the experiment, the juvenile daphnids were incubated at 22 ± 2 0 C, using a 16 h/8 h light/dark cycle (500-1000 lx). e mobility (viability) of the test organisms was observed after the 24-h exposure. e experiment was performed in triplicate. e daphnids were considered immobilized only if they did not swim during the 15 s which follow gentle agitation of the test and control solutions, even if they could still move their antennae. e percentage of viability (V (%)) of Daphnia magna was calculated according to the formula:  (2) 10.121 (2) 12.5635 (6) 8.6275 (13) 19.3092 (8) 4.2978 (7) 10.003 (2) 10.6121 (7) 22.8273 (17) 16.5371 (7) (5)) 1.351 (7) 1.350 (5) C2-C3 1.442 (5) (1.465 (7)) 1.445 (7) 1.453 (7) V(%) � N -Number of viability of Daphnia magna. e compounds median lethal concentration that kills 50% of the juvenile daphnids (LC 50 ) values were calculated from the dose-response equation determined by the least squares fit method.

Results and Discussion
e thiosemicarbazone HL has been previously described by Zeglis B.M. et al. [39] and has been obtained using the described method. Its yield is 89%. It was characterized by FT-IR, UV-Vis, 1 H NMR, and 13 C NMR spectroscopy. Also, its structure as well as structures of other two thiosemicarbazones with different position of pyridine nitrogen atom (HL a and HL b ) has been determined using X-ray diffraction analysis. Complexes 1-6 were synthesized by the interaction of ethanolic solution of 2-formylpyridine N 4allylthiosemicarbazone (HL) and zinc chloride (1), copper(II) salts (2-4) in a 1 : 1 molar ratio or nickel(II) (5) and cobalt(II) (6) salts in a 2 : 1 molar ratio. e composition of the thiosemicarbazone HL and complexes 1-6 has been confirmed using elemental analyses data.  2)) core is practically planar within 0.06Å in these ligands and the dihedral angle between given core and pyridine ring ranges from 6.5 to 10. Meanwhile, the HL, HL a , and HL b are nonplanar because the C 3 H 5 substituent in the thiosemicarbazone moiety, the dihedral angles between the mean planes of A, and (S(8)-S(10)) fragments lie in interval 46.4 ÷ 52.0°. e "shifting" of nitrogen atom N(4) in the pyridine rings has led to different hydrogen bonds (HB) architecture in the crystal structures of HL, HL a , and HL b . In the crystal, the ligand HL forms the zigzag chains along b-axis, where the molecules are linked by N(2)-H. . .N(4) hydrogen bonds, while in HL a the similar HB forms two types of chains, which are oriented along different diagonals of (110) plane (Table S2, Figures 2(a), 2(b) and 2(f)). In HL b , there are three independent molecules (L1, L2, L3) in the asymmetric unit cell. L2 forms layers parallel to (011) plane due to N2A-H. . .S1A and N3A-H. . .N4 hydrogen bonds, while L1 and L3 are joined into the chains along b-axis through the N3-H· · · N4B and N3B-H· · · N4 HB. e chains and layers are further alternated along a-axis (Table S2, Figures 2(c), 2(g) and 2(h)).

Study of Substances in Solid
In complexes 4 and 5, the ligand HL acts as mononegative and neutral tridentate around the metallic ions, respectively, through an SNN set of donor atoms. Deprotonation of N2 atom in 4 has led to the decrease of N2-C1 bond distance if compared with that in HL. e bond lengths S1-C1 in 4 and 5 are increased due to the coordination of suphur atoms to central metals and the maximal changing of this bond is observed in 4 which is equal to 1.734(5) ( Table 2). However, the composition of the coordination polyhedron of the central atom in these compounds is different. In the crystal, the copper complexes form the centrosymmetric dimers where the monomers are hold together by bridge acetate molecules. e central atom in 4 is five-coordinated in a distorted square-pyramidal coordination geometry. Its basal plane includes three donor atoms of the L and oxygen atom O(1) of coordinated acetate molecule (Figure 2(d), Table 2). e deviations of these atoms from their mean plane are within 0.017Å, while the Cu atom deviates from this plane by 0.1Å toward the apex of pyramid. is apex of the metal's coordination pyramid in 4 is occupied by oxygen atom O(1) of adjacent acetate molecule with a distance of 2.429(3)Å. e polyhedral volume is related to the distortion of coordination polyhedral and its value for the copper complex is 6.765Å 3 . e structural characterization of 5 revealed that the nickel atom in this complex is in a distorted octahedral environment, being surrounded by two HL ligands (Figure 1(e), Table 2). e octahedral volume of Ni complex is equal to 13.338Å 3 . In the crystal of 4, there are water molecules in the outer coordination sphere which are joined by HB. In turn, these water molecules bound the neutral centrosymmetric dimers forming the 3D hydrogen bonding networks where dimers are linked also by C2-H· · · O2 HB (Table S2, Figure 2(d)). In the crystal of 5, the outer sphere contains both water molecules and negative chlorine ions. ese chlorine ions form the hydrogen bonds with positively charged octahedral nickel complexes (Table S2, Figure 2(e)) via bifurcated N2-H· · · Cl and N3-H· · · Cl bonds.

Infrared Spectra.
e process of formation of complexes 1-6 leads to the appearance of a number of changes in the IR spectra in comparison with the IR spectrum of free ligand HL. A change in the position of the ](S � N pyr ), ](C � N 1 ) bands in the IR spectra of 1-6 as well as the disappearance of the ](S � S) band in the spectra of 1-4 and 6 is observed. In case of complex 5, the ](S � S) does not disappear, but is only shifted by 49 cm −1 . e ](S � N pyr ) band is shifted by 17-45 cm −1 toward lower wavenumbers comparing with the free ligand HL (1583 cm −1 ). e ](C � N 1 ) band is shifted by 14-39 cm −1 . One ](NH) band disappears in the IR spectra of complexes 1-4 and 6. All the above changes indicate [47] that the azomethine and pyridine nitrogen atoms as well as the sulphur atom of the thiosemicarbazone HL are involved in the coordination process. While in the complexes 1-4 and 6, the thiosemicarbazone is in the thiol deprotonated form and in the composition of complex 5 thiosemicarbazone remains in the non-deprotonated thione form.

Magnetic Study.
A magnetochemical study was carried out in order to determine the presence of interaction between copper atoms in complexes 2-4. While the effective moment of complex 4 (1.25 μ B ) is lower than the spin only value (1.73 μ B ), which indicates the presence of interaction between copper atoms that is consistent with the determined structure by X-ray diffraction analysis, the effective magnetic moments of complexes 2 and 3 exceed the spin only value (1.85 μ B and 1.81 μ B , respectively). ese values indicate that there is no interaction between copper atoms in complexes 2 and 3, so that they should possess a monomeric structure [48]. e magnetic moment of complex 6 was measured in order to determine the oxidation state of cobalt. e complex 6 turned out to be diamagnetic, which is possible only if cobalt atom is in the +3 oxidation state in an octahedral ligand environment [49]. So, cobalt(II) undergoes oxidation by atmospheric oxygen during the synthesis process. e magnetic moment value of nickel complex 5 (2.84 μ B ) corresponds to two unpaired electrons and octahedral ligand environment that fully corresponds to the structure of the given substance determined by X-ray diffraction analysis. e obtained values for copper(II) complexes in methanol solution are rather low but rapidly change to normal values for 1 : 1 electrolytes upon adding of DMSO in the methanol solution. It indicates that methanol does not cause a complete substitution of anions from the inner sphere with the solvent molecule, while DMSO molecules shift this equilibrium toward the removal of anions from the inner sphere. Nickel complex 5 represents 1 : 2 electrolytes (χ �158 Ω −1 cm 2 mol −1 ). e corresponding spectrum of zinc(II) complex (1) contains two maxima at 325 and 397 nm. e spectra of complexes 2-4 contain two maxima at 300 and 409 nm and a shoulder at ca 335 nm. e corresponding spectrum of nickel(II) complex contains two maxima at 320 and 419 nm, while the spectrum of cobalt(III) complex contains two maxima at 310 and 371 nm and one shoulder at 416 nm. us, by comparing the results of the UV-Vis spectra of the complexes and the initial thiosemicarbazone in the 0.1 mM DMSO solutions, it can be concluded that thiosemicarbazpone moiety remains coordinated to the corresponding 3d metal ions when dissolved in DMSO [51]. It is interesting that the copper complexes have almost identical spectra. is may indicate that in the process of dissolution in DMSO the anion of the acid residue regardless of its nature is replaced by a DMSO molecule and a complex cation containing a copper atom and an Lligand is formed [52].

UV-Vis
is is consistent with the change in the electrical conductivity of methanol solutions of copper complexes that occur upon the addition of DMSO.

EPR Study.
e EPR signatures of complexes 2-4 were recorded in a 9-GHz band and the spectra obtained at 0.2 mM in DMSO are shown in Figure 3. ey are all very similar with a classical pattern for Cu(II) complexes lying in a square-based environment. e g//factor equals 2.20, while the hyperfine splitting equals 185 ± 5 10 −4 cm −1 (dotted lines). e parameters are in line with a equatorial donor set of [2N,S,O] [53] and strongly reminiscent with those of a Cu(II) complex bound to three nitrogen atoms and one thiolate ligand [54]. In the perpendicular region, superhyperfine lines are observed that mirror the interaction of the Cu(II) center with nitrogen atoms form the ligand. While the spectra of 3 and 4 are virtually superimposable and independent of the recording concentration, that of 2 shows broader hyperfine lines that witness the presence of a second minor species, the proportion of which increases with the concentration ( Figure S8). ese observations indicate, in line with the UV-Vis data, that the exogenous ligand is freed in solution and replaced by a solvent molecule with the exception of Cl − in 2 that partly stays bound.

Antiproliferative Activity and In Vivo Toxicity.
e antiproliferative activity of thiosemicarbazones HL, HL a , and HL b has been tested toward human leukemia HL-60 cell line. As it can be seen from Table 3, the position of nitrogen atom in the pyridine moiety of studied thiosemicarbazones has a tremendous influence on their antiproliferative activity. While HL inhibits the growth and proliferation of HL-60 cells by ca 100% at 10 and 1 μM concentrations, and HL a and HL b practically do not affect their proliferation. Because changing of the position of nitrogen atom in pyridine moiety leads to a complete disappearance of antiproliferative properties, only 2-formylpyridine N 4 -allylthiosemicarbazone was used for the synthesis of complexes. Subsequently, antiproliferative activity of the thiosemicarbazone HL and  (Table 4). e nature of the central atom has a major influence on the antiproliferative activity of the corresponding complex. e activity of zinc(II) complex (1) toward HL-60 cells is practically equal to the activity of the free ligand HL. Meanwhile, the activity of cobalt complex (6) is significantly lower, but copper (2-4) and nickel (5) complexes are much more active compared to the activity of HL. Complexes 2-5 inhibit the proliferation of HL-60 cells by ca 100% at 10 and 1 μM concentration and also manifest inhibitory activity (24-55%) at 0.1 μM concentration.
It is known from the literature sources that copper(II) complexes in many cases manifest the best antiproliferative activity. In this series of complexes, the nickel complex 5 showed the highest inhibitory activity toward HL-60 cells.
e copper complexes 2-4 manifest the highest activity toward the studied series of cancer cell lines. Moreover, these complexes are superior in antiproliferative activity to doxorubicin that is used in medical practice. Also, zinc complex 1 manifests high activity toward RD cells which is approximately equal to the activity of the most active copper complex 4.
Antiproliferative activity of these compounds toward normal MDCK cell line has been studied in order to determine the selectivity of their anticancer activity. Nickel (5) and cobalt (6) complexes practically do not affect the growth and proliferation of this normal cell line. Meanwhile, complex 6 does not possess promising anticancer activity, and nickel complex 5 manifests the best activity toward human leukemia cells and its activity is highly selective with selectivity index (SI) more than 1000. Also, free ligand HL and zinc complex 1 have high values of SI toward HL-60 that significantly exceed the selectivity of doxorubicin used in medical practice.
As the obtained results of anticancer activity of the studied series of complexes and their selectivity are of both theoretical and practical interest, their toxicity has been studied in vivo on Daphnia magna (Table 5). e percentage of viability has been measured after 24 h of staying in the environment with different concentrations of studied substances. As it can be seen from Table 5

Antibacterial and Antifungal Activity.
Despite the fact that the role of the virus and bacteria in the development of cancer is underestimated, antimicrobial prophylaxis during cancer therapy is important for reducing some side effects and the risk of death due to infection [55,56].
It is known from the literature that thiosemicarbazones and biometal complexes with these ligands in many cases also exhibit antimicrobial activity. erefore, the antimicrobial and antifungal properties of the synthesized compounds have been studied. ree strains have been selected for this study: S. aureus (Gram-positive bacterial strain), E. coli (Gram-negative bacterial strain), and C. albicans (fungal strain). e obtained results are presented in Table 6 in the form of minimum inhibitory, bactericidal, and fungicidal concentrations. Table 6 also shows data on nystatin and furacillinum [11] to compare the activity of the synthesized compounds.
As it can be seen from the obtained data, zinc and copper complexes manifest the highest antimicrobial activity. While their activity toward Gram-positive bacterial strain is only 2-4 times higher than the activity of thiosemicarbazone HL, the activity toward Gram-negative bacterial strain already exceeds 4-16 times the corresponding activity of free ligand. Copper complexes 2-4 also manifest promising antifungal activity that exceeds the activity of nystatin. e nickel (5) and cobalt (6) complexes manifest much lower antimicrobial and antifungal activities comparing with free ligand HL, complexes 1-4, and standard substances.

Electronic Properties of Studied Compounds.
e electronic structures and maps of molecular electrostatic potential (MEP) of HL, HL a , HL b , 4, 5 compounds have been calculated by Density Functional eory (DFT) of Gaussian 16 suite of quantum chemical codes based on B3LYP level of theory [59]. e 6-31G and LANL2DZ basis sets were used for ligands and metal atoms, respectively. e spin-polarized calculations were performed for compound 4. e molecular electrostatic potential was generated through a constant value of electron density. e values of MEP were calculated with the help of Multiwfn software [60] using the wave functions generated by the Gaussian.

Molecular Electrostatic Potential.
e molecular electrostatic potential is used to study and predict the molecular reactive centers. MEP is related to the electron density,         which is descriptor for identifying nucleophilic and electrophilic attack sites as well as intermolecular hydrogen bonding [61]. e electrostatic potential is also applied to analyze the processes based on the ''recognition'' of one molecule by another, as in drug-receptor, and enzyme-substrate interactions [62,63]. e molecular electrostatic potential surfaces of all studied molecules were mapped, using crystallographically determined geometries, and presented in Figure 4 e MEPs increase in the order −31.7 kcalˑmol −1 � red < yellow < green < blue � 31.7 kcalˑmol −1 . e negative (red and yellow) regions of MEP are related to electrophilic reactivity and the positive (blue) regions to nucleophilic reactivity. In HL, HL a , and HL b , the negative electrostatic potential is located at the regions of nitrogen atoms of pyridine rings and in the vicinity of sulphur atoms, while the positive MEP covers the outer tips of the CH-and NH-groups in the thiosemicarbazone moieties. However, the values of MEP of N4 depend on its position in the pyridine rings, and the   In compound 4, the negative MEP covers the metallacycle (Cu1S1N1N2C1) and oxygen atoms of acetate molecule, and the positive region is over the methyl group. In compound 5, it is observed only the strong positive electrostatic potential surrounding over whole molecule due to the positive charge of this complex.

Analysis of Frontier Orbitals.
It is known that the structure-activity relationship models for understanding the biological activity and toxicity of molecules may be generated in terms of global and local reactivity descriptors within a conceptual density functional theory framework [64]. e global reactivity indexes of studied compounds were estimated based on the one-electron energies of the frontier molecular orbitals (FMO) [65,66]. e calculated molecular descriptors and the relationships between them for studied compounds as well as biological activities are summarized in Table 8. e high occupied molecular orbitals (HOMO) and lower unoccupied molecular orbitals (LUMO) are referred to as frontier molecular orbital and act as an electron donor and acceptor, respectively. e energy gap between the HOMO and LUMO is related to kinetic stability as well as the chemical reactivity of the compound, the molecules having large/small values of energy gap are known as hard/soft molecules, as well as HOMO-LUMO band gap supports the bioactive property of the molecule [67]. e contour plots of the ground state FMO are shown in Figure 5. e HOMOs of HL, HL a , and HL b are located mainly on the sulphur atoms with a small contribution of nitrogen atoms of the chains, while the LUMOs are distributed in thiosemicarbazone moieties and pyridine rings. In compound 4, the spin np HOMO is localized on the central atom and acetate molecule, and the spin down HOMO is located mainly on metallocycles with some contribution of oxygen atoms. e spin np LUMO is distributed on acetate molecule, metallocycles, and pyridine ring, whereas the spin down LUMO is located mainly on acetate molecule and copper atom. It is interesting to observe that both HOMO and LUMO of compound 5 are substantially distributed over the conjugation planes of compound 5 which include metallocycles and pyridine rings. Compound 5 has the lowest energetic gap which allows it to be the softest molecule among the studied ones. Compounds 4 and 5 with high-lying E HOMO and low-lying E LUMO are the best electron donor and the electron acceptor, respectively, in a given series of compounds. e absolute electronegativity (χ) and the absolute hardness (η) are related to the one-electron orbital energies of the HOMO and LUMO. It is known that the stability of molecules is related to hardness and electronegativity is the power of an atom in a molecule to attract electrons to itself. e chemical hardness value of compound 5 is lesser among all the molecules, and it is found to be more reactive than all the studied compounds. is compound possesses higher electronegativity value than all compounds, so it is the best electron acceptor. e electrophilicity reactivity index (ɷ) measures the stabilization in energy when the system becomes saturated by electrons coming from environment. A good electrophile is characterized by a high value of µ (chemical potential) and a low value of η (chemical hardness). e higher electrophilicity of 5 is displaying that this compound would be more susceptible for nucleophilic attack because the electrophilicity is generally proportional to electron-poorness. e analyses of dependencies between the global reactivity indexes of studied compounds and their biological activities revealed that the experimental data were related to the theoretical ones (Table 8). It was found that the energetic gaps and hardness correlate with antiproliferative activities. e values of ionization potential, electron affinity, electronegativity, and electrophilicity are related to antibacterial and antifungal activities. e energies of LUMO and HOMO as well as chemical potentials correlate with antiradical activity. However, the most significant values of molecular descriptors and corresponding correlations with biological activities were observed for compound 5.

Conclusion
ree thiosemicarbazones (2-, 3-, and 4-formylpyridine N 4allylthiosemicarbazones) and six complexes of Zn(II), Cu(II), Ni(II), and Co(III) have been prepared and characterized using different chemical, physical, and physicochemical methods. e structures of HL, HL a , HL b and complexes 4 and 5 have been determined using single-crystal X-ray diffraction analysis. e HL ligand remains in a nondeprotonated form in the composition of nickel(II) complex and deprotonates in the case of compounds 1-4 and 6. It was found that the position of nitrogen atom N(4) in the pyridine rings affects the different hydrogen bonds architecture in the crystal structures of HL, HL a , and HL b . e copper atom in 4 is five-coordinated in a distorted square-pyramidal coordination geometry and in the crystal this compound forms the centrosymmetric dimers where the monomers are hold together by bridge acetate molecules. In 5, the nickel atom is in a distorted octahedral environment being surrounded by two HL ligands. e anticancer activity of the synthesized zinc, copper, and nickel compounds (1-5) is of interest. ese complexes show high activity against HL-60, HeLa, BxPC-3, and RD cancer cell lines. Nickel complex 5 manifests high antiproliferative activity toward HL-60 cells and practically does not affect normal cells showing the highest selectivity toward this cell line exceeding the selectivity of doxorubicin used in medical practice, and also the activity of other substances of thiosemicarbazone class such as triapine [68] that undergoes preclinical and clinical testing as anticancer drug. It also shows the lowest in vivo toxicity. It makes substance 5 promising for preclinical and clinical trials as anticancer substance. e copper and zinc complexes (1-4) also manifest high selective activity toward cancer cells toward a wide range of studied cancer cells, but they have a higher inhibitory activity toward normal cell line and are more toxic.
In addition to high anticancer properties, zinc and copper complexes also exhibit high antimicrobial properties with MIC/MBC/MFC values in the range of 0.12-250 μg/mL). e molecular electrostatic potential and the global reactivity indexes were calculated at B3LYP level of theory and then discussed in terms of their biological activities. It was found that the different reactivity descriptors of studied compounds strongly correlate with their biological activity and such dependences provide information for structureactivity study to design the biologically important molecules that can be applied in the pharmacy and medicine.
Data Availability e data that supports the findings of this study are available in the supplementary material of this article.

Supplementary Materials
CCDC 2129910-2129914 contains the supplementary crystallographic data for HL, HL a , HL b , and complexes 4 and 5.