Death-associated protein kinase 1 (DAPK1), as a calcium/calmodulin (CaM) regulated serine/threonine kinase, functions in apoptotic and autophagy pathways and represents an interesting drug target for inflammatory bowel disease and Alzheimer's disease. The crystal structure of the DAPK1 catalytic domain and the autoregulatory domain (ARD) in complex with CaM provides an understanding of CaM-dependent regulation of DAPK1 activity. However, the molecular basis of how distinct Trp305 (W305Y and W305D) mutations in the ARD modulate different DAPK1 activities remains unknown. Here, we performed multiple, μs-length molecular dynamics (MD) simulations of the DAPK1-CaM complex in three different (wild-type, W305Y, and W305D) states. MD simulations showed that the overall structural complex did not change significantly in the wild-type and W305Y systems, but underwent obvious conformational alteration in the W305D system. Dynamical cross-correlation and principal component analyses revealed that the W305D mutation enhanced the anti-correlated motions between the DAPK1 and CaM and sampled a broader distribution of conformational space relative to the wild-type and W305Y systems. Structural and energetical analyses further exhibited that CaM binding was unfavored in response to the W305D mutation, resulting in the decreased binding of CaM to the W305D mutant. Furthermore, the hydrogen bonds and salt bridges responsible for the loss of CaM binding on the interface of the DAPK1-CaM complex were identified in the W305D mutant. This result may provide insights into the key role of Trp305 in the regulation of CaM-mediated DAPK1 activity.
Keywords: binging free energy calculation; calmodulin; death-associated protein kinase 1 (DAPK1); molecular dynamics simulations; principal component analysis.