Next Article in Journal
Current Status of Management and Outcome for Patients with Ewing Sarcoma
Next Article in Special Issue
Candidate Markers of Olaparib Response from Genomic Data Analyses of Human Cancer Cell Lines
Previous Article in Journal
Subsequent Malignant Neoplasms in Retinoblastoma Survivors
Previous Article in Special Issue
The Paradoxical Effect of PARP Inhibitor BGP-15 on Irinotecan-Induced Cachexia and Skeletal Muscle Dysfunction
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 6
10.3390/cancers13061201
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bioinformatic Analysis of the Nicotinamide Binding Site in Poly(ADP-Ribose) Polymerase Family Proteins

by
Garri Manasaryan
1,
Dmitry Suplatov
2,
Sergey Pushkarev
3,
Viktor Drobot
2,
Alexander Kuimov
2,
Vytas Švedas
2,3 and
Dmitry Nilov
2,*
1
Faculty of Medicine, Lomonosov Moscow State University, Lomonosov Ave. 27, bldg. 1, 119991 Moscow, Russia
2
Belozersky Institute of Physicochemical Biology, Lomonosov Moscow State University, Lenin Hills 1, bldg. 40, 119991 Moscow, Russia
3
Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, Lenin Hills 1, bldg. 73, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(6), 1201; https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13061201
Submission received: 2 February 2021 / Revised: 2 March 2021 / Accepted: 5 March 2021 / Published: 10 March 2021
(This article belongs to the Special Issue PARP Inhibitors for Cancer Therapy)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The PARP family consists of 17 proteins, and some of them are responsible for cancer cells’ viability. Much attention is therefore given to the search for chemical compounds with the ability to suppress distinct PARP family members (for example, PARP-5a and 5b). Here, we present the results of a family-wide bioinformatic analysis of an important functional region in the PARP structure and describe factors that can guide the design of highly selective compounds.

Abstract

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.

1. Introduction

Poly(ADP-Ribose) polymerase proteins (PARPs 1–16) catalyze the transfer of ADP-ribose from the nicotinamide adenine dinucleotide (NAD+) substrate to target proteins and are involved in many cellular functions (Table 1) [1,2,3,4,5,6,7]. In particular, the activity of the founding family member PARP-1 at DNA damaged sites recruits the base excision repair proteins XRCC1, DNA polymerase β, and DNA ligase III [8,9,10]. The most studied PARP family members are PARP-1 and 2 involved in DNA repair and PARP-5a and 5b (also known as tankyrases 1 and 2) regulating the Wnt signaling pathway. In the past years, aberrant Wnt signaling was found to be implicated in numerous malignancies, including gastrointestinal cancers, leukemia, and breast cancer [11,12].
Much attention is being paid to the search for inhibitors of distinct PARP family members because of their role in breast/ovarian cancers [29,30,31,32,33], gastrointestinal cancers, and many other, non-oncological diseases [34,35,36]. A number of PARP inhibitors contain an amide group attached to an aromatic ring or lactam group built in an aromatic ring system and mimic the nicotinamide (NA) moiety of NAD+. They suppress enzyme activity at nanomolar concentrations but lack sufficient selectivity [37,38,39], for example, isoquinolinone derivatives [40] are able to bind to PARPs 1, 3, and 5a. FDA-approved inhibitors of PARP-1/2 (olaparib, rucaparib, and niraparib) [41,42] can cause serious side effects likely related to the nonselective interaction with numerous NAD+-binding proteins, including other PARP family members [43,44,45,46].
The design of highly effective inhibitors is usually focused on a region of the PARP active site responsible for binding the NA group of NAD+. The architecture of the NA site provides strong interaction between the substrate/inhibitor amide group and Gly863 residue (PARP-1 numbering) via two hydrogen bonds [47,48,49,50,51]. In many crystal structures, PARP inhibitors form additional interactions with the other NA site residues: A hydrogen bond with Ser904, hydrophobic contact with Ala898, and π-stacking with Tyr907 [39,52,53]. The above residues are crucial for the effective binding in the PARP active site, and their mutations may modulate the affinity for inhibitors. Interactions with flexible loops around the NA site may also be important for molecular recognition. In this article, we present the results of a detailed family-wide analysis of the NA binding site in PARPs 1–16 that opens up new prospects for selective PARP inhibition.

2. Results

2.1. Multiple Sequence Alignment of PARPs 1–16

A multiple sequence alignment highlighting residues of the NA binding site in PARPs 1–16 is shown in Figure 1. Gly is the most prevalent residue at position 863, Ala at position 898, Ser at position 904, and Tyr at position 907 (PARP-1 numbering). Substitutions at these positions, marked in Table 2, may have a direct impact on the substrate and inhibitor binding. Notably, the key active site residue Gly863 is substituted only in catalytically inactive PARP family members, PARP-9 and PARP-13. Crystal structures of PARPs 1–3, 5a, 5b, 10, and 12–16 are available, which helps considerably in understanding how substitutions affect the binding site architecture, and the models of PARPs 4, 6–9, and 11 can be constructed from their sequences and structures of close homologs. In Section 2.2 and Section 2.3, interactions of PARPs mediated by amino acid residues at positions 863, 898, 904, and 907 in the multiple alignment are analyzed in detail using molecular dynamics (MD) and homology modeling.

2.2. Modeling of PARPs with Known Structures

