We perform a next to leading order QCD global analysis of nuclear deep inelastic scattering and Drell-Yan data using the convolution approach to parametrize nuclear parton densities. We find both a significant improvement in the agreement with data compared to previous extractions, and substantial differences in the scale dependence of nuclear effects compared to leading order analyses.