Ball bearings are one of the most critical components of rotating machines. They ensure shaft support and friction reduction, thus their malfunctioning directly affects the machine’s performance. As a consequence, it is necessary to monitor the health conditions of such a component to avoid major degradations which could permanently damage the entire machine. In this context, HMS (Health Monitoring Systems) and PHM (Prognosis and Health Monitoring) methodologies propose a wide range of algorithms for bearing diagnosis and prognosis. The present article proposes an end-to-end PHM approach for ball bearing RUL (Remaining Useful Life) estimation. The proposed methodology is composed of three main steps: HI (Health Indicator) construction, bearing diagnosis and RUL estimation. The HI is obtained by processing non-stationary vibration data with the MODWPT (Maximum Overlap Discrete Wavelet Packet Transform). After that, a degradation profile is defined and coupled with crack initiation and crack propagation fatigue models. Lastly, a MB-HMM (Hidden Markov Model) is trained to capture the bearing degradation dynamics. This latter model is used to estimate the current degradation state as well as the RUL. The obtained results show good RUL prediction capabilities. In particular, the fatigue models allowed a reduction of the ML (Machine Learning) model size, improving the algorithms training phase.