3. Discussion
Chemotherapy plays a pivotal role in the management of breast cancer. The chemotherapy regimens are built on anthracyclines, cyclophosphamides, taxanes, antimetabolites (5-fluorouracil, gemcitabine, capecitabine), navelbine [
49], targeted therapeutic agents as trastuzumab, pertuzumab and trastuzumab-emtansine, lapatinib [
50], endocrine therapy, including selective estrogen receptor modulators (SERMs), aromatase inhibitors and GNRH-analogs [
50] and novel therapeutic agents as PARP inhibitors [
51,
52,
53] or CDK4/6 (cyclin-dependent kinases) inhibitors [
54]. In this study we assessed those inhibitors that can be utilized in cell-based model systems, as their action does not require the activation in the liver, interaction with the immune system (as for humanized antibodies) or systemic endocrine loops (e.g., SERMs).
In this study we investigated whether there is an interaction between cytostatic bacterial metabolites and the above-mentioned chemotherapy agents. Numerous bacterial metabolites were identified, the majority of which have cytostatic properties [
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14]. The production of these metabolites decline in breast cancer patients, nevertheless, administration of minute quantities of these metabolites reduce the mitotic rate and the metastatic potential of the primary tumor [
4,
5,
6,
7].
Multiple chemotherapy agents were shown to modulate the composition of the microbiome [
55,
56], bacterial metabolism of chemotherapy agents were also evidenced [
32,
34,
35,
36,
37,
38,
39,
40] and the efficiency of humanized antibodies [
57] were also linked to compositional changes to the microbiome. These observations raised the possibility that bacterial metabolites may interfere with the cytostatic or cytotoxic effects of chemotherapy agents. Despite our negative observations on the interactions between bile acids and chemotherapy agents in pancreatic adenocarcinoma cells [
58,
59], in the current study we identified CAD and IPA that do interfere with PAC and 5FU. Neither metabolites displayed toxicity towards non-transformed cells in previous studies and can be applied in low concentrations [
5,
6].
CAD increased the Hill coefficient of 5FU suggesting a more collaborative binding of the drug molecules [
60] that has unknown pharmacological relevance. In contrast, IPA increased the IC
50 value of 5FU suggesting a lower efficiency that has negative pharmacological and, likely, clinical consequences. While in the case of PAC the IC
50 value was halved in the presence of IPA making this combination a potential. These findings also suggest that IPA may have adverse effects when PAC+5FU combinations are applied.
IPA is a bacterial metabolite that is the synthesized from tryptophan through deamination by tryptophanase (TnaA) [
5,
45]. A significant portion of tryptophan (4-6%) undergoes bacterial catabolism [
61]. Multiple studies have shown that disturbances to indole/tryptophan metabolism correlates with survival in breast cancer ([
6,
7,
62,
63], reviewed in [
6]). Our observations extend these by adding that higher IPA levels may support PAC responsiveness.
4. Materials and Methods
Chemicals
Bacterial metabolites (Cadaverine-CAD, cat # C8561; Indoxylsulfate-IS, cat # 13875; Indolepropionic acid-IPA, cat # 220027) were purchased from Sigma-Aldrich (St. Louis, MI, USA). All metabolites were dissolved in dimethyl-sulfoxide (DMSO) at a stock concentration of 100 mM. CAD
was used at concentrations of 0.8 μM, IS at 4 μM and IPA at 1 μM, corresponding to normal human serum concentrations of these metabolites [
41,
42,
43,
44,
45,
46]
.
Chemotherapy drugs, Irinotecan (IRI, cat # I1406), 5-fluorouracil (5FU, cat # F6627), Methotrexate (MTX, cat # PHR1396), Rucaparib (RUCA, cat # PZ0036) and Gemcitabine (GEM, cat # G6423) were from Sigma-Aldrich. The drugs IRI, 5-FU, MTX and RUCA were dissolved in DMSO at a stock concentration of 100 mM; GEM was dissolved in water at a stock concentration of 100 mM. Liposomal Encapsuled Doxorubicin (DOX-NP, cat # 300112) was purchased from Avanti and a stock solution of 50 mM was prepared. Paclitaxel (PAC, cat # A0451335) was from Thermo Fisher Scientific (Waltham, MA, USA) and 50 mM stock solution was prepared in DMSO.
Chemotherapy compounds were used at different concentrations as indicated in the figures.
Cell Line
The 4T1 breast cancer cell line was obtained from the American Type Culture Collection. Cells were cultured in RPMI-1640 medium (Sigma-Aldrich, cat # R5886) containing 10% fetal bovine serum (FBS), 1% penicillin/streptomycin, 2 mM L-glutamine and 1% pyruvate at 37°C in a humidified incubator with 5% CO2. Cells were regularly checked for Mycoplasma contamination.
MTT Assay
4T1 cells were plated in 96-well plates (1.5 x 103 cell/well). On the next day cells were treated with chemotherapy agents alone or in combination with bacterial metabolites for 48 h. After treatments cell numbers was determined using an MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assay. Briefly, cells were treated with MTT solution (0.5 mg/ml) and incubated at 37°C for 90 minutes. Then, culture medium was discarded and the formazan crystals were dissolved in DMSO. The absorbance was measured on a plate reader (Thermo Labsystems Multiskan MS, Walthman, MA, USA) at 540 nm. In the calculations, the absorbance values for the vehicle-treated cells were considered 1, and all treatment were expressed relative to 1.
Statistical Analysis
Each analyses were performed using GraphPad Prism 8 software. Experiments were repeated at least three times and results are presented as mean ± SEM values. Normal distribution of the values was tested using D’Agostino & Pearson normality test. Where appropriate, values were log-normalized or normalized using Box-Cox normalization method [
64]. Nonlinear regression was performed using the GraphPad “[Inhibitor] vs. response - Variable slope (four parameters)” utility, from which IC
50 and Hill slope values were obtained unless otherwise stated. Two-way analysis of variance test followed by Tukey’s honestly significant post hoc test were used for multiple comparisons.
Author Contributions
Conceptualization, E.M. and P.B.; methodology, E.M., P.B., E.J. and S.S; software, S.S., P.N and T.I.B.; validation, E.M. and P.B.; formal analysis, E.M., P.B. and S.S.; investigation, S.S., T.I.B. P.N. and E.J..; resources, E.M. and P.B.; data curation, S.S.; writing—original draft preparation, S.S., E.M. and P.B; writing—review and editing, E.M and P.B.; visualization, E.M.; supervision, E.M. and P.B.; project administration, E.M and S.S.; funding acquisition, E.M and P.B. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Cytostatic bacterial metabolites do not interfere with the cytostatic effect of doxorubicin. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with doxorubicin alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. Statistical difference between the inhibitory curves was determined using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. ## and ### indicate p < 0.01 and p < 0.001, respectively, DOX-treated vs. vehicle-treated cells. Abbreviations: CAD – cadaverine, DOX– doxorubicin, IPA – indolepropionic acid, IS – indoxylsulfate.
Figure 1.
Cytostatic bacterial metabolites do not interfere with the cytostatic effect of doxorubicin. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with doxorubicin alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. Statistical difference between the inhibitory curves was determined using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. ## and ### indicate p < 0.01 and p < 0.001, respectively, DOX-treated vs. vehicle-treated cells. Abbreviations: CAD – cadaverine, DOX– doxorubicin, IPA – indolepropionic acid, IS – indoxylsulfate.
Figure 2.
Cytostatic bacterial metabolites do not interfere with the cytostatic effect of gemcitabine. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with gemcitabine alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (three parameters)” utility) was performed on datasets to obtain IC50 values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values using the Shapiro-Wilk test. Dataset normality was achieved by the Box-Cox normalization method. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 values non-paired, two-sided t-test was applied. ## and ### indicate p < 0.01 and p < 0.001, respectively, GEM-treated vs. vehicle-treated cells. Abbreviations: CAD – cadaverine, GEM – gemcitabine, IPA – indolepropionic acid, IS – indoxylsulfate.
Figure 2.
Cytostatic bacterial metabolites do not interfere with the cytostatic effect of gemcitabine. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with gemcitabine alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (three parameters)” utility) was performed on datasets to obtain IC50 values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values using the Shapiro-Wilk test. Dataset normality was achieved by the Box-Cox normalization method. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 values non-paired, two-sided t-test was applied. ## and ### indicate p < 0.01 and p < 0.001, respectively, GEM-treated vs. vehicle-treated cells. Abbreviations: CAD – cadaverine, GEM – gemcitabine, IPA – indolepropionic acid, IS – indoxylsulfate.
Figure 3.
Cytostatic bacterial metabolites do not interfere with the cytostatic effect of irinotecan. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with irinotecan alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. ## and ### indicate p < 0.01 and p < 0.001, respectively, IRI-treated vs. vehicle-treated cells. Abbreviations: CAD – cadaverine, IPA – indolepropionic acid, IRI – irinotecan, IS – indoxylsulfate.
Figure 3.
Cytostatic bacterial metabolites do not interfere with the cytostatic effect of irinotecan. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with irinotecan alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. ## and ### indicate p < 0.01 and p < 0.001, respectively, IRI-treated vs. vehicle-treated cells. Abbreviations: CAD – cadaverine, IPA – indolepropionic acid, IRI – irinotecan, IS – indoxylsulfate.
Figure 4.
Cytostatic bacterial metabolites do not interfere with the cytostatic effect of methotrexate. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with methotrexate alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. ### indicate p < 0.001, respectively, MTX-treated vs. vehicle-treated cells. Abbreviations: CAD – cadaverine, IPA – indolepropionic acid, IS – indoxylsulfate, MTX – methotrexate.
Figure 4.
Cytostatic bacterial metabolites do not interfere with the cytostatic effect of methotrexate. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with methotrexate alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. ### indicate p < 0.001, respectively, MTX-treated vs. vehicle-treated cells. Abbreviations: CAD – cadaverine, IPA – indolepropionic acid, IS – indoxylsulfate, MTX – methotrexate.
Figure 5.
Cytostatic bacterial metabolites do not interfere with the cytostatic effect of rucaparib. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with rucaparib alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. ## and ### indicate p < 0.01 and p < 0.001, respectively, RUCA-treated vs. vehicle-treated cells. Abbreviations: CAD – cadaverine, IPA – indolepropionic acid, IS – indoxylsulfate, RUCA – rucaparib.
Figure 5.
Cytostatic bacterial metabolites do not interfere with the cytostatic effect of rucaparib. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with rucaparib alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. ## and ### indicate p < 0.01 and p < 0.001, respectively, RUCA-treated vs. vehicle-treated cells. Abbreviations: CAD – cadaverine, IPA – indolepropionic acid, IS – indoxylsulfate, RUCA – rucaparib.
Figure 6.
Cadaverine and indolepropionic acid interfere with the cytostatic effects of 5-fluorouracil. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with 5-fluorouracil alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. To achieve normal distribution, datasets were log-normalized. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. #, ##, and ### indicate p < 0.05, p < 0.01, and p < 0.001, respectively, 5FU-treated vs. vehicle-treated cells. * represents significance at p < 0.05 between the indicated groups. Abbreviations: CAD – cadaverine, IPA – indolepropionic acid, IS – indoxylsulfate, 5FU – 5-fluorouracil.
Figure 6.
Cadaverine and indolepropionic acid interfere with the cytostatic effects of 5-fluorouracil. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with 5-fluorouracil alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. To achieve normal distribution, datasets were log-normalized. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. #, ##, and ### indicate p < 0.05, p < 0.01, and p < 0.001, respectively, 5FU-treated vs. vehicle-treated cells. * represents significance at p < 0.05 between the indicated groups. Abbreviations: CAD – cadaverine, IPA – indolepropionic acid, IS – indoxylsulfate, 5FU – 5-fluorouracil.
Figure 7.
Indolpropionic acid improves the cytostatic effect of paclitaxel. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with paclitaxel alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. ### indicate p < 0.001, respectively, 5FU-treated vs. non-treated cells. ** represents significance at p < 0.01 between the indicated groups. Abbreviations: CAD – cadaverine, IPA – indolepropionic acid, IS – indoxylsulfate, PAC - paclitaxel.
Figure 7.
Indolpropionic acid improves the cytostatic effect of paclitaxel. 4T1 cells were plated in 96-well plates (1500 cells/well). Cells were treated with paclitaxel alone or in combination with CAD (0.8 µM), IS (4 µM) or IPA (1 µM) for 48 hours, then cell numbers were determined by MTT assay. Data are presented as means ± SEM, from at least three biological replicates. Individual assays were measured in quadruplicate or in triplicate. Values were normalized to vehicle-treated cells (absorbance is equal to 1). Nonlinear regression (Graphpad “[Inhibitor] vs. response (four parameters)” utility) was performed on datasets to obtain IC50 and Hill slope values. Normality was determined for the inhibitory curves using the D’Agostino and Pearson normality test, while for the IC50 values or the Hill slope values using the Shapiro-Wilk test. Statistical difference between the inhibitory curves was performed using a two-way ANOVA test, and all data points were compared with each other (in Tukey post hoc tests). For the comparison of the IC50 and Hill slope values non-paired, two-sided t-test was applied. ### indicate p < 0.001, respectively, 5FU-treated vs. non-treated cells. ** represents significance at p < 0.01 between the indicated groups. Abbreviations: CAD – cadaverine, IPA – indolepropionic acid, IS – indoxylsulfate, PAC - paclitaxel.