Semiemperical Investigation of the Postmenopausal Breast Cancer Treatment Potential of Xanthone Derivatives

Sajad Shahbazi1, Ananya Kuanar2, Deepak Reddy Gade3, Dattatreya Kar2, Anish Shrivastava2, Pavan Kunala4 and Manoj Kumar Mahto5* 1Department of Biotechnology, Panjab University, Chandigarh, India 2Center of Biotechnology, Siksha ‘O’ Anusandan University, Orissa, India 3Department of Pharmaceutical Chemistry, JNTU-OTRI, Anantapur, Andhra Pradesh, India 4Department of Biology, Arkansas State University, Jonesboro, Arkansas, USA 5Department of Bioinformatics, GVK Biosciences Private Limited, Andhra Pradesh, India


Introduction
From age's breast cancer, a type of ductal carcinoma is one of the most common cancers affecting female originating from the inner lining of the milk ducts of breast tissue [1]. Breast cancer alone comprises about 22.9% of all cancers in women. About 4.6 hundred thousand deaths were recorded worldwide in 2008 only due to breast cancer and it is about 13.7% of cancer deaths in women [2]. Many pathophysiological reasons like mutations in genes, inherited genetic predisposition, environmental, dietary and hormonal exposure leads to the development of breast cancer [3]. Many potential targets like Epidermal Growth Factor Receptor (EGFR), Vascular Endothelial Growth Factor Receptor (VEGFR), and enzymes like Aromatase, which catalyzes aromatization of androgens to estrogen, and also hormone receptors like estrogen receptor were identified for treatment of breast cancer [4]. The role of Aromatase in estrogen synthesis is depicted in Figure 1.
Aromatase Inhibitors (AIs), as a class of drugs in the treatment of postmenopausal women's breast and ovarian cancer, are generally categorized into two groups: irreversible steroidal inhibitors and nonsteroidal inhibitors. The former, by binding to the active site, deals with permanently deactivating it, but the latter, by reversible competition for Aromatase enzyme with androgen, hinders binding androgens to the Aromatase enzyme [5]. The reversible Androgen-Aromatase binding inhibition contributes to the inhibition of androgens to estrogens conversions, in turn, preventing the steroidogenesis. It considerably reduces the circulation and the intra-tissue concentrations of estrogens. The reversible is more efficient than the permanent inhibition of Aromatase in chemotherapy of breast cancer, due to maintain the estrogen levels in normal rates approximately 70, 20 and 400 pmol/L for estrone, estradiol, and estrone sulfate, respectively in plasma for naturally postmenopausal women [6]. Figure 2 shows the mode of action of Aromatase inhibitor in Estron synthesis pathway.
AIs can also be used to regulate Aromatase enzyme through other pathways and receptors such as, the modulation of Liver Receptor Homolog-1 (LRH-1) and orphan receptor. AIs regulate Aromatase in adipose tissue, testis, and granulose cells as well as contribute to over-expression of Aromatase in breast cancer patients [7]. Selective Aromatase modulators (SAMs) may be found based on the evidence for tissue-specific promoters of Aromatase expression [8].
With the clinical success of several synthetic Aromatase inhibitors (AIs) for the treatment of postmenopausal breast cancer, researchers have been investigating the potential of natural products as AIs [9]. The medicinal use of Natural Products (NPs) has a long history in treatment of various diseases. NPs have been utilized in the forms of herbal remedies, purified compounds, and precursor for combinatorial chemistry. The NPs can be considered for the inhibition of Aromatase due to the minimum side effects which can be as a result of natural product matrix components [9]. These components can be considered as the chemotherapeutic agents for the future clinical trials in breast cancer chemoprevention, treating postmenopausal breast cancer or preventing secondary recurrence of breast cancer [10]. It is also possible to synthetically modify of natural product scaffolds to enhance the inhibitory potential and improve the ADMET properties of such compounds [11].

