Lafora disease (LD) is an autosomal recessive, progressive form of myoclonus epilepsy which affects worldwide. LD occurs mainly in countries like southern Europe, northern Africa, South India, and in the Middle East. LD occurs with its onset mainly in teenagers and leads to decline and death within 2 to 10 years. The genes EPM2A and EPM2B are commonly involved in 90 % of LD cases. EPM2A codes for protein laforin which contains an amino terminal carbohydrate binding module (CBM) belonging to the CBM20 family and a carboxy terminal dual specificity phosphatase domain. Mutations in laforin are found to abolish glycogen binding and have been reported in wet lab methods. In order to investigate on structural insights on laforin mutation K81A, we performed molecular dynamics (MD) simulation studies for native and mutant protein. MD simulation results showed loss of stability due to mutation K87A which confirmed the structural reason for conformational changes observed in laforin. The conformational change of mutant laforin was confirmed by analysis using root mean square deviation, root mean square fluctuation, solvent accessibility surface area, radius of gyration, hydrogen bond, and principle component analysis. Our results identified that the flexibility of K87A mutated laforin structure, with replacement of acidic amino acid to aliphatic amino acid in functional CBM domain, have more impact in abolishing glycogen binding that favors LD.
Matched MeSH terms: Protein Tyrosine Phosphatases, Non-Receptor/chemistry*
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.
Matched MeSH terms: Protein Tyrosine Phosphatases, Non-Receptor/chemistry*