Next Article in Journal
Hansen Solubility Parameters Applied to the Extraction of Phytochemicals
Next Article in Special Issue
Eugenol: In Vitro and In Ovo Assessment to Explore Cytotoxic Effects on Osteosarcoma and Oropharyngeal Cancer Cells
Previous Article in Journal
Transcription Factor GmERF105 Negatively Regulates Salt Stress Tolerance in Arabidopsis thaliana
Previous Article in Special Issue
New Insights into the Latest Advancement in α-Amylase Inhibitors of Plant Origin with Anti-Diabetic Effects
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 16
10.3390/plants12163006
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An In Vitro Evaluation and Network Pharmacology Analysis of Prospective Anti-Prostate Cancer Activity from Perilla frutescens

by
Patrick Jay B. Garcia
1,2,†,
Steven Kuan-Hua Huang
3,4,5,†,
Kathlia A. De Castro-Cruz
1,
Rhoda B. Leron
1 and
Po-Wei Tsai
3,*
1
School of Chemical, Biological, and Materials Engineering and Sciences, Mapúa University, Intramuros, Manila 1002, Philippines
2
School of Graduate Studies, Mapúa University, Intramuros, Manila 1002, Philippines
3
Department of Medical Science Industries, College of Health Sciences, Chang Jung Christian University, Tainan 711, Taiwan
4
Division of Urology, Department of Surgery, Chi Mei Medical Center, Tainan 711, Taiwan
5
School of Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
*
Author to whom correspondence should be addressed.
†
These authors have equally contributed to this study.
Plants 2023, 12(16), 3006; https://doi.org/10.3390/plants12163006
Submission received: 26 July 2023 / Revised: 12 August 2023 / Accepted: 18 August 2023 / Published: 21 August 2023
(This article belongs to the Special Issue New Insights in the Research of Bioactive Compounds from Plant Origin with Nutraceutical and Pharmaceutical Potential II)
Browse Figures
Versions Notes

Abstract

:
Perilla frutescens (L.) Britt. is extensively cultivated in East Asia as a dietary vegetable, and nutraceuticals are reportedly rich in bioactive compounds, especially with anticancer activities. This study explored the in vitro cytotoxic effects of P. frutescens parts’ (stems, leaves, and seeds) extracts on prostate cancer cells (DU-145) and possible interactions of putative metabolites to related prostate cancer targets in silico. The ethanol extract of P. frutescens leaves was the most cytotoxic for the prostate cancer cells. From high-performance liquid chromatography analysis, rosmarinic acid was identified as the major metabolite in the leaf extracts. Network analysis revealed interactions from multiple affected targets and pathways of the metabolites. From gene ontology enrichment analysis, P. frutescens leaf metabolites could significantly affect 14 molecular functions and 12 biological processes in five cellular components. Four (4) KEGG pathways, including for prostate cancer, and six (6) Reactome pathways were shown to be significantly affected. The molecular simulation confirmed the interactions of relevant protein targets with key metabolites, including rosmarinic acid. This study could potentially lead to further exploration of P. frutescens leaves or their metabolites for prostate cancer treatment and prevention.

1. Introduction

Prostate cancer (PCa) affected 14.1% of men worldwide, particularly old-aged men, in 2020 [1]. In this report, it is the fifth most death-causing cancer disease. Most PCa cells are potentially metastatic and are sometimes detectable from the emergence of urethral blockage or hematuria [2]. These are usually late-stage manifestations that allude to appropriate monitoring, especially for high-risk individuals. The early detection of PCa is currently tested with prostate-specific antigen (PSA) levels in blood serum [3]. However, it is reported to have broad-ranged screening specificity and sensitivity, and PSA levels are increased for patients with benign prostate hyperplasia (BPH), which is usually a noncancerous prostate gland enlargement [4]. For this reason, additional tests such as biopsy or the detection of other biomarkers are required. Current non-metastatic hormone-sensitive PCa (nmHSPC) treatments are hormonal therapy, radiotherapy, immunotherapy, or castration/orchiectomy [3,5]. However, metastatic HSPC (mHSPC) and castrate-resistant PCa (mCRPC), which are advanced PCa, are unresponsive to such hormonal therapy and surgical procedures, which requires life-maintaining treatment and palliative care since they are generally incurable [5,6]. It is a wide-ranging disease that follows particular cell growth and progression mechanisms with complex molecular pathogenesis and etiology, which require multi-targeted treatments [5].
An increasing number of investigations are being conducted of PCa for traditional and complementary medicines, solely or in decoctions, as supplemental to such treatments or as a preventative measure [7,8,9]. These medicinal systems are commonly validated with non-physical mechanisms and alternative anatomical interpretations, which are viewed in obscurity and discord with biological mechanisms to some extent. Current efforts in traditional and complementary medicine research to uncover molecular mechanisms and therapeutic effects have increased the confidence of health consumers in this medicinal system [7,10,11,12]. Still, natural products have been a great source of biologically active molecules in conventional medicine, such as the chemotherapeutic drug Paclitaxel and Taxifolin from Taxus chinensis (Chinese yew) [8,11,13]. While plant materials may be low in concentration of these bioactive molecules compared to current medications, potent plant materials and extracts have been discovered to be indicative of excellent hit compounds [14,15]. Previously, it was reported that conjunctive use of traditional Chinese medicine (TCM), especially the Chai-Hu-Jia-Long-Gu-Mu-Li-Tang decoction, increased the survivability of metastatic PCa patients [16].
One of these traditional medicines is the Perilla frutescens (L.) Britt., which is an annual plant native and widely cultivated in eastern, southeastern, and southern Asia, mainly in subtropical climates, whereas it was introduced in northeastern America and southern Europe [17]. This plant has excellent nutraceutical and pharmaceutical value in TCM, and its leaves are a staple food in East Asian cuisine [18]. This plant has been traditionally used for but is not limited to athlete’s foot, generalized edema, cough, antibiotic, antidepressants, anti-lung, and colon cancer [19]. In fact, Jeong et al. [20] confirmed that ethanol extracts of this plant inhibited the metastatic ability of cancer cells, specifically breast and liver cancer cells, through Src kinase deactivation, thereby blocking the epithelial–mesenchymal transition (EMT) process. Also, a study by Lee et al. [21] uncovered the activity of isoegomaketone, an isolate in P. frutescens oils, as an inducer of apoptotic mechanisms on PCa cells from the activation of death receptors 4 and 5. Additionally, numerous secondary metabolites from plants’ adaptive mechanisms in different conditions could be a good source of compounds for multiple targets.
Network pharmacology is the current paradigm of molecular investigations of herb and herbal formulae in traditional medicines [22,23,24]. This computational approach includes the chemical diversity of TCMs into multiple targets through experimental data and predictive models aimed at advances in polypharmacology and multi-targeted drug discovery approach [23,24]. This network-driven approach overcomes the barriers of the conventional one-drug-one-target scheme and, therefore, is appropriate for complex diseases, mainly due to the diversity of cancer mechanisms. Based on the common targets between the disease and active compounds from the TCM herbals, this network-based method can characterize potential key molecular targets and mechanisms of action.
The abundance of compounds with anticancer activities in various parts of P. frutescens var. frutescens (PF), the green variety, has potential applications in PCa treatment or prevention [25]. Furthermore, this study focused on polar compounds, including phenolic and flavonoid compounds, which are abundant in water and ethanol extracts of P. frutescens and commonly have anticancer properties [26,27,28,29]. The cytotoxic effects of these extracts were explored with prostate adenocarcinoma (PRAD) DU-145 cell lines. The metabolites’ action was briefly studied with in silico network pharmacology. There is limited information about their potential anti-PCa activity; hence, this study addressed initial interests in the bioactivity of the plant’s extracts as supplemental to PCa treatment or its prevention. Also, understanding the molecular mechanism of the plant could potentially decipher drug-to-herb interaction as a necessary criterion for their adjuvant or concomitant use with conventional medicines.

2. Results

2.1. Total Phytochemical Content

The apparent phytochemical contents of the PF parts relevant to anticancer activity were assessed from total phenolics and flavonoid contents as shown in Table 1. The leaf extracts revealed superior amounts of phytochemicals compared to the other studied parts. For the total phenolics content with a Folin–Ciocâlteu assay, PF-L-W > PF-L-E > PF-SD-W > PF-S-E > PF-S-W > PF-SD-E, whereas for total flavonoid content with an AlCl3 assay, PF-L-W > PF-L-E > PF-SD-E > PF-S-E > PF-SD-W. PF leaves (PF-L) and water extracts were observed to obtain the highest amount of such phytochemicals.

2.2. Anti-Prostate Cancer Activity

Dose–response behavior was observed for all samples; however, the water extract of PF seeds was not significantly different from the negative control and had no dose–response behavior. In Figure 1, the 5.00 mg/mL dose of ethanolic PF leaf and seeds and aqueous leaf extracts surpassed the cytotoxicity of 5-fluorouracil (5FU). The highest half-maximal cytotoxic concentrations were observed from PF-L-E (1.3982 ± 0.3453 mg/mL) > PF-L-W (1.6224 ± 0.2926 mg/mL) > PF-SD-W (2.0612 ± 0.0327 mg/mL) > PF-SD-E (2.1631 ± 0.7595 mg/mL) > PF-S-E (3.5275 ± 0.4574 mg/mL). The PF part extract with the most cytotoxic potential of DU-145 cells is found in their leaves, specifically for the ethanol extract.

2.3. Detection of the Major Compound

Metabolomic analyses of PF-L show that its major metabolite is rosmarinic acid (RA), a phenolic acid [30,31]. For verification, the presence of rosmarinic acid in the crude PF-L extracts was evaluated with HPLC analysis. Peak identification revealed its retention time at tR = 29.82 min in Figure 2.

2.4. Network Pharmacology

