Next Article in Journal
Genomic Variation Underlying the Breeding Selection of Quinoa Varieties Longli-4 and CA3-1 in China
Next Article in Special Issue
From Nucleus to Organs: Insights of Aryl Hydrocarbon Receptor Molecular Mechanisms
Previous Article in Journal
Methyl Gallate Alleviates Acute Ulcerative Colitis by Modulating Gut Microbiota and Inhibiting TLR4/NF-κB Pathway
Previous Article in Special Issue
Crosstalk between CXCR4/ACKR3 and EGFR Signaling in Breast Cancer Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deciphering the Role and Signaling Pathways of PKCα in Luminal A Breast Cancer Cells

by
Emilio M. Serrano-López
1,2,†,
Teresa Coronado-Parra
1,3,†,
Consuelo Marín-Vicente
1,4,
Zoltan Szallasi
5,6,
Victoria Gómez-Abellán
1,7,
María José López-Andreo
1,8,
Marcos Gragera
1,9,
Juan C. Gómez-Fernández
1,2,
Rubén López-Nicolás
1,2,10,* and
Senena Corbalán-García
1,2,*
1
Department of Biochemistry and Molecular Biology A, Veterinary School, Universidad de Murcia, CEIR Campus Mare Nostrum (CMN), 30100 Murcia, Spain
2
Instituto Murciano de Investigación Biosanitaria IMIB-Arrixaca, El Palmar, 30120 Murcia, Spain
3
Microscopy Core Unit, Área Científica y Técnica de Investigación, Universidad de Murcia, 30100 Murcia, Spain
4
Cardiovascular Proteomics and Developmental Biology Program, Centro Nacional de Investigaciones Cardiovasculares (CNIC), 28029 Madrid, Spain
5
Computational Health Informatics Program, Boston Children’s Hospital, Harvard Medical School, Boston, MA 02115, USA
6
Department of Bioinformatics, Semmelweis University, H-1092 Budapest, Hungary
7
Department of Cellular Biology and Histology, Biology School, Universidad de Murcia, CEIR Campus Mare Nostrum (CMN), 30100 Murcia, Spain
8
Molecular Biology Unit, Área Científica y Técnica de Investigación, Universidad de Murcia, 30100 Murcia, Spain
9
Centro Nacional Biotecnología, Consejo Superior de Investigaciones Científicas, 28006 Madrid, Spain
10
Department of Bromatology and Nutrition, Veterinary School, Universidad de Murcia, CEIR Campus Mare Nostrum (CMN), 30100 Murcia, Spain
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(22), 14023; https://doi.org/10.3390/ijms232214023
Submission received: 18 September 2022 / Revised: 11 October 2022 / Accepted: 26 October 2022 / Published: 14 November 2022
(This article belongs to the Special Issue State-of-the-Art Biochemistry in Spain)

Abstract

:
Protein kinase C (PKC) comprises a family of highly related serine/threonine protein kinases involved in multiple signaling pathways, which control cell proliferation, survival, and differentiation. The role of PKCα in cancer has been studied for many years. However, it has been impossible to establish whether PKCα acts as an oncogene or a tumor suppressor. Here, we analyzed the importance of PKCα in cellular processes such as proliferation, migration, or apoptosis by inhibiting its gene expression in a luminal A breast cancer cell line (MCF-7). Differential expression analysis and phospho-kinase arrays of PKCα-KD vs. PKCα-WT MCF-7 cells identified an essential set of proteins and oncogenic kinases of the JAK/STAT and PI3K/AKT pathways that were down-regulated, whereas IGF1R, ERK1/2, and p53 were up-regulated. In addition, unexpected genes related to the interferon pathway appeared down-regulated, while PLC, ERBB4, or PDGFA displayed up-regulated. The integration of this information clearly showed us the usefulness of inhibiting a multifunctional kinase-like PKCα in the first step to control the tumor phenotype. Then allowing us to design a possible selection of specific inhibitors for the unexpected up-regulated pathways to further provide a second step of treatment to inhibit the proliferation and migration of MCF-7 cells. The results of this study suggest that PKCα plays an oncogenic role in this type of breast cancer model. In addition, it reveals the signaling mode of PKCα at both gene expression and kinase activation. In this way, a wide range of proteins can implement a new strategy to fine-tune the control of crucial functions in these cells and pave the way for designing targeted cancer therapies.

Graphical Abstract

1. Introduction

