Coagulation driven by supersonic turbulence is primarily an astrophysical problem because coagulation processes on Earth are normally associated with incompressible fluid flows at low Mach numbers, while dust aggregation in the interstellar medium for instance is an example of the opposite regime. We study coagulation of inertial particles in compressible turbulence using high-resolution direct and shock-capturing numerical simulations with a wide range of Mach numbers from nearly incompressible to moderately supersonic. The particle dynamics is simulated by representative particles and the effects on the size distribution and coagulation rate due to increasing Mach number is explored. We show that the time evolution of particle size distribution mainly depends on the compressibility (Mach number). We find that the average coagulation kernel C-ij scales linearly with the average Mach number M-rms multiplied by the combined size of the colliding particles, that is, < Cij > similar to <(a(i) + a(j))(3)> M-rms tau(-1)(eta), is proposed and can serve as a benchmark for future studies. We argue that the coagulation rate < R-c > is also enhanced by compressibility-induced compaction of particles.