Exploring the structural insights on human laforin mutation K87A in Lafora disease--a molecular dynamics study

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.