In Silico Prediction of Heliannuol A, B, C, D, and E Compounds on Estrogen Receptor Β Agonists
Total Page:16
File Type:pdf, Size:1020Kb
Indonesian Journal of Cancer Chemoprevention, February 2021 ISSN: 2088–0197 e-ISSN: 2355-8989 In Silico Prediction of Heliannuol A, B, C, D, and E Compounds on Estrogen Receptor β Agonists Roihatul Mutiah, Alif Firman Firdausy, Yen Yen Ari Indrawijaya, Hibbatullah* Department of Pharmacy, Faculty of Medical and Health Sciences, Maulana Malik Ibrahim State Islamic University of Malang, Indonesia Abstract Heliannuols has a benzoxepine ring that produces anticancer activity by the inhibition mechanism of phosphoinositide 3 kinases (PI3K). Heliannuols are a compound that can be found in the leaves of sunflower Helianthus( annuus L.). The purpose of this study is to predict interactions, toxicity, physicochemical, and pharmacokinetics of Heliannuol A, B, C, D, and E based in silico as candidate anticancer drugs. Estrogen receptor beta (ERβ) is a new potential therapy for glioma with an antiproliferative effect. Ligands agonist ERβ have the potential activity to inhibit the proliferation of glioma cells and the discovery of this ligand has opened new therapy through the ERβ to prolong survival in cancer patients. Prediction of physicochemical properties based on Lipinski rules and penetrate in the blood-brain barrier. Receptor validation shows that 2I0G(A) has a smaller RMSD value than 2I0G(B), receptor validation is valid if the RMSD value less than 2. The result of molecular docking shows that Heliannuols comply with Lipinski rules and have low toxicity. Heliannuols also have a similar amino acid with comparison drug (Erteberel), but the rerank score of Erteberel still lower than Heliannuols. Keywords: Helianthus annuus, Heliannuols, estrogen receptor β (ERβ), in silico, toxicity. INTRODUCTION using HPLC with hexane and ethyl acetate solvents (Macias, et al., 1994; Macias, et al., 2000; Macias, Sunflower plant (Helianthus annuus L.) has et al., 2002). The new anticancer treatment therapy many chemical compounds, especially in the leaves involves estrogen receptor beta-agonists (ERβ) of this plant containing monoterpene, diterpene, by the pathway of inhibiting cancer growth and alkaloids, phenols, and lactone sesquiterpenes anticancer therapy with estrogen receptors can (Ceccarini, 2004). Heliannuols are one of the reduce glioma cell proliferation (Sareddy, et al., sesquiterpenes which has a benzoxepine ring and 2012; Sareddy and Vadlamudi, 2015). anticancer activity produced by the inhibition ERβ has quite a different function mechanism of phosphoinositide 3 kinases (PI3K) from ERα, ERβ has antiproliferative action and (Ghantous, et al., 2010; Hefron, et al., 2011; Kuntala, et al., 2017; Ren, et al., 2003). Heliannuol A, B, C, D, and E were obtained from sunflower Submitted: September 29, 2020 Revised: February 9, 2021 leaves that had been immersed in water for 24 h and Accepted: February 15, 2021 extracted with chloroform (CH2Cl2) then observed *Corresponding author: [email protected] 37 Mutiah, et al., 2021 Indones. J. Cancer Chemoprevent., 12(1), 37-45 overexpression of ERβ can reduce proliferation in and Molegro Virtual Docker 6 (Lisensi expires on: cancer cells (Strom, et al., 2004; Hartman, et al., January 01, 2099). 2009; Nilsson, et al., 2011). This research focuses on estrogen receptor beta agonists because it has the Target and Template Selection potential as an antitumor therapy agent by activating The target of this research is the estrogen beta subtypes and preventing proliferative effects receptor beta 2I0G from Protein Data Bank (https:// on cancer. Ligands that are included as estrogen www.rcsb.org). The test compounds of this research receptor beta-agonists that have the potential activity are Heliannuol A, B, C, D, and E where the SMILES to inhibit the proliferation of glioma cells are DPN code from Chem Bio Draw 12.0 and the Erteberel (Diarylpropionitrile), MFF101, Liquiri-tigenin, and SMILES code from PubChem. LY500307 (comparison drug (Erteberel)) (Sareddy, et al., 2012; Sareddy and Vadlamudi, 2015). Prediction of Physicochemical, Erteberel is one of the agonists ERβ that has high Pharmacokinetic, and Toxicity Properties potential in reducing glioma cell proliferation and Prediction of physicochemical properties it can penetrate the blood-brain barrier (Sareddy, et using the pkCSM online tool and the Protox II al., 2016). The discovery of a specific agonists ERβ online tool by entering the SMILES code of the provided naturally has opened the development of compound. The prediction of physicochemical new therapies through the estrogen receptor beta to properties is based on Lipinski’s rules and the prolong survival in cancer patients (Sareddy and ability of compounds to penetrate the brain Vadlamudi, 2015). barrier membrane. Lipinski’s rules use to evaluate Based on the research by Sareddy, et al. compounds that have pharmacological activity with (2012) show that ERβ pathway activation is a predictable physical and chemical properties as potential therapeutic target for glioma because drug candidates for humans. While the prediction ERβ agonists are under clinical trials and are well of pharmacokinetic properties using the online tolerated with fewer side effects. The use of ERβ pkCSM tool is based on the absorption, distribution, as a therapeutic agent can be extended to clinical metabolism, and excretion parameters. Toxicity use and is predicted to be a new class of drugs prediction using the pkCSM online tool and the for treating glioma. Future studies examining the Protox II online tool based on several parameters mechanism of ERβ in cancer progression will be and toxicity classes according to the Globally useful for maximizing treatment using ERβ natural Harmonized System (GHS). ligands (Sareddy and Vadlamudi, 2015). The purpose of this study was to determine the activity Molecular Docking of Heliannuols on estrogen receptor beta and the Molecular docking using Molegro Virtual lack of information on estrogen receptor beta- Docker 6.0 to know the interaction between the agonists as therapeutic agents in cancer patients, test compound and the receptor. 2D structures of encouraging researchers to compile this study. Heliannuol A, B, C, D, and E were drawn using Chem Draw Ultra 12.0 and to know the SMILES MATERIALS AND METHODS code of the test compounds then 3D structures were made using Avogadro and energy minimization Software was performed. Validation of molecular docking This research uses Chem Bio Draw Ultra ERβ (2I0G) using Molegro Virtual Docker 6.0 and 12.0 application, Avogadro, pkCSM online tool, a docking simulation process was also performed. Protox II online tool (accessed on January 2020), Validation of molecular docking result is the Root 38 Indonesian Journal of Cancer Chemoprevention, February 2021 ISSN: 2088–0197 e-ISSN: 2355-8989 Mean Square Deviation (RMSD), the docking of compounds to penetrate the brain barrier mem- simulation process can be run if the RMSD value brane, and toxicity are shown in Table 1. was less than 2Å. The parameters used in docking The results of physicochemical prediction simulation include rerank score, bond distance, and show that all compounds comply with Lipinski type of interaction. rules, it can be predicted that all compound easy to absorption and has good permeability. The result RESULTS of TPSA showed that all compounds have a TPSA value of less than 80, it can be predicted that all Prediction of Physicochemical Properties compounds could penetrate the brain barrier mem- and Toxicity brane (Hughes, 2008). Heliannuol A, B, and C are Prediction of physicochemical properties in toxicity class 4 with Lethal Dose (LD50) 300-2000 using Lipinski rules of five with several parameters mg/kg. Heliannuol D and Erteberel are in toxicity includes molecular weight, the partition coefficient class 5 with LD50 2000-5000 mg/kg. Heliannuol A, of octanol/water (log P), Hydrogen Bond Donors B, C, D, E, and Erteberel were predicted to be non- (HBD), and Hydrogen Bond Acceptors (HBA). toxic in the AMES Mutagenic test, not toxic in the Topological Polar Surface Area (TPSA) also a Hepatotoxicity test, and did not cause skin irritation parameter to determine compound can penetrate in the Skin Sensitization test. the brain barrier membrane (Kelder, et al., 1999). The physicochemical test was used to predict that Prediction of Pharmacokinetic Properties the compound has good permeability, high absorp- Prediction of pharmacokinetic properties tion, and can penetrate the blood-brain barrier. The is based on the prediction of absorption, distribu- results of physicochemical properties, the ability tion, metabolism, and excretion (ADME) using the Table 1. Prediction of physicochemical and toxicity. Description: *pkCSM online tool; **Protox II online tool. pkCSM online tool. The results of the prediction of website (http://biosig.uni-melb.edu.au/pkcsm/the- pharmacokinetic properties are shown in Table 2. ory), the compound has a high skin permeability if The result from the prediction of phar- the Log value of Kp>-2.5 cm/h, and all compound macokinetic properties shows ADME of the com- are predicted have high skin permeability. Caco2 pound, in this research shows that all compound permeability is good if the Papp value>0.90 cm/s is predicted to have good intestinal absorption and all compounds are predicted to have good per- with the value is more than 80% (Chander, et al., meability. Heliannuol C and E do not include as 2017). According to the pkCSM online tool page Pgp substrates, then they are not predicted to be 39 Mutiah, et al., 2021 Indones. J. Cancer Chemoprevent., 12(1), 37-45 Table 2. Prediction of pharmacokinetic properties. Description: Vdss (Volume of Distribution at Steady State); BBB (Blood Brain Barrier); CNS (Central Nervous Sys- tem); CYP2D6 (Cytochrome P2D); CYP3A4 (Cytochrome 3A4); CYP1A2 (Cytochrome 1A2); CYP2C19 (Cyto- chrome 2C19); CYP2C9 (Cytochrome 2C9); Renal OCT2 (Renal Organic Cation Transporter 2). *The highest value of comparative drugs; **The highest value of the first test compound; +The second highest value of the test compound; ++Different result from other. removed from the target cell (Chakraborty and Ra- high if the Log Vdss value >0.45, all compound makrishnan, 2016). Heliannuol A, B, C, D, E, and has good distribution volume.