Cancer therapy is a strategic measure in inhibiting breast cancer stem cell (BCSC) pathways. Naringenin, a citrus flavonoid, was found to increase breast cancer cells’ sensitivity to chemotherapeutic agents. Bioinformatics study and 3D tumorsphere
in vitro
modeling in breast cancer (mammosphere) were used in this study, which aims to explore the potential therapeutic targets of naringenin (PTTNs) in inhibiting BCSCs. Bioinformatic analyses identified direct target proteins (DTPs), indirect target proteins (ITPs), naringenin-mediated proteins (NMPs), BCSC regulatory genes, and PTTNs. The PTTNs were further analyzed for gene ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, protein–protein interaction (PPI) networks, and hub protein selection. Mammospheres were cultured in serum-free media. The effects of naringenin were measured by MTT-based cytotoxicity, mammosphere forming potential (MFP), colony formation, scratch wound-healing assay, and flow cytometry-based cell cycle analyses and apoptosis assays. Gene expression analysis was performed using real-time quantitative polymerase chain reaction (q-RT PCR). Bioinformatics analysis revealed p53 and estrogen receptor alpha (ERα) as PTTNs, and KEGG pathway enrichment analysis revealed that TGF-ß and Wnt/ß-catenin pathways are regulated by PTTNs. Naringenin demonstrated cytotoxicity and inhibited mammosphere and colony formation, migration, and epithelial to mesenchymal transition in the mammosphere. The mRNA of tumor suppressors P53 and ERα were downregulated in the mammosphere, but were significantly upregulated upon naringenin treatment. By modulating the P53 and ERα mRNA, naringenin has the potential of inhibiting BCSCs. Further studies on the molecular mechanism and formulation of naringenin in BCSCs would be beneficial for its development as a BCSC-targeting drug.