We explore proton decay in a class of realistic supersymmetric flipped $SU(5)$ models supplemented by $U(1)_R$ symmetry which plays an essential role implementing hybrid inflation. Two distinct neutrino mass models, based on inverse seesaw and type I seesaw, are identified, with the latter arising from breaking nonrenormalizable superpotential terms. Depending model appropriate set intermediate...