Selection and preparation of protein
Aromatase structures with PDB ID-3EQM of obtained from the RCSB Protein Data Bank (http://www.rcsb.org/pdb/) with X-ray diffraction resolutions of 2.90 Å. Preparation of the retrieved protein was performed by using protein Preparation Wizard of Schrodinger explicit hydrogens and topological duplicates. The 2D structures of all compounds were presented in Table 1.

Receptor-Ligand interactions
Docking studies were performed by using Glide 5.5 module in Extra Precision (XP) mode [20][21][22] and the molecular mechanics/ Generalized Born Surface Area (MMGBSA) [23] for interaction of each ligand-protein complex has been calculated via Prime 2.1 application of Schrodinger Suite 2009. The results of docking and MMGBS are available in Table 1

Preparation of ligands
Preparation of ligands was done by using LigPrep 2.3 module of Schrodinger Suite 2009 using the OPLS force-field 2005 at biologically relevant PH. It performed by assigning the protonation states, including the disconnection of group I metals in simple salts, the deprotonation of strong acids and the protonation of strong bases, while adding   been detected by the pose of ligand presented in receptor crystal which was extracted from PDB structure.

Preparation of the receptor and ligands
Aromatase PDB structure has been downloaded from RCSB Protein Data Bank (http://www.rcsb.org/pdb/). Then via ProPrep wizard of Shrodinger 2009, the protein preparation was carried out by removing the disulphide and trisulphide bonds of a protein, all water molecules except those which bind through hydrogen bonds with the residues in the binding side of a protein and adding hydrogen's to the molecule to satisfy the valences of the molecule. It employed OPLS 2005 force-field with RMSD as 0.30. The binding site of proteins has * DG (ΔG bind) = G complex -(G protein + G ligand) where ΔG bind is Ligand binding energy in kcal/mol; XPG: Extra Precision Glide score, LPE = Ligand Potential Energy in kcal/mol units, HB-R = Hydrogen Bound Residue numbers, HB-D = hydrogen bound distance in Å.   body) [24]. The stereoisomer has generated at most 32 combinations per ligand. The potential energies of ligands are presented in Table 1.
According to the affinity values, binding energy and number of hydrogen bonds with the active site of receptor presented in Table 1

ADMET investigation
For designing a druggable molecule, it is necessary to check the drug ability of ligands in terms of oral and intestinal absorption level, the ability of ligands to distribute through the blood stream, level of metabolism, the ability of excretion from the body besides their toxicity profiling. For estimation of drug pharmacokinetics and metabolism, QikProp software investigates around 24 molecular descriptors, which determine the #star parameter. The recommended values for #star by Schrodinger is a range 0-5 for the computed property of a molecule out of the range for 95% of known drugs, including MW, dipole, IP, EA, SASA, FOSA, FISA, PISA, WPSA, PSA, volume, #rotor, donorHB, accptHB, glob, QPpolrz, QPlogPC16, QPlogPoct, QPlogPw, QPlogPo/w, logos, QPLogKhsa, QPlogBB, #metabol.
The likely oral availabilities of the compounds have been evaluated using [16] MW ≤ 500 Da, accptHB ≤ 10, donorHBD ≤ 5, QPlog Po/w ≤ 5, #rotor ≤ 10 [25] which are known as the Lipinski's criteria. The oral absorption or the likelihood of the oral availability is computed by using parameters including log S>-5.7, QPPCaco>22 nm/s and # primary metabolites<7 which are known as the Jorgensen's famous ''Rule of Three'' (ro3). The Jotgensen's ro3 beside human oral absorption percentage, the predicted qualitative human oral absorption, the predicted of human gut barrier absorption, QPPCaco, and CIlog S are used to predict the bioavailability [26]. The permeability prediction depends on the molecular properties such as the size, flexibility which depends on the #rotor [27], overall lipophilicity, shape, and the capacity to make hydrogen bonds.
The QPlogBB is used to predict the blood brain barrier permeability for each compound and accessibility of bio actives for central nervous system. It has contrary relation to the polarity of the compound. For CNS availability and activity of compounds, in addition of OPlogBB, the CNS activity, and QPPMDCK also have to be considered. The recommended range for CNS activity, QPlogBB and QPPMDCK are (-2-+2), (-3-1.2) and (25-500 nm/s) respectively.
The QplogKp is known as the skin permeability parameter. It is used to predict the penetration of drugs/compounds through the skin. For the maximum transdermal transport rates, Jm, predicted by the equation (1): Where K p is the skin permeability obtained from QPlogK p , MW is molecular weight and S is the aqueous solubility obtained from QPlogS and Jm is the prediction of the maximum transdermal transport rate in µg cm -2 hr -1 .
Another parameter which is used for prediction of availability of drugs for their target is Prediction of binding to human serum  albumin. QPlogKhsa is used to predict the ability of compounds to bind to the plasma proteins such as globulins, human serum albumin, lipoprotein and glycoprotein. The plasma proteins binding ability is directly influencing the drug efficacy, the distribution of drug through the blood stream and the availability of drugs for their target. The plasma-proteins binding tendency of a drug has an inverse relation to the target availability. Therefore, the less degree of plasma-protein binding is desirable for designing drug. The recommended range for QPlogKhsa is -1.5-1.5 for 95% of known drugs.
The #metab parameter is known as the number of likely metabolic reactions which is necessary for determining the level of accessibility of compounds to their target sites after entering into the blood stream. The recommended range of #metab is 1-8.
QPlogHERG is an important parameter to predict cardiac toxicity of compounds. It is used to predict the IC50value for blockage of Human Ether-a-go-go-Related Gene Potassium (HERG K + ) channel. HERG K + plays role in the electrical activity of the heart during systolic and diastolic periods by encoding the potassium ion (K + ) channel. This channel has also modulating function in the nervous system [28]. The Blockage of HERG K + channel is potentially toxic for the cardiac and nervous system [29]. The recommended range of IC50 values for blockage of HERG K + channels is QPlogKhsa >-5.
The ADME investigation of Sameathxanthone A showed that all parameters except #metab score10, CIQPlogS score 6.6 and QPlogeHERG score 5.9 were in the recommended range of druglikeness, pharmacokinetics and metabolism criteria. Among all, only one of Xanthone derivatives, 1-isomangostin, showed as same ADME profiling as Dihydroisocoumarin 1 and possessed all ADME parameters within the recommended range. However, 1-isomangostin showed very low the affinity value, around -1.37, for inhibition of Aromatse. Figure 6 depicted the comparison of ADME criteria between Dihydroisocoumarin 1 and Sameathxanthone A.

Toxicity
Some important parameters for toxicity investigation have been predicted via online TOPKAT approaches of Accelrys Environmental Chemistry and Toxicology Workbench. The carcinogenicity of compounds have been predicted based on structural similarity between compounds and structures available in both the FDA (U.S. Food and Drug Administration) and NTP (National Toxicology Program) databases for male and female of rat and mouse (FR, FM, MR and MM). According to the predicted carcinogenicity based on the FDA database, Sameathxanthone A is non-carcinogenic for both female and male mouse and male rat, but weak carcinogen for a female rat with probability 0.34, and based on NTP database, it is safe for both male and female rat and female mouse but carcinogen for male mouse with probability 0.66. It is not skin and ocular irritant, but strong skin sensitizer with probability 0.83. According to WOE (weight of evidence for rodent carcinogenicity), it showed no carcinogenic activity in rodents. Based on the DTP (developmental toxicity potential), it was toxic, and got a positive discriminant score with the probability 0.62 for the toxicity, but showed no genotoxicity or mutagenicity potential. The Rat Oral LD50 value evaluated as 1.32 gm/kg, which was within the Optimum Prediction Space (OPS) and indicated the higher safety of this compound. The half of effective concentration values for Daphina Magna model indicated that the EC50 value of Sameathxanthone A, 0.22 mg/l, was within the recommended range. The half of the lethal concentration value (LC50) was around 0.0001 g/l and the half of the tolerance doses for carcinogenicity by feeding for mouse and rat were 31.003 and 93.278 mg/kg_body_weight/day respectively. It is liable for degradability through aerobic-biodegradability.

Conclusion
In this paper, Schrodinger 2009 and online Topkat approaches of Accelry were employed for in silico investigation of the pharmacodynamics and pharmacokinetics of the 12 xanthone based molecules with target "Aromatase" in breast cancer. The affinity value of 12 Xanthone derivatives has been compared with each other and with the affinity value of Dihydroisocoumarin 1 (IUPAC ID: [(3R, 4R)-(-)-6-methoxy-1-oxo-3-pentyl-3,4-dihydro-1H-isochromen-4-yl acetate]) as a selective Aromatase inhibitor. The result indicated that Sameathxanthone A with -7.96 kcal/mol and the binding energy -38.46 possessed significantly higher docking score than the rest of Xanthones and Dihydroisocoumarin 1and it can be introduced as a candidate for breast cancer treatment. It showed a perfect toxicity profile, but still further study is necessary to improve the ADME properties of Sameathxanthone A.

Conflict of Interest
There is no conflict of interest among the authors.