Anticancer potential of Phoenix dactylifera L. seed extract in human cancer cells and pro-apoptotic effects mediated through caspase-3 dependent pathway in human breast cancer MDA-MB-231 cells: an in vitro and in silico investigation

Background Phoenix dactylifera L. has a diverse set of pharmacological properties due to its distinct phytochemical profile. The purpose of this study was to investigate the anticancer potential of Phoenix dactylifera seed extract (PDSE) in human breast cancer MDA-MB-231 and MCF-7 cells, as well as liver cancer HepG2 cells, and to investigate the anticancer efficacy in triple-negative MDA-MB-231 cells, followed by in silico validation of the molecular interaction between active components of PDSE and caspase-3, an apoptosis executioner protein . Methods In this study, human cancer cell lines were cultured and subsequently treated with 10 to 100 μg/mL of PDSE. MTT test was performed to determine the cell viability, MMP was measured using fluorescent probe JC-1, nuclear condensation was determined by Hoechst 33258 dye, Annexin V-FITC & PI staining and cell cycle analysis were evaluated through flow cytometer, and apoptotic markers were detected using western blotting. The bioactive agents in PDSE were identified using high-performance liquid chromatography (HPLC) analysis. The binding affinity was validated using molecular docking tools AutoDock Vina and iGEMDOCK v2.1. Results Cell viability data indicated that PDSE inhibited cell proliferation in both breast cancer cells and liver cancer cells. MDA-MB-231 cells showed maximum growth inhibition with an IC50 value of 85.86 μg/mL for PDSE. However, PDSE did not show any significant toxicity against the normal Vero cell line. PDSE induced MMP loss and formation of apoptotic bodies, enhanced late apoptosis at high doses and arrested cells in the S phase of cell cycle. PDSE activated the enzymatic activity of cleaved caspase-3 and caused the cleavage of poly-ADB ribose polymerase (PARP) protein. PDSE upregulated pro-apoptotic Bax protein markedly but no significant effect on tumor suppressor protein p53, while it downregulated the anti-apoptotic Bcl-2 protein expression. HPLC analysis showed the presence of rutin and quercetin bioactive flavonols in ethanolic extract of PDS. Interestingly, both active components revealed a strong binding interaction with amino acid residues of caspase-3 (PDB ID: 2XYP; Hetero 4-mer - A2B2) protein. Conclusion PDS could serve as a potential medicinal source for apoptotic cell death in human breast cancer cells and, thus, could be used as a promising and crucial candidate in anticancer drug development. This study warrants further in vivo research, followed by clinical investigation. Supplementary Information The online version contains supplementary material available at 10.1186/s12906-022-03533-0.