To model the NA binding site of PARPs 1–3, 5a, 5b, 10, and 12–16, we have selected representative crystal structures for each family member. Available structures listed in Table S1 were manually clustered into groups based on the similarity of the NA site conformation (outliers/minor conformations were excluded), and then a representative structure with the best resolution was selected from each cluster (Figure S1, Table S1). Single clusters were produced for all PARPs except PARP-5a, whose structures were divided into groups I and II differing in the D-loop conformation. The mobile D-loop of PARP-5a is located around the NA site [63] and may be involved in inhibitor binding. The D-loop residue Tyr1203 is oriented towards the NA site in conformation I, and Phe1208 in conformation II (Figure S2). The phenyl ring of the D-loop residue occupies a similar position in conformations I and II, indicating that the D-loop can mediate nonpolar interactions with NA site ligands in both possible conformational states.
The representative crystal structures of PARPs shown in Figure 2 provide a good starting point for modeling intermolecular interactions with NA mimics and probing the selectivity. We have chosen 7-methylguanine (7-MG) as a probe PARP inhibitor for the following reasons. (i) 7-MG is a promising competitive inhibitor of the founding family member PARP-1 [64,65,66,67], (ii) 7-MG is a small NA mimic that occupies only the NA binding region, forming crucial interactions with the Gly863, Ala898, Ser904, and Tyr907 residues [51,65], (iii) high-quality force field parameters of 7-MG are available and ready for use in molecular modeling [64]. The 7-MG molecule was docked into the NA binding site of PARPs, and then its position was refined using MD simulation and analyzed.
The overall structure of the catalytic domain is similar in PARPs 1–3 and includes an (ADP-ribosyl) transferase subdomain and a regulatory helical subdomain [68]. The residues of the NA binding site are conserved in these proteins and form interactions with 7-MG typical for PARP inhibitors: Two hydrogen bonds between Gly863/429/385 (PARP-1/2/3) and the lactam group, which act simultaneously as a donor and an acceptor, a hydrogen bond between Ser904/470/422 and the lactam oxygen, hydrophobic contact between Ala898/464/416 and the 7-methyl group, π-stacking between Tyr907/473/425 and the 7-MG fused rings (Figure 3, Figure S3 showing a close-up view of the NA site, and Table S2 providing mean distances and angles from 10 ns simulation). It is, therefore, not surprising that 7-MG inhibits PARP-1 and 2 with comparable potency [64].
In PARPs 5a and 5b, also known as tankyrases, the NA site residues are the same as in PARPs 1–3 and mediate similar interactions with 7-MG (Figure 4, Table S2), but additional intermolecular contacts may be formed due to the D-loop mobility. In the starting model of PARP-5a conformation I, the Tyr1203 side chain contacted with an aromatic nitrogen of 7-MG, but during the MD simulation, it moved away (without affecting the main chain coordinates). Meanwhile, the unfavorable contact with the Phe1208 side chain was persistent in conformation II (Figure 5). The starting D-loop position in the PARP-5b model was similar to conformation I of PARP-5a, and the corresponding Tyr1050 residue also moved away from 7-MG during the simulation. Both PARP-5a and 5b models indicate that D-loop interactions with donor/acceptor atoms of NA mimics are unfavorable, in accordance with our preliminary experimental data demonstrating the 7-MG selectivity for PARP-1 over PARP-5b (Table S3) [69]. However, more hydrophobic inhibitors may be selective against tankyrases due to additional contacts with Tyr1203/1050 or Phe1208. For example, the chloro substituent of an aromatic inhibitor interacts with the Tyr1050 side chain in the PARP-5b structure 4j1z [70].
In PARP-10 and 12, the NA site residues are also conserved (Figure 6). However, the hydrogen bond with the Ser604 residue in the active site of PARP-12 was characterized by an increased mean distance (Table S2) because its side chain periodically formed an alternative bond with the Asp600 main chain. In PARP-13 the crucial Gly residue is replaced with Ala788 whose methyl group is oriented away from the inhibitor and does not significantly affect the formation of hydrogen bonds between the main chain and lactam group. The Tyr residue is replaced with Asn830 in the NA site of PARP-13, which results in loss of the stacking interaction (Figure 7).
The NA site residues in PARP-14 and 15 are the same as in PARPs 1–3 (Figure 8). However, during the simulation of PARP-15 an additional hydrogen bond was formed between the aromatic nitrogen of 7-MG and the D-loop residue Gly558 (Figure 9). In PARP-16 one hydrogen bond with the lactam group is lost due to the Ser replacement by Ala190, while the replacement of another NA site residue Ala by Thr184 does not affect the ability to form favorable contacts with the inhibitor’s hydrophobic group (Figure 10, Table S2).

2.3. Modeling of PARPs with Unknown Structures

The NA binding site architecture in PARPs of unknown structure (4, 6–9, and 11) was predicted and studied using homology modeling (template structures are listed in Table S4). We consider such an approach reasonable because the functional sites of the protein tend to be more conserved in evolution (than the rest of the fold) and thus more accurately modeled [71]. The coordinates of 7-MG were transferred into the NA site of the obtained PARP models and then optimized. The modeled fold of PARP-4 was similar to its homologs PARPs 1–3, and the replacement of Ala with Ser479 did not affect significantly the NA site structure because the Ser479 side chain was oriented away from the inhibitor (Figure 11, Table S5). Similarly to PARP-4, PARP-6, and PARP-8 contain non-essential substitutions of Ala by Ser510 and Ser734, respectively (Figure 11).
In PARP-7 and 11, the replacement of Tyr by Phe575/Phe247 does not disrupt π-stacking because this interaction with the inhibitor can be mediated by a phenyl ring of both Tyr and Phe residues (Figure 12, Table S5). Lastly, the modeled PARP-9 protein is a PARP family member in which all the NA site residues are substituted (Figure 13). Gly is replaced by Gln706 whose side chain is oriented away from the inhibitor and does not affect the formation of hydrogen bonds between the main chain and lactam group. Ser is replaced by Leu745, which results in the loss of an additional bond with the lactam group. Similarly to PARP-16, Ala is replaced by Thr739, which does not significantly affect the hydrophobic interaction. The stacking interaction with the inhibitor is lost in PARP-9 due to the replacement of Tyr with Lys750.

3. Discussion

PARP inhibitors represent a promising new component of cancer chemotherapy, and a lot of attention is given to the design of compounds selective towards distinct PARP family members. In particular, PARP-1 and 2 are involved in the elimination of single- and double-strand breaks of DNA, and tankyrases (PARP-5a and 5b) modify the axin protein in the Wnt pathway to promote cell proliferation. Therefore, PARP-1/2 inhibitors can be effective against BRCA-deficient tumors [72,73], whilst tankyrase inhibitors against tumors with abnormal Wnt signaling [74,75,76,77]. Here, we present a detailed study of the NA site, crucial for substrate and inhibitor binding, in PARP family proteins. Bioinformatic analysis and molecular modeling with a probe 7-MG inhibitor allowed us to identify structural features of the NA site in PARPs 1–16 important for the selective binding. It should be noted that a recently published paper by Kam et al., also dedicated to an in silico family-wide analysis of PARPs, is mostly focused on new potential targets for inhibition beyond the NA site [78], and therefore, our findings do not conflict with or duplicate previously reported data.
The most prevalent NA site residue at position 863 (PARP-1 numbering) is Gly which forms key hydrogen bonds with the inhibitor’s lactam group in the obtained PARP models. The replacement of Gly in PARP-9 and 13 (with Gln and Ala, respectively) does not significantly affect the formation of these main-chain hydrogen bonds. Ala, Ser, or Thr at position 898 provide favorable contacts with the inhibitor’s hydrophobic group. Ser at position 904 forms an additional hydrogen bond with the lactam group, except PARP-9 and 16 where it is replaced by Leu and Ala, respectively. Tyr or Phe at position 907 form π-stacking interactions, except PARP-9 and 13 containing substitutions with Lys and Asn, respectively. Thus, substitutions at positions 904 and 907 presumably provide selectivity of NA mimics for PARPs 1–8, 10–12, 14, and 15 over PARPs 9, 13, and 16. The mobility of the D-loop around the NA site of PARPs 5a, 5b and 15 can result in additional interactions: hydrophobic (PARP-5a, 5b) or polar (PARP-15). This may be exploited in the design of selective inhibitors of PARP-1/2 or PARP-5a/5b: more hydrophobic NA mimics would tend to bind with PARP-5a and 5b, while more polar compounds, forming unfavorable contacts with D-loop, would preferentially target PARP-1 and 2. As surgery remains the primary modality of cure in cancers associated with aberrant Wnt signaling, additional targeted treatments with selective PARP-5a/5b inhibitors may be of great interest [79,80].

4. Materials and Methods

