Breast cancer is a heterogeneous disease. In clinical practice, tumors are classified as hormonal receptor positive, Her2 positive and triple negative tumors. In previous works, our group defined a new hormonal receptor positive subgroup, the TN-like subtype, which has a prognosis and a molecular profile more similar to triple negative tumors. In this study, proteomics and Bayesian networks were used to characterize protein relationships in 106 breast tumor samples. Components obtained by these methods had a clear functional structure. The analysis of these components suggested differences in processes such as metastasis or proliferation between breast cancer subtypes, including our new subtype TN-like. In addition, one of the components, mainly related with metastasis, had prognostic value in this cohort. Functional approaches allow to build hypotheses about regulatory mechanisms and to establish new relationships among proteins in the breast cancer context.Author SummaryBreast cancer classification in the clinical practice is defined by three biomarkers (estrogen receptor, progesterone receptor and HER2) into hormone receptor positive, HER2+ and triple negative breast cancer (TNBC). Our group recently described a new ER+ subtype with molecular characteristics and prognosis similar to TNBC. In this study we propose a mathematical method, the Bayesian networks, as a useful tool to study protein interactions and differential biological processes in breast cancer subtypes, characterizing differences in relevant processes such as proliferation or metastasis and associated them with patient prognosis.