Botulinum neurotoxins (BoNTs) are known as the most poisonous biological substances, and they are also used to treat a wide range of medical conditions as well as in the cosmetic applications. Recently, the complex structures of the BoNT/A receptor-binding domain (BoNT/A-RBD) and the synaptic vesicle protein 2C luminal domain (SV2C-LD) were determined by X-ray crystallography. In this article, the wild type (WT) and four mutants of the new structure are studied by molecular dynamics (MD) simulations. The differently decreased structural stabilities of the mutants relative to WT are shown to be consistent with the experimental data of binding affinities. The conformational changes of the five systems are explored by using principal component analysis (PCA) and free energy landscape (FEL) methods. Based on the calculation of interactions at the binding interface, we divide the interface between BoNT/A-RBD and SV2C-LD into two crucial binding regions. Through the comparison of WT and four mutants, we further propose the relationship between the conformational changes of BoNT/A-RBD:SV2C-LD and the interfacial interactions. This study would provide some new insights into the understanding of the dynamics and the interaction mechanism of BoNT/A-RBD:SV2C-LD.