Next Article in Journal
A New Laboratory Workflow Integrating the Free Light Chains Kappa Quotient into Routine CSF Analysis
Next Article in Special Issue
Proteome-Wide Detection and Annotation of Receptor Tyrosine Kinases (RTKs): RTK-PRED and the TyReK Database
Previous Article in Journal
Small Non-Coding RNAs as New Biomarkers to Evaluate the Quality of the Embryo in the IVF Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Evolutionary Perspective on the Origin, Conservation and Binding Partner Acquisition of Tankyrases

1
Faculty for Biochemistry and Molecular Medicine & Biocenter Oulu, University of Oulu, 90220 Oulu, Finland
2
Structural Bioinformatics Laboratory, Biochemistry, Faculty of Science and Engineering and InFLAMES Research Flagship Center, Åbo Akademi University, 20520 Turku, Finland
*
Author to whom correspondence should be addressed.
Biomolecules 2022, 12(11), 1688; https://doi.org/10.3390/biom12111688
Submission received: 29 August 2022 / Revised: 8 November 2022 / Accepted: 8 November 2022 / Published: 15 November 2022

Abstract

:
Tankyrases are poly-ADP-ribosyltransferases that regulate many crucial and diverse cellular processes in humans such as Wnt signaling, telomere homeostasis, mitotic spindle formation and glucose metabolism. While tankyrases are present in most animals, functional differences across species may exist. In this work, we confirm the widespread distribution of tankyrases throughout the branches of multicellular animal life and identify the single-celled choanoflagellates as earliest origin of tankyrases. We further show that the sequences and structural aspects of TNKSs are well-conserved even between distantly related species. We also experimentally characterized an anciently diverged tankyrase homolog from the sponge Amphimedon queenslandica and show that the basic functional aspects, such as poly-ADP-ribosylation activity and interaction with the canonical tankyrase binding peptide motif, are conserved. Conversely, the presence of tankyrase binding motifs in orthologs of confirmed interaction partners varies greatly between species, indicating that tankyrases may have different sets of interaction partners depending on the animal lineage. Overall, our analysis suggests a remarkable degree of conservation for tankyrases, and that their regulatory functions in cells have likely changed considerably throughout evolution.

1. Introduction

Tankyrases (TNKSs) are part of the ARTD enzyme family and they catalyze a sequential transfer of ADP-ribose from NAD+ to their protein substrates leaving them modified with poly-ADP-ribosyl (PAR) chains [1,2]. TNKSs regulate a wide variety of pathways and physiological processes, such as telomere maintenance [3,4], glucose metabolism [5,6], bone growth [7], mitotic spindle formation [8,9], Wnt/β-catenin signaling [10] and Hippo/YAP-signaling [11]. Humans and other vertebrates possess two highly similar TNKS paralogs: TNKS1 and TNKS2. It is yet unclear if these paralogs have different functions in the cell; however, they are at least redundant to some degree as indicated by knockout experiments in mice [12]. Among ARTD members, which are also referred to as PARPs [13], TNKSs have a unique domain architecture (Figure 1a). In addition to their catalytic ADP-ribosyltransferase (ART) domain, TNKSs have five N-terminal ankyrin-repeat (ANK) domains, which are also termed ankyrin repeat clusters (ARCs) and they are numbered ARC1 to ARC5. With the exception of ARC3, which may have a structural role [14,15,16], the ARCs bind interaction partners of TNKSs. These interaction partners bind via a small linear TNKS binding motif (TBM). TBMs have a canonical sequence Rxx[ACGP]xGxx, although also less constrained variations such as Rx(4)Gxx and Rx(5)Gxx were reported to bind to TNKSs [10,15,16,17]. A sterile alpha motif (SAM) domain is located between the ARC and ART domains and mediates the multimerization of TNKSs [18,19,20,21]. In contrast to TNKS2, TNKS1 has an additional and likely disordered N-terminal extension of unknown function which is rich in histidine, proline and serine residues and was accordingly termed an HPS region [22].
In the current model, PARylation of proteins binding to TNKSs leads to their ubiquitination and subsequent proteasomal degradation (Figure 1b). Additionally, TNKSs may act as scaffolding proteins by providing a platform for protein–protein interactions [19,21,26]. Importantly, the regulatory functions that TNKSs exert on different pathways in the cell are a consequence of specific binding partners interacting with TNKSs. Numerous proteins interacting with TNKSs have been reported in humans, although it has become clear that these interactions are conserved in varying degrees between species (Figure 1c). For example, TNKSs regulate telomere length homeostasis through the interaction with the telomeric protein TRF1 [4,27]. This interaction and thus the function of TNKSs in telomere homeostasis is absent in mice due to a missing TBM in the mouse TRF1 ortholog [23,28]. In contrast, the function of TNKSs as positive regulators of Wnt signaling through interaction and regulation of Axin levels appears well-conserved across species, and presence of the TBM sequence in Axin orthologs was reported even in invertebrates [10]. The regulatory function of TNKSs on Wnt signaling was shown in several species such as human [10], mouse [29], zebrafish [30] and fruit fly [31].
While many interaction partners of TNKSs have been confirmed in humans [1,32,33,34], the presence of TBMs in orthologs of other species was only shown in some cases and to varying degrees in terms of taxonomic coverage (Figure 1c). In this work, we examined the distribution of TNKSs across species, the conservation in terms of sequence and structure as well as the conservation of the presence of TBMs in orthologs of 21 human binding partners. We show that TNKSs are widely distributed throughout multicellular animals and we further identify single-celled choanoflagellates as the likely earliest origin of TNKSs. We demonstrate further that TNKSs display a high degree of conservation and show through in vitro experiments with an early diverged TNKS ortholog from the demosponge Amphimedon queenslandica that the basic molecular functions of TNKSs are also highly conserved. On the other hand, the interaction with different TNKS binding partners appears to be far less conserved, as the presence of TBMs in orthologs varies greatly among different metazoans.

2. Materials and Methods

2.1. Sequence Acquisition

