Theoretical Investigation by DFT and Molecular Docking of Synthesized Oxidovanadium(IV)-Based Imidazole Drug Complexes as Promising Anticancer Agents

Vanadium compounds have been set in various fields as anticancer, anti-diabetic, anti-parasitic, anti-viral, and anti-bacterial agents. This study reports the synthesis and structural characterization of oxidovanadium(IV)-based imidazole drug complexes by the elemental analyzer, molar conductance, magnetic moment, spectroscopic techniques, as well as thermal analysis. The obtained geometries were studied theoretically using density functional theory (DFT) under the B3LYP level. The DNA-binding nature of the ligands and their synthesized complexes has been studied by the electronic absorption titrations method. The biological studies were carried with in-vivo assays and the molecular docking method. The EPR spectra asserted the geometry around the vanadium center to be a square pyramid for metal complexes. The geometries have been confirmed using DFT under the B3LYP level. Moreover, the quantum parameters proposed promising bioactivity of the oxidovanadium(IV) complexes. The results of the DNA-binding revealed that the investigated complexes bind to DNA via non-covalent mode, and the intrinsic binding constant (Kb) value for the [VO(SO4)(MNZ)2] H2O complex was promising, which was 2.0 × 106 M−1. Additionally, the cytotoxic activity of the synthesized complexes exhibited good inhibition toward both hepatocellular carcinoma (HepG-2) and human breast cancer (HCF-7) cell lines. The results of molecular docking displayed good correlations with experimental cytotoxicity findings. Therefore, these findings suggest that our synthesized complexes can be introduced as effective anticancer agents.


Introduction
Vanadium (V) is one of the transition metals, with atomic number 23 and composed of two isotopes, 50 V and 51 V with atomic abundance of 0.25% and 99.75%, respectively. Vanadium is disseminated across the whole planet in small amounts and is present in the human body in trace amounts [1,2]. The ability of vanadium to influence numerous metabolic processes in the cell is due to possession of the same character traits as phosphate, which plays a vital role in medicinal applications [3]. Therefore, vanadium compounds have been determined in various fields as anticancer, anti-diabetic, anti-parasitic, anti-viral, anti-bacterial, anti-thrombotic, anti-hypertensive, anti-atherosclerotic, spermicidal, anti-human deficiency virus, and anti-tuberculous agents [4,5]. Although there are various drugs for cancer treatment, it remains the second most common cause of death worldwide [6,7]. The limitations of treatment include several side effects, the development of multidrug resistance, and high toxicity, which lead to a critical need to discover novel anticancer compounds with potent anticancer activity, less toxic than the platinum-based drugs, and less side effects [8]. The historical use of oxidovanadium(IV) complexes in various treatments is promising by the addition of new complexes that may be effective against cancer [9,10]. Thus, several studies have reported the cytostatic activity and ability to defeat tumor cell growth in vitro and in vivo of the inorganic compounds of vanadium (valence varied from III to V) coordinated with organic ligands [11].
Imidazole is a five-member ring molecule that has a wide range of biological activities and is an essential pharmacophore in drug discovery [12]. Several imidazole derivatives have been used as effective anticancer agents due to the remarkable pharmaceutical potencies to control cancer [13,14]. One of the significant properties of imidazole is easy binding to DNA or protein molecules via different interactions, such as electrostatic, intercalation, and groove binding [15,16].
In the search, we merged their benefits to obtain novel vanadyl-based drug complexes with a valuable anticancer potency. Three imidazole molecules, named Clotrimazole (CTNZ), Miconazole (MNZ), and Pantoprazole (PNZ), were reacted with oxidovanadium(IV) salt in ethanol solution and then characterized by various approaches. The molecular structures of all new complexes were optimized by theory level DFT/B3LYP to investigate the geometry properties and the quantum parameters. The biological efficiency of the metal complexes was tested by a DNA-binding study and molecular docking. Moreover, an in vitro cytotoxicity assay was carried out against two cell lines, hepatocellular carcinoma (HepG-2) and breast adenocarcinoma (MCF-7).

Materials and Instrumentation
The reagents and chemicals utilized in this research were purchased from various commercial sources with 100% purity (Fisher Scientific, CAD Middle East; Aldrich, Riyadh, Saudi Arabia) and used without purification. Melting points were measured by MPA100, an automated melting point system. The conductivity measurements were determined using a conductor TDA meter at room temperature in 10 −3 M DMSO solution. IR spectra were recorded using solid samples on a PerkinElmer Spectrum 100 instrument (Waltham, MA, USA), and peaks are reported over the range 4000-400 cm −1 . The electronic absorption spectra were determined in DMSO from the MultiSpec-1501 UV-Vis spectrophotometer. The EPR spectra for the complexes in solution were recorded at room temperature by a continuous-wave Bruker EMX PLUS spectrometer (Bruker BioSpin, Rheinstetten, Germany) and collected with Bruker Xenon software with DPPH standard (υ = 9.440103 GHz). TGA was obtained from the PerkinElmer TGA system, in the range of 25-900 • C and at a 10 • C min −1 heating rate under nitrogen atmosphere.