Putative metabolites of PF-L were collected from the literature as summarized in Table S1, excluding the nonpolar compounds [30,31,32,33]. There were 358 predicted targets of these metabolites after duplicate removal. From differentially expressed genes of PRAD versus normal prostate cells, 14 upregulated and 46 downregulated genes are associated with cytotoxicity of PF-L metabolites to DU-145 cells, shown in Figure 3a. The higher degree of node connections’ strength and amount in these network analyses indicates its importance. These nodes with higher connection degrees could potentially affect multiple targets simultaneously. In the metabolite-target network (MTN) for these common targets in Figure S1, simplified in Figure 3b, with at least 20 affected targets were tuberonic acid (37, 20 targets), scutellarein-7-O-glucuronide (32, 20 targets), 5’-gluco-pyranosyoxyjasmanic acid or tuberonic acid glucoside (1, 20 targets), and n-octanoylsucrose (23, 21 targets). The major compound, RA (28), in this MTN analysis had 17 common targets.
Interaction of the proteins from the protein-coding genes seen in Figure S2 revealed a highly significant network (p = 9.04 × 10−12) with 76 connections, an average node degree of 2.53, an average local clustering coefficient of 0.382, and a maximum interaction score of 0.50. An increased required interaction score of 0.90 revealed a significant network (p = 0.00033) with 19 significant interactions, an average node degree of 0.633, and an average local clustering coefficient of 0.298 (Figure 4). Henceforth, the top five (5) targets from topological analysis in Cytoscape with maximal clique centrality (MCC) and density of maximum neighborhood component (DMNC) algorithms were considered key targets for PCa cytotoxicity of PF-L metabolites. The MCC and DMNC algorithms were previously observed to best identify essential targets in a complex interactome [34]. In this analysis, the p59-Fyn proto-oncogene (FYN), platelet-derived growth factor receptors α and β (PDGFRA and PDGFRB), phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1), and integrin αVβ3 (ITGB3) were top-ranking. Based on the network in Figure S3, metabolites with at least three (3) affected PPI network-relevant targets were 5′-gluco-pyranosyoxyjasmanic acid (1, 3 targets), apigenin-7-O-diglucuronide (5, 3 targets), coumaric acid-4-O-glucoside (15, 3 targets), perillaldehyde (24, 3 targets), trans-p-menth-8-en-yl caffeate (36, 4 targets), and tuberonic acid (37, 5 targets). From these analyses, compounds 1 and 37 were defined for both MTN.
The gene ontology (GO) terms classified expressed properties from genes and their products according to their molecular functions (MF), involvement in biological processes (BP), and presence in cellular components (CC). In Figure S4, GO terms and pathways were ordered according to the adjusted p-value from the g:SCS algorithm. A higher p-value indicates that it has a significant PCa cytotoxicity effect from the metabolite action on protein-coding gene targets. According to the enriched GO terms, 14 MFs, 12 BPs, and 5 CCs were considered significant from 60 intersecting gene targets. Biological processes relevant to the PCa cytotoxic effect of PF-L metabolites are the responses to oxygen-containing compounds, protein phosphorylation, and cellular component disassembly, mainly distributed in the cytoplasm, cell surface, and periphery. These results reveal that the metabolites may influence the activity of ions, protein kinase, and identical protein binding.
The apparent mechanism of action of PF-L activity on PCa may be deciphered from significant pathways from KEGG and Reactome. Furthermore, pathway enrichment analysis showed four (4) KEGG pathways and six (6) Reactome pathways, which involved 28 targets, which are acetyl-CoA carboxylase 1 (ACACA), acetyl-CoA carboxylase 2 (ACACB), acetylcholinesterase (ACHE), disintegrin and metalloproteinase domain-containing protein 10 (ADAM10), aldose reductase (AKR1B1), 20α/3α-hydroxysteroid/dihydrodiol dehydrogenase (AKR1C1), 3α-hydroxysteroid 3-dehydrogenase (AKR1C2), tyrosine-protein kinase receptor UFO (AXL), calcium/calmodulin-dependent protein kinase kinase 2 (CAMKK2), M-phase inducer phosphatase 2 (CDC25B), fructose-1,6-bisphosphatase 1 (FBP1), tyrosine-protein kinase Fyn (FYN), non-lysosomal glucosylceramidase (GBA2), glutathione S-transferase P (GSTP1), cortisone reductase (HSD11B1), integrin αVβ3 (ITGB3), protein Mdm4 (MDM4), gelatinase A (MMP2), gelatinase B (MMP9), nuclear receptor 4A1 (NR4A1), platelet-derived growth factor receptor α (PDGFRA), platelet-derived growth factor receptor β (PDGFRB), phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1), membrane associated phospholipase A2 (PLA2G2A), protein kinase C-alpha (PRKCA), prostaglandin G/H synthase 2 (PTGS2), retinoic acid receptor RXR-β (RXRB), and stearoyl-CoA desaturase (SCD). The first two pathways were cell growth and proliferation-related pathways; there are three (3) cancer-related pathways, two (2) metabolism-related pathways, a cell segregation-related pathway, and an angiogenesis-related pathway. Based on the protein targets affected visualized in Figures S5–S14, PF-L metabolites have diverse targets where their activities are still in need of further validation and research.

2.5. Molecular Docking

For protein–ligand interaction validation, selected protein targets and metabolites from network pharmacology results were subjected to molecular docking. The binding site was selected with the volume of co-crystallized ligands for FYN (PDB code: 2DQ7), PIK3R1 (PDB code: 4JPS), ITGB3 (PDB code: 6MK0), and PDGFRA (PDB code: 6JOK). At the same time, the most significant volume from the grid-search and eraser algorithm was determined for PDGFRB (PDB code: 3MJG). Controls used were 5-fluorouracil (CID: 3385), [RGD-ChgE]-CONH2 [35], alpelisib (CID: 56649450), nintedanib (CID: 135423438), ponatinib (CID: 24826799), and sunitinib (CID: 5329102). The highest docking score from Table S5 for each ligand was selected, shown in Figure 5. RA (28) has the best docking score compared to 5FU. For metabolites in network pharmacology, all proteins were most stably docked with scutellarein-7-O-glucuronide (SG, 32). These compounds were more stable and spontaneously binding than the positive controls except for [RGD-ChgE]-CONH2. Prevalent drug–receptor interactions, such as hydrophobic contact, hydrogen bonding, and π-stacked interactions [36], were observed in all protein-28 complexes shown in Figure 6.

3. Discussion

PF has a wide array of pharmacologically active compounds historically used in TCM. Phytochemicals, such as polyphenols, are secondary metabolites of plants that are generally responsible for their defense mechanism against external stresses, such as infection, predators, environmental changes, and, to some extent, their growth [37]. Most polyphenols are anticancer agents targeting apoptotic and cell cycle pathways, along with cancer cell proliferation, tumorigenesis, angiogenesis, and metastatic abilities [26,27,28,29]. Plant extracts are an excellent source of these cancer-curative and preventative compounds due to their diverse targets and potentially synergistic action against cancer. Since PF-L had the most total phenolic and flavonoid content, its anti-PCa activity could be superior to other parts of PF.
Cancer is a multifactorial disease triggered through genetic alterations and instabilities, epigenetics, weakened immunosurveillance, infection, and so on [38]. Hallmarks of cancer pathogenesis revolve around the heightened activity of oncogenes (like growth-related genes responsible for tumor proliferation) and contrarywise for tumor suppressor genes [39]. A multi-targeted approach is appropriate since this disease takes advantage of each cell type’s specific cellular and metabolic machinery. The development of PCa at every stage has apparent changes in cellular mechanisms [40]. Henceforth, numerous in vitro models for PCa, such as the DU-145, PC-3, and LNCaP cell lines, are utilized for anti-prostate cancer activity research [41,42]. The prostate adenocarcinoma DU-145 cells, which are castration-resistant prostate adenocarcinoma from central nervous system (CNS) metastasis, were utilized for screening in this study [43]. While it is a negative-androgen receptor-expressing cell line, current reports reveal that DU-145 has positive androgen receptor expression in mediocre levels compared to androgen receptor-expressing cell lines [44]. Nevertheless, the phenolic and flavonoid-rich crude ethanol extracts of PF-L exhibited the lowest required concentration to become 50% cytotoxic to DU-145 cells. It indicates that PF-L has deficient activity relative to potent anticancer plant extracts [45,46,47]. This activity could signify that metabolites at low concentrations are more bioactive or that major metabolites require high concentrations. Deciphering the probable molecular targets of these metabolites are explored with network pharmacology and molecular docking analyses. This network-based bioinformatic analysis of metabolite targets and PRAD revealed 60 affected targets. The results show that compounds 1, 5, 15, 23, 24, 32, 36, and 37 have the most relevant targets on PCa cytotoxicity.
Previous studies revealed that many of the main compounds and such derivatives have been able to impact PCa cell proliferation and growth. It was known that apigenin targeted PI3K/Akt/Fox-O, β-catenin, and insulin-like growth factor-I signaling pathways, impacting PCa proliferation [48]. Also, apigenin derivatives have been known to induce cell death and cell cycle arrest [49]. Coumaric acid is known to have a similar effect, but it has targeted mitochondrial-related apoptosis, downregulated specific cyclin-dependent kinases, and decreased the expression of specific oncogenes [50,51]. Perillaldehyde has been shown to reduce bone metastasis of PC-3 prostate cancer cells with repression of the nuclear factor-κB (NF-κB) pathway and receptor activator of NF-κB ligand (RANKL) [52]. At the same time, scutellarin is considered a promising cancer treatment due to its multifactorial effect of cell cycle arrest, apoptosis induction, and angiogenic and metastatic reduction [53,54,55,56]. However, there are limited studies regarding their structure analogs with PCa. The major compound of PF-L, RA, was observed to have an excellent activity to DU-145, which could be the main-effect compound of its extracts [57,58,59,60].
In this study, affected targets from PRAD with PF-L were subjected to a GO term, KEGG, and Reactome pathway enrichment analysis, which uncovered prospective cytotoxic mechanisms to PCa in silico. Relevant processes that could be involved in this phenomenon from GO enrichment are the response to oxygen-containing compounds, protein phosphorylation, cellular component disassembly, and to a lesser extent, angiogenesis. From pathway enrichment, PF-L metabolites potentially regulate the AMPK signaling pathway, EGFR tyrosine kinase inhibitor resistance, miRNAs in cancer, prostate cancer, steroid and lipid metabolism, SCF-KIT and PIK3/Akt signaling pathways, EPH-ephrin mediated cell repulsion, and the VEGFA-VEGFR2 pathway. Collectively, these pathways were discovered to significantly affect the tumorigenesis, growth, progression, energy regulation, metastasis, and adhesion of PCa cells [61,62,63,64,65,66,67,68,69].
One of the relevant protein-coding genes from the PPI network is Fyn, under the Src family kinases. This gene is typically overexpressed in PCa, a signaling protein typically attributed to tumor cell proliferation and metastasis [70,71]. It regulates diverse biological functions, such as cell growth, proliferation, migration, and adhesion. However, overexpressed Fyn mediates the PIK3/Akt pathway to induce anti-apoptotic mechanisms and EMT for cell invasion and metastatic initiation, preventing the cell death pathways of normal cells [70,71,72]. It is a recent interest for cancer targets due to their functions in cancer cells and their drug resistance ability; however, the latter is effective in direct gene expression regulation. Indeed, previously reported immunochemical assays for Fyn kinase activity of RA and caffeic acid exhibited noncompetitive inhibition [73,74,75].
The molecular action of PF-L metabolites was also implicated with ITGB3, a transmembrane integrin receptor primarily responsible for cell adhesion, migration, and macrophage phagocytosis ability [76,77]. It was previously reported that the antagonism of integrin αVβ3 or vitronectin receptor revealed the most effective reduction in tumor cell angiogenesis and metastasis compared to the other classes of integrins [76,77,78,79]. Since it is overexpressed in highly metastatic neuroendocrinal metastasis of PCa, it is an antiangiogenic target, and it affects the migration, survival, plasticity, and metastatic ability of tumor cells as signal initiators [76,78]. However, previous results for endometrial cells showed that exposure to aqueous extracts of PF-S&L corresponded to increased integrin β3 expression and no observable changes in integrin αV [80]. This result is unfavorable for cancer angiogenesis and should be considered for further study.
Another PPI network-significant protein-coding gene with PF-L is the PIK3R1 from the PIK3/Akt pathway, which typically regulates cell death, proliferation, metabolism, and angiogenesis [81]. It is an obligate heterodimer composed of the p110α catalytic subunit (PIK3CA) and the p85α regulatory subunit (PIK3R1) [82]. This protein-coding gene is underexpressed and altered in both primary and metastatic PCa cells and has inverse negative feedback with androgen receptor (AR) [81,83]. Hence, the dynamic interplay of PIK3R1 and AR indicates that the current androgen-deprivation treatment could enable cell proliferation and survivability, suggesting vital consideration of combinatorial inhibition of the two targets [81]. For PF-L, the PI3K/Akt/NF-κB axis was regarded as the primary pathway of inhibited migration and invasion, induced apoptosis, and cell cycle arrest from RA [75,84,85].
Lastly, PF-L metabolites were shown to have significant interaction with PDGFRs, which are transmembrane receptor tyrosine kinases that can initiate various pathways with the SH2 domain-containing signal transduction molecules as inducers of cell growth and division, such as in the PI3K/Akt pathway [86]. Dimerization of these receptor monomers with PDGF ligands induces autophosphorylation [86,87,88]. The PDGFRα and β were implicated in PCa due to their heightened expression in PCa cells and oncogenic activity on its growth, angiogenesis, and metastatic potential [86,89,90,91]. In fact, ligand and dimerization-free PDGFRα are about as bone-metastatic compared to the complete form of the receptor for PCa, indicative of PDGF-independent mechanism on its metastatic ability [89,90]. In comparison, the high expression of PDGFRβ in prostate tumor stroma was associated with PCa aggressiveness and low patient survivability [92,93]. PDGF-induced proliferation of mesangial cells was inhibited upon exposure to RA [94].
Molecular docking was utilized to test whether the selected metabolites could interact with the key targets. This simulation revealed stable binding with respective proteins consistent with previous activities of their analogs. RA had the best activity, which could affirm the activity of PF-L on PCa. These results suggest that added structural features in the main compound may have increased their potency.
Anti-PCa activity from PF-L extracts can be promising, especially as a preventative measure and an additional diet for individuals with a high risk of PCa development. However, pharmacokinetics and toxicity with normal cells in vitro and in vivo tests are required to successfully evaluate PF-L’s potential as an anticancer agent for PCa or other cancer types. Also, the polar and nonpolar metabolites present in PF-L could be synergistic or antagonistic for the anti-PCa activity of PF and could be addressed from compound isolation and comparative explorations. At the same time, predicted mechanisms and actions from bioinformatic analyses should be further validated and investigated with gene expression studies. Meta-analysis of PCa patients and people at risk with their diet on such TCMs could provide an expansive vantage point of its potential as a cancer curative or preventative agent. Interestingly, other environmental conditions for the plant may lead to unique secondary metabolites, hence different activity.

