In Vitro and In Silico Analysis of the Anti-cancer Effects Eurycomanone and Eurycomalactone From Eurycoma longifolia.

Eurycomanone and eurycomalactone are known quassinoids present in the roots and stems of Eurycoma longifolia. These compounds had been reported to have cytotoxic effects, however, their mechanism of action in a few cancer cell lines have yet to be elucidated. This study was aimed to investigate the anti-cancer effects and mechanisms of action of eurycomanone and eurycomalactone in cervical (HeLa), colorectal (HT29) and ovarian (A2780) cancer cell lines via Sulforhodamine B assay. Their mechanism of cell death was evaluated based on Hoechst 33342 assay and in silico molecular docking toward DHFR and TNF-α as putative protein targets. Eurycomanone and eurycomalactone exhibited in vitro anti-cancer effects manifesting IC50 values of 4.58 ± 0.090 µM and 1.60 ± 0.12 µM (HeLa), 1.22 ± 0.11 µM and 2.21 ± 0.049 µM (HT-29), and 1.37 ± 0.13 µM and 2.46 ± 0.081 µM (A2780) respectively. They induced apoptotic cancer cell death in dose- and time-dependent manners. Both eurycomanone and eurycomalactone were also predicted to have good inhibitory potential as demonstrated by the docking into TNF-α with binding affinity of -8.83 and -7.51 kcal/mol, respectively, as well as into DHFR with binding affinity results of -8.05 and -8.87 kcal/mol, respectively. These results support the evidence of eurycomanone and eurycomalactone as anti-cancer agents via apoptotic cell death mechanism that could be associated with TNF-α and DHFR inhibition as among possible protein targets.

Keywords:

Subject: Medicine and Pharmacology - Medicine and Pharmacology

1. Introduction

Cancer is the second leading cause of morbidity and mortality globally, with colorectal, cervical and ovarian cancers among the top leading causes of cancer-related deaths worldwide [1]. Cancer death rates continue to rise globally from 7.6 million deaths in 2005 to 10 million deaths in 2020, which is nearly one in six deaths as reported by WHO [2]. Many types of cancers are asymptomatic during their early stages, which can make it difficult to detect them until they have advanced to a later stage. Late detection of these cancers normally leads to poor prognosis after treatments in many cancer patients. The main modality of treatment for treating advanced cancer is by using chemotherapy [3]. Currently, 90% of failures in the chemotherapy are during the invasion and metastasis of cancers related to drug resistance [4], and due to various toxicity effects [5]. Cancer is also consisted of more than 100 types due to their heterogeneity [6] but no single treatment or drug can treat all of these cancers, thus different drugs and treatment combinations are needed. Thus, the searches for new anti-cancer agents are still on-going around the world to overcome these issues. Natural products, particularly those derived from plants, continue to be important sources of developing new anti-cancer agents. The discovery of compounds such as taxanes and vinca alkaloids, which are respectively derived from the Pacific Yew tree and the Madagascar periwinkle plants that have been used to treat breast, leukemias, lymphomas and other various types of cancers, have motivated scientists around the world to search for anti-cancer agents from plants. From 1981 to 2019, it was estimated that 64.9% of either natural products, or synthetic variations from novel structures of natural products, had been developed into cancer drugs [7].

Eurycoma longifolia is a popular herbal plant species belonging to the Simaroubaceae family that is widely distributed in Southeast Asian countries. It is an unbranched tree that can grew up to 8 - 12 m tall that is crowned by an-umbrella-like rosette of leaves [8]. Its yellow roots are bitter in taste and famously consumed in capsulated form or beverage for health maintenance including as energy booster. The plant, also known locally as ‘tongkat Ali’ has been used traditionally to treat fever, wounds, ulcer, syphilitic sores, tonic after childbirth, to name a few [9]. E. longifolia especially for its roots had also been reported to have various medicinal properties including anti-bacterial, cytotoxic/ anti-cancer, anti-malarial, anti-ulcer, anti-parasitic and anti-pyretic [10]. The major compounds found in the roots of E. longifolia belonging to quassinoids (degraded triterpenes) are among the compounds known to contribute to various medicinal effects based on in vitro or/and in vivo studies. Eurycomanone had been reported to have anti-parasitic [11], anti-malaria [12] and anti-ulcer [13] properties, whilst, eurycomalactone had been reported to have anti-malaria [14], anti-inflammatory [15] and anti-viral [16].

Apoptosis (programmed cell death) plays critical role in removing damaged or abnormal cells from the body preventing them from dividing and growing uncontrollably by committing cell suicide. Apoptosis involves regulation and dysregulation of certain genes and proteins and these proteins often classified as pro-apoptotic and oncoproteins, respectively. However, when apoptosis is malfunctioned, these damaged or abnormal cells may continue to divide and multiply, leading to the development of tumor or cancer [17,18,19]. One of the strategies for finding new anti-cancer drug candidates is to search for compounds that target key proteins/genes in cancer cells that can re-induce the mechanism of apoptosis. Among the proteins which are involved in apoptosis directly include death receptors such as tumor necrosis factor alpha (TNF-α), Fas ligand, the caspases, Bcl-2 family proteins, to name a few. TNF-α is a multifunctional cytokine that plays a key role in cell proliferation, metabolism, differentiation, and survival [20,21,22,23,24]. It has also been identified as a potential therapeutic target for cancer treatment due to its ability to inhibit cancer development by inducing apoptosis [20,21,22,23,24].

In addition, there are also many proteins that are indirectly involved in the regulation apoptosis but also important in cell growth and proliferations such as dihydrofolate reductase (DHFR), heat shock proteins (HSPs), kinases and phosphatases, signal transduction proteins such as Ras and Raf. DHFR is one of the key proteins that involved in folate metabolism pathway that indirectly lead to apoptosis [25]. Folate metabolism is one of the well-known targets for cancer chemotherapy due to its role in nucleic acid synthesis [26,27,28]. It catalyzes the reduction of folate to tetrahydrofolate and maintains reduced cellular folate pools [26,27,28]. Inhibition of DHFR enzyme reduces the amount of tetrahydro-folate that is required for the synthesis of pyrimidines and purines, which are necessary for one-carbon metabolism reactions [29,30,31]. As a result, inhibition of DHFR can halt the synthesis of RNA and DNA, leading to the death of cancer cells via apoptosis [25,29,30,31].

As reviewed by Nurhanan et al. [32], eurycomanone was also reported to have anti-cancer effects against breast cancer, colon adenocarcinoma, fibrosarcoma, lung cancer cells and melanoma whereas eurycomalactone was reported to have anti-cancer effects against murine lymphocytic leukemia, epidermoid, melanoma, breast, lung, and colon cancer cells. Despite their reported anti-cancer effects, the mechanism of action of these compounds in several cancer cell lines such as ovary, colorectal and cervical has not been fully elucidated. Previously, 68 quassinoid analogues were designed and docked against DHFR via molecular docking [33]. Other quassinoids such as bruceantin and brusatol were also found to inhibit DHFR in P388 leukemia cells [34]. In another studies, brusatol was also reported to inhibit TNF-α [35]. In our previous proteomic analysis, we found that [33]. TNF-α was involved in apoptosis when treated with 9-methoxycanthin-6-one, a compound also isolated E. longifolia [32]. DHFR and TNF-α are also known to be involved in apoptosis when the cancer cells were treated with known chemotheraphy drug, methotrexate [36,37]. It is therefore suggested that DHFR and TNF-α are among attractive targets for the search and development of anti-cancer agents based on the apoptotic cell death mechanism [20,27,38]. Hence, the aim of these studies was to evaluate the in vitro anti-cancer effects and the mechanisms of action of eurycomanone and eurycomalactone via apoptosis. The interactions of both compounds with DHFR and TNF-α and their druglikeness activities were also investigated via in silico studies.

2. Results and Discussion

2.1. Isolation and characterization of eurycomalactone and eurycomanone

