Next Article in Journal
Illicium verum (Star Anise) and Trans-Anethole as Valuable Raw Materials for Medicinal and Cosmetic Applications
Next Article in Special Issue
Thiazolidinedione Derivatives: In Silico, In Vitro, In Vivo, Antioxidant and Anti-Diabetic Evaluation
Previous Article in Journal
Gas-Phase Chemistry of 1,1,2,3,3,4,4-Heptafluorobut-1-ene Initiated by Chlorine Atoms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

QSAR Modeling, Molecular Docking and Cytotoxic Evaluation for Novel Oxidovanadium(IV) Complexes as Colon Anticancer Agents

by
Fatimah Y. Alomari
1,
Abeer A. Sharfalddin
2,
Magda H. Abdellattif
3,
Doaa Domyati
4,
Amal S. Basaleh
2 and
Mostafa A. Hussien
2,5,*
1
Chemistry Department, College of Science, Imam Abdulrahman Bin Faisal University, P.O. Box 76971, Dammam 31441, Saudi Arabia
2
Department of Chemistry, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
3
Department of Chemistry, College of Science, Taif University, Al-Haweiah, P.O. Box 11099, Taif 21944, Saudi Arabia
4
Department of Chemistry, College of Science, University of Jeddah, P.O. Box 80327, Jeddah 21589, Saudi Arabia
5
Department of Chemistry, Faculty of Science, Port Said University, Port Said 42521, Egypt
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(3), 649; https://doi.org/10.3390/molecules27030649
Submission received: 7 December 2021 / Revised: 7 January 2022 / Accepted: 12 January 2022 / Published: 19 January 2022

Abstract

:
Four new drug-based oxidovanadium (IV) complexes were synthesized and characterized by various spectral techniques, including molar conductance, magnetic measurements, and thermogravimetric analysis. Moreover, optimal structures geometry for all syntheses was obtained by the Gaussian09 program via the DFT/B3LYP method and showed that all of the metal complexes adopted a square-pyramidal structure. The essential parameters, electrophilicity (ω) value and expression for the maximum charge that an electrophile molecule may accept (ΔNmax) showed the practical biological potency of [VO(CTZ)2] 2H2O. The complexes were also evaluated for their propensity to bind to DNA through UV–vis absorption titration. The result revealed a high binding ability of the [VO(CTZ)2] 2H2O complex with Kb = 1.40 × 10⁶ M−1. Furthermore, molecular docking was carried out to study the behavior of the VO (II) complexes towards colon cancer cell (3IG7) protein. A quantitative structure–activity relationship (QSAR) study was also implemented for the newly synthesized compounds. The results of validation indicate that the generated QSAR model possessed a high predictive power (R2 = 0.97). Within the investigated series, the [VO(CTZ)2] 2H2O complex showed the greatest potential the most selective compound comparing to the stander chemotherapy drug.

1. Introduction

Cancer is one of the leading causes of death, the fatal characteristic of which is uncontrollable and irregular growth. In cancer, some cells of the diseased part start to grow abnormally and can metastasize to a different part of the body [1]. Among the different types of cancer, colon cancer represents the fourth most common type of cancer and the fifth major cause of cancer-associated death [1,2]. Thus, there is a need to effect corrective measures to reduce these numbers, such as finding influential anticancer agents and achieving early screening of the cancer cell, which is dependent on the patient and the cancer type [3].
In the past few years, there has been increased interest in developing additional transition metal compounds as anticancer drugs due to their lower toxicity than platinum-based agents [4]. Among the transitional metals, vanadium complexes have been studied and have shown practical activity towards various tumor cell types [5]. This finding suggests that vanadium complexes have beneficial therapeutic properties that could be extended to treat common types of cancer [5,6,7,8]. The most practical role of these anti-cancer agents is the checking and correcting of the division of cancer cells before they continue their amorphous growth [9,10,11]. There are massive studies that have explored the introduction of the metal ion to medication drugs to enhance their efficiency for pharmacological applications, as well as for several other unique therapeutic opportunities [12,13,14]. The new metal-based drugs have lower toxicity and high lipophilicity to penetrate the cell, and greater absorbance and stability than the free drugs [15]. Oxidovanadium(IV) complexes of interest are given in the work of Biswal and his colleagues, who synthesized three novel water-soluble oxidovanadium(IV) compounds [16], and two of the obtained complexes showed significant anticancer activities against the human hepatic carcinoma cell line Hep3B, which could be considered as a candidate for further studies. In this light, the current work aims to conduct descriptive and predictive research concerning anticancer activity for a new series of oxidovanadium(IV) complexes. Four drugs have an important medicinal feature, such as antibacterial (sulfonamide), antihistamine (cetirizine), antithyroid (carbimazole), and anti-inflammatory (Lornoxicam) activities, and have different donating atoms to react with oxidovanadium(IV) and to test their efficiency.
One of the approaches to accelerating the drug design process through the compounds’ advantages is through computer-aided drug design (CADD) [17,18]. Two approaches that are often used in CADD are quantitative structure–activity relationships (QSAR) and molecular docking (MD). The basis of QSAR is classification and regression techniques that can be characterized as a strategy attempting to foretell the characteristics, reactivity, and activity of a fuzzy set of molecules, based on the investigation of an equation that connects the molecular configurations to their specific measured property and activity [19,20]. Currently, QSAR is one of the most important methodologies in pharmacy and chemistry [21,22] used to develop mathematical models to identify statistically noteworthy correlations between toxicological, biological, binary, categorical, and continuous properties, and chemical structure. Additionally, molecular docking is based on the binding characteristics of the target and ligand. It creates various potential adduct structures, clustered and classified together through software that scores them. Examples of the areas where molecular docking is utilized include the design of drugs and in understanding the molecular interaction of drugs [23,24,25].
Therefore, an in silico QSAR study was performed to rationalize the structure–activity relationship using a molecular operating environment (MOE) for the studied metal complexes. This method can reduce the experimental time and effort required to examine all of the investigated compounds and focus on the promising molecules to further investigate by in vitro study. Additionally, the molecular docking technique was also carried out for all oxidovanadium(IV) compounds to exhibit the binding modes and interactions between the compounds and the colon cancer protein. An in vitro assay was used to confirm the in-silico simulation results.

2. Results and Discussion

2.1. Physical Properties

All synthesized VO(II) complexes were stable at room temperature and virtually insoluble in water and most organic solvents, but were soluble in some solvents, including DMSO and DMF. The molar conductance values of [VO(SO4)(CBZ)] 8H2O, [VO(CTZ)2] 2H2O and [VO(SO4)(SCZ)] 7H2O complexes at 1 × 10−3 M in DMSO solution were 7.47, 8.85, and 17.72 Ω−1 mol−1 cm2, respectively, as given in Table 1. These low values indicate that the complexes were nonelectrolytes. At the same time, the [VO(LOR)2] SO4 complex gave a value of 69.85 Ω−1 cm2 mol−1, reflecting the formation of a cationic molecule or the existence of a sulfate group in the coordination sphere [26,27].
The proposed complex structures were confirmed by elemental analysis, conductometrics, and physical measurements. The elemental analyses and physical properties of the ligands and their vanadium complexes are summarized in Table 1 and the suggested structure is shown in Figure 1.

2.2. The Molar Ratio of the Complexes

The stoichiometric environments around the oxidovanadium(IV) complexes were identified using a molar ratio technique [28]. The absorbance values were plotted toward the ratio of [M]/([M]+[L]) and are shown in Figure S1. The reflection line upon increasing the ligand concentration was around 0.5, which indicates the formation of [VO(SO4)(CBZ)] 8H2O and [VO(SO4)2(SCZ)] 7H2O complexes with 1:1 (M:L). In contrast, [VO(CTZ)2] 2H2O and [VO(LOR)2] SO4 complexes gave a reflection line around 0.33, which suggested a molar ratio of 1:2 (M:L).

2.3. Infrared Spectra (IR)

