Previous studies demonstrated the efficiency of the Molecular Dynamics with excited Normal modes (MDeNM) method on the characterization of large structural changes at a low computational cost. We present here MDeNM-EMfit, an extension of the original method designed to the flexible fit of structures into cryo-EM maps. Here, instead of a uniform exploration of the collective motions described by normal modes, sampling is directed towards conformations with increased correlations with the experimental map. Future perspectives to improve the accuracy of fitting and speed of calculations are discussed in light of the results.