New properties of multiple harmonic sums modulo $p$ and $p$-analogues of Leshchiner's series