The IR spectra of the free ligands and their complexes are shown in Figure 2 and Figure S3. The main IR bands and their assignments are listed in Table 2. The IR spectra of complexes showed broad bands in the region of 3450–3540 cm−1, corresponding to the OH stretching vibrations of the lattice water molecules. The IR spectrum of the free CBZ ligand showed three bands for the ester group, at 1677, 1228, and 1096 cm−1 for the ν (C=O) and two ν (C–O) stretches, respectively. No notable change was seen in the ν(C=O) band, indicating that it did not coordinate with the vanadyl ion. Moreover, there was a noticeable shift in the ester oxygen bands to 1234 and 1107 cm−1, and these low frequencies confirmed the coordination of the ester oxygen with the vanadium ion [15]. The band at 1271 cm−1 related to the ν (C=S) vibration shifted to 1279 cm−1, which was a strong indication that the sulfur atom contributed to the complex was a donor [15]. A new band observed in the range 1310–1276 cm−1 was attributed to νas(SO4) and νs(SO4) bands.
In the [VO(CBZ)2] 2H2O complex, it is possible to observe the disappearance of the absorption band in the region around 1742 cm−1 assigned to the free carboxylic acid of CBZ, which indicated the formation of a bond between the vanadium ion and the carboxylate oxygen. The band corresponding to ν(C=O) was replaced by two bands at 1599 and 1409 cm−1, assigned to antisymmetric νasym(COO), and symmetric νsym(COO), respectively. The difference Δν [Δν = νasym(COO)–νsym(COO)] is a useful tool for determining the coordination mode of the carboxylate ligands [29]. The calculated value was in the range = 189 cm−1, which indicated a bidentate coordination mode for the cetirizine ligand [30,31].
The free LOR ligand showed bands in the ranges 1647 and 1591 cm−1, which were assigned to ν(C=O), and the ν(C=N) stretching frequency shifted to 1668 and 1621 cm−1, respectively, in the metal complex, indicating the involvement of the nitrogen of the pyridine ring and the oxygen of the amide group in the complex formation.
The IR spectrum of the free SCZ ligand showed a band corresponding to the ν(N-H) mode of the sulfonamide group and appeared at 3282 cm−1. This band shifted from the main position due to hydration or may have coordinated with metal [32,33]. The band at 1670 cm−1 assigned to ν(C=N) remained in complex, indicating the non-involvement of this group in the coordination sphere [34]. The ν(S-N) band at 1097 cm−1 shifted to 1105 cm−1, and hence confirmed the coordination of the metal center to the nitrogen of the sulfonamide. Two new bands observed in the [VO(SO4)(SCZ)] 7H2O spectrum at 1145–1011 cm−1 were assigned to νas(SO4) and νs(SO4) bands, in bidentate mode.
All four complexes exhibited a band at 955–960 cm−1, which can be ascribed to ν(V=O) stretching frequency. The band position of ν(V=O) agrees with the square-pyramidal structure [34]. Moreover, the bands in the regions of 560–580, 635–655, and 422 cm−1 may be assigned to (V–O), (V–N), and (V–S) stretching frequencies, respectively.

2.4. Electronic Spectroscopic and Magnetic Susceptibility

The UV–Vis spectroscopic analysis was performed for the ligands and their VO(II) complexes in DMSO in the range of 200–1000 nm and is shown in Figure 3 and Figure S4. The transitions in the free ligands attributed to π–π* and n–π* at the 260–280 and 311–373 nm regions in the UV spectra of the complexes shifted to low or high frequencies, confirming the coordination of the ligands to vanadyl ion. The observed weak-intensity bands in the range 850–900 nm in the metal complexes were assumed to metal-centered d–d transitions.
All synthesized complexes had magnetic moment values between 1.76 and 1.74 B.M. at 25 °C in accordance with the 3d1 electronic configuration of the VIVO2+ moiety. These values are compatible with the reported values for the square-pyramidal geometry [35]. The magnetic moment values and the electronic transition data for the compounds are presented in Table 3.

2.5. Electron Paramagnetic Resonance Spectra (EPR)

The X-band EPR spectra for the oxidovanadium(IV) complexes were acquired in DMF at room temperature. The stable isotope, 51V, of the vanadium (IV) ion carries electron and nuclear spins of S = 1/2 and I = 7/2, respectively. The EPR spectra of oxidovanadium(IV) complexes display eight absorptions due to the coupling of the vanadyl unpaired electron with the vanadium nuclear spin moment. These data suggest that just single mononuclear vanadium (IV) species containing a lone unpaired electron exist in the examined sample solution. The EPR parameters g and A values are presented in Table 4. The calculated values are consistent with a square-pyramidal structure [36,37]. Overall, the g values were less than those for the free electron: ge = 2.002. This phenomenon could be a consequence of the spin–orbit interplay of the ground state, dxy level, with the low-lying excited states [38]. The EPR spectra for oxidovanadium(IV) complexes are shown in Figure 4 and Figure S5.

2.6. Thermal Analysis (TGA)

Thermogravimetric studies of all complexes showed two decomposition stages. The obtained curves are displayed in Figure 5 and Figure S6. The initial step for [VO(SO4)(CBZ)] 8H2O occurring between a temperature of 25 °C and 218 °C corresponded to a loss of C4H7N2O3S and eight water molecules of hydration with weight loss 55.19% The second stage occurred within the temperature of 218 °C to 800 °C was related to loss of C5H3O2S with 27.55% mass loss, followed by the formation of VO2 as a residue.
The first decomposition step of [VO(CTZ)2] 2H2O complex was removed two lattice water molecules in between 55 and 100 °C with weight loss 4.09%. The organic decomposition began at ~120 °C with 86.46% The residue consisted of VO2 with a weight loss of 9.44%.
In contrast, the [VO(LOR)2] SO4 complex did not show a weight loss up to 119 °C corresponds to the absence of lattice water molecules. Therefore, the weight loss in the range 119–341 °C has corresponded to the loss of C7H6Cl2O10S2 with a total weight loss of 42.36%. At the same time, a weight loss of 48.26% within the temperature range 341–552 °C is assigned to losing C19H12N6OS3 and forming VO2; as a residue with a final weight of 9.38%.
The complex [VO(SO4)(SCZ)] 7H2O lost C6H8N2O4S2 and seven hydrated water molecules with a weight loss of 69.37%, then followed by losing 16.22% of the mass in a temperature range of 232–998 °C which assigned to C4HN2O. The final product is metal residue with a total mass loss amounts to 14.40%. The TGA data for the VO(II) complexes present in Table S2.

2.7. Kinetic Analysis for Thermogravimetric Data

The kinetic parameters were calculated using Coast–Redfern (CR) and Horowitz–Metzger (HM) methods. The kinetic parameters were obtained via equations presented in previous work [39] {Sharfalddin AA, 2021 #63} and showed in Table S3. The values of Ea indicate that the complexes are extremely stable. The values of ΔS are negative and indicated slower thermal decomposition than the typical reaction [40]. Positive values of ΔH and ΔG indicate that all decomposition steps are endothermic and non-spontaneous processes. The high values of ΔG for each stage are due to the increase in TΔS values from one step to the next. The CR and HM curves for the oxidovanadium(IV) complexes are represented in Table S4.

2.8. Molecular Modeling of Investigated Compounds

2.8.1. Geometry Description

The results obtained from elemental analysis, molar ratio, molar conductance, spectroscopic techniques, and thermal studies were used to build the input files for all synthesized oxidovanadium(IV) complexes and free ligands, running the jobs using the DFT method through the 6–31G and LANL2DZ basis sets. The best-optimized geometries and the bond length of V=O were visualized by Gauss View and are presented in Table 5 and Figure S7. The formed structure of the metal complexes is square-pyramidal geometry with negligible distortion. The bond length around the metal ion was in the range 1.95–2.60 Aº between the donating atoms and the vanadyl ion, which showed a range of small and medium ionic character (bond length > 2 Å) [41]. On the other hand, there was a direct correlation between the length of a metal–oxygen bond and its stretching frequency. Thus, we can directly determine the metal–oxygen bond length from the measurement of stretching frequencies, calculated by the below equation [42,43]:
v = 21,349 exp(−1.9176 R (Å))
where v is the stretching frequency for V=O in cm−1, and R is the V=O bond length in angstroms. The stretching frequencies and the respective metal−oxygen bond lengths have an inverse relationship. This means that a higher stretching frequency corresponds to a shorter metal−oxygen bond length. Therefore, it is possible to utilize IR spectra information to compute the V=O bond length.
A comparison of the experimental frequencies and the extracted theoretical data is presented in Table 6. The obtained experimental bond length by the above equation for the investigated compounds was in the range 1.603–1.616 Å and correlated with the theoretical data.
The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) represent the nucleophilicity of the molecule (electron donating) and the electrophilicity (electron accepting) of the molecule, respectively. Thus, the images and energy values of HOMO and LUMO were extracted in the ground state to study the electronic characteristics of molecular systems for the ligand and around the metal center. Table S5 compares the molecular orbital diagrams between the ligand and the oxidovanadium(IV) complexes. Interestingly, both HOMO and LUMO orbitals in the free CTZ were concentrated in the ((4-chlorophenyl)(phenyl)methyl)-4-methyl piperazine part. After complexing with a metal ion, the HOMO orbital appeared prolonged over a 1-((4-chlorophenyl) (phenyl)methyl)-4-methyl piperazine part of the complex, while the LUMO orbital was located around the central metal. In contrast, a free LOR ligand showed its HOMO over the whole molecule, except the pyridine ring, while the LUMO level presented in 2-chlorothiophene. These locations of the HOMO and LUMO were maintained after reaction with VO(II).

2.8.2. Global Reactivity Descriptors