Synthesis of VO(II) Complexes
Vanadyl sulfate (VOSO 4 .5H 2 O) was dissolved in 3 mL of distilled H 2 O and added dropwise to 25 mL of ethanolic solution of ligand and was refluxed for 2-3 h. The produced precipitates were filtered, collected, washed by hot ethanol, drops of diethyl ether were added, and they were dried in an oven at 50 • C.

Computational Details
Theoretical calculations of the investigated compounds were carried out by the Gaussian 09 program package [19], employing the density function theory method (DFT), and the hybrid density functional B3LYP in the gas phase. The basis set 6-31G was utilized for the free drug molecules while the metal complexes were performed with the correlationconsistent basis set LANL2DZ to accommodate the large size of the target complexes. Additionally, the MEP map was projected by applying the B3LYP/LANL2DZ level within DFT. The vibrational frequencies calculation was obtained by visualizing output files via Gauss View 6.0. All stationary points for optimized molecules are minima on the potential energy surface and do not have imaginary frequencies, and were scaled by factor values of 0.99 and 0.96 to correct for anharmonicity and the neglected aspect of electron correlation.

DNA Interaction Studies
The interaction of oxidovanadium(IV) complexes with Calf Thymus-DNA (CT-DNA, purchased from Sigma Aldrich) was examined in Tris-HCl buffer solution at a pH = 7.4 using UV-Vis absorption spectroscopy. The CT-DNA stock was stored in a refrigerator and used within one week. The solution of CT-DNA was confirmed to be free from proteins by checking the ratio of UV absorbance at 260 and 280 nm (A 260 /A 280 ), which was in the 1.8-1.9 range [20,21]. The CT-DNA concentration was determined from the absorbance calculated at λ max 260, with a molar extinction coefficient of 6600 M −1 cm −1 . Stock solutions of oxidovanadium(IV) complexes were dissolved in DMSO (1 × 10 −4 M) and suitably diluted with the corresponding buffer solution to the required concentrations. The intrinsic binding constant (K b ) of the complexes was calculated from the plot of [DNA]/(εa − εf) vs.
[DNA] using the following equation, where K b is given by the ratio of the slope to the intercept: [DNA] is the concentration of CT-DNA, εf is the extinction coefficient of the free compound, εa is equal to A obsd /[Compound], and εb is the extinction coefficient of the complexes.

