This work aims at studying the $MBH-Mdyn$ relation of a sample of quasars by constraining their host galaxy masses through a full kinematical modelling of the cold gas kinematics thus avoiding all possible biases and effects introduced by the rough virial estimates usually adopted so far. For this purpose we retrieved public observations of $72$ quasar host galaxies observed in $[CII]158\mu m$ or $CO$ transitions with the Atacama Large Millimeter Array (ALMA). We then selected those quasars whose line emission is spatially resolved and performed a kinematic analysis on ALMA observations. We estimated the dynamical mass of the systems by modelling the gas kinematics with a rotating disc taking into account geometrical and instrumental effects. Our dynamical mass estimates, combined with $MBH$ obtained from literature and our own new $CIV\lambda 1550$ observations have allowed us to investigate the $ M_{BH}/M_{dyn}$ in the early Universe. Overall we obtained a sample of $10$ quasars at $z\sim 2-7$ in which line emission is detected with high SNR ($>5-10$) and the gas kinematics is spatially resolved and dominated by ordered rotation. The estimated dynamical masses place $6$ out of $10$ quasars above the local relation yielding to a $MBH/Mdyn$ ratios $\sim 10\times$ higher than those estimated in low-$z$ galaxies. On the other hand we found that $4$ quasars at $z\sim 4-6$ have dynamical-to-BH mass ratios consistent with what is observed in early type galaxies in the local Universe.