Using details from these molecular orbital compositions and energy levels, global reactivity descriptors were calculated as the following: energy gap (∆E = Egap = ELUMO −EHOMO), absolute electronegativities (χ = −EHOMO + ELUMO/2), absolute hardness (η = ELUMO − EHOMO/2), chemical potentials (μ = −w), global softness (S = 1/2η), and global electrophilicity (ω = π2/2η) [39] and are shown in Table 7.
The energy gap, Egap, shows a practical guide for the reactivity of the compound [39]. Decreasing the level between the HOMO and LUMO orbitals led to lowering the energy gap Egap, thus enhancing the transition of the electrons and the softness of the complex. We can observe from Table 7 that complexing the ligand significantly minimized the energy gap of the complexes. The oxidovanadium(IV) of LOR has the lowest energy gap among the investigated complexes. Moreover, the global softness (S) has moderate values, indicating the minor characteristics of synthesized compounds. The electrophilicity (ω) is an expression of the ability to obtain electronic charge from other molecules in the environment, increasing the biological activity descriptor. The highest value appeared with the VO–CTZ complex (5.03 eV), followed by VO–SCZ and VO–CBZ, which reveal their biological potency. The ΔNmax is an expression of the maximum charge an electrophile may accept from the environment [44]. We could arrange the metal complexes depending on this value in descending order, VO–CTZ > VO–LOR > VO–SCZ > VO–CBZ, which represents the capacity of VO–CTZ to be an anticancer agent.

2.8.3. DNA Binding Study

A DNA-binding assay is an effective technique for investigating a new compound’s biological activity by studying its ability to interfere with DNA replication and transcription [45,46]. Absorption titration is a standard and easy approach to determine the interaction mode and the nature of DNA binding [47,48]. The spectral parameters for the interaction of the tested complexes with DNA are presented in Table 8, whereas the chromic absorption spectra are presented in Figure 6 and Figure S8.
All compounds showed an increasing hyperchromic shift with increasing CT–DNA concentration, except the LOR ligand, which showed a hypochromic shift. The hyperchromic effect implies that the ligand and its complex bind to the external of the major or minor grooves of DNA via electrostatic binding [49,50]. Moreover, hyperchromic observation indicates the alterations in the structure of the genetic material following engagement with the complex [45]. Otherwise, the hypochromic effect is an intercalative reaction between the aromatic part of the compound with DNA bases [51,52]. Obviously, the complex [VO(CTZ)2] 2H2O showed the highest Kb value (1.40 × 106), which implies that the compound has a better binding ability compared with other compounds.

2.8.4. Molecular Docking

Molecular docking analysis plays a significant part in investigating the complexation of a ligand at the active site of a receptor. Thus, molecular docking analysis was conducted to elucidate how the free ligands and the oxidovanadium(IV) complex interact with the investigated protein. The dock procedures were evaluated against colon cancer (3IG7), which was chosen based on previous research suggesting its inhabitation as colon cancer therapy [53,54]. The binding energy scores for free ligands and their synthesized oxidovanadium(IV) compounds are presented in Table 9. The oxidovanadium(IV) complexes initially exhibited more potent binding energy than their free ligands against 3IG7 protein. They indicated that the modified vanadyl complexes enhance the binding affinity of the free ligands toward the examined protein. At the same time, the [VO(CTZ)2] 2H2O complex has the highest binding with the lowest negative scoring of −9.81 KJ mol−1, which is related to the π-acceptor interaction between a six-membered ring containing a nitrogen atom in ASN 132 amino acid with an energy score of −1.6 kcal mol−1. The docking scores for all oxidovanadium(IV)) complexes and free ligands are shown in Table 10. The interactions between all complexes and 3IG7 proteins are represented in Table S6.

2.8.5. QSAR Study

The common applications of QSAR modeling include the evaluation of the structural characteristics of a compound, and the design and prediction of the activities of new compounds. Therefore, MOE software was applied to develop a QSAR model of oxidovanadium(IV) complexes. Several descriptors, including ASA_P, ASA_H, h_pkp, and dipole moment, were utilized to display the predicted IC50 values. These were close to the actual IC50 values for the training set, as presented in Figure S2. Out of the MLR analysis, the model exhibited a regression coefficient (R) of 0.9882 and coefficient of determination (R2) of 0.9766, thus explaining the 97.66% IC50 activity. The R and R2 figures portray the accuracy of the QSAR model.
The plot of the relationship between the predicted values obtained from MLR and the experimental data is shown in Figure 7. Some predicted and experimental values exhibit a close correlation, and residual (RE) values were close to zero, indicating that the predictions are reliable, as shown in Table S7. According to the results, the approach was recorded as the QSAR model and applied to predict the compound activities of the test data set. Based on the QSAR model, a fit was used to obtain the predicted activity (PRED IC50) data of the synthesized oxidovanadium(IV) complexes and their free ligands (testing set), and the results are presented in Table 11. The [VO(CTZ)2] 2H2O complex exhibited the best predicted activity (PRED IC50 = 1.45 µM), which indicated that the complex was more effective against colon cancer. On the other hand, the complexes [VO(SO4)(SCZ)] 7H2O and [VO(SO4)(CBZ)] 8H2O showed the lowest predicted activities against colon cancer. The other members showed a moderate predicted activity, which ranged between PRED IC50 values of 2.61 and 32.87 µM. Interestingly, the predicted IC50 values of [VO(SO4)(CBZ)] 8H2O, [VO(CTZ)2] 2H2O, and [VO(LOR)2] SO4 complexes were higher than their free ligands, while the other complexes were less active than their free ligands.

2.9. Anticancer Study

Based on the results obtained from the QSAR and docking studies, the [VO(CTZ)2] 2H2O complex was selected for examination against the human colon cancer cell line (HCT116) to evaluate its anticancer properties. The cytotoxicity of cisplatin against the HCT116 cell line was used to benchmark the anticancer activity of the oxidovanadium(IV) complex. The results are listed in Table 12. The IC50 for [VO(CTZ)2] 2H2O was 2.11 μM, indicating potent anticancer activity compared with the standard chemotherapy drug and with cisplatin, which has an IC50 of 44.13 μM. Furthermore, the obtained value strongly correlates with the predicted score from the QSAR equation (PRED IC50 = 1.45 µM). The selectivity of the compound showed low toxicity toward normal cell line (LLC-MK2) with an IC50 value of 649.8 µM which illustrated good selectivity and practical activity.
Comparing this result with previous work, the Schiff base oxidovanadium(IV) complex were tested against the colon cancer cell line (HCT116), with an IC50 = 2.39 µM. Hence, the [VO(CTZ)2] 2H2O has higher potency and could be considered as an anticancer agent [55].

3. Materials and Methods

3.1. Materials and Reagent

VOSO4 was supplied from Sigma Aldrich. Carbimazole (CBZ) and Sulfaclozine (SCZ) were obtained from Alibaba (China). Lornoxicam (LOR) and Cetirizine (CTZ) were obtained from Tabuk Pharmaceuticals (Tabuk, Saudi Arabia). Calf thymus DNA (CT-DNA) from SIGMA type 1A36 was used for the DNA studies. All the other solvents were obtained from Sigma-Aldrich (St. Louis, MO, USA) and were used as received.

3.2. Apparatuses

The element’s content (CHN) was observed using Perkin-Elmer 2400 CHN Elemental Analyzer instrument. Melting points were measured by MPA100–Automated melting point system. The conductivity measurements are determined using a conductor TDA meter at room temperature in 10−3 M DMSO solutions. IR spectra were recorded using solid samples on a PerkinElmer Spectrum 100 instrument (Waltham, MA, USA); peaks are reported over the range 4000–400 cm−1. The electronic absorption spectra were determined in the DMSO solution for the free ligands and their complexes using a MultiSpec-1501 UV–Vis spectrophotometer with 1 cm quartz cells; bands in the range 200–1100 nm. The X-band EPR spectra for all the complexes are recorded as in DMF solution at room temperature from the continuous wave Bruker EMX PLUS spectrometer (Bruker BioSpin, Rheinstetten, Germany) and collected with Bruker Xenon software with DPPH standard (υ = 9.440103 GHz). Thermal analysis was carried out for the complexes under a nitrogen atmosphere, a constant heating rate was applied by 10 °C min−1 at a temperature range from 25 to 1000 °C using the PerkinElmer TGA system.

3.3. Molar Ratio Method

A standard solution of the metal ion was consistent at 0.72 × 10−4 M and divided between seven volumetric flasks by pipette into the volumes of 0, 1, 2, 3, 4, 5, and 6 mL. Then, various ligand concentrations were added in 6, 5, 4, 3, 2, 1, and 0 mL volumes, as presented in Table S1. All measurements were performed in the 200 to 500 nm spectrum [28,56].

3.4. Synthesis of VO(II) Complexes