Docking Studies
Molecular docking investigations were performed using the Molecular Operating Environment (MOE) program. The crystal structures of breast cancer (PDB code 1HK7) binding to this protein present the inhibition of hsp90 alpha, which leads to the degradation of client proteins involved in the initiation of breast cancer pathogenesis, and human thymidylate synthase for hepatocellular carcinoma (PDB code 1JU6) proteins were retrieved from the protein data bank (https://www.rcsb.org/ accessed on 30 March 2021). All the water molecules around the proteins were removed, and hydrogen atoms were added. Then, the optimized geometry of each tested compound with the minimum binding energy was attained using the MMFF94x force-field [22,23]. The site finder in the MOE module was used to generate alpha-site spheres. The DFT-optimized structures of the compounds were used to generate the five best binding configurations with flexible molecular rotation. The hydrogen bonds that formed between proteins and the investigated compounds were used to rank the binding affinity, which was expressed as the free binding energy (S, kcal/mol).
The validation of our docking method was performed by docking the co-crystalline ligand with the same protocol and comparing the interaction between it and the downloaded one, and we obtained RMSD in the range of 1.37-1.52 Å.

Anticancer Studies
Oxidovanadium(IV) complexes were evaluated for their anticancer activity against hepatocellular carcinoma (HepG-2) and breast adenocarcinoma (MCF-7). In this assay, a 96-well plate with 100 pl/well was used to plate the cell culture (1 × 10 4 cells/mL). The tested complex was diluted into each well with calculated concentrations (3.125, 6.25, 12.50, 25.00, 50.00, and 100 ug/mL) and incubated. After 24 h of incubation, MTT ((3-14,5-dimethylthiazol-2-yl)-2,5-diphenylte trazolium bromide) liquid was added and the plate incubated for a further 24 h. Each well's optical density (OD) was measured spectrophotometrically, and the 50% concentration causing growth inhibitory action (IC 50 value) to tumor cells was estimated using the following formula: Cell viability % = OD value of experimental sample (mean) OD value of experimental control (mean) × 100 The acute oral toxicity study was performed as documented in [24]. Briefly, healthy BALB/C mice (

Stoichiometry and Physical Properties of the Prepared Complexes
The maximum absorbance wavelengths were measured relative to a blank sample.  Figure S1. The reflection line appeared at 0.5 upon increasing the ligand concentration for complex VO-CTNZ, which revealed that the metal:ligand combination ratio was 1:1. In contrast, the reflection line for complexes VO-MNZ and VO-PNZ was around 0.33, indicating a 1:2 M:L ratio.
The physical properties and the elemental analysis for the ligands and their complexes are listed in Table 1. All complexes are insoluble in water and most traditional organic solvents, but soluble in dimethylformamide (DMF) and dimethylsulfoxide (DMSO). Molar conductivity measurements of the obtained complexes were performed at 25 • C in DMSO at a concentration of 10 −3 M. The molar conductivity values for complexes VO-CTNZ and VO-MNZ ranged between 11.72 and 3.00 Ω −1 cm 2 mol −1 , respectively, thus revealing their non-electrolytic nature. Complex VO-PNZ afforded higher values, indicating the presence of sulfate outside the coordination sphere. The suggested structures of the complexes were established by CHN analysis, various spectroscopic methods, and thermal analysis. They displayed a square pyramidal geometry and were stable toward air and moisture ( Figure 1). organic solvents, but soluble in dimethylformamide (DMF) and dimethylsulfoxide (DMSO). Molar conductivity measurements of the obtained complexes were performed at 25 °C in DMSO at a concentration of 10 −3 M. The molar conductivity values for complexes VO-CTNZ and VO-MNZ ranged between 11.72 and 3.00 Ω −1 cm 2 mol −1 , respectively, thus revealing their non-electrolytic nature. Complex VO-PNZ afforded higher values, indicating the presence of sulfate outside the coordination sphere. The suggested structures of the complexes were established by CHN analysis, various spectroscopic methods, and thermal analysis. They displayed a square pyramidal geometry and were stable toward air and moisture ( Figure 1

Optimized Molecular Geometry
The results from the spectroscopic analysis (IR, UV-Vis, and EPR) were used to generate the input file to proceed with a computational study and to optimize the complexes' geometry in gas phase using the DFT method. The extracted optimized geometries from the log file are presented with the number and label atoms in Table S2. The comparison between the free ligand and the formed complexes showed elongation or compression after the reaction with metal ions. The bond lengths around the metal coordination for PNZ, C18-S19, C20-S19, and S19-O31 were longer than the free ligand, except C18-S19 was shorter with 0.04 Å A. The bond angles also showed a practical change after complexation (Table 2), where the bond angles decreased by 0.009-0.038° due to the electric repulsion after complexing. On the other hand, the changes in the MNZ and CTNZ were negligible after the reaction with oxidovanadium(IV) for both bond lengths and angles.

Optimized Molecular Geometry
The results from the spectroscopic analysis (IR, UV-Vis, and EPR) were used to generate the input file to proceed with a computational study and to optimize the complexes' geometry in gas phase using the DFT method. The extracted optimized geometries from the log file are presented with the number and label atoms in Table S2. The comparison between the free ligand and the formed complexes showed elongation or compression after the reaction with metal ions. The bond lengths around the metal coordination for PNZ, C18-S19, C20-S19, and S19-O31 were longer than the free ligand, except C18-S19 was shorter with 0.04 Å A. The bond angles also showed a practical change after complexation (Table 2), where the bond angles decreased by 0.009-0.038 • due to the electric repulsion after complexing. On the other hand, the changes in the MNZ and CTNZ were negligible after the reaction with oxidovanadium(IV) for both bond lengths and angles.
The new bonds between oxidovanadium(IV) and O and N sites in all investigated compounds had average bond lengths of 1.09-2.09, which are assumed to have an ionic character of the metal to oxygen and nitrogen, respectively [25], excluding the bond length of V-O45 in the CTNZ, which showed a covalent bond character with bond length >2 Å. Moreover, the angles around the individual central vanadium ion indicated that the compounds implemented a pyramidal geometry, as presented in Table S2, Supplementary Materials. The dihedral angles around the metal ions were between 175 • and −175 • , which are away from 0 • or 180 • , indicating that the ligand molecules and metal ions are not on the same plane.

Assignments of Vibrational Spectra
Vibration frequency analysis is an essential tool to investigate the binding coordination mode for metal complexes through peak shift observation [26]. The IR spectra were measured between 4000 and 400 cm −1 for the ligands and their oxidovanadium(IV) complexes. Moreover, the FT-IR analysis for the free and metal compounds was performed with DFT/B3LYP in the same range. The essential functional groups and their assignments were compared with theoretical calculation and listed in Table 3. The IR spectra of the ligands showed a characteristic peak between 1590 and 1494 cm -1 for the imidazole ν(C=N) group, which was shifted to higher frequencies in all complexes in the range of 1622-1593 cm -1 , confirming the coordination of imidazole nitrogen to the central vanadium ion. Further, the IR spectrum of free PNZ showed a band at 1043 cm −1 assigned to ν(S=O) of the side chain. This band appearing at 1027 cm −1 in VO-PNZ indicated the coordination of the O atom to the vanadium ion [27]. The broad peaks in the region from 3445 to 3315 cm −1 are due to the stretching vibrations of water molecules, ν(OH 2 ). Furthermore, new bands observed in the spectrum of complex VO-CTNZ at 874 and 799 cm −1 which might be assigned to δ r (H 2 O) at 874 cm −1 and δ w (H 2 O) at 799 cm −1 affirmed coordination with the water molecule [28]. The bands at 1091 and 1070 cm −1 in the spectrum of the free MNZ and PNZ drugs, respectively, corresponding to the stretching vibrations of the C-O-C group have shifted in the obtained coordination compounds to 1090 cm −1 for complex VO-MNZ and to 1068 cm −1 for complex VO-PNZ. Additionally, the bands relating to ν(SO 4 ) in the range 1289-1040 cm −1 support its bidentate coordination [29]. In all VO(II) complexes, a significant ν(V=O) band was recorded around ≈970 cm −1 . There were new bands that appeared in the complexes' spectra in the regions 588-525 and 451-438 cm −1 , which were attributed to the formation of ν(V-O) and ν(V-N) bonds, respectively. The calculated frequencies were well-connected with the experimental results.
The metal-oxygen bond length (V=O) was determined from the stretching frequencies and calculated by the below equation [30,31]:  Table 3, a real correlation between the experimental results with the obtained theoretical values can be seen, 1.60-1.59 Å. Figure 2 shows the comparison of the theoretical and experimental values of the IR spectra assignments of the free ligands and the synthesized metal complexes.

Frontier Molecular Orbitals (FMOs) Parameters
A practical method to describe the electronic characteristics of molecular arrangements and study the orbital compositions is the frontier molecular orbitals (FMO) [32]. These orbitals were generated from the fchk files and visualization for both the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) to study the molecular orbital compositions and energy levels of a molecule, and are presented in Figure 3. The ground state, HOMO, for the PNZ was mainly distributed over most of the ligand molecule, while it was concentrated in one of the dimethoxypyridin-2-yl methanethiols in the metal complex, while the LUMO was at the other one in the ligand molecule of the ratio 2:1.

Frontier Molecular Orbitals (FMOs) Parameters
A practical method to describe the electronic characteristics of molecular arrangements and study the orbital compositions is the frontier molecular orbitals (FMO) [32]. These orbitals were generated from the fchk files and visualization for both the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) to study the molecular orbital compositions and energy levels of a molecule, and are presented in Figure 3. The ground state, HOMO, for the PNZ was mainly distributed over most of the ligand molecule, while it was concentrated in one of the dimethoxypyridin-2yl methanethiols in the metal complex, while the LUMO was at the other one in the ligand molecule of the ratio 2:1. In contrast, MNZ showed a distribution of HOMO and LUMO over one of the coordination ligand molecules. The HOMO for CTNZ was located on the metal center and the donating atoms, which raise the ability of the metal complex to react as a donor or an accepter to react with the biological residue. The obtained values for the HOMO and In contrast, MNZ showed a distribution of HOMO and LUMO over one of the coordination ligand molecules. The HOMO for CTNZ was located on the metal center and the donating atoms, which raise the ability of the metal complex to react as a donor or an accepter to react with the biological residue. The obtained values for the HOMO and LUMO for the metal complexes and the free organic molecules have been used to calculate the energy gap, E g = E LUMO − E HOMO . This indicator is usually used to assess the reactivity of the molecules [33]. A small E g value indicates the reactivity of the complexes through to the softness character, while the large energy gap between HOMO and LUMO complicates the transition of the electron, and we can express this as hard compounds. Table 4 shows that the E g values were decreased after coordination to the metal ion, and CTNZ had the smallest value among the investigated compounds and was thus the most reactive molecule. For more details of the molecules' properties, we estimated the global reactivity descriptors, which are absolute electronegativities (χ), absolute hardness (η), chemical potentials (µ), global softness (S), and global electrophilicity (ω), using equations in our previous work [33] and illustrated in Table 4. The absolute softness is the reverse of the absolute hardness, which describes the sense to react with other molecules such as biological compounds, and we could arrange them according to the absolute softness (S) as follows: VO-CTNZ > VO-MNZ > VO-PNZ. Thus, VO-CTNZ complex had the maximum S value, indicating high anticipation of effective biological activity. The electrophilicity indicator (χ) has a positive value, but the electronic chemical potential (µ) has a negative value, which revealed the capability of the new complexes to capture electrons from the environment to decrease their energy and to stabilize the reaction. Therefore, the following order presents the most potent compound with a high electrophilicity indicator: VO-CTNZ > VO-MNZ > MNZ >VO-PNZ > CTNZ > PNZ.

