- Open Access
Beta 2 adrenergic receptor and mu opioid receptor interact to potentiate the aggressiveness of human breast cancer cell by activating the glycogen synthase kinase 3 signaling
Breast Cancer Research volume 24, Article number: 33 (2022)
Opioid and beta-adrenergic receptors are recently shown to cross talk via formation of receptor heterodimers to control the growth and proliferation of breast cancer cells. However, the underlying cell signaling mechanism remained unclear.
To determine the effect of the interaction of the two systems in breast cancer, we employed triple-negative breast cancer cell lines MDA-MB-231 and MDA-MB-468, CRISPR or chemical inhibition or activation of beta-adrenergic receptors (B2AR) and mu-opioid receptors (MOR) gene, and PCR array technology and studied aggressive tumor phenotype and signaling cascades.
We show here that in triple-negative breast cancer cells, the reduction in expression B2AR and MOR by genetic and pharmacological tools leads to a less aggressive phenotype of triple-negative breast cancer cells in vitro and in animal xenografts. Genomic analysis indicates the glycogen synthase kinase 3 (GSK3) pathway as a possible candidate messenger system involved in B2AR and MOR cross talk. GSK3 inactivation in MDA-MB-231 and MDA-MB-468 cells induced similar phenotypic changes as the inhibition of B2AR and/or MOR, while a GSK3 activation by wortmannin reversed the effects of B2AR and/or MOR knockdown on these cells. GSK3 inactivation also prevents B2AR agonist norepinephrine or MOR agonist DAMGO from affecting MDA-MB-231 and MDA-MB-468 cell proliferation.
These data confirm a role of B2AR and MOR interaction in the control of breast cancer cell growth and identify a possible role of the GSK3 signaling system in mediation of these two receptors’ cross talk. Screening for ligands targeting B2AR and MOR interaction and/or the GSK3 system may help to identify novel drugs for the prevention of triple-negative breast cancer cell growth and metastasis.
Breast cancer is the world’s leading cause of cancer mortality in women in the world [1,2,3,4,5]. Research has demonstrated a strong link between chronic stress and a weakened immune system, leaving the body prone to disease, like cancer [6,7,8,9]. Stress hormones, such as norepinephrine and epinephrine, acting primarily via β2-adrenergic receptors (B2AR), have been shown to promote tumor growth and differentiation, angiogenesis, migration, and metastasis [10,11,12]. Additionally, B2AR is highly expressed in different types of cancer, including breast cancer, and is associated with a poor prognosis [13, 14]. Also, recently B2AR blockers have been shown to inhibit metastatic spread and reduce cancer risk , providing clinical evidence for the involvement of B2AR in cancer progression. Body stress also affects the production of endogenous opioid peptides, which are critical to control stress homeostasis [16,17,18]. Opioid receptors are also found to be upregulated in several types of cancer [19, 20]. Administration of morphine, an opioid peptide which binds to μ-opioid receptors (MOR), has been shown to promote tumor growth, tumor apoptosis, angiogenesis, migration, and invasiveness [21,22,23,24], while a MOR antagonist, methylnaltrexone, has been shown to possess antiproliferative and antimetastatic properties in lung cancer . Hence, stress promotion of cancer might involve both B2AR and MOR.
Both MOR and B2AR are G-protein-coupled receptors (GPCRs) and are closely related receptor systems and coexist in many cells, including breast tumor cells [10, 11, 26]. Within the breast cancer cell lines, these receptors are present in both triple-negative breast cancer (TNBC) lines like MDA-MB-231 and MDA-MB-468 cells that represent the basal B and basal A subtypes of breast cancer, and non-TNBC lines like MCF-7 and T47D cells that represent the luminal A subtype of breast cancer [27,28,29]. GPCRs are known to heterodimerize with closely related members, resulting in the modulation of their functions [12, 13, 15]. Ligand binding studies have shown that opioidergic and adrenergic agents can bind to their complementary receptors and enhance the activity of other compound types through an allosteric mechanism . Physical interaction of B2AR and MOR has also been reported by co-immunoprecipitation studies in breast cancer cells , in CHO cells exogenously expressing both receptors , and in neuronal cells . These data support the notion that opioidergic and adrenergic systems may interact to modulate cellular processes. However, opioidergic and adrenergic systems interaction in breast cancer cells remain unclear. We recently have shown that propranolol (PRO), a beta blocker and naltrexone (NTX), an opiate antagonist, when combined, produce marked inhibitory effects on breast tumor growth and tumor mass while improving the survival rate . These antitumor effects resulted from a reduction in tumor cell proliferation, the induction of cellular apoptosis, and the prevention of epithelial–mesenchymal transition in the tumor. However, the underlying cell signaling mechanism remained unclear. We demonstrate in this study that B2AR and MOR interact to control triple-negative breast cancer (TNBC) cell growth and progression via activation of the GSK3 signaling system.
Materials and methods
The animal care and treatment were performed in accordance with institutional guidelines, and protocols were approved by the Rutgers Institutional Animal Care and Facilities Committee (PROTO 999900286) and complied with National Institutes of Health policy.
MDA-MB-231 and MDA-MB-468 cell line
Breast cancer cell line MDA-MB-231 and MDA-MB-468 was obtained from the ATCC® (reference HTB-26™ and HTB-132™; Gaithersburg, MD). The cell lines were cultivated in DMEM high glucose medium (Gibco; Gaithersburg, MD) with a final concentration of 10% fetal bovine serum (FBS) (Gibco) and 1% antibiotic–antimycotic (Gibco). Cells were trypsinized every 3 to 4 days, never exceeding 10 passages in all experiments, and maintained at 37 °C and 5% CO2.
The specific design of forward and reverse primers used to generate the gRNA by in vitro transcription (IVT) was done with the help of the Invitrogen™ GeneArt™ CRISPR Search and Design tool (http://www.thermofisher.com/crisprdesign). The CRISPR design tool identifies the top six CRISPR target sequences. Among them, three paired primers were selected for each target gene. (The different CRISPR sequences used are listed (Additional file 1: Table S1.) The full-length gRNA DNA template was generated by PCR using the previously forward and reverse overlapping oligonucleotides. These primers contain the target DNA sequence (CRISPR sequence) and the Tracr Fragment and T7 Primer Mix. They play the dual role of primer and template. The gRNA sequence has been generated by in vitro transcription using the TranscriptAid™ Enzyme Mix and then it has been purified using GeneJet™ RNA Purification Micro Columns. All the procedures were done using the GeneArt™ Precision gRNA Synthesis Kit (Invitrogen). Before the day of the transfection, MDA-MB-231 and MDA-MB-468 cells were trypsinized and seeded into a 24-well plate (3.2 × 105 cells/mL). After changing media, the transfection was done using Opti-MEM I reduced serum medium, GeneArt™ CRISPR Nuclease mRNA (Invitrogen, Thermo Fisher Scientific; Princeton, NJ), and Lipofectamine® CRISPRMax™ Reagent (Invitrogen). Three gRNAs were used to target each gene (Additional file 1: Table S1). Two controls (one media only control another Lipofectamine only control) were used. Cells were maintained for 48 h at 37 °C and 5% CO2.
Subcutaneous Xenograft experiments
T-cell-deficient, athymic nude (Crl: NIH-Foxn1rnu) female rats aged 21–28 days old were purchased from Charles River (Portage, MI, USA) and maintained in a pathogen-free condition with a 12-h light/dark cycle at our institute’s animal research facility. Animal care and treatment were performed in accordance with institutional guidelines, and protocols were approved by the Rutgers Institutional Animal Care and Facilities Committee and complied with National Institutes of Health policy.
MDA-MB 231 cells, at a final concentration of 1 × 107 cells/rat in 200-µl of PBS-50% Matrigel (BD Biosciences; San Jose, CA) mixture, were injected subcutaneously (SC) into the right flank of the athymic nude female rats (Crl: NIH-Foxn1rnu; Charles River; Portage, MI, USA). After tumors reached a diameter of approximately 50 mm3, the animals were randomly assigned to different treatment groups and injected s.c. daily with saline (control), NTX (10 mg/kg Sigma-Aldrich, St. Louis, MO, USA), PRO (10 mg/kg, Sigma-Aldrich, St. Louis, MO, USA), or in combination for 3 weeks. Tumor shrinkage/growth was measured in animals daily and animal weights were measured every other day. Animals were euthanized when the tumor reached 5000 mm3. Three dimensions of each tumor were measured using electronic calipers, and tumor volumes were calculated by the formula L x W2/2. The mean ± SEM of tumor volume was calculated weekly for each experimental group and presented. After the animals were sacrificed, the tumors were excised and were snap-frozen by immersion in liquid nitrogen and stored at − 80 °C until further use.
Cell growth assay
Experiments were performed in 96‐well plates over 72 h (1.5 × 105 cells/mL). Cell growth was estimated at 0, 24, 48, and 72 h by a colorimetric assay based on the conversion of tetrazolium dye (MTT) to a blue formazan product by live mitochondria. Optical density at 595 nm corresponding to solubilized formazan was read for each well on a Multiskan FC (Thermo Scientific Branchburg, NJ, USA). Data are presented as percent of control at 72 h.
Transwell (8 µm pore size insert) migration assay was performed in 24-well plates. Lower compartments contained DMEM (Gibco, Thermo Scientific, Branchburg, NJ, USA) with 0.5% FBS (Gibco, Thermo Scientific, Branchburg, NJ, USA) and 40 µg/mL collagen (Sigma-Aldrich, St. Louis, MO, USA). Upper compartments contained the cells in DMEM with 0.5% (5 × 105 cells/mL). Cells were incubated in the Transwell plate at 37 °C and 5% CO2 for 3 h. Then they were fixed on the lower side of the insert filter using 5% glutaraldehyde for 10 min and stained with crystal violet for 20 min. Cell density was measured using Image J analysis software (National Institutes of Health).
Colony formation assay
Experiments were performed in 6-well plates (103 cells/mL). Cells were incubated for 14 days at 37 °C and 5% CO2, with fresh media added every 3 days. After carefully washing with 1X phosphate-buffered saline (PBS), colonies were stained with crystal violet for 15 min. Cell density was measured using Image J analysis software (National Institutes of Health).
Cells were treated by SB-216763 (Sigma-Aldrich; Allentown, PA), wortmannin (Sigma-Aldrich, St. Louis, MO, USA), epinephrine ((Sigma-Aldrich, St. Louis, MO, USA)), and [D-Ala2, N-MePhe4, Gly-ol]-enkephalin (DAMGO) ((Sigma-Aldrich, St. Louis, MO, USA)). Experiments were performed in 96-well or 6-well plates (5 × 105 cells/mL).
Cancer stem cell gene qPCR array
The human cancer stem cell RT2 Profiler PCR array (Qiagen; Frederick, MD) profiles the expression of 84 genes linked to cancer stem cells (CSCs). The genes profiled with this array include CSC molecular markers and genes regulating CSC proliferation, self-renewal, or pluripotency and also those involved in CSC asymmetric cell division, migration, and metastasis. They are all listed in Additional file 1: Table S2. Synthesis of cDNA was done using the RT2 First Strand Kit (Qiagen; Frederick, MD) with 0.5 µg total RNA. Real-time PCR of the 96-well PCR array was performed using RT2 SYBR Green Mastermix on an ABI 7500 (Applied Biosystems™; Wilmington, DL). A set of controls present in the array enabled analysis using the ΔΔCT method using the Data Analysis Center (Qiagen; Frederick, MD). Fold regulation values were then used to design a heatmap with gplot package (heatmap2) in R Studio and to analyze canonical pathways with Ingenuity Pathway Analysis (IPA, Qiagen; Frederick, MD).
Real-time RT-PCR was carried out in 384-well plates in a final volume of 10 µL with 1:30 dilution cDNA (1 µg/µL) in Power SYBR® Green PCR Master Mix (Applied Biosystems™) using the comparative CT (ΔΔCT) standard run on a the QuantStudio™ 7 Flex Real-Time PCR System (Applied Biosystems™, Wilmington, DL). The primers were designed using Primer 3 (http://bioinfo.ut.ee/primer3/). Analysis of the results was done using QuantStudio™ Real-Time PCR (Applied Biosystems™, Wilmington, DL) software. The primers used are listed in Additional file 1: Table S3; four genes were used as housekeeping control genes: GAPDH, ATCB, B2M, and HPRT1.
For immunoblotting, about 50 µg of total protein from MDA-MB-231 or MDA-MB-468 cells were run in 12% SDS PAGE and transferred to PVDF membranes (GE Health Care, Piscataway, NY, USA) at 75 V for 2 h at 4 °C. The membranes were blocked in 5% non-fat dry milk-TBS-0.1% Tween 20 (TBST) at room temperature for 2 h and then incubated with primary antibody in the same buffer at 4 °C overnight. The following primary antibodies were used: antibodies against B2AR (1:1000, Abcam, Branford, CT, RRID: AB_2747383), MOR (1:1000, Cell Signaling, Danvers, MA, RRID: AB_177512) GSK3α (1:1000, Cell Signaling, Danvers, MA, RRID: RRID:AB_659836), GSK3β (1:1000, Cell Signaling, Danvers, MA, RRID: RRID:AB_10998934), and β-actin (1:10,000, Abcam, Branford, CT, RRID:AB_2737344). The membranes were washed for 10 min 6 times in TBST and then incubated with corresponding peroxidase conjugated secondary antibody (1:5000, Vector Labs, Burlingame, CA) at room temperature for 1 h. The membranes were washed in TBST for 10 min 6 times and incubated with ECL reagent (Thermo Fisher Scientific, and were developed on the film by autoradiography. The protein band intensities were determined by Image J analysis software (National Institutes of Health, Bethesda, MD) and normalized with corresponding β-actin band intensity.
Prism 7.0 (GraphPad Software) was used for statistical analysis. The data were expressed as means ± SEM values of three independent experiments. Significance of difference was analyzed using one-way analysis of variance (ANOVA) with post hoc analysis using the Newman–Keuls posttest. For all analyses, values of p < 0.05 were considered significant.
Inhibition of the expression of B2AR and MOR receptors alone or together reduces cell proliferation, cell migration, and colony formation of MDA-MB-231 and MDA-MB-468 cells
We have recently shown that B2AR and MOR interact to reduce the growth and proliferation of three human breast cancer cells . We first confirmed the interactive actions of B2AR and MOR on MDA-MB-231 and MDA-MB-468 cell functions. Two technological approaches were used to determine the role of B2AR and MOR in the regulation of growth, migration, and invasion properties of MDA-MB-231 and MDA-MB-468 cells: one involving CRISPR technology to knock down B2AR and MOR receptors and another involving the pharmacological blockade of these two receptors with the use of PRO and NTX.
In the CRISPR knockdown approach, 3 gRNAs were used for each gene and no treatment control and lipofectamine treatment was used as CRISPR negative control. Western blot assays were conducted to verify the effectiveness of the gene knockdown by CRISPR (Additional file 1: Fig. S1). Cell growth response to various treatments was evaluated at 24-h intervals for 72 h using MTT assay. We show that CRISPR knockdown of B2AR and MOR alone or in combination inhibited the growth of MDA-MB-231 (Fig. 1A) and MDA-MB-468 cells (Fig. 1G) at all time points as compared to those in control groups. The highest cell growth inhibition was seen in the B2AR and MOR combined knockdown group at all time points. The capacity of cells to migrate following B2AR, MOR, and B2AR + MOR knockdown was also evaluated using a Transwell membrane and collagen I as a chemoattractant molecule for 3 h (Fig. 1B, H, Additional file 1: Fig. S2A, B). The loss of expression of MOR or B2AR significantly reduced the capacity of MDA-MB-231 and MDA-MB-468 cells to migrate as compared to the controls. The lowest cell migration was seen when both receptors were inhibited. Finally, MOR and B2AR knockdown cells also show a lesser ability to form colonies compared to controls, with the double knockdown cells being the most impacted compared to the single knockdown (Fig. 1C, I, Additional file 1: Fig. S2C and S2D). These results suggest that both B2AR and MOR act individually or interactively to control MDA-MB-231 and MDA-MB-468 cell growth, migration, and invasion functions.
The effects of the treatment of 100 µM (IC50 dose; ) of PRO or NTX alone or in combination at 24-h intervals for 72 h on cell growth, colony formation, migration, and invasion were determined. Treatment with PRO and NTX alone or in combination significantly reduced cell growth response as determined by MTT assay (Fig. 1D, J). Interestingly, maximal growth inhibition was observed when NTX and PRO were combined, while only moderate inhibition was achieved by PRO and NTX when given alone as compared to the control group. The effects of these drugs on the capacity of cells to migrate were evaluated using a Transwell membrane and collagen I as a chemoattractant molecule (Fig. 1E, K, Additional file 1: Fig. S2A and S2B). PRO or NTX treatment alone decreased the migration of the MDA-MB-231 and MDA-MB-468 cells as compared to the controls. Combination treatment of these two drugs markedly decreased the capacity of cells to migrate as compared to the control. Next, the effectiveness of a 100 µM dose of PRO and NTX on colony formation of MDA-MB-231 and MDA-MB-468 was tested (Fig. 1F, L, Additional file 1: Fig. S2C and S2D). The results show that PRO and NTX moderately decreased the number of colonies alone. The combination of PRO + NTX produced the maximum inhibition of colony formation in these cells.
B2AR and MOR interaction may involve GSK3 signaling
The inhibition of B2AR or/and MOR receptors induced modifications of the MDA-MB-231 and MDA-MB-468 phenotype with a decrease in cell proliferation, cell migration, and colony formation; all three important functions are used by the cancer cell to maintain its growth and progression . Therefore, we determined if there was a common signaling mechanism these two receptors used to promote this cancer cell growth and progression. For this, we employed a cancer stem cell PCR array system which detects modulation of transcript abundance of sets of genes involved in regulation of cancer cell growth and progression (List of genes in Additional file 1: Table S2). We also employed both pharmacological and genetic approaches to knock down B2AR and MOR functions in MDA-MB-231 cells. Changes of the gene transcription profiles following knockdown of B2AR with or without MOR are shown in Fig. 2. Between 39 and 24 transcripts were found to be differentially expressed as shown in the heatmap for the gene knockdown (Fig. 2A) and the pharmacological blockade conditions (Fig. 2B). Upregulated genes are presented in red and downregulated genes are in green in the subset figures. Fold change values of the gene transcripts were then analyzed by Gene Ontology tools for different gene clusters and by IPA (Ingenuity Pathway Analysis) for canonical pathways. Significant changes were found in gene clusters involved in biological process such as cell movement, invasion, migration, adhesion, cell proliferation, apoptosis, and cell death of breast cancer cell lines (Fig. 2C) and the cell signaling pathways involved in these biological processes (Fig. 2D). The Venn diagrams, performed using the online Venny tool (https://bioinfogp.cnb.csic.es/tools/venny/), of the genes that were expressed in the three different conditions (B2AR inhibition, MOR inhibition, and both B2AR and MOR inhibition) identified several groups of genes modulated under a specific treatment condition, but MUC1 expression is only modulated in all three conditions in both gene knockdown (Fig. 2E) and pharmacological blockade studies (Fig. 2F).
Further analysis of PCR array gene data was performed using the IPA My Pathway tool, Molecular Activity Predictor (MAP), in order to identify the predictive pathways altered following genetic or pharmacological inhibition of B2AR with or without MOR. This analysis identified upregulation of the MUC1 gene . In addition, it identified downregulation of GSK3 signaling, which is known to negatively influence MUC1 and CTNNB1 expression  and positively affect STAT3, NFKB1, and NOTCH1 expression to control cell proliferation, apoptosis, colony formation, epithelial-mesenchymal transition, migration, and metastasis (Fig. 3A; ). Verification of the predicted gene expression changes was conducted using qPCR determination of GSK3 pathway related genes in MDA-MB-231 and MDA-MB-468 cells with knockdown of B2AR, MOR, or B2AR and MOR (achieved by both genetic and pharmacological approaches). As predicted, we found increased levels of MUC1 and CTNNB1 expression and decreased levels of GSK3, NOTCH1, STAT3, and NKB1 expression (Fig. 3B–E). These data led us to hypothesize that activation of B2AR and MOR receptors increases GSK3 signaling to control growth and progression of breast cancer cells.
To determine if the GSK3-MUC signaling pathway is also impacted by modulation of B2AR and/or MOR receptor activity in vivo, we studied the effects of PRO (10 mg/kg body weight; ), NTX (10 mg/kg body weight; ), and PRO + NTX on the growth of the tumor and on the expression of genes related to GSK3-MUC1 signaling in MDA-MB-231 cell rat xenografts. Like in vitro, in rat xenografts, the growth of tumors, as determined by the tumor volume, was reduced by both PRO and NTX. The tumor growth suppression effect of these two drugs was enhanced when they were combined (Fig. 4A–E). Also, like in vitro, we found increased levels of GSK3, MUC1, CTNNB1, and caspase 3 gene expression and decreased levels of GSK3, NOTCH1, STAT3, and NKB1 expression in MDA-MB-231 cell rat xenografts following B2AR inhibition by PRO and MOR inhibition by NTX. An additive effect of these drugs on GSK3 signaling molecules was also observed in MDA-MB-231 cell rat xenografts (Fig. 4F–M).
GSK3 may mediate B2AR and MOR regulated cell proliferation, cell migration, and colony formation in breast cancer cells
Because suppression of GSK3 signaling is identified as a putative mechanism involved in mediation of B2AR and MOR effects on breast cancer cell growth and progression, we first tested if GSK3 gene knockdown reduces MDA-MD-231 and MDA-MB-468 cell proliferation, migration, and colony formation in vitro. We employed both gene knockdown by CRISPR and pharmacological knockdown by a specific chemical inhibitor SD-216763. In the case of gene knockdown by CRISPR, Western blot data confirmed the effectiveness of the gene knockdown (Additional file 1: Fig. S3). GSK3 gene knockdown decreased the cell proliferation level for a period of 72 h as compared to control groups (Fig. 5A, G). The migration capacity of MDA-MD-231 and MDA-MB-468 cells was markedly decreased in GSK3 knockdown cells (Fig. 5B, H, Additional file 1: Fig. S2A and S2B). Additionally, the cells with GSK3 knockdown showed a lower number of colonies compared to the controls (Fig. 5C, I, Additional file 1: Fig. S2C and S2D). In the case of pharmaceutical knockdown by SD-216763, we also found this drug is effective in inhibition of cell proliferation (Fig. 5D, J) for a period of 72 h and determined that a dose of 29.2 µM produced the half-maximal effect or IC50. The IC50 dose of SD-216763 also decreased cell migration (Fig. 5E, K, Additional file 1: Fig. S2A and S2B) and reduced the number of colony formation (Fig. 5F, L, Additional file 1: Fig. S2C and S2D) in both MDA-MB-231 and MDA-MB-468 breast cancer cell lines.
To determine if GSK3 is involved in B2AR and MOR receptor signaling, we tested the effect of activation of one partner in reversing the effects of downregulation of the other partner on cell proliferation of MDA-MB-231 and MDA-MB-468 cells. Cell proliferation of untreated control, vehicle control (lipofectamine), B2AR knockdown, and MOR knockdown cells in the presence and absence of GSK3 activator wortmannin are compared and the data are presented in Fig. 6A, B. These data revealed that treatment with wortmannin dose-dependently increased % cell proliferation in control-treated cells and reversed the inhibitory effects of B2AR, MOR, and B2AR + MOR knockdown in MDA-MB-231 and MDA-MB-468 cells. We also tested the effects of a B2AR agonist epinephrine or a MOR agonist DAMGO on cell proliferation in GSK3 knockdown cells. Interestingly, B2AR agonist epinephrine at a dose of 100 nM increased the % cell proliferation in control cells but failed to alter the inhibitory effect of GSK3 knockdown in MDA-MB-231 cells (Fig. 6C) and MDA-MB-468 (Fig. 6D). Similarly, DAMGO at a dose of 100 nM increased the % cell proliferation in control cells and failed to reverse the inhibitory effect of GSK3 knockdown in MDA-MB-231 and MDA-MB-468 cells (Fig. 6E, F).
B2AR and MOR are G-protein-coupled receptors which are known to heterodimerize with closely related members, resulting in the modulation of their functions and signal transduction . Through the inhibition by knockdown of each of the receptors, we demonstrated that the loss of expression of either of these receptors reduces the growth and the capacity to migrate and to form colonies of MDA-MB-231 and MDA-MB-468 cells. The study also identifies synergistic effects following suppression of both receptors on these cell functions. Together, these data suggest the possibility of a role of a B2AR and MOR interaction in the regulation of growth and progression of MDA-MB-231 and MDA-MB-468 cells.
A role of B2AR and MOR interaction in the regulation of cancer cell functions is further highlighted by the demonstration of expression changes in subcellular molecules involved in the regulation of cell invasion, migration, apoptosis, and cell proliferation following suppression of B2AR and/or MOR activity both by genetic modifications using CRSPR and by pharmacological approaches. Our PCR array data also identified some critical signaling pathways. Of these pathways, some noticeable pathways are NANOG signaling, a key regulator of embryonic development and cellular reprogramming that possesses protumorigenic attributes and is broadly expressed in human cancers ; Notch signaling, commonly activated in cancer and plays a key role in the progress of cancer  and has been shown to be regulated by the norepinephrine-activated β2-AR–PKA–mTOR pathway ; Stat3 signaling, also constitutively activated in many cancers, which regulates cellular proliferation, invasion, migration, and angiogenesis that are critical for cancer metastasis [40, 41]; and also PTEN, Wnt/b-catenin, IL-6, p53, and PI3K/AKT signaling, known to be modified in many cancers. Some of the gene expression data are also in agreement with previous findings that MOR regulates cancer cell proliferation, migration, and EMT formation by using PI3K, Akt, and STAT3 signaling [23, 42], while B2AR actions on cancer cells involve P38MAPK signaling [11, 43]. Bioinformatic analysis of the gene expression data also identifies MUC1 as the signaling molecule affected by the inhibition of both B2AR and MOR. MUC1 has been shown to be overexpressed in a variety of epithelial cancers, including breast cancer [44, 45]. Several studies have indicated that MUC1 plays a critical role in the transcriptional regulation of genes associated with tumor invasion, metastasis, angiogenesis, proliferation, apoptosis, drug resistance, inflammation, and immune regulation [31, 41, 46]. In epithelial cells, disruption of the interaction between E-cadherin and β-catenin by overexpression of MUC1 may be important in the progression to carcinoma . MUC1 is shown to contribute to the malignant phenotype by preventing GSK3β-mediated phosphorylation of β-catenin and thereby stabilizes β-catenin . Hence, the interaction between MUC1 and β-catenin is tightly regulated by GSK3β, and MUC1 may transiently stabilize β-catenin before GSK3β induces dissociation of the MUC1-β-catenin complex . GSK3β decreases binding of β-catenin to MUC1 and stimulates the association of β-catenin and E-cadherin. We confirmed the involvement of GSK3 signaling by qRT-PCR determination of genes in B2AR and MOR knockdown cells by showing modification of the expression of different genes implicated in the pathways. Our data show a downregulation of the expression of GSK3 but an upregulation of MUC1and CTNNB1 (catenin- β], which is consistent with the view that GSK3 negatively influences the expression of MUC1 and β-catenin. In addition, we found downregulation of NOTCH1, STAT3, NF-κB, and upregulation of caspase 3. These data are consistent with several previously observed findings. For example, it has been found that a reduction in GSK-3 phosphorylation leads to a lower expression of Notch pathway members and reduces cancer cell growth . Also, activation of STAT3 has been shown to be dependent upon GSK in various cell systems . In addition, GSK3 has been shown to interact with NF-κB to control cellular apoptosis [49, 50]. The involvement of GSK3 signaling in the mediation of B2AR and MOR-regulated MDA-MB-231 and MDA-MB-468 cell growth and progression is further established from our data showing that the inhibition of GSK3 with chemical antagonist SD-216763 or using CRISPR technology decreased cell proliferation, cell migration, and colony formation in MDA-MB-231 and MDA-MB-468 cells. This shows that the inhibition of GSK3 in MDA-MB-231 and MDA-MB-468 cells induced the similar phenotypic changes as the inhibition of B2AR and/or MOR. In B2AR and/or MOR knockdown MDA-MB-231 and MDA-MB-468 cells, a treatment by wortmannin, agonist of GSK3, restored the proliferation rate observed in the control cells. In addition, our data indicated that activation of B2AR by norepinephrine or MOR by DAMGO has no effect on cell proliferation of GSK3 knockdown MDA-MB-231 and MDA-MB-468 cells. These data suggest that GSK3 signaling is involved in the mediation of B2AR and MOR regulated MDA-MB-231 and MDA-MB-468 cell growth and progression, possibly following dimerization of these two receptors.
Previously it has been shown that GSK3β signaling plays an important role in opioid-mediated apoptosis in breast cancer cells . MOR activation by morphine has also been shown to increase phosphorylation of GSK-3β and ERK in the heart . In the L6 skeletal muscle cells, it has been shown that adrenergic stimulation through β2 adrenoceptors, but not involving cyclic AMP or Gi, activates a PI3K pathway that stimulates GSK3 [53, 54]. Evidence from cultured cell studies also support interaction of B2AR with GSK3-Wnt signaling via G proteins . The data presented here add the role of GSK3 in mediation of B2AR and MOR regulated TNBC cell growth and progression.
Whether or not B2AR and MOR interact similarly to affects non-TNBC cell growth and migration is not addressed in this study. In vivo studies of B2AR signaling in animal models of breast tumors demonstrated an association with increased nodal involvement and the development of tumor metastasis . Also, B2AR antagonist PRO decreased proliferation, migration and invasion of TNBC and non-TNBC cells both in vitro and in vivo systems [11, 29, 57, 58]. In the case of MOR signaling, significant predominance of Ki-67 and μ-opioid receptor immuno-expression in the lymph nodes was observed in breast cancer patients with opioid medications . Also, MOR antagonist naloxone or naltrexone increased the cell death while decreased the cell metastasis in an animal model of TNBC [29, 60]. In contrast, it was found that a MOR agonist morphine significantly reduced the cell vitality, growth, and colony formation rate in MCF-7 cells, but naloxone along with morphine did not reverse these effects. This indicates that the inhibition of MCF-7 cell growth and proliferation by morphine could be its independent effect, not associated with opioid receptors . Further research is needed to establish if the interactive role of B2AR and MOR similarly operates to control cellular functions in both TNBC and non-TNBC tumor cells.
Overall, we demonstrated in this study a mutual activation of B2AR and MOR through heterodimerization. The joint activation of these G protein-coupled receptors induces the activation of the tyrosine-kinase effector GSK3 to induce cell proliferation, cell migration, and colony formation possibly via alteration of MUC1, NOTCH, and STAT3 expression resulting in a more aggressive tumor profile (Fig. 7). A better understanding of the molecular mechanism implicated in the aggressiveness of a tumor led by adrenergic and opioidergic activation present B2AR, MOR, and GSK3 as potential targets for the treatment of cancers such as TNBC that are resistant to current chemotherapies.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files].
Analysis of variance
Cancer stem cells
[D-Ala2, N-MePhe4, Gly-ol]-enkephalin
Dulbecco's Modified Eagle Medium
Glycogen synthase kinase 3
Ingenuity Pathway Analysis
Polymerase chain reaction
5% Non-fat dry milk-TBS-0.1% Tween 20
Triple-negative breast cancer
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.
Foulkes WD, Smith IE, Reis-Filho JS. Triple-negative breast cancer. N Engl J Med. 2010;363(20):1938–44.
Liedtke C, Mazouni C, Hess KR, André F, Tordai A, Mejia JA, et al. Response to neoadjuvant therapy and long-term survival in patients with triple-negative breast cancer. J Clin Oncol Off J Am Soc Clin Oncol. 2008;26(8):1275–81.
Anders CK, Carey LA. Biology, metastatic patterns, and treatment of patients with triple-negative breast cancer. Clin Breast Cancer. 2009;9(Suppl 2):S73-81.
von Minckwitz G, Untch M, Blohmer J-U, Costa SD, Eidtmann H, Fasching PA, et al. Definition and impact of pathologic complete response on prognosis after neoadjuvant chemotherapy in various intrinsic breast cancer subtypes. J Clin Oncol Off J Am Soc Clin Oncol. 2012;30(15):1796–804.
Chiriac V-F, Baban A, Dumireascu DL. Psychological stress and breast cancer incidence: a systematic review. Clujul Med. 2018;91(1):18–26.
Buccheri G. Depressive reactions to lung cancer are common and often followed by a poor outcome. Eur Respir J. 1998;11(1):173–8.
Watson M, Haviland JS, Greer S, Davidson J, Bliss JM. Influence of psychological response on survival in breast cancer: a population-based cohort study. Lancet Lond Engl. 1999;354(9187):1331–6.
Satin JR, Linden W, Phillips MJ. Depression as a predictor of disease progression and mortality in cancer patients. Cancer. 2009;115(22):5349–61.
Braadland PR, Grytli HH, Ramberg H, Katz B, Kellman R, Gauthier-Landry L, et al. Low β 2 -adrenergic receptor level may promote development of castration resistant prostate cancer and altered steroid metabolism. Oncotarget. 2015;7(2):1878–94.
Ouyang X, Zhu Z, Yang C, Wang L, Ding G, Jiang F. Epinephrine increases malignancy of breast cancer through p38 MAPK signaling pathway in depressive disorders. Int J Clin Exp Pathol. 2019;12(6):1932–46.
Zhao Y, Li W. Beta-adrenergic signaling on neuroendocrine differentiation, angiogenesis, and metastasis in prostate cancer progression. Asian J Androl. 2019;21(3):253–9.
Powe DG, Voss MJ, Habashy HO, Zänker KS, Green AR, Ellis IO, et al. Alpha- and beta-adrenergic receptor (AR) protein expression is associated with poor clinical outcome in breast cancer: an immunohistochemical study. Breast Cancer Res Treat. 2011;130(2):457–63.
Kurozumi S, Kaira K, Matsumoto H, Hirakata T, Yokobori T, Inoue K, et al. β2-Adrenergic receptor expression is associated with biomarkers of tumor immunity and predicts poor prognosis in estrogen receptor-negative breast cancer. Breast Cancer Res Treat. 2019;177(3):603–10.
Algazi M, Plu-Bureau G, Flahault A, Dondon M-G, Lê MG. Could treatments with beta-blockers be associated with a reduction in cancer risk? Rev Epidemiol Sante Publique. 2004;52(1):53–65.
Sarkar DK, Zhang C, Murugan S, Dokur M, Boyadjieva NI, Ortigüela M, et al. Transplantation of β-endorphin neurons into the hypothalamus promotes immune function and restricts the growth and metastasis of mammary carcinoma. Cancer Res. 2011;71(19):6282–91.
Zhang C, Murugan S, Boyadjieva N, Jabbar S, Shrivastava P, Sarkar DK. Beta-endorphin cell therapy for cancer prevention. Cancer Prev Res (Phila). 2015;8(1):56–67.
Sarkar DK, Murugan S, Zhang C, Boyadjieva N. Regulation of cancer progression by β-endorphin neuron. Cancer Res. 2012;72(4):836–40.
Valentino RJ, Volkow ND. Untangling the complexity of opioid receptor function. Neuropsychopharmacology. 2018;43(13):2514–20.
Gach K, Piestrzeniewicz M, Fichna J, Stefanska B, Szemraj J, Janecka A. Opioid-induced regulation of mu-opioid receptor gene expression in the MCF-7 breast cancer cell line. Biochem Cell Biol Biochim Biol Cell. 2008;86(3):217–26.
Gach K, Wyrębska A, Fichna J, Janecka A. The role of morphine in regulation of cancer cell growth. Naunyn Schmiedebergs Arch Pharmacol. 2011;384(3):221–30.
Lennon FE, Moss J, Singleton PA. The μ-opioid receptor in cancer progression: is there a direct effect? Anesthesiol J Am Soc Anesthesiol. 2012;116(4):940–5.
Lennon FE, Mirzapoiazova T, Mambetsariev B, Poroyko VA, Salgia R, Moss J, et al. The Mu opioid receptor promotes opioid and growth factor-induced proliferation, migration and Epithelial Mesenchymal Transition (EMT) in human lung cancer. PLoS ONE. 2014;9(3):e91577.
Nguyen J, Luk K, Vang D, Soto W, Vincent L, Robiner S, et al. Morphine stimulates cancer progression and mast cell activation and impairs survival in transgenic mice with breast cancer. Br J Anaesth. 2014;113(Suppl 1):i4-13.
Mathew B, Lennon FE, Siegler JH, Mirzapoiazova T, Mambetsariev N, Sammani S, et al. Novel role of the Mu opioid receptor in lung cancer progression: a laboratory study. Anesth Analg. 2011;112(3):558–67.
Baker JG, Hill SJ, Summers RJ. Evolution of β-blockers: from anti-anginal drugs to ligand-directed signalling. Trends Pharmacol Sci. 2011;32(4):227–34.
Gach K, Piestrzeniewicz M, Fichna J, Stefanska B, Szemraj J, Janecka A. Opioid-induced regulation of mu-opioid receptor gene expression in the MCF-7 breast cancer cell line. Biochem Cell Biol. 2008;86:217–26.
Amaro F, Silva D, Reguengo H, Oliveira JC, Quintas C, Vale N, et al. β-adrenoceptor activation in breast MCF-10A cells induces a pattern of catecholamine production similar to that of tumorigenic MCF-7 cells. Int J Mol Sci. 2020;21:7968.
Murugan S, Rousseau B, Sarkar DK. Beta 2 Adrenergic receptor antagonist propranolol and opioidergic receptor antagonist naltrexone produce synergistic effects on breast cancer growth prevention by acting on cancer cells and immune environment in a preclinical model of breast cancer. Cancers (Basel). 2021;13(19):4858.
Root-Bernstein R, Churchill B, Turke M, Subhramanyam UKT, Labahn J. Mutual enhancement of opioid and adrenergic receptors by combinations of opioids and adrenergic ligands Is reflected in molecular complementarity of ligands: drug development possibilities. Int J Mol Sci. 2019;20(17):4137.
Samoshkin A, Convertino M, Viet CT, Gauthier J, Li Fraine SP, Sharif-Naeini R, et al. Structural and functional interactions between six-transmembrane μ-opioid receptors and β2-adrenoreceptors modulate opioid signaling. Sci Rep. 2015;5:18198.
Sever R, Brugge JS. Signal transduction in cancer. Cold Spring Harb Perspect Med. 2015;5(4):a006098.
Kufe DW. Mucins in cancer: function, prognosis and therapy. Nat Rev Cancer. 2009;9:874–85.
Huang L, Chen D, Liu D, Yin L, Kharbanda S, Kufe D. MUC1 oncoprotein blocks glycogen synthase kinase 3beta-mediated phosphorylation and degradation of beta-catenin. Cancer Res. 2005;65(22):10413–22.
Maji S, Shriwas O, Samal SK, Priyadarshini M, Rath R, Panda S, Das Majumdar SK, Muduly DK, Dash R. STAT3- and GSK3β-mediated Mcl-1 regulation modulates TPF resistance in oral squamous cell carcinoma. Carcinogenesis. 2019;40(1):173–83.
Jordan BA, Trapaidze N, Gomes I, Nivarthi R, Devi LA. Oligomerization of opioid receptors with beta 2-adrenergic receptors: a role in trafficking and mitogen-activated protein kinase activation. Proc Natl Acad Sci U S A. 2001;98(1):343–8.
Jeter CR, Yang T, Wang J, Chao H-P, Tang DG. Concise Review: NANOG in cancer stem cells and tumor development: an update and outstanding questions. Stem Cells. 2015;33(8):2381–90.
Yuan X, Wu H, Xu H, Xiong H, Chu Q, Yu S, et al. Notch signaling: an emerging therapeutic target for cancer treatment. Cancer Lett. 2015;369(1):20–7.
Chen H, Liu D, Yang Z, Sun L, Deng Q, Yang S, et al. Adrenergic signaling promotes angiogenesis through endothelial cell–tumor cell crosstalk. Endocr Relat Cancer. 2014;21(5):783–95.
Kamran MZ, Patil P, Gude RP. Role of STAT3 in Cancer Metastasis and Translational Advances. Biomed Res Int. 2013;2013: 421821.
Yuen JWF, So IYW, Kam AYF, Wong YH. Regulation of STAT3 by mu-opioid receptors in human neuroblastoma SH-SY5Y cells. NeuroReport. 2004;15(9):1431–5.
Kharmate G, Rajput PS, Lin YC, Kumar U. Inhibition of tumor promoting signals by activation of SSTR2 and opioid receptors in human breast cancer cells. Cancer Cell Int. 2013;13(1):93.
Xie WY, He RH, Zhang J, He YJ, Wan Z, Zhou CF, Tang YJ, Li Z, Mcleod HL, Liu J. β-blockers inhibit the viability of breast cancer cells by regulating the ERK/COX-2 signaling pathway and the drug response is affected by ADRB2 single-nucleotide polymorphisms. Oncol Rep. 2019;41(1):341–50.
Nath S, Mukherjee P. MUC1: a multifaceted oncoprotein with a key role in cancer progression. Trends Mol Med. 2014;20(6):332–42.
Hattrup CL, Gendler SJ. MUC1 alters oncogenic events and transcription in human breast cancer cells. Breast Cancer Res. 2006;8:R37.
Rowse GJ. Delayed mammary tumor progression in Muc-1 null mice. J Biol Chem. 1995;270:30093–101.
Kunnimalaiyaan S, Gamblin TC, Kunnimalaiyaan M. Glycogen synthase kinase-3 inhibitor AR-A014418 suppresses pancreatic cancer cell growth via inhibition of GSK-3-mediated Notch1 expression. HPB (Oxford). 2015;17(9):770–6.
Beurel E, Jope RS. Differential regulation of STAT family members by glycogen synthase kinase-3. J Biol Chem. 2008;283(32):21934–44.
Kadry MO, Abdel-Megeed RM, El-Meliegy E, Abdel-Hamid AZ. Crosstalk between GSK-3, c-Fos, NFκB and TNF-α signaling pathways play an ambitious role in Chitosan Nanoparticles Cancer Therapy. Toxicol Rep. 2018;5:723–7.
Chen S, Sun KX, Liu BL, Zong ZH, Zhao Y. The role of glycogen synthase kinase-3β (GSK-3β) in endometrial carcinoma: a carcinogenesis, progression, prognosis, and target therapy marker. Oncotarget. 2016;7(19):27538–51.
Zhao M, Zhou G, Zhang Y, Chen T, Sun X, Stuart C, Hanley G, Li J, Zhang J, Yin D. Beta-arrestin2 inhibits opioid-induced breast cancer cell death through Akt and caspase-8 pathways. Neoplasma. 2009;56:108–13.
He SF, Jin SY, Yang W, Pan YL, Huang J, Zhang SJ, Zhang L, Zhang Y. Cardiac μ-opioid receptor contributes to opioid-induced cardioprotection in chronic heart failure. Br J Anaesth. 2018;121(1):26–37.
Zhou J, Lal H, Chen X, Shang X, Song J, Li Y, Kerkela R, Doble BW, MacAulay K, DeCaul M, Koch WJ, Farber J, Woodgett J, Gao E, Force T. GSK-3alpha directly regulates beta-adrenergic signaling and the response of the heart to hemodynamic stress in mice. J Clin Invest. 2010;120(7):2280–91.
Yamamoto DL, Hutchinson DS, Bengtsson T. Beta(2)-Adrenergic activation increases glycogen synthesis in L6 skeletal muscle cells through a signalling pathway independent of cyclic AMP. Diabetologia. 2007;50(1):158–67.
Jernigan KK, Cselenyi CS, Thorne CA, Hanson AJ, Tahinci E, Hajicek N, Oldham WM, Lee LA, Hamm HE, Hepler JR, Kozasa T, Linder ME, Lee E. Gbetagamma activates GSK3 to promote LRP6-mediated beta-catenin transcriptional activity. Sci Signal. 2010;3(121):ra37.
Sloan EK, Priceman SJ, Cox BF, Yu S, Pimentel MA, Tangkanangnukul V, et al. The sympathetic nervous system induces a metastatic switch in primary breast cancer. Cancer Res. 2010;70:7042–52.
Barron TI, Sharp L, Visvanathan K. Beta-adrenergic blocking drugs in breast cancer: a perspective review. Ther Adv Med Oncol. 2012;4:113–25.
Spini A, Roberto G, Gini R, Bartolini C, Bazzani L, Donnini S, et al. Evidence of β-blockers drug repurposing for the treatment of triple negative breast cancer: a systematic review. Neoplasma. 2019;66:963–70.
De Sousa AM, Dantas TS, Barros Silva PG, Martins CDS, Freire GE, Junior HLR, et al. Analysis of the immunoexpression of opioid receptors and their correlation with markers of angiogenesis, cell proliferation and apoptosis in breast cancer. Asian Pac J Cancer Prev. 2021;22:633–40.
Bimonte S, Barbieri A, Cascella M, Rea D, Palma G, Del Vecchio V, et al. The effects of naloxone on human breast cancer progression: in vitro and in vivo studies on MDA.MB231 cells. Onco Targets Ther. 2018;11:185–91.
Chen Y, Qin Y, Li L, Chen J, Zhang X, Xie Y. Morphine can inhibit the growth of breast cancer MCF-7 cells by arresting the cell cycle and inducing apoptosis. Biol Pharm Bull. 2017;40(10):1686–92.
We thank Stacey Pontoriero for editing the manuscript.
This research was funded by the National Institutes of Health, Grant Number R01 CA208632.
Ethics approval and consent to participate
All animal care and procedures were approved by the Rutgers Institutional Animal Care and Facilities Committee (PROTO 999900286) and complied with National Institutes of Health policy.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Rousseau, B., Murugan, S., Palagani, A. et al. Beta 2 adrenergic receptor and mu opioid receptor interact to potentiate the aggressiveness of human breast cancer cell by activating the glycogen synthase kinase 3 signaling. Breast Cancer Res 24, 33 (2022). https://doi.org/10.1186/s13058-022-01526-y
- Triple-negative breast cancer
- Growth and metastasis
- Beta-adrenergic receptor
- Mu-opioid receptor
- GSK3 signaling
- Gene knockdown
- Pharmacological agents