Amino acid sequences of the catalytic domain of PARPs 1–16 were obtained from the UniProt database [81]. To construct a multiple alignment of PARPs, various state-of-the-art methods (COBALT, PROMALS3D, Matt, Clustal Omega) were independently used, followed by manual expert evaluation, and the alignment by Clustal Omega [82,83,84] was found to be the most accurate due to correct superimposition of key NA site residues.
Crystal structures of the catalytic domain of PARPs 1–3, 5a, 5b, 10, and 12–16 were obtained from the Protein Data Bank [85], superimposed with Matt 1.0 [86], and manually clustered based on the NA binding site conformation. Molecular models of PARPs 1–3, 5a (conformations I and II), 5b, 10, and 12–16 were then constructed based on the selected representative crystal structures. The coordinates of missing residues of the PARP catalytic domain were predicted with Modeller 9.20 or transferred from other structures (Table S6) [87,88]. Next, the protein structures were optimized with AmberTools 15 and Amber 14 [89,90], according to the following protocol. Hydrogen atoms were added to the structure considering ionization of amino acid residues, and then it was solvated by 12 Å-thick layer of TIP3P water. Chloride or sodium ions were added to neutralize the net charge. The energy minimization was performed with positional restraints on heavy atoms of the protein (2500 steepest descent steps + 2500 conjugate gradient steps). The 7-MG molecule was docked into the active site of the PARP models with Lead Finder 1.1.16 [91,92]. The obtained PARP–7-MG complexes were re-optimized in two stages, one with positional restraints on the protein and inhibitor atoms (2500 steepest descent steps + 2500 conjugate gradient steps) and the other without restraints (5000 steepest descent steps + 5000 conjugate gradient steps). The system was then heated up from 0 to 300 K (50 ps, constant volume) and equilibrated at 300 K (500 ps, constant pressure). Lastly, 10 ns trajectory of equilibrium simulation was calculated and analyzed. Structures of PARPs 4, 6–9, and 11 were obtained using homology modeling. Close homologs of these family members were identified with a HHpred server [93,94] and used as templates for Modeller 9.20. 7-MG coordinates were transferred from the docking model of PARP-1 complex, and the obtained homology models of PARP–7-MG complexes were then energy minimized in two stages, as described above.
Control data for energy minimization and MD simulation are provided in Table S7. The ff14SB force field [95] was used to describe the protein with molecular mechanics, and recently developed parameters [64] were used to describe the 7-MG molecule. VMD 1.9.2 was used for the visualization of structures [96].

5. Conclusions

The present paper systematically describes the architecture of the NA binding site in 17 PARP family proteins (PARPs 1–4, 5a, 5b, 6–16) and can serve as a useful guide to estimate the selectivity of NA mimics towards distinct family members. Certain factors may lead to the selective inhibition: (i) Mutations in the NA site and (ii) D-loop mobility around the NA site. An important finding of our study is that only in tankyrases (PARP-5a and 5b) the mobile D-loop can form additional hydrophobic contacts with NA mimics, which provides opportunities for the development of highly selective tankyrase inhibitors as promising anticancer agents.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6694/13/6/1201/s1, Figure S1: Cluster of similar conformations of the NA binding site in PARP-1 crystal structures, Figure S2: Two possible conformations of the D-loop in crystal structures of PARP-5a, Figure S3: Interactions of 7-MG in the NA binding site of PARP-1 revealed by molecular modeling, Table S1: Crystal structures of PARPs used in the analysis of the NA binding site architecture, Table S2: Interactions between a probe inhibitor (7-MG) and NA site residues in PARPs revealed by 10-ns MD simulation, Table S3: Activity of PARP-1 and PARP-5b (tankyrase 2) at 7-MG concentration of 360 µM determined with an immunochemical assay, Table S4: PARPs of unknown structure and their close homologues, Table S5: Interactions between a probe inhibitor (7-MG) and NA site residues in PARPs revealed using homology modeling, Table S6: Missing residues in representative PARP structures, TableS7: Control data used for energy minimization and MD simulation of the PARP–7-MG complexes.

Author Contributions

Conceptualization and funding acquisition, D.N.; investigation, G.M., D.S., S.P., and V.D.; writing—original draft preparation, G.M. and D.N.; supervision and writing—review and editing, A.K. and V.Š. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation, grant number 19-74-10072.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

7-MG7-methylguanine
MDmolecular dynamics
NAnicotinamide
NAD+nicotinamide adenine dinucleotide
PARPpoly(ADP-ribose)polymerase

