A systematic molecular dynamics approach to the study of peptide Keap1-Nrf2 protein-protein interaction inhibitors and its application to p62 peptides.

Protein-protein interactions (PPIs) as drug targets have been gaining growing interest, though developing drug-like small molecule PPI inhibitors remains challenging. Peptide PPI inhibitors, which can provide informative data on the PPI interface, are good starting points to develop small molecule modulators. Computational methods combining molecular dynamics simulations and binding energy calculations could give both the structural and the energetic perspective of peptide PPI inhibitors. Herein, we set up a computational workflow to investigate Keap1-Nrf2 peptide PPI inhibitors and predict the activity of novel sequences. Furthermore, we applied this method to investigate p62 peptides as PPI inhibitors of Keap1-Nrf2 and explored the activity change induced by the phosphorylation of serine. Our results showed that because of the unfavorable solvation effects, the binding affinity of the phosphorylated p62 peptide is lower than the Nrf2 ETGE peptide. Our research results not only provide a useful method to investigate the Keap1-Nrf2 peptide inhibitors, but also give a good example to show how to incorporate computational methods into the study of peptide PPI inhibitors. Besides, applying this method to p62 peptides provides a detailed explanation for the expression of cytoprotective Nrf2 targets induced by p62 phosphorylation, which may benefit the further study of the crosstalk between the Keap1-Nrf2 pathway and p62-mediated selective autophagy.