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