In this article, we used a hybrid simulation method to sample the conformational space to characterize the structural dynamics and global motions of WT SARS-CoV-2 Mpro and 48 mutants, including several mutations that appear in P.1, B.1.1.7, B.1.351, B.1.525 and B.1.429+B.1.427 variants. Integrated Hybrid methods combining NMA and MD have been useful to study the correlation between the complex structural dynamics of macromolecules and their functioning mechanisms. Here, we applied this hybrid approach to elucidate the effects of mutation in the structural dynamics of SARS-CoV-2 Mpro, considering their flexibility, solvent accessible surface area analyses, global movements, and catalytic dyad distance. Furthermore, some mutants showed significant changes in their structural dynamics and conformation, which could lead to distinct functional properties.