All the metal complexes were prepared by adding a 3 mL solution of VOSO4.5H2O (1 mol) to 25 mL of ethanol to dissolve the free ligand, according to the stoichiometric amount extracted from the molar ratio experiment. The obtained mixture was refluxed for 2 h and the formed precipitate was filtered, collected, and then washed several times with diethyl ether, before being oven-dried overnight at 55 °C.
[VO(SO4)(CBZ)] 8H2O: Yield 63%, m.p. 125 °C. IR (KBr, cm−1): 3275 ν(OH2); 1310 νas(SO2); 1035 νs(SO2); 1234, 1107 ν(C-O); 1081 ν(C-Cl); 980 ν(V=O); 596 ν(V-O); 422 ν(V-S). Elemental analysis of C7H26N2O15S2V (%): Found: C 17.04; H 5.31; N 5.68; Anal. Calc. (%): C 17.25; H 5.32; N 5.60.
[VO(CTZ)2] 2H2O: Yield 89%, m.p. 134 °C. IR (cm−1): 3394 ν(OH2); 1599 νasym(COO); 1409 νsym(COO); 963 ν(V=O); 544 ν(V-O). Elemental analysis of C42H52Cl2N4O9V (%): Found: C 57.41; H 5.96; N 6.38; Anal. Calc. (%): C 57.43; H 5.97; N 6.41.
[VO(LOR)2] SO4: Yield 68%, m.p. 124 °C. IR (cm−1): 3348 ν(OH); 1668 ν(CONH); 1621 ν(C=N); 1325 νasym(SO2); 1034 νsym(SO2); 968 ν(V=O); 521 ν(V-O); 469 ν(V-N). Elemental analysis of C26H20Cl2N6O11S5V (%): Found: C 35.71; H 2.31; N 9.61; Anal. Calc. (%): C 35.73; H 2.30; N 9.54.
[VO(SO4)(SCZ)] 7H2O: Yield 61%, m.p. 209 °C. IR (cm−1): 3445 ν(OH2); 3282 ν(N-H); 1670 ν(C=N); 1271 νasym(SO2); 1033 νsym(SO2); 1105 ν(S-N); 968 ν(V=O); 589 ν(V-O); 424 ν(V-N). Elemental analysis of C10H23N4O14S2V (%): Found: C 20.93; H 4.04; N 9.76; Anal. Calc. (%): C 30.10; H 4.13; N 9.98.

3.5. Molecular Modeling Method

Theoretical calculations for the free ligands and their metal complexes were obtained using the Gaussian 09 package [57]. The geometry optimization was carried out by the density function theory method (DFT) with Becke’s three-parameter hybrid and Lee–Yang–Parr correlation functional method (B3LYP). The basis set (6–31G) was used for the organic molecules and the correlation-consistent basis set (LANL2DZ) as its preferred basis set for the metal complex [58]. The extracted data from the log files showed that the vibration frequency did not have any imaginary frequencies and that the molecule was a minimum on the potential energy surface. Gauss View software was used to visualize the fchk files and extract the ELUMO and EHOMO levels and energy orbitals to estimate the essential quantum parameters using the following equations: energy gap (Egap = ELUMO − EHOMO), absolute electronegativities (χ = −EHOMO +ELUMO/2), absolute hardness (ɳ = ELUMO − EHOMO/2), chemical potentials (μ = −χ), global softness (S = 1/2η), and global electrophilicity (ω = π2/2η) [41].

3.6. Molecular Docking Method