4. Materials and Methods

4.1. Preparation of Plant Extracts

Dried P. frutescens stem (PF-S), leaves (PF-L), and seeds (PF-SD) were acquired from a local TCM store located in Tainan City, Taiwan. The identification of P. frutescens was authenticated by Dr. Chia-Jung Lee, Ph.D. Program in Clinical Drug Development of Herbal Medicine, College of Pharmacy, Taipei Medical University. A voucher specimen was deposited as #2023-CJCU-PF-001 at the Department of Medical Sciences Industry at Chang Jung Christian University, Taiwan.
The samples were collected from manually separated leaves and stems, where the petioles were included with the latter. These samples were mechanically ground for surface area enlargement and particle size reduction for effective extraction. For ethanol extracts, dried samples (50 g) were immersed with one-liter 95% ethanol under a two-hour 65 °C water bath and total reflux. For water extracts, dried samples (50 g) were leached with one liter of distilled and deionized water in a traditional Chinese decoction pot until reduced to ~200 mL. The crude extracts were vacuum-filtered with 0.45 µm filter paper, concentrated by vacuum evaporation, and lyophilized to obtain a solvent-free extract stored in a −20 °C refrigerator for further analysis. The yield from the extraction is reported in Table 2.

4.2. Total Phytochemical Content Assay

The total phytochemical contents of PF extracts were assessed according to Tsai et al. [95] with slight modifications. A 1000 μg/mL stock solution of crude PF extracts was prepared for all assays. The chemicals used in the analyses are reagent grade. Standards were prepared from two-fold serial dilutions of 1000 μg/mL stock solution from 500 to 15.60 μg/mL. Results of the following assays were expressed in milligrams of standard equivalents per gram of crude extract. All measurements were taken in triplicates and color-corrected with blank solutions.

4.2.1. Total Phenolics Content (TPC) Assay

Gallic acid standards were prepared for TPC analysis. The standard and sample solutions (80 μL) were added with 400 μL of 0.2 N Folin–Ciocâlteu reagent. After five-minute equilibration, 320 μL of 7.5% (w/v) Na2CO3 was added. The mixture was incubated for 30 min at room temperature and transferred to a 96-well microarray plate. The absorbance of the mixtures was measured at 600 nm.

4.2.2. Total Flavonoids Content (TFC) Assay

Rutin standards were prepared for TFC analysis. The standard and sample solutions (500 μL) were reacted with 500 μL of 2.0% (w/v) AlCl3 reagent. The mixture was incubated for 15 min at room temperature and transferred to a 96-well microarray plate. The absorbance of the mixtures was measured at 430 nm.

4.3. Anti-Prostate Cancer Activity

4.3.1. Human Prostate Cancer Cell Line DU-145

The cell culture of DU-145 was adapted from Park et al. [96]. Briefly, human prostate cancer cell line DU-145 was obtained from Bioresource Collection and Research Center (BCRC, Taiwan). The DU-145 cells were cultured in Eagle’s Minimum Essential Medium (EMEM), containing 10% fetal bovine serum (FBS) and 1% penicillin-streptomycin. The cultured cells were kept at 37 °C in a humidified atmosphere containing 5% CO2. The cells were sub-cultured within two-day intervals after reaching 70–80% confluence.

4.3.2. Cell Treatment and Cell Viability with WST-1 Assay

DU-145 cells were cultured in 96-well microarray plates (2 × 104 cells/well) and incubated at 37 °C. After 24 h of incubation, the cells were treated with various concentrations (two-fold dilution from 5.00 to 0.625 mg/mL) of crude extracts for another 24 h. The positive control in this analysis was 1.00 mg/mL of 5-fluorouracil (5FU). Cell viability was determined using a WST-1 assay (AbcamTM) [97]. Briefly, 100 µL of a fresh medium was placed in each treated well with DU-145 cells. A total of 10 µL of WST-1 reagent was added to each well. The plate was incubated for 2 h at 37 °C. The absorbances were measured at 570 nm with background control. Measurements were performed in triplicates for statistical significance and reproducibility.

4.4. Chromatographic Analysis

The major compound was separated and identified in the extracts with liquid chromatography (LC-2050C 3D; Shimadzu Corporation, Japan) for verification through a Hypersil™ BDS C18 column (250 × 4.6 mm, 5 µm, Thermo Fischer Scientific Inc., Taiwan) kept at 25 °C. Chromatography conditions were adapted from Lee et al. [30] with slight modifications. A 5 μL of 10 mg/mL crude extracts dissolved in methanol was filtered with a 0.45 µm syringe filter and injected into the column. Chromatography runs were performed with mobile phases of 0.05% trifluoroacetic acid in water (eluent A) and methanol (eluent B). The implemented gradient program was 30% B (0–5 min); 30-50% B (5–20 min); 50–0% B (20-40 min); 90–100% B (40–45 min); and 100% B (45–50 min); then, it was equilibrated back to initial conditions for five minutes at a flow rate of 0.4 mL/min. Data acquisition was set until 45 min.

4.5. Data Treatment and Statistical Analysis

Microsoft Excel was used for data processing, and GraphPad Prism software was used for statistical analysis and data visualization. Experimental data are reported as mean ± standard deviation (SD). Multiple comparisons of means were conducted through one-way ANOVA and Dunnett’s pairwise comparison test. A statistical significance of p < 0.05 was set throughout the analysis.

4.6. Network Pharmacology

A brief prediction of the mechanism of action of PF with PCa cytotoxicity was performed with network analysis adapted from Jin et al. [98] and Peng et al. [99] with modifications. The methods and resulting data of this analysis were validated according to the guidelines of the World Federation of Chinese Medicine Societies [100].

4.6.1. Target Prediction and Identification

The experimental and predicted gene targets of previously identified putative metabolites from PF extracts were identified from SuperPred (prediction.charite.de, accessed in June 2023). This target prediction tool utilizes logistic regression and 2048-long Morgan fingerprints with 94.1% target prediction accuracy [101]. The default set cut-offs for predicted target probability and model accuracy were ≥ 50%. Differentially expressed genes between PRAD (ICD-11: 2C82.0) and normal prostate cells were determined from the Gene Expression Profiling Interactive Analysis (GEPIA2, gepia2.cancer-pku.cn, accessed in June 2023). Genomic datasets of the GEPIA2 web server utilized were from TCGA and GTEx isoform expression data [102]. Comparisons were performed with the LIMMA differential method. Statistically significant genes were selected with the set significance of p < 0.01 and |log2 FC|> 1. Gene nomenclatures were standardized into their official gene symbols using SynGO (www.syngoportal.org/convert, accessed in June 2023) [103]. Matched gene targets between the metabolites and PRAD were selected for further analysis with Venny 2.1 (bioinfogp.cnb.csic.es/tools/venny/, accessed in June 2023) and Intervene (asntech.shinyapps.io/intervene, accessed in June 2023) [104,105].

4.6.2. Network Construction

The identified matched protein-coding genes were imported to the STRING database 11.5 (string-db.org, accessed in June 2023) for the construction of a PPI network [106,107,108,109,110,111,112,113,114,115,116,117,118]. Furthermore, proteins with more than one (1) interaction and FDR stringency of 0.01 limited to the “Homo sapiens” species were considered significant. The PPI network is transferred CytoScape 3.9.1 for further analysis [119]. Protein relevance in the network was evaluated with topological analysis using MCC and DMNC algorithms from the cytoHubba plug-in in CytoScape [34]. The metabolite target network was constructed similarly and analyzed with the degree-ranking algorithm.

4.6.3. Enrichment Analysis

Functional profiling was performed with g:Profiler (biit.cs.ut.ee/gprofiler/gost, accessed in June 2023) for the exploration of the predicted pharmacological mechanism and signaling pathways involved in the action of metabolites to the selected targets [120,121]. Relevant results from databases of gene ontology (GO) terms [122,123], KEGG [124,125,126] and Reactome [127,128,129,130,131,132,133,134] pathways were characterized from a g:SCS significance threshold of less than 0.05 limited to the “Homo sapiens” species. Secondary data filtering of GO terms was automatically performed from g:Profiler with a simple greedy search algorithm. Data visualizations were performed in Python 3.11.

4.6.4. Molecular Docking Validation

