We analyze the contributions to the neutrinoless double \(\beta\) decay (\(0\nu\beta\beta\)-decay) coming from the Grand Unified Theory (GUT) constrained Minimal Supersymmetric Standard Model (MSSM) with trilinear R-parity breaking. We discuss the importance of two-nucleon and pion-exchange realizations of the quark-level \(0\nu\beta\beta\)-decay transitions. In this context, the questions of reliability of the calculated relevant nuclear matrix elements within the Renormalized Quasiparticle Random Phase Approximation (pn-RQRPA) for several medium and heavy open-shell nuclei are addressed. The importance of gluino and neutralino contributions to \(0\nu\beta\beta\)-decay is also analyzed. We review the present experiments and deduce limits on the trilinear R-parity breaking parameter \(\lambda_{111}'\) from the non-observability of \(0\nu\beta\beta\)-decay for different GUT constrained SUSY scenarios. In addition, a detailed study of limits on the MSSM parameter space coming from the \(B \to X_s \gamma\) processes by using the recent CLEO and OPAL results is performed. Some studies in respect to the future \(0\nu\beta\beta\)-decay project GENIUS are also presented.