Supersymmetric and superintegrable models of quantum mechanics with matrix potentials will be discussed.
In particular, a countable set of exactly solvable quantum mechanical systems will be presented which admit a dynamical symmetry with respect to algebra o(4). This algebra is generated by the Runge-Lenz vector generalized to the case
of arbitrary spin. The presented systems describe neutral particles with non-trivial multipole momenta. Their spectra can be found algebraically like in the case of Hydrogen atom. Solutions for the systems with spins 1/2 and 1 are presented explicitly, solutions for spin 3/2 are given in implicit form.