Molecular selectivity design of mitogen-inducible gene-derived phosphopeptides between oncogenic HER kinases

Mitogen-inducible gene (MIG) is a natural negative regulator of the oncogenic HER kinase signaling by binding at the activation interface of kinase domain to disrupt the kinase dimerization. In this study, we systematically examine the binding structures, dynamics and energetics of MIG region 2 to four HER kinases based on their crystal or modeled complex structures, and identify an 8-mer phosphopeptide segment pYpY from the core strand sequence of MIG region 2 as the binding hotspot of MIG protein to HER kinases. We demonstrate that the small pYpY phosphopeptide can partially restore the binding affinity of full-length MIG protein, but exhibit a moderate selectivity over different HER kinases (S = 2.3-fold). In addition, the two phosphotyrosine residues pTyr394 and pTyr395 play an essential role in MIG-HER binding; dephosphorylation of them would fully eliminate the binding capability. A machine evolution algorithm is used to optimize the wild-type pYpY phosphopeptide, aiming to simultaneously improve affinity for these kinases and to maximize the affinity gap between different kinases. Consequently, a population is computationally evolved as selective phosphopeptide candidates; the dissociation constants of four representatives with HER kinases are systematically determined using binding affinity analysis, from which their selectivity is derived. The designed pYpYp3 phosphopeptide possesses a high selectivity over different HER kinases (S = 4.8-fold) and satisfactory affinity profile to these kinase (K_D = 140-1000 μM). Structural analysis observes that the global binding modes of pYpYp3 to different kinases are roughly consistent, but its local conformation may vary considerably, thus conferring specificity to the phosphopeptide.