Threshold factorization of the Drell-Yan process at next-to-leading power