We consider energy solutions of the inhomogeneous parabolic p-Laplacien system ∂ t u − div(|∇u| p−2 ∇u) = −divg). We show in the case p ≥ 2 that if the right hand side g is locally in L ∞ (BMO), then u is locally in L ∞ (C 1 ), where C 1 is the 1-Hölder-Zygmund space. This is the borderline case of the Calderón-Zygmund theorey. We provide local quantitative estimates. We also show that finer properties of g are conserved by ∇u, e.g. Hölder continuity. Moreover, we prove a new decay for gradients of p-caloric solutions for all 2n n+2 < p < ∞.