References

  1. Amé, J.C.; Spenlehauer, C.; de Murcia, G. The PARP superfamily. Bioessays 2004, 26, 882–893. [Google Scholar] [CrossRef] [PubMed]
  2. Otto, H.; Reche, P.A.; Bazan, F.; Dittmar, K.; Haag, F.; Koch-Nolte, F. In silico characterization of the family of PARP-like poly(ADP-ribosyl)transferases (pARTs). BMC Genom. 2005, 6, 139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Vyas, S.; Chesarone-Cataldo, M.; Todorova, T.; Huang, Y.H.; Chang, P. A systematic analysis of the PARP protein family identifies new functions critical for cell physiology. Nat. Commun. 2013, 4, 2240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Bai, P. Biology of poly(ADP-Ribose) polymerases: The factotums of cell maintenance. Mol. Cell 2015, 58, 947–958. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Barkauskaite, E.; Jankevicius, G.; Ahel, I. Structures and mechanisms of enzymes employed in the synthesis and degradation of PARP-dependent protein ADP-ribosylation. Mol. Cell 2015, 58, 935–946. [Google Scholar] [CrossRef] [Green Version]
  6. Kuny, C.V.; Sullivan, C.S. Virus-host interactions and the ARTD/PARP family of enzymes. PLoS Pathog. 2016, 12, e1005453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Gupte, R.; Liu, Z.; Kraus, W.L. PARPs and ADP-ribosylation: Recent advances linking molecular functions to biological outcomes. Genes Dev. 2017, 31, 101–126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Masson, M.; Niedergang, C.; Schreiber, V.; Muller, S.; Menissier-de Murcia, J.; de Murcia, G. XRCC1 is specifically associated with poly(ADP-ribose) polymerase and negatively regulates its activity following DNA damage. Mol. Cell Biol. 1998, 18, 3563–3571. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Brem, R.; Hall, J. XRCC1 is required for DNA single-strand break repair in human cells. Nucleic Acids Res. 2005, 33, 2512–2520. [Google Scholar] [CrossRef] [Green Version]
  10. Alemasova, E.E.; Lavrik, O.I. Poly(ADP-ribosyl)ation by PARP1: Reaction mechanism and regulatory proteins. Nucleic Acids Res. 2019, 47, 3811–3827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Duchartre, Y.; Kim, Y.M.; Kahn, M. The Wnt signaling pathway in cancer. Crit. Rev. Oncol. Hematol. 2016, 99, 141–149. [Google Scholar] [CrossRef] [Green Version]
  12. Zhan, T.; Rindtorff, N.; Boutros, M. Wnt signaling in cancer. Oncogene 2017, 36, 1461–1473. [Google Scholar] [CrossRef] [PubMed]
  13. Schreiber, V.; Amé, J.C.; Dollé, P.; Schultz, I.; Rinaldi, B.; Fraulob, V.; Ménissier-de Murcia, J.; de Murcia, G. Poly(ADP-ribose) polymerase-2 (PARP-2) is required for efficient base excision DNA repair in association with PARP-1 and XRCC1. J. Biol. Chem. 2002, 277, 23028–23036. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Schreiber, V.; Dantzer, F.; Ame, J.C.; de Murcia, G. Poly(ADP-ribose): Novel functions for an old molecule. Nat. Rev. Mol. Cell Biol. 2006, 7, 517–528. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, Z.; Wang, F.; Tang, T.; Guo, C. The role of PARP1 in the DNA damage response and its application in tumor therapy. Front. Med. 2012, 6, 156–164. [Google Scholar] [CrossRef]
  16. Hassler, M.; Ladurner, A.G. Towards a structural understanding of PARP1 activation and related signalling ADP-ribosyl-transferases. Curr. Opin. Struct. Biol. 2012, 22, 721–729. [Google Scholar] [CrossRef]
  17. Zarkovic, G.; Belousova, E.A.; Talhaoui, I.; Saint-Pierre, C.; Kutuzov, M.M.; Matkarimov, B.T.; Biard, D.; Gasparutto, D.; Lavrik, O.I.; Ishchenko, A.A. Characterization of DNA ADP-ribosyltransferase activities of PARP2 and PARP3: New insights into DNA ADP-ribosylation. Nucleic Acids Res. 2018, 46, 2417–2431. [Google Scholar] [CrossRef] [PubMed]
  18. Grundy, G.J.; Polo, L.M.; Zeng, Z.; Rulten, S.L.; Hoch, N.C.; Paomephan, P.; Xu, Y.; Sweet, S.M.; Thorne, A.W.; Oliver, A.W.; et al. PARP3 is a sensor of nicked nucleosomes and monoribosylates histone H2B(Glu2). Nat. Commun. 2016, 7, 12404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Rodriguez-Vargas, J.M.; Nguekeu-Zebaze, L.; Dantzer, F. PARP3 comes to light as a prime target in cancer therapy. Cell Cycle 2019, 18, 1295–1301. [Google Scholar] [CrossRef]
  20. Huang, S.M.; Mishina, Y.M.; Liu, S.; Cheung, A.; Stegmeier, F.; Michaud, G.A.; Charlat, O.; Wiellette, E.; Zhang, Y.; Wiessner, S.; et al. Tankyrase inhibition stabilizes axin and antagonizes Wnt signalling. Nature 2009, 461, 614–620. [Google Scholar] [CrossRef] [PubMed]
  21. Mariotti, L.; Pollock, K.; Guettler, S. Regulation of Wnt/β-catenin signalling by tankyrase-dependent poly(ADP-ribosyl)ation and scaffolding. Br. J. Pharmacol. 2017, 174, 4611–4636. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Ferri, M.; Liscio, P.; Carotti, A.; Asciutti, S.; Sardella, R.; Macchiarulo, A.; Camaioni, E. Targeting Wnt-driven cancers: Discovery of novel tankyrase inhibitors. Eur. J. Med. Chem. 2017, 142, 506–522. [Google Scholar] [CrossRef] [PubMed]
  23. Kim, M.K. Novel insight into the function of tankyrase. Oncol. Lett. 2018, 16, 6895–6902. [Google Scholar] [CrossRef] [Green Version]
  24. Qi, G.; Kudo, Y.; Tang, B.; Liu, T.; Jin, S.; Liu, J.; Zuo, X.; Mi, S.; Shao, W.; Ma, X.; et al. PARP6 acts as a tumor suppressor via downregulating Survivin expression in colorectal cancer. Oncotarget 2016, 7, 18812–18824. [Google Scholar] [CrossRef] [Green Version]
  25. Zhao, Y.; Hu, X.; Wei, L.; Song, D.; Wang, J.; You, L.; Saiyin, H.; Li, Z.; Yu, W.; Yu, L.; et al. PARP10 suppresses tumor metastasis through regulation of Aurora A activity. Oncogene 2018, 37, 2921–2935. [Google Scholar] [CrossRef]
  26. Li, L.; Zhao, H.; Liu, P.; Li, C.; Quanquin, N.; Ji, X.; Sun, N.; Du, P.; Qin, C.F.; Lu, N.; et al. PARP12 suppresses Zika virus infection through PARP-dependent degradation of NS1 and NS3 viral proteins. Sci. Signal. 2018, 11, eaas9332. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Todorova, T.; Bock, F.J.; Chang, P. Poly(ADP-ribose) polymerase-13 and RNA regulation in immunity and cancer. Trends Mol. Med. 2015, 21, 373–384. [Google Scholar] [CrossRef] [Green Version]
  28. Barbarulo, A.; Iansante, V.; Chaidos, A.; Naresh, K.; Rahemtulla, A.; Franzoso, G.; Karadimitris, A.; Haskard, D.O.; Papa, S.; Bubici, C. Poly(ADP-ribose) polymerase family member 14 (PARP14) is a novel effector of the JNK2-dependent pro-survival signal in multiple myeloma. Oncogene 2013, 32, 4231–4242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Cepeda, V.; Fuertes, M.A.; Castilla, J.; Alonso, C.; Quevedo, C.; Soto, M.; Pérez, J.M. Poly(ADP-ribose) polymerase-1 (PARP-1) inhibitors in cancer chemotherapy. Recent Pat. Anticancer Drug Discov. 2006, 1, 39–53. [Google Scholar] [CrossRef]
  30. Martin, S.A.; Lord, C.J.; Ashworth, A. DNA repair deficiency as a therapeutic target in cancer. Curr. Opin. Genet. Dev. 2008, 18, 80–86. [Google Scholar] [CrossRef]
  31. Ferraris, D.V. Evolution of poly(ADP-ribose) polymerase-1 (PARP-1) inhibitors. From concept to clinic. J. Med. Chem. 2010, 53, 4561–4584. [Google Scholar] [CrossRef]
  32. Efremova, A.S.; Zakharenko, A.L.; Shram, S.I.; Kulikova, I.V.; Drenichev, M.S.; Sukhanova, M.V.; Khodyreva, S.N.; Myasoedov, N.F.; Lavrik, O.I.; Mikhailov, S.N. Disaccharide pyrimidine nucleosides and their derivatives: A novel group of cell-penetrating inhibitors of poly(ADP-ribose) polymerase 1. Nucleosides Nucleotides Nucleic Acids 2013, 32, 510–528. [Google Scholar] [CrossRef] [PubMed]
  33. Lord, C.J.; Tutt, A.N.; Ashworth, A. Synthetic lethality and cancer therapy: Lessons learned from the development of PARP inhibitors. Annu. Rev. Med. 2015, 66, 455–470. [Google Scholar] [CrossRef]
  34. Curtin, N.J.; Szabo, C. Therapeutic applications of PARP inhibitors: Anticancer therapy and beyond. Mol. Asp. Med. 2013, 34, 1217–1256. [Google Scholar] [CrossRef] [Green Version]
  35. Henning, R.J.; Bourgeois, M.; Harbison, R.D. Poly(ADP-ribose) polymerase (PARP) and PARP inhibitors: Mechanisms of action and role in cardiovascular disorders. Cardiovasc. Toxicol. 2018, 18, 493–506. [Google Scholar] [CrossRef] [PubMed]
  36. Berger, N.A.; Besson, V.C.; Boulares, A.H.; Bürkle, A.; Chiarugi, A.; Clark, R.S.; Curtin, N.J.; Cuzzocrea, S.; Dawson, T.M.; Dawson, V.L.; et al. Opportunities for the repurposing of PARP inhibitors for the therapy of non-oncological diseases. Br. J. Pharmacol. 2018, 175, 192–222. [Google Scholar] [CrossRef] [PubMed]
  37. Wahlberg, E.; Karlberg, T.; Kouznetsova, E.; Markova, N.; Macchiarulo, A.; Thorsell, A.G.; Pol, E.; Frostell, Å.; Ekblad, T.; Öncü, D.; et al. Family-wide chemical profiling and structural analysis of PARP and tankyrase inhibitors. Nat. Biotechnol. 2012, 30, 283–288. [Google Scholar] [CrossRef] [PubMed]
  38. Ekblad, T.; Camaioni, E.; Schüler, H.; Macchiarulo, A. PARP inhibitors: Polypharmacology versus selective inhibition. FEBS J. 2013, 280, 3563–3575. [Google Scholar] [CrossRef]
  39. Jain, P.G.; Patel, B.D. Medicinal chemistry approaches of poly ADP-Ribose polymerase 1 (PARP1) inhibitors as anticancer agents—A recent update. Eur. J. Med. Chem. 2019, 165, 198–215. [Google Scholar] [CrossRef]
  40. Chiarugi, A.; Meli, E.; Calvani, M.; Picca, R.; Baronti, R.; Camaioni, E.; Costantino, G.; Marinozzi, M.; Pellegrini-Giampietro, D.E.; Pellicciari, R.; et al. Novel isoquinolinone-derived inhibitors of poly(ADP-ribose) polymerase-1: Pharmacological characterization and neuroprotective effects in an in vitro model of cerebral ischemia. J. Pharmacol. Exp. Ther. 2003, 305, 943–949. [Google Scholar] [CrossRef]
  41. Mittica, G.; Ghisoni, E.; Giannone, G.; Genta, S.; Aglietta, M.; Sapino, A.; Valabrega, G. PARP inhibitors in ovarian cancer. Recent Pat. Anticancer Drug Discov. 2018, 13, 392–410. [Google Scholar] [CrossRef]
  42. Zimmer, A.S.; Gillard, M.; Lipkowitz, S.; Lee, J.M. Update on PARP inhibitors in breast cancer. Curr. Treat. Options Oncol. 2018, 19, 21. [Google Scholar] [CrossRef]
  43. Malyuchenko, N.V.; Kotova, E.Y.; Kulaeva, O.I.; Kirpichnikov, M.P.; Studitskiy, V.M. PARP1 inhibitors: Antitumor drug design. Acta Nat. 2015, 7, 27–37. [Google Scholar] [CrossRef]
  44. Sonnenblick, A.; de Azambuja, E.; Azim, H.A., Jr.; Piccart, M. An update on PARP inhibitors–moving to the adjuvant setting. Nat. Rev. Clin. Oncol. 2015, 12, 27–41. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, Y.Q.; Wang, P.Y.; Wang, Y.T.; Yang, G.F.; Zhang, A.; Miao, Z.H. An update on poly(ADP-ribose)polymerase-1 (PARP-1) inhibitors: Opportunities and challenges in cancer therapy. J. Med. Chem. 2016, 59, 9575–9598. [Google Scholar] [CrossRef]
  46. Ohmoto, A.; Yachida, S. Current status of poly(ADP-ribose) polymerase inhibitors and future directions. Onco Targets Ther. 2017, 10, 5195–5208. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Ruf, A.; Mennissier de Murcia, J.; de Murcia, G.; Schulz, G.E. Structure of the catalytic fragment of poly(AD-ribose) polymerase from chicken. Proc. Natl. Acad. Sci. USA 1996, 93, 7481–7485. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Jagtap, P.; Szabó, C. Poly(ADP-ribose) polymerase and the therapeutic effects of its inhibitors. Nat. Rev. Drug Discov. 2005, 4, 421–440. [Google Scholar] [CrossRef] [PubMed]
  49. Zakharenko, A.L.; Sukhanova, M.V.; Khodyreva, S.N.; Novikov, F.N.; Stroylov, V.S.; Nilov, D.K.; Chilov, G.G.; Švedas, V.K.; Lavrik, O.I. Improved procedure of the search for poly(ADP-Ribose) polymerase-1 potential inhibitors with the use of the molecular docking approach. Mol. Biol. (Mosc.) 2011, 45, 517–521. [Google Scholar] [CrossRef]
  50. Nilov, D.K.; Yashina, K.I.; Gushchina, I.V.; Zakharenko, A.L.; Sukhanova, M.V.; Lavrik, O.I.; Švedas, V.K. 2,5-Diketopiperazines: A new class of poly(ADP-ribose)polymerase inhibitors. Biochemistry (Mosc.) 2018, 83, 152–158. [Google Scholar] [CrossRef] [PubMed]
  51. Nilov, D.K.; Pushkarev, S.V.; Gushchina, I.V.; Manasaryan, G.A.; Kirsanov, K.I.; Švedas, V.K. Modeling of the enzyme–substrate complexes of human poly(ADP-ribose) polymerase 1. Biochemistry (Mosc.) 2020, 85, 99–107. [Google Scholar] [CrossRef]
  52. Ruf, A.; de Murcia, G.; Schulz, G.E. Inhibitor and NAD+ binding to poly(ADP-ribose) polymerase as derived from crystal structures and homology modeling. Biochemistry 1998, 37, 3893–3900. [Google Scholar] [CrossRef]
  53. Gushchina, I.V.; Polenova, A.M.; Suplatov, D.A.; Švedas, V.K.; Nilov, D.K. vsFilt: A tool to improve virtual screening by structural filtration of docking poses. J. Chem. Inf. Model. 2020, 60, 3692–3696. [Google Scholar] [CrossRef]
  54. Papeo, G.; Posteri, H.; Borghi, D.; Busel, A.A.; Caprera, F.; Casale, E.; Ciomei, M.; Cirla, A.; Corti, E.; D’Anello, M.; et al. Discovery of 2-[1-(4,4-Difluorocyclohexyl)piperidin-4-yl]-6-fluoro-3-oxo-2,3-dihydro-1H-isoindole-4-carboxamide (NMS-P118): A potent, orally available, and highly selective PARP-1 inhibitor for cancer therapy. J. Med. Chem. 2015, 58, 6875–6898. [Google Scholar]
  55. Lindgren, A.E.; Karlberg, T.; Thorsell, A.G.; Hesse, M.; Spjut, S.; Ekblad, T.; Andersson, C.D.; Pinto, A.F.; Weigelt, J.; Hottiger, M.O.; et al. PARP inhibitor with selectivity toward ADP-ribosyltransferase ARTD3/PARP3. ACS Chem. Biol. 2013, 8, 1698–1703. [Google Scholar] [CrossRef] [PubMed]
  56. Johannes, J.W.; Almeida, L.; Barlaam, B.; Boriack-Sjodin, P.A.; Casella, R.; Croft, R.A.; Dishington, A.P.; Gingipalli, L.; Gu, C.; Hawkins, J.L.; et al. Pyrimidinone nicotinamide mimetics as selective tankyrase and wnt pathway inhibitors suitable for in vivo pharmacology. ACS Med. Chem. Lett. 2015, 6, 254–259. [Google Scholar] [CrossRef] [Green Version]
  57. Hua, Z.; Bregman, H.; Buchanan, J.L.; Chakka, N.; Guzman-Perez, A.; Gunaydin, H.; Huang, X.; Gu, Y.; Berry, V.; Liu, J.; et al. Development of novel dual binders as potent, selective, and orally bioavailable tankyrase inhibitors. J. Med. Chem. 2013, 56, 10003–10015. [Google Scholar] [CrossRef] [PubMed]
  58. Nkizinkiko, Y.; Desantis, J.; Koivunen, J.; Haikarainen, T.; Murthy, S.; Sancineto, L.; Massari, S.; Ianni, F.; Obaji, E.; Loza, M.I.; et al. 2-Phenylquinazolinones as dual-activity tankyrase-kinase inhibitors. Sci. Rep. 2018, 8, 1680. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Thorsell, A.G.; Ekblad, T.; Karlberg, T.; Löw, M.; Pinto, A.F.; Trésaugues, L.; Moche, M.; Cohen, M.S.; Schüler, H. Structural basis for potency and promiscuity in poly(ADP-ribose) polymerase (PARP) and tankyrase inhibitors. J. Med. Chem. 2017, 60, 1262–1271. [Google Scholar] [CrossRef] [PubMed]
  60. Karlberg, T.; Klepsch, M.; Thorsell, A.G.; Andersson, C.D.; Linusson, A.; Schüler, H. Structural basis for lack of ADP-ribosyltransferase activity in poly(ADP-ribose) polymerase-13/zinc finger antiviral protein. J. Biol. Chem. 2015, 290, 7336–7344. [Google Scholar] [CrossRef] [Green Version]
  61. Andersson, C.D.; Karlberg, T.; Ekblad, T.; Lindgren, A.E.; Thorsell, A.G.; Spjut, S.; Uciechowska, U.; Niemiec, M.S.; Wittung-Stafshede, P.; Weigelt, J.; et al. Discovery of ligands for ADP-ribosyltransferases via docking-based virtual screening. J. Med. Chem. 2012, 55, 7706–7718. [Google Scholar] [CrossRef]
  62. Karlberg, T.; Thorsell, A.G.; Kallas, Å.; Schüler, H. Crystal structure of human ADP-ribose transferase ARTD15/PARP16 reveals a novel putative regulatory domain. J. Biol. Chem. 2012, 287, 24077–24081. [Google Scholar] [CrossRef] [Green Version]
  63. Haikarainen, T.; Krauss, S.; Lehtio, L. Tankyrases: Structure, function and therapeutic implications in cancer. Curr. Pharm. Des. 2014, 20, 6472–6488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Nilov, D.K.; Tararov, V.I.; Kulikov, A.V.; Zakharenko, A.L.; Gushchina, I.V.; Mikhailov, S.N.; Lavrik, O.I.; Švedas, V.K. Inhibition of poly(ADP-ribose) polymerase by nucleic acid metabolite 7-methylguanine. Acta Naturae 2016, 8, 108–115. [Google Scholar] [CrossRef] [PubMed]
  65. Nilov, D.; Maluchenko, N.; Kurgina, T.; Pushkarev, S.; Lys, A.; Kutuzov, M.; Gerasimova, N.; Feofanov, A.; Švedas, V.; Lavrik, O.; et al. Molecular mechanisms of PARP-1 inhibitor 7-methylguanine. Int. J. Mol. Sci. 2020, 21, 2159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Nilov, D.; Kirsanov, K.; Antoshina, E.; Maluchenko, N.; Feofanov, A.; Kurgina, T.; Zakharenko, A.; Khodyreva, S.; Gerasimova, N.; Studitsky, V.; et al. 7-Methylguanine: A natural DNA repair inhibitor and a promising anticancer compound. FEBS Open Bio 2018, 8. [Google Scholar] [CrossRef]
  67. Maluchenko, N.; Nilov, D.; Feofanov, A.; Lys, A.; Kutuzov, M.; Gerasimova, N.; Studitsky, V. 7-Methylguanine traps PARP-1 on nucleosomes: spFRET microscopy study. Microsc. Microanal. 2019, 25. [Google Scholar] [CrossRef] [Green Version]
  68. Langelier, M.F.; Eisemann, T.; Riccio, A.A.; Pascal, J.M. PARP family enzymes: Regulation and catalysis of the poly(ADP-ribose) posttranslational modification. Curr. Opin. Struct. Biol. 2018, 53, 187–198. [Google Scholar] [CrossRef]
  69. Kuimov, A.; Nilov, D. Personal Communication; Belozersky Institute of Physicochemical Biology, Lomonosov Moscow State University: Moscow, Russia, 2021. [Google Scholar]
  70. Larsson, E.A.; Jansson, A.; Ng, F.M.; Then, S.W.; Panicker, R.; Liu, B.; Sangthongpitag, K.; Pendharkar, V.; Tai, S.J.; Hill, J.; et al. Fragment-based ligand design of novel potent inhibitors of tankyrases. J. Med. Chem. 2013, 56, 4497–4508. [Google Scholar] [CrossRef] [PubMed]
  71. Baker, D.; Sali, A. Protein structure prediction and structural genomics. Science 2001, 294, 93–96. [Google Scholar] [CrossRef] [Green Version]
  72. Frampton, J.E. Olaparib: A review of its use as maintenance therapy in patients with ovarian cancer. BioDrugs 2015, 29, 143–150. [Google Scholar] [CrossRef]
  73. Keung, M.Y.T.; Wu, Y.; Vadgama, J.V. PARP inhibitors as a therapeutic agent for homologous recombination deficiency in breast cancers. J. Clin. Med. 2019, 8, 435. [Google Scholar] [CrossRef] [Green Version]
  74. Riffell, J.L.; Lord, C.J.; Ashworth, A. Tankyrase-targeted therapeutics: Expanding opportunities in the PARP family. Nat. Rev. Drug Discov. 2012, 11, 923–936. [Google Scholar] [CrossRef]
  75. Qiu, W.; Lam, R.; Voytyuk, O.; Romanov, V.; Gordon, R.; Gebremeskel, S.; Vodsedalek, J.; Thompson, C.; Beletskaya, I.; Battaile, K.P.; et al. Insights into the binding of PARP inhibitors to the catalytic domain of human tankyrase-2. Acta Crystallogr. D Biol. Crystallogr. 2014, 70, 2740–2753. [Google Scholar] [CrossRef]
  76. Thorvaldsen, T.E. Targeting tankyrase to fight WNT-dependent tumours. Basic Clin. Pharmacol. Toxicol. 2017, 121, 81–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Berishvili, V.P.; Kuimov, A.N.; Voronkov, A.E.; Radchenko, E.V.; Kumar, P.; Choonara, Y.E.; Pillay, V.; Kamal, A.; Palyulin, V.A. Discovery of novel tankyrase inhibitors through molecular docking-based virtual screening and molecular dynamics simulation studies. Molecules 2020, 25, 3171. [Google Scholar] [CrossRef] [PubMed]
  78. Kam, C.M.; Tauber, A.L.; Oosthuizen, D.L.; Levonis, S.M.; Schweiker, S.S. In silico family-wide profiling and 3D modelling of the poly(ADP-ribose) polymerase superfamily. Future Med. Chem. 2020, 12, 2105–2122. [Google Scholar] [CrossRef]
  79. Koushyar, S.; Powell, A.G.; Vincan, E.; Phesse, T.J. Targeting Wnt signaling for the treatment of gastric cancer. Int. J. Mol. Sci. 2020, 21, 3927. [Google Scholar] [CrossRef] [PubMed]
  80. Waaler, J.; Leenders, R.G.G.; Sowa, S.T.; Alam Brinch, S.; Lycke, M.; Nieczypor, P.; Aertssen, S.; Murthy, S.; Galera-Prat, A.; Damen, E.; et al. Preclinical lead optimization of a 1,2,4-triazole based tankyrase inhibitor. J. Med. Chem. 2020, 63, 6834–6846. [Google Scholar] [CrossRef]
  81. UniProt Consortium. UniProt: A worldwide hub of protein knowledge. Nucleic Acids Res. 2019, 47, D506–D515. [Google Scholar] [CrossRef] [Green Version]
  82. Sievers, F.; Higgins, D.G. Clustal Omega, accurate alignment of very large numbers of sequences. Methods Mol. Biol. 2014, 1079, 105–116. [Google Scholar]
  83. Madeira, F.; Park, Y.M.; Lee, J.; Buso, N.; Gur, T.; Madhusoodanan, N.; Basutkar, P.; Tivey, A.R.N.; Potter, S.C.; Finn, R.D.; et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019, 47, W636–W641. [Google Scholar] [CrossRef] [Green Version]
  84. Sievers, F.; Higgins, D.G. The Clustal Omega multiple alignment package. In Methods in Molecular Biology; Katoh, K., Ed.; Humana: New York, NY, USA, 2021; Volume 2231, pp. 3–16. [Google Scholar]
  85. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [Green Version]
  86. Menke, M.; Berger, B.; Cowen, L. Matt: Local flexibility aids protein multiple structure alignment. PLoS Comput. Biol. 2008, 4, e10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Sali, A.; Blundell, T.L. Comparative protein modelling by satisfaction of spatial restraints. J. Mol. Biol. 1993, 234, 779–815. [Google Scholar] [CrossRef]
  88. Webb, B.; Sali, A. Comparative protein structure modeling using MODELLER. Curr. Protoc. Bioinform. 2016, 54, 1–37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Case, D.A.; Berryman, J.T.; Betz, R.M.; Cerutti, D.S.; Cheatham, T.E., 3rd; Darden, T.A.; Duke, R.E.; Giese, T.J.; Gohlke, H.; Goetz, A.W.; et al. AMBER 2015; University of California: San Francisco, CA, USA, 2015; pp. 1–883. [Google Scholar]
  90. Salomon-Ferrer, R.; Case, D.A.; Walker, R.C. An overview of the Amber biomolecular simulation package. WIREs Comput. Mol. Sci. 2013, 3, 198–210. [Google Scholar] [CrossRef]
  91. Stroganov, O.V.; Novikov, F.N.; Stroylov, V.S.; Kulkov, V.; Chilov, G.G. Lead finder: An approach to improve accuracy of protein-ligand docking, binding energy estimation, and virtual screening. J. Chem. Inf. Model. 2008, 48, 2371–2385. [Google Scholar] [CrossRef]
  92. Novikov, F.N.; Stroylov, V.S.; Zeifman, A.A.; Stroganov, O.V.; Kulkov, V.; Chilov, G.G. Lead Finder docking and virtual screening evaluation with Astex and DUD test sets. J. Comput. Aided Mol. Des. 2012, 26, 725–735. [Google Scholar] [CrossRef]
  93. Zimmermann, L.; Stephens, A.; Nam, S.Z.; Rau, D.; Kübler, J.; Lozajic, M.; Gabler, F.; Söding, J.; Lupas, A.N.; Alva, V. A completely reimplemented MPI Bioinformatics Toolkit with a new HHpred server at its core. J. Mol. Biol. 2018, 430, 2237–2243. [Google Scholar] [CrossRef]
  94. Gabler, F.; Nam, S.Z.; Till, S.; Mirdita, M.; Steinegger, M.; Söding, J.; Lupas, A.N.; Alva, V. Protein sequence analysis using the MPI Bioinformatics Toolkit. Curr. Protoc. Bioinform. 2020, 72, e108. [Google Scholar] [CrossRef]
  95. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [Green Version]
  96. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
