A theoretical framework is developed for the evaluation of favorable residual stress profiles, suppressing fatigue damage initiation in rolling contact fatigue. Non-metallic inclusions at the microstructure of bearings are one of the most important reasons for fatigue damage initiation since they act as stress risers. In order to evaluate the stress state around such inclusions at the micro-scale, macroscopic stress histories are determined by Hertzian contact theory at different depths below the raceway for a typical roller bearing. These stress distributions are then used as far-field stresses for a micro-scale model accounting for single inclusions of different geometries and orientations. Eshelby's method is used to relate far-field and local stresses in the vicinity of inclusions. The von Mises stress criterion is then used as a conservative estimator of crack initiation due to micro-scale plasticity. The effect of compressive residual stresses added to the axial and circumferential normal stress components at different depths is analyzed. The von Mises stress field around different inclusions at different depths is investigated in order to determine the most critical case in terms of microscale plastic deformation. Finally, an optimization process is carried out in order to determine the residual stresses that minimize the maximum observed von Mises stress as a function of depth.