Structural Determinants of PARP1 Selectivity from Molecular Dynamics Analysis of PARP1 and PARP2 Complexes

Selective inhibition of poly(ADP-ribose) polymerase 1 (PARP1) may reduce the hematologic toxicity associated with dual PARP1/PARP2 inhibition. We performed molecular dynamics simulations for five selective inhibitors in complexes with PARP1 and PARP2, using three independent 50 ns runs per complex after docking and equilibration, followed by protein–ligand interaction fingerprint and statistical analyses. All complexes remained dynamically stable, with ligand root-mean-square deviation values generally within 0.3 nm. Comparative analysis identified three αF-helix residue pairs with nominally reduced interaction frequencies in PARP2: Asn767/Ala336, Leu769/Gly338, and Asp770/Asp339 (p < 0.05). After Benjamini–Hochberg correction for multiple comparisons, Leu769/Gly338 remained significant (q < 0.05), indicating that this pair represents the most statistically robust interaction difference within this region. Using palacaparib as the most selective inhibitor, these differences were associated with weakened or lost hydrophobic, van der Waals, and cation–π interactions in PARP2. Selective binding of modern PARP1 inhibitors appears to be associated with αF-helix-dependent interaction patterns, providing a mechanistic basis for the rational design of next-generation selective inhibitors with improved selectivity and potentially reduced toxicity.