Assume a system where there only two 1 kg solid iron balls floating in space. The two balls are touching each other, so the potential gravitational energy between them is 0. Now I move them 1000 m apart, now the potential gravitational energy between them is $\int_{2r}^{2r+1000} \frac{G}{x^2} dx$, which is approximately $1.07 \times 10^{-9}$ J (considering a radius derived from pure iron density). According to Einstein's $M = \dfrac{E}{c^2}$, the added mass to the system, which was originally 2 kg, is $1.19 \times 10^{-17}$ kg, which I am not sure it could be detected at all. But the fact is that we had a system, added energy to it, and its mass increased.
If instead of solid iron balls, we have sun-like stars (same mass and radius), and instead of 2, we have 90 billions, and the average distance between each two is $1.28 \times 10^{21}$ m (as in a spiral galaxy with 150 kly radius), and finally (and most importantly, as most of the potential gravitational energy is gained in close proximity) we assume the radius of two of the stars, if merged, is $\sqrt[3]{2}$ greater than the original radius (i.e. the volume is preserved), then we are separating 90 billions stars from each other, from a distance of 2 half-sphere barycenters (of the merged stars) to the average star distance in a spiral galaxy, (90 billions)² times.
Using that simple galaxy model, the potential gravity energy accounts for 99.9995 % of the total mass of the galaxy (my python code), spread diffusely through it.
So, logic says such mass exists, because the potential energy exists in the system, and judging from the numbers, it seems very relevant, but I have never heard anyone talking about it. Is such energy accounted for in the estimates of the total mass of visible matter in the galaxy? If not, was it ever considered as a possible explanation for the dark matter?