Eurycomanone (1), was isolated from E. longifolia as a white powder. Its molecular formula C₂₀H₂₄O₉ was established by positive-ion HREIMS [M+H]⁺409.3851,(calcd C₂₀H₂₄O₉ for m/z 408.403. Comparing its ¹H and ¹³C-NMR spectra (Table 1) with those of (2), we found that chemical shift of methyl protons for 4 and 8-CH₃ was increased to δ_H 1.81 and 2.00 respectively. The presence of ortho protons in C-14 and C-15 were displayed δ_H 3.25 and δ_H 5.25 respectively. 20 signals were displayed in its ¹³C-NMR spectrum including two carbonyl carbons at δ_C 197.89 (C-2) and δ_C 174.35 and four methylene carbons at δ_C 162.98 (C-3), 126.48 (C-4), 119.80 (C-13) and 108.71 (C-13’).

Eurycomalactone (2), was a white powder and its molecular formula was determined as C₁₉H₂₄O₆ by HREIMS at m/z 349.1647 [M+H]⁺, ¹H and ¹³C-NMR spectra. The ¹H-NMR spectrum of (1), showed signals at δ_H 6.13 (1H, s) corresponding to vinyl proton at C-3. Four signals resolved at δ_H 1.64, 1.96, 1.27 and 1.18 (3H, s) indicated methyl protons in compound (1). The ¹³C NMR peak assignments showed four methyl carbons resonated at δ_C 23.64, 21.97, 12.19 and 32.33. Meanwhile signals at δ 176.31, 197.41 and 205.56 attributed to carbonyl carbon (C-2, C-3 and C-15).

The known constituent (1) and (2) were identified by comparison of its spectral data ¹H and ¹³C NMR and MS with those reported in the literature [13,39,40]. The structure of eurycomanone (1) and eurycomalactone (2) are shown in Figure 1.

2.2. Percentage of Cells Viability and IC₅₀ values

The in vitro anti-cancer effects of eurycomanone and eurycomalactone were evaluated against ovarian (A2780), cervical (HeLa), colorectal (HT29) cancer cells, as well as normal cardiomyocyte (H9C2) and liver (WRL-68) cell lines. Cisplatin and methotrexate were also evaluated on its in vitro anti-cancer effects as for comparison studies. Both drugs had been reported to treat various type of cancers including ovarian, cervical and colorectal cancers [41,42,43]. Both eurycomanone and eurycomalactone gave significant in vitro anti-cancer effects, with IC₅₀ values ranging from 1.22 ± 0.11 µM to 4.58 ± 0.090 µM (for eurycomanone) and 1.60 ± 0.12 µM to 2.46 ± 0.081 µM (for eurycomalactone) as illustrated in Figure 2 and tabulated in Table 2. Cisplatin exerted comparable in vitro anti-cancer activities with the range of IC₅₀ between 1.38 ± 0.037 µM to 1.77 ± 0.018 µM, whereas, methotrexate exerted better in vitro anti-cancer activities with the range of IC₅₀ between 0.016 ± 0.00050 µM to 0.094 ± 0.0043 µM respectively (Figure 2, Table 2). Eurycomanone, eurycomalactone and methotrexate showed less cytotoxic effects in cardiomyocyte H9c2 normal cell line as opposed to liver WRL-68 normal cell line (Table 2). Clinically, cisplatin was reported to have several side effects including cardiotoxicity, nephrotoxicity and neurotoxicity [5]. Whereas, methotrexate was also reported to have adverse effects (e.g. hepatotoxicity, nephrotoxicity, gastrointestinal toxicity and death due to infections and hemorrhage) when given to cancer patients [44].

From our current studies, the in vitro anti-cancer analysis had shown that eurycomanone and eurycomalactone were in the potent range since the IC₅₀ values were less than 50 µM [45]. To the best of our knowledge, this is the first report that revealed eurycomalactone had in vitro anti-cancer effects against ovarian cancer cells. Other studies had also reported that eurycomanone and eurycomalactone gave in vitro anti-cancer effects against other cancer cell lines. Eurycomanone was reported to have anti-cancer effects in breast [46,47], colon, fibrosarcoma, lung and melanoma cancer cell lines [46] with IC₅₀ values ranging from 0.49-35 µM. Whereas, eurycomalactone was reported to have anti-cancer effects in murine lymphocytic leukemia (P388) and epidermoid (KB) [48], lung cancer (A-549), breast cancer (MCF-7) [47] and colon (26-L5), melanoma (B16-BL6), lung (LLC and A549) cancer cell lines [49] with IC₅₀ values ranging from 0.57-23.25 µM. However, the informations on the mechanisms of action of eurycomanone and eurycomalactone in killing the cancer cells are still lacking.

2.3. Apoptotic Effects of Eurycomanone and eurycomalactone via Hoechst 33342 Assay

Apoptotic effects of eurycomanone and eurycomalactone against A2780, HT-29 and HeLa cancer cell lines were analyzed after performing Hoechst 33342 assay. These cell lines were also treated with these compounds at different concentrations and incubation times to preliminary evaluate their pharmacodynamic effects in inducing apoptosis. Apoptosis is a programmed cell death in which the cells will undergo distinct changes in its morphology which normally started with chromatin condensation, DNA fragmentation, cell shrinkage, cell membrane blebbing and finally formation of apoptotic bodies. Hoechst assay is commonly used to stain the apoptotic cells in fluorescence blue [50].

As shown in Figure 3, Hoechst dye had stained the occurrence of chromatin condensation, dense chromatin at the periphery of the nucleus and apoptotic bodies upon treatment with eurycomanone, whereas non-apoptotic cells had spherical nuclei and evenly distributed chromatin. The percentage of apoptotic indices were shown in Table 3. Eurycomanone had significantly increased the percentage of apoptotic index in A2780, HT-29 and HeLa cells in a concentration- and time-dependent manners especially at 24 h and 48 h incubation times and at the highest concentration of eurycomanone tested.

As shown in Figure 4, eurycomalactone had also induced apoptosis A2780, HT-29 and HeLa cells in which different stages and hallmark signs of apoptosis were clearly shown such as chromatin condensation, dense chromatin at the periphery of the nucleus and formation of apoptotic bodies. The percentage of apoptotic indices were shown in Table 4. Eurycomalactone had significantly increased the percentage of apoptotic index in HT-29 and HeLa cells in a concentration- and time-dependent manners especially at the highest concentration of eurycomalactone tested.

2.4. Molecular Docking Analysis

Molecular docking is a widely used computational technique in modern drug design for predicting drug-receptor interactions and identifying potential inhibitors. In this study, we conducted a comparative analysis of the free binding energy and binding interactions of eurycomanone and eurycomalactone with the co-crystallized ligands in the active binding sites of TNF-α (PDB ID: 2AZ5) and DHFR (PDB ID: 5HQY) (refer to Table 5, Figure S1, and Figure 5). Our aim is to gain insights into the potential mechanisms of action of eurycomanone and eurycomalactone against TNF-α and DHFR. To validate the docking process, we first re-docked the co-crystallized (original) ligands as controls into the active sites of DHFR and TNF-α. The resulting score energies were -7.93 kcal/mol and -8.19 kcal/mol, respectively, with small root mean square deviations (RMSDs) of 0.98 Å for TNF-α and 0.62 Å for DHFR (Table 5). An RMSD value of ≤ 2.0 Å is typically considered acceptable in the literature [51,52,53,54,55]. Based on these results, we applied the same docking parameters to dock eurycomanone and eurycomalactone.

Figure S1(a, b) provides valuable insights into the interactions between co-crystalized ligands and the active binding sites of TNF-α and DHFR. In Figure S1a, the original ligand complexed into the active binding site of TNF-α (PDB ID: 2AZ5) is observed to form a single hydrogen bond with the tyrosine residue TYR151 at a distance of 2.45 Å. In addition, the ligand forms pi-sigma and pi-alkyl interactions with TYR119 and TYR59, which are known to be important in protein-ligand recognition and may contribute to the inhibitory activity of the compound. These findings are consistent with those reported in [36]. In Figure S1b, the original co-crystalized ligand in PDB ID: 2AZ5 adopts a bent conformation and forms four hydrogen bonds (H-bonds) with ILE7 (3.09 Å), GLU30 (2.59 Å), ASN64 (3.44 Å), and VAL115 (3.15 Å), with relatively weak to moderate interactions. The inhibitor also establishes hydrophobic contacts with several active site residues of DHFR, including ALA9, LEU22, PHE31, PHE34, and ILE60. These interactions likely contribute to the inhibitory activity of the co-crystalized ligands. A comprehensive understanding of these interactions can shed light on the mechanism of action of TNF-α and DHFR inhibitors and facilitate the development of more effective drugs targeting these enzymes.

Table 5 presents the free binding energies of eurycomanone and eurycomalactone, compared to the co-crystallized ligand, within the active binding site of TNF-α (PDB ID: 2AZ5). Remarkably, the calculated binding free energies for eurycomanone and eurycomalactone are -8.83 kcal/mol and -7.51 kcal/mol, respectively, which are substantially similar to that of the co-crystallized ligand. Figure 5(a,b) show detailed 2D and 3D molecular interaction analyses of eurycomanone and eurycomalactone with the active site of TNF-α. Eurycomanone is found to form three hydrogen bonds, with two of them formed between the hydroxyl group and carbonyl group of the 6-hydroxy-3-methylcyclohex-2-en-1-one ring with GLY121(A) at a distance of 2.24 Å and TYR151(B) at a distance of 1.92 Å. The third hydrogen bond is formed between the hydroxyl group of the 2-methylenecyclohex-3-en-1-ol ring and LEU120(A) at a distance of 2.02 Å. Eurycomanone also engages in pi-alkyl interactions with TYR59(A), TYR119(A), and TYR119(B) (See Figure 5a). In the same pattern, eurycomalactone interacts with the same residues through pi-alkyl interactions with TYR59(B), TYR119(B), and TYR151(B), and its carbonyl group of the 6-hydroxy-3-methylcyclohex-2-en-1-one ring forms a single hydrogen bond interaction with TYR119(A) at a distance of 2.01 Å. These detailed molecular interaction patterns and the Ki values suggest that both eurycomanone and eurycomalactone may serve as promising hit candidates against TNF-alpha, with their strong binding affinity likely translating to a potent inhibitory effect on TNF-α activity.

In the active binding site of DHFR (PDB ID: 5HQY), both eurycomanone and eurycomalactone exhibit free binding energies similar to the co-crystallized ligand. The calculated binding free energies for eurycomanone and eurycomalactone, which are -8.83 kcal/mol and -7.51 kcal/mol, respectively, suggest a promising potential for these compounds to bind to DHFR within its active binding site (Table 5).

Figure 5(c,d) illustrate the 2D and 3D molecular interactions of eurycomanone and eurycomalactone with DHFR. Eurycomanone is observed to establish two intermolecular hydrogen bonds with crucial residues in the active binding site of the enzyme (Figure 5c). The first hydrogen bond is formed between the carbonyl group of the 3-hydroxytetrahydro-2H-pyran-2-one ring and the guanidine group of ARG70 at a distance of 2.56 Å, while the second hydrogen bond is formed between the hydroxyl group of the 3-hydroxytetrahydro-2H-pyran-2-one at a distance of 1.67 Å. In the case of eurycomalactone, the 6-hydroxy-4-methylcyclohex-2-en-1-one also forms two intermolecular hydrogen bonds with TRP24 at a distance of 2.47 Å and with GLU30 at a distance of 1.91 Å (Figure 5d). Along with these hydrogen bonds, both eurycomanone and eurycomalactone are stabilized in the active site of DHFR through hydrophobic interactions with several residues, including ILE7, ALA9, LEU22, TRP24, PHE31, PHE34, and ILE60 for eurycomanone, and VAL8, ILE16, LEU22, PHE34, and VAL115 for eurycomalactone. These interactions bear substantial resemblance to those observed with the original co-crystallized ligand, indicating that eurycomanone and eurycomalactone may inhibit DHFR through similar mechanisms or contribute to the inhibitory activity of the compounds.

From these molecular docking analyses, eurycomanone and eurycomalactone had the potential to inhibit both TNF-α and DHFR. TNF-α and DHFR also had been reported to be involved in inhibiting cell proliferations via apoptosis by others, making them among potential targets for the treatment of cancer and other diseases [22,23,27,28]. Some anti-cancer drugs had been developed and targeted these TNF-α and DHFR as part of its mechanism in killing the cancer cells. For example, a few anti-cancer drugs that already approved by the FDA and others undergoing clinical trials include methotrexate [56,57], pemetrexed [58] and pyrimethamine [59]. Whereas, TNF-α is a cytokine involved in inflammation and targeted for the treatment of autoimmune, inflammatory disorders [20,28] and cancer [60]. A few anti-cancer drugs that targeted TNF- α include doxorubicin [61], melphalan [60] and pembrolizumab [62]. To the best of our knowledge, this is the first finding that eurycomanone and eurycomalactone were found to target TNF-α and DHFR via this molecular docking analysis. Nevertheless, further investigations shall be conducted to validate the inhibitory effects of eurycomanone and eurycomalactone against TNF-α and DHFR.

2.5. Lipinski’s Rule and ADMET of eurycomanone and eurycomalactone

It is known that Lipinski's Rule of Five (RoF) is utilised to evaluate the potential of a drug to be orally bioavailable [63]. The assessment relies on molecular properties, including molecular weight, number of hydrogen bond donors and acceptors, and lipophilicity [63]. It has become essential to predict the pharmacokinetic properties of a lead compound to assess its druggable potential before entering the drug development phase [64,65]. In this study, we assessed eurycomanone, eurycomalactone and methotrexate (control) using the ADMETlab 2.0 web service tool to evaluate RoF and their pharmacokinetic properties [66]. The results of the predicted values for eurycomanone, eurycomalactone, and methotrexate are presented in Table 6 and Table 7.

According to RoF, a compound is more likely to be orally active if it has no more than one violation of the following criteria: Log P is less than 5; molecular weight is less than 500 Da; hydrogen bond donor is less than 5; and hydrogen bond acceptor is less than 10 [63]. In Table 6, the predicted values of these properties for three compounds, eurycomanone, eurycomalactone, and methotrexate are presented. Eurycomanone has a molecular weight of 408.14 g/mol, 9 hydrogen bond acceptors, 4 hydrogen bond donors, and a partition coefficient (logP) of 0.215. Eurycomalactone has a molecular weight of 348.16 g/mol, 6 hydrogen bond acceptors, 1 hydrogen bond donor, and a logP of 0.655. Both compounds satisfy the Lipinski's Rule of Five, as their molecular weights are less than 500, their logP values are within the acceptable range of 0-3, and their hydrogen bond acceptor and donor counts are less than 10 and 5, respectively.

On the other hand, methotrexate has a molecular weight of 454.17 g/mol, 13 hydrogen bond acceptors, 7 hydrogen bond donors, and a logP of -2.747. Methotrexate violates the Lipinski's Rule of Five since it has more than 10 hydrogen bond acceptors and more than 5 hydrogen bond donors. Therefore, methotrexate might have poor oral bioavailability. The Lipinski's Rule of Five predicted that eurycomanone and eurycomalactone are orally bioavailable, while methotrexate might have poor oral bioavailability due to its molecular properties. However, it should be acknowledged that this rule alone is not sufficient to guarantee the efficacy or safety of a drug, as it overlooks other factors that may affect the pharmacokinetic and pharmacodynamic profiles of a compound.

To address these limitations, it is necessary to complement the Rule of Five with more comprehensive ADMET evaluations that evaluate the absorption, distribution, metabolism, excretion, and toxicity of drug candidates. In this study, we utilized the ADMETlab 2.0 web service tool to evaluate the pharmacokinetic properties of the compounds under investigation, including eurycomanone, eurycomalactone, and methotrexate (See Table 7), which were also analyzed using Lipinski's Rule of Five.

By combining the results of the Lipinski rule with the ADMET evaluation, we can obtain a more accurate and reliable prediction of the drug-likeness, pharmacokinetics and safety of the compounds, which is crucial for guiding the drug discovery and development process. Therefore, the combination of these two methods is highly recommended for identifying promising drug candidates with optimal pharmacokinetic properties, and this approach has the potential to improve the efficiency and success rate of the drug development process.

Caco-2 permeability is an important parameter when determining oral absorption and permeability in the early stages of drug development, and the ideal value is greater than −5.15 cm/s [54,67]. Eurycomanone and methotrexate showed low permeability, while eurycomalactone displayed higher permeability compared to the rest. An in vivo test on bioavailability studies had reported that eurycomanone is poorly bioavailable when given orally [68] which supported this in silico finding. Whereas methotrexate is reported to be consumed either orally or intravenously [69,70]. To the best of our knowledge, there are no in vitro or in vivo studies being reported on pharmacokinetic properties on eurycomalactone. Eurycomanone, eurycomalactone and methotrexate, all showed high levels of human intestinal absorption satisfaction (HIA% with a value of 30%), which is a crucial parameter associated with human intestinal absorption [54,71], affecting the way compounds pass through biological membranes under the influence of physicochemical properties. In the body, plasma protein binding plays a significant role in the dynamics of compounds [66,72]. This phenomenon is well-known as plasma protein binding percent (PPB%). The predicted data suggest a low protein binding potential for eurycomanone and eurycomalactone compared to the reference value of 90%. Drugs with low protein binding may have a high therapeutic index [66]. Another well-known parameter is the Blood-Brain Barrier (BBB), which facilitates the selective transfer of drug molecules between the blood and the brain parenchyma [54,66,73]. As predicted, eurycomanone and eurycomalactone were not able to cross the BBB as compared to methotrexate, which adds to the safety profile of these compounds.

A drug's fate can be determined by its metabolism, which is characterized by the enzymatic modification or degradation of its molecules according to its therapeutic response [74]. Cytochrome P450 (CYP) enzymes are essential for metabolism of drugs. Among 57 functional CYPs, the isoforms belonging to CYP1, CYP2 and CYP3 are responsible for the metabolism of 80% of clinical drugs [75]. Table 7 shows eurycomanone, eurycomalactone, and methotrexate showed no inhibition on CYP metabolism on other drug. An in vitro assay on eurycomanone that were tested on these cytochromes also had shown that there was no inhibition [76], thus indicating that eurycomanone may not interact with other drugs. The information on half-life and clearance of a drug candidate has important implications for preparing dosing regimen clinically. Half-life is defined as the time required for the concentration of a drug (typically in blood or plasma) to reduce to half of its initial value when the concentrations of the drug are in simple exponential (log−linear) decline [77]. Short half-time periods and low clearance rates were predicted from our study for eurycomanone, eurycomalactone, and methotrexate. As a leading cause of drug withdrawals from the market, drug-induced toxicity will remain a key concern for the development of novel molecules [52,54,78]. Eurycomanone and eurycomalactone showed low risk in exerting hepatotoxicity and carcinogenicity effects. Whereas, eurycomalactone had lower risk in exerting mutagenicity effect as compared to eurycomanone. Overall, eurycomanone and eurycomalactone showed good druglikeness and ADMET properties but further in vitro and in vivo studies need to be conducted to validate these predictions.

3. Materials and Methods

3.1. Isolation and Characterisation of Eurycomanone and Eurycomalactone

The dried roots of E. longifolia (900 g) was extracted in ethanol (9 L) using soxhlet extraction at 55-60°C for 18 h. The ethanol solution of extract was filtered and evaporated to yield the ethanol extract (46 g, 5.1%). The ethanol extract was fractionated using vacuum liquid chromatography on silica gel 60 (230-400 mesh) using n-hexane, n-hexane-ethyl acetate (EtOAc) (4:1, 3:2, 2:3 and 1:4) as mobile phase to obtain 4 fractions (FR1-FR4). A combined 2 fractions (FR1-FR2) (1.57 g, 0.17%) was re-fractionated by column chromatography (cc) on silica gel 60 (70-230 mesh) and eluted with n-hexane-EtOAc (9:1, 3:2, 2:3 and 1:9) yielded 4 sub-fractions (fr1-fr4). Subfractions fr2 (1.3 g) and fr3 (0.33g) were then re-chromatographed by cc to afford eurycomalactone (2) (7.5 mg, 0.001%) and eurycomanone (1) (25.1 mg, 0.002%) respectively.

3.2. Cell Culture and Treatments

The cancer cell lines used for this study were ovarian cancer (A2780), cervical cancer (HeLa) and colorectal cancer (HT29) as well as normal cardiomyocyte (H9C2) and liver (WRL-68) cell lines. All cell lines were purchased from American Type Culture Collections except A2780 from the European Collection of Authenticated Cell Cultures. These cells were sub-cultured in Dulbecco’s Modified Eagle’s medium (Sigma) supplemented with 10 % fetal bovine serum (Sigma, USA), 1% penicillin-streptomycin (Sigma, USA), 1% amphotericin B (Sigma, USA) and 1% gentamicin (Sigma, USA). The cells were seeded in each well of the 96 well plates and incubated in a humidified incubator at 37 °C and 5% carbon dioxide in air for 24 h. Each cell line was then treated with the eurycomanone and eurycomalactone at five different concentrations (0.08, 0.4, 2, 10 and 50 µM) in triplicate. Cisplatin and methotrexate (Sigma, USA), known chemo-drugs, were also treated on these cell lines at five different concentrations (0.08, 0.4, 2, 10, 50 µM for cisplatin, 0.0008, 0.004, 0.02, 0.1, 0.5 µM for methotrexate) as for the comparative studies. The treated cells were then incubated in the same incubator with the mentioned conditions for 72 h. The experiment was repeated at least three times.

3.3. Cells Viability Assay

Sulforhodamine B (SRB) assay [79,80] had been performed after the treated cells were incubated for 72 h. Briefly, 50 μL of ice cold tricholoroacetic acid (TCA) was added to each well and allowed to stand for 30 min at room temperature, followed by rinsing each well with tap water. Then, 100 μL of 0.4% SRB was added to each well to stain living cells for 30 min followed by a rinse with 1% acetic acid. Finally, 100 μL of Tris buffer was added to each well and the optical density (OD) of the treated and non-treated cells were read at 492 nm with a Magellan V.4 microtiter plate reader (Tecan, Switzerland). The percentage of cell viability was calculated based on (OD_492nm of the treated cells/ OD_492nm of the non-treated cells) × 100. The IC₅₀ values were determined from the dose–response curve of percentage of cell viability versus the concentration of the compounds (μM). Cells viability assay for each treatment was performed in triplicate in at least three independent experiments, and the IC₅₀ values are given as the mean ± SEM.

3.4. Apoptotic Hoechst 33342 Assay

A2780, HT-29 and HeLa cancer cell lines were chosen for evaluating the apoptotic effects of eurycomanone and eurycomalactone. The morphological changes of cells undergoing apoptosis were visualised under a fluoresecence microscope after performing the Hoechst 33342 assay. Briefly, 10 × 10⁴ cells were seeded for 500 µL/ well of four-Labtek® Chamber Slides (Thermo Fisher Scientific, Waltham, Massachusetts, USA) and incubated in 5% carbon dioxide in air for 24 h for cells’ attachment. The cells were treated with the concentrations of IC₅₀/5, IC₅₀ and IC₅₀ x5 values of eurycomanone and eurycomalactone in different wells, and further incubated for 6 h, 24 h and 48 h, respectively. Non-treated cells were also included in the experimental design, which acted as a negative control in this study. Following this, the media from each well were discarded and cells were fixed with 4% (w/v) paraformaldehyde for 30 min. The cells were then washed with cold phosphate buffer solution (PBS) prior to staining using Hoechst 33342 according to the method of [32]. Quantitative assessment of apoptotic cells were determined by counting apoptotic nuclei from five randomly chosen areas corroborating visual impression based on the morphology shown in Figure 3 and 4. The apoptotic index (number of apoptotic nuclei per total nuclei x 100) was expressed as percentage of mean ± S.E.M and analysed using ANOVA followed by Tukey’s multiple comparison test.

3.5. Molecular Docking simulation and ADMET predictions

The Protein Data Bank database [81] was used to retrieve the human crystal structures of dihydrofolate reductase (DHFR) (PDB ID: 5HQY) [82] and tumor necrosis factor alpha (TNF-α) (PDB ID: 2AZ5) [83]. All water molecules and heteroatoms have been eliminated using the Biovia Discovery Studio Visualizer (San Diego, CA, USA, 2019) [84]. With MODELLER 9.18 [85], all missing residues in the crystal structures have been added and refined. Using the PDB2PQR web service (https://pdb2pqr.poissonboltzmann.org/pdb2pqr, accessed on October 14th, 2022), additional treatments were performed on the crystal structures, including reconstructing missing atoms, assigning atomic charges, and radii using the SWANSON force field (AMBER ff99 charges with optimized radii) [53,54,55,86]. At pH 7.40, PROPKA3 [87], which is the most commonly used empirical pKa predictor, was used to figure out the protonation states of the ionizable groups in the crystal structures. After the protonation states were assigned, the proteins were uploaded to MolProbity (http://molprobity.biochem.duke.edu/, accessed October 14th, 2022) to correct bad contacts, hydrogen atom additions, and flipping of HIS, GLU, and ASN residues [53,54,55,88].

The co-crystalized ligands 5-methyl-6-[(2,3,4-trifluorophenyl)sulfanyl]thieno[2,3-d]pyrimidine-2,4-diamine) and 6,7-dimethyl-3-[[methyl-[2-[methyl-[[1-[3 (trifluorome-thyl)phenyl]indol-3-yl]methyl]amino]ethyl]amino]methyl]chromen-4-one) were taken from the crystal structures of 5HQY.PDB and 2AZ5.PDB, respectively. While the chemical structures of the derivatives were downloaded from the PubChem database [89], identification references were provided for eurycomanone (PubChem ID: 13936691) and eurycomalactone (PubChem ID: 441793). Then, the derivatives (eurycomanone and eurycomalactone) were subjected to energy minimization using the Molecular Mechanics 2 (MM2) force field by PerkinElmer Chem3D 17.1 (PerkinElmer, MA, USA) [52,53].

