Polymer melts undergoing large deformation by elongation are studied by molecular dynamics simulations of bead-spring chains in melts. By applying a primitive path analysis to strongly deformed polymer melts, the role of topological constraints in highly entangled polymer melts is investigated and quantified. We show that the overall, large scale conformations of the primitive paths (PPs) of stretched chains follow affine deformation while the number and the distribution of entanglement points along the PPs do not. Right after deformation, PPs of chains retract in both directions parallel and perpendicular to the elongation. Upon further relaxation we observe a long-lived clustering of entanglement points. Together with the delayed relaxation time this leads to a metastable inhomogeneous distribution of topological constraints in the melts.