Canonical sequences from human orthologs were obtained from UniProt. Sequences from non-human orthologs were obtained from NCBI non-redundant (NR) database using sequences of human orthologs as a query by BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi (accessed on 31 October 2022)) [35]. If no ortholog could be identified in protein databases, the tBLASTn implementation against the Genebank database (TSA and WGS) was used. Sequences of TNKSs from E. burgeri and M. leidyi were acquired from the Ensembl database (https://metazoa.ensembl.org/ (accessed on 20 April 2022)). Possible functional annotations for the TNKS-like protein from N. fowleri were examined in the Amoeba database (https://amoebadb.org/ (accessed on 20 April 2022)) under the gene identifier NF0130940.

2.2. Sequence Alignments and Construction of Phylogenetic Trees

Multiple sequence alignments were conducted with Clustal Omega [36] or the Clustal W implementation in the MEGA11 software [37]. Global pairwise sequence alignments were conducted with EMBOSS Needle and local pairwise alignments were conducted with EMBOSS Matcher. The construction of phylogenetic trees was conducted with the maximum likelihood and neighbor joining methods in MEGA11 [37] using the Jones–Taylor–Thornton (JTT) substitution model [38] and a bootstrap analysis with 1000 replicates was performed for each tree.

2.3. Structure Visualization and Prediction

Models of experimentally derived protein structures were obtained from the Protein Data Bank (PDB) [39]. The structure models were visualized in the PyMOL Molecular Graphics System (PyMOL, version 1.8.4.0). The prediction of protein structures from A. manzaensis and N. fowleri were conducted with RoseTTAFold [40] using only the protein sequences as input without additional constraints.

2.4. Position-Specific Sequence Conservation Analysis of Tankyrases

Sequences of TNKSs originating from a broad taxonomic range were sampled and obtained as described above. Sequences annotated as low-quality proteins were removed and the remaining sequences were inspected for obvious indications of assembly errors. A multiple sequence alignment from 113 TNKS sequences was used as input for residue specific conservation analysis by ConSurf [41]. JTT was determined to be the best evolutionary model for the data provided and was used by ConSurf to estimate the evolutionary relationships. ConSurf conservation scores were calculated based on the Bayesian method.

2.5. Identification of Potential TBMs in Orthologs of TNKS Binders

A list of 21 human TNKS binding partners was manually curated based on reports in the literature. Sequences for the orthologs of the human TNKS interaction partners were obtained as described above. A multiple sequence alignment was performed for each set of orthologs, and the presence of potential TBMs in respect to the human ortholog was examined for the proteins in each set. The different configurations Rxx[ACGP]xGxx, Rx(4)Gxx and Rx(5)Gxx were accepted as possible TBMs. Generally, only those TBMs that aligned with the position of the TBM in the human ortholog were considered potential TBMs; although not aligned, some possible TBMs were annotated due to presence in a likely disordered region in proximity to the TBM in human orthologs. For RNF146, all possible TBMs were located in the protein sequences of the orthologs. The positions of domains in the RNF146 orthologs were predicted using InterPro [42].

2.6. Molecular Cloning

The cloning of human TNKS1 constructs was conducted as previously described for CFP-ARC5 [43] or SAM-ART and ART domains [44]. Cloning of the YFP-REAGDGEE construct was conducted as previously described [43]. The E. coli codon-optimized constructs for A. queenslandica TNKS (AqTNKS: XP_019848937.1) were obtained as gBlock from IDT. The construct ARC5 (AqTNKS Glu658-Met811) was cloned into the pNIC-CFP expression vector (Addgene #173074). The AqTNKS constructs SAM-ART (Pro889-Thr1181) and ART (Thr968-Thr1181) were cloned into the pNIC-MBP vector, which was prepared as previously described [43]. Cloning was conducted using the SLIC method [45]. Briefly, the vectors were linearized and mixed with insert and T4 DNA polymerase. The reaction mixture was used to transform chemicompetent E. coli NEB 5α cells (New England BioLabs). Colonies were grown at 37 °C on LB agar plates containing 5% sucrose and 50 μg/mL kanamycin; genes in the vector encoding kanamycin-resistance and SacB enzyme [46,47] served as selection markers for successful transformation and vector linearization, respectively. All constructs were verified by sequencing of the insert regions.

2.7. Protein Expression

Chemicompetent E. coli BL21(DE3) cells were transformed with the expression vectors. For each construct, 500 mL Terrific Broth (TB) autoinduction media including trace elements (Formedium, Hunstanton, Norfolk, England) were supplemented with 8 g/l glycerol and 50 μg/mL kanamycin and inoculated with overnight preculture. The flasks were incubated at 37 °C (shaking) until an OD600 of 1 was reached. The incubation continued overnight at 18 °C (shaking). At the end of the incubation period, the cells were harvested by centrifugation at 5000× g, 4 °C for 30 min. The pellets were resuspended in lysis buffer (50 mM HEPES pH 7.5, 500 mM NaCl, 0.5 mM TCEP, 10 mM imidazole, 10% glycerol v/v) and used immediately for purification or frozen at −20 °C.

2.8. Protein Purification

Purification of human TNKS1 ART and TNKS1 SAM-ART [44] or CFP-ARC5 (TNKS1) and YFP-REAGDGEE [43] was conducted as previously described. Frozen cells were thawed first. The cells were supplemented with 0.1 mM Pefabloc SC (Roche) and 20 μg/mL DNase I (Roche) and lysed by sonication. The lysate was centrifuged (30,000× g, 4 °C, 30 min) and filtered (0.45 μm). As initial purification step, all the proteins were purified by immobilized metal affinity chromatography (IMAC) with a 5 mL HiTrap HP column equilibrated with lysis buffer and charged with Ni2+. Lysis buffer with higher concentrations of imidazole was used for washing steps (25 mM) and elution (300 mM).
The IMAC purification of the AqTNKS ART domain (MBP-tagged) was followed by buffer exchange with size-exclusion chromatography (SEC) buffer (20 mM HEPES pH 7.5, 350 mM NaCl, 0.5 mM TCEP, 10% v/v glycerol) using an Amicon Ultra-15 centrifugal filter (30 kDa MWCO). To remove MBP, the sample was incubated with TEV protease (1:30 molar ratio) for 48 h at 4 °C and was then run over a 5 mL MBPTrap HP column equilibrated with SEC buffer. SEC was performed as last step of the purification using an S75 16/600 size-exclusion chromatography column and SEC buffer.
For AqTNKS SAM-ART (MBP-tagged), the IMAC step was followed by MBP affinity chromatography using a 5 mL MBPTrap HP column equilibrated with SEC buffer. After loading, the column was washed with 5–10 column volumes of SEC buffer and eluted with SEC buffer containing 10 mM maltose.
For the AqTNKS ARC5 domain (CFP-tagged), the elution from IMAC was purified by SEC as above. To remove contamination by E. coli chaperone proteins, a successive IMAC purification was performed and 5 mM MgCl2 and 5 mM ATP disodium salt were included in the washing step [48]. Finally, a second purification by SEC was performed as above. At the end of each purification, all the proteins underwent concentration, they were aliquoted, flash frozen in liquid nitrogen and stored at −70 °C.

2.9. Activity Analysis by Western Blot

The constructs of TNKS1 or AqTNKS were prepared in absence or presence of a mixture of biotinylated-NAD+/NAD+ (1:10). Additionally, a control with PARG (0.2 mg/mL) was prepared. Final total NAD+-concentrations were 2 μM and 10 μM for TNKS1 and AqTNKS, respectively. The reactions were incubated for 1 h at room temperature. At the end of the incubation period, 10 μL of each reaction were run on SDS-PAGE (Mini-Protean TGX 4–20% gradient gel, Bio-Rad). The transfer of the proteins to a nitrocellulose membrane (Trans-Blot Turbo Mini 0.2 μm Nitrocellulose Transfer Pack, Bio-Rad) was performed using a Trans-Blot semi-dry system (Bio-Rad). Afterwards, the membrane was washed with TBS-T and subsequently stained with Ponceau S solution. The membranes were imaged (ChemiDoc Imaging System, Bio-Rad) and the staining was removed by washing the membrane in TBS-T. The membranes were blocked using 1× TBS with 1% casein (Bio-Rad) for 30 min. The membranes were incubated with streptavidin-HRP (PerkinElmer) diluted 1:10,000 in blocking buffer for 30 min. The membranes were washed with TBS-T buffer and HRP detection was conducted using ECL solution (Advansta). Following the imaging step, the membranes were washed in TBS-T containing 5% (w/v) skimmed milk powder for 15 min. The membranes were transferred to blocking solution including 0.1 μg/mL nanoluciferase-ALC1 [49] for 30 min. Finally, the membrane was washed in TBS-T and detection of nanoluciferase on the membrane was conducted with 1:500 NanoGlo (Promega) diluted in PBS buffer.

2.10. Analytical Size-Exclusion Chromatography

In separate independent runs, 0.5 mg of each protein sample (AqTNKS SAM-ART, AqTNKS ART, TNKS1 SAM-CAT and MBP) was loaded to a Superdex 200 Increase 10/300 GL column. The column was equilibrated in SEC buffer (20 mM HEPES pH 7.5, 350 mM NaCl, 0.5 mM TCEP, 10% glycerol) and samples were run at a flow rate of 0.5 mL/min at 4 °C. The absorbance at 280 nm was recorded. In each run, 1 mL fractions were collected. Fractions corresponding to the different peaks observed for the SAM-ART constructs were analyzed by SDS-PAGE.

2.11. Differential Scanning Fluorimetry

Samples of ART constructs (5 μM) from TNKS1 or AqTNKS were prepared in absence or presence of TNKS inhibitors (25 μM) IWR-1, G007-LK, OM-153, compound 40, XAV939, E7449 or NVP-TNKS656 [50] and 5× SYPRO orange dye (ThermoFisher Scientific). Each condition was prepared in 4 replicates with a final volume of 20 μL. The sample buffer was 10 mM Bis-Tris-Propane (pH 7.0), 3% PEG 20,000 (w/v), 0.01% Triton X-100 (v/v), 0.5 mM TCEP. The samples were transferred into a 96-well transparent qPCR plate and the measurements were performed using a C1000 CFX96 thermal cycler (Bio-Rad). The temperature was increased from 20 °C to 95 °C (1 °C/min) and data points for melting curves were recorded in 1 min intervals. The data analysis was performed with GraphPad Prism 9 using a nonlinear regression analysis (Boltzmann sigmoid equation) of normalized data.

2.12. Determination of the Dissociation Constants by FRET

The experiment was performed as previously described [43]. Briefly, CFP-fused ARC5 domains (100 nM) from TNKS1 or AqTNKS were mixed with YFP-REAGDGEE (0–1500 nM). The FRET emission for each condition was determined. The calculation of dissociation constants was conducted using a method described by Song et al. [51].

3. Results

3.1. Distribution and Origin of Tankyrases

Species that are frequently mentioned in this manuscript are shown in Table 1. We started our analysis by characterizing the distribution of tankyrases throughout the tree of life (Figure 2a). The presence of TNKSs was searched within major phylogenetic groups. An early study on the evolutionary history of PARP proteins by Citarelli et al. only found canonical TNKSs present in bilaterians, while Perina et al. later identified TNKSs also in cnidarians and sponges, concluding that TNKSs may be confined to metazoans [52,53]. Indeed, we found and confirmed the presence of at least one TNKS homolog in all vertebrates analyzed. While TNKS homologs are present in early diverged bilaterians such as the fruit fly Drosophila melanogaster and the mollusk Octopus bimaculoides, no TNKS homolog was found in the nematode Caenorhabditis elegans. The absence of TNKSs and other PARP homologs in C. elegans was reported previously, and the loss of these PARP genes was discussed [52]. It is worth mentioning that TNKSs appear to be present in other nematodes such as Brugia malayi and Trichinella spiralis.
Furthermore, our search confirmed that TNKSs are present in early diverging metazoan species such as the sponge Amphimedon queenslandica, the cnidarians Hydra vulgaris and Nematostella vectensis and the ctenophore Mnemiopsis leidyi. However, we could not identify a TNKS homolog in the placozoan Trichoplax adhaerens. While the exact taxonomic placement of placozoans is still debated, it is widely accepted that they diverged after the appearance of sponges [54,55], which indicates a loss of TNKS genes in this lineage similarly to C. elegans. In fact, it was suggested that many physiological traits such as muscles, nerves and guts were lost in placozoans in favor of a simpler lifestyle [54,56], which likely resulted in the loss of many genes.
It has been suggested that the gene duplication event leading to TNKS1 and TNKS2 occurred at some point during the evolution of fishes [53]. Our analysis shows that the two TNKS paralogs can be identified in representatives of bony fishes (Osteichthyes) and cartilaginous fishes (Chondrichthyes), but not in the jawless fishes Petromyzon marinus and Eptatretus burgeri, indicating that the duplication event took place during the early evolution of vertebrates.
To interrogate a possible pre-metazoan origin, we searched for the presence of TNKSs in representative subsets of earlier diverging eukaryotes. No TNKSs were identified in the major eukaryotic groups Archaeplastida and Holomycota, which include plants and fungi, respectively. We then searched for the presence of TNKSs in groups more closely related to metazoans. While we did not find TNKSs in Filasterea or Ichthyosporea, we identified sequences coding for TNKSs in transcriptome shotgun assembly (TSA) data from the choanoflagellates Salpingoeca kvevrii and Salpingoeca macrocollata. Choanoflagellates are the proposed sister group of metazoans and are studied as models for pre-metazoan evolution along with their ability to live in communities and differentiate into different cell types [57,58]. It is worth mentioning that we did not detect the presence of TNKSs in the model choanoflagellates Monosiga brevicollis and Salpingoeca rosetta, for which the genomes have been sequenced [58,59]. However, these model choanoflagellates are relatively closely related and do not sufficiently represent the genetic diversity found among choanoflagellates [60]. A phylogenetic tree constructed from the TNKS protein sequences of several metazoan species as well as the choanoflagellates S. kvevrii and S. macrocollata reproduced the expected phylogenetic relationships (Figure 2b and Figure S1). These results suggest an origin of TNKSs in choanoflagellates.
Although several proteins such as Shank1, ANKS1b and ASZ1 are known to comprise combinations of ANK and SAM domains (Figure 2c), this domain architecture appears unique to TNKSs among PARPs. We speculated that gene-fusion or exon-shuffling [61] between an ancestral PARP and an ANK-SAM protein could represent obvious mechanisms by which TNKSs may have originated. To see if tangible TNKS precursor proteins could be identified, we searched with sub-domain sequences from human TNKS1 against non-metazoan sequence databases. Surprisingly, we identified proteins with ANK domains that share 40–50% sequence identity with the ARC5 domain of TNKS1 in diverse phylogenetic groups including prokaryotes. However, the high similarity may be explained by the fact that many residues in ANK-folds are well-conserved irrespective of the domain’s function [62]. Indeed, closer inspection of the hit proteins such as the ANK domain from the archaeon Acidianus manzaensis (Figure 2c) showed that binding of TBM peptides in a similar mode as observed for TNKS ARC domains would not be possible (Figure S2). From a similar search using the TNKS1 ART domain as the query, we identified a yet uncharacterized protein from the pathogenic amoeba Naegleria fowleri that shares 37% sequence identity with the ART domain of TNKS1 (Figure 2c). Further analysis indicates that the TNKS-like ART domain from N. fowleri (NfTNKSL) contains a putative zinc-binding motif corresponding to the motif in the ART domain of TNKSs (Figure S3). It was thought that this motif was unique to TNKSs among PARPs [63], and we showed recently that it is important for the structural integrity of the TNKS ART domain [44]. The available sequence information of NfTNKSL indicates that it does not contain other domains. Considering the similarity, NfTNKSL likely shares a common ancestor with TNKSs and it may possibly resemble an ancestral version of TNKSs before acquisition of ANK and SAM domains.

3.2. Conservation of Tankyrase Sequence and Structure

The results above confirm the vast distribution of TNKSs among metazoan groups and additionally in choanoflagellates. We next aimed to characterize how well-conserved TNKSs are throughout evolution in terms of sequence, structure and function. The following analyses were performed in respect to metazoans with well-characterized genomes; choanoflagellates were omitted from the analyses as their sequence information originated from TSA and may thus be incomplete and less reliable.
First, we compared the sequence identity of human TNKS1/2 and other human ARTD family members with their orthologs from representative metazoan species as derived from multiple sequence alignments (Figure 3). For further comparison, the ubiquitous and well-conserved proteins cytochrome C [64] and seryl-tRNA synthetase [65] were included as controls. It was reported that some PARPs such as PARP4 and PARP9 have experienced rapid evolution likely due to positive selection in response to host–virus conflicts [66]. Indeed, PARP4, PARP9 and PARP12 share less than 50% sequence identity between the homologs of humans and the fish S. salar. In contrast, TNKSs showed the highest degree of conservation of any analyzed family member, displaying approximately 60% sequence identity between human TNKS1/2 and the distant TNKS ortholog found in the demosponge A. queenslandica. The degree of conservation appears comparable to that of cytochrome C and seryl-tRNA-synthetase, suggesting that TNKSs are highly conserved members of the ARTD family. Remarkably, the level of conservation for TNKSs within the metazoan lineage appears greater than that for PARP1, which is known to act as key-player in DNA damage response [67] and its orthologs are widely distributed throughout eukaryotes in protists, animals, fungi and plants [53].
To gain insight into the localization of potentially highly or poorly conserved regions in TNKSs, we plotted the number of different residues for each position from an alignment of 113 TNKS sequences from taxonomically diverse metazoans (Figure 4a). We additionally performed analysis by ConSurf, which provides normalized grades for the conservation at each position [41]. Because the N-terminal region in TNKS1 is poorly conserved among species, we used human TNKS2 as reference protein to avoid ambiguity due to insufficient alignments in this region.
This analysis shows strong conservation for the ARC domains 1–2 and 4–5, the SAM domain and the ART domain. The ARC3 domain clearly shows a lower degree of conservation in comparison to the other ARCs. ARC3 was reported to not bind TBMs and likely has a structural role [14,15,69]. Overall, the lowest degree of conservation is found in the regions linking the domains ARC3 to ARC4, ARC5 to SAM and SAM to ART as well as the N- and C-terminal regions. These regions are likely intrinsically disordered, thus allowing a higher freedom in residue composition. A study by Eisemann et al. investigated the structural features of the ARC region from TNKS1 and demonstrated that the linking region from ARC3 to ARC4 is indeed flexible, while the regions from ARC1 to ARC2 and ARC2 to ARC3 form rigid interfaces [14]. Although points of flexibility were also observed between ARC4 and ARC5 [14], our analysis shows a high level of conservation at their interface, possibly indicating that the flexibility does not originate from an intrinsically disordered linker region.
In order to further characterize the sequence conservation of TNKSs in respect to their structures, we mapped strongly conserved residues to TNKS domain structures. From an alignment of the TNKS sequences from 11 metazoan species with well-characterized genomic information (Figure 4b), we mapped only identical residues to the domain structures of ARCs, SAM and ART (Figure 4c–e and Figure S4). The ARC domains 1, 2, 4 and 5 show a large, continuous patch of highly conserved residues around the TBM-binding pocket (Figure 4c and Figure S4), indicating that the ability to bind the same peptide motif is conserved throughout TNKS evolution. Such a conserved patch is missing for ARC3 (Figure S4), which is in agreement with its inability to bind TBM peptides in human TNKSs. Indeed, inspection of ARC3 sequences from early diverging metazoans confirms the absence of crucial residues required for binding TBMs (Figure S5), suggesting that the ARC3 domain did not only recently lose the TBM-binding function. Sequences from the SAM domains have fewer conserved sites in common (Figure 4d); however, the residues at its interaction sites appear highly conserved. These interaction sites are termed “mid-loop” and “end-helix” and allow the TNKS SAM domains to multimerize via head-to-tail interactions [18,19,20]. The ART domain shows a high degree of conservation around the NAD+-binding site (Figure 4e). The binding mode of NAD+ is therefore likely the same across TNKSs from different species. The residues of the donor loop (D-loop) are likewise entirely conserved. The D-loop was shown to adapt different conformations and it partially occupies the NAD+-binding pocket in crystal structures of apo-TNKSs [1]; it was speculated to be involved in the activity regulation of TNKSs [70]. Overall, these results suggest a remarkable level of conservation among TNKSs and imply that the molecular functions across TNKSs are likely conserved.

3.3. The Basic Molecular Characteristics of TNKSs Are Conserved

The human paralogs TNKS1 and TNKS2 are the most studied TNKS proteins. As shown above, the sequences of TNKSs across all metazoans are highly conserved. While the distant TNKS homolog in drosophila was functionally characterized through in vivo studies in terms of cellular function [31,71,72,73], we wanted to confirm that the molecular functions of distantly related TNKS homologs are conserved using in vitro studies with purified proteins. For this, we recombinantly produced TNKS proteins of a basal metazoan, the demosponge A. queenslandica (AqTNKS), comprising the ART domain, the ART-SAM domains and the ARC5 domain. Corresponding constructs from human TNKS1 were tested for comparison. First, we tested the poly-ADP-ribosylation activity of the ART and SAM-ART constructs by auto-modification (Figure 5a). The proteins were incubated with biotinylated NAD+ and the incorporation of biotin into the PAR chains was detected in Western blot using HRP conjugated streptavidin. A smear corresponding to the PARylation of the proteins was clearly detected for the SAM-ART constructs of AqTNKS and TNKS1, but not for the isolated ART domains. Moreover, to detect PAR chains that originate from PARylation activity during the recombinant expression, we used the PAR-binder ALC1 fused to nanoluciferase [49]. Similarly, PAR chains were detected only for the SAM-ART constructs, and incubation with PAR hydrolyzing human PARG resulted in loss of the smears. These results are in agreement with previous reports showing a significant loss of catalytic activity of human TNKS1/2 in absence of the SAM domain [7,19,20,21,70]. Therefore, the role as poly-ADP-ribosyltransferases and the SAM-dependent activation mechanism are likely conserved among TNKSs. Moreover, size-exclusion experiments with AqTNKS constructs showed formation of multimeric species for the SAM-ART protein but not the ART domain, which is in agreement with the reported multimerization function of the SAM domain in human TNKSs (Figure S6) [18,19,20,21].
To further demonstrate the high similarity between the ART domains of AqTNKS and TNKS1, we tested the thermal stabilization with established potent and selective TNKS inhibitors that were recently benchmarked [50]. This included adenosine-site binders IWR-1, G007-LK, OM-153 and compound 40 [74,75,76,77], nicotinamide-site binders XAV939, E7449, AZ6102 [10,78,79] and the dual-site binder NVP-TNKS656 [80]. For all inhibitors, the melting temperatures for both TNKS1 and AqTNKS were similarly increased (Figure 5b), underlining the structural similarity between these ART domains.
Next, we tested the binding of a canonical TBM peptide to the ARC5 domain constructs of AqTNKS and TNKS1. Guettler et al. developed an optimized peptide (sequence: REAGDGEE) that possesses high binding affinity to the ARC domains of human TNKSs [15,43]. We expressed the ARC5 domain constructs and the TBM peptide as fusion proteins with the fluorescent proteins CFP and YFP, respectively. To study the interaction of the ARC5 constructs with the TBM peptide, we measured FRET emissions upon mixing of the constructs (Figure 5c). Both ARC5 constructs from TNKS1 and AqTNKS interact with the TBM fused to YFP with similar binding affinities, while they showed no interaction with YFP lacking the TBM. In combination with the high conservation of the residues lining the ARC5 domain, these results suggest that the ability of TNKSs to bind the same TBM sequences is highly conserved.

3.4. Acquisition of Tankyrase Binding Partners throughout Evolution

Experimental data supporting bona fide interaction for many of the reported TNKS binding partners in humans exists, although only limited information is available for the interaction of proteins with TNKS orthologs in other metazoans [31]. Considering the strong conservation of TNKSs and the ability of AqTNKS to bind a peptide with the canonical TBM, it is reasonable to think that TNKSs throughout the metazoan lineage can bind this motif. The conservation of the TBMs in orthologs from binding partners of human TNKSs has been reported for selected binders, covering often only a very limited taxonomic breadth [10,23,24,25].
To gain insight into the conservation of TNKS binding partners, we selected a set of 21 experimentally confirmed TNKS interaction partners (Table 2) and determined the presence of TBMs in the orthologs of several metazoan species (Figure 6). Only motifs were considered that aligned with the position of TBMs from the human orthologs. In addition to the canonical TBM Rxx[ACGP]xGxx, we also marked the presence of the less constrained motifs Rx(4)Gxx and Rx(5)Gxx, which were reported for some TNKS interaction partners in humans [10,16,17], although the target may require multiple of such motifs to achieve sufficient avidity.
Overall, the analysis revealed that the conservation of corresponding TBMs compared to the human orthologs varies greatly for different proteins. A few of the analyzed orthologs seem to have acquired the TBMs relatively recently, such as GDP-mannose 4,6-dehydratase (GMD) or the telomere length regulator TRF1. Both do not harbor corresponding TBMs in orthologs of mice or in any of the earlier-diverging metazoans, with the exception of a possible TBM in the GMD ortholog from P. marinus. Interestingly, GMD is an otherwise highly conserved enzyme important in basic metabolic pathways, and orthologs of GMD are present in prokaryotes and eukaryotes [81]. While it was reported that GMD is not a direct target of PARylation by TNKSs and moreover appears to inhibit the catalytic activity of TNKSs [24], the exact function of the GMD-TNKS interaction is not well-understood at this point. In contrast, the interaction of TNKSs with Axin1 (Axin) appears highly conserved, as possible TBMs were found in orthologs from all species used for the analysis. Axin is a crucial component in the Wnt signaling pathway and interacts in humans with two TBMs with TNKSs; however, we omitted the non-canonical second motif from the analysis as the binding rules for it are not clear [16]. Earlier work by Feng et al. also confirmed interaction of TNKS and Axin in D. melanogaster [31]. TBMs in Golgin-45 were identified in orthologs from many of the analyzed species with exception of A. queenslandica. Golgin-45 is a Golgi-associated protein required for Golgi structure and protein transport. Zhang et al. identified Golgin-45 as a target for PARylation by TNKSs [82]. More recently, it was also shown to be responsible for TNKS1 localization to the Golgi apparatus during interphase [83].
Table 2. Overview of tankyrase interaction partners used for the analysis. The required arginine and glycine in each TBM is shown in bold and underlined. The residue number of the arginine in the first position of the TBM from the human ortholog is shown.
Table 2. Overview of tankyrase interaction partners used for the analysis. The required arginine and glycine in each TBM is shown in bold and underlined. The residue number of the arginine in the first position of the TBM from the human ortholog is shown.
TNKS Interaction PartnerRole or Cellular FunctionTBM in Human OrthologReference
AMOTHippo/YAP signaling77-RQEPQGQE[11]
ArpinCell migration213-REQGDGAE[84,85]
ATG9APexophagy233-RLPGLGEA[26]
Axin1 *Wnt signaling22-RPPVPGEE[10,14,16,31]
CASC3Pre-mRNA splicing146-RQSGDGQE[15,82]
Dicer1RNA interference656-RELPDGTF[86]
GMDGDP-mannose 4,6 dehydratase12-RGSGDGEM[24,87]
Golgin-45Golgi structure and protein maturation18-RGAGDGME[15,82,83]
IRAPVesicle trafficking96-RQSPDGAC[5,6,14,23]
Mcl-1Apoptosis regulation78-RPPPIGAE[15,88]
MERIT40DNA damage repair28-RSNPEGAE
48-RSEGEGEA
[15,89]
NKD2Wnt signaling16-RESPEGDS[86]
Notch2Notch signaling1726-RREPVGQD[86]
NuMAMitotic spindle assembly1743-RTQPDGTS[15,23,90]
PEX14Peroxisome homeostasis310-RMEVQGEE
350-RRGGDGQI
[26]
SH3BP2Osteoclast formation415-RSPPDGQS[7,15]
SH3BP5Epithelial lumen formation269-RGCGVGAE
368-RSECSGAS
[91]
SOX9Chondrocyte differentiation257-RPLPEGGR
271-RDVDIGEL
[92]
StriatinDendritic Ca2+ signaling302-RSAGDGTD[15,93]
TAB182Cytoskeletal maintenance1508-RPQPDGEA[15,23]
TRF1Telomere maintenance13-RGCADGRD[2,23,87]
USP25Ubiquitin protease1049-RTPADGR[15,25]
* The second non-canonical TBM from Axin1 is not shown here.
In our view, TNKSs are not likely to bind the orthologs that show an absence of corresponding TBMs, although it cannot be ruled out that some of them may still interact with TNKSs, for example due to the presence of non-canonical binding motifs or because of TBMs situated in different positions compared to the binding partner in human. Moreover, it is striking that many of the human TNKS binding partners appear to entirely lack corresponding orthologs in several species. For example, we could not identify orthologs of TAB182 in X. laevis or any earlier-diverging metazoans. Orthologs of AMOT are present in all analyzed vertebrates, while no orthologs could be identified in any of the invertebrates. Notably, all identified orthologs of AMOT possess a TBM corresponding to the human ortholog, indicating that an interaction with TNKSs may have evolved not long after the emergence of AMOT.
Due to the many TBMs present in the C-terminus of the prominent TNKS binding partner RNF146 [17], a different analysis was performed for this protein (Figure S7). RNF146 is an E3 ubiquitin ligase, and it is associated with PAR-dependent proteasomal degradation of TNKS target proteins [82,94,95,96,97]. All possible TBMs were mapped to the sequences of RNF146 orthologs, showing that many of these proteins possess possible TBMs in different numbers, configurations and positions compared to human RNF146. While it is difficult to predict which of these orthologs may interact with TNKSs, a recent study showed that in D. melanogaster, TNKS PARylation-dependent proteasomal degradation is mediated by RNF146 [98], although a direct interaction with TNKSs was not investigated.
The results obtained here lead to the overall picture that TNKSs themselves are highly conserved, while the interactions with binding partners seem to be far more variable throughout evolution as suggested in part by the varying presence of TBMs in orthologs and also because several of the binding partner orthologs appear to be simply not present in many species.

4. Discussion

4.1. Origin and Distribution of Tankyrases

The evolutionary transition from a single-celled urchoanozoan (ancestor of choanoflagellates and metazoans) to metazoans occurred together with the innovation and loss of gene families [60,99]. Proteins that originated during this transition and that were retained in the metazoan lineage are of interest for investigations about the evolutionary origin of metazoans, as these proteins may have been important factors for allowing the development of multicellular animals.
TNKSs seem to appear only sparsely in choanoflagellates but are vastly distributed throughout metazoans, which suggests that TNKSs may be important for multicellular development. However, the presence of TNKSs seems not to be obligatory for metazoan life, as evidenced by the reported loss of TNKS in C. elegans [52,53] and possibly in placozoans as shown in this work.
As part of the analysis, we identified a yet uncharacterized ARTD family member in the pathogenic amoeba N. fowleri, here termed NfTNKSL. Based on the sequence information available, NfTNKSL contains only a single ART domain, which shows more similarity to the ART domain of TNKSs than that of any other ARTD family member. It is tempting to speculate that NfTNKSL may resemble an ancestral TNKS precursor, which later obtained ANK and SAM domains through domain-shuffling, a common mechanism for the innovation of novel protein functions in the stem animal lineage [100,101,102].

4.2. Conservation of Tankyrases

We have shown that the sequences of TNKSs are highly conserved throughout metazoans, displaying similar levels of conservation to fundamentally important proteins such as cytochrome C. As it was expected, the degree of conservation is not equal throughout the TNKS sequence but corresponds to the structure of the protein. For example, regions thought to be flexible linkers appear less conserved, although some of them are likely important for the function of TNKSs [14].
Mapping of conserved regions to the TNKS structures showed a very high degree of conservation particularly around the TBM binding site of ARCs and the NAD+ binding site of the ART domain. Functional in vitro studies comparing human TNKS1 with AqTNKS from A. queenslandica, a sponge that diverged from other metazoans over 600 million years ago [103], confirmed that the basic molecular functions are conserved. We demonstrated that AqTNKS is a poly-ADP-ribosyltransferase that has strongly increased activity in the presence of the SAM domain, and that the ARC5 domain binds a canonical TBM with similar affinity compared to human TNKS1. The ability of AqTNKS to bind structurally diverse TNKS inhibitors that display high selectivity and nanomolar binding affinity against human TNKSs further highlights the high degree of conservation.

4.3. Acquisition of Tankyrase Binding Partners

From an analysis of the presence of TBMs in orthologs of human TNKS binding partners, the number of orthologs with corresponding TBMs decreased according to the species’ divergence from humans. Our analysis indicates that new TBMs formed throughout different animal lineages. This may leave the impression of a trend favoring accumulation of TNKS binding partners as the complexity of animal species increases; however, it is important to consider that the loss of TBMs may have occurred as well. The vast majority of TNKS binding partners were only experimentally confirmed in humans, which limits our analysis to a possible emergence of TBMs in orthologs of these binding partners; a comprehensive analysis for the loss in TBMs would require multiple experimentally verified TNKS binding partners in early diverging metazoan lineages. From the current results, we cannot exclude the possibility that the number of TNKS interaction partners for example in early diverging groups such as sponges or arthropods may be similar to the number found in humans.
Hub proteins that bind small linear motifs (SLiMs) were recently suggested to be grouped under the term of “linear motif-binding hubs” (LMB-hubs) [104]. With TBMs being SLiMs, and considering the functions of TNKSs, they certainly can be classified as LMB-hubs. Evolutionary analyses previously performed on other LMB-hub proteins established striking parallels also observed here for TNKSs: While the hub proteins are generally highly conserved, the presence of interacting SLiMs is often poorly retained in the orthologs of binding partners [105,106]. This may be the case because SLiMs can rapidly evolve de novo or similarly be lost again through small changes in amino acid sequences [107]. SLiMs are located on intrinsically disordered protein regions such as the protein termini, which are often far less mutationally constrained than structured regions—allowing rapid changes in sequence composition [108]. Furthermore, SLiM-hub interactions may allow a level of fine-adjustment for the regulation of binding partners through change of the SLiM sequence, resulting in changes in binding affinity to the hub protein [107].
Indeed, it has been reported that differences in the sequence of TBMs result in varied binding affinities to TNKSs [15,87]. A study by Eisemann et al. recently demonstrated that mutations in the TBM of TRF1 reduced the binding affinity and also resulted in lowered PARylation of TRF1 by TNKS1 in vitro. On an evolutionary level, the TNKS-TBM interaction may allow for TNKS binding partners to adapt their TBMs and thereby adjust the degree by which they are targeted by TNKS-mediated PARylation. This adjustment of binding may not only happen through changes in the TBM sequences but could also take place through evolution of additional TBMs in the TNKS binding partners or via their dimerization, which creates a complex relationship in terms of multivalent binding to the multimeric TNKS scaffolds.

4.4. Implications for the Development of Therapeutics Targeting Tankyrases

The cellular or physiological functions of TNKSs and the biological outcomes thereof can be attributed to the diverse assortment of binding partners from many different pathways that are regulated by TNKSs. While the structure and molecular functions were assessed to be highly conserved among TNKSs, the cellular functions of TNKSs may vary greatly depending on the species as different sets of TNKS binding partners may be regulated by TNKSs. For example, in D. melanogaster, we could identify only 11 orthologs of the 21 analyzed human TNKS binding partners, and only 4 of these orthologs harbor corresponding TBMs. Another example can be found in TRF1, a nuclear protein repressing telomere elongation which interacts through a TBM with TNKSs in humans. While it was reported that TNKSs have an important role in telomere maintenance through regulation of TRF1 in humans [2,3,4,109,110,111], the TBM is absent from the TRF1 ortholog in mice, where TNKSs were shown not to participate in the regulation of telomere length in multiple studies [28,112,113]. Although the attenuation of telomere elongation through inhibition of TNKSs could be a potential avenue for cancer therapeutics [114], mouse models would not adequately represent the potential efficacy of such drugs. On the other hand, inhibition of TNKSs was shown to reduce Wnt signaling activity through blocking the PARylation of Axin [10], which has a highly conserved TBM across species. Treatment with TNKS inhibitors leads to reduced proliferation of cancer cells in Wnt-dependent tumors [115] and mice are frequently used as model organisms in preclinical studies.
For the evaluation of TNKS inhibitors as therapeutic drugs, it should be considered that the cellular functions of TNKSs may differ greatly in humans and model organisms, which may limit for example how well the perceived efficacy or toxicity translates to humans.

5. Conclusions

The main points derived from this work are the following: First, we showed that TNKSs are distributed widely throughout metazoans, although they seem to have a pre-metazoan origin as indicated by their presence in some choanoflagellates. Second, our analyses indicate that TNKSs are remarkably conserved in terms of sequence, structure and molecular function. Last, the conservation of TBMs in TNKS binding partners and their orthologs appears to vary greatly, indicating a dynamic gain and possibly a loss of TBMs throughout evolution. Moreover, many orthologs of human TNKS binding partners appear to be missing in several of the species analyzed, which further contributes to the notion that TNKSs are likely regulating different sets of processes in different species.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biom12111688/s1, Table S1: Protein sequence identifiers from the analysis of tankyrase distribution and origin. Related to Figure 2; Figure S1: Phylogenetic analysis of tankyrases. Related to Figure 2; Figure S2: Predicted structure of the TNKS-like ankyrin-repeat domain from the archaeon Acidianus manzaensis. Related to Figure 2; Figure S3: Predicted structure of the TNKS-like ART domain from the amoeba Naegleria fowleri. Related to Figure 2. Figure S4: Conservation of TNKS ARC domain residues mapped to the structures. Related to Figure 4; Figure S5: The ARC3 domains in human and early diverging tankyrases are unlikely to bind TBM peptides. Related to Figure 4; Figure S6: Size-exclusion chromatography of TNKS1 and AqTNKS constructs. Related to Figure 5; Figure S7: Distribution of possible tankyrase binding motifs in orthologs of RNF146. Related to Figure 6; Data S1: Sequences of TNKSs for the construction of the phylogenetic trees. The sequences were truncated. TNKS sequences start from the ARC3 domain to match the likely incomplete N-terminal sequence from S. macrocollata. Related to Figure 2b; Data S2: Sequences from orthologs of PARP proteins, cytochrome C and seryl-tRNA synthetase in separate files. Files contain multiple sequence alignments (Clustal Omega) and corresponding percent identity matrices. Related to Figure 3; Data S3: Sequences used as input for multiple sequence alignment and analysis of position-specific conservation by ConSurf. Related to Figure 4a; Data S4: Position-specific amino acid variety and ConSurf grade mapped to the reference sequence of human TNKS2. Related to Figure 4a. Data S5: Sequences of TNKSs from multiple species and multiple sequence alignment of these sequences. Residues that were completely conserved were mapped to the domain structures of human TNKSs. Related to Figure 4b–e, Figure S4 and S5; Data S6: Files containing raw data of DSF and FRET experiments. Related to Figure 5b,c; Data S7: Sequences from orthologs of TNKS binding partners and multiple sequence alignments. The potential and corresponding TBMs are highlighted in the sequences and the multiple sequence alignments. Each set of orthologs is contained in a separate file. Related to Figure 6 and Figure S7.

Author Contributions

Conceptualization, L.L. and S.T.S.; validation, S.T.S., C.B. and A.G.-P.; formal analysis, S.T.S. and M.S.J.; investigation, S.T.S., C.B. and A.G.-P.; resources, L.L.; data curation, S.T.S. and M.S.J.; writing—original draft preparation, S.T.S.; writing—review and editing, S.T.S., C.B., A.G.-P., M.S.J. and L.L.; visualization, S.T.S. and C.B.; supervision, A.G.-P., M.S.J. and L.L.; project administration, L.L.; funding acquisition, L.L. All authors have read and agreed to the published version of the manuscript.

Funding

The work was funded by the Jane and Aatos Erkko Foundation.

Data Availability Statement

The data presented in this study are available in manuscript and in the Supplementary Materials.

Acknowledgments

Biocenter Oulu Structural Biology core facility, member of Biocenter Finland, Instruct-ERIC Centre Finland and FINStruct, as well as of “Proteomics and Protein Analysis” and Sequencing core facilities are gratefully acknowledged. We also acknowledge the support of the Biocenter Finland Bioinformatics platform node at Åbo Akademi University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Haikarainen, T.; Krauss, S.; Lehtiö, L. Tankyrases: Structure, Function and Therapeutic Implications in Cancer. Curr. Pharm. Des. 2014, 20, 6472–6488. [Google Scholar] [CrossRef] [Green Version]
  2. Smith, S.; Giriat, I.; Schmitt, A.; de Lange, T. Tankyrase, a Poly(ADP-Ribose) Polymerase at Human Telomeres. Science 1998, 282, 1484–1487. [Google Scholar] [CrossRef] [Green Version]
  3. Smith, S.; de Lange, T. Tankyrase Promotes Telomere Elongation in Human Cells. Curr. Biol. 2000, 10, 1299–1302. [Google Scholar] [CrossRef] [Green Version]
  4. Cook, B.D.; Dynek, J.N.; Chang, W.; Shostak, G.; Smith, S. Role for the Related Poly(ADP-Ribose) Polymerases Tankyrase 1 and 2 at Human Telomeres. Mol. Cell. Biol. 2002, 22, 332–342. [Google Scholar] [CrossRef] [Green Version]
  5. Yeh, T.-Y.J.; Sbodio, J.I.; Tsun, Z.-Y.; Luo, B.; Chi, N.-W. Insulin-Stimulated Exocytosis of GLUT4 Is Enhanced by IRAP and Its Partner Tankyrase. Biochem. J. 2007, 402, 279–290. [Google Scholar] [CrossRef] [Green Version]
  6. Chi, N.-W.; Lodish, H.F. Tankyrase Is a Golgi-Associated Mitogen-Activated Protein Kinase Substrate That Interacts with IRAP in GLUT4 Vesicles. J. Biol. Chem. 2000, 275, 38437–38444. [Google Scholar] [CrossRef] [Green Version]
  7. Levaot, N.; Voytyuk, O.; Dimitriou, I.; Sircoulomb, F.; Chandrakumar, A.; Deckert, M.; Krzyzanowski, P.M.; Scotter, A.; Gu, S.; Janmohamed, S.; et al. Loss of Tankyrase-Mediated Destruction of 3BP2 Is the Underlying Pathogenic Mechanism of Cherubism. Cell 2011, 147, 1324–1339. [Google Scholar] [CrossRef] [Green Version]
  8. Chang, P.; Coughlin, M.; Mitchison, T.J. Tankyrase-1 Polymerization of Poly(ADP-Ribose) Is Required for Spindle Structure and Function. Nat. Cell Biol. 2005, 7, 1133–1139. [Google Scholar] [CrossRef]
  9. Chang, P.; Coughlin, M.; Mitchison, T.J. Interaction between Poly(ADP-Ribose) and NuMA Contributes to Mitotic Spindle Pole Assembly. Mol. Biol. Cell 2009, 20, 4575–4585. [Google Scholar] [CrossRef] [Green Version]
  10. Huang, S.-M.A.; 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]
  11. Wang, W.; Li, N.; Li, X.; Tran, M.K.; Han, X.; Chen, J. Tankyrase Inhibitors Target YAP by Stabilizing Angiomotin Family Proteins. Cell Rep. 2015, 13, 524–532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Chiang, Y.J.; Hsiao, S.J.; Yver, D.; Cushman, S.W.; Tessarollo, L.; Smith, S.; Hodes, R.J. Tankyrase 1 and Tankyrase 2 Are Essential but Redundant for Mouse Embryonic Development. PLOS ONE 2008, 3, e2639. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Lüscher, B.; Ahel, I.; Altmeyer, M.; Ashworth, A.; Bai, P.; Chang, P.; Cohen, M.; Corda, D.; Dantzer, F.; Daugherty, M.D.; et al. ADP-ribosyltransferases, an Update on Function and Nomenclature. FEBS J. 2021. ahead of print. [Google Scholar] [CrossRef] [PubMed]
  14. Eisemann, T.; McCauley, M.; Langelier, M.-F.; Gupta, K.; Roy, S.; Van Duyne, G.D.; Pascal, J.M. Tankyrase-1 Ankyrin Repeats Form an Adaptable Binding Platform for Targets of ADP-Ribose Modification. Structure 2016, 24, 1679–1692. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Guettler, S.; LaRose, J.; Petsalaki, E.; Gish, G.; Scotter, A.; Pawson, T.; Rottapel, R.; Sicheri, F. Structural Basis and Sequence Rules for Substrate Recognition by Tankyrase Explain the Basis for Cherubism Disease. Cell 2011, 147, 1340–1354. [Google Scholar] [CrossRef] [Green Version]
  16. Morrone, S.; Cheng, Z.; Moon, R.T.; Cong, F.; Xu, W. Crystal Structure of a Tankyrase-Axin Complex and Its Implications for Axin Turnover and Tankyrase Substrate Recruitment. Proc. Natl. Acad. Sci. USA 2012, 109, 1500–1505. [Google Scholar] [CrossRef] [Green Version]
  17. DaRosa, P.A.; Klevit, R.E.; Xu, W. Structural Basis for Tankyrase-RNF146 Interaction Reveals Noncanonical Tankyrase-Binding Motifs. Protein Sci. 2018, 27, 1057–1067. [Google Scholar] [CrossRef] [Green Version]
  18. DaRosa, P.A.; Ovchinnikov, S.; Xu, W.; Klevit, R.E. Structural Insights into SAM Domain-Mediated Tankyrase Oligomerization. Protein Sci. 2016, 25, 1744–1752. [Google Scholar] [CrossRef] [Green Version]
  19. Mariotti, L.; Templeton, C.M.; Ranes, M.; Paracuellos, P.; Cronin, N.; Beuron, F.; Morris, E.; Guettler, S. Tankyrase Requires SAM Domain-Dependent Polymerization to Support Wnt-β-Catenin Signaling. Mol. Cell 2016, 63, 498–513. [Google Scholar] [CrossRef] [Green Version]
  20. Riccio, A.A.; McCauley, M.; Langelier, M.-F.; Pascal, J.M. Tankyrase Sterile α Motif Domain Polymerization Is Required for Its Role in Wnt Signaling. Structure 2016, 24, 1573–1581. [Google Scholar] [CrossRef]
  21. Rycker, M.D.; Price, C.M. Tankyrase Polymerization Is Controlled by Its Sterile Alpha Motif and Poly(ADP-Ribose) Polymerase Domains. Mol. Cell. Biol. 2004, 24, 9802–9812. [Google Scholar] [CrossRef] [Green Version]
  22. Kaminker, P.G.; Kim, S.-H.; Taylor, R.D.; Zebarjadian, Y.; Funk, W.D.; Morin, G.B.; Yaswen, P.; Campisi, J. TANK2, a New TRF1-Associated Poly(ADP-Ribose) Polymerase, Causes Rapid Induction of Cell Death upon Overexpression. J. Biol. Chem. 2001, 276, 35891–35899. [Google Scholar] [CrossRef] [Green Version]
  23. Sbodio, J.I.; Chi, N.-W. Identification of a Tankyrase-Binding Motif Shared by IRAP, TAB182, and Human TRF1 but Not Mouse TRF1. J. Biol. Chem. 2002, 277, 31887–31892. [Google Scholar] [CrossRef] [Green Version]
  24. Bisht, K.K.; Dudognon, C.; Chang, W.G.; Sokol, E.S.; Ramirez, A.; Smith, S. GDP-Mannose-4,6-Dehydratase Is a Cytosolic Partner of Tankyrase 1 That Inhibits Its Poly(ADP-Ribose) Polymerase Activity. Mol. Cell. Biol. 2012, 32, 3044–3053. [Google Scholar] [CrossRef] [Green Version]
  25. Xu, D.; Liu, J.; Fu, T.; Shan, B.; Qian, L.; Pan, L.; Yuan, J. USP25 Regulates Wnt Signaling by Controlling the Stability of Tankyrases. Genes Dev. 2017, 31, 1024–1035. [Google Scholar] [CrossRef] [Green Version]
  26. Li, X.; Han, H.; Zhou, M.-T.; Yang, B.; Ta, A.P.; Li, N.; Chen, J.; Wang, W. Proteomic Analysis of the Human Tankyrase Protein Interaction Network Reveals Its Role in Pexophagy. Cell Rep. 2017, 20, 737–749. [Google Scholar] [CrossRef] [Green Version]
  27. Smith, S.; de Lange, T. Cell Cycle Dependent Localization of the Telomeric PARP, Tankyrase, to Nuclear Pore Complexes and Centrosomes. J. Cell Sci. 1999, 112, 3649–3656. [Google Scholar] [CrossRef]
  28. Muramatsu, Y.; Ohishi, T.; Sakamoto, M.; Tsuruo, T.; Seimiya, H. Cross-Species Difference in Telomeric Function of Tankyrase 1. Cancer Sci. 2007, 98, 850–857. [Google Scholar] [CrossRef]
  29. Qian, L.; Mahaffey, J.P.; Alcorn, H.L.; Anderson, K.V. Tissue-Specific Roles of Axin2 in the Inhibition and Activation of Wnt Signaling in the Mouse Embryo. Proc. Natl. Acad. Sci. USA 2011, 108, 8692–8697. [Google Scholar] [CrossRef] [Green Version]
  30. Kuusela, S.; Wang, H.; Wasik, A.A.; Suleiman, H.; Lehtonen, S. Tankyrase Inhibition Aggravates Kidney Injury in the Absence of CD2AP. Cell Death Dis. 2016, 7, e2302. [Google Scholar] [CrossRef]
  31. Feng, Y.; Li, X.; Ray, L.; Song, H.; Qu, J.; Lin, S.; Lin, X. The Drosophila Tankyrase Regulates Wg Signaling Depending on the Concentration of Daxin. Cell. Signal. 2014, 26, 1717–1724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Hsiao, S.J.; Smith, S. Tankyrase Function at Telomeres, Spindle Poles, and Beyond. Biochimie 2008, 90, 83–92. [Google Scholar] [CrossRef] [PubMed]
  33. Kim, M.K. Novel Insight into the Function of Tankyrase. Oncol. Lett. 2018, 16, 6895–6902. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Zamudio-Martinez, E.; Herrera-Campos, A.B.; Muñoz, A.; Rodríguez-Vargas, J.M.; Oliver, F.J. Tankyrases as Modulators of Pro-Tumoral Functions: Molecular Insights and Therapeutic Opportunities. J. Exp. Clin. Cancer Res. CR 2021, 40, 144. [Google Scholar] [CrossRef] [PubMed]
  35. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic Local Alignment Search Tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  36. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; McWilliam, H.; Remmert, M.; Söding, J.; et al. Fast, Scalable Generation of High-quality Protein Multiple Sequence Alignments Using Clustal Omega. Mol. Syst. Biol. 2011, 7, 539. [Google Scholar] [CrossRef]
  37. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef]
  38. Jones, D.T.; Taylor, W.R.; Thornton, J.M. The Rapid Generation of Mutation Data Matrices from Protein Sequences. Bioinformatics 1992, 8, 275–282. [Google Scholar] [CrossRef]
  39. Berman, H.M. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [Green Version]
  40. Baek, M.; DiMaio, F.; Anishchenko, I.; Dauparas, J.; Ovchinnikov, S.; Lee, G.R.; Wang, J.; Cong, Q.; Kinch, L.N.; Schaeffer, R.D.; et al. Accurate Prediction of Protein Structures and Interactions Using a Three-Track Neural Network. Science 2021, 373, 871–876. [Google Scholar] [CrossRef]
  41. Ashkenazy, H.; Abadi, S.; Martz, E.; Chay, O.; Mayrose, I.; Pupko, T.; Ben-Tal, N. ConSurf 2016: An Improved Methodology to Estimate and Visualize Evolutionary Conservation in Macromolecules. Nucleic Acids Res. 2016, 44, W344–W350. [Google Scholar] [CrossRef] [Green Version]
  42. Blum, M.; Chang, H.-Y.; Chuguransky, S.; Grego, T.; Kandasaamy, S.; Mitchell, A.; Nuka, G.; Paysan-Lafosse, T.; Qureshi, M.; Raj, S.; et al. The InterPro Protein Families and Domains Database: 20 Years On. Nucleic Acids Res. 2021, 49, D344–D354. [Google Scholar] [CrossRef]
  43. Sowa, S.T.; Vela-Rodríguez, C.; Galera-Prat, A.; Cázares-Olivera, M.; Prunskaite-Hyyryläinen, R.; Ignatev, A.; Lehtiö, L. A FRET-Based High-Throughput Screening Platform for the Discovery of Chemical Probes Targeting the Scaffolding Functions of Human Tankyrases. Sci. Rep. 2020, 10, 12357. [Google Scholar] [CrossRef]
  44. Sowa, S.T.; Lehtiö, L. The Zinc-Binding Motif in Tankyrases Is Required for the Structural Integrity of the Catalytic ADP-Ribosyltransferase Domain. Open Biol. 2022, 12, 210365. [Google Scholar] [CrossRef]
  45. Jeong, J.-Y.; Yim, H.-S.; Ryu, J.-Y.; Lee, H.S.; Lee, J.-H.; Seen, D.-S.; Kang, S.G. One-Step Sequence- and Ligation-Independent Cloning as a Rapid and Versatile Cloning Method for Functional Genomics Studies. Appl. Environ. Microbiol. 2012, 78, 5440–5443. [Google Scholar] [CrossRef] [Green Version]
  46. Hynes, M.F.; Quandt, J.; O’Connell, M.P.; Pühler, A. Direct Selection for Curing and Deletion of Rhizobium Plasmids Using Transposons Carrying the Bacillus Subtilis SacB Gene. Gene 1989, 78, 111–120. [Google Scholar] [CrossRef]
  47. Reyrat, J.-M.; Pelicic, V.; Gicquel, B.; Rappuoli, R. Counterselectable Markers: Untapped Tools for Bacterial Genetics and Pathogenesis. Infect. Immun. 1998, 66, 4011–4017. [Google Scholar] [CrossRef]
  48. Rial, D.V.; Ceccarelli, E.A. Removal of DnaK Contamination during Fusion Protein Purifications. Protein Expr. Purif. 2002, 25, 503–507. [Google Scholar] [CrossRef]
  49. Sowa, S.T.; Galera-Prat, A.; Wazir, S.; Alanen, H.I.; Maksimainen, M.M.; Lehtiö, L. A Molecular Toolbox for ADP-Ribosyl Binding Proteins. Cell Rep. Methods 2021, 1, 100121. [Google Scholar] [CrossRef]
  50. Brinch, S.A.; Amundsen-Isaksen, E.; Espada, S.; Hammarstrom, C.; Aizenshtadt, A.; Olsen, P.A.; Holmen, L.; Hoyem, M.; Scholz, H.; Grodeland, G.; et al. The Tankyrase Inhibitor OM-153 Demonstrates Anti-Tumor Efficacy and a Therapeutic Window in Mouse Models. Cancer Res. Commun. 2022, 2, 233–245. [Google Scholar] [CrossRef]
  51. Song, Y.; Rodgers, V.G.J.; Schultz, J.S.; Liao, J. Protein Interaction Affinity Determination by Quantitative FRET Technology. Biotechnol. Bioeng. 2012, 109, 2875–2883. [Google Scholar] [CrossRef] [PubMed]
  52. Citarelli, M.; Teotia, S.; Lamb, R.S. Evolutionary History of the Poly(ADP-Ribose) Polymerase Gene Family in Eukaryotes. BMC Evol. Biol. 2010, 10, 308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Perina, D.; Mikoč, A.; Ahel, J.; Ćetković, H.; Žaja, R.; Ahel, I. Distribution of Protein Poly(ADP-Ribosyl)Ation Systems across All Domains of Life. DNA Repair 2014, 23, 4–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Laumer, C.E.; Gruber-Vodicka, H.; Hadfield, M.G.; Pearse, V.B.; Riesgo, A.; Marioni, J.C.; Giribet, G. Support for a Clade of Placozoa and Cnidaria in Genes with Minimal Compositional Bias. eLife 2018, 7, e36278. [Google Scholar] [CrossRef] [PubMed]
  55. Srivastava, M.; Begovic, E.; Chapman, J.; Putnam, N.H.; Hellsten, U.; Kawashima, T.; Kuo, A.; Mitros, T.; Salamov, A.; Carpenter, M.L.; et al. The Trichoplax Genome and the Nature of Placozoans. Nature 2008, 454, 955–960. [Google Scholar] [CrossRef] [Green Version]
  56. Miller, D.J.; Ball, E.E. Animal Evolution: Trichoplax, Trees, and Taxonomic Turmoil. Curr. Biol. 2008, 18, R1003–R1005. [Google Scholar] [CrossRef] [Green Version]
  57. Carr, M.; Leadbeater, B.S.C.; Hassan, R.; Nelson, M.; Baldauf, S.L. Molecular Phylogeny of Choanoflagellates, the Sister Group to Metazoa. Proc. Natl. Acad. Sci. USA 2008, 105, 16641–16646. [Google Scholar] [CrossRef] [Green Version]
  58. King, N.; Westbrook, M.J.; Young, S.L.; Kuo, A.; Abedin, M.; Chapman, J.; Fairclough, S.; Hellsten, U.; Isogai, Y.; Letunic, I.; et al. The Genome of the Choanoflagellate Monosiga Brevicollis and the Origin of Metazoans. Nature 2008, 451, 783–788. [Google Scholar] [CrossRef] [Green Version]
  59. Fairclough, S.R.; Chen, Z.; Kramer, E.; Zeng, Q.; Young, S.; Robertson, H.M.; Begovic, E.; Richter, D.J.; Russ, C.; Westbrook, M.J.; et al. Premetazoan Genome Evolution and the Regulation of Cell Differentiation in the Choanoflagellate Salpingoeca Rosetta. Genome Biol. 2013, 14, R15. [Google Scholar] [CrossRef] [Green Version]
  60. Richter, D.J.; Fozouni, P.; Eisen, M.B.; King, N. Gene Family Innovation, Conservation and Loss on the Animal Stem Lineage. eLife 2018, 7, e34226. [Google Scholar] [CrossRef]
  61. Nagy, A.; Patthy, L. Reassessing Domain Architecture Evolution of Metazoan Proteins: The Contribution of Different Evolutionary Mechanisms. Genes 2011, 2, 578–598. [Google Scholar] [CrossRef] [Green Version]
  62. Mosavi, L.K.; Cammett, T.J.; Desrosiers, D.C.; Peng, Z.-Y. The Ankyrin Repeat as Molecular Architecture for Protein Recognition. Protein Sci. Publ. Protein Soc. 2004, 13, 1435–1448. [Google Scholar] [CrossRef] [Green Version]
  63. Lehtiö, L.; Collins, R.; van den Berg, S.; Johansson, A.; Dahlgren, L.-G.; Hammarström, M.; Helleday, T.; Holmberg-Schiavone, L.; Karlberg, T.; Weigelt, J. Zinc Binding Catalytic Domain of Human Tankyrase 1. J. Mol. Biol. 2008, 379, 136–145. [Google Scholar] [CrossRef] [Green Version]
  64. Baba, M.L.; Darga, L.L.; Goodman, M.; Czelusniak, J. Evolution of Cytochromec Investigated by the Maximum Parsimony Method. J. Mol. Evol. 1981, 17, 197–213. [Google Scholar] [CrossRef]
  65. Härtlein, M.; Cusack, S. Structure, Function and Evolution of Seryl-TRNA Synthetases: Implications for the Evolution of Aminoacyl-TRNA Synthetases and the Genetic Code. J. Mol. Evol. 1995, 40, 519–530. [Google Scholar] [CrossRef]
  66. Daugherty, M.D.; Young, J.M.; Kerns, J.A.; Malik, H.S. Rapid Evolution of PARP Genes Suggests a Broad Role for ADP-Ribosylation in Host-Virus Conflicts. PLoS Genet. 2014, 10. [Google Scholar] [CrossRef] [Green Version]
  67. Pascal, J.M. The Comings and Goings of PARP-1 in Response to DNA Damage. DNA Repair 2018, 71, 177–182. [Google Scholar] [CrossRef]
  68. Haikarainen, T.; Narwal, M.; Joensuu, P.; Lehtiö, L. Evaluation and Structural Basis for the Inhibition of Tankyrases by PARP Inhibitors. ACS Med. Chem. Lett. 2014, 5, 18–22. [Google Scholar] [CrossRef] [Green Version]
  69. Seimiya, H.; Muramatsu, Y.; Smith, S.; Tsuruo, T. Functional Subdomain in the Ankyrin Domain of Tankyrase 1 Required for Poly(ADP-Ribosyl)Ation of TRF1 and Telomere Elongation. Mol. Cell. Biol. 2004, 24, 1944–1955. [Google Scholar] [CrossRef] [Green Version]
  70. Fan, C.; Yarravarapu, N.; Chen, H.; Kulak, O.; Dasari, P.; Herbert, J.; Yamaguchi, K.; Lum, L.; Zhang, X. Regulation of Tankyrase Activity by a Catalytic Domain Dimer Interface. Biochem. Biophys. Res. Commun. 2018, 503, 1780–1785. [Google Scholar] [CrossRef]
  71. Cho-Park, P.F.; Steller, H. Proteasome Regulation by ADP-Ribosylation. Cell 2013, 153, 614–627. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Feng, Y.; Li, Z.; Lv, L.; Du, A.; Lin, Z.; Ye, X.; Lin, Y.; Lin, X. Tankyrase Regulates Apoptosis by Activating JNK Signaling in Drosophila. Biochem. Biophys. Res. Commun. 2018, 503, 2234–2239. [Google Scholar] [CrossRef] [PubMed]
  73. Wang, Z.; Tian, A.; Benchabane, H.; Tacchelly-Benites, O.; Yang, E.; Nojima, H.; Ahmed, Y. The ADP-Ribose Polymerase Tankyrase Regulates Adult Intestinal Stem Cell Proliferation during Homeostasis in Drosophila. Development 2016, 143, 1710–1720. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Bregman, H.; Chakka, N.; Guzman-Perez, A.; Gunaydin, H.; Gu, Y.; Huang, X.; Berry, V.; Liu, J.; Teffera, Y.; Huang, L.; et al. Discovery of Novel, Induced-Pocket Binding Oxazolidinones as Potent, Selective, and Orally Bioavailable Tankyrase Inhibitors. J. Med. Chem. 2013, 56, 4320–4342. [Google Scholar] [CrossRef]
  75. Chen, B.; Dodge, M.E.; Tang, W.; Lu, J.; Ma, Z.; Fan, C.-W.; Wei, S.; Hao, W.; Kilgore, J.; Williams, N.S.; et al. Small Molecule-Mediated Disruption of Wnt-Dependent Signaling in Tissue Regeneration and Cancer. Nat. Chem. Biol. 2009, 5, 100–107. [Google Scholar] [CrossRef] [Green Version]
  76. Leenders, R.G.G.; Brinch, S.A.; Sowa, S.T.; Amundsen-Isaksen, E.; Galera-Prat, A.; Murthy, S.; Aertssen, S.; Smits, J.N.; Nieczypor, P.; Damen, E.; et al. Development of a 1,2,4-Triazole-Based Lead Tankyrase Inhibitor: Part II. J. Med. Chem. 2021, 64, 17936–17949. [Google Scholar] [CrossRef]
  77. Voronkov, A.; Holsworth, D.D.; Waaler, J.; Wilson, S.R.; Ekblad, B.; Perdreau-Dahl, H.; Dinh, H.; Drewes, G.; Hopf, C.; Morth, J.P.; et al. Structural Basis and SAR for G007-LK, a Lead Stage 1,2,4-Triazole Based Specific Tankyrase 1/2 Inhibitor. J. Med. Chem. 2013, 56, 3012–3023. [Google Scholar] [CrossRef]
  78. 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]
  79. McGonigle, S.; Chen, Z.; Wu, J.; Chang, P.; Kolber-Simonds, D.; Ackermann, K.; Twine, N.C.; Shie, J.-L.; Miu, J.T.; Huang, K.-C.; et al. E7449: A Dual Inhibitor of PARP1/2 and Tankyrase1/2 Inhibits Growth of DNA Repair Deficient Tumors and Antagonizes Wnt Signaling. Oncotarget 2015, 6, 41307–41323. [Google Scholar] [CrossRef]
  80. Shultz, M.D.; Cheung, A.K.; Kirby, C.A.; Firestone, B.; Fan, J.; Chen, C.H.-T.; Chen, Z.; Chin, D.N.; DiPietro, L.; Fazal, A.; et al. Identification of NVP-TNKS656: The Use of Structure–Efficiency Relationships To Generate a Highly Potent, Selective, and Orally Active Tankyrase Inhibitor. J. Med. Chem. 2013, 56, 6495–6511. [Google Scholar] [CrossRef]
  81. Somoza, J.R.; Menon, S.; Schmidt, H.; Joseph-McCarthy, D.; Dessen, A.; Stahl, M.L.; Somers, W.S.; Sullivan, F.X. Structural and Kinetic Analysis of Escherichia Coli GDP-Mannose 4,6 Dehydratase Provides Insights into the Enzyme’s Catalytic Mechanism and Regulation by GDP-Fucose. Structure 2000, 8, 123–135. [Google Scholar] [CrossRef] [Green Version]
  82. Zhang, Y.; Liu, S.; Mickanin, C.; Feng, Y.; Charlat, O.; Michaud, G.A.; Schirle, M.; Shi, X.; Hild, M.; Bauer, A.; et al. RNF146 Is a Poly(ADP-Ribose)-Directed E3 Ligase That Regulates Axin Degradation and Wnt Signalling. Nat. Cell Biol. 2011, 13, 623–629. [Google Scholar] [CrossRef] [PubMed]
  83. Yue, X.; Tiwari, N.; Zhu, L.; Ngo, H.D.T.; Lim, J.-M.; Gim, B.; Jing, S.; Wang, Y.; Qian, Y.; Lee, I. Tankyrase-1-Mediated Degradation of Golgin45 Regulates Glycosyltransferase Trafficking and Protein Glycosylation in Rab2-GTP-Dependent Manner. Commun. Biol. 2021, 4, 1370. [Google Scholar] [CrossRef] [PubMed]
  84. Fetics, S.; Thureau, A.; Campanacci, V.; Aumont-Nicaise, M.; Dang, I.; Gautreau, A.; Pérez, J.; Cherfils, J. Hybrid Structural Analysis of the Arp2/3 Regulator Arpin Identifies Its Acidic Tail as a Primary Binding Epitope. Structure 2016, 24, 252–260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Simanov, G.; Dang, I.; Fokin, A.I.; Oguievetskaia, K.; Campanacci, V.; Cherfils, J.; Gautreau, A.M. Arpin Regulates Migration Persistence by Interacting with Both Tankyrases and the Arp2/3 Complex. Int. J. Mol. Sci. 2021, 22, 4115. [Google Scholar] [CrossRef]
  86. Bhardwaj, A.; Yang, Y.; Ueberheide, B.; Smith, S. Whole Proteome Analysis of Human Tankyrase Knockout Cells Reveals Targets of Tankyrase-Mediated Degradation. Nat. Commun. 2017, 8, 1–13. [Google Scholar] [CrossRef] [Green Version]
  87. Eisemann, T.; Langelier, M.-F.; Pascal, J.M. Structural and Functional Analysis of Parameters Governing Tankyrase-1 Interaction with Telomeric Repeat-Binding Factor 1 and GDP-Mannose 4,6-Dehydratase. J. Biol. Chem. 2019, 294, 14574–14590. [Google Scholar] [CrossRef]
  88. Bae, J.; Donigian, J.R.; Hsueh, A.J.W. Tankyrase 1 Interacts with Mcl-1 Proteins and Inhibits Their Regulation of Apoptosis. J. Biol. Chem. 2003, 278, 5195–5204. [Google Scholar] [CrossRef] [Green Version]
  89. Okamoto, K.; Ohishi, T.; Kuroiwa, M.; Iemura, S.; Natsume, T.; Seimiya, H. MERIT40-Dependent Recruitment of Tankyrase to Damaged DNA and Its Implication for Cell Sensitivity to DNA-Damaging Anticancer Drugs. Oncotarget 2018, 9, 35844–35855. [Google Scholar] [CrossRef] [Green Version]
  90. Chang, W.; Dynek, J.N.; Smith, S. NuMA Is a Major Acceptor of Poly(ADP-Ribosyl)Ation by Tankyrase 1 in Mitosis. Biochem. J. 2005, 391, 177–184. [Google Scholar] [CrossRef]
  91. Chandrakumar, A.A.; Coyaud, É.; Marshall, C.B.; Ikura, M.; Raught, B.; Rottapel, R. Tankyrase Regulates Epithelial Lumen Formation via Suppression of Rab11 GEFs. J. Cell Biol. 2021, 220, e202008037. [Google Scholar] [CrossRef]
  92. Kim, S.; Han, S.; Kim, Y.; Kim, H.-S.; Gu, Y.-R.; Kang, D.; Cho, Y.; Kim, H.; Lee, J.; Seo, Y.; et al. Tankyrase Inhibition Preserves Osteoarthritic Cartilage by Coordinating Cartilage Matrix Anabolism via Effects on SOX9 PARylation. Nat. Commun. 2019, 10, 4898. [Google Scholar] [CrossRef] [Green Version]
  93. Lahav-Ariel, L.; Caspi, M.; Nadar-Ponniah, P.T.; Zelikson, N.; Hofmann, I.; Hanson, K.K.; Franke, W.W.; Sklan, E.H.; Avraham, K.B.; Rosin-Arbesfeld, R. Striatin Is a Novel Modulator of Cell Adhesion. FASEB J. 2019, 33, 4729–4740. [Google Scholar] [CrossRef]
  94. Zhou, Z.; Chan, C.H.; Xiao, Z.; Tan, E. Ring Finger Protein 146/Iduna Is a Poly(ADP-Ribose) Polymer Binding and PARsylation Dependent E3 Ubiquitin Ligase. Cell Adhes. Migr. 2011, 5, 463–471. [Google Scholar] [CrossRef]
  95. Callow, M.G.; Tran, H.; Phu, L.; Lau, T.; Lee, J.; Sandoval, W.N.; Liu, P.S.; Bheddah, S.; Tao, J.; Lill, J.R.; et al. Ubiquitin Ligase RNF146 Regulates Tankyrase and Axin to Promote Wnt Signaling. PloS ONE 2011, 6, e22595. [Google Scholar] [CrossRef] [Green Version]
  96. Wang, Z.; Michaud, G.A.; Cheng, Z.; Zhang, Y.; Hinds, T.R.; Fan, E.; Cong, F.; Xu, W. Recognition of the Iso-ADP-Ribose Moiety in Poly(ADP-Ribose) by WWE Domains Suggests a General Mechanism for Poly(ADP-Ribosyl)Ation-Dependent Ubiquitination. Genes Dev. 2012, 26, 235–240. [Google Scholar] [CrossRef] [Green Version]
  97. DaRosa, P.A.; Wang, Z.; Jiang, X.; Pruneda, J.N.; Cong, F.; Klevit, R.E.; Xu, W. Allosteric Activation of the RNF146 Ubiquitin Ligase by a Poly(ADP-Ribosyl)Ation Signal. Nature 2015, 517, 223–226. [Google Scholar] [CrossRef] [Green Version]
  98. Wang, Z.; Tacchelly-Benites, O.; Noble, G.P.; Johnson, M.K.; Gagné, J.-P.; Poirier, G.G.; Ahmed, Y. A Context-Dependent Role for the RNF146 Ubiquitin Ligase in Wingless/Wnt Signaling in Drosophila. Genetics 2019, 211, 913–923. [Google Scholar] [CrossRef]
  99. Paps, J.; Holland, P.W.H. Reconstruction of the Ancestral Metazoan Genome Reveals an Increase in Genomic Novelty. Nat. Commun. 2018, 9, 1730. [Google Scholar] [CrossRef] [Green Version]
  100. Ekman, D.; Björklund, Å.K.; Elofsson, A. Quantification of the Elevated Rate of Domain Rearrangements in Metazoa. J. Mol. Biol. 2007, 372, 1337–1348. [Google Scholar] [CrossRef]
  101. López-Escardó, D.; Grau-Bové, X.; Guillaumet-Adkins, A.; Gut, M.; Sieracki, M.E.; Ruiz-Trillo, I. Reconstruction of Protein Domain Evolution Using Single-Cell Amplified Genomes of Uncultured Choanoflagellates Sheds Light on the Origin of Animals. Philos. Trans. R. Soc. B Biol. Sci. 2019, 374, 20190088. [Google Scholar] [CrossRef] [PubMed]
  102. Tordai, H.; Nagy, A.; Farkas, K.; Bányai, L.; Patthy, L. Modules, Multidomain Proteins and Organismic Complexity: Mobile Domains. FEBS J. 2005, 272, 5064–5078. [Google Scholar] [CrossRef] [PubMed]
  103. Srivastava, M.; Simakov, O.; Chapman, J.; Fahey, B.; Gauthier, M.E.A.; Mitros, T.; Richards, G.S.; Conaco, C.; Dacre, M.; Hellsten, U.; et al. The Amphimedon Queenslandica Genome and the Evolution of Animal Complexity. Nature 2010, 466, 720–726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Jespersen, N.; Barbar, E. Emerging Features of Linear Motif-Binding Hub Proteins. Trends Biochem. Sci. 2020, 45, 375–384. [Google Scholar] [CrossRef]
  105. Davey, N.E.; Van Roey, K.; Weatheritt, R.J.; Toedt, G.; Uyar, B.; Altenberg, B.; Budd, A.; Diella, F.; Dinkel, H.; Gibson, T.J. Attributes of Short Linear Motifs. Mol BioSyst 2012, 8, 268–281. [Google Scholar] [CrossRef]
  106. Van Roey, K.; Uyar, B.; Weatheritt, R.J.; Dinkel, H.; Seiler, M.; Budd, A.; Gibson, T.J.; Davey, N.E. Short Linear Motifs: Ubiquitous and Functionally Diverse Protein Interaction Modules Directing Cell Regulation. Chem. Rev. 2014, 114, 6733–6778. [Google Scholar] [CrossRef] [PubMed]
  107. Davey, N.E.; Cyert, M.S.; Moses, A.M. Short Linear Motifs - Ex Nihilo Evolution of Protein Regulation. Cell Commun. Signal. CCS 2015, 13, 43. [Google Scholar] [CrossRef] [Green Version]
  108. van der Lee, R.; Buljan, M.; Lang, B.; Weatheritt, R.J.; Daughdrill, G.W.; Dunker, A.K.; Fuxreiter, M.; Gough, J.; Gsponer, J.; Jones, D.T.; et al. Classification of Intrinsically Disordered Regions and Proteins. Chem. Rev. 2014, 114, 6589–6631. [Google Scholar] [CrossRef]
  109. Dynek, J.N.; Smith, S. Resolution of Sister Telomere Association Is Required for Progression Through Mitosis. Science 2004, 304, 97–100. [Google Scholar] [CrossRef]
  110. Ha, G.-H.; Kim, H.-S.; Go, H.; Lee, H.; Seimiya, H.; Chung, D.H.; Lee, C.-W. Tankyrase-1 Function at Telomeres and during Mitosis Is Regulated by Polo-like Kinase-1-Mediated Phosphorylation. Cell Death Differ. 2012, 19, 321–332. [Google Scholar] [CrossRef]
  111. Rippmann, J.F.; Damm, K.; Schnapp, A. Functional Characterization of the Poly(ADP-Ribose) Polymerase Activity of Tankyrase 1, a Potential Regulator of Telomere Length. J. Mol. Biol. 2002, 323, 217–224. [Google Scholar] [CrossRef]
  112. Donigian, J.R.; de Lange, T. The Role of the Poly(ADP-Ribose) Polymerase Tankyrase1 in Telomere Length Control by the TRF1 Component of the Shelterin Complex. J. Biol. Chem. 2007, 282, 22662–22667. [Google Scholar] [CrossRef] [Green Version]
  113. Hsiao, S.J.; Poitras, M.F.; Cook, B.D.; Liu, Y.; Smith, S. Tankyrase 2 Poly(ADP-Ribose) Polymerase Domain-Deleted Mice Exhibit Growth Defects but Have Normal Telomere Length and Capping. Mol. Cell. Biol. 2006, 26, 2044–2054. [Google Scholar] [CrossRef] [Green Version]
  114. Seimiya, H.; Muramatsu, Y.; Ohishi, T.; Tsuruo, T. Tankyrase 1 as a Target for Telomere-Directed Molecular Cancer Therapeutics. Cancer Cell 2005, 7, 25–37. [Google Scholar] [CrossRef] [Green Version]
  115. Thorvaldsen, T.E. Targeting Tankyrase to Fight WNT-Dependent Tumours. Basic Clin. Pharmacol. Toxicol. 2017, 121, 81–88. [Google Scholar] [CrossRef]