Background
Cancer is an ensemble of diseases that develop over a long stretch of time, causing millions of mortalities across the globe [1]. Owing to its extremely aggressive nature and poor prognosis and survival rate, cancer remains an important public health issue worldwide. Amongst various cancers, breast and liver cancer are of great concern globally, accounting for 2.26 and 2.21 million new cases, respectively in 2020 [2]. Recent cancer statistics have suggested a significant rise in the number of breast cancer patients, indicating breast cancer (11.7%), as a most prevalent type of cancer diagnosed every year, followed by lung (11.4%), colorectal (10.0%), prostate (7.3%), and stomach (5.6%) cancers [3]. As per the statistic, 1 out of 4 cancer patients are diagnosed with breast cancer, causing 1 out of 6 cancer related mortalities [4]. Among Indian females, breast cancer ranks as number one cancer, having a mortality rate of 12.7 per 100,000 women [5]. The incidence of liver and breast cancers has increased in many parts of the world, notably in India, North and South America, as well as in most European countries.
Natural products have emerged as a benchmark in the process of discovery and development of novel drugs, particularly for anticancer and anti-infectious agents [6]. It is noteworthy that about fifty percent of anti-cancer drugs and therapies are either derived from natural sources or natural products [7]. A range of herbal derived components such as alkaloids, polyphenols, flavonoids, terpenoids and polysaccharides are being used against various cancers and communicable diseases [8,9]. Various therapeutic strategies and interventions are available for these diseases; however, most of them remain incurable because of drug resistance. In addition, some natural products are not free from side effects. Therefore, the search and discovery of novel plant-derived products are essential to improve the impact and outcome of plantbased therapeutic agents.
Ajwa dates are the fruits of the native tree Phoenix dactylifera L. of Saudi Arabia. For thousands of years in Egypt and the Middle East, the tree has been enlisted for its diverse medicinal uses in Pharmacopeia [10]. All parts of the plant, particularly fruits, seeds, leaves, flowers and roots, are utilized in Ayurveda-based formulations and have a significant role as a Rasayana i.e. rejuvenation medicine [11]. The therapeutic effects of date palm have been well-documented [12]. Ajwa dates are cultivated only in Saudi Arabia/Al-Madinah Al-Munawara and be effective against several types of diseases [13]. Ajwa date seeds have been utilized in Egyptian traditional medicine to treat a variety of infectious disorders diabetes and cancer [14,15]. Traditionally and historically, Ajwa dates are used mainly for their antioxidant, anti-inflammatory and hepatoprotective, anticancer and cardiac function improvement effect [16,17].
Primary breast cancer is confined to the milk glands or ducts of the breast. However, in some cases, primary breast cancer becomes aggressive and spreads to distant organs such as the liver causing metastasis [18]. Therefore, in the present study, cancer cell lines from two different tissues viz. breast (MDA-MB-231 and MCF-7), and liver (HepG2) have been used to investigate the anticancer potential of Phoenix dactylifera seed extract (PDSE). African green monkey normal kidney epithelial cell line (Vero) has been widely used in various toxicological, virology and pharmacological research domains, particularly as a sensitive model for toxicity assays of compounds of various nature, either chemical or microbial toxins [19][20][21][22]. Thus, to test the toxicity of PDSE, Vero cell line was used as a control to mimic normal mammalian cells. HPLC characterization was carried out to detect the presence of active phytoconstituents and rutin and quercetin were found in the extract. In silico computational analyses revealed that these active ingredients of PDSE displayed excellent binding interaction with caspase-3 protein. Thus, the present study indicates that PDSE might be utilized as an adjunct to the mainline of cancer treatment or be developed as a therapeutic anti-cancer agent in breast cancer therapy.

HPLC analysis of PDSE
An Waters 515 HPLC Pump system (Milford, USA) equipped with a W2998 PDA detector, a pump control module, a Waters column temperature controller, and an empower chromatography workstation was utilized to characterize PDSE, as reported previously [23]. An X BRIDGE C18 5 m, 4.6 × 250 mm reverse-phase column with gradient elution as the mobile phase was used for chromatographic analysis. The mobile phase was a gradient of water (Solvent A) and acetonitrile (Solvent B). Rutin and quercetin were used as standards. HPLC was examined at 257 nm to obtain real-time chromatograms of both standards and PDSE.

Cell line and culture
Human triple-negative MDA-MB-231 breast cancer, human ER and PR positive MCF-7 breast cancer, human liver HepG2 carcinoma , and normal kidney epithelial Vero cell lines were purchased from the NCCS, Pune, India. Cells were cultured in DMEM:F12 (1:1) medium in 25 cc tissue culture flasks in an incubator (Thermo Scientific, USA) as reported previously [24].

MTT assay
The anti-proliferative activity of PDSE was done by MTT assay using an established protocol [25]. Briefly, all three cancer cell lines viz. MDA-MB-231, MCF-7, HepG2 and one normal kidney epithelial cell line Vero were seeded at a density of 1 × 10 4 cells/mL in 96-well microtiter culture plates overnight. PDSE was dissolved in culture media for stock preparation and diluted in the same media at various concentrations ranging from 10, 25, 50, 75, and 100 μg/mL for treatment of cultured cells over a 24-and 48-h period. At the end of the incubation period, absorbance values were read using MTT dye (5 mg/mL stock) through an ELISA plate reader (Biorad-PW41, USA) at 550 nm. IC 50 values were calculated using GraphPad Prism software.

Nuclear condensation assay
The apoptosis-inducing effect of PDSE was assessed at two effective doses of 50 and 100 μg/mL. Hoechst 33258 staining was used to assess nuclear condensation, as described previously [23]. In brief, cells fixation was performed using 4% paraformaldehyde following exposure with 50 and 100 µg/mL of PDSE for 48 h. After permeabilization and then staining with Hoechst 33258 (5 μg/ mL) dye, cells were photographed under an inverted fluorescence microscope (Zeiss AxioVert 135, US).

