We describe the self-adjoint realizations of the operator $H:=-i\alpha\cdot \nabla + m\beta + \mathbb V(x)$, for $m\in\mathbb R $, and $\mathbb V(x)= |x|^{-1} ( \nu \mathbb{I}_4 +\mu \beta -i \lambda \alpha\cdot{x}/{|x|}\,\beta)$, ...