We introduce a novel real-time formulation of lattice NRQCD designed for simulations in the background of an highly occupied gluon field. By evolving quarks in the background of a dynamically evolving gluon field we computed the time-evolution of heavy-quarkonium spectral functions as well as the static and for finitely heavy quarks generalised potential. We conclude that the back reaction of the quarks is necessary for any binding process. Here we discuss the methodology, our results and the origin of the absence of a binding process.