Figure 1. Tankyrase domain organization, function and binding partners. (a) Schematic representation of the domain architecture from human TNKS1 and TNKS2. The sequence identities between domains of TNKS1 and TNKS2 are shown in percentages. TNKS1 has a Histidine, Proline and Serine-rich region (HPS) at the N-terminus, which is not present in TNKS2. (b) A model of TNKS function. TNKSs form multimers mediated by the SAM domain and interact with proteins that contain TNKS binding motifs (TBMs), which bind to the ARC1, 2, 4 and 5 domains of TNKSs. Many binding partners are PARylated by TNKSs, which leads to their subsequent ubiquitination and proteasomal degradation. Additionally, TNKSs may also act as scaffolding proteins by providing a platform for protein–protein interactions. (c) Examples of tankyrase binding partners TRF1, Axin1, GMD and USP25. The TBM sequence for each binding partner is shown. Crucial residues for the interaction with TNKS ARC domains are highlighted in blue. Previously reported presence or absence of the TBM in orthologs from other species is shown for TRF1 [23], Axin1 [10], GMD [24] and USP25 [25].
Figure 1. Tankyrase domain organization, function and binding partners. (a) Schematic representation of the domain architecture from human TNKS1 and TNKS2. The sequence identities between domains of TNKS1 and TNKS2 are shown in percentages. TNKS1 has a Histidine, Proline and Serine-rich region (HPS) at the N-terminus, which is not present in TNKS2. (b) A model of TNKS function. TNKSs form multimers mediated by the SAM domain and interact with proteins that contain TNKS binding motifs (TBMs), which bind to the ARC1, 2, 4 and 5 domains of TNKSs. Many binding partners are PARylated by TNKSs, which leads to their subsequent ubiquitination and proteasomal degradation. Additionally, TNKSs may also act as scaffolding proteins by providing a platform for protein–protein interactions. (c) Examples of tankyrase binding partners TRF1, Axin1, GMD and USP25. The TBM sequence for each binding partner is shown. Crucial residues for the interaction with TNKS ARC domains are highlighted in blue. Previously reported presence or absence of the TBM in orthologs from other species is shown for TRF1 [23], Axin1 [10], GMD [24] and USP25 [25].
Biomolecules 12 01688 g001
Figure 2. Distribution and possible origin of tankyrases. (a) The presence or absence of TNKSs in phylogenetic groups was determined from searches of various protein and genomic databases. Phylogenetic relationships among groups are schematically depicted. Representative species are shown for metazoans and selected groups and the presence of a single TNKS or both TNKS1/2 paralogs is indicated. Information with common names for the phylogenetic groups and species displayed here together with identifiers for the protein sequences can be found in Table S1. Species in which TNKSs were identified only from TSA data are marked (*). LUCA = Last universal common ancestor. (b) Phylogenetic reconstruction from TNKS sequences using the maximum likelihood method. Percentages for bootstrap support values are shown for each node. Bootstrap tests were conducted with 1000 bootstrap replicates. The analysis was performed with truncated TNKS sequences starting from the ARC3 domain to match the likely incomplete sequence from S. macrocollata. In Figure S1, the reconstruction was alternatively performed with the neighbor joining method. (c) Schematic representation of proteins with domain architectures related to TNKSs. Sequence identity percentages to the ARC5 or ART domain of human TNKS1 were derived from local sequence alignments. Shank1, ANKS1b and ASZ1 are human proteins that contain both ANK and SAM domains. The ANK protein from Acidianus manzaensis (AmANK) shares a high sequence identity with ARC5 of TNKS1. The number of ankyrin repeat units is indicated for each ANK domain. PARP7, PARP10 and PARP15 are human representatives from the clade 3 of the ARTD family, which is most closely related to clade 4 (TNKSs) [52]. The TNKS-like protein from the amoeba N. fowleri (NfTNKSL) is described in more detail in Figure S3.
Figure 2. Distribution and possible origin of tankyrases. (a) The presence or absence of TNKSs in phylogenetic groups was determined from searches of various protein and genomic databases. Phylogenetic relationships among groups are schematically depicted. Representative species are shown for metazoans and selected groups and the presence of a single TNKS or both TNKS1/2 paralogs is indicated. Information with common names for the phylogenetic groups and species displayed here together with identifiers for the protein sequences can be found in Table S1. Species in which TNKSs were identified only from TSA data are marked (*). LUCA = Last universal common ancestor. (b) Phylogenetic reconstruction from TNKS sequences using the maximum likelihood method. Percentages for bootstrap support values are shown for each node. Bootstrap tests were conducted with 1000 bootstrap replicates. The analysis was performed with truncated TNKS sequences starting from the ARC3 domain to match the likely incomplete sequence from S. macrocollata. In Figure S1, the reconstruction was alternatively performed with the neighbor joining method. (c) Schematic representation of proteins with domain architectures related to TNKSs. Sequence identity percentages to the ARC5 or ART domain of human TNKS1 were derived from local sequence alignments. Shank1, ANKS1b and ASZ1 are human proteins that contain both ANK and SAM domains. The ANK protein from Acidianus manzaensis (AmANK) shares a high sequence identity with ARC5 of TNKS1. The number of ankyrin repeat units is indicated for each ANK domain. PARP7, PARP10 and PARP15 are human representatives from the clade 3 of the ARTD family, which is most closely related to clade 4 (TNKSs) [52]. The TNKS-like protein from the amoeba N. fowleri (NfTNKSL) is described in more detail in Figure S3.
Biomolecules 12 01688 g002
Figure 3. Sequence conservation of tankyrases and representative PARP proteins. For the different proteins analyzed, the graph shows the sequence identity shared by the human ortholog with orthologs from different metazoan species. The identities were determined from multiple sequence alignments of the orthologs from each set of ARTD family members. In species that only encode a single TNKS paralog, the sequence identity shared with human TNKS1 or TNKS2 was determined against the same TNKS ortholog. The conserved proteins cytochrome C and seryl-tRNA synthase were used for comparison.
Figure 3. Sequence conservation of tankyrases and representative PARP proteins. For the different proteins analyzed, the graph shows the sequence identity shared by the human ortholog with orthologs from different metazoan species. The identities were determined from multiple sequence alignments of the orthologs from each set of ARTD family members. In species that only encode a single TNKS paralog, the sequence identity shared with human TNKS1 or TNKS2 was determined against the same TNKS ortholog. The conserved proteins cytochrome C and seryl-tRNA synthase were used for comparison.
Biomolecules 12 01688 g003
Figure 4. Conservation of tankyrase residues linked to the protein structure. (a) Residue-specific conservation plot. The number of different residue identities and ConSurf grade was plotted for each residue from an alignment of 113 TNKS sequences. The ConSurf grade is a normalized measure for the conservation and ranges from 1 (low conservation) to 9 (high conservation). The three likely flexible linker regions are designated L1, L2 and L3. The projected domain boundaries were assigned based on existing structural information of individual domains from human TNKS1 and TNKS2. A scale corresponding to 100 amino acid residues is shown. (b) A representative section of a multiple sequence alignment from the TNKS sequences of 11 metazoan species. This section covers a part of the ART domain sequence. Residues highlighted in gray are different compared to human TNKS1. Only residues identical in every TNKS sequence were mapped to the structures and are shown in green: (c) Structure of the ARC5 domain from human TNKS1 in complex with the TBM peptide from USP25 (orange) (PDB: 5GP7, [25]). (d) Structure of the SAM domain from human TNKS1 (PDB: 5KNI, [20]). (e) Structure of the ART domain from human TNKS2 in complex with the NAD+-mimicking compound EB-47 (orange) (PDB: 4BJ9, [68]).
Figure 4. Conservation of tankyrase residues linked to the protein structure. (a) Residue-specific conservation plot. The number of different residue identities and ConSurf grade was plotted for each residue from an alignment of 113 TNKS sequences. The ConSurf grade is a normalized measure for the conservation and ranges from 1 (low conservation) to 9 (high conservation). The three likely flexible linker regions are designated L1, L2 and L3. The projected domain boundaries were assigned based on existing structural information of individual domains from human TNKS1 and TNKS2. A scale corresponding to 100 amino acid residues is shown. (b) A representative section of a multiple sequence alignment from the TNKS sequences of 11 metazoan species. This section covers a part of the ART domain sequence. Residues highlighted in gray are different compared to human TNKS1. Only residues identical in every TNKS sequence were mapped to the structures and are shown in green: (c) Structure of the ARC5 domain from human TNKS1 in complex with the TBM peptide from USP25 (orange) (PDB: 5GP7, [25]). (d) Structure of the SAM domain from human TNKS1 (PDB: 5KNI, [20]). (e) Structure of the ART domain from human TNKS2 in complex with the NAD+-mimicking compound EB-47 (orange) (PDB: 4BJ9, [68]).
Biomolecules 12 01688 g004
Figure 5. Experimental validation of the tankyrase molecular functions from Amphimedon queenslandica. (a) Western blot of ART and SAM-ART constructs from human TNKS1 and A. queenslandica TNKS (AqTNKS). SAM-ART constructs contain an N-terminal maltose-binding protein (MBP) tag. TNKS constructs (2 μM) were incubated for 1 h at room temperature in presence or absence of a mixture containing biotin-NAD+/NAD+ (1:10). The total NAD+ concentration was 2 μM for TNKS1 and 10 μM for AqTNKS due to a lower activity observed in AqTNKS. PAR with incorporated biotin was detected with streptavidin conjugated to horseradish peroxidase (SA-HRP), while total PAR modification independent of biotin was detected using nanoluciferase fused to ALC1 (Nluc-ALC1). PARG (0.2 mg/mL) was added to verify a removal of PAR chains. (b) Thermal stabilization of TNKS ART domain constructs by selective TNKS inhibitors. Differential scanning fluorimetry (DSF) was performed after mixing TNKS ART domain constructs (5 μM) and TNKS inhibitors (25 μM). A control without inhibitor was included. The absolute melting temperatures (Tm) are shown on the left and data shown are mean ± standard deviation with number of replicates n = 4; the melting temperature relative to the respective controls calculated from the mean values are shown on the right. (c) FRET-based determination of the dissociation constant KD for ARC5 constructs with a TBM peptide. CFP-fused ARC5 constructs (100 nM) from TNKS1 or AqTNKS were mixed with different concentrations of YFP fused to the TBM peptide (REAGDGEE). As control, the CFP-ARC5 constructs were mixed with YFP without the TBM peptide. The fluorescence emissions from FRET (EmFRET) were determined by a method described by Song et al. [51]. Data shown are mean ± standard deviation with number of replicates n = 4.
Figure 5. Experimental validation of the tankyrase molecular functions from Amphimedon queenslandica. (a) Western blot of ART and SAM-ART constructs from human TNKS1 and A. queenslandica TNKS (AqTNKS). SAM-ART constructs contain an N-terminal maltose-binding protein (MBP) tag. TNKS constructs (2 μM) were incubated for 1 h at room temperature in presence or absence of a mixture containing biotin-NAD+/NAD+ (1:10). The total NAD+ concentration was 2 μM for TNKS1 and 10 μM for AqTNKS due to a lower activity observed in AqTNKS. PAR with incorporated biotin was detected with streptavidin conjugated to horseradish peroxidase (SA-HRP), while total PAR modification independent of biotin was detected using nanoluciferase fused to ALC1 (Nluc-ALC1). PARG (0.2 mg/mL) was added to verify a removal of PAR chains. (b) Thermal stabilization of TNKS ART domain constructs by selective TNKS inhibitors. Differential scanning fluorimetry (DSF) was performed after mixing TNKS ART domain constructs (5 μM) and TNKS inhibitors (25 μM). A control without inhibitor was included. The absolute melting temperatures (Tm) are shown on the left and data shown are mean ± standard deviation with number of replicates n = 4; the melting temperature relative to the respective controls calculated from the mean values are shown on the right. (c) FRET-based determination of the dissociation constant KD for ARC5 constructs with a TBM peptide. CFP-fused ARC5 constructs (100 nM) from TNKS1 or AqTNKS were mixed with different concentrations of YFP fused to the TBM peptide (REAGDGEE). As control, the CFP-ARC5 constructs were mixed with YFP without the TBM peptide. The fluorescence emissions from FRET (EmFRET) were determined by a method described by Song et al. [51]. Data shown are mean ± standard deviation with number of replicates n = 4.
Biomolecules 12 01688 g005
Figure 6. Identification of corresponding tankyrase binding motifs in orthologs of confirmed interaction partners. Orthologs of 21 TNKS interaction partners (TIPs) were examined for presence or absence of possible TBM sequences. Three configurations Rxx[ACGP]xGxx, Rx(4)Gxx and Rx(5)Gxx were considered as TBM. TBMs were only considered if the motif aligned with the TBM in the human ortholog in a multiple sequence alignment, although in some cases also TBMs that did not align but were found in nearby positions are shown (see legend). * The non-canonical second motif from Axin1 was not included in this analysis.
Figure 6. Identification of corresponding tankyrase binding motifs in orthologs of confirmed interaction partners. Orthologs of 21 TNKS interaction partners (TIPs) were examined for presence or absence of possible TBM sequences. Three configurations Rxx[ACGP]xGxx, Rx(4)Gxx and Rx(5)Gxx were considered as TBM. TBMs were only considered if the motif aligned with the TBM in the human ortholog in a multiple sequence alignment, although in some cases also TBMs that did not align but were found in nearby positions are shown (see legend). * The non-canonical second motif from Axin1 was not included in this analysis.
Biomolecules 12 01688 g006
Table 1. Species frequently referred to in this work.
Table 1. Species frequently referred to in this work.
SpeciesCommon NameGroupNotes
Homo sapiensModern humanApesMost studies of TNKS and its interaction partners were carried out in human model systems.
Pongo abeliiSumatran orangutanApesExtant ape species that is phylogenetically most distant from human.
Mus musculusHouse mouseMammalsMost widely used mammalian model organism.
Gallus gallusRed junglefowlBirdsIncludes the domesticated chicken.
Xenopus laevisAfrican clawed frogAmphibiansModel amphibian.
Salmo salarAtlantic salmonBony fishesRepresentative of fishes with high-quality genome sequence data.
Petromyzon marinusSea lampreyJawless fishesModel for vertebrate evolution. An early diverging vertebrate with only one TNKS instead of two paralogs (TNKS1/2).
Branchiostoma belcheriBelcher’s lanceletLanceletsClosely related to vertebrates. Model organism for vertebrate evolution.
Drosophila melanogasterFruit flyInsectsPopular model organism in genetics and developmental biology.
Caenorhabiditis elegansC. elegansRoundwormsPopular model organism for studying neuronal development. Gene-loss of TNKS and other PARPs was reported.
Nematostella vectensisStarlet sea anemoneCnidariansModel cnidarian.
Hydra vulgarisFresh-water polypCnidariansModel organism in regenerative biology.
Amphimedon queenslandicaN/ASpongesModel for metazoan evolution and developmental biology. Functions from its TNKS ortholog were experimentally characterized in this work.
Salpingoeca kvevriiN/AChoanoflagellatesThe most diverged TNKS orthologs identified in this work are from these choanoflagellates.
Salpingoeca macrocollata
Naegleria fowleriBrain-eating amoebaPercolozoaNon-metazoan eukaryote that encodes an ART domain with high similarity to that of TNKSs.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sowa, S.T.; Bosetti, C.; Galera-Prat, A.; Johnson, M.S.; Lehtiö, L. An Evolutionary Perspective on the Origin, Conservation and Binding Partner Acquisition of Tankyrases. Biomolecules 2022, 12, 1688. https://doi.org/10.3390/biom12111688

AMA Style

Sowa ST, Bosetti C, Galera-Prat A, Johnson MS, Lehtiö L. An Evolutionary Perspective on the Origin, Conservation and Binding Partner Acquisition of Tankyrases. Biomolecules. 2022; 12(11):1688. https://doi.org/10.3390/biom12111688

Chicago/Turabian Style

Sowa, Sven T., Chiara Bosetti, Albert Galera-Prat, Mark S. Johnson, and Lari Lehtiö. 2022. "An Evolutionary Perspective on the Origin, Conservation and Binding Partner Acquisition of Tankyrases" Biomolecules 12, no. 11: 1688. https://doi.org/10.3390/biom12111688

APA Style

Sowa, S. T., Bosetti, C., Galera-Prat, A., Johnson, M. S., & Lehtiö, L. (2022). An Evolutionary Perspective on the Origin, Conservation and Binding Partner Acquisition of Tankyrases. Biomolecules, 12(11), 1688. https://doi.org/10.3390/biom12111688

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