Virtual element methods for parabolic problems on polygonal meshes