Cancers 13 01201 g001 550
Figure 1. Multiple sequence alignment of PARPs 1–16 with Clustal Omega. The residues 863, 898, 904, and 907 of the NA binding site (PARP-1 numbering) are marked with an asterisk. PARP family members with known structures are marked in green, and PARPs with unknown structures in red.
Figure 1. Multiple sequence alignment of PARPs 1–16 with Clustal Omega. The residues 863, 898, 904, and 907 of the NA binding site (PARP-1 numbering) are marked with an asterisk. PARP family members with known structures are marked in green, and PARPs with unknown structures in red.
Cancers 13 01201 g001
Cancers 13 01201 g002 550
Figure 2. Structure of the NA binding site in PARPs 1–16. Underlined PARP structures (1–3, 5a, 5b, 10, 12–16) were obtained from the Protein Data Bank, and other structures (4, 6–9, 11) were obtained by homology modeling. Amino acid substitutions in the NA site are marked in red. The NA site architecture is quite similar in PARPs 1–3, 5a, 5b, 10, 12, 14, and 15 due to the lack of substitutions. The coordinates of the 7-MG molecule (colored by atom type) are transferred from the docking model of PARP-1–7-MG complex.
Figure 2. Structure of the NA binding site in PARPs 1–16. Underlined PARP structures (1–3, 5a, 5b, 10, 12–16) were obtained from the Protein Data Bank, and other structures (4, 6–9, 11) were obtained by homology modeling. Amino acid substitutions in the NA site are marked in red. The NA site architecture is quite similar in PARPs 1–3, 5a, 5b, 10, 12, 14, and 15 due to the lack of substitutions. The coordinates of the 7-MG molecule (colored by atom type) are transferred from the docking model of PARP-1–7-MG complex.
Cancers 13 01201 g002
Cancers 13 01201 g003 550
Figure 3. 7-MG position in the NA binding site of PARPs 1–3 obtained using MD modeling. The NA site residues are colored in green.
Figure 3. 7-MG position in the NA binding site of PARPs 1–3 obtained using MD modeling. The NA site residues are colored in green.
Cancers 13 01201 g003
Cancers 13 01201 g004 550
Figure 4. 7-MG position in the NA binding site of PARP-5a (conformations I and II) and 5b obtained using MD modeling. The NA site residues are colored in green.
Figure 4. 7-MG position in the NA binding site of PARP-5a (conformations I and II) and 5b obtained using MD modeling. The NA site residues are colored in green.
Cancers 13 01201 g004
Cancers 13 01201 g005 550
Figure 5. D-loop conformation in PARP-5a and its interaction with the NA binding site revealed by MD modeling. The NA site residues are colored in green, D-loop conformation I is shown in blue, conformation II in yellow.
Figure 5. D-loop conformation in PARP-5a and its interaction with the NA binding site revealed by MD modeling. The NA site residues are colored in green, D-loop conformation I is shown in blue, conformation II in yellow.
Cancers 13 01201 g005
Cancers 13 01201 g006 550
Figure 6. 7-MG position in the NA binding site of PARP-10 and 12 obtained using MD modeling. The NA site residues are colored in green.
Figure 6. 7-MG position in the NA binding site of PARP-10 and 12 obtained using MD modeling. The NA site residues are colored in green.
Cancers 13 01201 g006
Cancers 13 01201 g007 550
Figure 7. 7-MG position in the NA binding site of PARP-13 obtained using MD modeling. The conserved NA site residues are colored in green, and unique residues in red.
Figure 7. 7-MG position in the NA binding site of PARP-13 obtained using MD modeling. The conserved NA site residues are colored in green, and unique residues in red.
Cancers 13 01201 g007
Cancers 13 01201 g008 550
Figure 8. 7-MG position in the NA binding site of PARP-14 and 15 obtained using MD modeling. The NA site residues are colored in green.
Figure 8. 7-MG position in the NA binding site of PARP-14 and 15 obtained using MD modeling. The NA site residues are colored in green.
Cancers 13 01201 g008
Cancers 13 01201 g009 550
Figure 9. D-loop conformation in PARP-15 and its interaction with 7-MG. The NA binding site is colored in green, and D-loop in blue.
Figure 9. D-loop conformation in PARP-15 and its interaction with 7-MG. The NA binding site is colored in green, and D-loop in blue.
Cancers 13 01201 g009
Cancers 13 01201 g010 550
Figure 10. 7-MG position in the NA binding site of PARP-16 obtained using MD modeling. The conserved NA site residues are colored in green, and unique residues in red. The Cγ atom of Thr184 forms a hydrophobic contact with the 7-MG methyl group.
Figure 10. 7-MG position in the NA binding site of PARP-16 obtained using MD modeling. The conserved NA site residues are colored in green, and unique residues in red. The Cγ atom of Thr184 forms a hydrophobic contact with the 7-MG methyl group.
Cancers 13 01201 g010
Cancers 13 01201 g011 550
Figure 11. 7-MG position in the NA binding site of PARPs 4, 6, and 8 obtained using homology modeling. The conserved NA site residues are colored in green and unique residues in red. The Cβ atom of Ser479/510/734 forms a hydrophobic contact with the 7-MG methyl group.
Figure 11. 7-MG position in the NA binding site of PARPs 4, 6, and 8 obtained using homology modeling. The conserved NA site residues are colored in green and unique residues in red. The Cβ atom of Ser479/510/734 forms a hydrophobic contact with the 7-MG methyl group.
Cancers 13 01201 g011
Cancers 13 01201 g012 550
Figure 12. 7-MG position in the NA binding site of PARP-7 and 11 obtained using homology modeling. The conserved NA site residues are colored in green and unique residues in red.
Figure 12. 7-MG position in the NA binding site of PARP-7 and 11 obtained using homology modeling. The conserved NA site residues are colored in green and unique residues in red.
Cancers 13 01201 g012
Cancers 13 01201 g013 550
Figure 13. 7-MG position in the NA binding site of PARP-9 obtained using homology modeling. The Gln706 main chain forms hydrogen bonds with the lactam group. The Cγ atom of Thr739 forms a hydrophobic contact with the 7-MG methyl group.
Figure 13. 7-MG position in the NA binding site of PARP-9 obtained using homology modeling. The Gln706 main chain forms hydrogen bonds with the lactam group. The Cγ atom of Thr739 forms a hydrophobic contact with the 7-MG methyl group.
Cancers 13 01201 g013
Table
Table 1. Catalytic activity and cellular functions of PARP family proteins.
Table 1. Catalytic activity and cellular functions of PARP family proteins.
PARPActivityFunction
1poly 1DNA repair enzymes regulating transcription and eliminating single- and double-strand breaks of DNA [13,14,15,16,17]
2poly
3mono 2DNA repair [18,19]
4monopoorly studied
5apolytankyrase enzymes regulating the Wnt signaling pathway and controlling cellular proliferation and differentiation [20,21,22,23]
5bpoly
6monotumor suppressor regulating Survivin expression [24]
7monopoorly studied
8monopoorly studied
9inactivepoorly studied
10monoregulates Aurora A and suppresses tumor metastasis [25]
11monopoorly studied
12monosuppresses Zika virus infection [26]
13inactiveregulates the cellular response to stress [27]
14monosuppresses kinase-mediated apoptosis [28]
15monopoorly studied
16monopoorly studied
1 Poly(ADP-ribosyl)ation. 2 Mono(ADP-ribosyl)ation.
Table
Table 2. Residues of the NA binding site in PARPs 1–16. Amino acid replacements are marked in gray.
Table 2. Residues of the NA binding site in PARPs 1–16. Amino acid replacements are marked in gray.
PARPResidues of the NA SitePDB ID 1
1Gly863Ala898Ser904Tyr9074zzz [54]
2Gly429Ala464Ser470Tyr4734zzx [54]
3Gly385Ala416Ser422Tyr4254gv2 [55]
4Gly439Ser479Ser485Tyr488-
5aGly1185Ala1215Ser1221Tyr12244w6e/4msg [56,57]
5bGly1032Ala1062Ser1068Tyr10715nwg [58]
6Gly474Ser510Ser516Tyr519-
7Gly533Ala566Ser572Phe575-
8Gly698Ser734Ser740Tyr743-
9Gln706Thr739Leu745Lys750-
10Gly888Ala921Ser927Tyr9325lx6 [59]
11Gly205Ala238Ser244Phe247-
12Gly565Ala598Ser604Tyr6072pqf [60]
13Ala788Ala821Ser827Asn8302x5y [60]
14Gly1602Ala1635Ser1641Tyr16463smj [37]
15Gly538Ala571Ser577Tyr5824f0e [61]
16Gly153Thr184Ala190Tyr1934f0d [62]
1 Representative crystal structures.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Manasaryan, G.; Suplatov, D.; Pushkarev, S.; Drobot, V.; Kuimov, A.; Švedas, V.; Nilov, D. Bioinformatic Analysis of the Nicotinamide Binding Site in Poly(ADP-Ribose) Polymerase Family Proteins. Cancers 2021, 13, 1201. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13061201

AMA Style

Manasaryan G, Suplatov D, Pushkarev S, Drobot V, Kuimov A, Švedas V, Nilov D. Bioinformatic Analysis of the Nicotinamide Binding Site in Poly(ADP-Ribose) Polymerase Family Proteins. Cancers. 2021; 13(6):1201. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13061201

Chicago/Turabian Style

Manasaryan, Garri, Dmitry Suplatov, Sergey Pushkarev, Viktor Drobot, Alexander Kuimov, Vytas Švedas, and Dmitry Nilov. 2021. "Bioinformatic Analysis of the Nicotinamide Binding Site in Poly(ADP-Ribose) Polymerase Family Proteins" Cancers 13, no. 6: 1201. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13061201

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop