The important force is the electromagnetic, gravity is very weak and you can ignore it in these systems.
For only one electron or one antiproton, the orbitals have the same shape and classification. The only difference is that the size depends on the mass, so the orbital to put the antiproton are much smaller than the orbitals to put the electrons. Note that something similar happens with muons that have an intermediate mass and the orbital to put them have an intermediate size. https://en.wikipedia.org/wiki/Bohr_radius
It's more difficult when you have many electrons or muons or antiprotons. (I guess nobody had measured a system with many muons or antiprotons.) If you have many electrons, the problem is that you must calculate the attraction of the nuclei and the repulsion of the other electrons, so the orbitals change. In particular, the filling rules https://en.wikipedia.org/wiki/Electron_configuration#Atoms:_... don't follow the energies of the orbitals of an isolated electron.
In a system with many muons or antiprotons, all of them will be closer and the repulsion will be bigger, and I expect weird filling rules.
Electrons: 1s_up, 1s_down, 2s_up, 2s_down, ...
¿¿¿Antiprotons: 1s_up, 2s_up, 1s_down, 2s_down, ... ???
It is possible to calculate the filling order for muons and antiprotons numerically, but I'm too lazy do do the calculation now and also the standard programs [1] [2] have a lot of hidden assumptions to make the calculations with electrons more efficient and I'm not sure how difficult is to tweak the entry files to calculate this fast enough. [3]
[1] https://en.wikipedia.org/wiki/Gaussian_(software)
[2] https://en.wikipedia.org/wiki/PSI_(computational_chemistry)
[3] Perhaps it's implemented and I just need to RTFM, but otherwise it looks too straightforward for a PhD thesis, but it may be a nice undergraduate thesis.