Bangladesh Journal of Pharmacology Research Article Computational drug designing of fungal pigments as potential aroma- tase inhibitors BJP Introduction Monascus pigments biosynthzied from fungi of genus Monascus, typically comprised six major valuable bio- active polyketide azaphilone pigments; ankaflavin (yellow); rubropunctatin (orange), monascin, rubro- punctamine (purple red), monascorubrin, and mona- scorubramine, (Lee et al., 2006; Mapari et al., 2010; Lee et al., 2012; Figure 1). Breast cancer (BC) is one of the most widespread and multifarious disease that accounts in women approxi- mately 22.9% of all cancers in world. The means and molecular biology of breast cancer has not been com- pletely investigated (Standish et al., 2008) and limita- tions of existing therapies prompt investigators to carry on penetrating for more effective molecules (Sotiriou et al., 2006; Sotiriou and Pusztai, 2009). Estrogen hormones are caught up in the growth of breast cancer which is most frequent in post-meno- pausal women (Maiti et al., 2007). Aromatase is a cytochrome P450 enzyme playing a significant role in the translation of androgen to estrogen (Jongen et al., 2005). Aromatase inhibitors can stop the creation of estrogen, as a result can orderly reduces the develop- ment of breast cancer cells. Therefore, aromatase inhi- bitors have become attractive therapy (Chomcheon et al., 2009; Sureram et al., 2012). Although various compounds have been studied for their aromatase inhibition potential but complexity of enzyme structure lemmatizes the use of these com- pounds. An alternative strategy to design new inhibi- tors relies on ligand-based virtual screening by using different compounds. Therefore in this study, we have used natural poly- ketide pigments of Monascus fungus for virtual screen- ing. Our study would be supportive in adding new lead molecules and drug targets to smooth the progress of the diagnosis and the management of breast cancers. A Journal of the Bangladesh Pharmacological Society (BDPS) Bangladesh J Pharmacol 2014; 9: 575-79 Journal homepage: www.banglajol.info Abstracted/indexed in Academic Search Complete, Agroforestry Abstracts, Asia Journals Online, Bangladesh Journals Online, Biological Abstracts, BIOSIS Previews, CAB Abstracts, Current Abstracts, Directory of Open Access Journals, EMBASE/Excerpta Medica, Global Health, Google Scholar, HINARI (WHO), International Pharmaceutical Abstracts, Open J-gate, Science Citation Index Expanded, SCOPUS and Social Sciences Citation Index ISSN: 1991-0088 Abstract The existing aromatase inhibitors produced unwelcome effects impose the discovery of novel drugs with privileged selectivity, a reduced amount of toxi -city and humanizing potency. In this study, we illuminate the binding mode of polyketide azaphilanoid pigments monascin, ankaflavin, monascorubrin and monascorubramine isolated from Monascus fungus to the aromatase by molecular docking. The 3-dimensional structure of aromatase enzyme (PDB: 4KQ8) was obtained from the Protein Data Bank. PatchDock docking software was used to analyze structural complexes of the aromatase with monascus pigments. Comparatively, the AutoGrid model presented the most briskly constructive binding mode of monascin to aromatase. Docked energies in kcal/mol are: monascin:-13.2; monascorubramine:-12.8, monascorubrin:-12.3; ankaflavin:-10.5. These outcomes exposed these ligands could be potential drugs to treat hormone dependent breast cancer. Article Info Received: 13 September 2014 Accepted: 24 October 2014 Available Online: 7 November 2014 DOI: 10.3329/bjp.v9i4.20435 Cite this article: Fatima N, Kalsoom A, Mumtaz A, Muhammad SA. Computational drug designing of fungal pigments as po- tent ia l aromatase inhibi tors. Bangladesh J Pharmacol. 2014; 9: 575- 79. This work is licensed under a Creative Commons Attribution 3.0 License. You are free to copy, distribute and perform the work. You must attribute the work in the manner specified by the author or licensor. Computational drug designing of fungal pigments as potential aromatase inhibitors Nighat Fatima1, Abida Kalsoom2, Amara Mumtaz2 and Syed Aun Muhammad3 1Department of Pharmacy, COMSATS Institute of Information Technology Abbottabad 22 060, Pakistan; 2Department of Chemistry, COMSATS Institute of Information Technology Abbottabad 22 060, Pakistan; 3Institute of Molecular Biology and Biotechnology, Bahauddin Zakariya University Multan, Pakistan. Materials and Methods Accession of target protein: The 3-dimensional structure of aromatase (cytochrome p-450) enzyme (PDB: 4KQ8) was downloaded from the RCSB protein Data Bank (http://www.rcsb.org/pdb/home/home.do). Ligand selection: The chemical structures of monascin, monascorubrin, monascurubramine and ankaflavin were obtained from PubChem database. structures were drawn by using ChemBioDraw and MOL2 format of these ligands were changed to PDB file using Open Babel tool earlier to upload onto PachDock virtual software. Target and ligand optimization: For docking analysis, PDB coordinates of the target protein and ligands molecules were optimized by Drug Discovery Studio version 3.0 software and UCSF Chimera tool respectively. These coordinates had minimum energy and stable confor- mation. Analysis of target active binding sites: The active sites are the coordinates of the ligand in the original target protein grids and these active binding sites of target protein were investigated using the MetaPocket 2.0 virtual tools (Huang, 2009). Molecular docking analysis: A computational approach, ligand-target docking was adopted to investigate the aromatase structural complexes (cytochrome p-450) as target with chosen compounds (ligands) in order to comprehend the structural basis of this particular protein target. To begin with, Platinum software web server was used for the investigation of protein–ligand attraction of complexes and to determine their hydro- philic or hydrophobic properties (Pyrkov et al., 2009). As a final point, these complexes were subjected to docking by using PatchDock virtual docking software (Schneidman-Duhovny et al., 2005). The energy of inter- action of ligands with the target enzyme is depicted as “grid point”. At every step of the simulation, the ligand -protein interaction of energy was monitored by atomic affinity potentials computed on a grid. The remaining parameters were set as default. Results and Discussion Molecular docking is an important source which depicts helpful inquires regarding drug-receptor inter- actions and is repeatedly used to forecast the binding orientation of small molecule behaving as potential candidates for drug with their target protein to guess their potential activity and affinity towards their target (Vijesh et al., 2013). Therefore molecular docking studies of biologically active fungal pigments provided superior perceptive of the drug-receptor interaction. The minimum binding energy indicated that the aroma- tase protein (target enzyme) was successfully docked with ligands molecules (Table I). The possible binding modes of ligands at aromatase active sites have been shown in Figure 2. Aromatase protein residues Arg, Met, Phe, Ser, Tyr, Pro, Asn, Glu, Val, Lys, Gln, Ala, Ile was formed H-bond with ligands molecules. Monascin showed relatively good binding affinity (-13.2 kcal/ mol) as compared to other ligands. The division of hydrophilic or hydrophobic properties of aromatase and its binding site anticipated against the surface of the ligands, and their complementarity was inquired. Less hydrophilicity of aromatic groups is in common practice . These stacking interactions are based to grade the molecular docking to investigate the drug- target interaction. Taking into account the hydrophobic properties, molecular hydrophobicity potential was 576 Bangladesh J Pharmacol 2014; 9: 575-579 Figure 1: Different monascus pigements used as aromatase inhibitors (A) Monascin, (B) Ankaflavin, (C) Monascorubrin, (D) Monascorubramine calculated for these interacting molecules (aromatase protein target and ligands) (Muhammad et al., 2014) which were the bases for the study of progress of molecular dynamics run base on hydrophobicity clustering on the surface of membrane also shows the membrane-molecule interaction. The docking of ligands-aromatase target came up that all the computationally forecast lowest energy complex- es of aromatase are stabilized by stacking interactions and intermolecular hydrogen bonding. It was also originated that A, SA, OA, HD, N are the ligand atoms are responsible for in docking with the enzyme. The AutoGrid model came up with the consequences that the most energetically acceptable binding mode of monascin to aromatase. The monascin as ligand was docked into the generated combined grids, the RMSD from native pose and the binding energy were investi- gated. It is found that the weight averaged grids per- forms excellently. Base on RMSD values it was found that the ligand showed significant interaction with tar- get proteins as compared to other compounds. Beside RMSD clustering, PatchDock has deliberate the binding free energies of these interactive molecules to discover the finest binding mode. The intended final docked energies for monascin was -13.2 kcal/mol, for anka- flavin -10.5 kcal/mol (Figure 3), for monascorubrin-12.3 kcal/mol and -12.8 kcal/mol monascorubramine were observed (Figure 4). Docking results came up with the statement that these ligands molecules can interact with aromatase protein target more significantly and can help in the study of cancer research. Bangladesh J Pharmacol 2014; 9: 567-574 577 Figure 2: Potential binding sites of aromatase (cytochrome p- 450) with binding site center X: 78.117, Y: 45.497, Z: 58.722 indicating amino acids residues: Arg, Met, Phe, Ser, Tyr, Pro, Asn, Glu, Val, Lys, Gln, Ala, Ile (A) Binding sites for ligands (B) maximized region A B Table I Energy and RMSD values obtained during docking analysis of ligands molecules and aromatase enzyme as target protein SL. No. Complex Binding energy RMSD/UBa RMSD/LBa 1 Aromatase_Monascin -13.2 0 0 2 Aromatase_Monascin -12.8 1.1 1.5 3 Aromatase_Monascin -12.3 2.2 1.8 4 Aromatase_Monascin -11.2 2.9 1.9 1 Aromatase_Ankaflavin -10.5 0 0 2 Aromatase_Ankaflavin -9.4 2.3 1.3 3 Aromatase_Ankaflavin -8.7 2.9 1.8 4 Aromatase_Ankaflavin -8.1 3.7 2.0 1 Aromatase_Monascorubrin -12.3 0 0 2 Aromatase_Monascorubrin -11.4 2.2 1.6 3 Aromatase_Monascorubrin -9.3 3.1 1.9 4 Aromatase_Monascorubrin -7.9 3.9 2.1 1 Aromatase_Monascorubramine -12.8 0 0 2 Aromatase_Monascorubramine -12.1 2.5 1.9 3 Aromatase_Monascorubramine -10.5 3.4 2.1 4 Aromatase_Monascorubramine -9.9 4.1 3.0 aRMSD/UB: Root mean square deviation/upper bond; RMSD/LB: root mean square deviation/lower bond Conclusion Docking studies of the monascus, ankaflavin, monas- corubrin and monascorubramine with aromatase enzyme showed that these natural compounds are good ligands which dock well with aromatase target. References Akihisa T, Tokuda H, Ukiya M, Kiyota A, Yasukawa K, Sakamoto N, Kimura Y, Suzuki T, Takayasu J, Nishino H. Anti-tumor-initiating effects of monascin, an azaphilonoid pigment from the extract of Monascus pilosus fermented (red- mold rice). Chem Biodivers. 2005; 2: 1305-07. Chomcheon P, Wiyakrutta S, Sriubolmas N, Ngamrojana- vanich N, Kengtong S, Mahidol C, Ruchirawat S, Kittakoop P. Aromatase inhibitory, radical scavenging, and antioxi- dant activities of depsidones and diaryl ethers from the endophytic fungus Corynespora cassiicola L36. Phytoche- mistry. 2009; 70: 407-13. Huang B. MetaPocket: A meta approach to improve protein ligand binding site prediction, OMICS. 2009; 13: 325-30. Jongen VHWM, Thijssen JHH, Hollema H, Donker GH, Santema JG, Van der zee AGJ, Heineman MJ. Is aromatase cytochrome P450 involved in the pathogenesis of endometrioid endometrial cancer? Int J Gynecol Cancer. 2005; 15: 529-36. Lee BH, Hsu WH, Chang YY, Kuo HF, Hsu YW, Pan TM. Ankaflavin: A natural novel PPARγ agonist up-regulates Nrf2 to attenuate methylglyoxal-induced diabetes in vivo. Free Radic Biol Med. 2012; 53: 2008-16. Lee CL, Wang JJ, Kuo SL, Pan TM. Monascus fermentation of dioscorea for increasing the production of cholesterol- lowering agent-moacolin K and antiinlammation agent- monascin. Appl. Microbiol Biotecehnol. 2006; 72: 1254-62. Maiti A, Cuendet M, Croy VL, Endringer DC, Pezzuto JM, Cushman M. Synthesis and biological evaluation of (+/-)- abyssinone II and its analogues as aromatase inhibitors for chemoprevention of breast cancer. J Med Chem. 2007; 50: 2799-806. Mapari SA, Thrane U, Meyer AS. Fungal polyketide azaphilone pigments as future natural food colorants? Trends in Biotechnol. 2010; 28: 300-07. Muhammad SA, Ali A, Ismail T, Zafar R, Ilyas U, Ahmad J. In 578 Bangladesh J Pharmacol 2014; 9: 575-579 Figure 3: Visuals of docking procedure obtained using PatchDock virtual tool (A) The AutoGrid dimensions between monascin and aromatase (cytochrome p-450) are: grid center X: 19.1414, Y: 21.2513, Z: 11.3412 with dimension (Angstrom) X: Y: Z: 25.000 (B) maximized region presenting confirmation and pose of monascin ligand (C) The AutoGrid dimensions between ankaflavin and aromatase (cytochrome p-450) are: grid center X: 20.0211, Y: 23.2021, Z: 12.1598 with dimension (Angstrom) X: Y: Z: 25.000 (D) maximized region presenting confirmation and pose of ankaflavin ligand A B C D silico study of anti-carcinogenic lysyl oxidase-like 2 inhibitors. Com Biol Chem. 2014; 51: 71-82. Pyrkov TV, Chugunov AO, Krylov NA, Nolde DE, Efremov RG. PLAT-INUM: A web tool for analysis of hydrophobic/ hydrophilic organization of biomolecular complexes. Bioinformatics. 2009; 25: 1201-02. Schneidman-Duhovny D, Inbar Y, Nussinov R, Wolfson HJ. PatchDock and SymmDock: Servers for rigid and symmetric docking. Nucl Acids Res. 2005; 33: W363-67. Sotiriou C, Pusztai L. Gene-expression signatures in breast cancer. N Engl J Med. 2009; 360: 790-800. Sotiriou C, Wirapati P, Loi S, Harris A, Fox S, Smeds J, Nordgren H, Farmer P, Praz V, Haibe-Kains B, Desmedt C, Larsimont D, Cardoso F, Peterse H, Nuyten D, Buyse M, Van de Vijver MJ, Bergh J, Piccart M, Delorenzi M. Gene expression profiling in breast cancer: Understanding the molecular basis of histologic grade to improve prognosis. J Natl Cancer Inst. 2006; 98: 262-72. Standish LJ, Sweet ES, Novack J, Wenner CA, Bridge C, Nelson A, Martzen M, Torkelson C. Breast cancer and the immune system. J Soc Integr Oncol. 2008; 6: 158-68. Sureram S, Wiyakrutta S, Ngamrojanavanich N, Mahidol C, Ruchirawat S, Kittakoop P. Depsidones, aromatase inhibi- tors and radical scavenging agents from the marine- derived fungus Aspergillus unguis CRI282-03. Planta Med. 2012; 78: 582-88. Author Info Nighat Fatima (Principal contact) e-mail: nighatmrl@yahoo.com Bangladesh J Pharmacol 2014; 9: 575-579 579 Figure 4: Visuals of docking procedure obtained using PatchDock virtual tool (A) The AutoGrid dimensions between Monascoru- brin and aromatase (cytochrome p-450) are: grid center X: 18.1001, Y: 22.2412, Z: 10.1201 with dimension (Angstrom) X: Y: Z: 25.000 (B) maximized region presenting confirmation and pose of Monascorubrin ligand (C) The AutoGrid dimensions between Monascorubramine and aromatase (cytochrome p-450) are: grid center X: 19.0011, Y: 21.1231, Z: 11.0251 with dimension (Angstrom) X: Y: Z: 25.000 (D) maximized region presenting confirmation and pose of Monascorubramine ligand A B C D A B C D: Conclusion: DatePrinted: This article was downloaded by you on: Oct 05, 2018