Abstract

The anticancer study of nitrogen-chelating ligands can be of tremendous help in choosing ligands for the anticancer metal complexes design especially with ruthenium(II). The inhibitory anticancer activities of some nitrogen-chelating ligands containing bis-pyrazole, bipyridine, and phenanthroline were studied using experimental screening against cancer cell and theoretical docking methods. In vitro anticancer activities showed compound 11 as the most promising inhibitor, and the computational docking further indicates its strong inhibitory activities towards some cancer-related receptors. Among the twenty-one modelled ligands, pyrazole-based compounds 7, 11, and 15 are the most promising inhibitors against the selected receptors followed by 18 and 21 which are derivatives of pyridine and phenanthroline, respectively. The presence of the carboxylic unit in the top five ligands that displayed stronger inhibitory activities against the selected receptors is an indication that the formation of noncovalent interactions such as hydrogen bonding and a strong electron-withdrawing group in these compounds are very important for their receptor interactions. The thermodynamic properties, the polarizabilities, and the LUMO energy of the compounds are in the same patterns as the observed inhibitory activities.

1. Introduction

Nitrogen-chelating ligands such as bis-pyrazole (pz), bipyridine (bpyr), and phenanthroline (phn) derivatives are being used in the design of several metal complexes for ranges of applications from biological to nanostructured materials, especially the development of anticancer complexes of ruthenium [119]. Other studies have shown that these nitrogen-chelating compounds play significant roles in the experimentally observed anticancer activities [20] that might be ascribed to the electronic interactions between the metal centre and the π-electrons in rings [2123]. However, despite the use of these nitrogen-chelating ligands in the design of anticancer metal complexes, there has been very few or no attention on the anticancer activities of these ligands individually. To this end, we have selected some derivatives of common nitrogen-chelating ligands as shown in Figure 1 to study their individual anticancer activities with particular focus on their interactions with cancer-related receptors using various docking methods. We also report the experimental in vitro anticancer activities of some of these ligands.

The selected receptors used for the docking studies are carbonic anhydrase II (CA-II), cathepsins B (Cat B) [24], two different DNAs (DNA-1 [25] and DNA-2), DNA gyrase (Gyrase) [26], histone deacetylase7 (HDAC7) [27], histone protein in the nucleosome core particle (HIS) [28], BRAF kinase (Kinase) [29], recombinant human albumin (rHA) [30], ribonucleotide reductase (RNR) [31], topoisomerase II (Top II) [32], thioredoxin reductase (TrxR) [33], and thymidylate synthase (TS) [34]. These receptors play significant roles in cancer growth and are thus a unique target in cancer therapy. For instance, rHA plays a significant role in the pharmacokinetic availability, bioavailability, and toxicology [35] and helps either in delivery of metal-based anticancer drugs to their cellular targets or in deactivating them even before reaching the target(s) [36]. The RNR is responsible for the synthesis of DNA from the corresponding building blocks of RNA [37]. Top II plays a key role in relaxing supercoiled DNA for replication and transcription in the absence of inhibitors [38], while the presence of inhibitors forms a stable complex with the enzyme and keeps it from DNA cleavage [39]. Thioredoxin reductase (TrxR) regulates the cellular reduction/oxidation (redox) status [40, 41]. Thymidylate synthase (TS) is a critical enzyme in maintaining a balanced supply of deoxynucleotides required for DNA synthesis and repair [42]. It is a target of chemotherapy to test the vulnerability of cancer cells to the inhibition of TMP synthesis [37]. Gyrase was considered to establish possible dual roles for the ligands as potential antibacterial agents.

2. Synthesis and Structural Elucidation of the Ligands

The synthesis and careful structural elucidation of these ligands have been reported in our previous works [4345].

3. Experimental Methods

The in vitro anticancer activities of the ligands against the cancer cell line HT29 and normal cell line KMST were examined using the MTT colorimetric assay. All the ligands before their docking to the receptors were first optimized with DFT functional PBEPBE [46] and the basis set 6-31 + G(d,p) for all atoms using the Gaussian 09 package [47]. The docking analyses were carried out using Molegro [48] and Vina and AutoDock [49] packages. The docking of each ligand against the receptors was done five times: twice in Molegro first using the quantum Mulliken atomic charges for all atoms of each ligand (subsequently referred to as Molegro-QC) and second using the predicted atomic charges from the Molegro package (referred to as Molegro). One time in the Vina package using the predicted atomic charges from AutoDock tools and twice in AutoDock using both the QM charges (AutoDock-QC) and package-predicted atomic charges (AutoDock) has done in Molegro.