Breast cancer is the most commonly diagnosed cancer and the leading cause of cancer death among women worldwide [1]. Nowadays, breast cancer is a pathological disease that requires ongoing studies to obtain early prognostic and better treatments. Several signaling proteins are promising targets to confront breast cancer malignancy, and Protein Kinase C (PKC) family could be a potential candidate for designing new therapies [2].
Although PKC was identified 40 years ago by Nishizuka et al., like proteolysis-activated kinases [3], it was not until the discovery of PKC activation by tumor promotor phorbol ester that they turned into a new target in cancer studies [4]. This family has more than 60,000 citations in PubMed, and more than 10,000 citations are related to cancer. PKC family consists of at least ten members encoded by nine genes in mammals, which play essential roles in multiple signaling pathways such as cell proliferation, survival, and differentiation [5]. Notably, it has been impossible to establish whether PKCα isoform acts as an oncogene or tumor suppressor.
Mammary gland differentiation processes regulate the localization and expression of PKCs, and the overexpression of some PKCs in breast cancer has also been described [6]. Moreover, PKCs regulate several processes related to this disease, such as apoptotic and mitogenic signals [7]. Several reports indicate high PKC expression levels in this type of pathology, indicating the use of PKC as a therapeutic target, especially in the case of classical PKCα [8,9]. Approximately 8% of 1084 patients with invasive breast carcinoma showed alterations in PKCα, and 6% of them (67 cases) correspond to amplifications of this gene (cBioportal: https://www.cbioportal.org/) (accessed on 1 May 2022) [10,11].
The role of PKCα in tumor growth and its involvement in tumor progression is well known [12,13,14]. Particularly in breast cancer, the overexpression of PKCα confers a more aggressive phenotype. For example, cell lines that are Estrogen Receptor (ER)-positive become ER-negative in hormone-dependent breast cancer cells [15,16,17]. MDA-MB-231 cells, an example of triple-negative breast cancer cells, also exhibit high expression of PKCα, providing them with an enormous proliferative, migratory, and invasive capacity [18]. In addition, PKCα has also been found in the epithelial–mesenchymal transition process and confers high invasive and motility capacity to breast cancer cells [7,18,19].
To date, a wide variety of PKC inhibitors exist, and, despite numerous clinical trials in cancer, most of them have failed as therapeutic tools [20]. Since the scientific community understood how vital the phosphorylation events are in different signaling pathways, the activity to identify their dysregulation as hallmarks of cancer processes and the search for specific inhibitors has not ceased [21]. Imatinib is the first drug developed targeting a specific protein kinase (ABL tyrosine kinase). This kinase, expressed in nearly all cases of chronic myeloid leukemia (CML), paved the way for the design of new medicine based on the clinical development of kinase inhibitors [22]. Most kinase-targeted cancer therapies use multiple kinase inhibitors selected after screening and detecting critical markers in the patient’s tumor [23]. Many kinases have been considered potential targets for anticancer therapy [24]. Currently, some of the most important therapies inhibit the signaling of EGFR [25], PI3K/AKT/mTOR [26,27], CDK4/6 [28], or AURKA [29].
This strategy provides many examples of beneficial clinical applications but also faces many challenges and limitations due to the resistance that cancer cells develop to chemotherapy or the different response of patients to the same therapy. Therefore, the new therapies seek to combine the simultaneous inhibition of several kinases acting in the same pathway (vertical inhibition) or complementary pathways (horizontal inhibition) [30]. Recently, interest in developing so-called “personalized medicine” is leading to the creation of computational models combining multiple drugs to attack different targets and, at different times, improve patient treatment [31].
All this knowledge and experience led us to think about the possibility of systematically using a two-step strategy to inhibit PKCα expression. First, to slow down the tumor process and then to survey the signaling pathways affected. Second, to use the knowledge obtained to apply a set of specific inhibitors to modulate the activity of these cancer cells. Thus, we studied the effect of inhibiting PKCα expression using siRNA on cellular processes such as proliferation, migration, or apoptosis in a luminal A breast cancer cell line (MCF-7). Here, we explore the differential gene expression and kinase arrays to compare PKCα-Knock-Down (PKCα-KD) vs. PKCα-Wild Type (PKCα-WT) MCF-7 cells. Metascape functional enrichment analysis was used to ascertain the most critical pathways that were down- and up-regulated. Together with a phospho-kinase array, we identified a set of down-modulated signaling pathways that correlated well with controlling the growth and migration capacity observed in these cells. Analysis of the up-regulated pathways helped us select a set of targeted inhibitors to test whether their effect was synergistic with that of PKCα. We found that the inhibition of the PLC or tyrosine kinase pathways strongly affects proliferation and migration, further counteracting the tumorigenic phenotype.

2. Results

2.1. Knock-Down of PKCα in MCF-7 Cells Affects Their Tumorigenic Capacity

We wanted to elucidate whether the reduction in PKCα protein level influenced the tumor phenotype of this breast cancer cell line. Therefore, PKCα was knocked down in the MCF-7 breast cancer cell line using specific siRNA. We used a Western blot to measure the protein expression level with a specific anti-PKCα antibody (Figure 1A). The results show inhibition of PKCα expression of more than 55% from day 3 to day 8, reaching the maximum inhibition (70%) between days 4 and 6. During the same period (8 days), PKCα-KD cells showed a 48% reduction in proliferation on day 6 that persisted till day 8, indicating the role of PKCα in the proliferation process (Figure 1B).
As tumor expansion capacity, migration was also measured in PKCα-WT and compared to PKCα-KD MCF-7 cells by wound healing assay (Figure S1A). A 61% reduction in cell migration was observed in PKCα-KD cells after 48 h, indicating that downstream phosphorylation by PKCα might be a critical step in reorganizing the cell cytoskeleton. PKCα-KD showed a slight increase in apoptosis (Figure S1B). Overall, the results indicate that PKCα depletion reduces the tumor phenotype mainly by decreasing the proliferation and migration of this breast cancer model.

2.2. Transcriptomic Analysis of PKCα Knock-Down MCF-7 Cells Reveals Essential Alterations in Signaling Pathways

We interrogated the molecular pathways altered in MCF-7 breast cancer cells by decreasing PKCα protein expression. Thus, transcriptomic analysis was studied, and the mRNA expression transcripts of PKCα-WT and PKCα-KD MCF-7 cells were compared. A total of 13.461 unique mRNA transcripts were analyzed and the datasets were represented with the help of a volcano plot (Figure 2A). More than 5.000 probes were significantly different between the two cell types, displayed above the horizontal line (p < 0.05). When comparing PKCα-KD vs. PKCα-WT, 292 probes showed 1.5-fold up-regulation, whereas 248 probes showed 1.5-fold down-regulation (Figure 2A). Several genes distributed at different positions of the volcano plot were used for microarray validation by RT-qPCR (Figure 2A,B). The down-regulated PRKCA, EGFR, and ITGβ6 genes, as well as up-regulated PRKAR2B and unaffected SPAM1, were validated (Figure 2B and Table S1). Some genes (ERBB4, PDGFA, and PLCB4) showed low expression differences when validated by RT-qPCR, probably due to differences in cDNA synthesis and amplification methods compared to microarray.
We used Metascape (https://metascape.org/) accessed on 24 March 2022 [32] to perform the functional enrichment analysis of differentially expressed genes (DEGs) to better understand the functions and metabolic pathways after the knock-down of PKCα in MCF-7 cells. We studied up- and down-regulated DEGs separated into two sets for a more comprehensive analysis. As shown in Figure 3, the up-regulated genes group showed that the ten most overrepresented biological processes were: cellular response to hormone stimulus (GO:0032870), protein phosphorylation (GO:0006468), regulation of cellular response to stress (GO:0080135), intracellular protein transport (GO:0006886), regulation of cellular localization (GO:0060341), pathways in cancer (HSA05200), neuron projection development (GO:0031175), RHO GTPase cycle (R-HAS-9012999), regulation of growth (GO:0040008), and signal transduction by p53 class mediator (GO:0072331). The eleven most overrepresented biological processes in down-regulated genes were: Interferon Signaling (R-HSA-913531), regulation of response to biotic stimulus (GO:002831), antiviral mechanism by IFN-stimulated genes (R-HSA-1169410), negative regulation of cell population proliferation (GO:0008285), immune response to tuberculosis (WP4197), regulation of peptidase activity (GO:0052547), ER-phagosome pathway (R-HAS-1236974), response to wounding (GO:0009611), non-genomic actions of 1,25 dihydroxy vitamin D3 (WP4341), modulation by symbiont of entry into the host (GO:0052372), and signaling by receptor tyrosine kinases (R-HAS-9006934).

2.3. PKCα Controls the Phosphorylation of Essential Ser/Thr and Tyr Kinases and Transcription Factors

To obtain further insights into the kinases affected by the lack of PKCα, we performed an indicative assay using a human phospho-kinase antibody array. We compared PKCα-KD vs. PKCα-WT MCF-7 cell lysates (Figure S2). The array showed 14 proteins with lower phosphorylation levels in the PKCα-KD cells (Figure 4). These included PDGFr and AKT isoforms, a group of tyrosine kinases such as Lck, Fgr, or Hck, different isoforms of the STAT and β-catenin transcription factors, as well as DNA damage-related proteins Chk-2 and p27.
Although most proteins experienced decreased phosphorylation, surprisingly, we found that ERK1/2 and p53 (S392) were phosphorylated at a higher level when PKCα was down-regulated by siRNA treatment (Figure 4).
To obtain a complete overview of the complex cellular signaling depending on PKCα in MCF-7 cells, we analyzed the interactions and connections between the differentially expressed genes and the kinases affected in the absence of PKCα. We analyzed the protein–protein interaction (PPI) network using the STRING database [33] and Cytoscape software [34] (Figure 5). The PPI network consisted of 502 nodes and 1139 edges. The Cytoscape tool MCODE [35] identified the highly interconnected proteins. Two clusters appeared: Cluster 1 (22 nodes and 222 edges) and Cluster 2 (17 nodes and 51 edges) (Figure 5A). Cluster 1 is composed of proteins whose expression decreased after PKCα inhibition. The functional analysis of both clusters using the ClueGO tool [36] (Gene Ontology (Biological process) and KEGG database) revealed that cluster 1 is related to immune system regulation through interferon signaling and antigen processing and preparation via MHC I (Figure 5B). Cluster 2 showed a complex interaction between proteins and kinases. The main pathways were related to the epidermal growth factor signaling and the regulation of glucose or nucleocytoplasmic transport in the cell. KEGG functions such as proteoglycans in cancer or thyroid hormone signaling pathways (Figure 5C) also appeared.

2.4. Potential Targeted Therapy by Using Specific Inhibitors of Key Signaling Pathways

This work indicates that suppression of PKCα induces the down-regulation of critical oncogenic pathways. However, this study also shows unexpected up-regulation of crucial genes and kinases (e.g., ERBB4, PLCβ4, PRKAR2B, and PDGFA) whose roles in cancer development and progression are crucial (Figure 3). Therefore, we propose a dual strategy involving drugs that inhibit some of the up-regulated signaling pathways in PKCα-KD cells as a mode of targeted therapy to augment the inhibitory effect already exerted by PKCα suppression.
Thus, we used four different commercial inhibitors: U73122 (PLC and PLA2 families), KT5720 (PKA, catalytic domain), BMS 599626 (ErbB receptor family, catalytic domain), and Imatinib (ABL, KIT, and PDGF tyrosine kinases). We measured their effect on proliferation, migration, and apoptosis in PKCα-KD vs. PKCα-WT MCF-7 cells (Figure 6).
Inhibition of the PLC family by U73122 impairs the hydrolysis reaction of PIP2 into IP3 and DAG [37]. The effect of U73122 on PKCα-KD vs. PKCα-WT MCF-7 cells was measured using concentrations at 1 and 10 μM (Figure 6A). The results showed more significant inhibition of the proliferation rate in PKCα-KD cells (88%) compared to 49% in the same cells in the absence of the drug and 42% in the presence of PKCα and drug on day 8 (Figure 6A; Table S3). U73122 treatment also produced 89% inhibition in the migration of PKCα-KD cells at 10 µM concentration compared to the 53% produced in the control cells at 72 h (Figure 6B; Table S4). U73122 only promoted significant apoptosis at 20 µM in PKCα-KD cells (Table S5).
Treatment with KT5720 at 0.1 and 1 μM, to inhibit the catalytic domain of PKA, produced a slight decrease in the proliferation of PKCα-KD in addition to the absence of PKCα (Figure 6C; Table S3). The addition of 1 μM KT7520 reduced the ability of cells to migrate (89%) compared to 53% in the absence of the drug after 72 h (Figure 6D; Table S4). However, the KT5720 treatment did not affect apoptosis (Table S5).
BMS 599626 is a selective and efficacious inhibitor of the ErbB receptor family and is used at a wide range of concentrations (0.1, 1, 5, and 10 µM). The proliferative activity of PKCα-KD cells decreased to 92% compared to 45% in the absence of the drug. However, the drug also induced a robust inhibitory effect (82%) in PKCα-WT cells, which implies only an extra 10% inhibition by combining the two strategies (Figure 6E; Table S3). Combining the inhibition of PKCα and BMS 599626, we obtained 80% inhibition compared to control cells treated with the same drug concentration (Figure 6F; Table S4). Regarding apoptosis, combinatory treatment hardly improved the results of each treatment; although, a high concentration of 20μM showed a slight difference between PKCα-WT and PKCα-KD cells (Table S5).
We used Imatinib at two different concentrations: 1 and 10 µM. We observed that imatinib (10 µM) decreased the proliferation rate of PKCα-WT MCF-7 cells above 47% (Figure 6G; Table S3). When this treatment was applied to PKCα-KD cells, the proliferation inhibition was 85% compared to the control cells (Figure 6G; Table S3). Moreover, the effect of Imatinib on the migration capacity hardly decreased in PKCα-WT cells. On the contrary, a substantial inhibitory migration effect of 68% appeared in PKCα-KD cells (Figure 6H; Table S4). Imatinib did not exhibit a significant effect on apoptosis (Table S5).
In conclusion, the results indicate that inhibitors of the PLC and tyrosine kinase pathways are more appropriate to control the proliferation of PKCα-KD MCF-7 cells and to control migration, inhibitors of the PLC, PKA, and EGFR pathways seem to be more appropriate. Other more potent drugs that modulate the PLC and EGFR pathways in that direction should be applied to induce apoptosis.

3. Discussion

PKCα and other isoforms are potential targets for cancer therapy, and it is challenging to determine the specific role of each isoform for different cancer types. PKCα is considered a tumor promoter, but in recent years, reports show that it can also function as a tumor suppressor, and its activity regulates cellular processes such as apoptosis, proliferation, or migration [14,38,39,40]. This duality in PKCα activity reveals the importance of studying each protein in its specific biological context to understand its role. Our study focused on characterizing the role of PKCα in a specific subtype of breast cancer (Luminal A breast cancer) using MCF-7 cells as a model. The suppression of its protein expression using siRNA has allowed us to identify a set of cellular processes that could be affected in its absence and whose knowledge could allow us to modulate the activity of this cell type [41,42,43]. In this work, we have obtained an overview of the signaling pathways affected in MCF-7 cells after PKCα inhibition by combining the results of DEGs with the activation status of kinases and the application of the MCODE tool (major gene clusters) on the PPI network.
Following PKCα inhibition, we found that interferon signaling was one of the clusters affected. MHC class I isoforms (HLA-A, HLA-B, HLA-C, and HLA-F) are involved in antigen presentation to cytotoxic CD8+ T cells and activate the adaptive immune system response [44,45]. Down-regulation of MHC class I has been described in many cancers as a strategy of tumor cells to avoid the immune system response [46]. Another group of down-regulated proteins participating in the interferon signaling pathway is the family of 2′-5′-oligoadenylate synthetases (OAS), consisting of OAS1, OAS2, OAS3, and OAS-like (OASL). High expression of these proteins correlates with better overall survival in basal-like breast cancer [47]. However, high expression of OAS1 and OAS3 was significantly associated with worse overall survival in Luminal A breast cancer [48]. Further studies will elucidate whether the low expression of these genes for a long time can cause resistance to therapies or evasion of the immune system. The results indicate that, in the short term, the tumor phenotype decreases in its ability to proliferate and migrate (Figure 1).
Proteins involved in wound response were down-regulated, and this would explain the low migration rate observed in wound healing assays of PKCα-KD MCF-7 cells. This group includes ITGA2 and ITGB4, KRT6A, or ELK3, among others (Table S2a). The role of integrins in cell adhesion and their influence on cancer progression is well characterized [49]. Keratins are also diagnostic tumor markers; their low expression is a good indicator of the progression of the tumor [50]. At the same time, ELK3 expression correlates well with cell migration and invasion in the triple-negative breast cancer cell model (MDA-MB-231) [51], and the low expression of these proteins would favor the control of the cancer phenotype in MCF-7 cells.
One of the most common activated pathways in cancer processes is the phosphoinositide 3-kinase (PI3K)-AKT pathway, controlling cellular growth, survival, and key metabolic processes [52]. The PI3K-AKT pathway has been studied in breast cancer because it appears as one of the most affected signaling pathways by mutations, amplifications, and/or deregulation of proteins [53,54]. Our results showed that inhibition of PKCα affected PI3K-AKT signaling by controlling the expression of their regulatory subunits (PI3KR1 and PI3KR2). These essential elements mediate the activity of the catalytic subunit by stabilizing it, inhibiting it, or allowing its interaction with other downstream elements [55,56]. The overexpression of PI3KR1 [57,58] and the inhibition of PI3KR2 [59] directly affect the catalytic capacity of PI3K, decreasing the PI3K-AKT signaling and tumor progression. The PI3K inhibition may result in the down-activation of AKT isoforms (AKT1 and AKT2), suggesting that PKCα positively regulates the PI3K-AKT signaling pathway.
Down-regulation of PKCα affected the expression of several proteins related to signaling by receptor tyrosine kinases (ERLIN2, CDK2AP2, and TNS3, among others). Most of them promoted breast cancer cell survival [60]; in addition, the kinase array also showed four tyrosine kinases with a lower level of phosphorylation, thus reinforcing the idea that PKCα also controls these critical pathways.
The STAT family is a significant element in this complex signaling [61,62]. Our data suggest that PKCα directly or indirectly controls the down-regulation of STAT1 expression and phosphorylation in MCF-7 cells. Several reports showed that PKC isoforms activated these transcription factors in other cell lines [63,64]. The role of the STAT family in cancer is not straightforward, and depending on the cellular context, they function as oncogenes or tumor suppressors [65,66]. Elevated STAT1 levels are associated with therapeutic resistance in ER+ breast cancer cells [67]. Therefore, we could consider suppressing PKCα expression to a beneficial effect in this breast cancer cell line. The same situation occurs with other transcription factors such as SOX2 and SNAI2 that were also down-regulated. Both have been described as essential mediators for tumor progression and proliferation and provide endocrine resistance to conventional anticancer therapy [68,69].
Most of the signaling pathways inhibited due to PKCα down-regulation indicated that this kinase controls key signaling hubs in MCF-7 cells. Its inhibition correlates well with the results observed in proliferation and migration assays (Figure 1). The only pathways that alerted our attention were interferon-related since their downregulation correlates with disease progression and poor prognosis.
Analysis of the pathways regulated by PKCα inhibition also showed a high expression of the insulin-like growth factor type 1 receptor (IGF1R) that correlates with a favorable prognosis in hormone receptor-positive breast cancer [70].
SMAD5 appeared up-regulated but not Aurora kinase A (AURKA), which cooperates to promote chemoresistance in breast cancer therapies [71]. In addition, we identified Aurora kinase binding protein (AURKAIP1) up-regulated. This protein is a direct activator of AURKA degradation, and thus, it is responsible for the pathway inactivation [29]. All these results correlate well with those obtained in our work regarding growth inhibition and cell migration.
ERK1 and ERK2 are master regulators of different pathways, such as cell proliferation, survival, growth, metabolism, migration, and differentiation [72]. Different stimuli trigger the ERK signaling through receptor tyrosine kinases (RTKs), G-protein-coupled receptors, or integrins, in addition to small GTPases Ras and Rap [73]. MAPK kinase signaling is a major culprit in breast cancer. This pathway is also responsible for tumorigenesis and cancer progression, and there are different drugs to inhibit their signaling [74]. However, many studies have shown good prognostic features associated with increased activation of MAPKs in ER+ breast cancers [75]. Numerous studies suggest that PKCα activates ERK signaling [76,77,78,79]. However, our data showed that ERK1/2 experienced increased activation after PKCα inactivation and correlated with decreased proliferation and migration capacity of MCF-7 cells. These results could support the idea that ERK activation mediated by PKCα inhibition has an antitumor effect in ER+ breast cancer.
One of the most important tumor suppressors in cancer research is p53 due to its role in inducing growth arrest or apoptosis after DNA damage or alteration in the cell cycle process, besides controlling other cellular mechanisms such as metabolism [80]. Many processes in the cellular context contribute to determining the decision between cellular survival or death, where p53 acts take the final jury’s verdict [81,82]. Our results showed an increase in p53 phosphorylation and up-regulation of several proteins related to this pathway after PKCα inhibition, supporting our previous results where PKCα-KD MCF-7 cells showed less capability to proliferate compared to PKCα-WT MCF-7 cells. p53 activation could be determined indirectly through AKT activity modulation [83] or direct activation mediated by a lack of PKCα [84,85].
A large variety of unexpected genes (ESR1, ERBB4, PLCβ4, PDGFA, PRKAR2B, IGF1R, SMAD5, or MAPK3K2, among many others, see Table S2b) appeared up-regulated. This fact is an essential finding after applying a PKCα inhibition strategy since, in the short term, the signaling balance indicated the inhibition of the progression toward the tumor phenotype.
The estrogen receptor ESR1 dimerizes and translocates to the nucleus upon ligand binding. It controls the transcription of many genes involved in cell cycle regulation, DNA replication, cell differentiation, or apoptosis. ESR1 also remains in the cytoplasm to interact with other proteins, such as receptors or kinases, participating in various signaling pathways [86]. Overexpression of the ESR1 gene potentially reduces the sensitivity of breast cancer cells to endocrine therapies, leading to disease progression and metastasis [87]. Studies showed that PI3K inhibition increases ER activity by regulating the histone methyltransferase KMT2D activity [88,89]. Likely, the observed inhibitory effect on the PI3K pathway due to the lack of PKCα promotes this unexpected compensatory effect.
Another up-regulated gene involved in protein phosphorylation and cancer pathways was ERBB4, while EGFR was down-regulated. EGFR is a significant target in cancer, as it is frequently mutated and/or overexpressed. Numerous drugs inhibit its activity [90]. Seven ligands connected to EGFR activators induce specific cellular responses and intracellular trafficking (Ras/MAPK pathway, the PI3K/AKT pathway, and the phospholipase C (PLC)/protein kinase C (PKC) signaling) [91]. In addition, other kinase-independent functions result from the heterodimerization of EGFR with other ERBB family members or interaction with kinases/phosphatases that trigger its signaling [92]. ERBB4 differs from the rest of the HER family members by undergoing alternative splicing and has four different isoforms depending on where the splicing site locates: the extracellular juxtamembrane region (JMa or JMb) or the cytosolic C-terminus, close to the tyrosine kinase domain (Cyt1 and Cyt2) [93]. Upon interaction with ligand, ERBB4 is activated and begins to phosphorylate its substrates, triggering different signaling pathways such as Ras-MEK-ERK or PI3K-AKT or promoting gene expression through its nuclear translocation [94]. Although ERBB4 plays an essential role in different types of cancer, its role as an oncogene or tumor suppressor depends on the cellular context, and it is still under study [95]. Also important is the endosomal recycling and degradation pathway that acts as the primary regulator of HER family expression and activation by controlling the availability of these receptors in the cell membrane and other cellular compartments [96]. The microarray results indicate decreased expression of the endosomal machinery. Numerous studies have described the role of PKC isoforms in regulating the HER family by controlling their travel from the membrane to endocytic vesicles. This traffic determines what vesicle proportion will be recycled and which one will degrade in the lysosome [97,98,99,100,101,102]. The results of proliferation and migration experiments indicate that the balance between EGFR and ERBB4 expression shifts towards inhibiting the tumor phenotype of MCF-7 cells.
The DEGs analysis also detected overexpression of the PRKAR2B gene encoding for the regulatory subunit of PKA. This result appears as a favorable prognostic marker in pancreatic and colorectal cancer and unfavorable in thyroid cancer. On the contrary, other work has shown that PRKAR2B depletion induces resistance to apoptosis, and in that sense, overexpression of the gene would be desirable [103,104].
Several genes, such as PDGFA or PLCβ4 exhibited increased expression and are usually associated with poor prognosis in patients with solid tumors [105,106].
These latter results were significant in designing the rationale for selecting inhibitory drugs (Figure 6) to counteract the over-expression of these genes to push the balance toward a more inhibited tumor phenotype. This strategy demonstrated that inhibiting PKCα in the first step makes it possible to condition the MCF-7 cells into a low tumor profile due to the inhibition of critical cancer pathways (Figure 3 and Figure 4). In the second step, the specific strategy to inhibit unexpected up-regulated pathways allows us to further hinder this phenotype, with PLC and tyrosine kinases being the most important for proliferation and PLC, PKA, and EGFR being the most important for migration (Figure 6).
In conclusion, we have shed light on the role of PKCα in luminal A breast cancer cells. The integration of this information has allowed us to advance in two significant areas in cancer therapy. First, inhibiting a kinase that is not the leading cause of the tumorigenic effect is an appropriate strategy to curb the tumor phenotype. Second, to fully control the process, we can use differential expression and kinase activation to establish a personalized treatment based on the specific knowledge of the proteins affected. Implicitly, this research seeds the idea of a workflow for targeted therapy. Shortly, if perfected, it might be applied to a clinical performance protocol by testing in vitro this dual strategy on tumor cells extracted from the patient to precisely tailor the treatment.

4. Materials and Methods

Cell culture. MCF-7 breast cancer cells were purchased from ATCC and were cultured in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% FBS, 50 uds/mL penicillin, and streptomycin, 110 mg/mL pyruvate and 2 mM glutamine, in a 7.5% CO2 incubator at 37 °C.
siRNA transfection. PKCα expression was inhibited in MCF-7 cells using siRNA (Invitrogen, Waltham, MA, USA). In addition, a negative siRNA control was used in this experiment. The sequence used to inhibit PKCα expression was: 5′-UCCAAACGGGCUUUCAGAUCCUUAU -3′ (PKCα-KD). The negative control siRNA sequence was 5′-AGGUAGUGUAAUCGCCUUGUCGCCUUG-3′ (PKCα-WT).
MCF-7 cells were transfected with 1 µM siRNA (control non-target or inhibitory siRNA) using Mirus Ingenio® Electroporation Solution (Mirus Bio LLC, Madison, WI, USA) by electroporation. Then, 3 × 106 cells were resuspended in 300 μL of electroporation solution and added to 0.4 cm electroporation cuvettes (Bio-Rad, Hercules, CA, USA). Cells were electroporated in a GenePulser (Bio-Rad, Hercules, CA, USA) with two 900 V/0.4 ms square wave pulses and immediately placed on a fresh medium.
Western blot. Cells were lysed in lysis buffer (50 mM Tris HCl pH 7.5, 150 mM NaCl, 1% Nonidet P-40 (Roche Diagnostic, GmbH, Germany), 10% glycerol, complete commercial protease inhibitor, Mini, EDTA-free Protease Inhibitor Tablets and phosphatases inhibitor PhosSTOP (Sigma-Aldrich Chemistry, S.A., Madrid, Spain). After cellular lysis, cells were passed through a syringe (10×) and centrifugated at 14,000 rpm for 15 min at 4 °C. Electrophoresis loading buffer (200 mM Tris HCl pH 6.8, 250 mM DTT, 5% SDS, 37.5% glycerol, and 0.015% bromophenol blue) was added to supernatants and analyzed into polyacrylamide gel (SDS-PAGE). Proteins were separated and transferred to a nitrocellulose membrane (Hybond ECL (GE)). The membrane was blocked using blocking buffer (2% bovine serum albumin diluted in TBST (Tris-HCl 20 mM pH 7.5, NaCl 150 mM and Tween-20 0.1% (v/v)) for 2 h at room temperature. Then, primary antibodies were incubated O/N at 4 °C, and horse peroxidase-conjugated secondary antibodies were incubated for one hour at room temperature. WesternBright ECL-HRP substrate commercial kit (Advansta) was used to reveal the bands. The optical densities of Western blot bands analyzed with Fiji software [107] (https://imagej.net/software/fiji/downloads) calculated the percentage of PKCα inhibition and normalized by β-actin and using the level of each day PKCα expression as 100%. The bar plot shows the average percentage of inhibition per day ± standard deviation. Primary antibodies used were: anti-PKCα ab32376 (Abcam), anti-GAPDH ab9485 (Abcam and anti-β actin ab8227 (Abcam). Secondary antibodies used were: Goat F(ab’)2 Anti-Mouse IgG—F(ab’)2 (HRP), pre-adsorbed (ab5887) and Goat F(ab’)2 Anti-Rabbit IgG F(ab’)2 (HRP) preadsorbed (ab6112) (Abcam).
Compounds. U73122 and KT5720 were purchased from Sigma-Aldrich (St. Louis, MO, USA); BMS 599626 was purchased from Selleckchem (Houston, TX, USA), and Imatinib was purchased from Novartis (Basel, Switzerland). Compounds were dissolved in sterile DMSO, and they were used at the following concentrations: U73122 (1–10 µM), KT5720 (0.1–1 µM), BMS 599626 (0.1–10 µM), and Imatinib (1–10 µM).
Wound-healing assay. A wound-healing assay was used to evaluate the migration ability. Control and siRNA-transfected MCF-7 cells were seeded in 24-well plates to achieve a 90–95% confluence on the day of maximum inhibitory. The day before the assay, cell medium was replaced by low concentration FBS medium to synchronize the cellular cycle. The plates were washed with PBS after making a scratch in each well using a sterile pipette tip. The cultures were photographed until the wound was closed. In addition, different compound (previously described) was added to evaluate their effect on cell migration capability. Quantitative analysis of the wound area was carried out using Fiji [107] software from three independent wound-healing experiments.
Proliferation assay. Proliferation assay was assessed by fluorometric quantification of DNA using CyQUANT® NF Cell Proliferation Assay kit (Invitrogen, Waltham, MA, USA), according to the manufacturer’s instructions. Briefly, control and siRNA transfected cells were seeded in 96-well plates (1 × 104 cells) and allowed to attach for 24 h. Control and siRNA MCF-7 cell medium were replaced by medium containing drugs when appropriate. Then, the proliferation assay was measured for 8 days, where each condition was tested in three wells and three independent experiments at least.
Apoptosis assay. The percentage of apoptotic cells was determined by flow cytometry using the Vybrant apoptosis assay kit #4 (YO-PRO-1/propidium iodide; Molecular Probes/Invitrogen). Briefly, cells were washed twice with ice-cold PBS and resuspended in PBS containing YO-PRO-1/propidium iodide. Apoptotic cells were identified by flow cytometry after incubation for 20 min. Briefly, control and siRNA transfected cells were seeded in 96-well plates (1 × 104 cells) and allowed to attach for 24 h. Control and siRNA MCF-7 cell medium were replaced by medium containing drugs where appropriate. Cells were treated until the day of maximum inhibitory. Then, apoptosis assay was measured, where each condition was tested in three wells and three independent experiments at least.
RNA extraction. RNA was extracted from cultured cells with the QIAGEN RNeasy pPlusMini Kit according to the manufacturer’s protocol. MCF-7 cells were transfected with control and inhibitory siRNA and RNA was extracted on the day of maximum PKC α expression inhibition. RNA concentration and purity were measured on a NanoDrop ND-1000 spectrophotometer (NanoDropTech, Wilmington, DE, USA).
Microarray. RNA isolated from PKCα-KD and PKCα-WT MCF-7 cells on the day of maximum inhibitory was analyzed using GeneChip® Probe Array HG-U133 A2 (Affymetrix) according to the manufacturer’s instructions. Briefly, previously isolated RNA was amplificated and biotinylated using MessageAmp™ II-Biotin Enhanced Kit (ThermoFisher Scientific) according to the manufacturer’s protocol. Then, new biotinylated RNA was fragmented and hybridized with GeneChip® Probe Array for 16 h at 45 °C. Chips were scanned using Affymetrix® GeneChip® Scanner 3000. Triplicates of each condition were performed. Microarray quality control was performed with PCA analysis (m) after data normalization with the RMA algorithm. Results were analyzed using Partek Genomic Suite (Partek Incorporated, St. Luis, MO, USA) software to determine genes with a significative different expression between both conditions (Table S6).
Microarray data analysis. The transcripts selected using the above procedures were analyzed using Metascape web-based portal [32]. Metascape is a powerful web-based tool for gene annotation and gene list enrichment analysis that incorporates updated ontology databases such as KEGG Pathway, GO Biological Processes, Reactome Gene Sets, Canonical Pathways, and CORUM. To perform the analysis, all genes in the genome proceeded from the enrichment background. Terms with a p-value <0.01, a minimum count of 3, and an enrichment factor >1.5 are collected and grouped into clusters based on their membership similarities. In addition, the nodes that share the same cluster-ID are typically close to each other. The terminology used to name the transcripts is that approved by the HGNC.
RT-qPCR. RNA was harvested from transfected MCF-7 cells as described before and reverse transcribed using iScript Reverse Transcription Supermix kit (Bio-RadBio-RadLaboratories, Hercules, CA, USA) following the manufacturer’s instructions. The RT-qPCR was performed from cDNA templates using the SYBR Green PCR master mix and a 7500 fast real-time PCR system (Applied Biosystems). The transcript levels were determined after normalization against GAPDH and SPAM1, using the REST© 2022 software. All RT-qPCR experiments were repeated with at least three biologically independent replicates. Table S7 shows the list of the primers used in this study.
Protein–Protein Interaction (PPI) network analysis. In the present study, Protein–Protein Interaction (PPI) network of Differential Expression Genes (DEGs) from PKCα knock-down MCF-7 cells was constructed using the online database Search Tool for the Retrieval of Interacting Genes (STRING; http://string-db.org) (version 11.5) [33]. An interaction with a combined score >0.7 was considered statistically significant. The thicknesses of those edges were associated with the combined score. The STRING interaction network was imported into Cytoscape software (version 3.9.1) [34]. Cytoscape is an open-source bioinformatics software platform for visualizing molecular interaction networks and further analysis. Cytoscape app MCODE [35] was used to visualize the significant nodes and partition the network into different modules with degree cut-off = 2, cluster finding = haircut, node score cut-off = 0.2, K-score = 2, and maximum depth = 100, respectively. Functional enrichment analysis of identified clusters was performed by using ClueGo (version 2.5.8) [36] plugin in Cytoscape environment. The parameters used for ClueGo analysis were a two-sided (Enrichment/Depletion) tests based on the hypergeometric distribution for enrichment analysis. The p value  <  0.05 was corrected by Bonferroni step-down correction method. Only Gene Ontology (Biological Processes) and KEGG databases were used in the analysis. Terms that passed the p value threshold (p  <  0.05) were considered significantly enriched.
Proteome profiler human phospho-kinase array. The relative levels of protein phosphorylation were determined by proteome profiler human phospho-kinase array kit (Cat#: ARY003B, R&D Systems, Minneapolis, MN, USA). Protein extraction was performed as mentioned in Western blot assay. Briefly, 400 μg of each MCF-7 cell lysate (PKCα-KD and PKCα-WT) were incubated with the microarrays pre-coated with antibodies against 43 kinase phosphorylation sites at 4 °C overnight. After washing, the microarrays were incubated with biotin-labeled antibodies at RT for 2 h. Then microarrays were incubated with HRP-conjugated streptavidin at RT for 30 min on a rocking platform. Chemiluminescence and Density of spot pixels were analyzed using ImageJ software. Each membrane has its controls described in Figure S2, and the intensities are normalized with respect to the 100% signal (reference spot). Background signal from each spot was subtracted using the average signal from the negative control spots. Note that each kinase is tested in two independent spots, and several kinases of the same family are represented. Only one biological sample was used in this assay.
Statistical analysis. Statistical data analysis and visualization were performed with Estimation Statistics database (www.estimationstats.com) [108]. All data points were presented in a Cumming estimation plot to display the underlying distribution. The raw data are plotted on the upper axis. On the lower axis, mean differences are plotted as bootstrap sampling distributions (5000). Each mean difference is depicted as a dot. Each 95% confidence interval is indicated by the ends of the vertical error bars. An additional permutation test was also applied to confirm these results. The p values to accept/reject the null hypothesis of no differences between PKCα-WT and PKCα-KD cells are indicated below each comparison (CI 95%).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232214023/s1.

Author Contributions

Conceptualization, S.C.-G.; methodology and software R.L.-N., S.C.-G., E.M.S.-L., M.J.L.-A. and T.C.-P.; validation, R.L.-N., C.M.-V., T.C.-P., V.G.-A., M.G. and Z.S.; formal analysis, R.L.-N., T.C.-P. and E.M.S.-L.; investigation, R.L.-N., T.C.-P., C.M.-V., V.G.-A., M.G., M.J.L.-A. and J.C.G.-F.; writing—original draft preparation, E.M.S.-L., T.C.-P. and S.C.-G.; writing—review and editing, S.C.-G., J.C.G.-F., Z.S., C.M.-V., E.M.S.-L., T.C.-P. and R.L.-N.; supervision, S.C.-G.; project administration and funding acquisition, S.C.-G. and J.C.G.-F. All authors have read and agreed to the published version of the manuscript.

Funding

Work in Murcia was supported by grants BFU2017-87222-P (MICINN, Spain-FEDER) to S.C.-G. and J.C.G.-F. and Fundación Séneca Region de Murcia 20885/PI/18 to S.C.-G.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Microarray data are available in the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-12374.

Acknowledgments

We would like to dedicate this work to the memory of Juan Carmelo Gomez-Fernandez, our mentor in the Biomembranes group at the University of Murcia, to whom we will always be grateful for his guidance and support in our continuous search for advancing scientific knowledge in this field. In addition, we would like to thank the Cell Culture, Molecular Biology, Image, and Microscopy Core Facilities at the Área Científica y Técnica de Investigación (ACTI) at the Universidad de Murcia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lei, S.; Zheng, R.; Zhang, S.; Wang, S.; Chen, R.; Sun, K.; Zeng, H.; Zhou, J.; Wei, W. Global Patterns of Breast Cancer Incidence and Mortality: A Population-Based Cancer Registry Data Analysis from 2000 to 2020. Cancer Commun. 2021, 41, 1183–1194. [Google Scholar] [CrossRef] [PubMed]
  2. Lønne, G.K.; Cornmark, L.; Zahirovic, I.O.; Landberg, G.; Jirström, K.; Larsson, C. PKCα Expression Is a Marker for Breast Cancer Aggressiveness. Mol. Cancer 2010, 9, 76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Inoue, M.; Kishimoto, A.; Takai, Y.; Nishizuka, Y. Studies on a Cyclic Nucleotide-Independent Protein Kinase and Its Proenzyme in Mammalian Tissues. II. Proenzyme and Its Activation by Calcium-Dependent Protease from Rat Brain. J. Biol. Chem. 1977, 252, 7610–7616. [Google Scholar] [CrossRef]
  4. Castagna, M.; Takai, Y.; Kaibuchi, K.; Sano, K.; Kikkawa, U.; Nishizuka, Y. Direct Activation of Calcium-Activated, Phospholipid-Dependent Protein Kinase by Tumor-Promoting Phorbol Esters. J. Biol. Chem. 1982, 257, 7847–7851. [Google Scholar] [CrossRef]
  5. Cooke, M.; Magimaidas, A.; Casado-Medrano, V.; Kazanietz, M.G. Protein Kinase C in Cancer: The Top Five Unanswered Questions. Mol. Carcinog. 2017, 56, 1531–1542. [Google Scholar] [CrossRef]
  6. Masso-Welch, P.A.; Verstovsek, G.; Ip, M.M. Alterations in the Expression and Localization of Protein Kinase C Isoforms during Mammary Gland Differentiation. Eur. J. Cell Biol. 1999, 78, 497–510. [Google Scholar] [CrossRef]
  7. Urtreger, A.J.; Kazanietz, M.G.; Bal De Kier Joffé, E.D. Contribution of Individual PKC Isoforms to Breast Cancer Progression. IUBMB Life 2012, 64, 18–26. [Google Scholar] [CrossRef]
  8. Lahn, M.; Köhler, G.; Sundell, K.; Su, C.; Li, S.; Paterson, B.M.; Bumol, T.F. Protein Kinase C Alpha Expression in Breast and Ovarian Cancer. Oncology 2004, 67, 1–10. [Google Scholar] [CrossRef]
  9. Tan, M.; Li, P.; Sun, M.; Yin, G.; Yu, D. Upregulation and Activation of PKC Alpha by ErbB2 through Src Promotes Breast Cancer Cell Invasion That Can Be Blocked by Combined Treatment with PKC Alpha and Src Inhibitors. Oncogene 2006, 25, 3286–3295. [Google Scholar] [CrossRef] [Green Version]
  10. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the CBioPortal. Sci. Signal. 2013, 6, pl1. [Google Scholar] [CrossRef]
  11. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The CBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Griner, E.M.; Kazanietz, M.G. Protein Kinase C and Other Diacylglycerol Effectors in Cancer. Nat. Rev. Cancer 2007, 7, 281–294. [Google Scholar] [CrossRef] [PubMed]
  13. Kang, J.-H. Protein Kinase C (PKC) Isozymes and Cancer. New J. Sci. 2014, 8, 168–173. [Google Scholar] [CrossRef] [Green Version]
  14. Newton, A.C.; Brognard, J. Reversing the Paradigm: Protein Kinase C as a Tumor Suppressor. Trends Pharmacol. Sci. 2017, 38, 438–447. [Google Scholar] [CrossRef] [Green Version]
  15. Tonetti, D.A.; Chisamore, M.J.; Grdina, W.; Schurz, H.; Jordan, V.C. Stable Transfection of Protein Kinase C Alpha CDNA in Hormone-Dependent Breast Cancer Cell Lines. Br. J. Cancer 2000, 83, 782. [Google Scholar] [CrossRef] [Green Version]
  16. Tonetti, D.A.; Morrow, M.; Kidwai, N.; Gupta, A.; Badve, S. Elevated Protein Kinase C Alpha Expression May Be Predictive of Tamoxifen Treatment Failure. Br. J. Cancer 2003, 88, 1400. [Google Scholar] [CrossRef] [Green Version]
  17. Correction: Novel Antitumor Effect of Estradiol in Athymic Mice Injected with a T47D Breast Cancer Cell Line Overexpressing Protein Kinase C Alpha. Clin. Cancer Res. 2016, 22, 4756. [CrossRef] [Green Version]
  18. Hsu, Y.-H.; Yao, J.; Chan, L.-C.; Wu, T.-J.; Hsu, J.L.; Fang, Y.-F.; Wei, Y.; Wu, Y.; Huang, W.-C.; Liu, C.-L.; et al. Definition of PKC-, CDK6, and MET as Therapeutic Targets in Triple-Negative Breast Cancer. Cancer Res. 2014, 74, 4822–4835. [Google Scholar] [CrossRef] [Green Version]
  19. Tam, W.L.; Lu, H.; Buikhuisen, J.; Soh, B.S.; Lim, E.; Reinhardt, F.; Wu, Z.J.; Krall, J.A.; Bierie, B.; Guo, W.; et al. Protein Kinase C α Is a Central Signaling Node and Therapeutic Target for Breast Cancer Stem Cells. Cancer Cell 2013, 24, 347–364. [Google Scholar] [CrossRef] [Green Version]
  20. Kawano, T.; Inokuchi, J.; Eto, M.; Murata, M.; Kang, J.H. Activators and Inhibitors of Protein Kinase C (PKC): Their Applications in Clinical Trials. Pharmaceutics 2021, 13, 1748. [Google Scholar] [CrossRef]
  21. Kannaiyan, R.; Mahadevan, D. A Comprehensive Review of Protein Kinase Inhibitors for Cancer Therapy. Expert Rev. Anticancer Ther. 2018, 18, 1249–1270. [Google Scholar] [CrossRef] [PubMed]
  22. Cohen, P.; Cross, D.; Jänne, P.A. Kinase Drug Discovery 20 Years after Imatinib: Progress and Future Directions. Nat. Rev. Drug Discov. 2021, 20, 551–569. [Google Scholar] [CrossRef] [PubMed]
  23. Bhullar, K.S.; Lagarón, N.O.; Mcgowan, E.M.; Parmar, I.; Jha, A.; Hubbard, B.P.; Rupasinghe, H.P.V. Kinase-Targeted Cancer Therapies: Progress, Challenges and Future Directions. Mol. Cancer 2018, 17, 48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Yamaoka, T.; Kusumoto, S.; Ando, K.; Ohba, M.; Ohmori, T. Receptor Tyrosine Kinase-Targeted Cancer Therapy. Int. J. Mol. Sci. 2018, 19, 3491. [Google Scholar] [CrossRef] [Green Version]
  25. Alharbi, K.S.; Javed Shaikh, M.A.; Afzal, O.; Alfawaz Altamimi, A.S.; Almalki, W.H.; Alzarea, S.I.; Kazmi, I.; Al-Abbasi, F.A.; Singh, S.K.; Dua, K.; et al. An Overview of Epithelial Growth Factor Receptor (EGFR) Inhibitors in Cancer Therapy. Chem. Biol. Interact. 2022, 366, 110108. [Google Scholar] [CrossRef]
  26. Alzahrani, A.S. PI3K/Akt/MTOR Inhibitors in Cancer: At the Bench and Bedside. Semin. Cancer Biol. 2019, 59, 125–132. [Google Scholar] [CrossRef]
  27. Anwar, S.; Mohammad, T.; Shamsi, A.; Queen, A.; Parveen, S.; Luqman, S.; Hasan, G.M.; Alamry, K.A.; Azum, N.; Asiri, A.M.; et al. Discovery of Hordenine as a Potential Inhibitor of Pyruvate Dehydrogenase Kinase 3: Implication in Lung Cancer Therapy. Biomedicines 2020, 8, 119. [Google Scholar] [CrossRef]
  28. Fassl, A.; Geng, Y.; Sicinski, P. CDK4 and CDK6 Kinases: From Basic Science to Cancer Therapy. Science 2022, 375, eabc1495. [Google Scholar] [CrossRef]
  29. Du, R.; Huang, C.; Liu, K.; Li, X.; Dong, Z. Targeting AURKA in Cancer: Molecular Mechanisms and Opportunities for Cancer Therapy. Mol. Cancer 2021, 20, 1–27. [Google Scholar] [CrossRef]
  30. Yesilkanal, A.E.; Johnson, G.L.; Ramos, A.F.; Rosner, M.R. New Strategies for Targeting Kinase Networks in Cancer. J. Biol. Chem. 2021, 297, 101128. [Google Scholar] [CrossRef]
  31. Moradi Kashkooli, F.; Soltani, M. Evaluation of Solid Tumor Response to Sequential Treatment Cycles via a New Computational Hybrid Approach. Sci. Rep. 2021, 11, 21475. [Google Scholar] [CrossRef] [PubMed]
  32. Zhou, Y.; Zhou, B.; Pache, L.; Chang, M.; Khodabakhshi, A.H.; Tanaseichuk, O.; Benner, C.; Chanda, S.K. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat. Commun. 2019, 10, 1523. [Google Scholar] [CrossRef] [PubMed]
  33. 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]
  34. 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]
  35. Bader, G.D.; Hogue, C.W.V. An Automated Method for Finding Molecular Complexes in Large Protein Interaction Networks. BMC Bioinform. 2003, 4, 2. [Google Scholar] [CrossRef] [Green Version]
  36. Bindea, G.; Mlecnik, B.; Hackl, H.; Charoentong, P.; Tosolini, M.; Kirilovsky, A.; Fridman, W.H.; Pagès, F.; Trajanoski, Z.; Galon, J. ClueGO: A Cytoscape Plug-in to Decipher Functionally Grouped Gene Ontology and Pathway Annotation Networks. Bioinformatics 2009, 25, 1091–1093. [Google Scholar] [CrossRef] [Green Version]
  37. Burgdorf, C.; Schäfer, U.; Richardt, G.; Kurz, T. U73122, an Aminosteroid Phospholipase C Inhibitor, Is a Potent Inhibitor of Cardiac Phospholipase D by a PIP2-Dependent Mechanism. J. Cardiovasc. Pharmacol. 2010, 55, 555–559. [Google Scholar] [CrossRef]
  38. Parker, P.J.; Brown, S.J.; Calleja, V.; Chakravarty, P.; Cobbaut, M.; Linch, M.; Marshall, J.J.T.; Martini, S.; McDonald, N.Q.; Soliman, T.; et al. Equivocal, Explicit and Emergent Actions of PKC Isoforms in Cancer. Nat. Rev. Cancer 2021, 21, 51–63. [Google Scholar] [CrossRef]
  39. Black, A.R.; Black, J.D. The Complexities of PKCα Signaling in Cancer. Adv. Biol. Regul. 2021, 80, 100769. [Google Scholar] [CrossRef]
  40. Isakov, N. Protein Kinase C (PKC) Isoforms in Cancer, Tumor Promotion and Tumor Suppression. Semin. Cancer Biol. 2018, 48, 36–52. [Google Scholar] [CrossRef]
  41. Caiazza, F.; Elliott, L.; Fennelly, D.; Sheahan, K.; Doherty, G.A.; Ryan, E.J. Targeting EGFR in Metastatic Colorectal Cancer beyond the Limitations of KRAS Status: Alternative Biomarkers and Therapeutic Strategies. Biomark. Med. 2015, 9, 363–375. [Google Scholar] [CrossRef] [PubMed]
  42. Caiazza, F.; Harvey, B.J.; Thomas, W. Cytosolic Phospholipase A 2 Activation Correlates with HER2 Overexpression and Mediates Estrogen-Dependent Breast Cancer Cell Growth. Mol. Endocrinol. 2010, 24, 953–968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Cubas, R.; Zhang, S.; Li, M.; Chen, C.; Yao, Q. Trop2 Expression Contributes to Tumor Pathogenesis by Activating the ERK MAPK Pathway. Mol. Cancer 2010, 9, 253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Wieczorek, M.; Abualrous, E.T.; Sticht, J.; Álvaro-Benito, M.; Stolzenberg, S.; Noé, F.; Freund, C. Major Histocompatibility Complex (MHC) Class I and MHC Class II Proteins: Conformational Plasticity in Antigen Presentation. Front. Immunol. 2017, 8, 292. [Google Scholar] [CrossRef] [Green Version]
  45. Pishesha, N.; Harmand, T.J.; Ploegh, H.L. A Guide to Antigen Processing and Presentation. Nat. Rev. Immunol. 2022, 26, 1–14. [Google Scholar] [CrossRef]
  46. Dersh, D.; Hollý, J.; Yewdell, J.W. A Few Good Peptides: MHC Class I-Based Cancer Immunosurveillance and Immunoevasion. Nat. Rev. Immunol. 2021, 21, 116–128. [Google Scholar] [CrossRef]
  47. Zhao, J.; Zhang, X.; Cheng, M.; Guan, J.; Gai, J.; Fu, L.; Zhang, R.; Du, T.; Li, Q. Expression of IFN-Induced 2′-5′-Oligoadenylate Synthetases Correlates with Immune Infiltration, Revealing Potential Targets and New Biomarkers for Basal-like Breast Cancer Prognosis. Int. Immunopharmacol. 2020, 88, 106916. [Google Scholar] [CrossRef]
  48. Zhang, Y.; Yu, C. Prognostic Characterization of OAS1/OAS2/OAS3/OASL in Breast Cancer. BMC Cancer 2020, 20, 575. [Google Scholar] [CrossRef]
  49. Cooper, J.; Giancotti, F.G. Cancer Cell Review Integrin Signaling in Cancer: Mechanotransduction, Stemness, Epithelial Plasticity, and Therapeutic Resistance. Cancer Cell 2019, 35, 347–367. [Google Scholar] [CrossRef]
  50. Karantza, V. Keratins in Health and Cancer: More than Mere Epithelial Cell Markers. Oncogene 2010, 30, 127–138. [Google Scholar] [CrossRef]
  51. Heo, S.H.; Lee, J.Y.; Yang, K.M.; Park, K.S. ELK3 Expression Correlates with Cell Migration, Invasion, and Membrane Type 1-Matrix Metalloproteinase Expression in MDA-MB-231 Breast Cancer Cells. Gene Expr. 2015, 16, 197–203. [Google Scholar] [CrossRef] [PubMed]
  52. Hoxhaj, G.; Manning, B.D. The PI3K–AKT Network at the Interface of Oncogenic Signalling and Cancer Metabolism. Nat. Rev. Cancer 2020, 20, 74–88. [Google Scholar] [CrossRef] [PubMed]
  53. Miricescu, D.; Totan, A.; Stanescu-Spinu, I.I.; Badoiu, S.C.; Stefani, C.; Greabu, M. PI3K/AKT/MTOR Signaling Pathway in Breast Cancer: From Molecular Landscape to Clinical Aspects. Int. J. Mol. Sci. 2021, 22, 173. [Google Scholar] [CrossRef] [PubMed]
  54. Costa, R.L.B.; Han, H.S.; Gradishar, W.J. Targeting the PI3K/AKT/MTOR Pathway in Triple-Negative Breast Cancer: A Review. Breast Cancer Res. Treat. 2018, 169, 397–406. [Google Scholar] [CrossRef]
  55. Thorpe, L.M.; Yuzugullu, H.; Zhao, J.J. PI3K in Cancer: Divergent Roles of Isoforms, Modes of Activation and Therapeutic Targeting. Nat. Rev. Cancer 2015, 15, 7–24. [Google Scholar] [CrossRef] [Green Version]
  56. Rathinaswamy, M.K.; Burke, J.E. Advances in Biological Regulation Class I Phosphoinositide 3-Kinase (PI3K) Regulatory Subunits and Their Roles in Signaling and Disease. Adv. Biol. Regul. 2020, 75, 100657. [Google Scholar] [CrossRef]
  57. Taniguchi, C.M.; Winnay, J.; Kondo, T.; Bronson, R.T.; Guimaraes, A.R.; Alemán, J.O.; Luo, J.; Stephanopoulos, G.; Weissleder, R.; Cantley, L.C.; et al. The Phosphoinositide 3-Kinase Regulatory Subunit P85α Can Exert Tumor Suppressor Properties through Negative Regulation of Growth Factor Signaling. Cancer Res. 2010, 70, 5305–5315. [Google Scholar] [CrossRef] [Green Version]
  58. Cizkova, M.; Vacher, S.; Meseure, D.; Trassard, M.; Susini, A.; Mlcuchova, D.; Callens, C.; Rouleau, E.; Spyratos, F.; Lidereau, R.; et al. PIK3R1 Underexpression Is an Independent Prognostic Marker in Breast Cancer. BMC Cancer 2013, 13, 545. [Google Scholar] [CrossRef] [Green Version]
  59. Cortés, I.; Sánchez-Ruíz, J.; Zuluaga, S.; Calvanese, V.; Marqués, M.; Hernández, C.; Rivera, T.; Kremer, L.; González-García, A.; Carrera, A.C. P85β Phosphoinositide 3-Kinase Subunit Regulates Tumor Progression. Proc. Natl. Acad. Sci. USA 2012, 109, 11318–11323. [Google Scholar] [CrossRef] [Green Version]
  60. Wang, G.; Liu, G.; Wang, X.; Sethi, S.; Ali-Fehmi, R.; Abrams, J.; Zheng, Z.; Zhang, K.; Ethier, S.; Yang, Z.Q. ERLIN2 Promotes Breast Cancer Cell Survival by Modulating Endoplasmic Reticulum Stress Pathways. BMC Cancer 2012, 12, 225. [Google Scholar] [CrossRef] [Green Version]
  61. Ivashkiv, L.B. IFNγ: Signalling, Epigenetics and Roles in Immunity, Metabolism, Disease and Cancer Immunotherapy. Nat. Rev. Immunol. 2018, 18, 545–558. [Google Scholar] [CrossRef] [PubMed]
  62. Barrat, F.J.; Crow, M.K.; Ivashkiv, L.B. Interferon Target-Gene Expression and Epigenomic Signatures in Health and Disease. Nat. Immunol. 2019, 20, 1574–1583. [Google Scholar] [CrossRef] [PubMed]
  63. Van der Vorst, E.P.C.; Theodorou, K.; Wu, Y.; Hoeksema, M.A.; Goossens, P.; Bursill, C.A.; Aliyev, T.; Huitema, L.F.A.; Tas, S.W.; Wolfs, I.M.J.; et al. High-Density Lipoproteins Exert Pro-Inflammatory Effects on Macrophages via Passive Cholesterol Depletion and PKC-NF-ΚB/STAT1-IRF1 Signaling. Cell Metab. 2017, 25, 197–207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Zheng, X.; Li, A.S.; Zheng, H.; Zhao, D.; Guan, D.; Zou, H. Different Associations of CD45 Isoforms with STAT3, PKC and ERK Regulate IL-6-Induced Proliferation in Myeloma. PLoS ONE 2015, 10, e0119780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Verhoeven, Y.; Tilborghs, S.; Jacobs, J.; De Waele, J.; Quatannens, D.; Deben, C.; Prenen, H.; Pauwels, P.; Trinh, X.B.; Wouters, A.; et al. The Potential and Controversy of Targeting STAT Family Members in Cancer. Semin. Cancer Biol. 2020, 60, 41–56. [Google Scholar] [CrossRef]
  66. Loh, C.Y.; Arya, A.; Naema, A.F.; Wong, W.F.; Sethi, G.; Looi, C.Y. Signal Transducer and Activator of Transcription (STATs) Proteins in Cancer and Inflammation: Functions and Therapeutic Implication. Front. Oncol. 2019, 9, 48. [Google Scholar] [CrossRef] [Green Version]
  67. Rah, B.; Rather, R.A.; Bhat, G.R.; Baba, A.B.; Mushtaq, I.; Farooq, M.; Yousuf, T.; Dar, S.B.; Parveen, S.; Hassan, R.; et al. JAK/STAT Signaling: Molecular Targets, Therapeutic Opportunities, and Limitations of Targeted Inhibitions in Solid Malignancies. Front. Pharmacol. 2022, 13, 821344. [Google Scholar] [CrossRef] [PubMed]
  68. Dey, A.; Kundu, M.; Das, S.; Jena, B.C.; Mandal, M. Understanding the Function and Regulation of Sox2 for Its Therapeutic Potential in Breast Cancer. Biochim. Biophys. Acta Rev. Cancer 2022, 1877, 188692. [Google Scholar] [CrossRef] [PubMed]
  69. Alves, C.L.; Elias, D.; Lyng, M.B.; Bak, M.; Ditzel, H.J. SNAI2 Upregulation Is Associated with an Aggressive Phenotype in Fulvestrant-Resistant Breast Cancer Cells and Is an Indicator of Poor Response to Endocrine Therapy in Estrogen Receptor-Positive Metastatic Breast Cancer. Breast Cancer Res. 2018, 20, 60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Obr, A.E.; Kumar, S.; Chang, Y.J.; Bulatowicz, J.J.; Barnes, B.J.; Birge, R.B.; Lazzarino, D.A.; Gallagher, E.; Leroith, D.; Wood, T.L. Insulin-like Growth Factor Receptor Signaling in Breast Tumor Epithelium Protects Cells from Endoplasmic Reticulum Stress and Regulates the Tumor Microenvironment 11 Medical and Health Sciences 1112 Oncology and Carcinogenesis. Breast Cancer Res. 2018, 20, 138. [Google Scholar] [CrossRef]
  71. Opyrchal, M.; Gil, M.; Salisbury, J.L.; Goetz, M.P.; Suman, V.; Degnim, A.; McCubrey, J.; Haddad, T.; Iankov, I.; Kurokawa, C.B.; et al. Molecular Targeting of the Aurora-A/SMAD5 Oncogenic Axis Restores Chemosensitivity in Human Breast Cancer Cells. Oncotarget 2017, 8, 91803–91816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Lavoie, H.; Gagnon, J.; Therrien, M. ERK Signalling: A Master Regulator of Cell Behaviour, Life and Fate. Nat. Rev. Mol. Cell Biol. 2020, 21, 607–632. [Google Scholar] [CrossRef] [PubMed]
  73. Kyriakis, J.M.; Avruch, J. MAP Kinase Pathways. In Compendium of Inflammatory Diseases; Springer: Berlin/Heidelberg, Germany, 2016; pp. 892–908. [Google Scholar] [CrossRef]
  74. Khojasteh Poor, F.; Keivan, M.; Ramazii, M.; Ghaedrahmati, F.; Anbiyaiee, A.; Panahandeh, S.; Khoshnam, S.E.; Farzaneh, M. Mini Review: The FDA-Approved Prescription Drugs That Target the MAPK Signaling Pathway in Women with Breast Cancer. Breast Dis. 2021, 40, 51–62. [Google Scholar] [CrossRef] [PubMed]
  75. Ahmad, D.A.J.; Negm, O.H.; Alabdullah, M.L.; Mirza, S.; Hamed, M.R.; Band, V.; Green, A.R.; Ellis, I.O.; Rakha, E.A. Clinicopathological and Prognostic Significance of Mitogen-Activated Protein Kinases (MAPK) in Breast Cancers. Breast Cancer Res. Treat. 2016, 159, 457–467. [Google Scholar] [CrossRef]
  76. Lin, R.; Bao, X.; Wang, H.; Zhu, S.; Liu, Z.; Chen, Q.; Ai, K.; Shi, B. TRPM2 Promotes Pancreatic Cancer by PKC/MAPK Pathway. Cell Death Dis. 2021, 12, 585. [Google Scholar] [CrossRef]
  77. Kim, S.J.; Pham, T.H.; Bak, Y.; Ryu, H.W.; Oh, S.R.; Yoon, D.Y. Orientin Inhibits Invasion by Suppressing MMP-9 and IL-8 Expression via the PKCα/ ERK/AP-1/STAT3-Mediated Signaling Pathways in TPA-Treated MCF-7 Breast Cancer Cells. Phytomedicine 2018, 50, 35–42. [Google Scholar] [CrossRef]
  78. Huang, C.; Cao, Z.; Ma, J.; Shen, Y.; Bu, Y.; Khoshaba, R.; Shi, G.; Huang, D.; Liao, D.F.; Ji, H.; et al. AKR1B10 Activates Diacylglycerol (DAG) Second Messenger in Breast Cancer Cells. Mol. Carcinog. 2018, 57, 1300–1310. [Google Scholar] [CrossRef]
  79. McCubrey, J.A.; Steelman, L.S.; Chappell, W.H.; Abrams, S.L.; Wong, E.W.T.; Chang, F.; Lehmann, B.; Terrian, D.M.; Milella, M.; Tafuri, A.; et al. Roles of the Raf/MEK/ERK Pathway in Cell Growth, Malignant Transformation and Drug Resistance. Biochim. Biophys. Acta-Mol. Cell Res. 2007, 1773, 1263–1284. [Google Scholar] [CrossRef] [Green Version]
  80. Kruiswijk, F.; Labuschagne, C.F.; Vousden, K.H. P53 in Survival, Death and Metabolic Health: A Lifeguard with a Licence to Kill. Nat. Rev. Mol. Cell Biol. 2015, 16, 393–405. [Google Scholar] [CrossRef]
  81. Liu, J.; Zhang, C.; Hu, W.; Feng, Z. Tumor Suppressor P53 and Metabolism. J. Mol. Cell Biol. 2019, 11, 284–292. [Google Scholar] [CrossRef] [Green Version]
  82. Roos, W.P.; Thomas, A.D.; Kaina, B. DNA Damage and the Balance between Survival and Death in Cancer Biology. Nat. Rev. Cancer 2016, 16, 20–33. [Google Scholar] [CrossRef] [PubMed]
  83. Boehme, K.A.; Kulikov, R.; Blattner, C. P53 Stabilization in Response to DNA Damage Requires Akt/PKB and DNA-PK. Proc. Natl. Acad. Sci. USA 2008, 105, 7785–7790. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Wu, T.T.; Hsieh, Y.H.; Hsieh, Y.S.; Liu, J.Y. Reduction of PKCα Decreases Cell Proliferation, Migration, and Invasion of Human Malignant Hepatocellular Carcinoma. J. Cell. Biochem. 2008, 103, 9–20. [Google Scholar] [CrossRef] [PubMed]
  85. Smith, S.D.; Enge, M.; Bao, W.; Thullberg, M.; Costa, T.D.F.; Olofsson, H.; Gashi, B.; Selivanova, G.; Strömblad, S. Protein Kinase Cα (PKCα) Regulates P53 Localization and Melanoma Cell Survival Downstream of Integrin Av in Three-Dimensional Collagen and in Vivo. J. Biol. Chem. 2012, 287, 29336–29347. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Brufsky, A.M.; Dickler, M.N. Estrogen Receptor-Positive Breast Cancer: Exploiting Signaling Pathways Implicated in Endocrine Resistance. Oncologist 2018, 23, 528–539. [Google Scholar] [CrossRef] [Green Version]
  87. Lei, J.T.; Gou, X.; Seker, S.; Ellis, M.J. ESR1 Alterations and Metastasis in Estrogen Receptor Positive Breast Cancer. J. Cancer Metastasis Treat. 2019, 5, 38. [Google Scholar] [CrossRef] [Green Version]
  88. Toska, E.; Osmanbeyoglu, H.U.; Castel, P.; Chan, C.; Hendrickson, R.C.; Elkabets, M.; Dickler, M.N.; Scaltriti, M.; Leslie, C.S.; Armstrong, S.A.; et al. PI3K Pathway Regulates ER-Dependent Transcription in Breast Cancer through the Epigenetic Regulator KMT2D. Science 2017, 355, 1324–1330. [Google Scholar] [CrossRef] [Green Version]
  89. Bosch, A.; Li, Z.; Bergamaschi, A.; Ellis, H.; Toska, E.; Prat, A.; Tao, J.J.; Spratt, D.E.; Viola-Villegas, N.T.; Castel, P.; et al. PI3K Inhibition Results in Enhanced Estrogen Receptor Function and Dependence in Hormone Receptor-Positive Breast Cancer. Sci. Transl. Med. 2015, 7, 283ra51. [Google Scholar] [CrossRef] [Green Version]
  90. Guardiola, S.; Varese, M.; Sánchez-Navarro, M.; Giralt, E. A Third Shot at EGFR: New Opportunities in Cancer Therapy. Trends Pharmacol. Sci. 2019, 40, 941–955. [Google Scholar] [CrossRef]
  91. Harris, R.C.; Chung, E.; Coffey, R.J. EGF Receptor Ligands. Exp. Cell Res. 2003, 284, 2–13. [Google Scholar] [CrossRef]
  92. Sigismund, S.; Avanzato, D.; Lanzetti, L. Emerging Functions of the EGFR in Cancer. Mol. Oncol. 2018, 12, 3–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. El-Gamal, M.I.; Mewafi, N.H.; Abdelmotteleb, N.E.; Emara, M.A.; Tarazi, H.; Sbenati, R.M.; Madkour, M.M.; Zaraei, S.O.; Shahin, A.I.; Anbar, H.S. A Review of HER4 (ErbB4) Kinase, Its Impact on Cancer, and Its Inhibitors. Molecules 2021, 26, 7376. [Google Scholar] [CrossRef] [PubMed]
  94. Segers, V.F.M.; Dugaucquier, L.; Feyen, E.; Shakeri, H.; De Keulenaer, G.W. The Role of ErbB4 in Cancer. Cell. Oncol. 2020, 43, 335–352. [Google Scholar] [CrossRef]
  95. Lucas, L.M.; Dwivedi, V.; Senfeld, J.I.; Cullum, R.L.; Mill, C.P.; Piazza, J.T.; Bryant, I.N.; Cook, L.J.; Miller, S.T.; Lott, J.H.; et al. The Yin and Yang of ERBB4: Tumor Suppressor and Oncoprotein. Pharmacol. Rev. 2022, 74, 18–47. [Google Scholar] [CrossRef] [PubMed]
  96. Avraham, R.; Yarden, Y. Feedback Regulation of EGFR Signalling: Decision Making by Early and Delayed Loops. Nat. Rev. Mol. Cell Biol. 2011, 12, 104–117. [Google Scholar] [CrossRef]
  97. Fernández-Ulibarri, I.; Vilella, M.; Lázaro-Diéguez, F.; Sarri, E.; Martínez, S.E.; Jiménez, N.; Claro, E.; Mérida, I.; Burger, K.N.J.; Egea, G. Diacylglycerol Is Required for the Formation of COPI Vesicles in the Golgi-to-ER Transport Pathway. Mol. Biol. Cell 2007, 18, 3250–3263. [Google Scholar] [CrossRef] [Green Version]
  98. Lonic, A.; Gehling, F.; Belle, L.; Li, X.; Schieber, N.L.; Nguyen, E.V.; Goodall, G.J.; Parton, R.G.; Daly, R.J.; Khew-Goodall, Y. Phosphorylation of PKCδ by FER Tips the Balance from EGFR Degradation to Recycling. J. Cell Biol. 2021, 220, e201902073. [Google Scholar] [CrossRef]
  99. Liu, M.; Idkowiak-Baldys, J.; Roddy, P.L.; Baldys, A.; Raymond, J.; Clarke, C.J.; Hannun, Y.A. Sustained Activation of Protein Kinase C Induces Delayed Phosphorylation and Regulates the Fate of Epidermal Growth Factor Receptor. PLoS ONE 2013, 8, e80721. [Google Scholar] [CrossRef] [Green Version]
  100. Heckman, C.A.; Biswas, T.; Dimick, D.M.; Cayer, M.L. Activated Protein Kinase C (PKC) Is Persistently Trafficked with Epidermal Growth Factor (EGF) Receptor. Biomolecules 2020, 10, 1288. [Google Scholar] [CrossRef]
  101. Bao, J.; Alroy, I.; Waterman, H.; Schejter, E.D.; Brodie, C.; Gruenberg, J.; Yarden, Y. Threonine Phosphorylation Diverts Internalized Epidermal Growth Factor Receptors from a Degradative Pathway to the Recycling Endosome. J. Biol. Chem. 2000, 275, 26178–26186. [Google Scholar] [CrossRef] [Green Version]
  102. Dietrich, M.; Malik, M.S.; Skeie, M.; Bertelsen, V.; Stang, E. Protein Kinase C Regulates ErbB3 Turnover. Exp. Cell Res. 2019, 382, 111473. [Google Scholar] [CrossRef] [PubMed]
  103. Basso, F.; Rocchetti, F.; Rodriguez, S.; Nesterova, M.; Cormier, F.; Stratakis, C.A.; Ragazzon, B.; Bertherat, J.; Rizk-Rabin, M. Comparison of the Effects of PRKAR1A and PRKAR2B Depletion on Signaling Pathways, Cell Growth, and Cell Cycle Control of Adrenocortical Cells. Horm. Metab. Res. 2014, 46, 883. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Saidy, B.; Kotecha, S.; Butler, A.; Rakha, E.A.; Ellis, I.O.; Green, A.R.; Martin, S.G.; Storr, S.J. PP1, PKA and DARPP-32 in Breast Cancer: A Retrospective Assessment of Protein and MRNA Expression. J. Cell. Mol. Med. 2021, 25, 5015. [Google Scholar] [CrossRef] [PubMed]
  105. Li, C.-F.; Liu, T.-T.; Chuang, I.-C.; Chen, Y.-Y.; Fang, F.-M.; Chan, T.-C.; Li, W.-S.; Huang, H.-Y. PLCB4 Copy Gain and PLCß4 Overexpression in Primary Gastrointestinal Stromal Tumors: Integrative Characterization of a Lipid-Catabolizing Enzyme Associated with Worse Disease-Free Survival Network; PLCB4, Phospholipase C-SS4. Oncotarget 2017, 8, 19997–20010. [Google Scholar] [CrossRef] [Green Version]
  106. Han, N.; Zhang, Y.Y.; Zhang, Z.M.; Zhang, F.; Zeng, T.Y.; Zhang, Y.B.; Zhao, W.C. High Expression of PDGFA Predicts Poor Prognosis of Esophageal Squamous Cell Carcinoma. Medicine 2021, 100, e25932. [Google Scholar] [CrossRef]
  107. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 Years of Image Analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef]
  108. Ho, J.; Tumkaya, T.; Aryal, S.; Choi, H.; Claridge-Chang, A. Moving beyond P Values: Data Analysis with Estimation Graphics. Nat. Methods 2019, 16, 565–566. [Google Scholar] [CrossRef]
Figure 1. Knock-down of PKCα in MCF-7 cells affects their proliferation capacity: (A) MCF-7 cells were transfected with siRNA specific for PKCα. Protein expression was assessed by WB relative to β-actin expression at the days indicated. (B) Proliferation rates of PKCα-WT and PKCα-KD MCF-7 cells were measured by fluorometric DNA quantitation. The mean difference for two comparisons against the PKCα-WT and PKCα-KD cells are shown in the Cumming estimation plots (see Materials and Methods for a more extensive explanation).
Figure 1. Knock-down of PKCα in MCF-7 cells affects their proliferation capacity: (A) MCF-7 cells were transfected with siRNA specific for PKCα. Protein expression was assessed by WB relative to β-actin expression at the days indicated. (B) Proliferation rates of PKCα-WT and PKCα-KD MCF-7 cells were measured by fluorometric DNA quantitation. The mean difference for two comparisons against the PKCα-WT and PKCα-KD cells are shown in the Cumming estimation plots (see Materials and Methods for a more extensive explanation).
Ijms 23 14023 g001
Figure 2. Transcriptomic analysis: (A) Volcano plot displaying the expression of probe sets in PKCα-KD vs. PKCα-WT MCF-7 cells. Log2 fold changes and their corresponding p-values of all genes in the microarray were taken to construct the volcano plot. Genes up-regulated with more than 1.5-fold change with a p-value <0.05 are depicted in blue dots, and those down-regulated with identical fold change and p-value are in red dots. All other genes in the array are not significantly altered and are represented in grey dots. (B) Validation of significant up-and down-regulated genes by qPCR.
Figure 2. Transcriptomic analysis: (A) Volcano plot displaying the expression of probe sets in PKCα-KD vs. PKCα-WT MCF-7 cells. Log2 fold changes and their corresponding p-values of all genes in the microarray were taken to construct the volcano plot. Genes up-regulated with more than 1.5-fold change with a p-value <0.05 are depicted in blue dots, and those down-regulated with identical fold change and p-value are in red dots. All other genes in the array are not significantly altered and are represented in grey dots. (B) Validation of significant up-and down-regulated genes by qPCR.
Ijms 23 14023 g002
Figure 3. Functional enrichment analysis with Metascape of DEGs obtained after comparing PKCα-KD vs. PKCα-WT MCF-7 cells. Network of enriched terms colored by cluster-ID. DEGs on the left represent the up-regulated genes and on the right represent the down-regulated genes. Here, we show the five top clusters with enriched representative terms, and the numbers in brackets indicate the (Log10 (p-value)). Table S2a,b show the genes included in each cluster.
Figure 3. Functional enrichment analysis with Metascape of DEGs obtained after comparing PKCα-KD vs. PKCα-WT MCF-7 cells. Network of enriched terms colored by cluster-ID. DEGs on the left represent the up-regulated genes and on the right represent the down-regulated genes. Here, we show the five top clusters with enriched representative terms, and the numbers in brackets indicate the (Log10 (p-value)). Table S2a,b show the genes included in each cluster.
Ijms 23 14023 g003
Figure 4. Human phospho-kinase array comparing PKCα-KD vs. PKCα-WT MCF-7 cells. The histogram shows the fold changes of the kinases affected by the down-regulation of the expression of PKCα in MCF-7 cells. Red bars represent phosphorylated proteins showing a decrease in phosphorylation levels, whereas blue bars show an increase.
Figure 4. Human phospho-kinase array comparing PKCα-KD vs. PKCα-WT MCF-7 cells. The histogram shows the fold changes of the kinases affected by the down-regulation of the expression of PKCα in MCF-7 cells. Red bars represent phosphorylated proteins showing a decrease in phosphorylation levels, whereas blue bars show an increase.
Ijms 23 14023 g004
Figure 5. Critical signaling pathway analysis: (A) Protein–protein interaction network of DEGs using the STRING database [33]. DEGs obtained from the microarray analysis were colored in dark red when down-expressed, whereas up-expressed DEGs were colored in dark blue. Down- and up-activated kinases obtained from the phospho-kinase array were depicted in light red and light blue, respectively. The MCODE tool was used to identify the two main protein clusters present in the protein-protein interaction network: Cluster 1 (22 nodes and 222 edges) (B) and Cluster 2 (17 nodes and 51 edges) (C).
Figure 5. Critical signaling pathway analysis: (A) Protein–protein interaction network of DEGs using the STRING database [33]. DEGs obtained from the microarray analysis were colored in dark red when down-expressed, whereas up-expressed DEGs were colored in dark blue. Down- and up-activated kinases obtained from the phospho-kinase array were depicted in light red and light blue, respectively. The MCODE tool was used to identify the two main protein clusters present in the protein-protein interaction network: Cluster 1 (22 nodes and 222 edges) (B) and Cluster 2 (17 nodes and 51 edges) (C).
Ijms 23 14023 g005
Figure 6. Effect of additional inhibition of the key signaling pathways identified in the arrays. The proliferation and migration capabilities of PKCα-WT and PKCα-KD MCF-7 cells were measured under different treatments: (A,B), U73122; (C,D), KT5720; (E,F), BMS 599626; (G,H), Imatinib on PKCα-WT and PKCα-KD MCF-7 cells. A heat map-type graph was used with the dark blue color corresponding to the lowest value for the DNA fold change or migration % and the dark red color to the highest value. The effect of each treatment on PKCα-KD vs. PKCα-WT MCF-7 cells at day 8 (proliferation assay) or 72 h (migration assay) is shown at the top of the heat map as the percentage of inhibition on the ability to proliferate or migrate.
Figure 6. Effect of additional inhibition of the key signaling pathways identified in the arrays. The proliferation and migration capabilities of PKCα-WT and PKCα-KD MCF-7 cells were measured under different treatments: (A,B), U73122; (C,D), KT5720; (E,F), BMS 599626; (G,H), Imatinib on PKCα-WT and PKCα-KD MCF-7 cells. A heat map-type graph was used with the dark blue color corresponding to the lowest value for the DNA fold change or migration % and the dark red color to the highest value. The effect of each treatment on PKCα-KD vs. PKCα-WT MCF-7 cells at day 8 (proliferation assay) or 72 h (migration assay) is shown at the top of the heat map as the percentage of inhibition on the ability to proliferate or migrate.
Ijms 23 14023 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Serrano-López, E.M.; Coronado-Parra, T.; Marín-Vicente, C.; Szallasi, Z.; Gómez-Abellán, V.; López-Andreo, M.J.; Gragera, M.; Gómez-Fernández, J.C.; López-Nicolás, R.; Corbalán-García, S. Deciphering the Role and Signaling Pathways of PKCα in Luminal A Breast Cancer Cells. Int. J. Mol. Sci. 2022, 23, 14023. https://doi.org/10.3390/ijms232214023

AMA Style

Serrano-López EM, Coronado-Parra T, Marín-Vicente C, Szallasi Z, Gómez-Abellán V, López-Andreo MJ, Gragera M, Gómez-Fernández JC, López-Nicolás R, Corbalán-García S. Deciphering the Role and Signaling Pathways of PKCα in Luminal A Breast Cancer Cells. International Journal of Molecular Sciences. 2022; 23(22):14023. https://doi.org/10.3390/ijms232214023

Chicago/Turabian Style

Serrano-López, Emilio M., Teresa Coronado-Parra, Consuelo Marín-Vicente, Zoltan Szallasi, Victoria Gómez-Abellán, María José López-Andreo, Marcos Gragera, Juan C. Gómez-Fernández, Rubén López-Nicolás, and Senena Corbalán-García. 2022. "Deciphering the Role and Signaling Pathways of PKCα in Luminal A Breast Cancer Cells" International Journal of Molecular Sciences 23, no. 22: 14023. https://doi.org/10.3390/ijms232214023

APA Style

Serrano-López, E. M., Coronado-Parra, T., Marín-Vicente, C., Szallasi, Z., Gómez-Abellán, V., López-Andreo, M. J., Gragera, M., Gómez-Fernández, J. C., López-Nicolás, R., & Corbalán-García, S. (2022). Deciphering the Role and Signaling Pathways of PKCα in Luminal A Breast Cancer Cells. International Journal of Molecular Sciences, 23(22), 14023. https://doi.org/10.3390/ijms232214023

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