In this study, we employed molecular docking simulations to investigate the affinity of the chemical functional groups of ligands with significant amino acids for binding to the active sites of TNF-α and DHFR [90]. The inhibition constant (Ki) is a critical measure of the affinity of a ligand, such as a drug, for its target receptor. Ki is defined as the concentration of the ligand that occupies 50% of the receptors in a population [91,92,93]. Pharmacologists frequently use Ki values to assess the potency of a ligand and its potential for interacting with a receptor [91,92,93]. A lower Ki value indicates that the ligand has a higher affinity for the receptor and thus a greater potential for receptor activation or inhibition [91,92,93]. This analysis may provide insights into the molecular interactions and binding affinities of these compounds with the proteins.

AutoDockTools 1.5.6 (The Scripps Research Institute, La Jolla, CA, USA [90]was used to add polar hydrogens and Kollman charges to the proteins, while Gasteiger charges were assigned to the co-crystallized ligands and the derivatives, and all were saved in PDBQT format. In addition, the flexibility of the ligand (active rotatable bonds) was preserved. The parameters were set as follows: Grid box size for DHFR is 40*40*40, grid spacing is 0.375, and coordinates are x = -4.39, y = 17.11, and z = 25.02. TNF-α has a grid box size of 40*40*40, a grid spacing of 0.375, and coordinates of x = -19.16, y = 74.45, and z = 33.83. These coordinates were centered on the active binding sites and were saved in grid parameter files (GPFs). Proteins were set as rigid, while ligands were set as flexible. The number of genetic algorithm runs was 100, the number of population sizes was 150, the maximum number of evaluations was 2,500,000 (medium), and the maximum number of generations was 27,000. The Lamarckian genetic algorithm was selected to accomplish this process, and the remaining parameters were retained as defaults and saved in docking parameter files (DPFs). The docking processes were simulated using AutoDock 4.2 [94]. The Biovia Discovery Studio Visualizer [84], which allows for both 2D and 3D visualization, was used to visualize the molecular interactions between the ligands and the active site of proteins.

ADMETlab 2.0 web service tool was used to predict the druglike-ness via Lipinski’s rule of five (RoF) and pharmacokinetic properties (https://admetmesh.scbdd.com/service/evaluation/cal (accessed on October 28, 2022)), which anticipates mutagenicity (Ames test), carcinogenicity, permeability of the BBB, as well as absorption of human intestinal protein and plasma protein binding [66]. The two-dimensional chemical structures of eurycomanone and eurycomalactone were converted to the SMILES format and submitted to ADMETlab 2.0 to obtain the results. Methotrexate, an anti-cancer drug with antifolate property was included in this pharmacokinetic analysis for comparison studies.

4. Conclusions

In conclusion, both eurycomanone and eurycomalactone have shown strong in vitro anti-cancer activities and their mechanisms of action was via apoptosis in killing the cancer cells. The molecular docking study showed that both compounds can target DHFR and TNF-α by binding to their active sites and forming hydrogen bonds with the key residues. The compounds also exhibited hydrophobic interactions that enhanced their stability. In addition, eurycomanone and eurycomalactone may be suitable to be developed as drug candidates because they both obeyed the Lipinski’s rule. These two compounds are also predicted to be safe, implying that these compounds could serve as direct inhibitors against the target proteins. However, further molecular dynamics, in vitro and immuno-assays should be conducted to confirm the predicted values.

Supplementary Materials

The following supporting information of binding interactions of eurycomanone and eurycomalactone with the co-crystallized ligands in the active binding sites of TNF-α (PDB ID: 2AZ5) and DHFR (PDB ID: 5HQY) for eurycomanone and eurycomalactone can be downloaded at: www.mdpi.com/xxx/s1.

Author Contributions

For research articles with several authors, a short paragraph specifying their individual contributions must be provided. The following statements should be used “Conceptualization, N.M.Y. and H.A.W..; methodology, N.M.Y. and H.A.W.; software, H.A.W..; validation, N.M.Y, M.G.A.-T.; formal analysis, N.M.Y.; investigation, N.M.Y.; M.G.A.-T.; resources, M.H.J.; data curation, N.M.Y.; writing—original draft preparation, N.M.Y.; writing—review and editing, N.M.Y.; H.A.W.; N.J.S; visualization, N.J.S.; M.G.A.-T.; supervision, N.M.Y.; H.A.W.; project administration, N.M.Y.; H.A.W.; funding acquisition, N.M.Y. All authors have read and agreed to the published version of the manuscript.”

Funding

This research was funded by Public Services Department, Malaysia (304/PFARMASI/6501298/J116) and the government of Malaysia under 12th Malaysian plan (grant number 24-01070-2001).

Data Availability Statement

All relevant data already included in the text.

Acknowledgments

The authors are grateful to Ruzana Rabuzin for assisting in the maintenance of cell cultures.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
World Health Organization (WHO) WHO Report on Cancer: Setting Priorities, Investing Wisely and Providing Care for All. 2020.
Ashdown, M.L.; Robinson, A.P.; Yatomi-Clarke, S.L.; Ashdown, M.L.; Allison, A.; Abbott, D.; Markovic, S.N.; Coventry, B.J. Chemotherapy for Late-Stage Cancer Patients: Meta-Analysis of Complete Response Rates. F1000 Research 2015, 4, 1–14. [Google Scholar] [CrossRef] [PubMed]
Dinić, J.; Efferth, T.; García-Sosa, A.T.; Grahovac, J.; Padrón, J.M.; Pajeva, I.; Rizzolio, F.; Saponara, S.; Spengler, G.; Tsakovska, I. Repurposing Old Drugs to Fight Multidrug Resistant Cancers. Drug Resistance Updates 2020, 52, 100713. [Google Scholar] [CrossRef]
El-Awady, E.S.E.; Moustafa, Y.M.; Abo-Elmatty, D.M.; Radwan, A. Cisplatin-Induced Cardiotoxicity: Mechanisms and Cardioprotective Strategies. Eur J Pharmacol 2011, 650, 335–341. [Google Scholar] [CrossRef]
Heppner, G.H. Tumor Heterogeneity. Cancer Research 1984, 44, 2259–2265. [Google Scholar]
Newman, D.J.; Cragg, G.M. Natural Products as Sources of New Drugs over the Nearly Four Decades from 01/1981 to 09/2019. J Nat Prod 2020, 83, 770–803. [Google Scholar] [CrossRef]
Kulip, J.; Wong, K.M. Simaroubaceae. In Tree Flora of Sabah and Sarawak Volume One; Soepadmo, E., Wong, K.M., Eds.; Ampang Press Sdn. Bhd.: Kuala Lumpur, 1995; Vol. 1, pp. 421–422. [Google Scholar]
Burkill, L.H. A Dictionary of the Economic Products of the Malay Peninsula; Ministry of Agriculture and Cooperatives, Kuala Lumpur., 1966; Vol. I & II.
Rehman, S.U.; Choe, K.; Yoo, H.H. Review on a Traditional Herbal Medicine, Eurycoma Longifolia Jack (Tongkat Ali): Its Traditional Uses, Chemistry, Evidence-Based Pharmacology and Toxicology. Molecules 2016, 21, 331. [Google Scholar] [CrossRef] [PubMed]
Girish, S.; Kumar, S.; Aminudin, N.; Hashim, N.M. Comparison of Apoptotic Responses in Blastocystis Sp. upon Treatment with Tongkat Ali and Metronidazole. Scientific Reports 2021, 11, 1–12. [Google Scholar] [CrossRef]
Low, B.S.; Teh, C.H.; Yuen, K.H.; Chan, K.L. Physicochemical Effects of the Major Quassinoids in a Standardized Eurycoma longifolia Extract (Fr 2) on the Bioavailability and Pharmacokinetic Properties, and Their Implications for Oral Antimalarial Activity. Nat Prod Commun 2011, 6, 337–341. [Google Scholar] [CrossRef]
Tada, H.; Yasuda, F.; Otani, K.; Doteuchi, M.; Ishihara, Y.; Shiro, M. New Antiulcer Quassinoids from Eurycoma longifolia. Eur J Med Chem 1991, 26, 345–349. [Google Scholar] [CrossRef]
Chan, K.L.; O’Neill, M.J.; Phillipson, J.D.; Warhurst, D.C. Plants as Sources of Antimalarial Drugs. Part 3. Eurycoma longifolia. Planta Med 1986, NO. 2, 105–107. [Google Scholar] [CrossRef]
Tran, T.V.A.; Malainer, C.; Schwaiger, S.; Atanasov, A.G.; Heiss, E.H.; Dirsch, V.M.; Stuppner, H. NF-ΚB Inhibitors from Eurycoma longifolia. J Nat Prod 2014, 77, 483–488. [Google Scholar] [CrossRef] [PubMed]
Choonong, R.; Ruangdachsuwan, S.; Churod, T.; Palabodeewat, S.; Punyahathaikul, S.; Juntarapornchai, S.; Ketsuwan, K.; Komaikul, J.; Masrinoul, P.; Kitisripanya, T.; et al. Evaluating the in Vitro Efficacy of Quassinoids from Eurycoma longifolia and Eurycoma harmandiana against Common Cold Human Coronavirus OC43 and SARS-CoV-2 Using In-Cell Enzyme-Linked Immunosorbent Assay. J Nat Prod 2022. [Google Scholar] [CrossRef] [PubMed]
Haiat, S.; Billard, C.; Quiney, C.; Ajchenbaum-Cymbalista, F.; Kolb, J.P. Role of BAFF and APRIL in Human B-Cell Chronic Lymphocytic Leukaemia. Immunology 2006, 118, 281–292. [Google Scholar] [CrossRef] [PubMed]
Lim, B.; Greer, Y.; Lipkowitz, S.; Takebe, N. Novel Apoptosis-Inducing Agents for the Treatment of Cancer, a New Arsenal in the Toolbox. Cancers 2019, 11, 1087. [Google Scholar] [CrossRef] [PubMed]
Niknejad, H.; Khayat-Khoei, M.; Peirovi, H.; Abolghasemi, H. Human Amniotic Epithelial Cells Induce Apoptosis of Cancer Cells: A New Anti-Tumor Therapeutic Strategy. Cytotherapy 2014, 16, 33–40. [Google Scholar] [CrossRef] [PubMed]
van Horssen, R.; ten Hagen, T.L.M.; Eggermont, A.M.M. TNF-α in Cancer Treatment: Molecular Insights, Antitumor Effects, and Clinical Utility. Oncologist 2006, 11, 397–408. [Google Scholar] [CrossRef]
Balkwill, F. Tumour Necrosis Factor and Cancer. Nature Reviews Cancer 2009, 9, 361–371. [Google Scholar] [CrossRef]
Zidi, I.; Mestiri, S.; Bartegi, A.; Amor, N. Ben TNF-α and Its Inhibitors in Cancer. Medical Oncology 2009, 27, 185–198. [Google Scholar] [CrossRef]
Zia, K.; Ashraf, S.; Jabeen, A.; Saeed, M.; Nur-e-Alam, M.; Ahmed, S.; Al-Rehaily, A.J.; Ul-Haq, Z. Identification of Potential TNF-α Inhibitors: From in Silico to in Vitro Studies. Scientific Reports 2020, 10, 1–9. [Google Scholar] [CrossRef]
Marriott, J.B.; Westby, M.; Dalgleish, A.G. Therapeutic Potential of TNF-α Inhibitors Old and New. Drug Discov Today 1997, 2, 273–282. [Google Scholar] [CrossRef]
Neradil, J.; Pavlasova, G.; Veselska, R. New Mechanisms for an Old Drug; DHFR-and Non-DHFR-Mediated Effects of Methotrexate in Cancer Cells. Klinicka Oncologie 2012, 25, 2S87–2S92. [Google Scholar]
Srivastava, V.; Kumar, A.; Mishra, B.N.; Siddiqi, M.I. Molecular Docking Studies on DMDP Derivatives as Human DHFR Inhibitors. Bioinformation 2008, 3, 188. [Google Scholar] [CrossRef] [PubMed]
Raimondi, M.V.; Randazzo, O.; Franca, M. La; Barone, G.; Vignoni, E.; Rossi, D.; Collina, S. DHFR Inhibitors: Reading the Past for Discovering Novel Anticancer Agents. Molecules 2019, 24, 1140. [Google Scholar] [CrossRef] [PubMed]
Wróbel, A.; Drozdowska, D. Recent Design and Structure-Activity Relationship Studies on the Modifications of DHFR Inhibitors as Anticancer Agents. Curr Med Chem 2019, 28, 910–939. [Google Scholar] [CrossRef]
Yang, C.; Zhang, J.; Liao, M.; Yang, Y.; Wang, Y.; Yuan, Y.; Ouyang, L. Folate-Mediated One-Carbon Metabolism: A Targeting Strategy in Cancer Therapy. Drug Discov Today 2021, 26, 817–825. [Google Scholar] [CrossRef]
Newman, A.C.; Maddocks, O.D.K. One-Carbon Metabolism in Cancer. British Journal of Cancer 2017, 116, 1499–1504. [Google Scholar] [CrossRef]
Shuvalov, O.; Petukhov, A.; Daks, A.; Fedorova, O.; Vasileva, E.; Barlev, N.A. One-Carbon Metabolism and Nucleotide Biosynthesis as Attractive Targets for Anticancer Therapy. Oncotarget 2017, 8, 23977. [Google Scholar] [CrossRef]
Nurhanan, M.Y.; Nor Datiakma, M.A.; Muhammad Haffiz, J.; Sui Kiong, L.; Nor Hasnida, H.; Nor Jannah, S. The In Vitro Anti-Cancer Activities and Mechanisms of Action of 9-Methoxycanthin-6-One from Eurycoma longifolia in Selected Cancer Cell Lines. Molecules 2022, 27, 585. [Google Scholar]
Shailima Vardhini, R.D. Insilico Analysis of Protein-Ligand Docking of DHFR (Dihydro Folate Reductase) and Quassinoids. Int J Comput Appl 2013, 62, 14–19. [Google Scholar]
Hall, I.H.; Liou, Y.F.; Okano, M.; Lee, K.H. Antitumor Agents XLVI: In Vitro Effects of Esters of Brusatol, Bisbrusatol, and Related Compounds on Nucleic Acid and Protein Synthesis of P-388 Lymphocytic Leukemia Cells. J Pharm Sci 1982, 71, 345–348. [Google Scholar] [CrossRef] [PubMed]
Zhou, J.; Wang, T.; Dou, Y.; Huang, Y.; Qu, C.; Gao, J.; Huang, Z.; Xie, Y.; Huang, P.; Lin, Z.; et al. Brusatol Ameliorates 2, 4, 6-Trinitrobenzenesulfonic Acid-Induced Experimental Colitis in Rats: Involvement of NF-ΚB Pathway and NLRP3 Inflammasome. Int Immunopharmacol 2018, 64, 264–274. [Google Scholar] [CrossRef] [PubMed]
Olsen, N.J.; Spurlock, C.F.; Aune, T.M. Methotrexate Induces Production of IL-1 and IL-6 in the Monocytic Cell Line U937. Arthritis Res Ther 2014, 16, 1–8. [Google Scholar] [CrossRef] [PubMed]
Esim, O.; Oztuna, A.; Sarper, M.; Hascicek, C. Chitosan-Coated Bovine Serum Albumin Nanocarriers Mediate Efficient Delivery of Methotrexate in Breast Cancer Therapeutics. J Drug Deliv Sci Technol 2022, 77, 103906. [Google Scholar] [CrossRef]
Balkwill, F. TNF-α in Promotion and Progression of Cancer. Cancer and Metastasis Reviews 2006, 25, 409–416. [Google Scholar] [CrossRef] [PubMed]
Lahrita, L.; Hirosawa, R.; Kato, E.; Kawabata, J. Isolation and Lipolytic Activity of Eurycomanone and Its Epoxy Derivative from Eurycoma longifolia. Bioorg Med Chem 2017, 25, 4829–4834. [Google Scholar] [CrossRef] [PubMed]
Miyake, K.; Tezuka, Y.; Awale, S.; Li, F.; Kadota, S. Quassinoids from Eurycoma Longifolia. J Nat Prod 2009, 72, 2135–2140. [Google Scholar] [CrossRef]
Thigpen, J.T.; Bertelsen, K.; Eisenhauer, E.A.; Hacker, N.E.; Lund, B.; Sessa, C. Long-Term Follow-up of Patients with Advanced Ovarian Carcinoma Treated with Chemotherapy. Annals of Oncology 1993, 4, 35–40. [Google Scholar] [CrossRef]
Go, R.S.; Adjei, A.A. Review of the Comparative Pharmacology and Clinical Activity of Cisplatin and Carboplatin. Journal of Clinical Oncology 1999, 17, 409–422. [Google Scholar] [CrossRef]
Köhne-Wömpner, C.H.; Schmoll, H.J.; Harstrick, A.; Rustum, Y.M. Chemotherapeutic Strategies in Metastatic Colorectal Cancer: An Overview of Current Clinical Trials. Semin Oncol 1992, 19, 105–125. [Google Scholar]
Malaviya, A.N.; Sharma, A.; Agarwal, D.; Kapoor, S.; Garg, S.; Sawhney, S. Low-Dose and High-Dose Methotrexate Are Two Different Drugs in Practical Terms. Int J Rheum Dis 2010, 13, 288–293. [Google Scholar] [CrossRef] [PubMed]
Boik, J.C.; Newman, R.A. Structure-Activity Models of Oral Clearance, Cytotoxicity, and LD50: A Screen for Promising Anticancer Compounds. BMC Pharmacol 2008, 8, 1–13. [Google Scholar] [CrossRef] [PubMed]
Kardono, L.B.S.; Angerhofer, C.K.; Tsauri, S.; Padmawinata, K.; Pezzuto, J.M.; Douglas Kinghorn, A. Cytotoxic and Antimalarial Constituents of the Roots of Eurycoma longifolia. J Nat Prod 1991, 54, 1360–1367. [Google Scholar] [CrossRef] [PubMed]
Kuo, P.C.; Damu, A.G.; Lee, K.H.; Wu, T.S. Cytotoxic and Antimalarial Constituents from the Roots of Eurycoma Longifolia. Bioorg Med Chem 2004, 12, 537–544. [Google Scholar] [CrossRef] [PubMed]
Itokawa, H.; Kishi, E.; Morita, H.; Takeya, K. Cytotoxic Quassinoids and Tirucallane-Type Triterpenes from the Woods of Eurycoma Longifolia. Chem Pharm Bull (Tokyo) 1992, 40, 1053–1055. [Google Scholar] [CrossRef]
Miyake, K.; Li, F.; Tezuka, Y.; Awale, S.; Kadota, S. Cytotoxic Activity of Quassinoids from Eurycoma Longifolia. Nat Prod Commun 2010, 5, 1009–1012. [Google Scholar] [PubMed]
Elmore, S. Apoptosis: A Review of Programmed Cell Death. Toxicol Pathol 2007, 35, 495–516. [Google Scholar] [CrossRef]
Kollman, P.A.; Massova, I.; Reyes, C.; Kuhn, B.; Huo, S.; Chong, L.; Lee, M.; Lee, T.; Duan, Y.; Wang, W.; et al. Calculating Structures and Free Energies of Complex Molecules: Combining Molecular Mechanics and Continuum Models. Acc Chem Res 2000, 33, 889–897. [Google Scholar] [CrossRef]
Al-Thiabat, M.G.; Saqallah, F.G.; Gazzali, A.M.; Mohtar, N.; Yap, B.K.; Choong, Y.S.; Wahab, H.A. Heterocyclic Substitutions Greatly Improve Affinity and Stability of Folic Acid towards FRα. an In Silico Insight. Molecules 2021, 26, 1079. [Google Scholar] [CrossRef]
Al-Thiabat, M.G.; Gazzali, A.M.; Mohtar, N.; Murugaiyah, V.; Kamarulzaman, E.E.; Yap, B.K.; Rahman, N.A.; Othman, R.; Wahab, H.A. Conjugated β-Cyclodextrin Enhances the Affinity of Folic Acid towards FRα: Molecular Dynamics Study. Molecules 2021, 26, 5304. [Google Scholar] [CrossRef]
Larue, L.; Kenzhebayeva, B.; Al-Thiabat, M.G.; Jouan–Hureaux, V.; Mohd–Gazzali, A.; Wahab, H.A.; Boura, C.; Yeligbayeva, G.; Nakan, U.; Frochot, C.; et al. TLyp–1: A Peptide Suitable to Target NRP–1 Receptor. Bioorg Chem 2023, 130, 106200. [Google Scholar] [CrossRef] [PubMed]
Amir Rawa, M.S.; Al-Thiabat, M.G.; Nogawa, T.; Futamura, Y.; Okano, A.; Wahab, H.A. Naturally Occurring 8ß,13ß-Kaur-15-En-17-al and Anti-Malarial Activity from Podocarpus polystachyus Leaves. Pharmaceuticals 2022, 15, 902. [Google Scholar] [CrossRef] [PubMed]
U.S. Food & Drug Administration Pediatric Oncology Drug Approvals | FDA Available online:. Available online: https://www.fda.gov/about-fda/oncology-center-excellence/pediatric-oncology-drug-approvals (accessed on 17 March 2023).
U.S. National Library of Medicine Safety and Efficacy of Oral Methotrexate Tablets Combined with Immunotherapy During Radiotherapy for Unresectable/Metastatic Solid Tumors: A Single-Center, Prospective Study Available online:. Available online: https://clinicaltrials.gov/ct2/show/NCT05522582?term=methotrexate&cond=cancer&draw=4&rank=19 (accessed on 17 March 2023).
U.S. Food & Drug Administration Drug Approval Package: Alimta (Pemetrexed) NDA #021677 Available online:. Available online: https://www.accessdata.fda.gov/drugsatfda_docs/nda/2004/021677s000_alimtatoc.cfm (accessed on 16 March 2023).
U.S. National Library of Medicine Early Phase I Window of Opportunity Biomarker and Safety Trial to Test Pyrimethamine as an Inhibitor of NRF2 in HPV-Negative, Locally Advanced Head and Neck Squamous Cell Carcinoma Available online:. Available online: https://clinicaltrials.gov/ct2/show/NCT05678348?term=DHFR+inhibitor+clinical+trials&draw=2&rank=1 (accessed on 17 March 2023).
U.S. National Library of Medicine Melphalan With or Without Tumor Necrosis Factor in Treating Patients 202 With Locally Advanced Melanoma of the Arm or Leg. Available online: https://clinicaltrials.gov/ct2/show/NCT00003789 (accessed on 16 March 2023).
Hayslip, J.; Dressler, E. V.; Weiss, H.; Taylor, T.J.; Chambers, M.; Noel, T.; Miriyala, S.; Keeney, J.T.R.; Ren, X.; Sultana, R.; et al. Plasma TNF-α and Soluble TNF Receptor Levels after Doxorubicin with or without Co-Administration of Mesna-A Randomized, Cross-Over Clinical Study. PLoS One 2015, 10. [Google Scholar] [CrossRef] [PubMed]
U.S. National Library of Medicine Oncolytic Adenovirus TILT-123 With Pembrolizumab as Treatment for Ovarian Cancer. Available online: https://clinicaltrials.gov/ct2/show/NCT05271318?term=TNF-%CE%B1+clinical+trials+on+cancer&draw=3&rank=11 (accessed on 16 March 2023).
Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and Computational Approaches to Estimate Solubility and Permeability in Drug Discovery and Development Settings. Adv Drug Deliv Rev 1997, 23, 3–25. [Google Scholar] [CrossRef]
van de Waterbeemd, H.; Gifford, E. ADMET in Silico Modelling: Towards Prediction Paradise? Nat Rev Drug Discov 2003, 2, 192–204. [Google Scholar] [CrossRef] [PubMed]
Li Di Edward Kerns, B.H.; York, N.; Diego, S. Drug-Like Properties Concepts, Structure, Design, and Methods from ADME to Toxicity Optimization; 2016; ISBN 9780128010761.
Xiong, G.; Wu, Z.; Yi, J.; Fu, L.; Yang, Z.; Hsieh, C.; Yin, M.; Zeng, X.; Wu, C.; Lu, A.; et al. ADMETlab 2.0: An Integrated Online Platform for Accurate and Comprehensive Predictions of ADMET Properties. Nucleic Acids Res 2021, 49, W5–W14. [Google Scholar] [CrossRef]
Wang, N.N.; Dong, J.; Deng, Y.H.; Zhu, M.F.; Wen, M.; Yao, Z.J.; Lu, A.P.; Wang, J.B.; Cao, D.S. ADME Properties Evaluation in Drug Discovery: Prediction of Caco-2 Cell Permeability Using a Combination of NSGA-II and Boosting. J Chem Inf Model 2016, 56, 763–773. [Google Scholar] [CrossRef]
Low, B.S.; Ng, B.H.; Choy, W.P.; Yuen, K.H.; Chan, K.L. Bioavailability and Pharmacokinetic Studies of Eurycomanone from Eurycoma longifolia. Planta Med 2005, 71, 803–807. [Google Scholar] [CrossRef]
Furst, D.E.; Koehnke, R.; Burmeister, L.F.; Kohler, J.; Cargill, I. Increasing Methotrexate Effect with Increasing Dose in the Treatment of Resistant Rheumatoid Arthritis. J Rheumatol 1989, 16, 313–320. [Google Scholar]
Colom, H.; Farré, R.; Soy, D.; Peraire, C.; Cendros, J.M.; Pardo, N.; Torrent, M.; Domenech, J.; Mangues, M.A. Population Pharmacokinetics of High-Dose Methotrexate after Intravenous Administration in Pediatric Patients with Osteosarcoma. Ther Drug Monit 2009, 31, 76–85. [Google Scholar] [CrossRef]
Wang, N.N.; Huang, C.; Dong, J.; Yao, Z.J.; Zhu, M.F.; Deng, Z.K.; Lv, B.; Lu, A.P.; Chen, A.F.; Cao, D.S. Predicting Human Intestinal Absorption with Modified Random Forest Approach: A Comprehensive Evaluation of Molecular Representation, Unbalanced Data, and Applicability Domain Issues. RSC Adv 2017, 7, 19007–19018. [Google Scholar] [CrossRef]
Clark, D.E. In Silico Prediction of Blood–Brain Barrier Permeation. Drug Discov Today 2003, 8, 927–933. [Google Scholar] [CrossRef] [PubMed]
Monika; Sharma, A.; Aggarwal, V.; Sharma, M.; Dhingra, N. Lantadenes Targeting NF-ΚB in Cancer: Molecular Docking and ADMET Predictions Article In. Int J Pharma Bio Sci 2021.
Kazmi, S.R.; Jun, R.; Yu, M.S.; Jung, C.; Na, D. In Silico Approaches and Tools for the Prediction of Drug Metabolism and Fate: A Review. Comput Biol Med 2019, 106, 54–64. [Google Scholar] [CrossRef] [PubMed]
Ingelman-Sundberg, M. Human Drug Metabolising Cytochrome P450 Enzymes: Properties and Polymorphisms. Naunyn Schmiedebergs Arch Pharmacol 2004, 369, 89–104. [Google Scholar] [CrossRef] [PubMed]
Pan, Y.; Tiong, K.H.; Abd-Rashid, B.A.; Ismail, Z.; Ismail, R.; Mak, J.W.; Ong, C.E. Effect of Eurycomanone on Cytochrome P450 Isoforms CYP1A2, CYP2A6, CYP2C8, CYP2C9, CYP2C19, CYP2E1 and CYP3A4 in Vitro. J Nat Med 2014, 68, 402–406. [Google Scholar] [CrossRef] [PubMed]
Smith, D.A.; Beaumont, K.; Maurer, T.S.; Di, L. Relevance of Half-Life in Drug Design. J Med Chem 2018, 61, 4273–4282. [Google Scholar] [CrossRef] [PubMed]
Yang, H.; Sun, L.; Li, W.; Liu, G.; Tang, Y. In Silico Prediction of Chemical Toxicity for Drug Design Using Machine Learning Methods and Structural Alerts. Front Chem 2018, 6, 30. [Google Scholar] [CrossRef]
Skehan, P.; Storeng, R.; Scudiero, D.; Monks, A.; Mcmahon, J.; Vistica, D.; Warren, J.T.; Bokesch, H.; Kenney, S.; Boyd, M.R. New Colorimetric Cytotoxicity Assay for Anticancer-Drug Screening. JNCI: Journal of the National Cancer Institute 1990, 82, 1107–1112. [Google Scholar] [CrossRef]
Nurhanan, M.; Nor Azah, M.; Zunoliza, A.; Siti Humeriah, A.; Siti Syarifah, M.; Nor Hayati, A. In Vitro Anticancer Activity and High-Performance Liquid Chromatography Profiles of Aquilaria Subintegra Fruit and Seed Extracts. Journal of Tropical Forest Science 2017, 29, 208–214. [Google Scholar]
Berman, H.; Henrick, K.; Nakamura, H. Announcing the Worldwide Protein Data Bank. Nature Structural & Molecular Biology 2003, 10, 980. [Google Scholar]
Cody, V.; Gangjee, A. RCSB PDB - 5HPB: Human Dihydrofolate Reductase Complex with NADPH and 5-Methyl-6-(Phenylthio-4’trifluoromethyl)Thieno[2,3-d]Pyrimidine-2,4-Diamine Available online:. Available online: https://www.rcsb.org/structure/5hpb (accessed on 26 November 2022).
He, M.M.; Smith, A.S.; Oslob, J.D.; Flanagan, W.M.; Braisted, A.C.; Whitty, A.; Cancilla, M.T.; Wang, J.; Lugovskoy, A.A.; Yoburn, J.C.; et al. Medicine: Small-Molecule Inhibition of TNF-α. Science (1979) 2005, 310, 1022–1025. [Google Scholar]
Biovia, D.S. Discovery Studio Visualizer. Available online: https://discover.3ds.com/discovery-studio-visualizer-download (accessed on 26 November 2022).
Eswar, N.; Eramian, D.; Webb, B.; Shen, M.Y.; Sali, A. Protein Structure Modeling with MODELLER. Methods Mol Biol 2008, 426, 145–159. [Google Scholar] [PubMed]
Dolinsky, T.J.; Czodrowski, P.; Li, H.; Nielsen, J.E.; Jensen, J.H.; Klebe, G.; Baker, N.A. PDB2PQR: Expanding and Upgrading Automated Preparation of Biomolecular Structures for Molecular Simulations. Nucleic Acids Res 2007, 35, W522–W525. [Google Scholar] [CrossRef] [PubMed]
Olsson, M.H.M.; SØndergaard, C.R.; Rostkowski, M.; Jensen, J.H. PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical p K a Predictions. J Chem Theory Comput 2011, 7, 525–537. [Google Scholar] [CrossRef] [PubMed]
Williams, C.J.; Headd, J.J.; Moriarty, N.W.; Prisant, M.G.; Videau, L.L.; Deis, L.N.; Verma, V.; Keedy, D.A.; Hintze, B.J.; Chen, V.B.; et al. MolProbity: More and Better Reference Data for Improved All-Atom Structure Validation. Protein Science 2018, 27, 293–315. [Google Scholar] [CrossRef] [PubMed]
Kim, S.; Thiessen, P.A.; Bolton, E.E.; Chen, J.; Fu, G.; Gindulyte, A.; Han, L.; He, J.; He, S.; Shoemaker, B.A.; et al. PubChem Substance and Compound Databases. Nucleic Acids Res 2016, 44, D1202–D1213. [Google Scholar] [CrossRef]
Morris, G.M.; Ruth, H.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated Docking with Selective Receptor Flexibility. J Comput Chem 2009, 30, 2785–2791. [Google Scholar] [CrossRef]
Rowinska-Zyrek, M.; Salerno, M.; Kozlowski, H. Neurodegenerative Diseases – Understanding Their Molecular Bases and Progress in the Development of Potential Treatments. Coord Chem Rev 2015, 284, 298–312. [Google Scholar] [CrossRef]
Rowland, M.; Tozer, T.N.; Derendorf, H.; Hochhaus, G. Chapter 1. Therapeutic Relevance. Clinical Pharmacokinetics and Pharmacodynamics: Concepts and Applications; 4th ed.; Wolters Kluwer/Lippincott Williams & Wilkins: Philadelphia, 2011; ISBN 9780781750097. [Google Scholar]
Shargel, L.; Andrew, B.; Wu-Pong, S. Applied Biopharmaceutics & Pharmacokinetics ; Appleton & Lange Stamford, 1999; Vol. 264.
Forli, S.; Olson, A.J. A Force Field with Discrete Displaceable Waters and Desolvation Entropy for Hydrated Ligand Docking. J Med Chem 2012, 55, 623–638. [Google Scholar] [CrossRef]

Figure 1. 2D structures of a) eurycomanone and b) eurycomalactone.

