The zero modes method is applied in order to get action of the monodromy matrix entries onto off-shell Bethe vectors in quantum integrable models associated with $U_q(mathfrak{gl}_N)$-invariant $R$-matrices. The action formulas allow to get recurrence relations for off-shell Bethe vectors and for highest coefficients of the Bethe vectors scalar product.