The default parameters were used in the Vina docking package but with little modifications in Molegro and AutoDock dockings. The scoring function used in Molegro was MolDock because it takes care of the hydrogen bonding, intermolecular protein ligand, and intramolecular ligand interactions and has been successfully applying for molecular docking [50]. The maximum interaction was set to 2500 instead of the default value of 1500, and the population number was increased from the default value of 50 to 100. In using AutoDock, the number of grid points in x-, y-, and z-axes was set to 60 × 60 × 60 with each point separated by 0.375 Å. The Lamarckian genetic algorithm was chosen based on its efficiency and reliability in comparison with others like simulated annealing (SA) and generic genetic algorithm (GA) methods in AutoDock [51, 52]. The maximum number of energy evaluations was set to 2,500,000 for each of the 20 independent runs, a maximum number of 27,000 GA operations were generated on a single population of 100 individuals, and step sizes of 2 Å for translation and 50° for rotation were chosen.

All the graphical representations of the docking results are prepared using the package Chimera [53].

4. Results and Discussion

Twenty-one ligands were modelled which comprise fifteen models of bis-pyrazole, three models of bipyridine (bpyr), and three models of phenanthroline (phn) (Figure 1). Six of the modelled compounds were screened to determine their in vitro anticancer activities, and their results are presented in Table 1. The in vitro activities of these ligands were tested against the cancer cell line HT29 and the normal cell line KMST using the MTT colorimetric assay. The results clearly showed that the six ligands pose no threat to the normal cell line as their inhibitory activities (IC50) in KMST are found to be greater than 50 μm as shown in Table 1. Ligands 11 and 14 are the most potent (<6.25 mm) against the cancer cell line HT29, while ligand 12 is the least active. Also, the same ligand 11 was found among the three best inhibitors of the selected receptors from the docking results (Table 2).

4.1. The Binding Site Predictions of the Compound from the Docking Methods

The results obtained from the different docking methods and packages: Molegro-QC (grey), Molegro (green), Vina (brown), AutoDock-QC (yellow), and AutoDock (cyan), are shown in Figure 2. The features of the interactions of the ligands with the receptors using different docking methods are similar. In most of the receptors, the same binding sites were located by docking methods, and very similar conformational orientations of the ligands in the binding sites were predicted with some found to overlap each other (Figure 2). Ligands 11 and 15 have similar binding orientation in CA-II, Cat B, and HDAC7 based on the results obtained from the binding site interactions. Compound 11 is of interest because it shows promising anticancer activities in vitro, and ligand 15 is predicted as the best inhibitor of many of the receptors according to the results obtained from Molegro-QC and Molegro. The interaction of 11 with CA-II shows that all the docking methods gave similar orientation for the compounds except Molegro-QC (grey) that is slightly different in orientation (Figure 2).

The orientation obtained from Molegro-QC (grey) gives the best interacting energy compared to others (Table 2). Also in the interaction of 15 with CA-II, all of the methods locate the same binding site and give similar orientation of the ligand except Vina (brown) which locates a different binding site compared to others. The Molegro-QC and Molegro orientations of 15 are superimposed and likewise AutoDock and AutoDock-QC. All the five methods except Vina also locate the same binding site for the interaction of 11 with Cat B, and AutoDock prediction was found to be superimposed with that of Molegro-QC and Molegro. The results obtained from AutoDock predict different binding sites for the interaction of 15 with Cat B, while all other methods predict the same site for its binding. Another receptor of interest is HDAC7 in which strong inhibitory activities are displayed by 11 and 15. HDAC7 was predicted as one of the most targeted receptors for many of the ligands according to the results obtained from Molegro-QC and Molegro. Vina locates a completely different binding site for the interaction of 11 with HDAC7 when compared to the rest of the methods. The binding site orientation of 11 from all the methods besides Vina is very similar to that of Molegro-QC completely superimposed with that of AutoDock. All methods predicted a similar binding site orientation of the ligand 15, of which the Molegro-QC and Molegro orientations are found to be superimposed and in close orientation with those of Vina and AutoDock-QC, which are different from those of AutoDock. There is a significant difference between the results from AutoDock-QC and AutoDock methods in terms of the features of the ligand interactions with receptors just because of changes in the accuracy of the atomic charges. The correlation of the results of AutoDock-QC and AutoDock ranges from −0.32 to 0.46, while that of Molegro-QC and Molegro ranges from 0.96 to 1.00. It is obvious that the accuracy of atomic charges plays a very strong role in the determination of the ligand interaction with the receptors especially when using AutoDock. The results show that the reason for the inconsistency from Vina, AutoDock-DC, and AutoDock could be ascribed to different binding sites which were predicted for many of the ligand interactions with the receptors, while high consistency and similarity were obtained in Molegro-QC and Molegro because they both predicted the same binding site as cocrystallized inhibitors of the receptors. In addition, the differences in the ligand-receptor inhibitory energies obtained from AutoDock methods for many of the ligands are within the standard error margin of ∼2.177 kcal/mol [43, 44] that makes the order obtained from AutoDock unreliable.

