Underlying mechanisms of cyclic peptide inhibitors interrupting the interaction of CK2α/CK2β: comparative molecular dynamics simulation studies

Protein–protein interactions (PPIs) are fundamental to all biological processes. Recently, the CK2β-derived cyclic peptide Pc has been demonstrated to efficiently antagonize the CK2α/CK2β interaction and strongly affect the phosphorylation of CK2β-dependent CK2 substrate specificity. The binding affinity of Pc to CK2α is destroyed to different extents by two single-point mutations of Tyr188 to Ala (Y188A) and Phe190 to Ala (F190A), which exert negative effects on the inhibitory activity (IC50) of Pc against the CK2α/CK2β interaction from 3.0 μM to 54.0 μM and ≫100 μM, respectively. However, the structural influences of Y188A and F190A mutations on the CK2α–Pc complex remain unclear. In this study, comparative molecular dynamics (MD) simulations, principal component analysis (PCA), domain cross-correlation map (DCCM) analysis and energy calculations were performed on wild type (WT), Y188A mutant, and F190A mutant systems. The results revealed that ordered communications between hydrophobic and polar interactions were essential for CK2α–Pc binding in the WT system. In addition to the loss of the hydrogen bond between Gln36 of CK2α and Gly189 of Pc in the two mutants, the improper recognition mechanisms occurred through different pathways. These pathways included the weakened hydrophobic interactions in the Y188A mutant as well as decreased polar and hydrophobic interactions in the F190A mutant. The energy analysis results qualitatively elucidated the instability of the two mutants and energetic contributions of the key residues. This study not only revealed the structural mechanisms for the decreased binding affinity of Y188A and F190A mutant CK2α–Pc complexes, but also provided valuable clues for the rational design of CK2α/CK2β subunit interaction inhibitors with high affinity and specificity.