The propagation of heavy quarks (HQs) in the quark-gluon plasma (QGP) is described by means of a full Boltzmann transport approach. The nonperturbative dynamics and the interaction between HQs and the bulk is taken into account by means of a Quasi-Particle Model. Including the intense electromagnetic and vortical fields, we discuss their impact on the directed flow of neutral D mesons at RHIC energy, clarifying the role of this observable in giving information on QGP transport properties. We also show a novel study of the correlations between different flow coefficients of D mesons and soft hadrons at LHC energy within a coalescence plus fragmentation hadronization scheme and including event-by-event initial state fluctuations. This investigations can put further constraints on HQ transport coefficients towards a solid determination in comparison to the lattice QCD calculations.