J Biomol Struct Dyn. 2026 Jun 02.
1-35
Phosphoinositide 3-kinase (PI3K) isoform selectivity remains a central challenge in cancer therapeutics, as highly conserved ATP-binding domains hinder the development of specific inhibitors while broad inhibition leads to severe toxicities. Here, we introduce a computational framework that integrates molecular dynamics (MD), residue-network analysis, and hydration thermodynamics to define the molecular determinants of PI3K isoform selectivity. Across 100-ns all-atom simulations of class I PI3K isoforms (α, β, γ, δ), we identify distinct 'dynamic signatures' rather than static pocket differences. PI3Kα and PI3Kβ exhibit the greatest structural stability, whereas PI3Kδ shows the largest conformational deviation, and PI3Kγ displays exceptional plasticity, reflecting a unique mechanism of dynamic regulation. Network analysis highlights isoform-specific supercritical hubs (SER838 in α, SER841 in β, SER815 in δ, and ASP849 in γ) that are predicted to act as bottlenecks of allosteric communication. Importantly, druggability emerges independently of pocket size; although PI3Kβ has the largest cavity, PI3Kγ is predicted to be the most druggable, consistent with favorable water displacement around hydrophobic rims and thermodynamically challenging hydration near polar hubs. Hydrogen bond residence time analysis independently validated the hydration site classifications, revealing a strong correlation (R2 = 0.91) between water residence times and hydration site free energies across all 20 supercritical residues. Together, these findings suggest that isoform selectivity may be primarily governed by dynamic ensembles, allosteric pathways, and water-mediated interactions, establishing a computational hypothesis framework to guide the future rational design of isoform-selective PI3K inhibitors that target dynamic rather than purely static features.
Keywords: PI3K isoform selectivity; allosteric communication; conformational dynamics; molecular dynamics simulations; selective inhibitors