Structural dynamics of non-synonymous SNPs of histone methyltransferase EZH2 involved in Weaver Syndrome

Enhancer of zeste homolog 2 (EZH2) is a histone H3 lysine 27 methyltransferases. Non-synonymous SNPs (nsSNPs) in the ezh2 gene may cause Weaver Syndrome that is a prominent and rare congenital disorder. Several EZH2 genetic variants have been characterized and reported, but there is no information available on their structural dynamics. Our study employs an in silico approach for structural and functional characterization of EZH2 nsSNPs. We identified 19 EZH2 nsSNPs majorly associated with Weaver Syndrome, among which four SET-domain nsSNPs including V621M, A677T, R679C, and H689Y significantly affected EZH2 structure. We conducted the triplicate of 100 ns of molecular dynamics simulations (MDS) to compare the dynamic interaction behaviors between wild-type and mutants of EZH2 with their substrate, H3K27me0 peptide. The simulation results revealed that the mutants had higher levels of structural variations as evidenced by secondary structure, density, distance plots, and principal component analysis. Compared to EZH2-WT, the mutants A677T and H689Y have shown lower binding energy with H3K27me0 peptide due to the denaturing of 310 helixes as exemplified by MM/PBSA calculations and secondary structure analysis. Our analysis shed light on mechanisms of structural variations of EZH2 nsSNPs and lay a stone to develop mutant-based therapeutic strategies for the design of target-specific scaffolds against Weaver Syndrome.