The molecular modeling and visualization software BIOVIA Discovery Studio was employed for the docking analysis. Three-dimensional (3D) structure-data files (SDFs) of ligand molecules were collected from PubChem (chem.ncbi.nlm.nih.gov, accessed in June 2023). The ligands were prepared with the ligand preparation (at pH 7.5 ± 1.0) and minimization protocols in Discovery Studio. Experimentally elucidated structures of protein targets were collected in PDB file format from the Protein Data Bank (PDB, www.rcsb.org, accessed in June 2023). Afterward, water molecules and irrelevant heteroatoms were removed, and polar hydrogens were added to the protein structure. The proteins were prepared with the protein preparation protocol (at pH 7.4 and 0.145 M ionic strength) in Discovery Studio. These proteins were minimized until an energy gradient of 0.01 with CHARMm force field and the Momany-Rone partial charge estimation method. The defined binding sites were determined from the co-crystallized ligands volume or eraser algorithm in Discovery Studio [135]. The docking runs were performed using the grid-based CDOCKER protocol [136]. Receptor-ligand binding site visualizations were performed in PyMOL™ [137].

5. Conclusions

Most PF extracts follow dose–response behavior for PCa cell cytotoxicity with DU-145 cells, where its ethanolic leaf extracts had the most cytotoxic activity. While characterized as having minimal activity as an anti-prostate cancer agent, the plant may have preventative potential for PCa development and progression. Regardless, bioinformatic analyses revealed that many minor metabolites and the major metabolite, RA, could have contributed to the cytotoxic potential of the mentioned extract, affecting the determined key targets FYN, PDGFRA, PDGFRB, PIK3R1, and ITGB3. Molecular docking of the compounds to the key targets consistently verified the protein–ligand interactions from network pharmacology results. Hence, the cytotoxicity of the leaf extract can be postulated from cell death induction and cell cycle arrest from action to the key targets. More advanced studies, such as in vivo trials, gene expression, or clinical trials, are recommended to verify activities and mechanisms of action from PF-L and their metabolites to PCa.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12163006/s1, Figure S1: Metabolite-target network of P. frutescens leaves against PRAD; Figure S2: Protein-protein interaction network of intersected genes set; Figure S3: Metabolite-target network of PIN-relevant proteins; Figure S4: GO term enrichment analysis by (a) molecular function, (b) biological process, and (c) cellular component; (d) KEGG and Reactome pathway enrichment analysis; Figure S5: Affected proteins in the AMPK signaling pathway from the KEGG database (has04152); Figure S6: Affected proteins in the EGFR tyrosine kinase inhibitor resistance pathway from the KEGG database (hsa01521); Figure S7: Affected proteins from protein-translated miRNAs in cancer from the KEGG database (hsa05206); Figure S8: Affected proteins in the prostate cancer pathway from the KEGG database (hsa05215); Figure S9: Affected proteins in the steroid metabolism pathway from the Reactome database (R-HSA-8957322); Figure S10: Affected proteins in the SCF-KIT signaling pathway from the Reactome database (R-HSA-1433557); Figure S11: Affected proteins in the lipid metabolism pathway from the Reactome database (R-HAS-556833); Figure S12: Affected proteins in the PI3K/Akt signaling pathway in cancers from the Reactome database (R-HAS-2219528); Figure S13: Affected proteins in the EPH-ephrin mediated cell repulsion pathway from the Reactome database (R-HSA-2682334); Figure S14: Affected proteins in the VEGFA-VEGFR2 pathway from the Reactome database (R-HAS-194138); Table S1: Summary of putative compounds of P. frutescens leaf extracts from metabolomic studies; Table S2: First-pass relevant gene ontology (GO) terms; Table S3: Second-pass relevant gene ontology (GO) terms; Table S4: Identified protein targets from protein-protein interaction network analysis; and, Table S5: CDOCKER energy (in kcal/mol) of protein-ligand complexes.

Author Contributions

P.-W.T. and S.K.-H.H.: conceptualization, resources, project administration, and funding acquisition; P.J.B.G., S.K.-H.H. and P.-W.T.: methodology; P.J.B.G. and S.K.-H.H.: software, validation, formal analysis, investigation, data curation, writing—original draft preparation, and visualization; P.-W.T., S.K.-H.H., K.A.D.C.-C. and R.B.L.: writing—review and editing, and supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the Division of Urology, Department of Surgery, Chi Mei Medical Center (CMFHR111140), National Science and Technology Council, Taiwan (NSTC 112-2221-E-309-001), for funding this project.

Data Availability Statement

The data and results presented in this study are available on request from the first author and corresponding author.

Acknowledgments

