We present a detailed study of the spatially resolved thermodynamic and hydrostatic mass profiles of the five most massive clusters detected at z~1 via the Sunyaev-Zel'dovich effect. These objects represent an ideal laboratory to test our models in a mass regime where structure formation is driven mainly by gravity. We present a method to study these objects that optimally exploits information from XMM-Newton and Chandra observations. The combination of Chandra’s excellent spatial resolution and XMM-Newton’s photon collecting power allows us to spatially resolve the profiles from the core to the outskirts, for the first time in such objects.  Evolution properties are investigated by comparison with the REXCESS  local galaxy cluster sample. Finally, we discuss the current limitations of this method in the context of  joint analysis of future Chandra and XMM large programs and, more generally, of  multi-wavelength efforts to study high redshift objects.