<p>Supersymmetry beyond the TeV scale offers several theoretical and phenomenological advantages, such as accommodating the observed Higgs mass and alleviating the flavor and CP problems. However, flavor and CP observables still impose stringent constraints even at the PeV scale, motivating a systematic study of flavor symmetries in this regime. In this work, we investigate nucleon decay induced by dimension-five operators in supersymmetric standard models and examine how flavor symmetries, particularly of the Froggatt-Nielsen type, can suppress these operators. We perform a Bayesian analysis combining flavor, CP, and proton-decay observables to quantify the allowed parameter space and identify characteristic predictions. Our results demonstrate that a multi-messenger approach, integrating flavor, CP, and baryon-number-violating observables, is essential for probing the underlying structure of supersymmetry beyond the TeV scale.</p>