The PARP family consists of 17 members with diverse functions, including those related to cancer cells’ viability. Several PARP inhibitors are of great interest as innovative anticancer drugs, but they have low selectivity towards distinct PARP family members and exert serious adverse effects. We describe a family-wide study of the nicotinamide (NA) binding site, an important functional region in the PARP structure, using comparative bioinformatic analysis and molecular modeling. Mutations in the NA site and D-loop mobility around the NA site were identified as factors that can guide the design of selective PARP inhibitors. Our findings are of particular importance for the development of novel tankyrase (PARPs 5a and 5b) inhibitors for cancer therapy.

Poly(ADP-ribose) polymerase 1 (PARP-1) is a key DNA repair enzyme and an important target in cancer treatment. Conventional methods of studying the reaction mechanism of PARP-1 have limitations because of the complex structure of PARP-1 substrates; however, the necessary data can be obtained by molecular modeling. In this work, a molecular dynamics model for the PARP-1 enzyme-substrate complex containing NAD⁺ molecule and the end of the poly(ADP-ribose) chain in the form of ADP molecule was obtained for the first time. Interactions with the active site residues have been characterized where Gly863, Lys903, Glu988 play a crucial role, and the S_N1-like mechanism for the enzymatic ADP-ribosylation reaction has been proposed. Models of PARP-1 complexes with more sophisticated two-unit fragments of the growing polymer chain as well as competitive inhibitors 3-aminobenzamide and 7-methylguanine have been obtained by molecular docking.