Modulation of localized solutions in an inhomogeneous saturable nonlinear Schrödinger equation