Measurement of intracellular ROS
Flow cytometry technique was used to assess the intracellular ROS levels using DCFH-DA dye as described previously [23]. After 12 h exposure of PDSE, cell were washed with PBS and incubated in PBS containing 10 μM DCFH-DA dye at 37 °C for 20 min. The cells were then washed twice with PBS and subjected to flow cytometry analysis (FACS Canto II flow cytometer, BD Biosciences, USA).

Evaluation of mitochondrial membrane potential (MMP, ΔΨ m )
The fluorescent probe JC-1 was used to assess MMP changes as reported previously [23]. Briefly, cells were treated with 50 and 100 μg/mL of PDSE for 48 h. For flow cytometry analysis, cells were incubated with JC-1 dye at a final concentration of 2 μM for 30 min in the dark. After washing with PBS twice, cells were resuspended in 500 μL PBS and analyzed using flow cytometry.

Determination of apoptosis by Annexin V-FITC & PI double stain
Flow cytometry was used to quantify apoptotic cells using an Annexin V-FITC Apoptosis Kit (BioVision, USA) manufacturer's protocol. In brief, cells at a density of 1 × 10 6 cells/mL were incubated for 48 h with PDSE at concentrations of 50 and 100 μg/mL. Cells were then harvested, resuspended in binding buffer and stained for 15 min at 25 °C in the dark with 2 μL Annexin V-FITC and 2 μL PI. Flow cytometry was used to assess the apoptotic index.

Analysis of cellular DNA content
Cells at density 1 × 10 6 cells/mL were used to treat with 50 and 100 μg/mL concentrations of PDSE into a 6-well plate for 48 h. Flow cytometry was used to examine different stages of the cell cycle and the contents of cellular DNA, as described previously [26].

Western blot analysis
The western blotting of PDSE treated and untreated cells was carried out as per a previously published method [27]. Briefly, cell lysates were prepared in ice-cold RIPA lysis buffer containing protease and phosphatase inhibitor cocktail (Thermo Scientific). Protein sample (30 μg each) was resolved on a 10-15% SDS-PAGE gel, transferred to the nitrocellulose membrane (Millipore). Finally, the immunodetection was done using enhanced chemiluminescence (Millipore) as per manufacturer's instructions. The full-length blots were cut before antibody hybridization and each section was incubated with primary antibody individually. Image J software (version 1.43, NIH, USA) was used to quantify the relative abundance of each band against housekeeping β-actin protein.

Molecular docking analysis through AutoDock Vina and iGEMDOCK v2.1
The binding interaction(s) of PDSE active components viz. rutin and quercetin with the apoptosis executioner protein caspase-3 was performed using AutoDock 4.2 and iGEMDOCK v2.1 [24]. PubChem database was used to access the 3D structures of rutin and quercetin components with PubChem CID: 5280805 and 5280343, respectively. Energy minimization of phytocomponents was done by ChemBio3D Ultra 14.0, with Force Field type MM2. The 3D X-ray diffraction crystal structures of caspase-3 protein (PDB ID: 2XYP; Hetero 4-mer -A2B2) were downloaded from RCSB Protein Data Bank. Complete PDB structure with no mutation and resolution 1.86 Å was selected for molecular docking study. Before docking, the refinement procedure was carried out by the addition of missing atoms to the residues, addition of polar hydrogen atoms and Kollman charges, removal of crystallographic water molecules, and external and irrelevant ligands and ions from the protein. AutoDock Vina reduces the computational effort required for binding pocket predictions. This tool was used to run docking simulations and to generate ten ligand-receptor complex conformations that were then ranked based on binding energy. iGEMDOCK can be used for post-screening analysis and predicting pharmacological interactions from screening compounds. In case of iGEMDOCK v2.1, genetic algorithm parameters were as follows: population size = 200, generations = 0 and number of solutions = 2. The best fitted was then selected displaying total binding energy in the form of van der Waals forces (VDW), hydrogen bond (HB) and electrostatic interactions (EI). AccelrysBiovia Discovery Studio 2017 R2 (Biovia, San Diego, CA, USA) was used to visualize the best docking sites and poses from both docking simulations [24,28].