4.2. The Inhibitory Activities

The ligands 7, 11, and 15 have the best inhibitory activities towards many of the receptors according to the results obtained from Molegro-QC and Molegro (Table 2 and Figure 3). Besides the first three bis-pyrazole ligands, next promising ligands in interaction with the selected receptors are 18 and 21 which are bpyr and phn compounds. The common feature possessed by all the promising ligands is carboxylic acid moieties which are predicted to enhance the noncovalent interactions such as hydrogen bonding and play the role of electron withdrawing. Since the results from the two methods Molegro-QC and Molegro are highly correlated (0.96 to 1.00), Molegro-QC is then used for the analysis of the receptor interaction of ligands 7, 11, and 15. Although in many of the ligand interactions with the receptors, AutoDock and Vina gave similar features in their interaction, the ranking of the ligands according to their inhibitory interaction is vague. It has been pointed out that, to have a good ranking in AutoDock, the error values of ∼2.177 kcal/mol [43, 44] have to be considered, and this makes AutoDock unsuitable for ranking of these ligands because the differences in their interacting energy are lower.

The three best inhibitors among all the modelled compounds are the bis-pyrazole-based compounds 7, 11, and 15 (Figure 3). The binding site interactions of these three ligands with each of the receptors are shown in Figure 4, while the interactions of each of the compounds 7, 11, and 15 with each of the receptors are shown in supplementary Figure S1. In most of the receptors, the three compounds 7 (grey), 11 (green), and 15 (yellow) are found to be clustering around the same point and sharing common angles of interaction with the receptors CA-II, Cat B, DNA-2, Gyrase, HDAC7, HIS, Kinase, rHA, RNR, Top II, and TS (Figure 4). The few exceptions to these common interaction features are found in the receptors DNA-1 and TrxR (Figure 4). A different binding site was predicted for the interaction of 11 with DNA-1 where 7 and 15 have inner groove interaction but the interaction of 11 is out of groove. In addition, the binding site predicted for 11 is separated from the one predicted for 7 and 15 in the receptor TrxR. In all of the receptor interactions especially for 7 and 15, the predicted binding sites are the same as the binding site for their cocrystallized inhibitors except in the receptor DNA-2 where the three ligands bind to a groove outside of the binding surface of the cocrystallized cis-platin.

The receptor interactions of the three most promising compounds 7, 11, and 15 follow a different order based on their receptor preferences. The order of their inhibitory activities in HDAC7, Kinase, rHA, RNR, TrxR, and TS is 15 > 7 > 11, while they follow the order of 7 > 5 > 11 in their interaction with CA-II, Cat B, and Top. The interaction of 11 was found to be better than that of 7 with the DNA-1, Gyrase, and HIS in the order of 15 > 11 > 7, and also, the inhibitory activity of 11 against DNA-2 is found to be higher than that of 7 and 15 in the order of 11 > 15 > 7. The presence of more than one carboxylic unit in 7 and 15 compared to 11 results in a greater number of HB interactions (Table 3), which could also be responsible for their greater inhibitory activities compared to 11 in many of the receptors. The number of the HB interactions of the ligand 11 in most of the receptors ranges from 0 to 2 except in the receptors Top II, TrxR, and TS where their binding site residues support more hydrogen bond interactions with the ligands (Table 3 and Figure S1).