Figure 2. Dose-response curves of (a) eurycomanone and (b) eurycomalactone when tested against ovarian (A2780), cervical (HeLa), colorectal cancer (HT29) cancer cell lines and, cardiomyocyte (H9C2) and liver (WRL-68) normal cell lines.

Figure 3. Visualization of (A) ovarian (A2780), (B) colorectal (HT29) and (C) cervical (HeLa) cancer cells when treated with eurycomanone at different concentrations (IC₅₀/ 5, IC₅₀, IC₅₀x 5) and incubation times (6, 24 and 48 h) with Hoechst 33342 staining. Arrows indicate different stages of apoptosis: (1) extensive chromatin condensation, (2) highly compact chromatin at the periphery of the nucleus and (3) apoptotic bodies. Magnification: x200.

Figure 4. Visualization of (A) ovarian (A2780), (B) colorectal (HT29) and (C) cervical (HeLa) cancer cells when treated with eurycomanone at different concentrations (IC₅₀/ 5, IC₅₀, IC₅₀x 5) and incubation times (6, 24 and 48 h) with Hoechst 33342 staining. Arrows indicate different stages of apoptosis: (1) extensive chromatin condensation, (2) highly compact chromatin at the periphery of the nucleus and (3) apoptotic bodies. Magnification: x200.

