Transverse momentum dependent (TMD) parton distributions obtained from the parton branching (PB) method are combined with next-to-leading-order (NLO) calculations of Drell-Yan (DY) production. We apply the MC@NLO method for the hard process calculation and matching with the PB TMDs. We compute predictions for the transverse momentum, rapidity, and ϕ* spectra of Z-bosons. We find that the theoretical uncertainties of the predictions are dominated by the renormalization and factorization scale dependence, while the impact of TMD uncertainties is moderate. The theoretical predictions agree well, within uncertainties, with measurements at the Large Hadron Collider (LHC). In particular, we study the region of lowest transverse momenta at the LHC, and comment on its sensitivity to nonperturbative TMD contributions.