Molecular docking study utilized the MOE 2019.0102 program [59]. All optimized drugs and their metal complexes from the density functional theory method (DFT) were used in the docking process. Crystal configurations of the HCT116 (PDB code = 3ig7) were obtained from the Protein Data Bank (http://www.rcsb.org/pdb/ (accessed on 1 December 2021)) [60]. Any attached water or ligand cofactors, together with the water moieties in proximity to the protein, were deleted; a process followed by the correction of the protein structure, and then the addition of hydrogen atoms. An MMFF94x force field was used to allocate parameters and charges. Following the MOE synthesis of alpha-site spheres, the configured complex models were docked on the leading active site using MOE’s DOCK module [61,62,63]. A London dG scoring function was used for dock scoring within the MOE software, enhanced with the use of unconnected refinement techniques, and the superlative ten binding positions were directed for analysis to attain the optimum score. A database browser was recruited to perform a comparison of the docking poses with the ligand within the co-crystallized configuration and to generate their RMSD. The RMSD and the native ligand’s method of engagement within the receptor’s configuration were utilized as a standard docked model. The same protocol was utilized to dock the co-crystalline compound“N-{1-[cis-3-(acetylamino)cyclobutyl]-1H-imidazol-4-yl}-2-(4-methoxyphenyl)acetamide” which downloaded with the 3IG7 protein and treated as positive control while for negative control we used four different amino acids. For all docking result we get same docking site with the positive control and −7.30 kcal/mol with hydrogen bond of 2.5 Å. Where the negative control compounds do not show any interaction with very small docking score Tables S8 and S9.

3.7. QSAR Technique

3.7.1. Data Set Used

A QSAR assessment was conducted on a set of 100 compounds to examine their colon anticancer activity. The compounds were grouped differently for testing and training. The literature shows that the imidazole derivative molecules and vanadyl complexes can be utilized to counter colon cancer activity. The structure and biological activities of the compounds obtained in the experiment are represented in Figure S2. Moreover, the synthesized vanadyl complexes and their ligands were used as test data. The training data set obtained from the biological activity (IC50) was used as the dependent variable, linearly correlated with the independent variables.

3.7.2. Physiochemical Descriptors

A descriptor is a mathematical argument describing molecular properties according to biological activity and structure. The physicochemical descriptors calculated were total polar surface carrier (ASA_P), total hydrophobic surface (ASA_H), pKb of the reaction to add a proton at pH 7 (h_pKb), and dipole moment (dipole).

3.7.3. Development of Experimental Data Model

The QSAR model was formulated based on anti-cancer activities (IC50) (as the dependent variable), while descriptors were used as model fields. Multi-linear regression (MLR) analysis was applied to obtain a training data set. MLR depicts a direct linkage between the independent variable X, which represents the molecular descriptors, and the dependent variable Y, which represents the biological activity. In an MLR approach, the average value of Y is dependent on X. Below is an MLR equation that includes more than a single independent variable (descriptors), while maintaining one response variable (activity).
Y = bx 1 + bx 2 + bx 3 + bx 4 + C
where Y is the dependent variable, b represents the regression coefficients of a varying x (independent variable), and C represents a regression constant.

3.8. Biological Applications

3.8.1. CT-DNA Interaction

All studies with CT–DNA were carried out in Tris-HCl buffer solution (pH 7.4). The UV absorbance ratio (A260/A280) was 1.8–1.9, confirming that the DNA was free of protein. Electronic absorption spectra were obtained with a constant compound concentration at varying CT–DNA concentrations. The Benesi–Hildebrand equation (Equation (2)) was used to calculate the effective binding constants (Kb) [45].
A 0 A A 0 = Ԑ G Ԑ H _ G Ԑ G + Ԑ G Ԑ H _ G Ԑ G × 1 K [ DNA ]
where K is the association/binding constant; A0 and A are the absorbances of the drug and its complex with DNA, respectively; and Ԑ G and Ԑ H _ G are the absorption coefficient of the drug and the drug–DNA complex, respectively. The intercept to slope ratios from the graphic representation of A 0 /( A A 0 ) vs. 1/[DNA] yielded the association constant. Graphs depicting [DNA] vs. [DNA]/ Ԑ a Ԑ f , where Ԑ a (or Ԑ G ) represents the absorption coefficient of the drug and Ԑ f (or Ԑ H _ G ) is the absorption coefficient of the drug–DNA complex.

3.8.2. In Vitro Antitumor Activity

HCT116 colon cancer cell line was purchased from the Egyptian Holding Company for Biological Products & Vaccines (VACSERA), Giza, Egypt, and then maintained in the tissue culture unit (Pharmacy College Taif University, Taif, Saudi Arabia). The cells were cultured in RPMI-1640 medium, supplemented with 10% heat-inactivated FBS, 50 units/mL of penicillin, and 50 mg/mL of streptomycin, and maintained at 37° in a humidified atmosphere containing 5% CO2. The cells were maintained as a “monolayer culture” by serial subculturing. Cell culture reagents were obtained from Lonza (Basel, Switzerland). Cytotoxicity was determined using the SRB method as previously described by Skehan et al. Exponentially growing cells were collected using 0.25% Trypsin–EDTA and seeded in 96-well plates at 1000–2000 cells/well in RPMI-1640-supplemented medium. After 24 h, cells were incubated for 72 h with various concentrations of the tested compounds. Following 72 h treatment, the cells were fixed with 10% trichloroacetic acid for 1h at 4 ºC. Next, wells were stained for 10 min at room temperature with 0.4% SRB (sulforhodamine B) dissolved in 1% acetic acid. Finally, the plates were air-dried for 24 h, and the dye was solubilized with Tris-HCl for 5 min on a shaker at 1600 rpm. Each well’s optical density (OD) was measured spectrophotometrically at 564 nm with an ELISA microplate reader (ChroMate-4300, FL, USA). The IC50 values were calculated according to the Boltzmann sigmoidal concentration–response curve equation using the nonlinear regression fitting models. The cytotoxic activity of the investigated compound was tested against a normal cell line aiming to get evidence about the selectivity of the compound.

4. Conclusions

In this study, four oxidovanadium(IV) complexes have been synthesized and characterized by several spectroscopic methods. The IR spectra suggested that the ligands have bidentate coordination mode to the vanadium ion. In addition, the molar conductance, EPR, magnetic moment values, and electronic data support a square pyramidal structure for all complexes.
Computational studies were applied to demonstrate the optimization geometry and essential quantum parameters and confirm their biological efficiency. The metal complexes showed a square-pyramidal geometry arrangement around the metal ion which was agreeable with experimental results. The form coordination bonds length presented strength bonded between oxidovanadium(IV) and the investigated ligand to form stable complexes. The bioactivity study was varied and involved DNA-binding, molecular docking, QSAR, and cytotoxicity analyses. DNA binding results revealed two behaviors for the synthesized oxidovanadium(IV) complexes with an increase in CT-DNA concentration, hyperchromic with electrostatic or grooves binding modes and hypochromic signifying an intercalation binding mode. Molecular docking results showed that the [VO(CTZ)2] 2H2O complex exhibited significant interaction with colon cancer (3IG7) Protein with good selectivity.
QSAR study has given significant information on biological activity by using the MLR method. QSAR model showed a good correlation between the predicted and the experimentally observed inhibitory activities. Based on the molecular docking and QSAR results, [VO(CTZ)2] 2H2O was selected and tested for its inhibitory activity against colon cancer cell line (HCT116). The selected complex showed higher anticancer activity than the standard cisplatin chemotherapy drug.

Supplementary Materials

Table S1: The experimental data of oxidovanadium(IV) complexes, Table S2: Thermogravimetric analysis data for the oxidovanadium(IV) complexes, Table S3: Kinetic parameters of thermal decomposition steps for oxidovanadium(IV) complexes, Table S4: Coats-Redfern (CR) and Horowitz–Metzger (HM) of the oxidovanadium(IV) complexes, Table S5: HOMO and LUMO plots of the free ligands and their metal complexes using DFT/B3LY, Table S6: Docking interaction calculation of all oxidovanadium(IV) complexes with colon cancer protein (3IG7), Table S7: The experimental and the predicted IC50 values of some of the training set against HCT116 cell, Table S8: Docking score and energies of Co-crystalline compound (positive control) and four amino acids (negative control), Table S9: 2D and 3D Docking interaction of the Co-crystalline compound (positive control) and four amino acids (negative control), Figure S1: Molar ratio method plots of oxidovanadium(IV) complexes, Figure S2: Some of the structures of imidazole and oxidovanadium(IV) compounds (training) that were used in the QSAR with their activity against HCT 116, Figure S3:The FT-IR spectra of free ligands and their oxidovanadium(IV) complex, Figure S4: UV-vis Spectra of free ligands and their oxidovanadium(IV) complexes (inset shows the d-d transition the metal complex), Figure S5: EPR spectra of oxidovanadium(IV) complexes in DMF and in room temperature, Figure S6: The TGA curves of synthesized oxidovanadium(IV) complex, Figure S7: The optimized geometry of the free ligands gas phase using DFT method LANL2DZ basis sets, Figure S8: Absorption spectra of free ligands and their complex in the presence of increasing DNA concentration.

Author Contributions

Synthesis and characterizations, F.Y.A.; Supervision and writing original draft preparation, A.S.B. and D.D.; Theoretical calculation, review, and editing; A.A.S.; biological studies, interpretations, M.H.A.; Supervision and the conception, M.A.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All relevant data are within the manuscript and its Supporting information.

Acknowledgments

The authors are thankful to the King Fahad Center for Medical Research for providing the necessary laboratory facility. The simulation in this work was performed at King Abdulaziz University’s High-Performance Computing Center (Aziz Supercomputer) (http://hpc.kau.edu.sa (accessed on 30 July 2021)). The authors also thank Taif University Researchers Supporting Project Number TURSP-2020/91, Taif University, Saudi Arabia.

Conflicts of Interest

The authors declare no competing financial interest and the authors received no specific funding for this work.

Sample Availability

Samples of the compounds available from the authors.

References

  1. Ma, L.; Xu, G.B.; Tang, X.; Zhang, C.; Zhao, W.; Wang, J.; Chen, H. Anti-cancer potential of polysaccharide extracted from hawthorn (Crataegus.) on human colon cancer cell line HCT116 via cell cycle arrest and apoptosis. J. Funct. Foods 2020, 64, 103677. [Google Scholar] [CrossRef]
  2. Sun, Y.; Mironova, V.; Chen, Y.; Lundh, E.P.; Zhang, Q.; Cai, Y.; Vasiliou, V.; Zhang, Y.; Garcia-Milian, R.; Khan, S.A. Molecular pathway analysis indicates a distinct metabolic phenotype in women with right-sided colon cancer. Transl. Oncol. 2020, 13, 42–56. [Google Scholar] [CrossRef]
  3. Ngameni, B.; Cedric, K.; Mbaveng, A.T.; Erdoğan, M.; Simo, I.; Kuete, V.; Daştan, A. Design, synthesis, characterization, and anticancer activity of a novel series of O-substituted chalcone derivatives. Bioorg. Med. Chem. Lett. 2021, 35, 127827. [Google Scholar] [CrossRef]
  4. Nardon, C.; Boscutti, G.; Fregona, D. Beyond platinums: Gold complexes as anticancer agents. Anticancer Res. 2014, 34, 487–492. [Google Scholar] [PubMed]
  5. Evangelou, A.M. Vanadium in cancer treatment. Crit. Rev. Oncol. Hematol. 2002, 42, 249–265. [Google Scholar] [CrossRef]
  6. Sgarbossa, S.; Diana, E.; Marabello, D.; Deagostino, A.; Cadamuro, S.; Barge, A.; Laurenti, E.; Gallicchio, M.; Boscaro, V.; Ghibaudi, E. Synthesis, characterization and cell viability test of six vanadyl complexes with acetylacetonate derivatives. J. Inorg. Biochem. 2013, 128, 26–37. [Google Scholar] [CrossRef] [PubMed]
  7. Sochor, V.; McCaffrey, V.; Rabquer, B. Vanadium complexes inhibit the growth of HT-29 and MCF-7 cancer cell lines. FASEB J. 2015, 29, 978.7. [Google Scholar] [CrossRef]
  8. Kostova, I. Titanium and vanadium complexes as anticancer agents. Anti Cancer Agents Med. Chem. 2009, 9, 827–842. [Google Scholar] [CrossRef]
  9. Banerjee, A.; Dash, S.P.; Mohanty, M.; Sanna, D.; Sciortino, G.; Ugone, V.; Garribba, E.; Reuter, H.; Kaminsky, W.; Dinda, R. Chemistry of mixed-ligand oxidovanadium (IV) complexes of aroylhydrazones incorporating quinoline derivatives: Study of solution behavior, theoretical evaluation and protein/DNA interaction. J. Inorg. Biochem. 2019, 199, 110786. [Google Scholar] [CrossRef]
  10. Crans, D.C.; Yang, L.; Haase, A.; Yang, X. Health benefits of vanadium and its potential as an anticancer agent. Met. Ions Life Sci. 2018, 18, 251–279. [Google Scholar]
  11. Naglah, A.M.; Refat, M.S.; Al-Omar, M.A.; Bhat, M.A.; AlKahtani, H.M.; Al-Wasidi, A.S. Synthesis of a vanadyl (IV) folate complex for the treatment of diabetes: Spectroscopic, structural, and biological characterization. Drug Des. Dev. Ther. 2019, 13, 1409. [Google Scholar] [CrossRef] [Green Version]
  12. Warra, A. Transition metal complexes and their application in drugs and cosmetics-a Review. J. Chem. Pharm. Res. 2011, 3, 951–958. [Google Scholar]
  13. Hariprasath, K.; Deepthi, B.; Babu, I.S.; Venkatesh, P.; Sharfudeen, S.; Soumya, V. Metal complexes in drug research—A review. J. Chem. Pharm. Res. 2010, 2, 496–499. [Google Scholar]
  14. Sharfalddin, A.A.; Emwas, A.-H.M.; Jaremko, M.; Abdulala, M.A. Practical and Computational studies of Bivalence Metal Complexes of Sulfaclozine and Biological Studies. Front. Chem. 2021. [Google Scholar] [CrossRef]
  15. Sharfalddin, A.A.; Hussien, M.A. Bivalence Metal Complexes of Antithyroid Drug Carbimazole; Synthesis, Characterization, Computational simulation, and Biological Studies. J. Mol. Struct. 2020, 1228, 129725. [Google Scholar] [CrossRef]
  16. Biswal, D.; Pramanik, N.R.; Chakrabarti, S.; Drew, M.G.; Acharya, K.; Chandra, S. Syntheses, crystal structures, DFT calculations, protein interaction and anticancer activities of water soluble dipicolinic acid-imidazole based oxidovanadium (iv) complexes. Dalton Trans. 2017, 46, 16682–16702. [Google Scholar] [CrossRef] [PubMed]
  17. Baig, M.H.; Ahmad, K.; Rabbani, G.; Danishuddin, M.; Choi, I. Computer aided drug design and its application to the development of potential drugs for neurodegenerative disorders. Curr. Neuropharmacol. 2018, 16, 740–748. [Google Scholar] [CrossRef] [PubMed]
  18. Sliwoski, G.; Kothiwale, S.; Meiler, J.; Lowe, E.W. Computational methods in drug discovery. Pharmacol. Rev. 2014, 66, 334–395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Adeniji, S.E.; Uba, S.; Uzairu, A. Multi-linear regression model, molecular binding interactions and ligand-based design of some prominent compounds against Mycobacterium tuberculosis. Netw. Model. Anal. Health Inform. Bioinform. 2020, 9, 8. [Google Scholar] [CrossRef]
  20. Ponzoni, I.; Sebastián-Pérez, V.; Martínez, M.J.; Roca, C.; De la Cruz Pérez, C.; Cravero, F.; Vazquez, G.E.; Páez, J.A.; Díaz, M.F.; Campillo, N.E. QSAR classification models for predicting the activity of inhibitors of beta-secretase (BACE1) associated with Alzheimer’s disease. Sci. Rep. 2019, 9, 9102. [Google Scholar] [CrossRef] [Green Version]
  21. Neves, B.J.; Braga, R.C.; Melo-Filho, C.C.; Moreira-Filho, J.T.; Muratov, E.N.; Andrade, C.H. QSAR-based virtual screening: Advances and applications in drug discovery. Front. Pharmacol. 2018, 9, 1275. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Verma, J.; Khedkar, V.M.; Coutinho, E.C. 3D-QSAR in drug design-A review. Curr. Top. Med. Chem. 2010, 10, 95–115. [Google Scholar] [CrossRef]
  23. Fan, J.; Fu, A.; Zhang, L. Progress in molecular docking. Quant. Biol. 2019, 7, 83–89. [Google Scholar] [CrossRef] [Green Version]
  24. Salmaso, V.; Moro, S. Bridging molecular docking to molecular dynamics in exploring ligand-protein recognition process: An overview. Front. Pharmacol. 2018, 9, 923. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Torres, P.H.; Sodero, A.C.; Jofily, P.; Silva, F.P., Jr. Key topics in molecular docking for drug design. Int. J. Mol. Sci. 2019, 20, 4574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Ali, I.; Wani, W.A.; Saleem, K. Empirical formulae to molecular structures of metal complexes by molar conductance. Synth. React. Inorg. Met. Org. Nano Met. Chem. 2013, 43, 1162–1170. [Google Scholar] [CrossRef]
  27. Geary, W.J. The use of conductivity measurements in organic solvents for the characterisation of coordination compounds. Coord. Chem. Rev. 1971, 7, 81–122. [Google Scholar] [CrossRef]
  28. Tirmizi, S.A.; Wattoo, F.; Sarwar, S.; Anwar, W.; Wattoo, F.H.; Memon, A.N.; Iqbal, J. Spectrophotometric study of stability constants of famotidine-Cu (II) complex at different temperatures. Arab. J. Sci. Eng. 2009, 34, 43–48. [Google Scholar]
  29. Ashafaq, M.; Khalid, M.; Raizada, M.; Ahmad, M.S.; Khan, M.S.; Shahid, M.; Ahmad, M. A Zn-based fluorescent coordination polymer as bifunctional sensor: Sensitive and selective aqueous-phase detection of picric acid and heavy metal ion. J. Inorg. Organomet. Polym. Mater. 2020, 30, 4496–4509. [Google Scholar] [CrossRef]
  30. Iqbal, M.; Karim, A.; Ullah, I.; Haleem, M.A.; Ali, S.; Tahir, M.N.; Abbas, S.M. Synthesis, characterization, structural description, TGA, micellization behavior, DNA-binding and antioxidant activity of mono-, di-and tri-nuclear Cu (II) and Zn (II) carboxylate complexes. J. Coord. Chem. 2021, 74, 762–778. [Google Scholar] [CrossRef]
  31. Iftikhar, S.H.; Gilani, S.R.; Taj, B.M.; Raheel, A.; Termizi, S.A.; Al-Shakban, M.; Ali, H.M. Design, synthesis and biological evaluation of organotin (IV) complexes of flumequine and cetirizine. J. Serb. Chem. Soc. 2018, 83, 425–437. [Google Scholar] [CrossRef] [Green Version]
  32. Sun, B.; Jin, Q.; Tan, L.; Wu, P.; Yan, F. Trace of the interesting “V”-shaped Dynamic Mechanism of Interactions between Water and Ionic Liquids. J. Phys. Chem. B 2008, 112, 14251–14259. [Google Scholar] [CrossRef] [PubMed]
  33. Rocha, M.; Piro, O.E.; Echeverría, G.A.; Pastoriza, A.C.; Sgariglia, M.A.; Soberón, J.R.; Gil, D.M. Co (II), Ni (II) and Cu (II) ternary complexes with sulfadiazine and dimethylformamide: Synthesis, spectroscopic characterization, crystallographic study and antibacterial activity. J. Mol. Struct. 2019, 1176, 605–613. [Google Scholar] [CrossRef]
  34. Jone Kirubavathy, S.; Velmurugan, R.; Tamilarasan, B.; Karvembu, R.; Bhuvanesh, N.; Chitra, S. Synthesis, Characterization, Single-Crystal XRD, and Biological Evaluation of Nickel (II) Salen Sulfadiazine Complex. Synth. React. Inorg. Met. Org. Nano Met. Chem. 2016, 46, 1751–1758. [Google Scholar] [CrossRef]
  35. Borge, V.V.; Patil, R.M. Stability of Vanadium Chalcone Complexes. In Stability and Applications of Coordination Compounds; IntechOpen: London, UK, 2020. [Google Scholar]
  36. Santoni, G.; Rehder, D. Structural models for the reduced form of vanadate-dependent peroxidases: Vanadyl complexes with bidentate chiral Schiff base ligands. J. Inorg. Biochem. 2004, 98, 758–764. [Google Scholar] [CrossRef]
  37. Yamaguchi, T.; Watanabe, S.; Matsumura, Y.; Tokuoka, Y.; Yokoyama, A. Oxovanadium complexes with quinoline and pyridinone ligands: Syntheses of the complexes and effect of alkyl chains on their apoptosis-inducing activity in leukemia cells. Bioorg. Med. Chem. 2012, 20, 3058–3064. [Google Scholar] [CrossRef]
  38. Mahmoud, M.A.; Zaitone, S.A.; Ammar, A.M.; Sallam, S.A. Synthesis, spectral, thermal and insulin-enhancing properties of oxovanadium (IV) complexes of metformin Schiff-bases. J. Therm. Anal. Calorim. 2017, 128, 957–969. [Google Scholar] [CrossRef]
  39. Sharfalddin, A.A.; Emwas, A.-H.; Jaremko, M.; Hussien, M.A. Complexation of uranyl (UO2)2+ with bidentate ligands: XRD, spectroscopic, computational, and biological studies. PLoS ONE 2021, 16, e0256186. [Google Scholar] [CrossRef]
  40. Poletto, M.; Zattera, A.J.; Santana, R.M. Effect of natural oils on the thermal stability and degradation kinetics of recycled polypropylene wood flour composites. Polym. Compos. 2014, 35, 1935–1942. [Google Scholar] [CrossRef]
  41. Sharfalddin, A.A.; Emwas, A.H.; Jaremko, M.; Hussien, M.A. Synthesis and Theoretical Calculations of Metal–Antibiotic Chelation with Thiamphenicol; In vitro DNA and HSA Binding, Molecular Docking, and Cytotoxic Studies. New J. Chem. 2021, 45, 9598–9613. [Google Scholar] [CrossRef]
  42. Hardcastle, F.D.; Wachs, I.E. Determination of vanadium-oxygen bond distances and bond orders by Raman spectroscopy. J. Phys. Chem. 1991, 95, 5031–5041. [Google Scholar] [CrossRef]
  43. Thalluri, S.M.; Martinez Suarez, C.; Hussain, M.; Hernandez, S.; Virga, A.; Saracco, G.; Russo, N. Evaluation of the parameters affecting the visible-light-induced photocatalytic activity of monoclinic BiVO4 for water oxidation. Ind. Eng. Chem. Res. 2013, 52, 17414–17418. [Google Scholar] [CrossRef]
  44. Sahu, V.; Sharma, P.; Kumar, A. Impact of Global and Local Reactivity Descriptors on the Hetero-Diels-Alder Reaction of Enaminothione with Various Electrophiles. J. Chil. Chem. Soc. 2014, 59, 2327–2334. [Google Scholar] [CrossRef] [Green Version]
  45. Sirajuddin, M.; Ali, S.; Badshah, A. Drug–DNA interactions and their study by UV–Visible, fluorescence spectroscopies and cyclic voltametry. J. Photochem. Photobiol. B Biol. 2013, 124, 1–19. [Google Scholar] [CrossRef] [PubMed]
  46. Ponya Utthra, P.; Kumaravel, G.; Senthilkumar, R.; Raman, N. Heteroleptic Schiff base complexes containing terpyridine as chemical nucleases and their biological potential: A study of DNA binding and cleaving, antimicrobial and cytotoxic tendencies. Appl. Organomet. Chem. 2017, 31, e3629. [Google Scholar] [CrossRef]
  47. Howsaui, H.B.; Basaleh, A.S.; Abdellattif, M.H.; Hassan, W.M.; Hussien, M.A. Synthesis, Structural Investigations, Molecular Docking, and Anticancer Activity of Some Novel Schiff Bases and Their Uranyl Complexes. Biomolecules 2021, 11, 1138. [Google Scholar] [CrossRef]
  48. Phadte, A.A.; Banerjee, S.; Mate, N.A.; Banerjee, A. Spectroscopic and viscometric determination of DNA-binding modes of some bioactive dibenzodioxins and phenazines. Biochem. Biophys. Rep. 2019, 18, 100629. [Google Scholar] [CrossRef] [PubMed]
  49. Kholod, Y.; Hoag, E.; Muratore, K.; Kosenkov, D. Computer-Aided Drug Discovery: Molecular Docking of Diminazene Ligands to DNA Minor Groove. J. Chem. Educ. 2018, 95, 882–887. [Google Scholar] [CrossRef]
  50. Pratviel, G.; Bernadou, J.; Meunier, B. DNA and RNA cleavage by metal complexes. In Advances in Inorganic Chemistry; Elsevier: Amsterdam, The Netherlands, 1998; Volume 45, pp. 251–312. [Google Scholar]
  51. Bao, Z.; Lai, D.; Shen, P.; Yu, M.; Kumar, R.; Liu, Y.; Chen, Z.; Liang, H. A new samarium (III) complex of liriodenine: Synthesis, crystal structure, antitumor activity, and DNA binding study. Z. Anorg. Allg. Chem. 2019, 645, 570–579. [Google Scholar] [CrossRef]
  52. Sankarganesh, M.; Vijay Solomon, R.; Dhaveethu Raja, J. Platinum complex with pyrimidine- and morpholine-based ligand: Synthesis, spectroscopic, DFT, TDDFT, catalytic reduction, in vitro anticancer, antioxidant, antimicrobial, DNA binding and molecular modeling studies. J. Biomol. Struct. Dyn. 2021, 39, 1055–1067. [Google Scholar] [CrossRef] [PubMed]
  53. Abdel-Rhman, M.H.; Hussien, M.A.; Mahmoud, H.M.; Hosny, N.M. Synthesis, characterization, molecular docking and cytotoxicity studies on N-benzyl-2-isonicotinoylhydrazine-1-carbothioamide and its metal complexes. J. Mol. Struct. 2019, 1196, 417–428. [Google Scholar] [CrossRef]
  54. Alamri, M.A.; Al-Jahdali, M.; Al-Radadi, N.S.; Hussien, M.A. Biological activity evaluation and computational study of novel triazene derivatives containing benzothiazole rings. J. Mol. Struct. 2020, 1227, 129507. [Google Scholar] [CrossRef]
  55. Sinha, A.; Banerjee, K.; Banerjee, A.; Sarkar, A.; Ahir, M.; Adhikary, A.; Chatterjee, M.; Choudhuri, S.K. Induction of apoptosis in human colorectal cancer cell line, HCT-116 by a vanadium-Schiff base complex. Biomed. Pharmacother. 2017, 92, 509–518. [Google Scholar] [CrossRef]
  56. Hussien, M.; Salama, H. Spectrophotometric study of stability constants of Semicarbazone-Cu(II) complex at different temperatures. Pharma Chem. 2016, 8, 44–47. [Google Scholar]
  57. Frisch, M.; Trucks, G.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. Gaussian 09; Revision d. 01; Gaussian, Inc.: Wallingford, CT, USA, 2009; Volume 201. [Google Scholar]
  58. Hay, P.J.; Wadt, W.R. Ab initio effective core potentials for molecular calculations. Potentials for K to Au including the outermost core orbitals. J. Chem. Phys. 1985, 82, 299–310. [Google Scholar] [CrossRef]
  59. MOE (The Molecular Operating Environment). Available online: https://www.chemcomp.com/Products.htm (accessed on 28 December 2021).
  60. Singh, S.P.; Nongalleima, K.; Singh, N.I.; Doley, P.; Singh, C.B.; Singh, T.R.; Sahoo, D. Zerumbone reduces proliferation of HCT116 colon cancer cells by inhibition of TNF-alpha. Sci. Rep. 2018, 8, 4090. [Google Scholar] [CrossRef] [Green Version]
  61. Abdellattif MH, H.M.a.A.E. New approaches of 4-aryl-2-hydrazinothiazole derivatives synthesis, molecular docking, and biological evaluations. Int. J. Pharm. Sci. Res. 2017, 9, 5060–5078. [Google Scholar] [CrossRef]
  62. Al-Khathami, N.D.; Al-Rashdi, K.S.; Babgi, B.A.; Hussien, M.A.; Nadeem Arshad, M.; Eltayeb, N.E.; Elsilk, S.E.; Lasri, J.; Basaleh, A.S.; Al-Jahdali, M. Spectroscopic and biological properties of platinum complexes derived from 2-pyridyl Schiff bases. J. Saudi Chem. Soc. 2019, 23, 903–915. [Google Scholar] [CrossRef]
  63. Mashat, K.H.; Babgi, B.A.; Hussien, M.A.; Nadeem Arshad, M.; Abdellattif, M.H. Synthesis, structures, DNA-binding and anticancer activities of some copper(I)-phosphine complexes. Polyhedron 2019, 158, 164–172. [Google Scholar] [CrossRef]
Figure 1. The suggested structures for oxidovanadium(IV) complexes.
Figure 1. The suggested structures for oxidovanadium(IV) complexes.
Molecules 27 00649 g001
Figure 2. A comparison of the IR spectra of free CTZ ligand and its oxidovanadium(IV) complex in the range 4000–500 cm−1.
Figure 2. A comparison of the IR spectra of free CTZ ligand and its oxidovanadium(IV) complex in the range 4000–500 cm−1.
Molecules 27 00649 g002
Figure 3. Electronic absorption spectrum of free CTZ ligand and its metal complex (inset shows the d-d transition the metal complex).
Figure 3. Electronic absorption spectrum of free CTZ ligand and its metal complex (inset shows the d-d transition the metal complex).
Molecules 27 00649 g003
Figure 4. EPR spectrum of [VO(CTZ)2] 2H2O in DMF at room temperature.
Figure 4. EPR spectrum of [VO(CTZ)2] 2H2O in DMF at room temperature.
Molecules 27 00649 g004
Figure 5. The TGA curve decomposition in the range 25–800 °C of [VO(CTZ)2] 2H2O complex.
Figure 5. The TGA curve decomposition in the range 25–800 °C of [VO(CTZ)2] 2H2O complex.
Molecules 27 00649 g005
Figure 6. Absorption spectra of CTZ and its metal complex in the presence of increasing DNA concentration (arrow indicates changes with increasing DNA concentration).
Figure 6. Absorption spectra of CTZ and its metal complex in the presence of increasing DNA concentration (arrow indicates changes with increasing DNA concentration).
Molecules 27 00649 g006
Figure 7. MLR correlation between Predicted anticancer activities of “IC50” and experimental values of training set.
Figure 7. MLR correlation between Predicted anticancer activities of “IC50” and experimental values of training set.
Molecules 27 00649 g007
Table 1. The effective analytical and physical data of ligands and their oxidovanadium (IV) complexes.
Table 1. The effective analytical and physical data of ligands and their oxidovanadium (IV) complexes.
ComplexM.wt. (g mol−1)ColorM.P.
(°C)
% Found (calc.)Λm
(Ω−1 mol−1 cm2)
CHN
CBZ185.2White11345.155.4115.041
[VO(SO4)(CBZ)] 8H2O493.3Green12517.04
(17.25)
5.31
(5.32)
5.68
(5.60)
7.47
CTZ461.8White14054.625.896.075.23
[VO(CTZ)2] 2H2O878.7Heavy blue13457.41
(57.43)
5.96
(5.97)
6.38
(6.41)
8.85
LOR371.8Yellow22542.002.7111.304.00
[VO(LOR)2] SO4906.6Brown12435.71
(35.73)
2.31
(2.30)
9.61
(9.54)
68.85
SCZ250.1White13042.193.1919.681.3
[VO(SO4)(SCZ)] 7H2O573.8Green20920.93
(30.10)
4.04
(4.13)
9.76
(9.98)
17.50
Table 2. The main bands in IR spectra in the range 4000–500 cm−1 for ligands and their oxidovanadium(IV) complexes.
Table 2. The main bands in IR spectra in the range 4000–500 cm−1 for ligands and their oxidovanadium(IV) complexes.
Compoundν(C=O)ν(C=N)ν(C=S)ν(S=O)ν(V-O)ν(V-N)ν(V-S)νasym, νsym (SO4)
CBZ1677-1271-----
VO(SO4)(CBZ)] 8H2O1677-1279-596422-1310–1035
CTZ1739------ -
[VO(CTZ)2] 2H2O1599, 1409---544 -
LOR16471591-----
[VO(LOR)2] SO416681621--521496
SCZ-1670-1271---
[VO(SO4)(SCZ)] 7H2O-1671-1033589-4241280–1011
Table 3. Electronic spectra data of ligand and their oxidovanadium(IV) complexes and their magnetic moment values.
Table 3. Electronic spectra data of ligand and their oxidovanadium(IV) complexes and their magnetic moment values.
Compoundπ–π*n–π*µeff.
(B.M.)
d-d (broad)
CBZ266- -
[VO(SO4)(CBZ)] 8H2O273-1.74850
CTZ262280 -
[VO(CTZ)2] 2H2O2603311.75900
LOR274394 -
[VO(LOR)2] SO42633111.76850
SCZ262275 -
[VO(SO4)(SCZ)] 7H2O2823731.75870
Table 4. EPR parameter in DMF at room temperature of oxidovanadium(IV) complexes.
Table 4. EPR parameter in DMF at room temperature of oxidovanadium(IV) complexes.
ComplexgA
[VO(SO4)(CBZ)] 8H2O1.99296
[VO(CTZ)2] 2H2O1.980106
[VO(LOR)2] SO41.982105
[VO(SO4)(SCZ)] 7H2O1.984108
Table 5. The optimized geometry and V=O bond length of the metal complexes in the gas phase using DFT method LANL2DZ basis sets.
Table 5. The optimized geometry and V=O bond length of the metal complexes in the gas phase using DFT method LANL2DZ basis sets.
CompoundOptimized Geometry
[VO(SO4)(CBZ)] 8H2O Molecules 27 00649 i001
[VO(CTZ)2] 2H2O Molecules 27 00649 i002
[VO(LOR)2]SO4 Molecules 27 00649 i003
[VO(SO4)(SCZ)] 7H2O Molecules 27 00649 i004
Table 6. Bond lengths and stretching frequencies for V=O Bonds in synthesized oxidovanadium(IV) complexes, experimental, and theoretical values.
Table 6. Bond lengths and stretching frequencies for V=O Bonds in synthesized oxidovanadium(IV) complexes, experimental, and theoretical values.
Complexν(V=O)V=O Radius
ExpDFTExpDFT
[VO(SO4)(CBZ)] 8H2O980 1037 1.6071.599
[VO(CTZ)2] 2H2O963 1068 1.616 1.579
[VO(LOR)2]SO4968 1075 1.6131.576
[VO(SO4)(SCZ)] 7H2O987 1047 1.603 1.596
Table 7. EHOMO, ELUMO, Eg, and the ligand’s quantum parameters and complexes using B3LYP/LANL2DZ.
Table 7. EHOMO, ELUMO, Eg, and the ligand’s quantum parameters and complexes using B3LYP/LANL2DZ.
CompoundHUMOLUMO∆ExησPiσSωΔN Max
CBZ−5.55−1.004.553.282.280.44−3.280.441.142.361.44
VO-CBZ−7.09−3.273.835.181.910.52−5.180.960.962.592.71
SCZ−6.75−2.783.974.771.990.50−4.770.990.992.382.40
VO-SCZ−7.41−3.643.775.531.890.53−5.530.940.942.762.93
LOR−3.96−1.991.972.980.981.02−2.980.490.491.493.03
VO-LOR−3.93−2.741.193.340.591.68−3.340.300.301.675.62
CTZ−6.59−4.681.915.640.951.05−5.640.480.482.825.91
VO-CTZ−10.26−8.841.4210.051.210.83−10.050.610.615.038.31
Table 8. The spectral parameters for the interaction of oxidovanadium(IV) complexes and their ligands with DNA.
Table 8. The spectral parameters for the interaction of oxidovanadium(IV) complexes and their ligands with DNA.
CompoundKb (M−1)λmax Free (nm)λmax Bound (nm)Type of Chromism
CBZ8.33 × 105261255Hyperchromic
[VO(SO4)(CBZ)] 8H2O5.00 × 105265256Hyperchromic
CTZ1.00 × 106263263Hyperchromic
[VO(CTZ)2] 2H2O 1.40 × 106265259Hyperchromic
LOR8.33 × 105395376Hypochromic
[VO(LOR)2] SO41.20 × 106281274Hyperchromic
SCZ1.00 × 106265257Hyperchromic
[VO(SO4)(SCZ)] 7H2O1.25 × 106262256Hyperchromic
Table 9. Docking score and energies of ligands and their oxidovanadium(IV) complexes with 3IG7 protein.
Table 9. Docking score and energies of ligands and their oxidovanadium(IV) complexes with 3IG7 protein.
CompoundSrmsd_RefineE_confE_PlaceE_Refine
CBZ−5.130.95−16.39−50.34−27.37
[VO(SO4)(CBZ)] 8H2O−5.952.88−624.55−49.81−15.08
CTZ−7.601.41129.88−92.44−28.38
[VO(CTZ)2] 2H2O−9.811.86−658.30−79.50−29.23
LOR−6.581.042.82−87.23−36.53
[VO(LOR)2] SO4−6.832.35−358.53−1.13−37.58
SCZ−5.831.93−40.53−57.04−31.67
[VO(SO4)(SCZ)] 7H2O−6.901.16−797.18−79.19−9.94
Table 10. 2D and 3D Docking interaction of the oxidovanadium(IV) complexes with colon cancer protein (PDB code = 3IG7).
Table 10. 2D and 3D Docking interaction of the oxidovanadium(IV) complexes with colon cancer protein (PDB code = 3IG7).
Compound2D Snapshot3D Snapshot
CBZ Molecules 27 00649 i005 Molecules 27 00649 i006
[VO(SO4)(CBZ)] 8H2O Molecules 27 00649 i007 Molecules 27 00649 i008
CTZ Molecules 27 00649 i009 Molecules 27 00649 i010
[VO(CTZ)2] 2H2O Molecules 27 00649 i011 Molecules 27 00649 i012
LOR Molecules 27 00649 i013 Molecules 27 00649 i014
[VO(LOR)2]SO4 Molecules 27 00649 i015 Molecules 27 00649 i016
SCZ Molecules 27 00649 i017 Molecules 27 00649 i018
[VO(SO4)(SCZ)] 7H2O Molecules 27 00649 i019 Molecules 27 00649 i020
Molecules 27 00649 i021
Table 11. The predicted IC50 of all synthesized oxidovanadium(IV) complexes and their free ligands (testing set) against HCT116 cell.
Table 11. The predicted IC50 of all synthesized oxidovanadium(IV) complexes and their free ligands (testing set) against HCT116 cell.
CompoundPredicted IC50 (µM)
CBZ42.43
[VO(SO4)(CBZ)] 8H2O39.87
CTZ12.30
[VO(CTZ)2] 2H2O 1.45
LOR10.26
[VO(LOR)2] SO42.61
SCZ17.62
[VO(SO4)(SCZ)] H2O37.06
Table 12. Anticancer activity of [VO(CTZ)2] 2H2O compared to Cisplatin and the stander division value.
Table 12. Anticancer activity of [VO(CTZ)2] 2H2O compared to Cisplatin and the stander division value.
ComplexIC50 ± SD (µM)
[VO(CTZ)2] 2H2O2.11 ± 0.02
Cisplatin2.13 0.03
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alomari, F.Y.; Sharfalddin, A.A.; Abdellattif, M.H.; Domyati, D.; Basaleh, A.S.; Hussien, M.A. QSAR Modeling, Molecular Docking and Cytotoxic Evaluation for Novel Oxidovanadium(IV) Complexes as Colon Anticancer Agents. Molecules 2022, 27, 649. https://doi.org/10.3390/molecules27030649

AMA Style

Alomari FY, Sharfalddin AA, Abdellattif MH, Domyati D, Basaleh AS, Hussien MA. QSAR Modeling, Molecular Docking and Cytotoxic Evaluation for Novel Oxidovanadium(IV) Complexes as Colon Anticancer Agents. Molecules. 2022; 27(3):649. https://doi.org/10.3390/molecules27030649

Chicago/Turabian Style

Alomari, Fatimah Y., Abeer A. Sharfalddin, Magda H. Abdellattif, Doaa Domyati, Amal S. Basaleh, and Mostafa A. Hussien. 2022. "QSAR Modeling, Molecular Docking and Cytotoxic Evaluation for Novel Oxidovanadium(IV) Complexes as Colon Anticancer Agents" Molecules 27, no. 3: 649. https://doi.org/10.3390/molecules27030649

APA Style

Alomari, F. Y., Sharfalddin, A. A., Abdellattif, M. H., Domyati, D., Basaleh, A. S., & Hussien, M. A. (2022). QSAR Modeling, Molecular Docking and Cytotoxic Evaluation for Novel Oxidovanadium(IV) Complexes as Colon Anticancer Agents. Molecules, 27(3), 649. https://doi.org/10.3390/molecules27030649

Article Metrics

Back to TopTop