Nuclei do have a finite volume, and this correction is made in accurate computations. It must be pretty useless for chemistry, but on the measures of the energy of last ionisation for varied elements, you can observe it. Or as well on the energy of a pseudo-hydrogen atom made of one proton and one muon, as this "atom" is smaller hence more sensitive to the proton size (...and figures don't fit, that's one open research question presently, and nobody knows the consequences).
The usual model is to stray the nucleus' charge uniformly over its volume rather than over the nucleus' surface. Though, we might argue that since they repel an other, protons should concentrate a bit more at the nucleus' surface. I ignore if someone tried to model that.
The effect is that:
- Outside the nucleus, the field and the potential are the same as for a point charge.
- Inside the nucleus, the field is less steep because a fraction of the charge is behind you as you dive to the centre, so the potential is less deep when the charge is strayed.
- For the same electron orbital, the attraction energy would be less strong, which means that the adapted orbital is wider.
This correction is smaller than the relativistic effect of the electron's kinetic energy.
And by the way, I still want to know why the energy of the electrostatic interaction between the electron and the nucleus has no inertia. It should be about as big as the relativistic correction, hence much bigger than other corrections which are experimentally confirmed, but this inertia is never included in successful models.