Statistical analysis
Cell viability data were expressed as the mean ± SEM of at least three independent experiments. GraphPad Prism software was used for statistical analysis, which comprised one-way ANOVA and Dunnett's Multiple Comparison Test (Version 5.01). A p-value of less than 0.05 was determined statistically significant.  Table S1. The HPLC chromatographic analysis provided a fine separation of rutin and quercetin with R t value of 13.632 and 19.049 min, respectively at 257 nm in chromatograms (Fig. 1f ). The corresponding peak of rutin and quercetin in PDSE was found with R t of 13.401 and 19.573 min, respectively under similar conditions (Fig. 1e). This study revealed the presence of rutin and quercetin as active components in PDSE.  (Fig. 4c). The cell viability data suggested that PDSE treatment significantly reduced cancer cell growth in both doses-and time-dependent manner. However, PDSE did not exert any significant morphology variation and effect on survival of normal cell line Vero as represented in Fig. 5a

Apoptotic body formation by PDSE
Based on cell viability data, two effective doses 50 and 100 μg/mL were selected to study the chromatin condensation in different cancer cell lines at 48 h incubation period. As apparent from the photomicrograph (Fig. 2f ), PDSE at 50 μg/mL concentration increased the chromatin condensation in MDA-MB-231 cells as compared to control, however,

PDSE causes intracellular ROS production
ROS measurement through flow cytometry analysis revealed that the mean fluorescence intensity of DCFDA dye in MDA-MB-231 control cells was 185.47 which was decreased to 95.71 and 68.11 at 50 and 100 μg/mL of PDSE treatment, respectively indicating a decrease in intracellular ROS levels ( Fig. 6a and b).

Apoptosis quantification at early and late stage
Apoptosis Detection Kit (Annexin V-FITC) was used further to examine the early and late apoptotic cells. Untreated cells demonstrated negligible apoptosis and dead cells, while PDSE at 50 μg/mL promoted cell death by lowering the number of viable cells and enhancing the percentage of early (19.8%) and late (18.97%) apoptotic cells. The concentration 100 μg/mL promoted the proportion of early and late apoptotic cells to 21.5 and 19.61%, respectively ( Fig. 6c and d). This indicates that a high dose of PDSE induced early and late apoptosis in treated cells.

Decrease in MMP by PDSE
As evident from flow cytometry data (Fig. 7), PDSE treatment resulted in an increase in green fluorescence by 62.61% at 50 μg/mL and 78.15% at 100 μg/mL doses

Effects of PDSE on anti-apoptotic and pro-apoptotic proteins
To analyze the underlying mechanisms of PDSE-induced cell death, western blotting was performed for the expression analysis of key apoptotic proteins viz. proapoptotic Bax, tumor suppressor p53, anti-apoptotic Bcl-2, effector caspase-3 and cleaved PARP-1. Results showed that pro-apoptotic protein Bax was upregulated, but it was independent of tumor suppressor protein p53 while anti-apoptotic protein Bcl-2 was down-regulated  (Fig. 9). These results indicate that PDSE induces apoptosis through intrinsic apoptotic pathways in breast cancer cells.

Binding analysis of rutin and quercetin present in PDSE against caspase-3 protein
The docking tools AutoDock Vina and iGEMDOCK v2.1 were used to study the binding interaction of rutin and quercetin present in PDSE with a caspase-3 target protein. The docking results obtained from both docking tools were visualized by BIOVIA Discovery Studio software. AutoDock Vina is based on the statistical scoring function that replaces the semi-empirical free energy force field of AutoDock 4.2. AutoDock Vina provides improved prediction accuracy and speed, which is not only due to the simplification of the scoring function but also due to the capability of multi-threading in presence of multiple CPU cores. Figure 10a and b depict the molecular structure of rutin and quercetin, respectively; Fig. 10c represents the 3-D crystal structure of a caspase-3 protein; Figs. 10d and Fig. 10e depict the docking interaction of rutin and quercetin complexed with caspase-3 protein, respectively as analyzed by Auto-Dock Vina while Fig. 10f and g are destined for rutin and quercetin, respectively as analyzed by iGEMDOCK v2.1 tool. Table 1 represents the active constituents of PDSE with their chemical structure, binding energy, dissociation constant, best docking poses with amino acid residues contributing to the binding pocket of the caspase-3 protein. As shown in Table 1, rutin and quercetin phytoconstituents exhibited potent binding interaction with caspase-3 protein. As analyzed by AutoDock Vina, the binding energy of rutin (BE = − 9.1 kcal/mol) was found to be lower than the binding energy of quercetin (BE = − 7.6 kcal/mol), which means that the binding affinity of rutin with caspase-3 protein is greater than quercetin. These results were further confirmed from analysis using iGEMDOCK v2.1, which showed consistent results with both components. The result suggests that both ligands have different binding sites and hence possess different binding energies towards caspase-3 protein.

Discussion
Cancer is a foremost health problem that affects millions of people across the world. At present, chemotherapy is considered one of the most effective cancer treatment strategies. Although chemotherapy significantly improves the prognosis of cancer patients, but this approach is not free from side effects. Chemotherapeutic drugs can damage normal, healthy cells resulting in side effects like hair loss, anemia, sores, nerve, muscle as well as kidney and fertility problems [29]. To avoid side effects, most cancer patients adopt alternative therapy based on herbal medicines. Plant-based cancer treatment is considered a better treatment option because phytotherapeutic agents are natural, readily available, easily assimilated in the body, and have fewer side effects and toxicity. Medicinal plants are well known for their antioxidant and immunomodulatory properties, as well as their anticancer activities [9,30]. The seeds of Ajwa dates are abundant in minerals, vitamins, dietary fibers, phenolic compounds, and different flavonoids [17]. Evidence indicates that individuals with a high intake of dietary fiber and phenolic compounds have reduced incidence of colorectal, prostate, lung, breast and ovarian cancer(s) [31][32][33]. Previous studies have shown that date palm seeds induce immunity in broiler chickens [34], cause an increase in the paraoxonase and arylesterase activities in hypercholesterolemic rats [35] and enhance the endogenous insulin secretion in type 1 diabetic rats [36]. An earlier study has also found that a combination of black pepper and Ajwa seed extract normalizes glucose levels and liver enzymes aspartate transaminase, alanine transaminase and alkaline phosphatase activities in alloxan-induced diabetic rats [37]. The therapeutic potential of PDSE has been investigated against DNA damage induced by N-nitroso-N-methyl urea (NMU) in mice [38] and carbon tetrachloride (CCl4)-induced hepatotoxicity in rats [15]. The protective role of PDSE against gastric ulcers has also been investigated [39]. Two ex vivo studies have been carried out to evaluate the efficacy of date seed oil extract to prevent oxidative stress; a major contributor towards cancer development [40,41]. Moreover, a study has shown that acetone extract of date palm seeds is highly cytotoxic against human colorectal cancer cell line HCT-15 and has significant antibacterial activity against Bacillus cereus and Escherichia coli [42]. However, further experimental work needs to be performed to confirm the previous traditional applications of Ajwa date seeds in cancer treatment. The current study not only presents a novel approach to understand the anticancer activity but also augments the existing knowledge about the traditional use of the PDSE. Although previous studies have reported the antidiabetic, hypolipidemic and antioxidant properties of PDSE [43,44], the cytotoxic effects of PDSE against human breast cancer cell lines MDA-MB-231 and MCF-7 and human liver cancer cell line HepG2 remain to be investigated. This study also attempted the HPLC characterization of PDSE to find the major active component(s) that might be contributing to the anticancer potential of the extract. Moreover, in silico molecular docking analysis between active components, viz. rutin and quercetin with apoptosis executioner caspase-3 protein further validated the anticancer potential of PDSE. For a comprehensive summary, Fig. 11 summarizes the phytochemical analysis, several in vitro anticancer parameters, and in silico analysis. The cell viability data indicated that PDSE had a cytotoxic effect against MDA-MB-231 cells with IC 50 values of 101.6 and 85.86 μg/mL following 24 and 48 h cells exhibit an estrogen-independent state and do not express estrogen receptors and hence they are ideal models for chemotherapeutic studies, however, MCF-7 cells possess estrogen and progesterone receptors and hence they are suitable models for investigations on hormone therapy [45]. Based on this principle, it can be concluded that PDSE could be a better therapeutic agent for the growth inhibition of triple-negative breast cancer cell line MDA-MB-231. A previous study has reported that hydromethanolic extract of Ardisia crispa showed moderate cytotoxic effect against MCF-7 and a weak cytotoxic effect against MDA-MB-231 [46]. Conversely, another study has stated that hexane fraction of Acanthopanax sessiliflorus stem bark extract displayed more cytotoxicity against MDA-MB-231 cells compared to MCF-7 cells [47]. Interestingly, in agreement with this previous study, ethanolic PDSE was found to be more cytotoxic against MDA-MB-231 cells as compared to MCF-7 cells. In addition, PDSE did not exert any significant morphology variation and toxic effect on survival of normal cell line Vero as represented (Fig. 5). Furthermore, the present study involved an investigation into the mechanism(s) responsible for PDSE mediated cytotoxicity. There are two major pathways of cell death viz. apoptosis or necrosis. Results of nuclear condensation by Hoechst staining revealed that PDSE-treated cancer cells exhibited the  (Figs. 2f, 3d, and 4d). Previous studies have also reported the apoptotic effect of several plant extracts against different cancer cell lines [48][49][50].
In cancer cells, high levels of ROS have been detected due to increased metabolic and peroxisomal activities, mitochondrial dysfunction, increased receptor signaling, oncogenic activity, increased enzymatic activity of oxidases, lipoxygenases, cyclooxygenases and thymidine phosphorylase [51]. In this study, PDSE decreased the level of ROS in treated MDA-MB-231 cells, which indicates that PDSE induced cell death independent of ROS pathway. ROS-independent cell apoptosis pathway has rarely been reported in cancer cells and may act as a natural regulator of important signaling pathways in cells [52]. Therefore, it can be expected that ROS signaling does not influence the mechanisms of PDSE-induced cell death. Interestingly, in earlier research, several prospective compounds viz. metformin, quercetin, curcumin and vitamin C have been found to downregulate ROS in the cellular apoptotic process and some of them have even been demonstrated to promote apoptosis in cancer cells [53]. Further, based on initial study of nuclear condensation, MDA-MB-231 cells were stained with Annexin-V-FITC/PI double stain and examined by flow cytometry. Results showed that the percentage of live cells was decreased with a simultaneous increase in early and late apoptosis of the cells. This study suggests that high dose of PDSE increases the percentage of cells undergoing late apoptosis (Fig. 6c). Similarly, in a previous study, methanolic extract of Ajwa dates pulp has been reported to increase the percentage of MCF-7 cells undergoing late apoptosis [54].
Mitochondria is not only the ATP factory of energy but also help in the regulation of the membrane potential, apoptosis, calcium signaling and regulation of cellular metabolism [55]. The efficacy of cancer therapy can be improved by altering the cellular metabolism of tumor cells or promoting MMP decrease. Permeabilization of the inner mitochondrial membrane causes disruption of MMP and thus, a decline of MMP is associated with the opening of a mitochondrial permeability transition pore which results in the rupture of the outer mitochondrial membrane and release of various apoptosis factors such as Cyt c, Smac, Endo G into the cytoplasm, finally leading to cell apoptosis [56]. The present study revealed that PDSE declined the MMP level by changing the JC-1 fluorescent color from red to green, suggesting that depletion of MMP is associated with the increasing concentration of PDSE. Cell cycle analysis showed that PDSE treatment resulted in S phase arrest of cell cycle followed by apoptosis in MDA-MB-231 cancer cells (Fig. 5). This study confirmed that PDSE hampers the initiation of DNA replication and thus arrests MDA-MB-231 cells at the S phase. Apoptotic events also require permeabilization of the outer membrane of the mitochondria which is regulated by the Bcl-2 family proteins. Based on the role of Bcl-2 family proteins, this study was also formulated to analyze protein expression of p53, Bax, Bcl-2, cleaved caspase-3 and PARP-1 cleavage in PDSE-treated cancer cells. Western blot analysis revealed that PDSE increased the expression level of pro-apoptotic Bax and effector cleaved caspase-3, whereas PDSE downregulated the anti-apoptotic Bcl-2 protein level (Fig. 9). While, PDSE did not reduce the tumor suppressor p53 protein, indicating PDSE induced cell death via p53 independent pathway. Following activation, Bax and Bak form homooligomers deactivating Bcl-2 protein and contribute to pore formation which causes permeabilization of the outer mitochondrial membrane, leading to the release of mitochondrial inner membrane space contents, including Cyt c and Smac, into the cytosol [57]. These contents drive the activation of apoptotic effector caspases including caspase-3 that cleave and degrade the crucial PARP protein, a DNA repair enzyme leading to DNA breakage and cellular apoptosis [58]. The apoptosis executioner enzyme caspase-3 is the main proteolytic cascade involved in the apoptosis of both intrinsic and death receptor pathways [59]. Western blot data indicated that PDSE induced cellular apoptosis through intrinsic pathways in breast cancer cells.
PDSE was further characterized using HPLC to identify the phytoconstituents present therein. Rutin and quercetin were the principal phytochemicals identified in PDSE. Rutin and quercetin are the most widely distributed plant flavonoids and major constituents of the human diet [60,61]. Previous studies have confirmed the anticancer activities of both bioactive agents against different human cancer cell lines [62][63][64]. In our recent publications, the pulp extract of Ajwa dates has shown IC 50 values of 20.03 and 16.78 mg/mL at 24 and 48 h periods, respectively, against human liver cancer HepG2 cells and IC 50 values of 17.45 and 16.67 mg/mL at 24 and 48 h, respectively against TNBC MDA-MB-231 cells [23,65]. While the active components of Ajwa date pulp extract in the previous study was β-D-glucan (a polysaccharide), rutin and quercetin were the active phytoconstituents identified in PDSE. Lower doses of PDSE were required as compared to Ajwa dates pulp extract, which suggests that seed extract has comparatively more potential against cancer cells as compared to pulp extract of Ajwa dates. Furthermore, an in silico computational study was carried out using chemoinformatic tools to validate caspase-3 activation. As caspase-3 is the main executioner protein accountable for apoptosis and also shares many structural characteristics with other-known caspases [59], thus, caspase-3 was selected as a potential target protein for ligand-protein binding affinity. Moreover, caspase-3 is active over a broad pH range [66], which indicates that caspase-3 would be fully active under normal and apoptotic cell conditions. Caspase-3, (also called apopain) is synthesized in the cell in its zymogen form of 32 kDa, consisting of an N-terminal pro-domain followed by a large 17 kDa (p17) and small 12 kDa (p12) subunit linked to each other by an inter-subunit linker. Caspase-3 in its functional form is a heterotetramer; formed by hydrophobic interactions of four anti-parallel beta-sheets from p17 and two from p12 subunits. Beta-sheet interacts with another heterodimer resulting in a 12-stranded beta-sheet structure, around which alpha-helices are positioned. A previous study has shown that catalytic residues of caspase-3 consist of sulfhydryl group of Cys-163 and the imidazole ring of His-121 [67]. The large subunit p17 harbors the active site catalytic dyad residues and the small subunits contain most of the dimer interface and the allosteric site [68]. Interestingly, as revealed by in silico binding interaction data using AutoDock Vina and iGEMDOCK v2.1 tools, rutin present in PDSE did not interacted with any catalytic residue, while quercetin interacted with both catalytic residues His121 and Cys163 in addition to other amino acid residues in the binding pocket of caspase-3 ( Table 1). As is apparent from binding interactions of amino acid residues in both AutoDock Vina and iGEM-DOCK v2.1 analyses, the slight variations in interacting amino acid residues are because of the differences in the grid box generation and determination of binding pockets on the target protein [28].

Conclusion
The current study revealed the cytotoxic activity of PDSE against breast and liver cancer and apoptotic potential against TNBC MDA-MB-231 cells. This activity could be attributed to the synergistic effect of bioactive agents particularly rutin and quercetin present in PDSE. Moreover, in silico analysis confirmed the potential binding affinity of rutin and quercetin with amino acid residues of caspase-3 executioner protein. However, because caspase-3 is a downstream protein, the apoptotic effect of PDSE in upstream pathways cannot be rationalized by binding these bioactive compounds to caspase-3. As a result, further study is warranted for upstream targetprotein(s) docking studies, as well as the isolation of bioactive compounds from PDSE and their mechanisms of action in vitro and in vivo. These studies would lead to a more complete assessment of PDSE and its feasibility as a future anticancer drug candidate and an adjunct.