We consider the azimuthal $cos varphi$ and $cos 2varphi$ distributions and the Callan-Gross ratio $R={rm d}sigma_L/{rm d}sigma_T$ in heavy-quark pair electroproduction, $lNrightarrow l^{prime}Qbar{Q}X$, as probes of linearly polarized gluons in unpolarized nucleons. Our analysis shows that the azimuthal asymmetries and Callan-Gross ratio are predicted to be large and very sensitive to the contribution of the gluonic counterpart of the Boer-Mulders function, $h_{1}^{perp g}$, describing the linear polarization of gluons inside unpolarized nucleon. In particular, the maximum values of the azimuthal distributions vary from 0 to 1 depending on $h_{1}^{perp g}$. We conclude that future measurements of these quantities at the proposed EIC and LHeC colliders could clarify in details the proton spin decomposition puzzle.