{'fraction_non_alphanumeric': 0.03528624037495815, 'fraction_numerical': 0.022363575493806496, 'mean_word_length': 4.946844515229942, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 3, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Cheminformatics Application in the Phytochemical and Biological Study of Eucalyptus globulus L. Bark as a Potential Hepatoprotective Drug

Natural products are considered as a good source of antifibrotic agents, but identifying and isolating bioactive molecule(s) is still challenging. Fortunately, numerous computational techniques have evolved to save time and efforts in this field. The aim of the current study was to utilize several cheminformatics software to study the chemical and biological features of the bark of Eucalyptus globulus cultivated in Egypt. Sirius software, with the aid of online databases, was used to process liquid chromatography–mass spectrometry (LC–MS) chemical profiling and predict precise molecular formulae, chemical classes, and structures. Accordingly, 37 compounds were tentatively identified, including 15 reported here for the first time from this species. Also, the BioTransformer tool was successfully applied for in silico virtual study of the human metabolism of these compounds, and 1960 different products were obtained through various metabolic pathways. Finally, an electronic library of the identified compounds and their metabolites were developed and docked in silico against eight different protein targets that are involved in the liver fibrosis process. The results revealed that the extract may have a potential hepatoprotective effect through several mechanisms and that the metabolites have the highest binding affinities to the relevant enzymes than their parent compounds. The extract was found to show potent cytotoxic activity against the liver cancer cell lines HEPG2 and HUH-7, and its absorption was enhanced through nanoformulation, as proved using the ex vivo everted gut sac method.


INTRODUCTION
Eucalyptus globulus Labill. is an evergreen tree that belongs to the family Myrtaceae and is cultivated worldwide. 1 In Egypt, several Eucalyptus species are grown as line plantings and used for shade, construction timbers, and poles or fuelwood. In the Eucalyptus industry, the bark is considered as one of the main by products. Eucalyptus bark is considered as a good source of phenolic compounds with diverse biological activities. 2,3 Polyphenols have several applications in cosmetics, food, and pharmaceutical industries. This class of compounds is known to demonstrate antioxidant, antimicrobial, antidiabetic, anti-inflammatory, antihyperlipidemic, hepatoprotective, nephroprotective, cardioprotective and, anticancer effects. 4 Ingestion of polyphenol-rich extracts such as green tea drink was found to be beneficial in reducing the risk of cancer incidence and mortality. 5 However, its poor bioavailability, fast metabolism in the gastrointestinal tract (GIT) and the liver, and rapid elimination from the body hinder its utilization in the pharmaceutical industry. 6 Generally, the bioavailability of polyphenols could be enhanced by increasing their solubility and dissolution rate in the gastrointestinal fluid. Several techniques were used for improving solubility such as solid-lipid nanoparticles, nanostructured lipid carriers, nanoemulsion, liposomes, and self-micro emulsifying drug-delivery systems (SMEDDS). 7 Recently, several cheminformatics software have been developed to reduce time and cost in developing bioactive molecules. Several in silico tools for the prediction of absorption, distribution, metabolism, and toxicity were successively employed. 8 Furthermore, molecular docking studies were applied to predict the efficacy in terms of its binding to human target proteins and reduce efforts in in vitro and preclinical in vivo studies.
In the present study, the ethanolic extract of E. globulus cultivated in Egypt was chemically characterized by the spectrophotometric estimation of the total phenolics and flavonoid contents, profiling using ultrahigh-performance liquid chromatography−electrospray ionization−quadrupole time-of-flight tandem mass spectrometry (UPLC-ESI-QTOF-MS−MS) as well as its antioxidant activity. Consequently, the hepatoprotective effect of the extract was evaluated in light of recent cheminformatics software. In addition, the cytotoxic activities of the extract against two liver cancer cell lines (HEPG2 and HUH-7) were determined. SMEDDS was applied to enhance the intestinal absorption of the extract, which was assessed ex vivo using the everted gut sac method.
2.4. Chemical Characterization of the Extract. 2.4.1. Determinations of Total Phenol and Total Flavonoid Contents. The Folin−Ciocalteu method using gallic acid as a standard was adopted to determine the total phenol content (TPC) of EE, as described before. 9 The measurement was performed at 630 nm (for TPC) and 420 nm (for TFC). Total phenols were expressed as gallic acid equivalents (GAE) mg/g extract (dw) using a calibration curve of a freshly prepared gallic acid solution. The total flavonoid content (TFC) was quantified using the aluminum chloride colorimetric method, as reported by Herald and Gadgil. 10 The TFC was expressed as rutin equivalents (RE) mg /g extract (dw).
2.4.3. Tentative Identification of the Secondary Metabolites. The tentative identification of the metabolites of EE was carried out using Sirius software (ver. 4.7.4), which is available through https://bio.informatik.uni-jena.de/software/sirius/ as a free open-source software to predict the fragmentation and the molecular formulae. 11 The chemical structures were predicted by CSI: FingerID, 12 while CANOPUS was used for the prediction of natural product classes directly from MS/ MS. 13 2.4.4. Antioxidant Activities. 2.4.4.1. DPPH Assay. The DPPH assay method is based on the use and reduction of the free radical of DPPH. 14 The absorbance was measured at 517 nm in a microplate reader. The control and reference compound ascorbic acid were also measured. The experiment was carried out in triplicate, and the radical scavenging activity (Inhibition %) was calculated according to the following equation: The result was presented as IC 50 (μM).

ABTS Assay.
The antioxidant capacity of sample to scavenge ABTS (2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid)) free radicals was evaluated according to the methodology proposed by Faitanin and Gomes. 15 The absorbance was measured at 690 nm using a microplate reader (Tecan, USA). Data are represented as mean ± SD according to the following equation: The result was presented as μM ascorbic acid equivalent per mg extract (μM AAE/mg extract) according to the following calibration curve equation:  17 At the end of incubation time, the decrease in the produced color intensity was measured at 540 nm using the microplate reader (Tecan, USA). Data are represented as mean ± SD according to the following equation: Results of the extract tested were presented as μM ethylenediaminetetraacetic acid (EDTA) equivalent/mg extract with the calibration curve equation: 2.5. Virtual Human Metabolism Prediction. The human metabolism of the tentatively identified compounds was virtually predicted using the Biotransformer tool. 8 The metabolism prediction was carried out using the Superbio metabolic transformation option to cover biotransformation occurring both in human tissues as well as the gut microbiota. The Biotransformer tool is freely available on https:// biotransformer.ca/ 2.6. Molecular Docking. All the docking studies were carried out using protein−ligand complexes with crystal structures. Before the docking process, the protein structures were processed by AutoDockTools GUI of AutoDock 4.2 software. 19 All Autodock software are freely available on https://autodock.scripps.edu/resources/. The processing included adding hydrogen atoms, removing water molecules and the co-crystallized ligands from the protein−ligand complexes, and the regulation of the amide groups in the side chains of asparagine and glutamine to modify their connections with nearby residues and atoms. After that, the docking process was continued using AutoDock Vina. 20 The target human receptors selected in this study were the ligand-binding domains (LBDs) of ZAP-70, Kaep1, Angio-II-Type-1, Jak-2, FGFR1 kinase-1, EGFR, TGFβR1, and VEGFR-2. All the LBDs were retrieved from the RCSB Protein Data Bank.
2.7. Preparation of SMEDDS. The SMEDDS were prepared according to the method that was modified by Nematallah et al. 21 Capryol 90 was the oily phase, while Tween 80 was the surfactant, and Transcutol P was the cosurfactant with a ratio 12.5:53.5:34, respectively. EE was dissolved in the oil/surfactant mixture with the aid of sonication at 25°C for 5 min. Then, Transcutol P was added while swirling and stirring. The prepared SMEDDS was stored in tightly sealed glass vials at a temperature range of 2−8°C until further use.
2.8. Characterization of the Microemulsion. SMEDDS were reconstituted with 4 times the volume of distilled water to form a microemulsion for characterization. The mean diameter of the droplet (droplet size) and the size distribution (polydispersity index; PDI) of the formed microemulsion were measured using the dynamic light scattering technique at a scattering angle of 173°on Zetasizer Nano ZS (Malvern Instruments, UK). The shape and size of the droplets were further characterized using Jeol transmission electron microscopy (TEM) (JEM-1230, Japan), a microemulsion drop was placed on a copper-coated grid, and a piece of filter paper was used to remove excess microemulsion. The magnification was 100,000× with an operating voltage of 80 kV.
2.9. Ex vivo Analysis of Microemulsion Absorption. The everted gut sac method was used to test the enhancement of the absorption of EE in intestines with the aid of microemulsion. The method was slightly modified from the method previously developed by Schilling and Mitra 22 and modified by Tambe and Mokashi. 23 The buffer used was a modified Kreb's Ringer phosphate bicarbonate (KRPB) solution, pH 7.4. Male Sprague−Dawley rats weighing 200− 300 g were fasted for 16−20 h and anesthetized by an i.p. injection of 60 mg/kg pentobarbital. After making a midline incision in the abdomen, the small intestine was cut at two positions at 3 cm distal to the duodenum and at the ileocecal junction, carefully removed, and placed immediately in ice-cold saline solution. The intestines were segmented into 8−10 cm length segments, followed by eversion on a thin stainless-steel rod, rinsing with ice-cold saline solution, and then secured to the tip of a 1 mL disposable syringe barrel. The other end of the segment was ligated with silk thread and attached to the bottom of the everted gut sac set up. The everted gut sac was filled with a known volume (1−2 mL) of the modified KRPB buffer solution.
Each everted gut sac set up was filled with 25 mL of KRPB buffer. In one set up, a known amount of EE was suspended in the buffer, and the same amount of EE dissolved in the SMEDDS was added to the second setup. Each setup was left for 30 min at a temperature of 37°C with simple agitation; then, the buffer inside the everted gut sac was withdrawn and tested for their content of extract spectrophotometrically.
2.10. Cytotoxicity Assay. Human hepatocellular (HCC) carcinoma cell lines (HEPG2 and HUH-7) were obtained from VACSERA (Giza, Egypt). The cells were maintained in Dulbecco's modified Eagle medium (DMEM) [+L-glutamine] supplemented with 10% fetal bovine serum, 100 units/mL of penicillin, 100 μg/mL of streptomycin, and 0.25 μg/mL amphotericin B (Gibco) at 37°C in a humidified incubator containing 5% CO 2 . All cell culture procedures were performed in a class II hood and incubations were done inside the incubator at 37°C. Initial cytotoxicity screening of EE was investigated using the MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) colorimetric cell viability assay that assess the cells' metabolic activity. Cells were seeded in 96-well plates at a density of 15,000 cells/well and allowed to attach overnight. Next day, the cells were treated with a serial dilution of EE at concentrations of 500, 250, 125, 62.5, 31.25, 15.625, 7.8, and 3.9 μg/mL for 72 h. Then, after 72 h, the culture medium was replaced with 100 μL/well of 5 mg/mL MTT powder solution dissolved in complete DMEM medium. Cells were then incubated for 2 h at 37°C inside the incubator. Finally, the media was discarded, and the formed formazan crystals were dissolved in 100 μL/well DMSO and further incubated for 10 min at 37°C inside the incubator. The optical density was then measured spectrophotometrically at 570 nm, and % cell viability was calculated relative to untreated controls. The inhibitory concentration 50 (IC 50 ) of the extracts on the different cell lines was determined via GraphPad Prism software (Ver 7.00, 2016, USA) using the nonlinear regression analysis.