Figure 5. A 3D and 2D interaction analysis of docked models of eurycomanone (a) and eurycomalactone (b) with the TNF-α active binding site, and eurycomanone (c) and eurycomalactone (d) with the DHFR active binding site, was conducted using BIOVIA Discovery Studio Visualizer. The original crystal structure is depicted in pink for carbon (C), red for oxygen (O), cyan for fluorine (F), orange for sulphur (S), and blue for nitrogen (N). Eurycomanone is shown in green for C and red for O, and eurycomalactone is depicted in cyan for C and red for O.

Table 1. ¹H and ¹³C NMR spectral data of eurycomanone (1) and eurycomalactone (2) (400, 100MHz in CDCl₃).

Eurycomanone (1)			Eurycomalactone (2)
No	¹H δ ppm	¹³C δ ppm	No	¹H δ ppm	¹³C δ ppm
1	4.03 (1H, d, 8 Hz)	81.40	1	4.05 (1H, s)	81.27
2	-	197.89	2	-	197.41
3	6.14 (1H, d, 1.5 Hz)	162.98	3	6.13 (1H, s)	124.39
4	-	126.48	4	-	162.17
5	2.36 (1H, td, 2.4 Hz)	48.14	5	2.91 (1H, m)	49.38
6	2.08 (2H, m)	42.58	6	2.79 (2H, m)	36.21
7	-	199.32	7	-	205.56
8	-	53.02	8	-	51.16
9	2.02 (1H, t, 2.7, 13.3 Hz)	50.42	9	1.88 (1H, d, 3.2 Hz)	49.06
10	-	46.34	10	-	46.90
11	4.53 (1H, t, 6.8 Hz)	68.10	11	4.78 (1H, t, 4 Hz)	69.75
12	4.59 (1H, d, 8 Hz)	72.21	12	4.39 (1H, d, 4.4 Hz)	83.16
13	-	119.80	13	2.89 (1H, m)	32.33
14	3.25 (1H,d, 12.5 Hz)	76.28	14	3.02 (1H, m)	52.88
1516	5.25 (1H,t, 2.5, 18.7 Hz)-	79.81174.35	15	-	176.31
4-CH₃	1.81 (3H, s)	10.81	4-CH₃	1.64 (3H, s)	23.64
10-CH₃	2.00 (3H, s)	26.11	8-CH₃	1.96 (3H, s)	21.97
8-CH₂	2.07 (2H, m)	84.95	10-CH₃	1.27 (3H, s)	12.19
13’	7.6 (2H,s)	108.71	13-CH₃	1.18 (3H, d, 6.4 Hz)	32.33

