In order to do that you would need to find a gravitational equivalent for every term. For Coulomb interactions this seems doable at first, by using Newton’s law of universal gravitation for example.
The density would need to be replaced with some other type of density, like a density of the planet’s mass.
In order to apply Kohn-Sham, the interacting planets would need to be replaced with non-interacting “pseudo”-planets that are moving in an effective potential.
Still, at first sight it would seem weird to try to model the solar system like this. Quantum effects aren’t relevant on this scale, so there’s no uncertainty about the speed and position of the planet. Also, DFT searches for the ground state, which in this case would entail all planets moving as close as possible to the sun.