The first results on the complete next-to-next-to-leading order (NNLO) Quantum Chromodynamic (QCD) corrections to the production of di-leptons at hadron colliders in large extra dimension models with spin-2 particles are reported in this article. In particular, we have computed these corrections to the invariant mass distribution of the di-leptons taking into account all the partonic subprocesses that contribute at NNLO. In these models, spin-2 particles couple through the energy-momentum tensor of the Standard Model (SM) with the universal coupling strength. The tensorial nature of the interaction and the presence of both quark annihilation and gluon fusion channels at the Born level make it challenging computationally and interesting phenomenologically. We have demonstrated numerically the importance of our results at the Large Hadron Collider (LHC) energies. The two loop corrections contribute an additional 10\% to the total cross section. We find that the QCD corrections are not only large but also important to make the predictions stable under renormalisation and factorisation scale variations providing an opportunity to stringently constrain the parameters of the models with a spin-2 particle.