Table 2. The IC₅₀ values (μM) ± SEM (n=9) of eurycomanone, eurycomalactone, cisplatin and methotrexate tested in ovarian (A2780), cervical (HeLa), colorectal cancer (HT29) cancer cell lines, and cardiomyocyte (H9C2) and liver (WRL-68) normal cell lines.

Compound	A2780	HeLa	HT-29	H9C2	WRL-68
Eurycomanone	1.37 ± 0.13	4.58 ± 0.090	1.22 ± 0.11	>50	1.34 ± 0.046
Eurycomalactone	2.46 ± 0.081	1.60 ± 0.12	2.21 ± 0.049	7.00 ± 0.43	2.71 ± 0.042
Cisplatin	1.77 ± 0.018	1.54 ± 0.12	1.38 ± 0.037	14.07 ± 1.14	1.13 ± 0.098
Methotrexate	0.016 ± 0.00050	0.094 ± 0.0043	0.059 ± 0.0010	>50	0.015 ± 0.00041

Note: SEM (Standard Error of the Mean) < 5%.

Table 3. Apoptotic Indices of A2780, HT-29 and HeLa when being treated with eurycomanone.

Concentrations/ Incubation Time	IC₅₀/5	IC₅₀	IC₅₀ x 5
Cell Line: A2780
6 h	4.30 ± 0.34^{a/ x}	8.48 ± 0.60^{ab/ x}	13.01 ± 0.29^{ac,bc/ x}
24 h	7.11 ± 1.60^{a/ xy}	13.13 ± 1.30^{ab/ y}	28.90 ± 0.93^{ac,bc/ xy}
48 h	14.72 ± 0.59^{a/ xz,yz}	31.74 ± 3.19^{ab/ xz,yz}	100.00 ± 0.00^{ac,bc/ xz,yz}

