We study the propagation of heavy quarks (HQs) in the quark-gluon plasma (QGP) by means of a relativistic Boltzmann transport (RBT) approach. The non-perturbative interaction between heavy quarks and light quarks is described by means of a quasi-particle approach able to describe simultaneously the experimental data for the nuclear suppression factor $R_{rm AA}$ and the elliptic flow $v_2(p_T)$ of D mesons from RHIC to LHC energies. In the same framework we predict the B meson nuclear modification factor at LHC. Finally, we discuss the relevance of initial state fluctuations that allows to extend the analysis to high order anisotropic flows $v_n(p_T)$ as well as to investigate the role of QCD interaction in developing correlations between the light and the heavy flavour anisotropic flows.