In this letter we present the three-loop virtual corrections to the Higgs boson production in the gluon fusion channel where finite top quark mass effects are taken into account. We perform an asymptotic expansion and manage to evaluate five terms in the expansion parameter $M_H^2/M_t^2$. A good convergence is observed almost until $M_H\approx 2M_t$.