Molecular dynamics simulations and principal component analysis on human laforin mutation W32G and W32G/K87A

Mutations in human laforin lead to an autosomal neurodegenerative disorder Lafora disease. In N-terminal carbohydrate binding domain of laforin, two mutations W32G and K87A are reported as highly disease causing laforin mutants. Experimental studies reported that mutations are responsible for the abolishment of glycogen binding which is a critical function of laforin. Our current computational study focused on the role of conformational changes in human laforin structure due to existing single mutation W32G and prepared double mutation W32G/K87A related to loss of glycogen binding. We performed 10 ns molecular dynamics (MD) simulation studies in the Gromacs package for both mutations and analyzed the trajectories. From the results, the global properties like root mean square deviation, root mean square fluctuation, radius of gyration, solvent accessible surface area and hydrogen bonds showed structural changes in atomic level observed in W32G and W32G/K87A laforin mutants. The conformational change induced by mutants influenced the loss of the overall stability of the native laforin. Moreover, the change in overall motion of protein was analyzed by principal component analysis and results showed protein clusters expanded more than native and also change in direction in case of double mutant in conformational space. Overall, our report provides theoretical information on loss of structure-function relationship due to flexible nature of laforin mutants. In conclusion, comparative MD simulation studies support the experimental data on W32G and W32G/K87A related to the lafora disease mechanism on glycogen binding.