The authors sincerely appreciate the financial support from the Office of Research and Development, Chang Jung Christian University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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]
  2. Siddiqui, E.; Mumtaz, F.H.; Gelister, J. Understanding Prostate Cancer. J. R. Soc. Promot. Health 2004, 124, 219–221. [Google Scholar] [CrossRef] [PubMed]
  3. Litwin, M.S.; Tan, H.-J. The Diagnosis and Treatment of Prostate Cancer. JAMA 2017, 317, 2532–2542. [Google Scholar] [CrossRef] [PubMed]
  4. Harvey, P.; Basuita, A.; Endersby, D.; Curtis, B.; Iacovidou, A.; Walker, M. A Systematic Review of the Diagnostic Accuracy of Prostate Specific Antigen. BMC Urol. 2009, 9, 14. [Google Scholar] [CrossRef] [PubMed]
  5. Sandhu, S.; Moore, C.M.; Chiong, E.; Beltran, H.; Bristow, R.G.; Williams, S.G. Prostate Cancer. Lancet 2021, 398, 1075–1090. [Google Scholar] [CrossRef]
  6. Teo, M.Y.; Rathkopf, D.E.; Kantoff, P. Treatment of Advanced Prostate Cancer. Annu. Rev. Med. 2019, 70, 479–499. [Google Scholar] [CrossRef] [PubMed]
  7. Basati, G.; Ghanadi, P.; Abbaszadeh, S. A Review of the Most Important Natural Antioxidants and Effective Medicinal Plants in Traditional Medicine on Prostate Cancer and Its Disorders. J. Herbmed. Pharmacol. 2020, 9, 112–120. [Google Scholar] [CrossRef]
  8. Cassileth, B.R.; Deng, G. Complementary and Alternative Therapies for Cancer. Oncologist 2004, 9, 80–89. [Google Scholar] [CrossRef]
  9. Yin, S.-Y.; Wei, W.-C.; Jian, F.-Y.; Yang, N.-S. Therapeutic Applications of Herbal Medicines for Cancer Patients. Evid.-Based Complement. Altern. Med. 2013, 2013, 302426. [Google Scholar] [CrossRef]
  10. Wang, J.; Ding, R.; Ouyang, T.; Gao, H.; Kan, H.; Li, Y.; Hu, Q.; Yang, Y. Systematic Investigation of the Mechanism of Herbal Medicines for the Treatment of Prostate Cancer. Aging 2023, 15, 1004–1024. [Google Scholar] [CrossRef]
  11. Ghosh, S.; Hazra, J.; Pal, K.; Nelson, V.K.; Pal, M. Prostate Cancer: Therapeutic Prospect with Herbal Medicine. Curr. Res. Pharmacol. Drug Discov. 2021, 2, 100034. [Google Scholar] [CrossRef] [PubMed]
  12. Lin, Y.-H.; Chen, K.-K.; Chiu, J.-H. Coprescription of Chinese Herbal Medicine and Western Medications among Prostate Cancer Patients: A Population-Based Study in Taiwan. Evid.-Based Complement. Altern. Med. 2012, 2012, 147015. [Google Scholar] [CrossRef] [PubMed]
  13. Li, C.; Huo, C.; Zhang, M.; Shi, Q. Chemistry of Chinese Yew, Taxus chinensis Var. Mairei. Biochem. Syst. Ecol. 2008, 36, 266–282. [Google Scholar] [CrossRef]
  14. Yuan, H.; Ma, Q.; Ye, L.; Piao, G. The Traditional Medicine and Modern Medicine from Natural Products. Molecules 2016, 21, 559. [Google Scholar] [CrossRef]
  15. Nie, J.; Zhao, C.; Deng, L.; Chen, J.; Yu, B.; Wu, X.; Pang, P.; Chen, X. Efficacy of Traditional Chinese Medicine in Treating Cancer. Biomed. Rep. 2016, 4, 3–14. [Google Scholar] [CrossRef]
  16. Liu, J.-M.; Lin, P.-H.; Hsu, R.-J.; Chang, Y.-H.; Cheng, K.-C.; Pang, S.-T.; Lin, S.-K. Complementary Traditional Chinese Medicine Therapy Improves Survival in Patients with Metastatic Prostate Cancer. Medicine 2016, 95, e4475. [Google Scholar] [CrossRef]
  17. POWO. Plants of the World Online. Available online: http://www.plantsoftheworldonline.org/ (accessed on 26 May 2023).
  18. Yu, H.; Qiu, J.-F.; Ma, L.-J.; Hu, Y.-J.; Li, P.; Wan, J.-B. Phytochemical and Phytopharmacological Review of Perilla frutescens L. (Labiatae), a Traditional Edible-Medicinal Herb in China. Food Chem. Toxicol. 2017, 108, 375–391. [Google Scholar] [CrossRef]
  19. Wu, X.; Dong, S.; Chen, H.; Guo, M.; Sun, Z.; Luo, H. Perilla frutescens: A Traditional Medicine and Food Homologous Plant. Chin. Herb. Med. 2023, 15, 369–375. [Google Scholar] [CrossRef]
  20. Jeong, J.H.; Park, H.J.; Chi, G.Y.; Choi, Y.H.; Park, S.H. An Ethanol Extract of Perilla frutescens Leaves Suppresses Adrenergic Agonist-Induced Metastatic Ability of Cancer Cells by Inhibiting Src-Mediated EMT. Molecules 2023, 28, 3414. [Google Scholar] [CrossRef]
  21. Lee, J.H.; Cho, H.D.; Jeong, I.Y.; Lee, M.K.; Seo, K. Il Sensitization of Tumor Necrosis Factor-Related Apoptosis-Inducing Ligand (TRAIL)-Resistant Primary Prostate Cancer Cells by Isoegomaketone from Perilla frutescens. J. Nat. Prod. 2014, 77, 2438–2443. [Google Scholar] [CrossRef]
  22. Chandran, U.; Mehendale, N.; Patil, S.; Chaguturu, R.; Patwardhan, B. Network Pharmacology. In Innovative Approaches in Drug Discovery; Elsevier: Amsterdam, The Netherlands, 2017; pp. 127–164. [Google Scholar]
  23. Zhao, L.; Zhang, H.; Li, N.; Chen, J.; Xu, H.; Wang, Y.; Liang, Q. Network Pharmacology, a Promising Approach to Reveal the Pharmacology Mechanism of Chinese Medicine Formula. J. Ethnopharmacol. 2023, 309, 116306. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, X.; Wang, Z.-Y.; Zheng, J.-H.; Li, S. TCM Network Pharmacology: A New Trend towards Combining Computational, Experimental and Clinical Approaches. Chin. J. Nat. Med. 2021, 19, 1–11. [Google Scholar] [CrossRef] [PubMed]
  25. Ahmed, H. Ethnomedicinal, Phytochemical and Pharmacological Investigations of Perilla frutescens (L.) Britt. Molecules 2018, 24, 102. [Google Scholar] [CrossRef] [PubMed]
  26. Zhou, Y.; Zheng, J.; Li, Y.; Xu, D.-P.; Li, S.; Chen, Y.-M.; Li, H.-B. Natural Polyphenols for Prevention and Treatment of Cancer. Nutrients 2016, 8, 515. [Google Scholar] [CrossRef]
  27. Ullah, A.; Munir, S.; Badshah, S.L.; Khan, N.; Ghani, L.; Poulson, B.G.; Emwas, A.-H.; Jaremko, M. Important Flavonoids and Their Role as a Therapeutic Agent. Molecules 2020, 25, 5243. [Google Scholar] [CrossRef]
  28. Kopustinskiene, D.M.; Jakstas, V.; Savickas, A.; Bernatoniene, J. Flavonoids as Anticancer Agents. Nutrients 2020, 12, 457. [Google Scholar] [CrossRef] [PubMed]
  29. Bhosale, P.B.; Ha, S.E.; Vetrivel, P.; Kim, H.H.; Kim, S.M.; Kim, G.S. Functions of Polyphenols and Its Anticancer Properties in Biomedical Research: A Narrative Review. Transl. Cancer Res. 2020, 9, 7619–7631. [Google Scholar] [CrossRef]
  30. Lee, Y.H.; Kim, B.; Kim, S.; Kim, M.S.; Kim, H.; Hwang, S.R.; Kim, K.; Lee, J.H. Characterization of Metabolite Profiles from the Leaves of Green Perilla (Perilla frutescens) by Ultra High Performance Liquid Chromatography Coupled with Electrospray Ionization Quadrupole Time-of-Flight Mass Spectrometry and Screening for Their Antioxidant Properties. J. Food Drug Anal. 2017, 25, 776–788. [Google Scholar] [CrossRef]
  31. Chen, J.; Guo, L.; Yang, G.; Yang, A.; Zheng, Y.; Wang, L. Metabolomic Profiling of Developing Perilla Leaves Reveals the Best Harvest Time. Front. Plant. Sci. 2022, 13, 989755. [Google Scholar] [CrossRef]
  32. Zhang, J.L.; Yan, R.J.; Yu, N.; Zhang, X.; Chen, D.J.; Wu, T.; Xin, J.G. A New Caffeic Acid Tetramer from the Dracocephalum moldavica L. Nat. Prod. Res. 2017, 32, 370–373. [Google Scholar] [CrossRef]
  33. Fujita, T.; Terato, K.; Nakayama, M. Two Jasmonoid Glucosides and a Phenylvaleric Acid Glucoside from Perilla frutescens. Biosci. Biotechnol. Biochem. 1996, 60, 732–735. [Google Scholar] [CrossRef]
  34. Chin, C.-H.; Chen, S.-H.; Wu, H.-H.; Ho, C.-W.; Ko, M.-T.; Lin, C.-Y. CytoHubba: Identifying Hub Objects and Sub-Networks from Complex Interactome. BMC Syst. Biol. 2014, 8 (Suppl. S4), S11. [Google Scholar] [CrossRef] [PubMed]
  35. Tomassi, S.; D’Amore, V.M.; Di Leva, F.S.; Vannini, A.; Quilici, G.; Weinmüller, M.; Reichart, F.; Amato, J.; Romano, B.; Izzo, A.A.; et al. Halting the Spread of Herpes Simplex Virus-1: The Discovery of an Effective Dual αvβ6/αvβ8 Integrin Ligand. J. Med. Chem. 2021, 64, 6972–6984. [Google Scholar] [CrossRef] [PubMed]
  36. Ferreira de Freitas, R.; Schapira, M. A Systematic Analysis of Atomic Protein–Ligand Interactions in the PDB. Med. Chem. Commun. 2017, 8, 1970–1981. [Google Scholar] [CrossRef] [PubMed]
  37. Velu, G.; Palanichamy, V.; Rajan, A.P. Phytochemical and pharmacological importance of plant secondary metabolites in modern medicine. In Bioorganic Phase in Natural Food: An Overview; Springer International Publishing: Cham, Switzerland, 2018; pp. 135–156. [Google Scholar]
  38. Cullen, J.M.; Breen, M. An Overview of molecular cancer pathogenesis, prognosis, and diagnosis. In Tumors in Domestic Animals; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2016; pp. 1–26. [Google Scholar]
  39. Novikov, N.M.; Zolotaryova, S.Y.; Gautreau, A.M.; Denisov, E.V. Mutational Drivers of Cancer Cell Migration and Invasion. Br. J. Cancer 2021, 124, 102–114. [Google Scholar] [CrossRef] [PubMed]
  40. Testa, U.; Castelli, G.; Pelosi, E. Cellular and Molecular Mechanisms Underlying Prostate Cancer Development: Therapeutic Implications. Medicines 2019, 6, 82. [Google Scholar] [CrossRef]
  41. Moya, L.; Walpole, C.; Rae, F.; Srinivasan, S.; Seim, I.; Lai, J.; Nicol, D.; Williams, E.D.; Clements, J.A.; Batra, J. Characterisation of Cell Lines Derived from Prostate Cancer Patients with Localised Disease. Prostate Cancer Prostatic Dis. 2023. [Google Scholar] [CrossRef]
  42. Saranyutanon, S.; Deshmukh, S.K.; Dasgupta, S.; Pai, S.; Singh, S.; Singh, A.P. Cellular and Molecular Progression of Prostate Cancer: Models for Basic and Preclinical Research. Cancers 2020, 12, 2651. [Google Scholar] [CrossRef]
  43. Cunningham, D.; You, Z. In Vitro and in Vivo Model Systems Used in Prostate Cancer Research. J. Biol. Methods 2015, 2, e17. [Google Scholar] [CrossRef]
  44. Alimirah, F.; Chen, J.; Basrawala, Z.; Xin, H.; Choubey, D. DU-145 and PC-3 Human Prostate Cancer Cell Lines Express Androgen Receptor: Implications for the Androgen Receptor Functions and Regulation. FEBS Lett. 2006, 580, 2294–2300. [Google Scholar] [CrossRef]
  45. Foster, K.; Oyenihi, O.; Rademan, S.; Erhabor, J.; Matsabisa, M.; Barker, J.; Langat, M.K.; Kendal-Smith, A.; Asemota, H.; Delgoda, R. Selective Cytotoxic and Anti-Metastatic Activity in DU-145 Prostate Cancer Cells Induced by Annona muricata L. Bark Extract and Phytochemical, Annonacin. BMC Complement. Med. Ther. 2020, 20, 1–15. [Google Scholar] [CrossRef]
  46. Fithrotunnisa, Q.; Arsianti, A.; Kurniawan, G.; Qorina, F.; Tejaputri, N.A.; Azizah, N.N. In Vitro Cytotoxicity of Hibiscus Sabdariffa Linn Extracts on A549 Lung Cancer Cell Line. Pharmacogn. J. 2020, 12, 14–19. [Google Scholar] [CrossRef]
  47. Nordin, M.L.; Abdul Kadir, A.; Zakaria, Z.A.; Abdullah, R.; Abdullah, M.N.H. In Vitro Investigation of Cytotoxic and Antioxidative Activities of Ardisia crispa against Breast Cancer Cell Lines, MCF-7 and MDA-MB-231. BMC Complement Altern. Med. 2018, 18, 87. [Google Scholar] [CrossRef] [PubMed]
  48. Salehi, B.; Venditti, A.; Sharifi-Rad, M.; Kręgiel, D.; Sharifi-Rad, J.; Durazzo, A.; Lucarini, M.; Santini, A.; Souto, E.; Novellino, E.; et al. The Therapeutic Potential of Apigenin. Int. J. Mol. Sci. 2019, 20, 1305. [Google Scholar] [CrossRef]
  49. Ashrafizadeh, M.; Bakhoda, M.R.; Bahmanpour, Z.; Ilkhani, K.; Zarrabi, A.; Makvandi, P.; Khan, H.; Mazaheri, S.; Darvish, M.; Mirzaei, H. Apigenin as Tumor Suppressor in Cancers: Biotherapeutic Activity, Nanodelivery, and Mechanisms with Emphasis on Pancreatic Cancer. Front. Chem. 2020, 8, 829. [Google Scholar] [CrossRef] [PubMed]
  50. Tehami, W.; Nani, A.; Khan, N.A.; Hichami, A. New Insights into the Anticancer Effects of P-Coumaric Acid: Focus on Colorectal Cancer. Dose Response 2023, 21, 155932582211507. [Google Scholar] [CrossRef]
  51. Gutiérrez Mercado, Y.K.; Mateos Díaz, J.C.; Ojeda Hernández, D.D.; López Gonzalez, F.J.; Reza Zaldivar, E.E.; Hernández Sapiens, M.A.; Gómez Pinedo, U.A.; Estrada, R.S.; Macías Carballo, M.; Canales Aguirre, A.A. Ortho-Coumaric Acid Derivatives with Therapeutic Potential in a Three-Dimensional Culture of the Immortalised U-138 MG Glioblastoma Multiforme Cell Line. Neurol. Perspect. 2022, 2, S19–S30. [Google Scholar] [CrossRef]
  52. Lin, Z.; Huang, S.; LingHu, X.; Wang, Y.; Wang, B.; Zhong, S.; Xie, S.; Xu, X.; Yu, A.; Nagai, A.; et al. Perillaldehyde Inhibits Bone Metastasis and Receptor Activator of Nuclear Factor-ΚB Ligand (RANKL) Signaling-Induced Osteoclastogenesis in Prostate Cancer Cell Lines. Bioengineered 2022, 13, 2710–2719. [Google Scholar] [CrossRef]
  53. Cao, P.; Liu, B.; Du, F.; Li, D.; Wang, Y.; Yan, X.; Li, X.; Li, Y. Scutellarin Suppresses Proliferation and Promotes Apoptosis in A549 Lung Adenocarcinoma Cells via AKT/MTOR/4EBP1 and STAT3 Pathways. Thorac. Cancer 2019, 10, 492–500. [Google Scholar] [CrossRef]
  54. Gao, C.; Zhou, Y.; Jiang, Z.; Zhao, Y.; Zhang, D.; Cong, X.; Cao, R.; Li, H.; Tian, W. Cytotoxic and Chemosensitization Effects of Scutellarin from Traditional Chinese Herb Scutellaria altissima L. in Human Prostate Cancer Cells. Oncol. Rep. 2017, 38, 1491–1499. [Google Scholar] [CrossRef]
  55. Iqbal, J.; Abbasi, B.A.; Mahmood, T.; Kanwal, S.; Ali, B.; Shah, S.A.; Khalil, A.T. Plant-Derived Anticancer Agents: A Green Anticancer Approach. Asian Pac. J. Trop. Biomed. 2017, 7, 1129–1150. [Google Scholar] [CrossRef]
  56. Vesaghhamedani, S.; Mazloumi Kiapey, S.S.; Gowhari Shabgah, A.; Amiresmaili, S.; Jahanara, A.; Oveisee, M.; Shekarchi, A.; Gheibihayat, S.M.; Jadidi-Niaragh, F.; Gholizadeh Navashenaq, J. From Traditional Medicine to Modern Oncology: Scutellarin, a Promising Natural Compound in Cancer Treatment. Prog. Biophys. Mol. Biol. 2023, 180–181, 19–27. [Google Scholar] [CrossRef] [PubMed]
  57. Yesil-Celiktas, O.; Sevimli, C.; Bedir, E.; Vardar-Sukan, F. Inhibitory Effects of Rosemary Extracts, Carnosic Acid and Rosmarinic Acid on the Growth of Various Human Cancer Cell Lines. Plant. Foods Hum. Nutr. 2010, 65, 158–163. [Google Scholar] [CrossRef] [PubMed]
  58. Moore, J.; Yousef, M.; Tsiani, E. Anticancer Effects of Rosemary (Rosmarinus officinalis L.) Extract and Rosemary Extract Polyphenols. Nutrients 2016, 8, 731. [Google Scholar] [CrossRef]
  59. Anwar, S.; Shamsi, A.; Shahbaaz, M.; Queen, A.; Khan, P.; Hasan, G.M.; Islam, A.; Alajmi, M.F.; Hussain, A.; Ahmad, F.; et al. Rosmarinic Acid Exhibits Anticancer Effects via MARK4 Inhibition. Sci. Rep. 2020, 10, 10300. [Google Scholar] [CrossRef] [PubMed]
  60. Jaglanian, A.; Termini, D.; Tsiani, E. Rosemary (Rosmarinus officinalis L.) Extract Inhibits Prostate Cancer Cell Proliferation and Survival by Targeting Akt and MTOR. Biomed. Pharmacother. 2020, 131, 110717. [Google Scholar] [CrossRef]
  61. Park, H.U.; Suy, S.; Danner, M.; Dailey, V.; Zhang, Y.; Li, H.; Hyduke, D.R.; Collins, B.T.; Gagnon, G.; Kallakury, B.; et al. AMP-Activated Protein Kinase Promotes Human Prostate Cancer Cell Growth and Survival. Mol. Cancer Ther. 2009, 8, 733–741. [Google Scholar] [CrossRef]
  62. Chau, V.; Madan, R.A.; Aragon-Ching, J.B. Protein Kinase Inhibitors for the Treatment of Prostate Cancer. Expert Opin. Pharm. 2021, 22, 1889–1899. [Google Scholar] [CrossRef]
  63. Kim, J.-H.; Xu, C.; Keum, Y.-S.; Reddy, B.; Conney, A.; Kong, A.-N.T. Inhibition of EGFR Signaling in Human Prostate Cancer PC-3 Cells by Combination Treatment with β-Phenylethyl Isothiocyanate and Curcumin. Carcinogenesis 2006, 27, 475–482. [Google Scholar] [CrossRef]
  64. Jones, H.E.; Dutkowski, C.M.; Barrow, D.; Harper, M.E.; Wakeling, A.E.; Nicholson, R.I. New EGF-R Selective Tyrosine Kinase Inhibitor Reveals Variable Growth Responses in Prostate Carcinoma Cell Lines PC-3 and DU-145. Int. J. Cancer 1997, 71, 1010–1018. [Google Scholar] [CrossRef]
  65. Figueira, M.I.; Cardoso, H.J.; Correia, S.; Maia, C.J.; Socorro, S. The Stem Cell Factor (SCF)/c-KIT System in Carcinogenesis of Reproductive Tissues: What Does the Hormonal Regulation Tell Us? Cancer Lett. 2017, 405, 10–21. [Google Scholar] [CrossRef] [PubMed]
  66. Pungsrinont, T.; Kallenbach, J.; Baniahmad, A. Role of PI3K-AKT-MTOR Pathway as a Pro-Survival Signaling and Resistance-Mediating Mechanism to Therapy of Prostate Cancer. Int. J. Mol. Sci. 2021, 22, 11088. [Google Scholar] [CrossRef] [PubMed]
  67. Kakarla, M.; ChallaSivaKanaka, S.; Dufficy, M.F.; Gil, V.; Filipovich, Y.; Vickman, R.; Crawford, S.E.; Hayward, S.W.; Franco, O.E. Ephrin B Activate Src Family Kinases in Fibroblasts Inducing Stromal Remodeling in Prostate Cancer. Cancers 2022, 14, 2336. [Google Scholar] [CrossRef] [PubMed]
  68. Roberts, E.; Cossigny, D.A.F.; Quan, G.M.Y. The Role of Vascular Endothelial Growth Factor in Metastatic Prostate Cancer to the Skeleton. Prostate Cancer 2013, 2013, 1–8. [Google Scholar] [CrossRef]
  69. Alberto, M.; Yim, A.; Lawrentschuk, N.; Bolton, D. Dysfunctional Lipid Metabolism—The Basis for How Genetic Abnormalities Express the Phenotype of Aggressive Prostate Cancer. Cancers 2023, 15, 341. [Google Scholar] [CrossRef]
  70. Gururajan, M.; Cavassani, K.A.; Sievert, M.; Duan, P.; Lichterman, J.; Huang, J.-M.; Smith, B.; You, S.; Nandana, S.; Chu, G.C.-Y.; et al. SRC Family Kinase FYN Promotes the Neuroendocrine Phenotype and Visceral Metastasis in Advanced Prostate Cancer. Oncotarget 2015, 6, 44072–44083. [Google Scholar] [CrossRef]
  71. Saito, Y.D.; Jensen, A.R.; Salgia, R.; Posadas, E.M. Fyn: A Novel Molecular Target in Prostate Cancer. Cancer 2010, 116, 1629–1637. [Google Scholar] [CrossRef]
  72. Peng, S.; Fu, Y. FYN: Emerging Biological Roles and Potential Therapeutic Targets in Cancer. J. Transl. Med. 2023, 21, 84. [Google Scholar] [CrossRef]
  73. Jelić, D.; Mildner, B.; Koštrun, S.; Nujić, K.; Verbanac, D.; Čulić, O.; Antolović, R.; Brandt, W. Homology Modeling of Human Fyn Kinase Structure: Discovery of Rosmarinic Acid as a New Fyn Kinase Inhibitor and in Silico Study of Its Possible Binding Modes. J. Med. Chem. 2007, 50, 1090–1100. [Google Scholar] [CrossRef]
  74. Kang, N.J.; Lee, K.W.; Shin, B.J.; Jung, S.K.; Hwang, M.K.; Bode, A.M.; Heo, Y.-S.; Lee, H.J.; Dong, Z. Caffeic Acid, a Phenolic Phytochemical in Coffee, Directly Inhibits Fyn Kinase Activity and UVB-Induced COX-2 Expression. Carcinogenesis 2008, 30, 321–330. [Google Scholar] [CrossRef]
  75. Liu, Y.; Xu, X.; Tang, H.; Pan, Y.; Hu, B.; Huang, G. Rosmarinic Acid Inhibits Cell Proliferation, Migration, and Invasion and Induces Apoptosis in Human Glioma Cells. Int. J. Mol. Med. 2021, 47, 67. [Google Scholar] [CrossRef]
  76. Quaglia, F.; Krishn, S.R.; Wang, Y.; Goodrich, D.W.; McCue, P.; Kossenkov, A.V.; Mandigo, A.C.; Knudsen, K.E.; Weinreb, P.H.; Corey, E.; et al. Differential Expression of αVβ3 and αVβ6 Integrins in Prostate Cancer Progression. PLoS ONE 2021, 16, e0244985. [Google Scholar] [CrossRef] [PubMed]
  77. Pang, X.; He, X.; Qiu, Z.; Zhang, H.; Xie, R.; Liu, Z.; Gu, Y.; Zhao, N.; Xiang, Q.; Cui, Y. Targeting Integrin Pathways: Mechanisms and Advances in Therapy. Signal. Transduct. Target. Ther. 2023, 8, 1. [Google Scholar] [CrossRef] [PubMed]
  78. Liu, Z.; Wang, F.; Chen, X. Integrin αvβ3-Targeted Cancer Therapy. Drug Dev. Res. 2008, 69, 329–339. [Google Scholar] [CrossRef] [PubMed]
  79. Krishn, S.R.; Singh, A.; Bowler, N.; Duffy, A.N.; Friedman, A.; Fedele, C.; Kurtoglu, S.; Tripathi, S.K.; Wang, K.; Hawkins, A.; et al. Prostate Cancer Sheds the αvβ3 Integrin in Vivo through Exosomes. Matrix Biol. 2019, 77, 41–57. [Google Scholar] [CrossRef] [PubMed]
  80. Kim, E.-Y.; Choi, H.-J.; Chung, T.-W.; Choi, J.-Y.; Kim, H.S.; Jung, Y.-S.; Lee, S.-O.; Ha, K.-T. Water-Extracted Perilla frutescens Increases Endometrial Receptivity Though Leukemia Inhibitory Factor-Dependent Expression of Integrins. J. Pharm. Sci. 2016, 131, 259–266. [Google Scholar] [CrossRef]
  81. Munkley, J.; Livermore, K.E.; McClurg, U.L.; Kalna, G.; Knight, B.; McCullagh, P.; McGrath, J.; Crundwell, M.; Leung, H.Y.; Robson, C.N.; et al. The PI3K Regulatory Subunit Gene PIK3R1 Is under Direct Control of Androgens and Repressed in Prostate Cancer Cells. Oncoscience 2015, 2, 755–764. [Google Scholar] [CrossRef]
  82. Dornan, G.L.; Burke, J.E. Molecular Mechanisms of Human Disease Mediated by Oncogenic and Primary Immunodeficiency Mutations in Class IA Phosphoinositide 3-Kinases. Front. Immunol. 2018, 9, 575. [Google Scholar] [CrossRef]
  83. Chakraborty, G.; Nandakumar, S.; Hirani, R.; Nguyen, B.; Stopsack, K.H.; Kreitzer, C.; Rajanala, S.H.; Ghale, R.; Mazzu, Y.Z.; Pillarsetty, N.V.K.; et al. The Impact of PIK3R1 Mutations and Insulin–PI3K–Glycolytic Pathway Regulation in Prostate Cancer. Clin. Cancer Res. 2022, 28, 3603–3617. [Google Scholar] [CrossRef]
  84. Zhao, J.; Xu, L.; Jin, D.; Xin, Y.; Tian, L.; Wang, T.; Zhao, D.; Wang, Z.; Wang, J. Rosmarinic Acid and Related Dietary Supplements: Potential Applications in the Prevention and Treatment of Cancer. Biomolecules 2022, 12, 1410. [Google Scholar] [CrossRef]
  85. Wang, L.; Yang, H.; Wang, C.; Shi, X.; Li, K. Rosmarinic Acid Inhibits Proliferation and Invasion of Hepatocellular Carcinoma Cells SMMC 7721 via PI3K/AKT/MTOR Signal Pathway. Biomed. Pharmacother. 2019, 120, 109443. [Google Scholar] [CrossRef] [PubMed]
  86. Heldin, C.-H. Targeting the PDGF Signaling Pathway in Tumor Treatment. Cell Commun. Signal. 2013, 11, 97. [Google Scholar] [CrossRef] [PubMed]
  87. Heldin, C.H.; Ernlund, A.; Rorsman, C.; Rönnstrand, L. Dimerization of B-Type Platelet-Derived Growth Factor Receptors Occurs After Ligand Binding and Is Closely Associated with Receptor Kinase Activation. J. Biol. Chem. 1989, 264, 8905–8912. [Google Scholar] [CrossRef] [PubMed]
  88. Herren, B.; Rooney, B.; Weyer, K.A.; Iberg, N.; Schmid, G.; Pech, M. Dimerization of Extracellular Domains of Platelet-Derived Growth Factor Receptors. A Revised Model of Receptor-Ligand Interaction. J. Biol. Chem. 1993, 268, 15088–15095. [Google Scholar] [CrossRef]
  89. Russell, M.R.; Liu, Q.; Lei, H.; Kazlauskas, A.; Fatatis, A. The α-Receptor for Platelet-Derived Growth Factor Confers Bone-Metastatic Potential to Prostate Cancer Cells by Ligand- and Dimerization-Independent Mechanisms. Cancer Res. 2010, 70, 4195–4203. [Google Scholar] [CrossRef] [PubMed]
  90. Liu, Q.; Jernigan, D.; Zhang, Y.; Fatatis, A. Implication of Platelet-Derived Growth Factor Receptor Alpha in Prostate Cancer Skeletal Metastasis. Chin. J. Cancer 2011, 30, 612. [Google Scholar] [CrossRef]
  91. Hofer, M.D.; Fecko, A.; Shen, R.; Setlur, S.R.; Pienta, K.G.; Tomlins, S.A.; Chinnaiyan, A.M.; Rubin, M.A. Expression of the Platelet-Derived Growth Factor Receptor in Prostate Cancer and Treatment Implications with Tyrosine Kinase Inhibitors. Neoplasia 2004, 6, 503–512. [Google Scholar] [CrossRef]
  92. Wikström, P.; Bergström, S.H.; Josefsson, A.; Semenas, J.; Nordstrand, A.; Thysell, E.; Crnalic, S.; Widmark, A.; Karlsson, C.T.; Bergh, A. Epithelial and Stromal Characteristics of Primary Tumors Predict the Bone Metastatic Subtype of Prostate Cancer and Patient Survival after Androgen-Deprivation Therapy. Cancers 2022, 14, 5195. [Google Scholar] [CrossRef]
  93. Hägglöf, C.; Hammarsten, P.; Josefsson, A.; Stattin, P.; Paulsson, J.; Bergh, A.; Östman, A. Stromal PDGFRβ Expression in Prostate Tumors and Non-Malignant Prostate Tissue Predicts Prostate Cancer Survival. PLoS ONE 2010, 5, e10747. [Google Scholar] [CrossRef]
  94. Makino, T.; Ono, T.; Muso, E.; Yoshida, H.; Honda, G.; Sasayama, S. Inhibitory Effects of Rosmarinic Acid on the Proliferation of Cultured Murine Mesangial Cells. Nephrol. Dial. Transplant. 2000, 15, 1140–1145. [Google Scholar] [CrossRef]
  95. Tsai, P.-W.; Mailem, R.C.; Tayo, L.L.; Hsueh, C.-C.; Tseng, C.-C.; Chen, B.-Y. Interactive Network Pharmacology and Electrochemical Analysis Reveals Electron Transport-Mediating Characteristics of Chinese Medicine Formula Jing Guan Fang. J. Taiwan Inst. Chem. Eng. 2023, 147, 104898. [Google Scholar] [CrossRef]
  96. Park, J.S.; Kwon, J.K.; Kim, H.R.; Kim, H.J.; Kim, B.S.; Jung, J.Y. Farnesol Induces Apoptosis of DU145 Prostate Cancer Cells through the PI3K/Akt and MAPK Pathways. Int. J. Mol. Med. 2014, 33, 1169–1176. [Google Scholar] [CrossRef]
  97. WST-1 Assay Reagent—Cell Proliferation (Ready to Use) (Ab155902). Available online: https://www.abcam.com/products/assay-kits/wst-1-assay-reagent-cell-proliferation-ready-to-use-ab155902.html (accessed on 3 May 2023).
  98. Jin, J.; Chen, B.; Zhan, X.; Zhou, Z.; Liu, H.; Dong, Y. Network Pharmacology and Molecular Docking Study on the Mechanism of Colorectal Cancer Treatment Using Xiao-Chai-Hu-Tang. PLoS ONE 2021, 16, e0252508. [Google Scholar] [CrossRef]
  99. Peng, J.; Zhang, K.; Wang, L.; Peng, F.; Zhang, C.; Long, K.; Chen, J.; Zhou, X.; Gao, P.; Fan, G. Integrating Network Pharmacology and Molecular Docking to Explore the Potential Mechanism of Xinguan No. 3 in the Treatment of COVID-19. Open Chem. 2022, 20, 570–582. [Google Scholar] [CrossRef]
  100. Li, S. Network Pharmacology Evaluation Method Guidance—Draft. World J. Tradit Chin. Med. 2021, 7, 148. [Google Scholar] [CrossRef]
  101. Nickel, J.; Gohlke, B.-O.; Erehman, J.; Banerjee, P.; Rong, W.W.; Goede, A.; Dunkel, M.; Preissner, R. SuperPred: Update on Drug Classification and Target Prediction. Nucleic Acids Res. 2014, 42, W26–W31. [Google Scholar] [CrossRef]
  102. Tang, Z.; Kang, B.; Li, C.; Chen, T.; Zhang, Z. GEPIA2: An Enhanced Web Server for Large-Scale Expression Profiling and Interactive Analysis. Nucleic Acids Res. 2019, 47, W556–W560. [Google Scholar] [CrossRef]
  103. Koopmans, F.; van Nierop, P.; Andres-Alonso, M.; Byrnes, A.; Cijsouw, T.; Coba, M.P.; Cornelisse, L.N.; Farrell, R.J.; Goldschmidt, H.L.; Howrigan, D.P.; et al. SynGO: An Evidence-Based, Expert-Curated Knowledge Base for the Synapse. Neuron 2019, 103, 217–234.e4. [Google Scholar] [CrossRef]
  104. Khan, A.; Mathelier, A. Intervene: A Tool for Intersection and Visualization of Multiple Gene or Genomic Region Sets. BMC Bioinform. 2017, 18, 287. [Google Scholar] [CrossRef]
  105. Oliveros, J.C. Venny. An Interactive Tool for Comparing Lists with Venn’s Diagrams. Available online: https://bioinfogp.cnb.csic.es/tools/venny/index.html (accessed on 1 June 2023).
  106. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING V11: Protein-Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-Wide Experimental Datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef]
  107. Snel, B.; Lehmann, G.; Bork, P.; Huynen, M.A. STRING: A Web-Server to Retrieve and Display the Repeatedly Occurring Neighbourhood of a Gene. Nucleic Acids Res. 2000, 28, 3442–3444. [Google Scholar] [CrossRef] [PubMed]
  108. Von Mering, C.; Huynen, M.; Jaeggi, D.; Schmidt, S.; Bork, P.; Snel, B. STRING: A Database of Predicted Functional Associations between Proteins. Nucleic Acids Res. 2003, 31, 258–261. [Google Scholar] [CrossRef] [PubMed]
  109. Von Mering, C.; Jensen, L.J.; Snel, B.; Hooper, S.D.; Krupp, M.; Foglierini, M.; Jouffre, N.; Huynen, M.A.; Bork, P. STRING: Known and Predicted Protein-Protein Associations, Integrated and Transferred across Organisms. Nucleic Acids Res. 2005, 33, D433–D437. [Google Scholar] [CrossRef] [PubMed]
  110. Von Mering, C.; Jensen, L.J.; Kuhn, M.; Chaffron, S.; Doerks, T.; Krüger, B.; Snel, B.; Bork, P. STRING 7—Recent Developments in the Integration and Prediction of Protein Interactions. Nucleic Acids Res. 2007, 35, D358–D362. [Google Scholar] [CrossRef]
  111. Jensen, L.J.; Kuhn, M.; Stark, M.; Chaffron, S.; Creevey, C.; Muller, J.; Doerks, T.; Julien, P.; Roth, A.; Simonovic, M.; et al. STRING 8—A Global View on Proteins and Their Functional Interactions in 630 Organisms. Nucleic Acids Res. 2009, 37, D412–D416. [Google Scholar] [CrossRef]
  112. Szklarczyk, D.; Franceschini, A.; Kuhn, M.; Simonovic, M.; Roth, A.; Minguez, P.; Doerks, T.; Stark, M.; Muller, J.; Bork, P.; et al. The STRING Database in 2011: Functional Interaction Networks of Proteins, Globally Integrated and Scored. Nucleic Acids Res. 2011, 39, D561–D568. [Google Scholar] [CrossRef]
  113. Franceschini, A.; Szklarczyk, D.; Frankild, S.; Kuhn, M.; Simonovic, M.; Roth, A.; Lin, J.; Minguez, P.; Bork, P.; Von Mering, C.; et al. STRING v9.1: Protein-Protein Interaction Networks, with Increased Coverage and Integration. Nucleic Acids Res. 2013, 41, D808–D815. [Google Scholar] [CrossRef]
  114. Franceschini, A.; Lin, J.; Von Mering, C.; Jensen, L.J. SVD-Phy: Improved Prediction of Protein Functional Associations through Singular Value Decomposition of Phylogenetic Profiles. Bioinformatics 2016, 32, 1085–1087. [Google Scholar] [CrossRef]
  115. Szklarczyk, D.; Franceschini, A.; Wyder, S.; Forslund, K.; Heller, D.; Huerta-Cepas, J.; Simonovic, M.; Roth, A.; Santos, A.; Tsafou, K.P.; et al. STRING V10: Protein-Protein Interaction Networks, Integrated over the Tree of Life. Nucleic Acids Res. 2015, 43, D447–D452. [Google Scholar] [CrossRef]
  116. Szklarczyk, D.; Morris, J.H.; Cook, H.; Kuhn, M.; Wyder, S.; Simonovic, M.; Santos, A.; Doncheva, N.T.; Roth, A.; Bork, P.; et al. The STRING Database in 2017: Quality-Controlled Protein-Protein Association Networks, Made Broadly Accessible. Nucleic Acids Res. 2017, 45, D362–D368. [Google Scholar] [CrossRef]
  117. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING Database in 2021: Customizable Protein-Protein Networks, and Functional Characterization of User-Uploaded Gene/Measurement Sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef] [PubMed]
  118. Szklarczyk, D.; Kirsch, R.; Koutrouli, M.; Nastou, K.; Mehryary, F.; Hachilif, R.; Gable, A.L.; Fang, T.; Doncheva, N.T.; Pyysalo, S.; et al. The STRING Database in 2023: Protein-Protein Association Networks and Functional Enrichment Analyses for Any Sequenced Genome of Interest. Nucleic Acids Res. 2023, 51, D638–D646. [Google Scholar] [CrossRef]
  119. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  120. Kanehisa, M.; Furumichi, M.; Sato, Y.; Ishiguro-Watanabe, M.; Tanabe, M. KEGG: Integrating Viruses and Cellular Organisms. Nucleic Acids Res. 2021, 49, D545–D551. [Google Scholar] [CrossRef] [PubMed]
  121. Raudvere, U.; Kolberg, L.; Kuzmin, I.; Arak, T.; Adler, P.; Peterson, H.; Vilo, J. G:Profiler: A Web Server for Functional Enrichment Analysis and Conversions of Gene Lists (2019 Update). Nucleic Acids Res. 2019, 47, W191–W198. [Google Scholar] [CrossRef]
  122. Aleksander, S.A.; Balhoff, J.; Carbon, S.; Cherry, J.M.; Drabkin, H.J.; Ebert, D.; Feuermann, M.; Gaudet, P.; Harris, N.L.; Hill, D.P.; et al. The Gene Ontology Knowledgebase in 2023. Genetics 2023, 224, iyad031. [Google Scholar] [CrossRef]
  123. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene Ontology: Tool for the Unification of Biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef]
  124. Kanehisa, M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  125. Kanehisa, M. Toward Understanding the Origin and Evolution of Cellular Organisms. Protein Sci. 2019, 28, 1947–1951. [Google Scholar] [CrossRef]
  126. Kanehisa, M.; Furumichi, M.; Sato, Y.; Kawashima, M.; Ishiguro-Watanabe, M. KEGG for Taxonomy-Based Analysis of Pathways and Genomes. Nucleic Acids Res. 2023, 51, D587–D592. [Google Scholar] [CrossRef]
  127. Gillespie, M.; Jassal, B.; Stephan, R.; Milacic, M.; Rothfels, K.; Senff-Ribeiro, A.; Griss, J.; Sevilla, C.; Matthews, L.; Gong, C.; et al. The Reactome Pathway Knowledgebase 2022. Nucleic Acids Res. 2022, 50, D687–D692. [Google Scholar] [CrossRef]
  128. Griss, J.; Viteri, G.; Sidiropoulos, K.; Nguyen, V.; Fabregat, A.; Hermjakob, H. ReactomeGSA—Efficient Multi-Omics Comparative Pathway Analysis. Mol. Cell. Proteom. 2020, 19, 2115–2125. [Google Scholar] [CrossRef] [PubMed]
  129. Jassal, B.; Matthews, L.; Viteri, G.; Gong, C.; Lorente, P.; Fabregat, A.; Sidiropoulos, K.; Cook, J.; Gillespie, M.; Haw, R.; et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res. 2020, 48, D498–D503. [Google Scholar] [CrossRef] [PubMed]
  130. Fabregat, A.; Korninger, F.; Viteri, G.; Sidiropoulos, K.; Marin-Garcia, P.; Ping, P.; Wu, G.; Stein, L.; D’Eustachio, P.; Hermjakob, H. Reactome Graph Database: Efficient Access to Complex Pathway Data. PLoS Comput. Biol. 2018, 14, e1005968. [Google Scholar] [CrossRef] [PubMed]
  131. Fabregat, A.; Sidiropoulos, K.; Viteri, G.; Marin-Garcia, P.; Ping, P.; Stein, L.; D’Eustachio, P.; Hermjakob, H. Reactome Diagram Viewer: Data Structures and Strategies to Boost Performance. Bioinformatics 2018, 34, 1208–1214. [Google Scholar] [CrossRef]
  132. Sidiropoulos, K.; Viteri, G.; Sevilla, C.; Jupe, S.; Webber, M.; Orlic-Milacic, M.; Jassal, B.; May, B.; Shamovsky, V.; Duenas, C.; et al. Reactome Enhanced Pathway Visualization. Bioinformatics 2017, 33, 3461–3467. [Google Scholar] [CrossRef] [PubMed]
  133. Fabregat, A.; Sidiropoulos, K.; Viteri, G.; Forner, O.; Marin-Garcia, P.; Arnau, V.; D’Eustachio, P.; Stein, L.; Hermjakob, H. Reactome Pathway Analysis: A High-Performance in-Memory Approach. BMC Bioinform. 2017, 18, 142. [Google Scholar] [CrossRef]
  134. Wu, G.; Haw, R. Functional Interaction Network Construction and Analysis for Disease Discovery. In Protein Bioinformatics; Springer: Berlin/Heidelberg, Germany, 2017; pp. 235–253. [Google Scholar]
  135. Venkatachalam, C.M.; Jiang, X.; Oldfield, T.; Waldman, M. LigandFit: A Novel Method for the Shape-Directed Rapid Docking of Ligands to Protein Active Sites. J. Mol. Graph. Model. 2003, 21, 289–307. [Google Scholar] [CrossRef] [PubMed]
  136. Wu, G.; Robertson, D.H.; Brooks, C.L.; Vieth, M. Detailed Analysis of Grid-Based Molecular Docking: A Case Study of CDOCKER—A CHARMm-Based MD Docking Algorithm. J. Comput. Chem. 2003, 24, 1549–1562. [Google Scholar] [CrossRef]
  137. Schrödinger LLC. The PyMOL Molecular Graphics System, version 2.5.4 2015; Schrödinger LLC: New York, NY, USA, 2015. [Google Scholar]
