Instead of using the Coulomb potential, you would have to use the gravitational potential. I guess, that way, you could write down a Kohn-Sham equation for every planet.