4.3. Molecular Properties of the Compound

The molecular properties of the modelled compounds are shown in Table 4. Compounds 13 and 14 have the highest hyperpolarizability values, while 13 and 21 have the lowest band gap, but 14 and 21 have the highest dipole. The hyperpolarizability of the ligands follows the order of 13 > 14 > 12 > 21. These four ligands can be useful as building blocks for nonlinear optical materials, but only the ligand 21 appears among the best five inhibitors obtained from the docking results. In order to study the possible effects of changes in the molecular properties on the docking interactions with the receptors, the molecular properties of the compound correlated with the docking results obtained from the Molegro-QC method.

All the thermodynamic energy, zero energy, thermal energy, enthalpy, and the Gibbs free energy of the compound have a correlation range of 0.92 to 0.98 with the docking results. In addition, the CV and entropy highly correlate with the inhibitory values of the compounds in the ranges −0.90 to −0.96 for CV and −0.92 to −0.97 for entropy. Besides the thermodynamic properties, other properties that show high correlation with receptors’ inhibitory activities are their polarizabilities (range from 0.74 to 0.90) and polarizabilities W (range from 0.78 to 0.93). Other computed properties like LUMO (0.46 to 0.64), gap (0.32 to 0.60), Hyp (−0.20 to −0.40), and dipole (−0.22 to −0.37) have very low correlation (as shown in the parenthesis). The HOMO has the poorest correlation with the binding activities of the ligands which range from 0.00 to 0.18. The better correlation of the LUMO with the receptor interactions further supports the hypothesis that the lower the LUMO energy of inhibitors the easier the overlapping of it with the HOMO of the DNA [54]. However, other factors like the hydrophobicity and electronic effect of the ligands play significant roles in their inhibitory activities [54].

5. Conclusion

The inhibitory potentials of nitrogen donor ligands in their interaction with twelve cancer-related receptors were studied using docking methods. Twenty-one derivatives of nitrogen-chelating ligands consisting of fifteen pyrazole, three bipyridine, and three phenanthroline derivatives were modelled. In addition, the experimental in vitro anticancer activities of the six ligands 5, 6, 9, 11, 12, and 14 against the cancer cell line HT29 are discussed. The result of the in vitro study shows that the ligands 11 and 14 are the most active ones compared to other ligands. Among the two, only the ligand 11 showed promising inhibitory activities towards the selected receptors. In many of the receptors, different docking methods locate similar binding sites and similar ligand orientations. The order of the ligand interactions with the receptors using Vina, AutoDock-DC, and AutoDock is inconsistent because they locate different binding sites in many cases and their energy difference is also within the range of the standard error of AutoDock (∼2.177 kcal/mol), while the order from Molegro-QC and Molegro is highly consistent and similar.

The results from the Molegro-QC and Molegro showed that 7, 11, and 15 have better inhibitory activities towards many of the receptors and have similar angles of interaction with CA-II, Cat B, DNA-2, Gyrase, HDAC7, HIS, Kinase, rHA, RNR, Top II, and TS. This is an evidence that they have common residue interactions at the binding sites of each of the receptors. The receptor interactions of the three most promising ligands 7, 11, and 15 follow a different order based on their receptor preference. Apart from these three, 18 and 21 with bipyridine and phenanthroline moieties also show promising interaction with many of the receptors. The most common feature of the best five inhibitors is the carboxylic unit, which indicates that the carboxylic units enhance the binding site interaction of the ligands through formation of stronger HB interactions with the receptor residues. The thermodynamic properties have a high correlation with the docking results of the ligands’ inhibitory activities. Other molecular properties, which have a high correlation with the ligands’ inhibitory activities are polarizabilities. The LUMO of the ligands shows a good correlation with their receptor-binding interaction which further supports the hypothesis that the lower the LUMO energy of inhibitors the easier the overlapping of it with the HOMO of the DNA.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the financial support of Sasol and ESKOM TESP, South Africa. The authors also gratefully acknowledge the CHPC of the Republic of South Africa for providing the G09 and computing facilities.

Supplementary Materials

Figure S1: the binding site interaction of each of the three ligands 7 (grey), 11 (green), and 15 (yellow) (horizontally) with the receptors. (Supplementary Materials)