The interaction between barrier-to-autointegration factor dimer (BAF2) and LEM domain of emerin (EmLEM) was studied by molecular simulation methods. Nonspecific fragment of double-strand DNA molecule was docked with each chain of BAF2 by ZDOCK program. The model of DNA2:BAF2:EmLEM was thus constructed. The mutant Gly25Glu of BAF2 was manually constructed to explore the detailed effect of the mutation on the binding of BAF2 and EmLEM. It has been experimentally suggested that point mutation Gly25Glu can disturb the binding between BAF2 and EmLEM. Then, molecular dynamics (MD) simulations were performed on DNA2:BAF2(WT):EmLEM and DNA2:BAF2(MT):EmLEM complexes. 30ns trajectories revealed that the trajectory fluctuations of MT complex are more violent than that of the WT complex. Further, the binding free energy analysis showed that the electronegative residues Asp57, Glu61 and Asp65 from chain A, glu36 from chain B of BAF2 mainly contribute to interact with EmLEM. Besides, a stable π–π stack between trp62 and phe39 from BAF2(WT) chain B is destroyed by Glu25 in BAF2(MT). As a result, trp62 forms an interaction with glu25, and phe39 converts to strengthen affinity to EmLEM. On the other hand, Trp62 from chain A also forms a strong interaction with MT Glu25. Thus, with the docking of DNA, BAF2(MT) has higher affinity with EmLEM than BAF2(WT). 
Molecular simulation investigation on the interaction between barrier-to-autointegration factor dimer or its Gly25Glu mutant and LEM domain of emerin