We will explain the integrable systems structure on the space of Willmore tori arising from the Novikov-Veselov hierarchy. This structure can be used to gain a picture of the moduli space of all (constrained) Willmore tori and to compute many examples of such surfaces.