Computational drug designing of fungal pigments as potential aromatase inhibitors

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 potent ia l aromatase inhibi tors. Bangladesh J Pharmacol. 2014; 9: 57579. 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 Department of Pharmacy, COMSATS Institute of Information Technology Abbottabad 22 060, Pakistan; Department of Chemistry, COMSATS Institute of Information Technology Abbottabad 22 060, Pakistan; Institute of Molecular Biology and Biotechnology, Bahauddin Zakariya University Multan, Pakistan.

Breast cancer (BC) is one of the most widespread and multifarious disease that accounts in women approximately 22.9% of all cancers in world. The means and molecular biology of breast cancer has not been completely investigated (Standish et al., 2008) and limitations 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 development of breast cancer cells. Therefore, aromatase inhibitors 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 compounds. An alternative strategy to design new inhibitors relies on ligand-based virtual screening by using different compounds.
Therefore in this study, we have used natural polyketide pigments of Monascus fungus for virtual screening. 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.
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 conformation.
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 hydrophilic 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 interaction 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 interactions 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 aromatase 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 drugtarget interaction. Taking into account the hydrophobic properties, molecular hydrophobicity potential was 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 complexes 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 investigated. It is found that the weight averaged grids performs excellently. Base on RMSD values it was found that the ligand showed significant interaction with target 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 ankaflavin -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.

Conclusion
Docking studies of the monascus, ankaflavin, monascorubrin and monascorubramine with aromatase enzyme showed that these natural compounds are good ligands which dock well with aromatase target.

Financial Support
Self-funded