Shikimate kinase (SK) enzyme is a suitable target for antimicrobial drugs as it is present in pathogenic microorganisms and absent in mammals. A complete understanding of the functioning of this enzyme can unveil novel methods to inactivate it. To do this, a clear understanding of SK performance is needed. Previously, the chemical step of SK was studied in detail, but a study of longer-term scale simulations is still missing. In the present work, we performed molecular dynamics (MD) simulations in the µs time scale that allowed us to explore further regions of the SK energy landscape than previously. Simulations were conducted on the wild-type (WT) enzyme and the R116A and R116K mutants. We analyzed the dynamics of the enzymes through standard MD tools, and we found that the global motions in the mutants were perturbed. These motions can be linked to the observed undetectable binding affinity of the WT enzyme and the R116A and R116K mutants.