Binding selectivity studies of phosphoinositide 3-kinases using free energy calculations

Phosphoinositide 3-kinases (PI3Ks) and their phosphatidylinositol 3,4,5-triphosphate (PIP(3)) products regulate a variety of cellular processes. Of these, PI3Kα is an attractive target for anticancer drug design. Mutations in the PI3Kα kinase domain alter the mobility of the activation loop resulting in gain of function. We employed molecular dynamics (MD) simulations-based energetic analysis using molecular mechanics/generalized born surface area (MM/GBSA) for PI3Kα and -γ. MD simulations were carried out for PI3K models based on the RESP (restrained electrostatic potential) and quantum mechanics (QM)-polarized ligand docking (QPLD)-derived partial charges. Computational alanine scanning was also used to evaluate the contributions of key binding residues to ligand binding. Our results show that both RESP and QPLD charge models of PI3Kα and PI3Kγ provide similar performance in MD simulations. For example, the predicted RESP and QPLD free energies of -9.5 and -9.3 kcal/mol for LY294002 binding to PI3Kγ and -10.9 and -11.7 kcal/mol for wortmannin binding to PI3Kα are in good agreement with experimental values. A significant loss in binding free energy was observed when hydrophobic residues were mutated to alanine, suggesting that specific hydrophobic interactions are important to optimal ligand binding. MM/GBSA calculations suggested that residues Ser774, Gln859, and Ile932 of PI3Kα might be used to design H1047R mutant-specific ligands, whereas Lys890 of PI3Kγ can be used for ligand design targeting PI3Kγ.