In this paper, we establish a Hölder-type quantitative estimate of unique continuation for solutions to the heat equation with Coulomb potentials in either a bounded convex domain or a $C^2$-smooth bounded domain. The approach is based on the frequency function method, as well as some parabolic-type Hardy inequalities.