The resolution of Dyson-Schwinger equations leads to the freezing of the QCD running coupling (effective charge) in the infrared, which is best understood as a dynamical generation of a gluon mass function, giving rise to a momentum dependence which is free from infrared divergences. We calculate the interquark static potential for heavy mesons by assuming that it is given by a massive One Gluon Exchange interaction and compare with phenomenologyical fits inspired by lattice QCD. We apply these potential forms to the description of quarkonia and conclude that, even though some aspects of the confinement mechanism are absent in the Dyson-Schwinger formalism, the spectrum can be reasonably reproduced. We discuss possible explanations for this outcome.