In Silico Molecular Docking Analysis of Natural Pyridoacridines as Anticancer Agents
Docking studies are proved to be an essential tool that facilitates the structural diversity of natural products to be harnessed in an organized manner. In this study, pyridoacridines containing natural anticancer pigments were subjected to docking studies using Glide (Schrodinger). Investigations were carried out to find out the potential molecular targets for these selected pigments. The docking was carried out on different cancer macromolecules involved in different cell cycle pathways, that is, CDK-2, CDK-6, Bcl-2, VEGFR-2, IGF-1R kinase, and G-Quadruplexes. CDK-6 was found to be the most suitable anticancer target for the pyridoacridines. In addition, effectiveness of the study was further evaluated by performing docking of known inhibitors against their respective selected macromolecules. However, the results are preliminary and experimental evaluation will be carried out in near future.
Cancer is the deadliest disease and search is in progress to identify the potential anticancer drugs. The important aspect and reason for nonavailability of suitable anticancer drug lies in the fact that cancer displayed its action in human body through different pathways, comprising different cancer macromolecules, and it is very difficult for a single molecule/drug to inhibit all macromolecules at once. Some important examples of cancer macromolecules are B-cell lymphoma 2 (Bcl-2), vascular endothelial growth factor receptor 2 (VEGFR-2), cyclin-dependent protein kinase 6 (CDK-6), CDK-2, IGF-1R kinase (insulin-like growth factor 1 receptor), and G-Quadruplexes .
Some important examples of naturally derived anticancer compounds are Vincristine and Vinblastine from Catharanthus roseus, Paclitaxel from Taxus brevifolia, and Topotecan and Irinotecan from Camptotheca acuminata [2, 3]. Various pigments like varamine, violacein, amphimedine, fascaplysin, monascin, chinikomycin, and so forth were discovered with anticancer activity. Pigments constitute different chemical moieties like pyridine, quinoline, acridine, pyridoacridines, and so forth, which may be responsible for their biological effects [4–7]. Pyridoacridine (Figure 1) is one such moiety present in these pigments which might be responsible for anticancer effects of these pigments.
Different in vitro, in vivo, and computational methods were employed to assess the anticancer potential of drugs or chemicals. Among these methods, docking has been used widely in drug designing for cancer [8, 9]. It is very difficult to identify the role of PCNPs against every single cancer macromolecule. So, it seems appropriate to limit the present study to some selected cancer macromolecules. Role of these macromolecules is well studied by different scientists from time to time  and their inhibition justifies their role in anticancer potential. In the present study, various aspects behind anticancer potential of selected “Pyridoacridine containing natural anticancer pigments (PCNPs)” were assessed. The selected PCNPs were, namely, amphimedine, deoxyamphimedine, neoamphimedine (Xestospongia and Amphimedon species), meridine (Amphicarpa meridiana), and varamine A (Lissoclinum vareau) [10–12]. Additionally, known inhibitors (KIs), that is, SU9516, pyridostatin pentahydrochloride, Palbociclib, OSI-906, HA14-1, and AG879, were docked against their respective macromolecule, that is, CDK-2, G-Quadruplex, CDK-6, IGF-1R kinase, Bcl-2, and VEGFR-2, in order to compare results of the study [13–19].
2. Results and Discussion
In the present study, PCNPs were evaluated through molecular docking studies using Glide (in silico analysis). Initially, the structures of these molecules were generated and energy was minimized. As assessed from literature survey, CDK inhibitors are strongly recommended for cancer treatment as CDKs are not necessary for the cell cycle . Docking studies against CDK-2 revealed that all PCNPs (except neoamphimedine and varamine A) form H-bonds with CDK-2. Meridine interacted with CDK-2 by forming two hydrogen bonds (Leu83, 2.427; 1.928). All PCNPs were suitably docked on CDK-2 (−5.31 to −6.21) and CDK-6 (−2.9 to −9.5) macromolecules with promising Glide score (Figures 2 and 3).
On the other hand, docked PCNPs against G-Quadruplex macromolecule showed Glide score in the range of −1.007 to 0.49; the reason for this poor score may be a consequence of fewer H-bond interactions. Although the Glide score is comparatively poor, visual inspection of docked ligands confirmed their binding affinity towards G-Quadruplex macromolecule (Figure 4). Interaction of PCNPs with VEGFR-2 demonstrated the absence of any H-bonding interaction while in case of Bcl-2, varamine A surprisingly exhibits three H-bonds. Furthermore, PCNPs docked on IGF-1R kinase displayed an H-bonding with varamine A, meridine, and neoamphimedine with Glide score in the range of −6.7 to −4.9.
It was noticed that studied PCNPs have different Glide score, Van der Waals energy, and coulomb energy which is apparent due to the structural difference between these natural pigments of same class. This can be exemplified from the fact that Van der Waals energy value for CDK-6 and CDK-2 inhibitors was among the lowest for almost all PCNPs. Such lower Van der Waals energy value denoted the impact of hydrogen bonding property of these compounds during protein/enzyme interaction . It is interesting to note that PCNPs with an intact pyridine ring, namely, meridine, amphimedine, deoxyamphimedine, and neoamphimedine, have common binding site at Val101 of CDK-6 protein as compared to varamine A. Absence of intact pyridine ring in varamine A might be responsible for its poor Glide score. The Glide score represents the best fit for a ligand in the active site of the target macromolecule and it is summarized as follows: Glide score for PCNPs was in the range of −7.16 to −9.53 (CDK-6); −5.5 to −6.78 (IGF-1R kinase); −5.31 to −6.21 (CDK-2); −4.38 to −4.95 (Bcl-2); and −1 to 0.49 (G-Quadruplex). Secondly, a closer look at each PCNP structure under study discloses that substituting methyl group on nitrogen with hydrogen atom (in case of meridine) makes meridine more suitable to form hydrogen bond with the target macromolecule as methyl group is more bulky and less polar as compared to H-atoms.
Figure 4 represents meridine docking on G-Quadruplex macromolecule. These points indicated meridine effectiveness to be considered as prominent anticancer PCNP molecule, where potential attachment positions of meridine for testing substituents could improve the anticancer activity by increasing the corresponding affinity for the target. On the other hand, each PCNP under investigation forms H-bonds with CDK-6 macromolecule which makes CDK-6 the most suitable anticancer target for the pyridoacridines. Therefore, in the first place, all PCNPs might have an important role in CDK-6 inhibition .
Different bonding interactions are shown in Table 1 while from Table 2 it can be easily anticipated that meridine was best docked to CDK-6 with Glide energy of −44.78, Van der Waals energy of −41.76, and coulomb energy of −3.02. In a study conducted by Lu and coworkers, carbonyl group at 4th position of Fisetin (a CDK-6 inhibitor) formed H-bond with Val101; this is in accordance with present results where different PCNPs formed H-bond with Val101 of CDK-6 macromolecule while varamine A formed H-bond with Lys43, that is, a catalytic residue that undergoes a large conformational change during activation of CDK-6 by cyclin binding so it can be ascertained that H-bonding of varamine A with Lys43 residue might activate CDK-6 macromolecule .
2.1. Effectiveness of the Study
Effectiveness of the study was further evaluated by performing docking of KIs against their respective macromolecule which can be assessed by comparing the interactions of “PCNPs-macromolecule complex” with “KIs-macromolecule complex.” Interacting residues in both complexes are similar (like meridine-CDK-6/Palbociclib-CDK-6 have “Val101” while meridine-CDK-2/SU9516-CDK-2 have “Leu83”) which might be a reason for their comparative Glide score while in other cases different interacting residues of KIs and PCNPs could be a reason for their contradictory Glide score (Figure 5 and Tables 1 and 3).
E-pharmacophores of different PCNPs were generated in accordance with the docked macromolecule. Important pharmacophoric sites (Table 4 and Figure 6) in each PCNP were identified by this method to highlight the various structural aspects of PCNP responsible for docking on a particular macromolecule. For example, in case of meridine docked on CDK-6 protein, numbers of pharmacophoric sites identified by e-pharmacophore were 6 while in case of meridine docked on IGF-1R kinase, the number of sites changes to 5 which signifies the role of e-pharmacophore generation.
Each pharmacophore displayed the presence of pharmacophoric sites, that is, aromatic ring (R), hydrogen bond acceptor (A), hydrogen bond donor (D), and hydrophobic group (H), as per their activity/affinity towards that particular macromolecule. Similar results were obtained on generating e-pharmacophores of each docked PCNP on a specific macromolecule.
2.2. Rational Modification of Pyridoacridine Moiety
All selected PCNPs were evaluated on their specific features and were subsequently ranked as per their Glide score. From docking results, some important structural aspects which could be helpful in designing new pyridoacridine derivatives are summarized as follows:(1)Presence of an intact pyridine ring is necessary in a PCNP molecule, as varamine A with open pyridine ring displayed poor Glide score when docked on CDK-6 molecule.(2)Poor Glide score was observed with substitution on nitrogen heteroatom in pyridinone ring. As meridine has no substitution at nitrogen heteroatom, its Glide score is better (Glide score: −9.53) compared to other PCNPs (except varamine with open pyridinone ring) having substitution at nitrogen heteroatom (Glide score range: −7.66 to −2.95). So, order of affinity towards CDK-6 macromolecule is as follows:(a)Meridine > neoamphimedine > deoxyamphimedine > amphimedine > varamine A(3)Meridine exhibits maximum number of pharmacophoric sites compared to any other PCNP when docked on any cancer macromolecule except Bcl-2. Similarly, varamine A exhibited maximum number of pharmacophoric sites when docked on G-Quadruplex.(4)Presence of either hydrogen bond acceptors or donors seems to be an important structural aspect. Careful substitution of electron withdrawing or electron donating groups in pyridine ring may improve the affinity/activity of the ligand.(5)As carbonyl group makes H-bonding in most of the cases, its presence seems to be essential in case of PCNP lead molecule.(6)It was also observed that more H-bonding between ligand-macromolecule complexes makes the complex more stable, so use of functional groups which are known for their H-bonding could improve the affinity and as a result could improve the anticancer action.
Based on these points, it can be asserted that appropriate modifications at meridine nucleus may provide novel anticancer derivatives.
3. Materials and Methods
Docking studies of PCNP ligands and KIs of selected macromolecules (Figure 7) were performed using Glide, version 5.6, 2010 , against the selected cancer macromolecules, where alteration of expression for each macromolecule corresponds to a different anticancer mechanism.
3.1. Ligand Preparation
It is the initial step for molecular docking studies. LigPrep module, version 2.4, 2010 , was used for geometrical refining of chemical structures (drawn in Maestro module) of PCNP and KIs ligands of selected macromolecules. LigPrep is intended to set up premium 3D structures with accurate chiralities. Original states of ionization were retained; tautomers and conformations were generated by the Monte Carlo method as implemented in MacroModel version 9.8, 2010 , using OPLS-2005 force field. The generated conformers were subsequently minimized using truncated Newton conjugate gradient (TNCG) minimization up to 500 iterations. The conformers with an energy difference of 30 kcal/mol as compared to the global energy minimum conformer were retained. The conformational searches were carried out for aqueous solution using the generalized born/solvent accessible surface (GB/SA) continuum solvation model [25, 26].
3.2. Protein Preparation
Protein preparation wizard of Maestro software was used for protein preparation. The protein structures, namely, 1DI8 (CDK-2), 1L1H (G-Quadruplex), 1XO2 (CDK-6), 2OJ9 (IGF-1R kinase), 2O2F (Bcl-2), and 2OH4 (VEGFR-2), were taken from Protein Data Bank [20, 27–32]. The selected chains were edited for missing hydrogens and for assigning proper bond orders. The H-bonds were optimized using sample orientations. All the polar hydrogens were displayed. Finally, the protein structure was minimized to the default Root Mean Square Deviation (RMSD) value of 0.30.
3.3. Receptor Grid Generation
From the defined receptor, the cocrystallized ligand was separated from the active site of receptor chain. The atoms were of size equal to Van der Waals radii of 1.0 Ǻ while the partial atomic charge was less than 0.25 defaults. The active site represents an enclosing box at the centroid of the workspace ligand. Following this protocol, a grid centered on the ligand was generated using the default Glide settings. All ligands were docked into this grid structure.
3.4. Molecular Docking Analysis
On a defined receptor grid, flexible docking was performed using the extra precision (XP) feature of Glide module, version 5.6, 2010 . The constraints to defined ligand-receptor interactions were not set. The structure output format was set to pose viewer file so as to view the output of the resulting docking studies from pose viewer.
3.5. Generation of E-Pharmacophore
E-pharmacophore was generated for each of the KIs as well as PCNP docked on different cancer macromolecules. The process started automatically with refined ligands where pharmacophoric sites were generated with PHASE, Schrodinger, LLC, New York . PHASE uses number of features like positive ionizable group (P), hydrogen bond acceptor (A), aromatic ring (R), hydrogen bond donor (D), negative ionizable group (N), and hydrophobe (H). An energetic value was assigned to each pharmacophoric site on the basis of Glide XP descriptors.
3.6. Selection of the Best-Scored Pose
The best docking poses for the PCNPs were selected primarily by considering the docking scores but values of different energies, number of H-bonds, and visual inspection of all docking poses in Maestro (Schrodinger, USA) were also taken into account. Interaction energy between protein and ligand can be related to binding affinities. Different criteria were laid down to select best-docked structure for each ligand. Then, rankings were derived by directly using the Glide GScore.
From the study, two questions arise: first, which macromolecule could be responsible for anticancer mechanism of PCNPs? And second one is, which PCNP could be considered as a “lead molecule” of the future? The answer to the first question is, since almost all PCNP ligands have good Glide score except varamine A, CDK-6 can be considered as the most suitable macromolecule through which PCNPs exhibit their anticancer mechanism while in second case, it was “meridine” which bears the best Glide score for each macromolecule except varamine A. Glide score of a particular ligand describes its effectiveness against the desired target. Docking results clearly judged meridine (best Glide score, maximum H-bonds, and promising Van der Waals and coulomb energy between the ligand and receptor) to be a promising PCNP and it could be considered as potential “lead molecule” for the development of novel anticancer drugs. Docking results of KIs against selected macromolecules further validate our assertion behind suitability of “meridine as a lead molecule,” as meridine and Palbociclib both interact with the same residue (i.e., Val101) of CDK-6 and exhibit comparative Glide scores of −9.533 and −8.066, respectively.
From the docking results, it is quite clear that nitrogen atom substitution of the pyridine ring of pyridoacridines results in lesser Glide score which corresponds to poor binding interaction with the target macromolecule. Furthermore, the results indicate higher affinity of pyridoacridines towards CDK-6 macromolecule, which might be considered as the potential mechanism for pyridoacridines. On the other hand, among the studied pyridoacridines, meridine comes with highest Glide score against each target macromolecule which makes it an appropriate “lead compound.” Some distinguished features of meridine are the following: (a) nitrogen atom of the pyridine group, (b) the presence of two carbonyl groups, and (c) presence of nitrogen (donor) alongside carbonyl group (acceptor) might decrease electron density in the region; all these features promote H-bonding with the active site of the target macromolecule which results in the high Glide score. Presence of bulkier methyl group on nitrogen atom of pyridine ring in case of amphimedine, deoxyamphimedine, and neoamphimedine with absence of pyridine ring in varamine A disfavors hydrogen bonding with the corresponding target macromolecule. However the results are preliminary and certainly need experimental confirmation, which will be conducted in near future via molecular biology studies but considering all these structural aspects and Glide score, meridine may possibly be a first choice as a PCNP lead molecule which could be exploited to design anticancer agents of future.
Binding models of PCNPs and KIs against different cancer macromolecules are shown in Figures S1 and S2, respectively, in Supplementary Material available online at http://dx.doi.org/10.1155/2016/5409387.
The authors declare that there is no conflict of interests regarding the publication of this article.
Mr. Vikas Sharma is thankful to All India Council for Technical Education, India, for providing National Doctoral fellowship for Ph.D. research work (vide letter no. F. No. 14/AICTE/RIFD/NDF(Policy-I)/01/2012-13). Director of Institute of Pharmaceutical Sciences, Kurukshetra University, Kurukshetra, Haryana 136 119, India, is duly acknowledged for providing necessary research facilities. The authors sincerely acknowledge Schrodinger LLC., Portland, USA (https://www.schrodinger.com/), for providing software for research purpose.
Binding models displays the interactions between PCNPs/KIs against different selected cancer macromolecules.
V. M. Srinivasan and C. S. Devi, “Novel anticancer compounds from marine actinomycetes: a review,” Journal of Pharmacy Research, vol. 4, pp. 1285–1287, 2011.View at: Google Scholar
R. Ali, Z. Mirza, G. M. D. Ashraf et al., “New anticancer agents: recent developments in tumor therapy,” Anticancer Research, vol. 32, no. 7, pp. 2999–3005, 2012.View at: Google Scholar
S. Tabassum, M. Zaki, M. Afzal, and F. Arjmand, “Synthesis and characterization of Cu(II)-based anticancer chemotherapeutic agent targeting topoisomerase Iα: in vitro DNA binding, pBR322 cleavage, molecular docking studies and cytotoxicity against human cancer cell lines,” European Journal of Medicinal Chemistry, vol. 74, pp. 509–523, 2014.View at: Publisher Site | Google Scholar
S. N. e Zahra, N. A. Khattak, and A. Mir, “Comparative modeling and docking studies of p16ink4/cyclin D1/Rb pathway genes in lung cancer revealed functionally interactive residue of RB1 and its functional partner E2F1,” Theoretical Biology & Medical Modelling, vol. 10, pp. 1–9, 2013.View at: Publisher Site | Google Scholar
F. J. Schmitz, S. K. Agarwal, S. P. Gunasekera, P. G. Schmidt, and J. N. Shoolery, “Amphimedine, new aromatic alkaloid from a pacific sponge, Amphimedon sp. Carbon connectivity determination from natural abundance carbon-13-carbon-13 coupling constants,” Journal of the American Chemical Society, vol. 105, no. 14, pp. 4835–4836, 1983.View at: Publisher Site | Google Scholar
J.-L. Wang, D. Liu, Z.-J. Zhang et al., “Structure-based discovery of an organic compound that binds Bcl-2 protein and induces apoptosis of tumor cells,” Proceedings of the National Academy of Sciences of the United States of America, vol. 97, no. 13, pp. 7124–7129, 2000.View at: Publisher Site | Google Scholar
Glide, Version 5.6, Schrodinger. LLC, New York, NY, USA, 2010.
LigPrep, Version 2.4, Schrodinger. LLC, New York, NY, USA, 2010.
MacroModel, Version 9.8, Schrodinger. LLC, New York, NY, USA, 2010.
U. Kaushik, V. Sharma, and V. Kumar, “Computation of pharmacophore models for the prediction of mitogen-activated protein kinase activated protein kinase-2 inhibitory activity of pyrrolopyridines,” Medicinal Chemistry Research, vol. 21, no. 11, pp. 3777–3784, 2012.View at: Publisher Site | Google Scholar
L. Shewchuk, A. Hassell, B. Wisely et al., “Binding mode of the 4-anilinoquinazoline class of protein kinase inhibitor: X-ray crystallographic studies of 4-anilinoquinazolines bound to cyclin-dependent kinase 2 and p38 kinase,” Journal of Medicinal Chemistry, vol. 43, no. 1, pp. 133–138, 2000.View at: Publisher Site | Google Scholar
U. Velaparthi, M. Wittman, P. Liu et al., “Discovery and initial SAR of 3-(1H-benzo[d]imidazol-2-yl)pyridin-2(1H)-ones as inhibitors of insulin-like growth factor 1-receptor (IGF-1R),” Bioorganic and Medicinal Chemistry Letters, vol. 17, no. 8, pp. 2317–2321, 2007.View at: Publisher Site | Google Scholar
PHASE, Version 3.0, Schrodinger, LLC, NY, 2008.