Cell Line: HT-29
6 h	4.08 ± 0.81^{a/ x}	5.24 ± 0.17^{b/ x}	5.97 ± 0.35^{ac/ x}
24 h	4.18 ± 0.21^{a/ y}	6.30 ± 1.01^{b/ y}	10.95 ± 0.71^{ac,bc/ y}
48 h	4.24 ± 0.20^{a/ z}	14.80 ± 0.56^{ab/ xz,yz}	26.20 ± 1.48^{ac,bc/ xz,yz}

Cell Line: HeLa
6 h	8.42 ± 0.20^{a/ x}	19.98 ± 0.76^{ab/ x}	41.37 ± 0.24^{ac,bc/ x}
24 h	13.76 ± 1.26^{a/ xy}	52.57 ± 1.40^{ab/ y}	94.56 ± 1.16^{ac,bc/ xy}
48 h	16.79 ± 0.92^{a/ xz,yz}	61.16 ± 0.63^{ab/ xz,yz}	100.00 ± 0.00^{ac,bc/ xz,yz}

Statistical analysis was performed using one-way analysis of variance (ANOVA) followed by Tukey’s multiple comparison test using Graphpad Prism version 7. ^a,b,c letters represent significant differences within the groups when compared to IC₅₀/5 (dose-dependent experiment) while ^x,y,z represent significant differences between the groups when compared to 6 h of treatment (time-dependent experiment).