RESULTS AND DISCUSSION
3.1. Extraction and Determination of Total Contents of Phenolics and Flavonoids. In order to prepare an extract of E. globulus bark rich in phenolic compounds, the powdered bark was extracted with 70% ethanol and the extract was loaded onto a column of diaion-HP20. Elution started with water to remove salts, carbohydrates, or water-soluble impurities. Then, elution with ethanol yielded an extract rich in phenolics (EE). The total phenolics (TPC) and total flavonoid (TFC) contents of the EE were found to be 436.70 ± 34.63 μg GAE/mg extract and 7.32 ± 0.55 μg RE/mg extract, respectively. TPC of EE was found to be similar to that reported for the aqueous methanolic extract of E. globulus bark grown in Portugal, 24 and more than that reported in different types of extracts of E. globulus bark grown in Spain and Chile. 3, 25 3.2. UPLC-ESI-QTOF-MS−MS Analysis of the Extract. For the first time, the extract of the bark of E. globulus growing in Egypt was analyzed using the high-resolution LC−MS−MS technique in the negative-ion mode ( Figure 1). Tentative identification of the metabolites was carried out on Sirius open-source software utilizing CSI: Finger ID and CANOPUS integrated tools. 11 Thirty-seven compounds were tentatively identified, including 12 phenolic acids derivatives, 4 flavonoids, 11 triterpene derivatives, and 5 fatty acids (Table 1)  relation between antioxidant activity and TPC that was previously reported for extracts of other plants. 42 The promising results of the antioxidant activity and richness of the extract with polyphenolic compounds encouraged us to utilize this extract as a hepatoprotective drug as it is well known that the main mechanism of the action of polyphenols as hepatoprotective drugs is mediated by decreasing the oxidation stress level in the liver. 43,44 3.4. Virtual Human Metabolism Prediction. One of the main objectives of the current study is to prepare a nanoformulation that has the ability to enhance the bioavailability of the E. globulus extract as a hepatoprotective/antifibrotic drug. There are two approaches of the nanoformulations: the first one is the nanoformulation that protects the active principles from the human metabolic degradation and deliver the active drug to its site of action, while the second approach is the nanoformulations that enhance the solubility and absorption of the drug, leaving it susceptible to metabolic activation after absorption. The decision is mainly based on the activity of the parent compounds versus their metabolic products. First, all the tentatively identified compounds from the LC− MS analysis were searched in the literature for their reported stereoisomers that result in the preparation of an electronic library of 47 different parent compounds. Subsequently, a virtual in silico metabolism was carried out for all those parent compounds using the Biotransformer tool. 8 The metabolism covered all the pathways including gut microbiota, EC-based, phase I and phase II metabolism. The process yielded an electronic library of 3093 metabolic products, and the number decreased to 1960 different metabolic products after the removal of the similar compounds (see the Supporting Information for full electronic library (.mol file) of the compounds and their metabolic products).
To validate the output of the virtual metabolism process, some of the metabolic products were compared with the reported metabolism processes in the literature. For example, from the 39 reported metabolic products of epicatechin, 37 compounds were successfully predicted using the Biotransformer tool. 8 Similarly, all the reported urolithin derivatives and their glucuronide glycosides in a study that was carried out by Tomas−Barberan, Garcia−Villalba 45 were found to be predicted successfully.
3.5. Molecular Docking. Molecular docking studies were carried out to evaluate the liver antifibrotic ability of E. globulus bark extract. The studies included both the parent-identified compounds as well as their metabolic products. The multitargeted approach focusing on different pathways is the most promising therapeutic strategy against fibrotic diseases. 46 Totally, 2008 molecules were docked against eight human protein targets that are involved in the liver fibrosis process. The results demonstrated a promising binding ability (Table  3). Collectively, the metabolic products showed binding affinities higher than their parent compounds as one of the parent compounds showed the highest binding affinity to a single target, while the metabolites showed the highest affinity to the other seven targets.
3.5.1. Janus Kinase-2 Enzyme (JAK-2). JAK-2 in a nonreceptor tyrosine kinase enzyme, and it has the ability to activate the signal transducer and activator of transcription (STAT). In many studies, it is proved that JAK2/STAT3 signaling promotes fibrosis, angiogenesis, and inflammation in many diseases. 47 Therefore, targeting and inhibition of the JAK-2 enzyme may have a promising therapeutic for liver fibrosis. 48 From the docking studies, it was found that the compound with the highest binding affinity with JAK-2 (PDB ID 3KCK) was a urolithin glucuronide derivative with a score of −9.8, forming five hydrogen bonds with GLU-898, MET-929, LEU-932, and ASP-994 ( Figure 2A). This binding is nearly better than the binding of the cocrystallized ligand 3chloro-4-(4H-3,4,7-triazadibenzo[cd,f ]azulen-6-yl) (binding score −9.7) with only three hydrogen bond interactions with the active site on GLU-898, LEU-932, and PHE-995. The parent compound with the highest binding affinity was found to be the flavonoid astilbin (16) with score −9.2. This urolithin glucuronide derivative is a metabolic product that could be produced by several metabolic pathways from several phenolic products identified in the E. globulus extract such as hexahydroxy diphenoyl hexoside (1), bis [hexahydroxy diphenoyl] hexoside (3), galloyl bis [hexahydroxy diphenoyl] hexoside (7), ellagic acid glucuronide (9), ellagic acid hexoside (10), and ellagic acid (15), which comprise around 20% of the total UPLC chromatogram area, and the predicted biotransformation pathway is illustrated in the Supporting Information ( Figure S1).
3.5.2. Angiotensin II Receptor Type 1 (ACE). The renin− angiotensin system is frequently activated in patients with chronic liver diseases. Angiotensin II has been suggested to play an important role in liver fibrogenesis. It induces hepatic stellate cell proliferation and upregulates the transforming growth factor beta-1 expression. 48 Both a parent compound and a metabolic product showed much higher binding affinities when compared to the standard inhibitor drug olmesartan (score −8.6) and the control drug silybin B (score −7.9). The highest binding affinity with ACE (PDB ID 4ZUD) was for another urolithin glucuronide derivative (score −9.4), which is another metabolic product from the same mentioned parent compounds in JAK-2 ( Figure S2). It forms seven hydrogen bonds with four amino-acid residues TYR-35, TYR-92, SER-105, and ARG-167, and it also has Pi-Pi interactions with residues TRP-84 and TYR-87 ( Figure 2B).
3.5.3. Zeta-Chain-Associated Protein Kinase 70 (ZAP Kinase). One of the major signaling pathways that control the liver homeostasis, pathophysiology and regulation of the inflammation−fibrosis−cancer axis is nuclear factor-kappa B (NF-κB). 49 The ZAP kinase enzyme mediates the activation of NF-κB pathway through tyrosine phosphorylation. Thus, the inhibition of ZAP kinase could have a promising effect in the cessation of the liver fibrosis process by the indirect deactivation of NF-κB. 49 From the tested compounds, the compound with the highest binding affinity with ZAP kinase (PDB ID 1U59) is the metabolite ellagic acid glucuronide with score −10.0. Although the binding score of this metabolite is less than that of the cocrystallized ligand staurosporine (score −11.3), it forms a strong interaction with the active site through four hydrogen bonds with residues CYS-346, ALA-417, GLY-418, and ASP-479 ( Figure 2C). This metabolite could be produced from the parent compounds ellagic acid (15), ellagic acid hexoside (10), and ellagic acid glucuronide (9) ( Figure S3). Although there is an ellagic acid glucuronide (9) in the parent compounds, its binding score (−9.2) was less than the metabolite, this could be due to different orientations of the glucuronide moiety in the active site. 3.5.4. Transforming Growth Factor Beta Receptor I (TGFβR-1). In most of the chronic liver diseases, the level of TGF-β increases markedly, activating the transformation of hepatic stellate cells to myofibroblasts and increased hepatocyte cell death, which causes liver fibrosis. 50 Therefore, the inhibition of TGF-β or blocking its downstream signaling pathway is one of the approaches to prevent liver fibrosis. The highest binding affinity with TGFβR-1 (PDB ID 2X7O) from the tested compounds was found to be astilbin with a score of −10.4, which is the only parent compound that shows higher binding affinity than the metabolites against all the eight target proteins in the current study. There are also two flavonoid glycosides from the metabolites that showed a relatively promising binding affinity with the active site of TGFβR-1 with a score of −10.1. These results are promising when compared with the binding score of the cocrystallized ligand (3Z)-Nethyl-N-methyl-2-oxo-3-[phenyl-[[4-(piperidin-1-ylmethyl)phenyl]amino]methylidene]-1H-indole-6-carboxamide (−10.5). Astilbin forms three hydrogen bonds with residues LYS-232 and LYS-337. It also interacts by 4 pi-alkyl and 2 pisigma interactions with several amino-acid residues in the active site of TGFβR-1 ( Figure 2D).
3.5.5. Kelch-Like ECH-Associated Protein 1 (Kaep1). Nrf2 protein is a key factor in the expression of the antioxidant defense system; phase II detoxifies enzymes by translocation into the nucleus and binding with the antioxidant responsive element. 51 Normally, Nrf2 is bound to Kaep1 protein, and it is activated by dissociation from Kaep1 by either Nrf2 phosphorylation or Kaep1 modification. 52 Finding a ligand with high binding affinity to Kaep1 will cause its modification and activation of Nrf2. Thus, targeting Kaep1 is considered beneficial, especially against liver diseases. One of the phenolic metabolites, which is a urolithin derivative, showed a very strong binding with Kaep1 (PDB ID 3VNG) with a score of −8.4 when compared with the binding score of the cocrystallized ligand 2-(3-((3-(5-(furan-2-yl)-1,3,4-oxadiazol-2yl)ureido)methyl)phenoxy)acetic acid (−7.5). This metabolite forms five hydrogen bonds with TYR-13, SER-42, ARG-59, and ASN-61 residues ( Figure 2E). This metabolite is a result of metabolism of several parent phenolic compounds identified in E. globulus extract such as compounds 1, 3, 7, 9, 10, and 15 ( Figure S4).
Collectively, it could be concluded from the docking studies that the metabolism products of the compounds identified in E. globulus extract are biologically active and have higher binding affinities to the protein targets under investigation in the current study. Furthermore, by the comparison of the binding affinities of these compounds with those of the cocrystallized ligand of each protein and silybin B as a reference hepatoprotective drug, it could be deduced that E. globulus bark extract has a promising liver antifibrotic activity and it could be included in further biological studies to prove its effect.
3.6. Preparation and Characterization of the Microemulsion. Based on the previous results of virtual metabolism and molecular docking, SMEDDS was chosen as an approach to enhance the bioavailability of EE. The EE was dissolved in a microemulsion system composed of Capryol 90, Tween 80, Transcutol P, and water. The size of the droplets loaded with E. globulus extract was measured using a zeta sizer and found to be 14.90 ± 1.2 nm with PDI 0.255 ± 0.032 ( Figure 3A). Loading EE in the microemulsion has nearly no effect on the droplet size as the size of the void droplets was 18.28 ± 1.4 nm with PDI 0.331 ± 0.041 ( Figure 3B). Additionally, the droplets were characterized by TEM, and it was found that the droplets have a spherical shape with uniform distribution with sizes ranging from 5 to 7 nm ( Figure 3C).
In previous studies, nanoemulsions of E. globulus essential oil were prepared using grape seed oil and a hydrophilic surfactant to develop an efficient transdermal nanocarrier for the alternative analgesic therapy. 56 Also, urea formaldehyde nanocapsules loaded by a hydroalcoholic extract of E. globulus leaf extract were developed as a natural nanopesticide for pest control. 57 To the best of our knowledge, the application of SMEDDS to enhance the bioavailability of E. globulus bark extract was reported here for the first time.
3.7. Ex Vivo Analysis of the Microemulsion Absorption. To study the effect of the microemulsion on the intestinal absorption of the EE, an everted gut sac technique was carried out. To determine the quantity of the EE, a collective UV−vis spectrophotometric calibration curve of EE dissolved in the microemulsion was established. The samples were measured spectrophotometrically against a blank composed of void microemulsion at wavelength (λ max ) 278 nm. After 30 min of the experiment, only 9.5% of the extract suspended in water was transferred into the intestinal sac, while 24% of the extract dissolved in the microemulsion was absorbed successfully in the intestinal sac. This result indicates that SMEDDS enhanced the absorption of the E. globulus bark extract by nearly 2.5-fold.
The everted gut sac method is a known ex vivo technique to study the effect of several parameters on the absorption of drugs and its different formulations. 58 It was also used to study the absorption mechanisms of some herbal extracts such as Schisandra chinensis and Inula cappa. 59, 60 3.8. Cytotoxic Activity on Liver Carcinoma Cell Lines. MTT cytotoxicity assay results revealed that the EE possesses strong cytotoxic potential on two liver cancer cell lines ( Figure  4). On the hepatic carcinoma cell lines, the extract was more potent on the HEPG2 cell line with IC 50 of 53.14 μg/mL compared to 153 μg/mL on HUH-7 cells.
Our results revealed a strong anticancer potential of the E. globulus bark extract on both tested cell lines with dosedependent reduction in cell viability and calculated IC 50 of less than 300 μg/mL. The most potent effect was observed on HEPG2 cells, followed by HUH-7.
Even though the anticancer potential of E. globulus essential oils has been reported before, the potency of its bark extract on different cancer cell lines is still under investigation. E. globulus essential oil demonstrated the cytotoxic effect on colon cancer cell line SW48 and liver cancer cell line HepG2, normal human embryonic kidney cells HEK293t, and was toxic to the normal skin fibroblasts tested. 61 Another study by Teixeira et al. investigated the anticancer potential of the decoction extract of aerial leafy parts of E. globulus on colorectal (HCT-15), pancreatic (PANC-1), and nonsmall cell lung cancer (NCI-H460) cell lines with the highest potency observed on NCI-H460 cell lines. The cell cytotoxicity effect was attributed to the G0/G1 cell cycle arrest and an increase in the expression of tumor suppression proteins (e.g. p53, p21, and cyclin D1 proteins). 62 On the other hand, the bark extract of E. globulus was reported to cause a dose-and time-dependent reduction in the cell viability of MDA-MB-231 cell lines. 63 In line with the published reports, our findings further confirmed that E. globulus bark extract could be further developed for anticancer therapies.

CONCLUSIONS
The present findings clearly demonstrate that the integration of in silico-predicted biological activity profiles of the parent compounds and their metabolites was crucial in understanding the antifibrotic activity of the bark extract of E. globulus cultivated in Egypt. The bark extract was found to be a promising antifibrotic agent, and the predicted metabolites were more active than their parent compounds. Also, SMEDDS was effective in developing a more bioavailable formula with better metabolism and absorption. However, to the best of our knowledge, an investigation that combines the in silico estimates of drug metabolism and predictions of the biological activity profiles for both the parent drug and its metabolites has never been conducted. The results clearly demonstrate that the integration of in silico-predicted biological activity profiles of the parent drug substance and its metabolites improves the computer-aided assessments.
A full electronic library of the tentatively identified compounds of the ethanolic extract of E. globulus bark and their metabolic products after the removal of the repeated molecules (MOL) Virtual metabolic pathways of the compounds with the highest binding affinities to the protein targets under investigation (Figures S1−S7) (PDF)