Maybe the apparent contradiction can be solved by actually making a multipole expansion for this model. Perhaps the monopole shift is still positive, but the higher order corrections together lead to the overall negative contribution and lowering of the energy.
Otherwise, maybe it has to do with the shape of the nucleus: in the model there are only 2 point charges at a fixed distance from each other instead of a true volume of nuclear charge.