Table 4. Percentage of apoptotic indices of A2780, HT-29 and HeLa when being treated with eurycomalactone.

Concentrations/ Incubation Time (h)	IC₅₀/5	IC₅₀	IC₅₀ x 5
Cell Line: A2780
6 h	3.84 ± 0.10^{a/ x}	4.26 ± 0.64^{b/ x}	15.81 ± 0.38^{ac,bc/ x}
24 h	4.12 ± 0.16^{a/ y}	6.29 ± 0.21^{ab/ y}	Cells died and completely detached^{ac,bc/ xy}
48 h	5.16 ± 0.065^{a/ xz,yz}	7.93 ± 2.28^{b/ xz}	Cells died and completely detached^{ac,bc/ xz}

Cell Line: HT-29
6 h	3.79 ± 0.45^{a/ x}	4.51 ± 0.23^{b/ x}	6.52 ± 1.22^{ac,bc/ x}
24 h	7.71 ± 0.51^{a/ xy}	12.41 ± 0.77^{ab/ xy}	20.21 ± 1.52^{ac,bc/ xy}
48 h	8.86 ± 0.68^{a/ xz}	14.28 ± 0.84^{ab/ xz,yz}	100.00 ± 0.00^{ac,bc/ xz,yz}

Cell Line: HeLa
6 h	5.45 ± 0.23^{a/ x}	17.67 ± 0.77^{ab/ x}	31.87 ± 2.19^{ac,bc/ x}
24 h	8.94 ± 0.21^{a/ xy}	32.00 ± 1.57^{ab/ xy}	62.20 ± 1.35^{ac,bc/ xy}
48 h	14.78 ± 0.12^{a/ xz,yz}	36.71 ± 1.19^{ab/ xz,yz}	100.00 ± 0.00^{ac,bc/ xz,yz}

Table 5. The free binding energies (ΔGbind) and inhibition constants (Ki) for eurycomanone, eurycomalactone, and the co-crystallized ligands of TNF-α (PDB ID: 2AZ5) and DHFR (PDB ID: 5HQY) were calculated into the active binding sites using AutoDock 4.2.

Compound	TNF-α		DHFR
Compound	*ΔG_bind (kcal/mol)	K_i (Micromolar uM)	*ΔG_bind (kcal/mol)	K_i (Micromolar uM)
Eurycomanone	-8.83	0.34	-8.05	1.25
Eurycomalactone	-7.51	3.11	-8.87	0.32
*124037103	*****	*****	-8.19	0.99
*5327044	-7.93	1.53	*****	*****

*ΔG_bind: Gibbs free binding energy (kcal/mol), *Ki: Inhibition constant (uM). *124037103: The PubChem ID of the co-crystalized ligands in DHFR (PDB ID: 5HQY). *5327044: The PubChem ID of the co-crystalized ligands in TNF-α (PDB ID: 2AZ5).

Table 6. Lipinski’s Rule of Five (RoF) predicted values for eurycomanone, eurycomalactone, and methotrexate.

Compounds	*M.W (g/mol)	*Hacc	*Hdon	*logP
Eurycomanone	408.14	9	4	0.215
Eurycomalactone	348.16	6	1	0.655
Methotrexate	454.17	13	7	-2.747

*M.W: Molecular weight (g/mol), *logP: Partition coefficient (Lipophilicity), *Hacc: Hydrogen bond acceptor, *Hdon: Hydrogen bond donor. *Lipinski Rule: M.W ≤ 500; logP ≤ 5; Hacc ≤ 10; Hdon ≤ 5, an orally active drug has no more than one violation of these criteria.

Table 7. Predicted ADMET properties of eurycomanone, eurycomalactone, and methotrexate using ADMETlab 2.0.

Property	Model Name	Predicted Value			Comment
Property	Model Name	Eurycomanone	Eurycomalactone	Methotrexate	Comment
Absorption	Papp (Caco-2 Permeability) cm/s	-5.54	-5.01	-6.73	* Papp ideal value is > −5.15 cm/s
Absorption	HIA (Human Intestinal Absorption) %	5.78	4.01	3.70	* HIA idea value is < 30%
Distribution	*PPB (Plasma Protein Binding) %	52.15	52.66	55.23	* PPB ideal value is < 90%
Distribution	Cross BBB (Blood Brain Barrier)	No	Yes	No
Metabolism	CYP1A2 substrate	No	No	No
	CYP2C19 substrate	No	No	No
	CYP2C9 substrate	No	No	No
	CYP2D6 substrate	No	No	No
	CYP1A2 inhibitor	No	No	No
	CYP2C19 inhibitor	No	No	No
	CYP2C9 inhibitor	No	No	No
	CYP3A4 inhibitor	No	No	No
Excretion	*CL (Clearance Rate) mL/min/kg	1.75	2.31	2.52	● High: CL >15 mL/min/kg ● Moderate: CL 5-15 mL/min/kg ● Low: CL <5 mL/min/kg
Excretion	T ½ (Half Lifetime) hr	0.03	0.18	0.89	● Long half-life: >3h ● Short half-life: <3h
Toxicity	H-HT (Human Hepatotoxicity)	+	+	+++	+ Low risk to be toxic. ++ Moderate risk to be toxic. +++ High risk to be toxic.
	AMES (Ames Mutagenicity)	++	+	+
	Carcinogenicity	+	+	++

*PBB: Plasma Protein Binding (PPB) Optimal: < 90%. Drugs with high protein-bound may have a low therapeutic index. *CL (Clearance Rate) mL/min/kg: high: CL >15 mL/min/kg, moderate: CL 5-15 mL/min/kg, and low: CL <5 mL/min/kg.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer