Gradient estimate for solutions of a class of nonlinear elliptic equations below the duality exponent