Plants 12 03006 g001aPlants 12 03006 g001b
Figure 1. Concentration–cytotoxicity curves of PF parts crude extracts, i.e., ethanol extracts of (a) seeds, (b) leaves, (c) seeds, and water extracts of (d) stems, (e) leaves, and (f) seeds. The negative (−) control for ethanol extract contains the media and added DMSO, whereas water extract only includes the media. Comparisons with the negative control ** p < 0.01, *** p < 0.001, and ns are not significant according to Dunnett’s pairwise comparison test.
Figure 1. Concentration–cytotoxicity curves of PF parts crude extracts, i.e., ethanol extracts of (a) seeds, (b) leaves, (c) seeds, and water extracts of (d) stems, (e) leaves, and (f) seeds. The negative (−) control for ethanol extract contains the media and added DMSO, whereas water extract only includes the media. Comparisons with the negative control ** p < 0.01, *** p < 0.001, and ns are not significant according to Dunnett’s pairwise comparison test.
Plants 12 03006 g001aPlants 12 03006 g001b
Plants 12 03006 g002
Figure 2. HPLC chromatogram of (a) 500 µg/Ml rosmarinic acid standard, (b) ethanol extract, and (c) water extract of PF-L.
Figure 2. HPLC chromatogram of (a) 500 µg/Ml rosmarinic acid standard, (b) ethanol extract, and (c) water extract of PF-L.
Plants 12 03006 g002
Plants 12 03006 g003
Figure 3. (a) Intersected genes set from predicted targets of PF-L metabolites and expressed/overexpressed genes of PRAD compared to normal prostate cells where DR are downregulated genes and UR are upregulated genes; (b) key metabolites (including RA) in the MTN of PF-L against PRAD.
Figure 3. (a) Intersected genes set from predicted targets of PF-L metabolites and expressed/overexpressed genes of PRAD compared to normal prostate cells where DR are downregulated genes and UR are upregulated genes; (b) key metabolites (including RA) in the MTN of PF-L against PRAD.
Plants 12 03006 g003
Plants 12 03006 g004
Figure 4. PPI network of intersected genes set with interaction score of at least 0.90 where red-labeled proteins are top-ranking in interactivity through MCC and DMNC analysis.
Figure 4. PPI network of intersected genes set with interaction score of at least 0.90 where red-labeled proteins are top-ranking in interactivity through MCC and DMNC analysis.
Plants 12 03006 g004
Plants 12 03006 g005
Figure 5. CDOCKER energy of protein–ligand interaction (in kcal/mol) where lower energy (darker color) indicates good binding affinity. Controls’ abbreviations: 5FU, 5-fluorouracil; RGD, [RGD-ChgE]-CONH2; ALP, alpelisib; NIN, nintedanib; PON, ponatinib; SUN, sunitinib.
Figure 5. CDOCKER energy of protein–ligand interaction (in kcal/mol) where lower energy (darker color) indicates good binding affinity. Controls’ abbreviations: 5FU, 5-fluorouracil; RGD, [RGD-ChgE]-CONH2; ALP, alpelisib; NIN, nintedanib; PON, ponatinib; SUN, sunitinib.
Plants 12 03006 g005
Plants 12 03006 g006
Figure 6. (a–e) Protein–ligand binding site and (a′–e′) intermolecular interactions of docked RA (28) with (a,a′) FYN, (b,b′) PDGFRB, (c,c′) PIK3R1, (d,d′) ITGB3, and (e,e′) PDGFRA.
Figure 6. (a–e) Protein–ligand binding site and (a′–e′) intermolecular interactions of docked RA (28) with (a,a′) FYN, (b,b′) PDGFRB, (c,c′) PIK3R1, (d,d′) ITGB3, and (e,e′) PDGFRA.
Plants 12 03006 g006
Table 1. Total phytochemical content of P. frutescens extracts.
Table 1. Total phytochemical content of P. frutescens extracts.
ExtractTPC (mg GAE/g CE)TFC (mg RE/g CE)
PF-S-E32.1602 ± 0.544518.8969 ± 0.9848
PF-S-W30.1641 ± 1.0134N.D.
PF-L-E83.1263 ± 0.943142.6199 ± 1.9120
PF-L-W87.2611 ± 1.084154.3026 ± 1.9616
PF-SD-E22.2391 ± 0.205827.3490 ± 0.8821
PF-SD-W44.6299 ± 0.471511.6966 ± 0.3097
Calibration Curvey = 5.6109x − 0.0192
R2 = 0.9990		y = 5.5923x + 0.0280
R2 = 0.9994
Abbreviations: TPC, total phenolics content; TFC, total flavonoids content; GAE, gallic acid equivalents; CE, crude extract; RE, rutin equivalents; S, stem; L, leaf; SD, seed; E, ethanol extract; W, water extract; N.D., not detected.
Table 2. Extraction yield (g CE/100 g PM) of lyophilized P. frutescens crude extracts.
Table 2. Extraction yield (g CE/100 g PM) of lyophilized P. frutescens crude extracts.
Plant Material
SolventStemLeafSeed
Ethanol1.88566.53103.2080
Water4.906015.56864.0512
Abbreviations: CE, crude extract; PM, dried plant material.
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.

Share and Cite

MDPI and ACS Style

Garcia, P.J.B.; Huang, S.K.-H.; De Castro-Cruz, K.A.; Leron, R.B.; Tsai, P.-W. An In Vitro Evaluation and Network Pharmacology Analysis of Prospective Anti-Prostate Cancer Activity from Perilla frutescens. Plants 2023, 12, 3006. https://doi.org/10.3390/plants12163006

AMA Style

Garcia PJB, Huang SK-H, De Castro-Cruz KA, Leron RB, Tsai P-W. An In Vitro Evaluation and Network Pharmacology Analysis of Prospective Anti-Prostate Cancer Activity from Perilla frutescens. Plants. 2023; 12(16):3006. https://doi.org/10.3390/plants12163006

Chicago/Turabian Style

Garcia, Patrick Jay B., Steven Kuan-Hua Huang, Kathlia A. De Castro-Cruz, Rhoda B. Leron, and Po-Wei Tsai. 2023. "An In Vitro Evaluation and Network Pharmacology Analysis of Prospective Anti-Prostate Cancer Activity from Perilla frutescens" Plants 12, no. 16: 3006. https://doi.org/10.3390/plants12163006

APA Style

Garcia, P. J. B., Huang, S. K.-H., De Castro-Cruz, K. A., Leron, R. B., & Tsai, P.-W. (2023). An In Vitro Evaluation and Network Pharmacology Analysis of Prospective Anti-Prostate Cancer Activity from Perilla frutescens. Plants, 12(16), 3006. https://doi.org/10.3390/plants12163006

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop