Molecular Docking and Drug-Likeness for the Identification
Total Page:16
File Type:pdf, Size:1020Kb
Biomedical & Pharmacology Journal, September 2018. Vol. 11(3), p. 1301-1307 Molecular Docking and Drug-Likeness for the Identification of Inhibitory Action of Acetogenins from Annona muricata as Potential Anticancer against Hypoxia Inducible Factor 1 Alpha Supri I. Handayani1, Rahmiati1, Lisnawati Rahmadi1, Rosmalena2 and Vivitri D. Prasasty3* 1Department of Anatomical Pathology, Faculty of Medicine, University of Indonesia, Jalan Salemba Raya 6, Jakarta 10430, Indonesia. 2Department of Medical Chemistry, Faculty of Medicine, University of Indonesia, Jalan Salemba Raya 6, Jakarta 10430, Indonesia. 3Faculty of Biotechnology, Atma Jaya Catholic University of Indonesia Jalan Jenderal Sudirman 51, Jakarta 12930, Indonesia. *Corresponding author E-mail: [email protected] http://dx.doi.org/10.13005/bpj/1492 (Received: 03 August 2018; accepted: 30 August 2018) Hypoxia inducible factor 1 alpha (HIF-1a) regulates cell growth and differentiation which is implicated in human cancers. HIF-1a activates its cascade carcinogenesis mechanism in cancer cells. It is well-understood that signaling is initiated by HIF-1± receptor. Overexpression of HIF-1a is associated with several different human cancers, including breast cancer, lung cancer and colon cancer. Thus, HIF-1a becomes potential target of therapeutic approach in developing HIF-1a inhibitors. The aim of this research is to investigate potential inhibitors which are known as Acetogenins (AGEs) isolated from Annona muricata against HIF-1a. In order to achieve this goal, chemical structures of all compounds were retrieved from PubChem database. Molecular docking was performed by AutoDock Vina program and the resulting binding modes were analyzed with AutoDock Tools program. Among all the compounds, murihexocin A showed the best binding modes compared to other two inhibitors based on the lowest binding energies (LBE = -7.9 kcal/mol) as high as gefitinib. This was indicating that murihexocin A has favorable interaction with the essential amino acid residues at catalytic site of HIF-1a. Drug-likeness calculation of AGEs were also performed. These in silico results could be beneficial as a compound model for further studies in-vitro and in-vivo. Keywords: HIF-1a, Annona muricata, inhibitor, acetogenins, molecular docking, drug-likeness. Hypoxia inducible factor-1 (HIF-1) is a sequence 52 -RCGTG-32 (which R is purine) at transcription factor that controls the expression of the response elements of hipoxia to target genes gene that involved in tumorigenesis and metastases 4. The transcription process of various genes are of malignancies 1, 2. The HIF-1a level could be activated by HIF-1, including glycolytic enzymes, enhanced within breast oncogénesis, and it closely gluconeogenesis, mediating glucose transporters, related with other tumor biomarkers 3. HIF-1a growth factors, high-energy phosphate metabolism, plays an important role in binding the consensus heme metabolism, iron transport, erythropoiesis, This is an Open Access article licensed under a Creative Commons license: Attribution 4.0 International (CC-BY). Published by Oriental Scientific Publishing Company © 2018 1302 HANDAYANI et al., Biomed. & Pharmacol. J, Vol. 11(3), 1301-1307 (2018) synthesis of nitric oxide, and regulation of with HIF-1a through molecular docking methods. vasomotor. Therefore, HIF-1 possibly promotes Computational methodologies have become a the tumor cell viability in hypoxic circumstances crucial component in drug discovery program, 5, 6. which involves identification to lead optimization. Hypoxia induces tumor cell proliferation, Molecular modeling is one of the methodologies metastasis, and the rate of cell apoptosis 7. primarily used as hit identification tool when only Moreover, HIF-1 is considered as a starting point of structure of target and its active or binding site are angiogenic process in tumor cells by transcription available 17. Docking method is an energy-based activation of cancer-related gene, such as vascular scoring function which identifies the energetically endothelial growth factor (VEGF) gene 8. The level most favorable ligand conformation that binds to of HIF-1a escalates the pathological stage which the target 18. is higher in poorly differentiated lesions than in well-differentiated lesions 3. The enhanced levels of METHOD HIF-1a are strongly related with high proliferation and enhanced ER as well as VEGF expressions 9. Protein structure preparation Therefore, the high level of HIF-1a has potency in The amino acid sequence of HIF-1a associating with further massive tumors 3. (Entry PDB code : 4z1v) was retrieved from RSCB Natural bioactive compounds which Protein Database 19. The attached ligand in the are derived from plants have been used for protein structure was removed from the binding maintaining health and remedies in many years. site and saved to a new file format: pdbqt. The The phytochemical constituents in plants have been Gasteiger charges and the solvation condition a critical pipeline for the discovery of bioactive were added to the protein structure using the substances in pharmaceutical field 10. A. muricata AutoDockTool 20. or Graviola has been greatly pressumed to have Ligand structures preparation valuable natural products that play important role Ligands which are AGEs consisting of in affecting anticancer activity 11. A. muricata eight 3D structures of natural bioactive compounds leaves have been used to investigate of numerous originally belong to A. muricata and one anticancer numbers of human diseases, including cancers 10. drug for molecular docking experiments and their The highly constituents screening are most possibly conformational energy were minimized by using affected by its major bioactive components MMFF94 force field. 8 molecules of AGEs. The known as annonaceous Acetogenins (AGEs) 12. molecule structures are retrieved from PubChem Many studies reported that isolated AGEs from database (Fig. 3). The structures were scored different extracts of the plants have significant based on their physicochemical properties under antiproliferative effects against various cancer Chemicalize (ChemAxon) and Molsoft platforms 21, cell lines 10. However, some of these studies have 22. These physicochemical properties are important defined the staple mechanism of action. Recent for developing drug candidate in every stages from in-vitro studies showed inhibition action of ethyl design to pre-clinical study. acetate extract from A. muricata leaves combating Drug likeness analysis of A. muricata bioactive lung cancer cells (A549) and colon cancer cells compounds (HCT-116 and HT-29) 10, 13. The leaf extract was 3D structures of Cinchona alkaloids capable of inducing colon carcinoma and lung were analyzed using a program based on the cancer cells apoptosis by way of mitocondrial physicochemical properties, Molsoft Drug - route. This antiproliferative effect was associated Likeness. Determination of physicochemical with cell cycle involved in the G1 phase. Moreover, properties is important in the development of drug the migration and invasion of colon cancer cells candidates in all stages ranging from study design were significantly halted by the leaf extract 14-16. through pre- clinical trials 22. The aim of this research is to determine Molecular docking of HIF-1a and A. muricata the inhibition mechanism by bioactive compounds bioactive compounds of A. muricata interact with HIF-1a. To study The three dimension structure of protein the binding interactions of bioactive compounds and ligands were prepared in pdb format. Molecular HANDAYANI et al., Biomed. & Pharmacol. J, Vol. 11(3), 1301-1307 (2018) 1303 docking simulation was run by Autodock Vina RESULTS (Vina, The Scripps Institute) 23. The AutodockTools (ADL) was utilized in minimizing energy and Bioactive compounds isolated from A. adding the partial charges of polar hydrogens of muricata , including annomuricin A, annomuricin receptor (protein). The ligands were prepared with B, annomuricin C, annomuricin E, annomutacin, flexible torsion angles and the protein was prepared murihexocin A, murihexocin B, murihexocin C, in a static (rigid) form. Furthermore, protein and and gefitinib (Figure 1) were docked into binding ligands were kept in pdbqt formats which suitable pocket of HIF-1a. for docking simulation. The affinity binding were The lowest binding energy (LBE) to the calculated as total intermolecular energies (kcal/ target protein was murihexocin A (-7.9 kcal/mol). mol) which involved hydrogen bond, Van Der The binding interactions between A. muricata Walls force, desolvation and electrostatic energies. bioactive compounds with of HIF-1a binding On the other hand, the apropriate torsion angles pocket residues were analysed as shown in of ligand is also induced as internal ligand energy. Table 1. The docking program evaluated the lowest binding From docking result, murihexocin A energy (LBE) to obtain the best binding mode. The compound interacted hydrophobically with Thr- Root-Mean-Square Deviation (RMSD) which less 183, Trp-296, Tyr-102, His-199, Gln-147, Gln-203, than 2.0 Å was scored during running docking Glu-202, Ser-184, Pro-235, Gln-239, Tyr-103, and program. Ile-281 in the binding pocket of HIF-1a. Fig. 1. AGEs compound of A. muricata leaves: a. annomuricin A; b. annomuricin B; c. annomuricin C; d. annomuricin E; e. annomutacin; f. murihexocin A; g. murihexocin B; h. murihexocin C; i. gefitinib 1304 HANDAYANI et al., Biomed. & Pharmacol. J, Vol. 11(3), 1301-1307 (2018) The molecular interaction