In high-energy heavy-ion collisions, the energy density profile of the produced quark-gluon plasma and its space-time dynamics are sensitive to the shape and radial profiles of the nuclei, described by the collective nuclear structure parameters including quadrupole deformation β2, octupole deformation β3, radius R0 and surface diffuseness a. Using hydrodynamic model simulations, we find a general scaling relation between these parameters and a large class of experimental observables such as elliptic flow v2, triangular flow v3 and particle multiplicity distribution p(N ch ). In particular, we show that the ratio of these observables between two isobar collision systems depends only on the differences of these parameters. Using this scaling relation, we show how the nuclear structure parameters of 96 44 Ru and 96 40 Zr conspire to produce the non-monotonic centrality dependence of ratios of v2, v3 and p(N ch ) between 96 44 Ru+ 96 44 Ru and 96 40 Zr+ 96 40 Zr collisions, in agreement with measurements by the STAR Collaboration. This scaling approach towards heavy-ion observables demonstrates that isobar collisions is a precision tool to probe the shape and radial structures, including the neutron skin, of the atomic nuclei across energy scales.