Not Keplerian, because it is not a conic-section. It is not even explained by Newtonian gravity. In contrast, Kepler's laws are explained by newtonian gravity.
The lowest orbital-energy from Keplerian orbit is circular. And the orbits of stars are observed to be approximately circular. Hence:
$$
\frac{mv^2}{r} = \frac{GMm}{r^2} \quad\Longrightarrow\quad
v = \sqrt{\frac{GM}{r}}
$$
So, a circular keplerian orbit would imply speed dependent of the distance from the star to the center, proportional to $r^{-1/2}$. However, this is not observed. The observation seems to indicate that exists a certain "independence" of the distance from the center. Therefore, the orbits are not keplerian.
Since gas and dust is observed, then this must not be the problem. To fix this small problem, one possible solution is to postulate $M(r)\propto r$. Since this is not observed, then must be some non-observed kind of matter, some dark matter.
Another solution is to say the force is not $ F = \frac{GMm}{r^2}$, and then invent a force that works in this case as well: MOND, which you've pointed out.