Molecular Electrostatic Potential Maps (MEPs)
The calculated MEPs were obtained to provide the molecule electron distribution over the compounds, which provided insight into the electrophilic and nucleophilic attack sites. Molecular electrostatic potential surface maps of the metal complexes are presented in Figure 4. The maps have a red color which presents the electron-rich zone and the nucleophile site, while the blue areas represent the positive electrophile region. The VO-CTNZ and VO-MNZ showed that the SO 4 group had the highest electronic density in the molecules, while the rest of the molecules were mostly in the blue area and reacted as electrophile. Interestingly, the blue color was distributed over the whole of the VO-PNZ complex, which increased the tendency of the molecule toward electrophilic attack with the negative phosphate group in the biological residues.

Electronic Spectroscopic and Magnetic Susceptibility
Electronic spectral data of all the ligands and metal complexes are listed in Table 5 and the spectra are shown in Figure 5, measured in DMSO solution (10 − 3 to 10 −4 M) in the

Electronic Spectroscopic and Magnetic Susceptibility
Electronic spectral data of all the ligands and metal complexes are listed in Table 5 and the spectra are shown in Figure 5, measured in DMSO solution (10 −3 to 10 −4 M) in the range of 200 to 1000 nm at room temperature. The complexes contain bands in the ranges of 257-273 and 280-299 nm, which are ascribed to π-π* transitions that originate mainly from the aromatic rings and n→π* transitions that could be assigned to the C=N chromophore. These bands shifted to higher wavelengths after the reaction with metal ions. Electronic spectral data of all the ligands and metal complexes ar and the spectra are shown in Figure 5, measured in DMSO solution (10 range of 200 to 1000 nm at room temperature. The complexes contain ba of 257-273 and 280-299 nm, which are ascribed to π-π* transitions tha from the aromatic rings and n→π* transitions that could be assigned to phore. These bands shifted to higher wavelengths after the reaction wit  Low-intensity peaks observed at higher wavelengths around 800 nm in VO-CTNZ and VO-MNZ are due to low-energy d-d transitions [34]. However, the electronic spectrum of VO-PNZ shows a shoulder band at 317 nm and a low-intensity band at 472 nm that may be attributed to LMCT, and the spin-allowed d-d electronic transitions, respectively [34]. Moreover, the effective magnetic moment values (µ eff ) were calculated from the measured magnetic susceptibility using the following equation: where T is the absolute temperature and χ m is the molar susceptibility constant [35,36]. Complexes VO-CTNZ, VO-MNZ, and VO-PNZ had magnetic moment values of 1.75, 1.75, and 1.76 Bohr's magneton (BM), respectively. The obtained values indicated that the VO(II) complexes have normal paramagnetic properties, with one unpaired electron and square pyramidal geometry around the central metal ion [37,38].

EPR Spectroscopy
EPR spectra of the VO(II) complexes were recorded in DMF at ambient temperature, and the spectral parameters are summarized in Table 6. The spectra exhibit eight lines, which are due to hyperfine splitting arising from the interaction of the unpaired electron with a 51 V nucleus having the nuclear spin I = 7/2 ( Figure 6). This confirms the presence of a single cation as the metallic center in the complex. The measured values of g and A are in good agreement for a square pyramidal structure [39,40].  Figure 6. EPR spectra of oxidovanadium(IV) complexes.

Thermal Analysis
The thermogravimetric analysis (TGA) was measured from 25 to 900 °C to co

Thermal Analysis
The thermogravimetric analysis (TGA) was measured from 25 to 900 • C to confirm the stability of our complexes and the molecular structure (Figure 7), and the results were summarized in Table S3. The TG curve  [VO(SO4)(MNZ)2] H2O complex was thermally decomposed in one step within the temperature range from 25 to 800 °C, accompanied by a weight loss of 80.35% (calc. 80.98%), and corresponding to the loss of C27H28Cl8N4O5S and half-hydrated water molecules, forming VO2 and 9 carbon atoms as a final product with a mass loss of 19.65% (calc. 19.01%).
The decomposition of [VO(PNZ)2]SO4·2H2O occurred in two steps. The first step was within a temperature range of 25-109 °C and corresponds to the loss of two moles of hydrated water with a mass loss of 3.88% (calc. 3.73%). The second step was within a temperature range of 109-800 °C and equals the loss of the organic moiety, C32H30F4N6O11S3, with a weight loss of 87.59% (calc. 87.68%). The final product is VO2 residue, with a total mass loss of 8.52% (calc. 8.59%).

Kinetic Studies of the TGA
Recently, the Coats-Redfern (CR) [41] and Horowitz-Metzger (HZ) [42] methods have been used effectively around the world to depict thermal decomposition kinetics. The entropy (ΔS), enthalpy (∆H), and Gibbs free energy (ΔG) were calculated using previous equations [43]. According to the kinetic parameters presented in Table S4, Supple- [VO(SO 4 )(MNZ) 2 ] H 2 O complex was thermally decomposed in one step within the temperature range from 25 to 800 • C, accompanied by a weight loss of 80.35% (calc. 80.98%), and corresponding to the loss of C 27 H 28 Cl 8 N 4 O 5 S and half-hydrated water molecules, forming VO 2 and 9 carbon atoms as a final product with a mass loss of 19.65% (calc. 19.01%).
The decomposition of [VO(PNZ) 2 ]SO 4 ·2H 2 O occurred in two steps. The first step was within a temperature range of 25-109 • C and corresponds to the loss of two moles of hydrated water with a mass loss of 3.88% (calc. 3.73%). The second step was within a temperature range of 109-800 • C and equals the loss of the organic moiety, C 32 H 30 F 4 N 6 O 11 S 3 , with a weight loss of 87.59% (calc. 87.68%). The final product is VO 2 residue, with a total mass loss of 8.52% (calc. 8.59%).

Kinetic Studies of the TGA
Recently, the Coats-Redfern (CR) [41] and Horowitz-Metzger (HZ) [42] methods have been used effectively around the world to depict thermal decomposition kinetics. The entropy (∆S), enthalpy (∆H), and Gibbs free energy (∆G) were calculated using previous equations [43]. According to the kinetic parameters presented in Table S4, Supplementary Materials, for each step of decomposition, large positive values of E a indicate that the complexes are extremely stable [9,44]. The values of ∆S for the complexes were negative, which indicates that the thermal decomposition was slower than the formation reaction, or that there was a greater degree of balance within the activated complex than in the reactants [45]. The positive ∆H values suggest that all the decomposition steps were endothermic in nature. The high values of ∆G for each stage are due to the increase in T∆S values from one step to the next. Moreover, the correlation coefficient values demonstrate that the two methods used provide a good linear function for the thermal decomposition stages of the complexes. The Coats-Redfern and Horowitz-Metzger plots of all complexes are represented in Table S5, Supplementary Materials.

DNA-Binding Study
The electronic absorption spectra were evaluated for the ligands and their oxidovanadium(IV) complexes by monitoring the changes in the maximum absorption bands in the range 262-301 nm assigned to the π→π * transitions of the aromatic rings (Figures 8 and S3). The UV-Vis absorbance of the free ligands and VO-PNZ showed increasing molar absorptivity (hyperchromic) with increasing CT-DNA concentration. Hyperchromism occurs due to the breakage of the secondary structure of DNA after the interaction with compounds and DNA, which derogates the interaction forces maintaining the DNA double-helix structure. Thus, the solution had many bases in free helix form and did not form hydrogen bonds with complementary bases, causing heightened absorbance [46,47]. The obtained results clearly show a groove-binding mode of binding of the free ligands and VO-PNZ complex to DNA. On the other hand, the absorption bands of VO-CTNZ and VO-MNZ decreased (hypochromic) with the rising amount of CT-DNA. Hypochromism is an intercalative reaction between the aromatic part of the compound with the base-pair π electrons of the DNA, resulting in a decrease of the distance between this and the π electron in the aromatic ring [48]. On further addition, there was a notable blue shift with hypochromism, indicating that VO-CTNZ and VO-MNZ have potentially used intercalation to bond with the DNA. The determined binding constants (K b ) for the complexes were ordered as follows: VO-MNZ > VO-PNZ > VO-CTNZ. The calculated K b and absorption data of the complexes are presented in Table S6

DNA-Binding Study
The electronic absorption spectra were evaluated for the ligands and their oxidovanadium(IV) complexes by monitoring the changes in the maximum absorption bands in the range 262-301 nm assigned to the π→π * transitions of the aromatic rings (Figures 8 and  S3). The UV-Vis absorbance of the free ligands and VO-PNZ showed increasing molar absorptivity (hyperchromic) with increasing CT-DNA concentration. Hyperchromism occurs due to the breakage of the secondary structure of DNA after the interaction with compounds and DNA, which derogates the interaction forces maintaining the DNA double-helix structure. Thus, the solution had many bases in free helix form and did not form hydrogen bonds with complementary bases, causing heightened absorbance [46,47]. The obtained results clearly show a groove-binding mode of binding of the free ligands and VO-PNZ complex to DNA. On the other hand, the absorption bands of VO-CTNZ and VO-MNZ decreased (hypochromic) with the rising amount of CT-DNA. Hypochromism is an intercalative reaction between the aromatic part of the compound with the base-pair π electrons of the DNA, resulting in a decrease of the distance between this and the π electron in the aromatic ring [48]. On further addition, there was a notable blue shift with hypochromism, indicating that VO-CTNZ and VO-MNZ have potentially used intercalation to bond with the DNA. The determined binding constants (Kb) for the complexes were ordered as follows: VO-MNZ > VO-PNZ > VO-CTNZ. The calculated Kb and absorption data of the complexes are presented in Table S6, Supplementary Materials.

Molecular Docking
Molecular docking helps us to understand the possible modes of interaction and binding affinities of new compounds as targeted drug molecules, which may offer better insight into how well they exert bioactivity against the target [45]. The docked complex conformations were evaluated according to the docking scores and binding interactions of the ligands and oxidovanadium(IV) complexes with the hepatocellular carcinoma (ID: 1JU6) and breast cancer (ID: 1HK7) proteins, which were selected according to the literature and previous studies [49]. The docking study results indicated that the ligands and synthesized complexes involve non-covalent interactions, such as π-cation, π-π, ionic, and hydrogen bonding interactions, as presented in Table S7. All the studied compounds exhibited excellent docking score values with the proteins (1JU6 and 1HK7). In general, the more negative the energy value (S), the greater the stability and the stronger the binding between the docking sites and protein receptors. The binding free energy docking scores (S) and root-mean-square deviations (RMSD) of atomic positions are presented in Tables S8 and S9, Supplementary Materials. The estimated binding scores of the free ligands CTNZ, MNZ, and PNZ were calculated as −6.66, −8.09, and −8.59 kcal/mol, respectively, for the 1JU6 protein, while the binding score values of the synthesized complexes VO-CTNZ, VO-MNZ, and VO-PNZ were −7.23, −11.09, and −11.30 kcal/mol, respectively. The score values against the 1JU6 protein toward the free ligand were −6.81, −7.11, and −7.11 kcal/mol, respectively, and −9.85, −7.64, and −11.10 kcal/mol, respectively, for their metal complexes. Notably, the results revealed high docking scores for the complexes compared to those for the free ligands. This indicates that the modified metal complexes enhance the binding affinity of the free ligands. In terms of the coordination pattern, the VO-CTNZ complex formed a π-π bond and four hydrogen bonds with the following amino acid residues: oxygen 45 with CYS 195, oxygen 46 with HIS 256, oxygen 48 with ARG 50, oxygen 49 with SER 216, and the benzene ring with TRP 109, with bond distances of 2.80-3.53 Å. In the 1HK7 protein, the VO-CTNZ complex formed H-acceptor bonds via oxygen 44 and oxygen 49 atoms with TRP 109 and ARG 344 amino acid residues, with bond distances of 3.00-3.42 Å. Additionally, in the 1JU6 protein, the main contributions to the activity of the VO-MNZ complex were the interactions with the amino acids CYS 195, ARG 50, ASN 112, ASP 218, SER 216, and LEU through hydrogen and ionic bonds, with bond distances of 3.10-4.11 Å. Moreover, the imidazole ring displayed π-π bonding with the PHE 140 amino acid residue, and π-H bonding with the ARG 51 amino acid residue in the 1HK7 protein, with bond distances of 3.73-3.80 Å. The VO-PNZ complex displayed the best binding scores of the docked complexes against both proteins (−11.30 and −11.10 kcal/mol), forming four hydrogen bonds with ARG 78, MET 311, ILE 108, and PHE 225 amino acids for the 1JU6 protein, with bond distances of 3.39-4.40 Å. An additional H-acceptor bond was observed between VO-PNZ and TYR 337 in the 1HK7 protein, with a bond distance of 2.63 Å. All the 3D and 2D interactions for the two targets are illustrated in Tables 7 and S10.

Cytotoxicity Study
Our aim in this study was to evaluate the anticancer activity of newly synthesized oxidovanadium(IV) complexes against two different human tumor cell lines: HepG-2 and MCF-7 cells. We have chosen MCF-7 and HepG-2 cells because they are considered the most common types of carcinomas. The calculated IC 50 values are shown in Table 8 and a graphical statistical analysis is represented in Figure 9. Their performance was compared to an anticancer chemotherapy drug, cisplatin, as a standard reference. Initially, the results were consistent with the docking predictions, revealing that all the synthesized vanadyl complexes displayed excellent activities against these two cell lines. The IC 50 (µg/mL) values indicate the following order of toxicity towards the HepG-2 cell line: VO-PNZ > VO-CTNZ > VO-MN. According to the results obtained for the HCF-7 cell line, the order of toxicity is as follows: VO-PNZ > VO-MNZ > VO-CTNZ. The molar concentration (µM) showed that the VO-PNZ potency is double the inhibition of VO-CTNZ toward HepG-2. Despite that VO-PNZ has a good effect against the HCF-7 cell line, with 2.93 µg/mL, the obtained µM value of VO-MNZ was better, with 0.05 µM. Morphological observations of the HepG-2 and MCF-7 cell lines showed noteworthy changes in their morphology and decreased cell density, as shown in Figures 10, S4 [50], a chemical substance with an LD 50 within the range of 5000-15,000 mg kg −1 is regarded as practically nontoxic, slightly toxic, moderately toxic, and highly toxic, respectively. The results indicated that both complexes fall within the category indicating a slightly toxic substance. Based on the results, it may be concluded that the synthesized oxidovanadium(IV) complexes exhibit promising anticancer activity against both cell lines [51].

Conclusions
New oxidovanadium(IV) complexes were synthesized and characterized by various techniques. From the elemental analysis, electronic spectra, magnetic moment values, molar conductance, and EPR, the newly synthesized complexes were found to have a square pyramidal geometry around the central metal ion, with a molar ratio of the metal to the ligand of 1:1 for VO-CTNZ, and a 1:2 (M:L) ratio for VO-MNZ and VO-PNZ. The DFT study confirmed this geometry structure and reveled that the complexes were more reactive than the free ligand. The VO-CTNZ showed the highest softness among the prepared complexes, which could indicate favorable biological activity. The positive charge around the VO-PNZ presented by MPE indicated a practical tendency to react toward biological molecules. A kinetic study was conducted of the thermal decomposition using the Coats-Redfern and Horowitz-Metzger methods, suggesting that all complexes were non-spontaneous, and thereby that the complexes were thermally stable. DNA-binding features were investigated by using spectrophotometric titrations. The binding ability of the mentioned complexes followed the order: VO-MNZ > VO-PNZ > VO-CTNZ. Anticancer activity was tested for all compounds against hepatocellular carcinoma and human breast cancer cell lines, in which all complexes exhibited a significant inhibition, with the highest cytotoxicity value for the VO-PNZ. Additionally, the in vitro cytotoxicity of compounds was supported by docking studies.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/molecules27092796/s1, Table S1: Experimental data by molar ratio method; Figure S1: The molar ratio for VO(II) complexes; Table S2: The optimized geometry and numbering system for the free ligands and Oxidovanadium(IV) complexes; Table S4: Kinetic parameters of thermal decomposition steps for VO(II) complexes; Table S5: Coats-Redfern (CR) and Horowitz-Metzger (HM) of the VO(II) complexes; Figure S2: The UV-Vis spectra of free ligands; Figure S3: Absorption spectra of free ligands and their complex in the presence of increasing DNA concentration; Table S6: The spectral parameters for the interaction of VO(IV) complexes and their ligands with DNA; Table S7: interaction table between the compounds and 1JU6 "Human Thymidylate Synthase" for Hepatocellular carcinoma protein and Breast cancer (1HK7) proteins; Table S8: Docking score and energy of the compounds and 1JU6 "Human Thymidylate Synthase" for Hepatocellular carcinoma protein; Table S9: Docking score and energy of the compounds and 1Hk7 Breast cancer protein; Table S10: 2D and 3D Docking interaction of the VO(II) complexes with colon cancer protein "Human Thymidylate Synthase" for Hepatocellular carcinoma protein (PDB code = 1JU6) and Breast cancer (PDB code = 1HK7) proteins; Figure