Molecular simulation investigation on the interaction between barrier-to-autointegration factor or its Gly25Glu mutant and DNA

Mono-phosphorylation at Ser4 of barrier-to-autointegration factor (Banf1) significantly reduces its DNA binding capability by inducing critical changes in its local conformation and DNA binding surface

Barrier-to-autointegration factor (Banf1) is a small DNA-bridging protein. The binding status of Banf1 to DNA is regulated by its N-terminal phosphorylation and dephosphorylation, which plays a critical role in cell proliferation. Banf1 can be phosphorylated at Ser4 into mono-phosphorylated Banf1, which is further phosphorylated at Thr3 to form di-phosphorylated Banf1. It was observed decades ago that mono-phosphorylated Banf1 cannot bind to DNA. However, the underlying molecular- and atomic-level mechanisms remain unclear. A clear understanding of these mechanisms will aid in interfering with the cell proliferation process for better global health. Herein, we explored the detailed atomic bases of unphosphorylated Banf1-DNA binding and how mono- and di-phosphorylation of Banf1 impair these atomic bases to eliminate its DNA-binding capability, followed by exploring the DNA-binding capability of mono- and di-phosphorylation Banf1, using comprehensive and systematic molecular modelling and molecular dynamics simulations. This work presented in detail the residue-level binding energies, hydrogen bonds and water bridges between Banf1 and DNA, some of which have not been reported. Moreover, we revealed that mono-phosphorylation of Banf1 causes its N-terminal secondary structure changes, which in turn induce significant changes in Banf1's DNA binding surface, thus eliminating its DNA-binding capability. At the atomic level, we also uncovered the alterations in interactions due to the induction of mono-phosphorylation that result in the N-terminal secondary structure changes of Banf1. Additionally, our modelling showed that phosphorylated Banf1 with their dominant N-terminal secondary structures bind to DNA with a significantly lower affinity and the docked binding pose are not stable in MD simulations. These findings help future studies in predicting effect of mutations in Banf1 on its DNA-binding capability and open a novel avenue for the development of therapeutics such as cancer drugs, targeting cell proliferation by inducing conformational changes in Banf1's N-terminal domain.

Molecular simulation investigation on the interaction between barrier-to-autointegration factor dimer or its Gly25Glu mutant and LEM domain of emerin

The interaction between barrier-to-autointegration factor dimer (BAF2) and LEM domain of emerin (Em^LEM) 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:Em^LEM 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 Em^LEM. It has been experimentally suggested that point mutation Gly25Glu can disturb the binding between BAF2 and Em^LEM. Then, molecular dynamics (MD) simulations were performed on DNA2:BAF2(WT):Em^LEM and DNA2:BAF2(MT):Em^LEM 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 Em^LEM. 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 Em^LEM. 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 Em^LEM than BAF2(WT).